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.