-
Improving methods for determining the subsurface structure of sunspots from
their seismic signature requires a better understanding of the interaction of
waves with magnetic field concentrations. We aim to quantify the impact of
changes in the internal structure of sunspots on local helioseismic signals. We
have numerically simulated the propagation of a stochastic wave field through
sunspot models with different properties, accounting for changes in the Wilson
depression between 250 and 550 km and in the photospheric umbral magnetic field
between 1500 and 3500 G. The results show that travel-time shifts at
frequencies above approximately 3.50 mHz (depending on the phase-speed filter)
are insensitive to the magnetic field strength. The travel time of these waves
is determined exclusively by the Wilson depression and sound-speed
perturbation. The travel time of waves with lower frequencies is affected by
the direct effect of the magnetic field, although photospheric field strengths
below 1500 G do not leave a significant trace on the travel-time measurements.
These results could potentially be used to develop simplified travel-time
inversion methods.
-
Context. Local helioseismology has detected spatially extended converging
surface flows into solar active regions. These play an important role in
flux-transport models of the solar dynamo.
Aims. We aim to validate the existence of the inflows by deriving horizontal
flow velocities around active regions with local correlation tracking of
granulation.
Methods. We generate a six-year long-time series of full-disk maps of the
horizontal velocity at the solar surface by tracking granules in continuum
intensity images provided by the Helioseismic and Magnetic Imager (HMI) onboard
the Solar Dynamics Observatory (SDO).
Results. On average, active regions are surrounded by inflows extending up to
10 deg from the center of the active region of magnitudes of 20-30 m/s,
reaching locally up to 40 m/s, which is in agreement with results from local
helioseismology. By computing an ensemble average consisting of 243 individual
active regions, we show that the inflows are not azimuthally symmetric but
converge predominantly towards the trailing polarity of the active region with
respect to the longitudinally and temporally averaged flow field.
-
While the velocity fluctuations of supergranulation dominate the spectrum of
solar convection at the solar surface, very little is known about the
fluctuations in other physical quantities like temperature or density at
supergranulation scale. Using SDO/HMI observations, we characterize the
intensity contrast of solar supergranulation at the solar surface. We identify
the positions of ${\sim}10^4$ outflow and inflow regions at supergranulation
scales, from which we construct average flow maps and co-aligned intensity and
magnetic field maps. In the average outflow center, the maximum intensity
contrast is $(7.8\pm0.6)\times10^{-4}$ (there is no corresponding feature in
the line-of-sight magnetic field). This corresponds to a temperature
perturbation of about $1.1\pm0.1$ K, in agreement with previous studies. We
discover an east-west anisotropy, with a slightly deeper intensity minimum east
of the outflow center. The evolution is asymmetric in time: the intensity
excess is larger 8 hours before the reference time (the time of maximum
outflow), while it has almost disappeared 8 hours after the reference time. In
the average inflow region, the intensity contrast mostly follows the magnetic
field distribution, except for an east-west anisotropic component that
dominates 8 hours before the reference time. We suggest that the east-west
anisotropy in the intensity is related to the wave-like properties of
supergranulation.
-
Wave propagation through sunspots involves conversion between waves of
acoustic and magnetic character. In addition, the thermal structure of sunspots
is very different than that of the quiet Sun. As a consequence, the
interpretation of local helioseismic measurements of sunspots has long been a
challenge. With the aim of understanding these measurements, we carry out
numerical simulations of wave propagation through sunspots. Helioseismic
holography measurements made from the resulting simulated wavefields show
qualitative agreement with observations of real sunspots. We use additional
numerical experiments to determine, separately, the influence of the thermal
structure of the sunspot and the direct effect of the sunspot magnetic field.
We use the ray approximation to show that the travel-time shifts in the thermal
(non-magnetic) sunspot model are primarily produced by changes in the wave path
due to the Wilson depression rather than variations in the wave speed. This
shows that inversions for the subsurface structure of sunspots must account for
local changes in the density. In some ranges of horizontal phase speed and
frequency there is agreement (within the noise level in the simulations)
between the travel times measured in the full magnetic sunspot model and the
thermal model. If this conclusion proves to be robust for a wide range of
models, it would suggest a path towards inversions for sunspot structure.
-
Context. Local correlation tracking of granulation (LCT) is an important
method for measuring horizontal flows in the photosphere. This method exhibits
a systematic error that looks like a flow converging towards disk center, also
known as the shrinking-Sun effect.
Aims. We aim at studying the nature of the shrinking-Sun effect for continuum
intensity data and at deriving a simple model that can explain its origin.
Methods. We derived LCT flow maps by running the local correlation tracking
code FLCT on tracked and remapped continuum intensity maps provided by the
Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory.
We also computed flow maps from synthetic continuum images generated from
STAGGER code simulations of solar surface convection. We investigated the
origin of the shrinking-Sun effect by generating an average granule from
synthetic data from the simulations.
Results. The LCT flow maps derived from HMI and from the simulations exhibit
a shrinking-Sun effect of comparable magnitude. The origin of this effect is
related to the apparent asymmetry of granulation originating from radiative
transfer effects when observing with a viewing angle inclined from vertical.
This causes, in combination with the expansion of the granules, an apparent
motion towards disk center.
-
Flow vorticity is a fundamental property of turbulent convection in rotating
systems. Solar supergranules exhibit a preferred sense of rotation, which
depends on the hemisphere. This is due to the Coriolis force acting on the
diverging horizontal flows. We aim to spatially resolve the vertical flow
vorticity of the average supergranule at different latitudes, both for outflow
and inflow regions. To measure the vertical vorticity, we use two independent
techniques: time-distance helioseismology (TD) and local correlation tracking
of granules in intensity images (LCT) using data from the Helioseismic and
Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Both maps
are corrected for center-to-limb systematic errors. We find that 8-h TD and LCT
maps of vertical vorticity are highly correlated at large spatial scales.
Associated with the average supergranule outflow, we find tangential (vortical)
flows that reach about 10 m/s in the clockwise direction at 40{\deg} latitude.
In average inflow regions, the tangential flow reaches the same magnitude, but
in the anti-clockwise direction. These tangential velocities are much smaller
than the radial (diverging) flow component (300 m/s for the average outflow and
200 m/s for the average inflow). The results for TD and LCT as measured from
HMI are in excellent agreement for latitudes between $-$60{\deg} and 60{\deg}.
From HMI LCT, we measure the vorticity peak of the average supergranule to have
a full width at half maximum of about 13 Mm for outflows and 8 Mm for inflows.
This is larger than the spatial resolution of the LCT measurements (about 3
Mm). On the other hand, the vorticity peak in outflows is about half the value
measured at inflows (e.g. 4/(10^6 s) clockwise compared to 8/(10^6 s)
anti-clockwise at 40{\deg} latitude). Results from MDI/SOHO obtained in 2010
are biased compared to the HMI/SDO results for the same period.
-
Supergranules in the quiet Sun are outlined by a web-like structure of
enhanced magnetic field strength, the so-called magnetic network. We aim to map
the magnetic network field around the average supergranule near disk center. We
use observations of the line-of-sight component of the magnetic field from the
Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory
(SDO). The average supergranule is constructed by coaligning and averaging over
3000 individual supergranules. We determine the positions of the supergranules
with an image segmentation algorithm that we apply on maps of the horizontal
flow divergence measured using time-distance helioseismology. In the center of
the average supergranule the magnetic (intranetwork) field is weaker by about
2.2 Gauss than the background value (3.5 Gauss), whereas it is enhanced in the
surrounding ring of horizontal inflows (by about 0.6 Gauss on average). We find
that this network field is significantly stronger west (prograde) of the
average supergranule than in the east (by about 0.3 Gauss). With time-distance
helioseismology, we find a similar anisotropy. The observed anisotropy of the
magnetic field adds to the mysterious dynamical properties of solar
supergranulation.
-
The helioseismic properties of the wave scattering generated by monolithic
and spaghetti sunspots are analyzed by means of numerical simulations. In these
computations, an incident f or p1 mode travels through the sunspot model, which
produces absorption and phase shift of the waves. The scattering is studied by
inspecting the wavefield, computing travel-time shifts, and performing
Fourier-Hankel analysis. The comparison between the results obtained for both
sunspot models reveals that the differences in the absorption coefficient can
be detected above noise level. The spaghetti model produces an steep increase
of the phase shift with the degree of the mode at short wavelengths, while
mode-mixing is more efficient for the monolithic model. These results provide a
clue for what to look for in solar observations to discern the constitution of
sunspots between the proposed monolithic and spaghetti models.
-
PLATO 2.0 has recently been selected for ESA's M3 launch opportunity
(2022/24). Providing accurate key planet parameters (radius, mass, density and
age) in statistical numbers, it addresses fundamental questions such as: How do
planetary systems form and evolve? Are there other systems with planets like
ours, including potentially habitable planets? The PLATO 2.0 instrument
consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2
with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large
photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars
in wide fields to detect and characterize planets down to Earth-size by
photometric transits, whose masses can then be determined by ground-based
radial-velocity follow-up measurements. Asteroseismology will be performed for
these bright stars to obtain highly accurate stellar parameters, including
masses and ages. The combination of bright targets and asteroseismology results
in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet
radii, masses and ages, respectively. The planned baseline observing strategy
includes two long pointings (2-3 years) to detect and bulk characterize planets
reaching into the habitable zone (HZ) of solar-like stars and an additional
step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will
observe up to 1,000,000 stars and detect and characterize hundreds of small
planets, and thousands of planets in the Neptune to gas giant regime out to the
HZ. It will therefore provide the first large-scale catalogue of bulk
characterized planets with accurate radii, masses, mean densities and ages.
This catalogue will include terrestrial planets at intermediate orbital
distances, where surface temperatures are moderate. Coverage of this parameter
range with statistical numbers of bulk characterized planets is unique to PLATO
2.0.
-
We comment on the potential for continuing asteroseismology of solar-type and
red-giant stars in a 2-wheel Kepler Mission. Our main conclusion is that by
targeting stars in the ecliptic it should be possible to perform high-quality
asteroseismology, as long as favorable scenarios for 2-wheel pointing
performance are met. Targeting the ecliptic would potentially facilitate unique
science that was not possible in the nominal Mission, notably from the study of
clusters that are significantly brighter than those in the Kepler field. Our
conclusions are based on predictions of 2-wheel observations made by a space
photometry simulator, with information provided by the Kepler Project used as
input to describe the degraded pointing scenarios. We find that elevated levels
of frequency-dependent noise, consistent with the above scenarios, would have a
significant negative impact on our ability to continue asteroseismic studies of
solar-like oscillators in the Kepler field. However, the situation may be much
more optimistic for observations in the ecliptic, provided that pointing resets
of the spacecraft during regular desaturations of the two functioning reaction
wheels are accurate at the < 1 arcsec level. This would make it possible to
apply a post-hoc analysis that would recover most of the lost photometric
precision. Without this post-hoc correction---and the accurate re-pointing it
requires---the performance would probably be as poor as in the Kepler-field
case. Critical to our conclusions for both fields is the assumed level of
pointing noise (in the short-term jitter and the longer-term drift). We suggest
that further tests will be needed to clarify our results once more detail and
data on the expected pointing performance becomes available, and we offer our
assistance in this work.
-
We use numerial simulations to study the absorption and phase shift of
surface-gravity waves caused by groups of magnetic flux tubes. The dependence
of the scattering coefficients with the distance between the tubes and their
positions is analyzed for several cases with two or three flux tubes embedded
in a quiet Sun atmosphere. The results are compared with those obtained
neglecting completely or partially multiple scattering effects. We show that
multiple scattering has a significant impact on the absorption measurements and
tends to reduce the phase shift. We also consider more general cases of
ensembles of randomly distributed flux tubes, and we have evaluated the effects
on the scattering measurements of changing the number of tubes included in the
bundle and the average distance between flux tubes. We find that for the
longest wavelength incoming waves multiple scattering enhances the absorption,
and its efficiency increases with the number of flux tubes and the reduction of
the distance between them.
-
The scattering of f-modes by magnetic tubes is analyzed using
three-dimensional numerical simulations. An f-mode wave packet is propagated
through a solar atmosphere embedded with three different flux tube models which
differ in radius and total magnetic flux. A quiet Sun simulation without a tube
present is also performed as a reference. Waves are excited inside the flux
tube and propagate along the field lines, and jacket modes are generated in the
surroundings of the flux tube, carrying 40% as much energy as the tube modes.
The resulting scattered wave is mainly an f-mode composed of a mixture of m=0
and m=+/-1 modes. The amplitude of the scattered wave approximately scales with
the magnetic flux. A small amount of power is scattered into the p_1-mode. We
have evaluated the absorption and phase shift from a Fourier-Hankel
decomposition of the photospheric vertical velocities. They are compared with
the results obtained from the emsemble average of 3400 small magnetic elements
observed in high-resolution MDI Doppler datacubes. The comparison shows that
the observed dependence of the phase shift with wavenumber can be matched
reasonably well with the simulated flux tube model. The observed variation of
the phase-shifts with the azimuthal order $m$ appears to depend on details of
the ensemble averaging, including possible motions of the magnetic elements and
asymmetrically shaped elements.
-
While sunspots are easily observed at the solar surface, determining their
subsurface structure is not trivial. There are two main hypotheses for the
subsurface structure of sunspots: the monolithic model and the cluster model.
Local helioseismology is the only means by which we can investigate
subphotospheric structure. However, as current linear inversion techniques do
not yet allow helioseismology to probe the internal structure with sufficient
confidence to distinguish between the monolith and cluster models, the
development of physically realistic sunspot models are a priority for
helioseismologists. This is because they are not only important indicators of
the variety of physical effects that may influence helioseismic inferences in
active regions, but they also enable detailed assessments of the validity of
helioseismic interpretations through numerical forward modeling. In this paper,
we provide a critical review of the existing sunspot models and an overview of
numerical methods employed to model wave propagation through model sunspots. We
then carry out an helioseismic analysis of the sunspot in Active Region 9787
and address the serious inconsistencies uncovered by
\citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find
that this sunspot is most probably associated with a shallow, positive
wave-speed perturbation (unlike the traditional two-layer model) and that
travel-time measurements are consistent with a horizontal outflow in the
surrounding moat.
-
We perform helioseismic holography to assess the noise in p-mode travel-time
shifts which would form the basis of inferences of large-scale flows throughout
the solar convection zone. We also derive the expected travel times from a
parameterized return (equatorward) flow component of the meridional circulation
at the base of the convection zone from forward models under the assumption of
the ray and Born approximations. From estimates of the signal-to-noise ratio
for measurements focused near the base of the convection zone, we conclude that
the helioseismic detection of the deep meridional flow including the return
component may not be possible using data spanning an interval less than a solar
cycle.
-
With the aim of studying magnetic effects in time-distance helioseismology,
we use the first-order Born approximation to compute the scattering of acoustic
plane waves by a magnetic cylinder embedded in a uniform medium. We show, by
comparison with the exact solution, that the travel-time shifts computed in the
Born approximation are everywhere valid to first order in the ratio of the
magnetic to the gas pressures. We also show that, for arbitrary magnetic field
strength, the Born approximation is not valid in the limit where the radius of
the magnetic cylinder tends to zero.
-
Wave travel-time shifts in the vicinity of sunspots are typically interpreted
as arising predominantly from magnetic fields, flows, and local changes in
sound speed. We show here that the suppression of granulation related wave
sources in a sunspot can also contribute significantly to these travel-time
shifts, and in some cases, an asymmetry between in and outgoing wave travel
times. The tight connection between the physical interpretation of travel times
and source-distribution homogeneity is confirmed. Statistically significant
travel-time shifts are recovered upon numerically simulating wave propagation
in the presence of a localized decrease in source strength. We also demonstrate
that these time shifts are relatively sensitive to the modal damping rates;
thus we are only able to place bounds on the magnitude of this effect. We see a
systematic reduction of 10-15 seconds in $p$-mode mean travel times at short
distances ($\sim 6.2$ Mm) that could be misinterpreted as arising from a
shallow (thickness of 1.5 Mm) increase ($\sim$ 4%) in the sound speed. At
larger travel distances ($\sim 24$ Mm) a 6-13 s difference between the ingoing
and outgoing wave travel times is observed; this could mistakenly be
interpreted as being caused by flows.
-
We calculate the scattering effects associated with the interaction of a
surface gravity or {\it f} mode with a thin magnetic flux tube embedded in a
realistically stratified medium. We find that the dominant scattered wave is an
$f$ mode with amplitude and phase of 1.17% and around $49^\circ$ relative to
the incident wave, compared to the values of 0.13% and $40^\circ$ estimated
from observations. The extent of scattering into high-order acoustic $p$ modes
is too weak to be accurately characterized. We recover the result that the
degree of scattering is enhanced as (a) the frequency of the incident wave
increases and (b) the flux tube becomes magnetically dominated.
-
We perform helioseismic holography on realistic solar convection simulations
and compare the observed travel-time perturbations with the expected travel
times from the horizontal flows in the simulations computed from forward models
under the assumption of the Born approximation. We demonstrate reasonable
agreement between the observed and model travel times which reinforces the
validity of helioseismic holography in the detection of subsurface horizontal
flows. From the variation of the signal-to-noise ratio with depth, we conclude
that the helioseismic detection of individual flow structures with spatial
scales of supergranulation or smaller is not possible for depths below about 5
Mm below the surface over time scales less than a day. Approximately half of
the observed signal originates within the first 2 Mm below the surface. A
consequence of this is a rapid decrease (and reversal in some cases) of the
travel-time perturbations with depth due to the contribution to the
measurements of oppositely directed surface flows in neighboring convective
cells. This confirms an earlier interpretation of similar effects reported from
observations.