Context: Solar sectoral Rossby waves have only recently been unambiguously identified in Helioseimsic and Magnetic Imager (HMI) and Michelson Doppler Imager (MDI) maps of flows near the solar surface. So far this has not been done with the Global Oscillation Network Group (GONG) ground-based observations, which have different noise properties. Aims: We utilize 17 years of GONG++ data, to identify and characterize solar Rossby waves using ring-diagram helioseismology. We compare directly with HMI ring-diagram analysis. Methods: Maps of the radial vorticity are obtained for flows within the top 2 Mm of the surface for 17 years of GONG++. The data is corrected for systematic effects including the annual periodicity related to the $B_0$ angle. We then compute the Fourier components of the radial vorticity of the flows in the co-rotating frame. We perform the same analysis on the HMI data that overlap in time. Results: We find that the solar Rossby waves have measurable amplitudes in the GONG++ sectoral power spectra for azimuthal orders between $m=3$ and $m=15$. The measured mode characteristics (frequencies, lifetimes and amplitudes) from GONG++ are consistent with the HMI measurements in the overlap period from 2010 to 2018 for $m\le9$. For higher-$m$ modes the amplitudes and frequencies agree within two sigmas. The signal-to-noise ratio of modes in GONG++ power spectra is comparable to HMI for $8\le m\le11$, but is lower by a factor of two for other modes. Conclusions: The GONG++ data provide a long and uniform data set to study solar global-scale Rossby waves from 2001.

The inversion of ring fit parameters to obtain subsurface flow maps in ring-diagram analysis for 8 years of SDO observations is computationally expensive, requiring ~3200 CPU hours. In this paper we apply machine learning techniques to the inversion in order to speed up calculations. Specifically, we train a predictor for subsurface flows using the mode fit parameters and the previous inversion results, to replace future inversion requirements. We utilize Artificial Neural Networks as a supervised learning method for predicting the flows in 15 degree ring tiles. To demonstrate that the machine learning results still contain the subtle signatures key to local helioseismic studies, we use the machine learning results to study the recently discovered solar equatorial Rossby waves. The Artificial Neural Network is computationally efficient, able to make future flow predictions of an entire Carrington rotation in a matter of seconds, which is much faster than the current ~31 CPU hours. Initial training of the networks requires ~3 CPU hours. The trained Artificial Neural Network can achieve a root mean-square error equal to approximately half that reported for the velocity inversions, demonstrating the accuracy of the machine learning (and perhaps the overestimation of the original errors from the ring-diagram pipeline). We find the signature of equatorial Rossby waves in the machine learning flows covering six years of data, demonstrating that small-amplitude signals are maintained. The recovery of Rossby waves in the machine learning flow maps can be achieved with only one Carrington rotation (27.275 days) of training data. We have shown that machine learning can be applied to, and perform more efficiently than the current ring-diagram inversion. The computation burden of the machine learning includes 3 CPU hours for initial training, then around 0.0001 CPU hours for future predictions.

The interpretation of helioseismic measurements, such as wave travel-time, is based on the computation of kernels that give the sensitivity of the measurements to localized changes in the solar interior. These are computed using the ray or the Born approximation. The Born approximation is preferable as it takes finite-wavelength effects into account, but can be computationally expensive. We propose a fast algorithm to compute travel-time sensitivity kernels under the assumption that the background solar medium is spherically symmetric. Kernels are typically expressed as products of Green's functions that depend upon depth, latitude and longitude. Here, we compute the spherical harmonic decomposition of the kernels and show that the integrals in latitude and longitude can be performed analytically. In particular, the integrals of the product of three associated Legendre polynomials can be computed thanks to the algorithm of Dong and Lemus (2002). The computations are fast and accurate and only require the knowledge of the Green's function where the source is at the pole. The computation time is reduced by two orders of magnitude compared to other recent computational frameworks. This new method allows for flexible and computationally efficient calculations of a large number of kernels, required in addressing key helioseismic problems. For example, the computation of all the kernels required for meridional flow inversion takes less than two hours on 100 cores.

The temporal covariance between seismic waves measured at two locations on the solar surface is the fundamental observable in time-distance helioseismology. Above the acoustic cut-off frequency ($\sim$5.3~mHz), waves are not trapped in the solar interior and the covariance function can be used to probe the upper atmosphere. We wish to implement appropriate radiative boundary conditions for computing the propagation of high-frequency waves in the solar atmosphere. We consider the radiative boundary conditions recently developed by Barucq et al. (2017) for atmospheres in which sound-speed is constant and density decreases exponentially with radius. We compute the cross-covariance function using a finite element method in spherical geometry and in the frequency domain. The ratio between first- and second-skip amplitudes in the time-distance diagram is used as a diagnostic to compare boundary conditions and to compare with observations. We find that a boundary condition applied 500 km above the photosphere and derived under the approximation of small angles of incidence accurately reproduces the `infinite atmosphere' solution for high-frequency waves. When the radiative boundary condition is applied 2 Mm above the photosphere, we find that the choice of atmospheric model affects the time-distance diagram. In particular, the time-distance diagram exhibits double-ridge structure when using a VAL atmospheric model.

Local helioseismology has so far relied on semi-analytical methods to compute the spatial sensitivity of wave travel times to perturbations in the solar interior. These methods are cumbersome and lack flexibility. Here we propose a convenient framework for numerically solving the forward problem of time-distance helioseismology in the frequency domain. The fundamental quantity to be computed is the cross-covariance of the seismic wavefield. We choose sources of wave excitation that enable us to relate the cross-covariance of the oscillations to the Green's function in a straightforward manner. We illustrate the method by considering the 3D acoustic wave equation in an axisymmetric reference solar model, ignoring the effects of gravity on the waves. The symmetry of the background model around the rotation axis implies that the Green's function can be written as a sum of longitudinal Fourier modes, leading to a set of independent 2D problems. We use a high-order finite-element method to solve the 2D wave equation in frequency space. The computation is `embarrassingly parallel', with each frequency and each azimuthal order solved independently on a computer cluster. We compute travel-time sensitivity kernels in spherical geometry for flows, sound speed, and density perturbations under the first Born approximation. Convergence tests show that travel times can be computed with a numerical precision better than one millisecond, as required by the most precise travel-time measurements. The method presented here is computationally efficient and will be used to interpret travel-time measurements in order to infer, e.g., the large-scale meridional flow in the solar convection zone. It allows the implementation of (full-waveform) iterative inversions, whereby the axisymmetric background model is updated at each iteration.

The use of acoustic holography in the high-frequency $p$-mode spectrum can resolve the source distributions of enhanced acoustic emissions within halo structures surrounding active regions. In doing so, statistical methods can then be applied to ascertain relationships with the magnetic field. This is the focus of this study. The mechanism responsible for the detected enhancement of acoustic sources around solar active regions has not yet been explained. Furthermore the relationship between the magnetic field and enhanced acoustic emission has not yet been comprehensively examined. We have used vector magnetograms from the \Helioseismic and Magnetic Imager (HMI) on-board the Solar Dynamics Observatory (SDO) to image the magnetic-field properties in the halo. We have studied the acoustic morphology of an active region, with a complex halo and "glories," and we have linked some acoustic properties to the magnetic-field configuration. In particular, we find that acoustic sources are significantly enhanced in regions of intermediate field strength with inclinations no different from the distributions found in the quiet Sun. Additionally we have identified a transition region between the active region and the halo, in which the acoustic source power is hindered by inclined fields of intermediate field strength. Finally, we have compared the results of acoustic emission maps, calculated from holography, and the commonly used local acoustic maps, finding that the two types of maps have similar properties with respect to the magnetic field but lack spatial correlation when examining the highest-powered regions.

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.

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.

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.

Multi-frequency power maps of the local acoustic oscillations show acoustic enhancements (\acoustic-power halos") at high frequencies surrounding large active region. Computational seismic holography reveals a high-frequency \acoustic-emission halo", or \seismic glory" surrounding large active regions. In this study, we have applied computational seismic holography to map the seismic seismic source density surrounding AR 11429. Studies of HMI/SDO Doppler data, shows that the "acoustic halos" and the "seismic glories" are prominent at high frequencies 5-8 mHz. We investigate morphological properties of acoustic-power and acoustic emission halos around an active region to see if they are spatially correlated. Details about the local magnetic field from vectormagnetograms of AR 11429 are included. We identify a 15" region of seismic deficit power (dark moat) shielding the white-light boundary of the active region. The size of the seismic moat is related to region of intermediate magnetic field strength. The acoustic moat is circled by the halo of enhanced seismic amplitude as well as enhanced seismic emission. Overall, the results suggest that features are related. However, if we narrow the frequency band to 5.5-6.5 mHz, we find that the seismic source density dominates over the local acoustic power, suggesting the existence of sources that emit more energy downward into the solar interior than upward toward the solar surface.

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.