Shock waves are simulated passing through the Alfv\'en-acoustic equipartition layer in a stratified isothermal magneto-atmosphere. The recent ray-theoretic calculations of N\'u\~nez (2019) predicted smoothing of the shock through this layer, causing both the fast and slow components to emerge as continuous waves. However, it is found that the partial mode conversion expected from linear theory for oblique incidence of the shock on the magnetic field is accompanied by a smoothing of the slow-shock only, whilst the fast-shock persists. Explanations are presented based on MHD mode conversion and shock theory.

Coupling between fast magneto-acoustic and Alfv\'en waves can be observe in fully ionized plasmas mediated by stratification and 3D geometrical effects. In Paper I, Cally & Khomenko (2015) have shown that in a weakly ionized plasma, such as the solar photosphere and chromosphere, the Hall current introduces a new coupling mechanism. The present study extends the results from Paper I to the case of warm plasma. We report on numerical experiments where mode transformation is studied using quasi-realistic stratification in thermodynamic parameters resembling the solar atmosphere. This redresses the limitation of the cold plasma approximation assumed in Paper I, in particular allowing the complete process of coupling between fast and slow magneto-acoustic modes and subsequent coupling of the fast mode to the Alfv\'en mode through the Hall current. Our results confirm the efficacy of the mechanism proposed in Paper I for the solar case. We observe that the efficiency of the transformation is a sensitive function of the angle between the wave propagation direction and the magnetic field, and of the wave frequency. The efficiency increases when the field direction and the wave direction are aligned for increasing wave frequencies. After scaling our results to typical solar values, the maximum amplitude of the transformed Alfv\'en waves, for a frequency of 1 Hz, corresponds to an energy flux (measured above the height of peak Hall coupling) of $\sim10^3$ $\rm W\,m^{-2}$, based on an amplitude of 500 $\rm m\,s^{-1}$ at $\beta=1$, which is sufficient to play a major role in both quiet and active region coronal heating.

Previous helioseismology of sunspots has been sensitive to both the structural and magnetic aspects of sunspot structure. We aim to develop a technique that is insensitive to the magnetic component so the two aspects can be more readily separated. We study waves reflected almost vertically from the underside of a sunspot. Time-distance helioseismology was used to measure travel times for the waves. Ray theory and a detailed sunspot model were used to calculate travel times for comparison. It is shown that these large distance waves are insensitive to the magnetic field in the sunspot. The largest travel time differences for any solar phenomena are observed. With sufficient modeling effort, these should lead to better understanding of sunspot structure.

We analyze data from the Helioseismic Magnetic Imager (HMI) and the Atmospheric Imaging Assembly (AIA) instruments on board the Solar Dynamics Observatory (SDO) to characterize the spatio-temporal acoustic power distribution in active regions as a function of the height in the solar atmosphere. For this, we use Doppler velocity and continuum intensity observed using the magnetically sensitive line at 6173 \AA\ as well as intensity at 1600 \AA\ and 1700 \AA. We focus on the power enhancements seen around AR 11330 as a function of wave frequency, magnetic field strength, field inclination and observation height. We find that acoustic halos occur above the acoustic cutoff frequency and extends up to 10 mHz in HMI Doppler and AIA 1700 \AA\ observations. Halos are also found to be strong functions of magnetic field and their inclination angle. We further calculate and examine the spatially averaged relative phases and cross-coherence spectra and find different wave characteristics at different heights.

Fast sausage modes in solar magnetic coronal loops are only fully contained in unrealistically short dense loops. Otherwise they are leaky, losing energy to their surrounds as outgoing waves. This causes any oscillation to decay exponentially in time. Simultaneous observations of both period and decay rate therefore reveal the eigenfrequency of the observed mode, and potentially insight into the tubes' nonuniform internal structure. In this article, a global spectral description of the oscillations is presented that results in an implicit matrix eigenvalue equation where the eigenvalues are associated predominantly with the diagonal terms of the matrix. The off-diagonal terms vanish identically if the tube is uniform. A linearized perturbation approach, applied with respect to a uniform reference model, is developed that makes the eigenvalues explicit. The implicit eigenvalue problem is easily solved numerically though, and it is shown that knowledge of the real and imaginary parts of the eigenfrequency is sufficient to determine the width and density contrast of a boundary layer over which the tubes' enhanced internal densities drop to ambient values. Linearized density kernels are developed that show sensitivity only to the extreme outside of the loops for radial fundamental modes, especially for small density enhancements, with no sensitivity to the core. Higher radial harmonics do show some internal sensitivity, but these will be more difficult to observe. Only kink modes are sensitive to the tube centres. {Variation in internal and external Alfv\'en speed along the loop is shown to have little effect on the fundamental dimensionless eigenfrequency, though the associated eigenfunction becomes more compact at the loop apex as stratification increases, or may even displace from the apex.

A simple model of a periodic ensemble of closely packed flux tubes sitting atop a vertically stratified layer reveals that an incident fast wave from below preferentially converts almost immediately to Alfv\'en waves in the flux tubes, with kink waves restricted to at most a very few Fourier modes. This suggests that observations of coronal kink modes in such structured systems may greatly underestimate the net wave energy flux being transported into and through the corona, much of which may reside in harder-to-observe Alfv\'en waves. The processes of mode conversion/resonant absorption and Alfv\'en phase mixing are implicated. It is suggested that the Sun's internal p-mode field -- the 5-minute oscillations -- may contribute substantially to the process by supplying incident fast waves in the chromosphere that scatter and mode-convert in the tube ensemble.

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.

We use the exact solutions for magnetoacoustic waves in a two dimensional isothermal atmosphere with uniform inclined magnetic field to calculate the wave reflection, transmission, and conversion of slow and fast waves incident from above ($z=\infty$). This is relevant to the question of whether waves excited by flares in the solar atmosphere can penetrate the Alfv\'en/acoustic equipartition layer (which we identify as the canopy) to reach the photosphere with sufficient energy to create sunquakes. It is found that slow waves above the acoustic cutoff frequency efficiently penetrate (transmit) as acoustic (fast) waves if directed at a small attack angle to the magnetic field, with the rest converting to magnetic (slow) waves, in accord with Generalized Ray Theory. This may help explain the compact nature of seismic sources of sunquakes identified using seismic holography. The incident slow waves can also efficiently transmit at low frequency in inclined field due to the reduction in acoustic cutoff frequency (ramp effect). Incident fast (magnetic) "waves" from infinity with specified nonzero horizontal wavenumber are necessarily evanescent, but can carry energy to the equipartition level by tunnelling. It is found that this can then efficiently convert to acoustic (fast) energy that can again reach the photosphere as a travelling wave. Overall, there appear to be ample avenues for substantial compressive wave energy to penetrate the canopy and impact the photosphere.

The photospheric temperature minimum in the Sun and solar-like stars is very weakly ionized, with ionization fraction $f$ as low as $10^{-4}$. In galactic star forming regions, $f$ can be $10^{-10}$ or lower. Under these circumstances, the Hall current can couple low frequency Alfv\'en and magneto\-acoustic waves via the dimensionless Hall parameter $\epsilon=\omega/\Omega_\text{i}f$, where $\omega$ is the wave frequency and $\Omega_\text{i}$ is the mean ion gyrofrequency. This is analysed in the context of a cold (zero-$\beta$) plasma, and in less detail for a warm plasma. It is found that Hall coupling preferentially occurs where the wave vector is nearly field-aligned. In these circumstances, Hall coupling in theory produces a continual oscillation between fast and Alfv\'en modes as the wave passes through the weakly ionized region. At low frequencies (mHz), characteristic of solar and stellar normal modes, $\epsilon$ is probably too small for more than a fraction of one oscillation to occur. On the other hand, the effect may be significant at the far higher frequencies (Hz) associated with magnetic reconnection events. In another context, characteristic parameters for star forming gas clouds suggest that $\mathcal{O}(1)$ or more full oscillations may occur in one cloud crossing. This mechanism is not expected to be effective in sunspots, due to their high ion gyrofrequencies and Alfv\'en speeds, since the net effect depends inversely on both and therefore inverse quadratically on field strength.

After more than half a century of community support related to the science of "solar activity'', IAU's Commission 10 was formally discontinued in 2015, to be succeeded by C.E2 with the same area of responsibility. On this occasion, we look back at the growth of the scientific disciplines involved around the world over almost a full century. Solar activity and fields of research looking into the related physics of the heliosphere continue to be vibrant and growing, with currently over 2,000 refereed publications appearing per year from over 4,000 unique authors, publishing in dozens of distinct journals and meeting in dozens of workshops and conferences each year. The size of the rapidly growing community and of the observational and computational data volumes, along with the multitude of connections into other branches of astrophysics, pose significant challenges; aspects of these challenges are beginning to be addressed through, among others, the development of new systems of literature reviews, machine-searchable archives for data and publications, and virtual observatories. As customary in these reports, we highlight some of the research topics that have seen particular interest over the most recent triennium, specifically active-region magnetic fields, coronal thermal structure, coronal seismology, flares and eruptions, and the variability of solar activity on long time scales. We close with a collection of developments, discoveries, and surprises that illustrate the range and dynamics of the discipline.

Our previous semi-analytic treatment of f- and p-mode multiple scattering from ensembles of thin flux tubes (Hanson and Cally, Astrophys. J. 781, 125; 791, 129, 2014) is extended by allowing both sausage and kink waves to freely escape at the top of the model using a radiative boundary condition there. As expected, this additional avenue of escape, supplementing downward loss into the deep solar interior, results in substantially greater absorption of incident f- and p-modes. However, less intuitively, it also yields mildly to substantially smaller phase shifts in waves emerging from the ensemble. This may have implications for the interpretation of seismic data for solar plage regions, and in particular their small measured phase shifts.

We present a technique to construct a spectropolarimetrically accurate magneto-hydrostatic model of a large-scale solar magnetic field concentration, mimicking a sunspot. Using the constructed model we perform a simulation of acoustic wave propagation, conversion and absorption in the solar interior and photosphere with the sunspot embedded into it. With the $6173\mathrm{\AA}$ magnetically sensitive photospheric absorption line of neutral iron, we calculate observable quantities such as continuum intensities, Doppler velocities, as well as full Stokes vector for the simulation at various positions at the solar disk, and analyse the influence of non-locality of radiative transport in the solar photosphere on helioseismic measurements. Bisector shapes were used to perform multi-height observations. The differences in acoustic power at different heights within the line formation region at different positions at the solar disk were simulated and characterised. An increase in acoustic power in the simulated observations of the sunspot umbra away from the solar disk centre was confirmed as the slow magneto-acoustic wave.

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.

Motivated by the observational results of Braun (1995), we extend the model of Hanson & Cally (2014) to address the effect of multiple scattering of f and p-modes by an ensemble of thin vertical magnetic flux tubes in the surface layers of the Sun. As in observational Hankel analysis we measure the scatter and phase shift from an incident cylindrical wave in a coordinate system roughly centred in the core of the ensemble. It is demonstrated that, although thin flux tubes are unable to interact with high order fluting modes individually, they can indirectly absorb energy from these waves through the scatters of kink and sausage components. It is also shown how the distribution of absorption and phase shift across the azimuthal order m depends strongly on the tube position, as well as on the individual tube characteristics. This is the first analytical study into an ensembles multiple scattering regime, that is embedded within a stratified atmosphere.

The Solar Orbiter mission, to be launched in July 2017, will carry a suite of remote sensing and in-situ instruments, including the Polarimetric and Helioseismic Imager (PHI). PHI will deliver high-cadence images of the Sun in intensity and Doppler velocity suitable for carrying out novel helioseismic studies. The orbit of the Solar Orbiter spacecraft will reach a solar latitude of up to 21 deg (up to 34 deg by the end of the extended mission) and thus will enable the first local helioseismology studies of the polar regions. Here we consider an array of science objectives to be addressed by helioseismology within the baseline telemetry allocation (51 Gbit per orbit, current baseline) and within the science observing windows (baseline 3 x 10 days per orbit). A particularly important objective is the measurement of large-scale flows at high latitudes (rotation and meridional flow), which are largely unknown but play an important role in flux transport dynamos. The full range of Earth-Sun-spacecraft angles provided by the orbit will enable helioseismology from two vantage points by combining PHI with another instrument: stereoscopic helioseismology will allow the study of the deep solar interior and a better understanding of the physics of solar oscillations in both quiet Sun and sunspots. We have used a model of the PHI instrument to study its performance for helioseismology applications. As input we used a 6 hr time-series of realistic solar magneto-convection simulation (Stagger code) and the SPINOR radiative transfer code to synthesize the observables. The simulated power spectra of solar oscillations show that the instrument is suitable for helioseismology. In particular, the specified point spread function, image jitter, and photon noise are no obstacle to a successful mission.

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.

We expand on recent studies to analytically model the behavior of two thin flux tubes interacting through the near- and acoustic far-field. The multiple scattering that occurs between the pair alters the absorption and phase of the outgoing wave, when compared to non-interacting tubes. We have included both the sausage and kink scatter produced by the pair. It is shown that the sausage mode's contribution to the scattered wave field is significant, and plays an equally important role in the multiple scattering regime. A disparity between recent numerical results and analytical studies, in particular the lack of symmetry between the two kink modes, is addressed. This symmetry break is found to be caused by an incorrect solution for the near-field modes.

Using advanced numerical magneto-hydrodynamic simulations of the magnetised solar photosphere, including non-grey radiative transport and a non-ideal equation of state, we analyse plasma motions in photospheric magnetic vortices. We demonstrate that apparent vortex-like motions in photospheric magnetic field concentrations do not exhibit "tornado"-like behaviour or a "bath-tub" effect. While at each time instance the velocity field lines in the upper layers of the solar photosphere show swirls, the test particles moving with the time-dependent velocity field do not demonstrate such structures. Instead, they move in a wave-like fashion with rapidly changing and oscillating velocity field, determined mainly by magnetic tension in the magnetised intergranular downflows. Using time-distance diagrams, we identify horizontal motions in the magnetic flux tubes as torsional Alfv\'en perturbations propagating along the nearly vertical magnetic field lines with local Alfv\'en speed.

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.

The role of the acoustic continuum associated with compact sources in the Sun's interior wave field is explored for a simple polytropic model. The continuum produces a near-field acoustic structure -- the so-called `acoustic jacket' -- that cannot be represented by a superposition of discrete normal modes. Particular attention is paid to monochromatic point sources of various frequency and depth, and to the surface velocity power that results, both in the discrete f- and p-mode spectrum and in the continuum. It is shown that a major effect of the continuum is to heal the surface wave field produced by compact sources, and therefore to hide them from view. It is found that the continuous spectrum is not a significant contributor to observable inter-ridge seismic power.

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.

Magneto-atmospheres with Alfv\'en speed [a] that increases monotonically with height are often used to model the solar atmosphere, at least out to several solar radii. A common example involves uniform vertical or inclined magnetic field in an isothermal atmosphere, for which the Alfv\'en speed is exponential. We address the issue of internal reflection in such atmospheres, both for time-harmonic and for transient waves. It is found that a mathematical boundary condition may be devised that corresponds to perfect absorption at infinity, and, using this, that many atmospheres where a(x) is analytic and unbounded present no internal reflection of harmonic Alfv\'en waves. However, except for certain special cases, such solutions are accompanied by a wake, which may be thought of as a kind of reflection. For the initial-value problem where a harmonic source is suddenly switched on (and optionally off), there is also an associated transient that normally decays with time as O(t-1) or O(t-1 ln t), depending on the phase of the driver. Unlike the steady-state harmonic solutions, the transient does reflect weakly. Alfv\'en waves in the solar corona driven by a finite-duration train of p-modes are expected to leave such transients.

Alfv\'en waves may be difficult to excite at the photosphere due to low ionization fraction and suffer near-total reflection at the transition region (TR). Yet they are ubiquitous in the corona and heliosphere. To overcome these difficulties, we show that they may instead be generated high in the chromosphere by conversion from reflecting fast magnetohydrodynamic waves, and that Alfv\'enic transition region reflection is greatly reduced if the fast reflection point is within a few scale heights of the TR. The influence of mode conversion on the phase of the reflected fast wave is also explored. This phase can potentially be misinterpreted as a travel speed perturbation, with implications for the practical seismic probing of active regions.

We study the conversion of fast magneto-acoustic waves to Alfven waves by means of 2.5D numerical simulations in a sunspot-like magnetic configuration. A fast, essentially acoustic, wave of a given frequency and wave number is generated below the surface and propagates upward though the Alfven/acoustic equipartition layer where it splits into upgoing slow (acoustic) and fast (magnetic) waves. The fast wave quickly reflects off the steep Alfven speed gradient, but around and above this reflection height it partially converts to Alfven waves, depending on the local relative inclinations of the background magnetic field and the wavevector. To measure the efficiency of this conversion to Alfven waves we calculate acoustic and magnetic energy fluxes. The particular amplitude and phase relations between the magnetic field and velocity oscillations help us to demonstrate that the waves produced are indeed Alfven waves. We find that the conversion to Alfven waves is particularly important for strongly inclined fields like those existing in sunspot penumbrae. Equally important is the magnetic field orientation with respect to the vertical plane of wave propagation, which we refer to as "field azimuth". For field azimuth less than 90 degrees the generated Alfven waves continue upwards, but above 90 degrees downgoing Alfven waves are preferentially produced. This yields negative Alfven energy flux for azimuths between 90 and 180 degrees. Alfven energy fluxes may be comparable to or exceed acoustic fluxes, depending upon geometry, though computational exigencies limit their magnitude in our simulations.

Motivated by observations that only a very few stars have been found to have antisolar differential rotation, much weaker in amplitude than that of the Sun, we analyze the stability of antisolar and solar type latitudinal differential rotations in the tachoclines of typical F,G and K stars. We employ two three-dimensional thin shell models, one for a Boussinesq but non-hydrostatic system, the other for a hydrostatic but non-Boussinesq system. We find that, in general, the combination of toroidal field band and differential rotation is more unstable, and unstable for lower toroidal fields, for antisolar than for solar-type differential rotation. In the antisolar case, the instability is always found to weaken the differential rotation, even if the primary energy source for the instability is the magnetic field. This favors surface antisolar differential rotations in stars being weaker than solar types, if the instability in the tachocline is felt at the surface of the star. This is most likely to happen in F stars, whose convection zones are much thinner than they are in G and K stars. This effect could help explain why the antisolar differential rotations that have been found are very weak compared to the Sun.

Modelling of adiabatic gravity wave propagation in the solar atmosphere showed that mode conversion to field guided acoustic waves or Alfv\'en waves was possible in the presence of highly inclined magnetic fields. This work aims to extend the previous adiabatic study, exploring the consequences of radiative damping on the propagation and mode conversion of gravity waves in the solar atmosphere. We model gravity waves in a VAL-C atmosphere, subject to a uniform, and arbitrarily orientated magnetic field, using the Newton cooling approximation for radiatively damped propagation. The results indicate that the mode conversion pathways identified in the adiabatic study are maintained in the presence of damping. The wave energy fluxes are highly sensitive to the form of the height dependence of the radiative damping time. While simulations starting from 0.2 Mm result in modest flux attenuation compared to the adiabatic results, short damping times expected in the low photosphere effectively suppress gravity waves in simulations starting at the base of the photosphere. It is difficult to reconcile our results and observations of propagating gravity waves with significant energy flux at photospheric heights unless they are generated in situ, and even then, why they are observed to be propagating as low as 70 km where gravity waves should be radiatively overdamped.

Alfv\'en waves may be generated via mode conversion from fast magneto-acoustic waves near their reflection level in the solar atmosphere, with implications both for coronal oscillations and for active region helio-seismology. In active regions this reflection typically occurs high enough that the Alfv\'en speed $a$ greatly exceeds the sound speed $c$, well above the $a=c$ level where the fast and slow modes interact. In order to focus on the fundamental characteristics of fast/Alfv\'en conversion, stripped of unnecessary detail, it is therefore useful to freeze out the slow mode by adopting the gravitationally stratified cold MHD model $c\to0$. This provides a benchmark for fast-to-Alfv\'en mode conversion in more complex atmospheres. Assuming a uniform inclined magnetic field and an exponential Alfv\'en speed profile with density scale height $h$, the Alfv\'en conversion coefficient depends on three variables only; the dimensionless transverse-to-the-stratification wavenumber $\kappa=kh$, the magnetic field inclination from the stratification direction $\theta$, and the polarization angle $\phi$ of the wavevector relative to the plane containing the stratification and magnetic field directions. We present an extensive exploration of mode conversion in this parameter space and conclude that near-total conversion to outward-propagating Alfv\'en waves typically occurs for small $\theta$ and large $\phi$ ($80^\circ$--$90^\circ$), though it is absent entirely when $\theta$ is exactly zero (vertical field). For wavenumbers of helioseismic interest, the conversion region is broad enough to encompass the whole chromosphere.

A fast-wave pulse in a simple, cold, inhomogeneous MHD model plasma is constructed by Fourier superposition over frequency of harmonic waves that are singular at their respective Alfven resonances. The pulse partially reflects before reaching the resonance layer, but also partially tunnels through to it to mode convert to an Alfven wave. The exact absorption/conversion coefficient for the pulse is shown to be given precisely by a function of transverse wavenumber tabulated in Paper I of this sequence, and to be independent of frequency and pulse width.

We study the coupling of magneto-acoustic waves to Alven waves using 2.5D numerical simulations. In our experiment, a fast magnetoacoustic wave of a given frequency and wavenumber is generated below the surface. The magnetic field in the domain is assumed homogeneous and inclined. The efficiency of the conversion to Alfven waves near the layer of equal acoustic and Alfven speeds is measured calculating their energy flux. The particular amplitude and phase relations between the oscillations of magnetic field and velocity help us to demonstrate that the waves produced after the transformation and reaching upper atmosphere are indeed Alfven waves. We find that the conversion from fast magneto-acoustic waves to Alfven waves is particularly important for the inclination and azimuth angles of the magnetic field between 55 and 65 degrees, with the maximum shifted to larger inclinations for lower frequency waves. The maximum Alfven flux transmitted to the upper atmosphere is about 2-3 times lower than the corresponding acoustic flux.

Resonant absorption and mode conversion are both extensively studied mechanisms for wave "absorption" in solar magnetohydrodynamics (MHD). But are they really distinct? We re-examine a well-known simple resonant absorption model in a cold MHD plasma that places the resonance inside an evanescent region. The normal mode solutions display the standard singular resonant features. However, these same normal modes may be used to construct a ray bundle which very clearly undergoes mode conversion to an Alfv\'en wave with no singularities. We therefore conclude that resonant absorption and mode conversion are in fact the same thing, at least for this model problem. The prime distinguishing characteristic that determines which of the two descriptions is most natural in a given circumstance is whether the converted wave can provide a net escape of energy from the conversion/absorption region of physical space. If it cannot, it is forced to run away in wavenumber space instead, thereby generating the arbitrarily small scales in situ that we recognize as fundamental to resonant absorption and phase mixing. On the other hand, if the converted wave takes net energy way, singularities do not develop, though phase mixing may still develop with distance as the wave recedes.

Using helioseismic holography strong evidence is presented that the phase (or equivalent travel-time) of helioseismic signatures in Dopplergrams within sunspots depend upon the line-of-sight angle in the plane containing the magnetic field and vertical directions. This is shown for the velocity signal in the penumbrae of two sunspots at 3, 4 and 5 mHz. Phase-sensitive holography demonstrates that they are significantly affected in a strong, moderately inclined magnetic field. This research indicates that the effects of the surface magnetic field are potentially very significant for local helioseismic analysis of active regions.

We numerically study the evolution of magnetic fields and fluid flows in a thin spherical shell. We take the initial field to be a latitudinally confined, predominantly toroidal flux tube. For purely toroidal, untwisted flux tubes, we recover previously known radial-shredding instabilities, and show further that in the nonlinear regime these instabilities can very effectively destroy the original field. For twisted flux tubes, including also a poloidal component, there are several possibilities, including the suppression of the radial-shredding instability, but also a more directly induced evolution, brought about because twisted flux tubes in general are not equilibrium solutions of the governing equations.

The detection of upward propagating internal gravity waves at the base of the Sun's chromosphere has recently been reported by Straus et al., who postulated that these may efficiently couple to Alfven waves in magnetic regions. This may be important in transporting energy to higher levels. Here we explore the propagation, reflection and mode conversion of linear gravity waves in a VAL C atmosphere, and find that even weak magnetic fields usually reflect gravity waves back downward as slow magnetoacoustic waves well before they reach the Alfven/acoustic equipartition height at which mode conversion might occur. However, for certain highly inclined magnetic field orientations in which the gravity waves manage to penetrate near or through the equipartition level, there can be substantial conversion to either or both upgoing Alfven and acoustic waves. Wave energy fluxes comparable to the chromospheric radiative losses are expected.

Phase perturbations due to inclined surface magnetic field of active region strength are calculated numerically in quiet Sun and simple sunspot models in order to estimate and compare the direct and indirect (thermal) effects of the fields on helioseismic waves. It is found that the largest direct effects occur in highly inclined field characteristic of penumbrae, and scale roughly linearly with magnetic field strength. The combined effects of sunspot magnetic and thermal anomalies typically yield negative travel-time perturbations in penumbrae. Travel-time shifts in umbrae depend on details of how the thermal and density structure differs from the quiet Sun. The combined shifts are generally not well approximated by the sum of the thermal and magnetic effects applied separately, except at low field strengths of around 1 kG or less, or if the thermal shift is small. A useful rule-of-thumb appears to be that travel-time perturbations in umbrae are predominantly thermal, whereas in penumbrae they are mostly magnetic.

Generalized Ray Theory (GRT) provides a simple description of MHD mode transmission and conversion between magnetoacoustic fast and slow waves and is directly applicable to solar active regions. Here it is tested in a simple two-dimensional, isothermal, gravitationally-stratified model with inclined magnetic field using previously published exact solutions and found to perform very well.

We study the near-field coupling of a pair of flux tubes embedded in a gravitationally stratified environment. The mutual induction of the near-field {\it jackets} of the two flux tubes can considerably alter the scattering properties of the system, resulting in sizable changes in the magnitudes of scattering coefficients and bizarre trends in the phases. The dominant length scale governing the induction zone turns out to be approximately half the horizontal wave length of the incident mode, a result that fits in quite pleasantly with extant theories of scattering. Higher-$\beta$ flux tubes are more strongly coupled than weaker ones, a consequence of the greater role that the near-field jacket modes play in the such tubes. We also comment on the importance of incorporating the effects of multiple scattering when studying the effects of mode absorption in plage and interpreting related scattering measurements. That the near-field plays such an important role in the scattering process lends encouragement to the eventual goal of observationally resolving sub-wavelength features of flux tubes using techniques of helioseismology.

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?

At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic waves and the local surface motion within two sunspots. These correlations are modeled assuming the surface motion is elliptical, and we explore the properties of the elliptical motion on the magnetic field inclination. We also demonstrate that the phase shift of the outward propagating waves is opposite to the phase shift of the inward propagating waves in stronger, more vertical fields, but similar to the inward phase shifts in weaker, more inclined fields.

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.

The efficacy of fast/slow MHD mode conversion in the surface layers of sunspots has been demonstrated over recent years using a number of modelling techniques, including ray theory, perturbation theory, differential eigensystem analysis, and direct numerical simulation. These show that significant energy may be transferred between the fast and slow modes in the neighbourhood of the equipartition layer where the Alfven and sound speeds coincide. However, most of the models so far have been two dimensional. In three dimensions the Alfven wave may couple to the magneto-acoustic waves with important implications for energy loss from helioseismic modes and for oscillations in the atmosphere above the spot. In this paper, we carry out a numerical ``scattering experiment'', placing an acoustic driver 4 Mm below the solar surface and monitoring the acoustic and Alfvenic wave energy flux high in an isothermal atmosphere placed above it. These calculations indeed show that energy conversion to upward travelling Alfven waves can be substantial, in many cases exceeding loss to slow (acoustic) waves. Typically, at penumbral magnetic field strengths, the strongest Alfven fluxes are produced when the field is inclined 30-40 degrees from the vertical, with the vertical plane of wave propagation offset from the vertical plane containing field lines by some 60-80 degrees.

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.

Aims: The primary objective of this study is to search for and identify wave modes within a sunspot penumbra. Methods: Infrared spectropolarimetric time series data are inverted using a model comprising two atmospheric components in each spatial pixel. Fourier phase difference analysis is performed on the line-of-sight velocities retrieved from both components to determine time delays between the velocity signals. In addition, the vertical separation between the signals in the two components is calculated from the Stokes velocity response functions. Results: The inversion yields two atmospheric components, one permeated by a nearly horizontal magnetic field, the other with a less-inclined magnetic field. Time delays between the oscillations in the two components in the frequency range 2.5-4.5 mHz are combined with speeds of atmospheric wave modes to determine wave travel distances. These are compared to expected path lengths obtained from response functions of the observed spectral lines in the different atmospheric components. Fast-mode (i.e., modified p-mode) waves exhibit the best agreement with the observations when propagating toward the sunspot at an angle ~50 degrees to the vertical.

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.