The well-observed acoustic halo is an enhancement in time-averaged Doppler velocity and intensity power with respect to quiet-sun values which is prominent for weak and highly inclined field around the penumbra of sunspots and active regions. We perform 3D linear wave modelling with realistic distributed acoustic sources in a MHS sunspot atmosphere and compare the resultant simulation enhancements with multi-height SDO observations of the phenomenon. We find that simulated halos are in good qualitative agreement with observations. We also provide further proof that the underlying process responsible for the halo is the refraction and return of fast magnetic waves which have undergone mode conversion at the critical $a=c$ atmospheric layer. In addition, we also find strong evidence that fast-Alfv\'en mode conversion plays a significant role in the structure of the halo, taking energy away from photospheric and chromospheric heights in the form of field-aligned Alfv\'en waves. This conversion process may explain the observed "dual-ring" halo structure at higher ($> 8 $ mHz) frequencies.

A crucial feature not widely accounted for in local helioseismology is that surface magnetic regions actually open a window from the interior into the solar atmosphere, and that the seismic waves leak through this window, reflect high in the atmosphere, and then re-enter the interior to rejoin the seismic wave field normally confined there. In a series of recent numerical studies using translation invariant atmospheres, we utilised a "directional time-distance helioseismology" measurement scheme to study the implications of the returning fast and Alfv\'en waves higher up in the solar atmosphere on the seismology at the photosphere (Cally & Moradi 2013; Moradi & Cally 2014). In this study, we extend our directional time-distance analysis to more realistic sunspot-like atmospheres to better understand the direct effects of the magnetic field on helioseismic travel-time measurements in sunspots. In line with our previous findings, we uncover a distinct frequency-dependant directional behaviour in the travel-time measurements, consistent with the signatures of MHD mode conversion. We found this to be the case regardless of the sunspot field strength or depth of its Wilson depression. We also isolated and analysed the direct contribution from purely thermal perturbations to the measured travel times, finding that waves propagating in the umbra are much more sensitive to the underlying thermal effects of the sunspot.

An enhancement in high-frequency acoustic power is commonly observed in the solar photosphere and chromosphere surrounding magnetic active regions. We perform 3D linear forward wave modelling with a simple wavelet pulse acoustic source to ascertain whether the formation of the acoustic halo is caused by MHD mode conversion through regions of moderate and inclined magnetic fields. This conversion type is most efficient when high frequency waves from below intersect magnetic field lines at a large angle. We find a strong relationship between halo formation and the equipartition surface at which the Alfv\'en speed $a$ matches the sound speed $c$, lending support to the theory that photospheric and chromospheric halo enhancement is due to the creation and subsequent reflection of magnetically dominated fast waves from essentially acoustic waves as they cross $a=c$. In simulations where we have capped $a$ such that waves are not permitted to refract after reaching the $a=c$ height, halos are non-existent, which suggests that the power enhancement is wholly dependent on returning fast waves. We also reproduce some of the observed halo properties, such as a dual 6 and 8 mHz enhancement structure and a spatial spreading of the halo with height.

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.

The rapid exponential increase in the Alfv\'en wave speed with height above the solar surface presents a serious challenge to physical modelling of the effects of magnetic fields on solar oscillations, as it introduces a significant CFL time-step constraint for explicit numerical codes. A common approach adopted in computational helioseismology, where long simulations in excess of 10 hours (hundreds of wave periods) are often required, is to cap the Alfv\'en wave speed by artificially modifying the momentum equation when the ratio between Lorentz and hydrodynamic forces becomes too large. However, recent studies have demonstrated that the Alfv\'en wave speed plays a critical role in the MHD mode conversion process, particularly in determining the reflection height of the upward propagating helioseismic fast wave. Using numerical simulations of helioseismic wave propagation in constant inclined (relative to the vertical) magnetic fields we demonstrate that the imposition of such artificial limiters significantly affects time-distance travel times unless the Alfv\'en wave-speed cap is chosen comfortably in excess of the horizontal phase speeds under investigation.

Active regions are open wounds in the Sun's surface. Seismic oscillations from the interior pass through them into the atmosphere, changing their nature in the process to fast and slow magneto-acoustic waves. The fast waves then partially reflect and partially mode convert to upgoing and downgoing Alfv\'en waves. The reflected fast and downgoing Alfv\'en waves then re-enter the interior through the active regions that spawned them, infecting the surface seismology with signatures of the atmosphere. Using numerical simulations of waves in uniform magnetic fields, we calculate the upward acoustic and Alfv\'enic losses in the atmosphere as functions of field inclination and wave orientation as well as the Time-Distance `travel time' perturbations, and show that they are related. Travel time perturbations relative to quiet Sun can exceed 40 seconds in 1 kG magnetic field. It is concluded that active region seismology is indeed significantly infected by waves leaving and re-entering the interior through magnetic wounds, with differing travel times depending on the orientation of the wave vector relative to the magnetic field. This presages a new directional-time-distance seismology.

We successfully utilise the SPARC code to model fast-Alfv\'en mode conversion in the region $c_A \gg c_S$ via 3-D MHD numerical simulations of helioseismic waves within constant inclined magnetic field configurations. This was achieved only after empirically modifying the background density and gravitational stratifications in the upper layers of our computational box, as opposed to imposing a traditional Lorentz Force limiter, to ensure a manageable timestep. We found that the latter approach inhibits the fast-Alfv\'en mode conversion process by severely damping the magnetic flux above the surface.

Helioseismology is the study of the variations in the internal structure and properties of the dynamics of the Sun from measurements of its surface oscillations. With the 2010 launch of the Solar Dynamics Observatory (SDO) we are undoubtedly approaching a new dawn for local helioseismology, as the extent and quality of raw surface oscillation data has never been better. However, advances in theory and modelling are still required to fully utilise these data, especially in magnetic active regions and sunspots, where the physics is poorly understood.

In this paper, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science, 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose to make additional corrections to the sound speed of Model S before stabilisation. No stabilised model can be perfectly solar-like, so we present three stabilised models with slightly different eigenmodes. The models are appropriate to study the f and p1 to p4 modes with spherical harmonic degrees in the range from 400 to 900. Background model CSM has a modified pressure gradient for stabilisation and has eigenfrequencies within 2% of Model S. Model CSM_A has an additional 10% increase in sound speed in the top 1 Mm resulting in eigenfrequencies within 2% of Model S and eigenfunctions that are, in comparison with CSM, closest to those of Model S. Model CSM_B has a 3% decrease in sound speed in the top 5 Mm resulting in eigenfrequencies within 1% of Model S and eigenfunctions that are only marginally adversely affected. These models are useful to study the interaction of solar waves with embedded three-dimensional heterogeneities, such as convective flows and model sunspots. We have also calculated the response of the stabilised models to excitation by random near-surface sources, using simulations of the propagation of linear waves. We find that the simulated power spectra of wave motion are in good agreement with an observed SOHO/MDI power spectrum. Overall, our convectively stabilised background models provide a good basis for quantitative numerical local helioseismology. The models are available for download from http://www.mps.mpg.de/projects/seismo/NA4/.

In sequel to Moradi et al. [2009, ApJ, 690, L72], we employ two established numerical forward models (a 3D ideal MHD solver and MHD ray theory) in conjunction with time-distance helioseismology to probe the lateral extent of wave-speed perturbations produced in regions of strong, near-surface magnetic fields. We continue our comparisons of forward modeling approaches by extending our previous surface-focused travel-time measurements with a common midpoint deep-focusing scheme that avoids the use of oscillation signals within the sunspot region. The idea is to also test MHD ray theory for possible application in future inverse methods.

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.

A data and computation center for helioseismology has been set up at the Max Planck Institute for Solar System Research in Germany to prepare for the SDO mission. Here we present the system infrastructure and the scientific aims of this project, which is funded through grants from the German Aerospace Center and the European Research Council.

We investigate the direct contribution of strong, sunspot-like magnetic fields to helioseismic wave travel-time shifts via two numerical forward models, a 3D ideal MHD solver and MHD ray theory. The simulated data cubes are analyzed using the traditional time-distance center-to-annulus measurement technique. We also isolate and analyze the direct contribution from purely thermal perturbations to the observed travel-time shifts, confirming some existing ideas and bring forth new ones: (i) that the observed travel-time shifts in the vicinity of sunspots are largely governed by MHD physics, (ii) the travel-time shifts are sensitively dependent on frequency and phase-speed filter parameters and the background power below the $p_1$ ridge, and finally, (iii) despite its seeming limitations, ray theory succeeds in capturing the essence of the travel-time variations as derived from the MHD simulations.

On 2005 January 15, the active region AR10720 produced an X1.2 solar flare that induced high levels of seismicity into the photospheric layers. The seismic source was detected using helioseismic holography and analysed in detail in Paper I. Egression power maps at 6 mHz with a 2 mHz bandwidth revealed a compact acoustic source strongly correlated with the footpoints of the coronal loop that hosted the flare. We present a magneto-seismic study of this active region in order to understand, for the first time, the magnetic topological structure of a coronal field that hosts an acoustically active solar flare. The accompanying analysis attempts to answer questions such as: Can the magnetic field act as a barrier and prevent seismic waves from spreading away from the focus of the sunquake? And, what is the most efficient magnetic structure that would facilitate the development of a strong seismic source in the photosphere?

One of the main uncertainties associated with time-distance measurements is how to interpret observations made in regions of strong magnetic field as isolating the effects of the magnetic field from thermal or sound-speed perturbations has proved to be quite complex. One possible way to decouple these effects is by using a 3D sunspot model, with a surrounding stratified, quiet-Sun atmosphere to model the magneto-acoustic ray propagation, and analyze the resulting ray travel-time perturbations that will directly account for the magnetic field. These artificial travel-time perturbation profiles have provided us with several related but distinct observations: i) that strong surface magnetic fields have a dual effect on helioseismic rays -- increasing their skip distance while at the same time speeding them up considerably compared to their quiet-Sun counterparts, ii) there is a clear and significant frequency dependence of both skip-distance and travel-time perturbations across the simulated sunspot radius, iii) the negative sign and magnitude of these perturbations appears to be directly related to the sunspot magnetic-field strength and inclination, iv) by ``switching off'' the magnetic field inside the sunspot, we are able to completely isolate the thermal component of the travel-time perturbations observed, which is seen to be both opposite in sign and much smaller in magnitude than those measured when the magnetic field is present. These results tend to suggest that purely thermal perturbations are unlikely to be the main effect seen in travel times through sunspots and that strong, near-surface magnetic fields may be directly and significantly altering the magnitude and lateral extent of sound-speed inversions of sunspots made by time-distance helioseismology.

On 10 March 2001 the active region NOAA 9368 produced an unusually impulsive solar flare in close proximity to the solar limb. This flare has previously been studied in great detail, with observations classifying it as a type 1 white-light flare with a very hard spectrum in hard X-rays. The flare was also associated with a type II radio burst and coronal mass ejection. The flare emission characteristics appeared to closely correspond with previous instances of seismic emission from acoustically active flares. Using standard local helioseismic methods, we identified the seismic signatures produced by the flare that, to date, is the least energetic (in soft X-rays) of the flares known to have generated a detectable acoustic transient. Holographic analysis of the flare shows a compact acoustic source strongly correlated with the impulsive hard X-ray, visible continuum, and radio emission. Time-distance diagrams of the seismic waves emanating from the flare region also show faint signatures, mainly in the eastern sector of the active region. The strong spatial coincidence between the seismic source and the impulsive visible continuum emission reinforces the theory that a substantial component of the seismic emission seen is a result of sudden heating of the low photosphere associated with the observed visible continuum emission. Furthermore, the low-altitude magnetic loop structure inferred from potential--field extrapolations in the flaring region suggests that there is a significant inverse correlation between the seismicity of a flare and the height of the magnetic loops that conduct the particle beams from the corona.

We carried out an electromagnetic acoustic analysis of the solar flare of 14 August 2004 in active region AR10656 from the radio to the hard X-ray spectrum. The flare was a GOES soft X-ray class M7.4 and produced a detectable sun quake, confirming earlier inferences that relatively low-energy flares may be able to generate sun quakes. We introduce the hypothesis that the seismicity of the active region is closely related to the heights of coronal magnetic loops that conduct high-energy particles from the flare. In the case of relatively short magnetic loops, chromospheric evaporation populates the loop interior with ionized gas relatively rapidly, expediting the scattering of remaining trapped high-energy electrons into the magnetic loss cone and their rapid precipitation into the chromosphere. This increases both the intensity and suddenness of the chromospheric heating, satisfying the basic conditions for an acoustic emission that penetrates into the solar interior.

We report the discovery of one of the most powerful sunquakes detected to date, produced by an X1.2-class solar flare in active region 10720 on 2005 January 15. We used helioseismic holography to image the source of seismic waves emitted into the solar interior from the site of the flare. Acoustic egression power maps at 3 and 6 mHz with a 2 mHz bandpass reveal a compact acoustic source strongly correlated with impulsive hard X-ray and visible-continuum emission along the penumbral neutral line separating the two major opposing umbrae in the $\delta$-configuration sunspot that predominates AR10720. The acoustic emission signatures were directly aligned with both hard X-ray and visible continuum emission that emanated during the flare. The visible continuum emission is estimated at $2.0 \times 10^{23}$ J, approximately 500 times the seismic emission of $\sim 4 \times 10^{20}$ J. The flare of 2005 January 15 exhibits the same close spatial alignment between the sources of the seismic emission and impulsive visible continuum emission as previous flares, reinforcing the hypothesis that the acoustic emission may be driven by heating of the low photosphere. However, it is a major exception in that there was no signature to indicate the inclusion of protons in the particle beams thought to supply the energy radiated by the flare. The continued strong coincidence between the sources of seismic emission and impulsive visible continuum emission in the case of a proton-deficient white-light flare lends substantial support to the ``back -- warming'' hypothesis, that the low photosphere is significantly heated by intense Balmer and Paschen continuum-edge radiation from the overlying chromosphere in white-light flares.