Emanating from the base of the Sun's corona, the solar wind fills the interplanetary medium with a magnetized stream of charged particles whose interaction with the Earth's magnetosphere has space-weather consequences such as geomagnetic storms. Accurately predicting the solar wind through measurements of the spatio-temporally evolving conditions in the solar atmosphere is important but remains an unsolved problem in heliophysics and space-weather research. In this work, we use deep learning for prediction of solar wind (SW) properties. We use Extreme Ultraviolet images of the solar corona from space based observations to predict the SW speed from the NASA OMNIWEB dataset, measured at Lagragian point 1. We evaluate our model against autoregressive and naive models, and find that our model outperforms the benchmark models, obtaining a best-fit correlation of 0.55 $\pm$ 0.03 with the observed data. Upon visualization and investigation of how the model uses data to make predictions, we find higher activation at the coronal holes for fast wind prediction ($\approx$ 3 to 4 days prior to prediction), and at the active regions for slow wind prediction. These trends bear an uncanny similarity to the influence of regions potentially being the sources of fast and slow wind, as reported in literature. This suggests that our model was able to learn some of the salient associations between coronal and solar wind structure without built-in physics knowledge. Such an approach may help us discover hitherto unknown relationships in heliophysics datasets.

Solar flares are explosions in the solar atmosphere that release intense bursts of short-wavelength radiation and are capable of producing severe space-weather consequences. Flares release free energy built up in coronal fields, which are rooted in active regions (ARs) on the photosphere, via magnetic reconnection. The exact processes that lead to reconnection are not fully known and therefore reliable forecasting of flares is challenging. Recently, photospheric magnetic-field data has been extensively analysed using machine learning (ML) and these studies suggest that flare-forecasting accuracy does not strongly depend on how long in advance flares are predicted (Bobra & Couvidat 2015; Raboonik et al. 2017; Huang et al. 2018). Here, we use ML to understand the evolution of AR magnetic fields before and after flares. We explicitly train convolutional neural networks (CNNs) to classify SDO/HMI line-of-sight magnetograms into ARs producing at least one M- or X-class flare or as nonflaring. We find that flaring ARs remain in flare-productive states -- marked by recall >60% with a peak of ~ 80% -- days before and after flares. We use occlusion maps and statistical analysis to show that the CNN pays attention to regions between the opposite polarities from ARs and the CNN output is dominantly decided by the total unsigned line-of-sight flux of ARs. Using synthetic bipole magnetograms, we find spurious dependencies of the CNN output on magnetogram dimensions for a given bipole size. Our results suggest that it is important to use CNN designs that eliminate such artifacts in CNN applications for processing magnetograms and, in general, solar image data.

Departures from standard spherically symmetric solar models, in the form of perturbations such as global and local-scale flows and structural asphericities, result in the splitting of eigenfrequencies in the observed spectrum of solar oscillations. Drawing from prevalent ideas in normal-mode coupling theory in geophysical literature, we devise a procedure that enables the computation of sensitivity kernels for general Lorentz stress fields in the Sun. Mode coupling due to any perturbation requires careful consideration of self- and cross-coupling of multiplets. Invoking the isolated-multiplet approximation allows for limiting the treatment to purely self-coupling, requiring significantly less computational resources. We identify the presence of such isolated multiplets under the effect of Lorentz stresses in the Sun. Currently, solar missions allow precise measurements of self-coupling of multiplets via "$a$-coefficients" and the cross-spectral correlation signal which enables the estimation of the "structure coefficients". We demonstrate the forward problem for both self-coupling ($a$-coefficients) and cross-coupling (structure coefficients). In doing so, we plot the self-coupling kernels and estimate $a$-coefficients arising from a combination of deep-toroidal and surface-dipolar axisymmetric fields. We also compute the structure coefficients for an arbitrary general magnetic field (real and solenoidal) and plot the corresponding "splitting function", a convenient way to visualize the splitting of multiplets under 3D internal perturbations. The results discussed in this paper pave the way to formally pose an inverse problem, and infer solar internal magnetic fields.

Rossby waves play an important role in mediating the angular momentum of rotating spherical fluids, creating weather on Earth and tuning exoplanet orbits in distant stellar systems (Ogilvie 2014). Their recent discovery in the solar convection zone provides an exciting opportunity to appreciate the detailed astrophysics of Rossby waves. Large-scale Rossby waves create subtle drifts in acoustic oscillations in the convection zone, which we measure using helioseismology to image properties of Rossby waves in the interior. We analyze 20 years of space-based observations, from 1999 to 2018, to measure Rossby-mode frequencies, line-widths, and amplitudes. Spatial leakage affects the measurements of normal model coupling and complicates the analysis of separating out specific harmonic degree and the azimuthal number of features on the Sun. Here we demonstrate a novel approach to overcome this difficulty and test it by performing synthetic tests. We find that the root-mean-square velocity of the modes is of the order of 0.5 m/s at the surface.

Supergranules are divergent 30-Mm sized cellular flows observed everywhere at the solar photosphere. Their place in the hierarchy of convective structures and their origin remain poorly understood (Rincon et al., 2018). Estimating supergranular depth is of particular interest since this may help point to the underlying physics. However, their subsurface velocity profiles have proven difficult to ascertain. Birch et al. (2006) had suggested that helioseismic inferences would benefit from an ensemble average over multiple realizations of supergranules due to the reduction in realization noise. Bhattacharya et al. (2017) used synthetic forward-modelled seismic wave travel times and demonstrated the potential of helioseismic inversions at recovering the flow profile of an average supergranule that is separable in the horizontal and vertical directions, although the premise of this calculation has since been challenged by Ferret (2019). In this work we avoid this assumption and carry out a validation test of helioseismic travel-time inversions starting from plausible synthetic non-separable profiles of an average supergranule. We compute seismic wave travel times and sensitivity kernels by simulating wave propagation through this background. We find that, while the ability to recover the exact profile degrades based on the number of parameters involved, we are nevertheless able to recover the peak depth of our models in a few iterations where the measurements are presumably above the noise cutoff. This represents an important step towards unraveling the physics behind supergranules, as we start appreciating the parameters that we may reliably infer from a time-distance helioseismic inversion.

Solar flares - bursts of high-energy radiation responsible for severe space-weather effects - are a consequence of the occasional destabilization of magnetic fields rooted in active regions (ARs). The complexity of AR evolution is a barrier to a comprehensive understanding of flaring processes and accurate prediction. Though machine learning (ML) has been used to improve flare predictions, the potential for revealing precursors and associated physics has been underexploited. Here, we train ML algorithms to classify between vector-magnetic-field observations from flaring ARs, producing at least one M-/X-class flare, and non-flaring ARs. Analysis of magnetic-field observations accurately classified by the machine presents statistical evidence for (1) ARs persisting in flare-productive states --- characterized by AR area --- for days, before and after M- and X-class flare events, (2) systematic pre-flare build-up of free energy in the form of electric currents, suggesting that associated subsurface magnetic field is twisted, (3) intensification of Maxwell stresses in the corona above newly emerging ARs, days before first flares. These results provide new insights into flare physics and improving flare forecasting.

Asteroseismology has undergone a profound transformation as a scientific field following the CoRoT and Kepler space missions. The latter is now yielding the first measurements of latitudinal differential rotation obtained directly from oscillation frequencies. Differential rotation is a fundamental mechanism of the stellar dynamo effect. Our goal is to measure the amount of differential rotation in the solar analogues 16 Cyg A and B, which are the components of a binary system. These stars are the brightest observed by Kepler and have therefore been extensively observed, with exquisite precision on their oscillation frequencies. We modelled the acoustic power spectrum of 16 Cyg A and B using a model that takes into account the contribution of differential rotation to the rotational frequency splitting. The estimation was carried out in a Bayesian setting. We then inverted these results to obtain the rotation profile of both stars under the assumption of a solar-like functional form. We observe that the magnitude of latitudinal differential rotation has a strong chance of being solar-like for both stars, their rotation rates being higher at the equator than at the pole. The measured latitudinal differential rotation, defined as the difference of rotation rate between the equator and the pole, is $320\pm269$ nHz and $440^{+363}_{-383}$ nHz for 16 Cyg A and B, respectively, confirming that the rotation rates of these stars are almost solar-like. Their equatorial rotation rates are $535\pm75$ nHz and $565_{-129}^{+150}$ nHz. Our results are in good agreement with measurements obtained from spectropolarimetry, spectroscopy, and photometry. We present the first conclusive measurement of latitudinal differential rotation for solar analogues. Their rotational profiles are very close to those of the Sun. These results depend weakly on the uncertainties of the stellar parameters.

Rossby waves play a fundamental role in angular momentum processes in rotating fluids. In addition to the potential to shed light on physical mechanisms operating in the solar convection zone, the recent detection of Rossby waves in the Sun (L\"optien et al. 2018; Liang et al. 2018) also serves as a means of comparison between different helioseismic methods. Time-distance helioseismology, ring-diagram analysis and other techniques have all proven successful in recovering the Rossby-wave dispersion relation from analyses of Helioseismic and Magnetic Imager data (HMI; Schou et al. 2012). In this article, we demonstrate that analyses of two years of HMI global-mode-oscillation data using the technique of normal-mode coupling also show signatures of Rossby waves. In addition to providing an independent means of inferring Rossby waves, this detection lends credence to the methodology of mode coupling and encourages a more complete exploration of its possibilities.

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 differentially rotating outer layers of stars are thought to play a role in driving their magnetic activity, but the underlying mechanisms that generate and sustain differential rotation are poorly understood. We report the measurement of latitudinal differential rotation in the convection zones of 40 Sun-like stars using asteroseismology. For the most significant detections, the stars' equators rotate approximately twice as fast as their mid-latitudes. The latitudinal shear inferred from asteroseismology is much larger than predictions from numerical simulations.

Accurate inference of solar meridional flow is of crucial importance for the understanding of solar dynamo process. Wave travel times, as measured on the surface, will change if the waves encounter perturbations e.g. in the sound speed or flows, as they propagate through the solar interior. Using functions called sensitivity kernels, we may image the underlying anomalies that cause measured shifts in travel times. The inference of large-scale structures e.g meridional circulation requires computing sensitivity kernels in spherical geometry. Mandal et al. (2017) have computed such spherical kernels in the limit of the first-Born approximation. In this work, we perform an inversion for meridional circulation using travel-time measurements obtained from 6 years of SDO/HMI data and those sensitivity kernels. We enforce mass conservation by inverting for a stream function. The number of free parameters is reduced by projecting the solution on to cubic B-splines in radius and derivatives of the Legendre-polynomial basis in latitude, thereby improving the condition number of the inverse problem. We validate our approach for synthetic observations before performing the actual inversion. The inversion suggests a single-cell profile with the return-flow occurring at depths below 0.78 $R_\odot$.

Normal modes are coupled by the presence of perturbations in the Sun, providing a novel and under-appreciated helioseismic technique with which to image the solar interior. The process of measuring coupling between normal modes is straightforward, much more so when compared with other prevalent helioseismic techniques. The theoretical framework to interpret these measurements is well developed with the caveat that it applies only in the case where the entire surface of the Sun is observed. In practice however, the limited visibility of the Sun and line-of-sight related effects diminish the resolution of the technique. Here, we compute realistic sensitivities of normal-mode coupling measurements to flows in the solar interior and describe how to mitigate the sometimes-overwhelming effect of leakage. The importance of being able to isolate individual spherical harmonics and observe the full Sun, to which future solar observatories may aspire, is thus highlighted in our results. In the latter part of the article, we describe the noise model for the variance of coupling coefficients, a critical component to the process of inference.

Oscillation properties are usually measured by fitting symmetric Lorentzian profiles to the power spectra of Sun-like stars. However the line profiles of solar oscillations have been observed to be asymmetrical for the Sun. The physical origin of this line asymmetry is not fully understood, although it should depend on the depth dependence of the source of wave excitation (convective turbulence) and details of the observable (velocity or intensity). For oscillations of the Sun, it has been shown that neglecting the asymmetry leads to systematic errors in the frequency determination. This could subsequently affects the results of seismic inferences of the solar internal structure. Using light curves from the {\it Kepler} spacecraft we have measured mode asymmetries in 43 stars. We confirm that neglecting the asymmetry leads to systematic errors that can exceed the $1\sigma$ confidence intervals for seismic observations longer than one year. Therefore, the application of an asymmetric Lorentzian profile is to be favoured to improve the accuracy of the internal stellar structure and stellar fundamental parameters. We also show that the asymmetry changes sign between cool Sun-like stars and hotter stars. This provides the best constraints to date on the location of the excitation sources across the Hertzsprung-Russel diagram.

We develop a helioseismic inversion algorithm that can be used to recover sub-surface vertical profiles of 2-dimensional supergranular flows from surface measurements of synthetic wave travel times. We carry out seismic wave-propagation simulations through a 2-dimensional section of a flow profile that resembles an averaged supergranule, and a starting model that has flows only at the surface. We assume that the wave measurements are entirely without realization noise for the purpose of our test. We expand the vertical profile of the supergranule stream function on a basis of B-splines. We iteratively update the B-spline coefficients of the supergranule model to reduce the travel-times differences observed between the two simulations. We carry out the exercise for four different vertical profiles peaking at different depths below the solar surface. We are able to accurately recover depth profiles of four supergranule models at depths up to $8-10\,\text{Mm}$ below the solar surface using $f-p_4$ modes, under the assumption that there is no realization noise. We are able to obtain the peak depth and the depth of the return flow for each model. A basis-resolved inversion performs significantly better than one where the flow field is inverted for at each point in the radial grid. This is an encouraging result and might act as a guide in developing more realistic inversion strategies that can be applied to supergranular flows in the Sun.

Understanding the governing mechanism of solar magnetism remains an outstanding challenge in astrophysics. Seismology is the most compelling technique with which to infer the internal properties of the Sun and stars. Waves in the Sun, nominally acoustic, are sensitive to the emergence and cyclical strengthening of magnetic field, evidenced by measured changes in resonant oscillation frequencies that are correlated with the solar cycle. The inference of internal Lorentz stresses from these measurements has the potential to significantly advance our appreciation of the dynamo. Indeed, seismological inverse theory for the Sun is well understood for perturbations in composition, thermal structure and flows but, is not fully developed for magnetism, owing to the complexity of the ideal magnetohydrodynamic (MHD) equation. Invoking first-Born perturbation theory to characterize departures from spherically symmetric hydrostatic models of the Sun and applying the notation of generalized spherical harmonics, we calculate sensitivity functions of seismic measurements to the general time-varying Lorentz stress tensor. We find that eigenstates of isotropic (i.e. acoustic only) background models are dominantly sensitive to isotropic deviations in the stress tensor and much more weakly so to anisotropic stresses (and therefore challenging to infer). The apple cannot fall far from the tree.

In this article, we derive and compute the sensitivity of measurements of coupling between normal modes of oscillation in the Sun to underlying flows. The theory is based on first-Born perturbation theory, and the analysis is carried out using the formalism described by \citet{lavely92}. Albeit tedious, we detail the derivation and compute the sensitivity of specific pairs of coupled normal modes to anomalies in the interior. Indeed, these kernels are critical for the accurate inference of convective flow amplitudes and large-scale circulations in the solar interior. We resolve some inconsistencies in the derivation of \citet{lavely92} and reformulate the fluid-continuity condition. We also derive and compute sound-speed kernels, paving the way for inverting for thermal anomalies alongside flows.

The inference of internal properties of the Sun from surface measurements of wave travel times is the goal of time-distance helioseismology. A critical step toward the accurate interpretation of travel-time shifts is the computation of sensitivity functions linking seismic measurements to internal structure. Here we calculate finite-frequency sensitivity kernels in spherical geometry for two-point travel-time measurements. We numerically build Green's function by solving for it at each frequency and spherical-harmonic degree and summing over all these pieces. These computations are performed in parallel ("embarrassingly"), thereby achieving significant speedup in wall-clock time. Kernels are calculated by invoking the first-order Born approximation connecting deviations in the wavefield to perturbations in the operator. Validated flow kernels are shown to produce travel times within 0.47% of the true value for uniform flows up to 750 m/s. We find that travel-time can be obtained with errors of 1 millisecond or less for flows having magnitudes similar to meridional circulation. Alongside flows, we also compute and validate sensitivity kernel for sound-speed perturbations. These accurate sensitivity kernels might improve the current inferences of sub-surface flows significantly.

The advent of space-based observatories such as CoRoT and Kepler has enabled the testing of our understanding of stellar evolution on thousands of stars. Evolutionary models typically require five input parameters, the mass, initial Helium abundance, initial metallicity, mixing- length (assumed to be constant over time), and the age to which the star must be evolved. Some of these parameters are also very useful in characterizing the associated planets and in studying galactic archaeology. How to obtain these parameters from observations rapidly and accurately, specifically in the context of surveys of thousands of stars, is an outstanding ques- tion, one that has eluded straightforward resolution. For a given star, we typically measure the effective temperature and surface metallicity spectroscopically and low-degree oscillation frequencies through space observatories. Here we demonstrate that statistical learning, using artificial neural networks, is successful in determining the evolutionary parameters based on spectroscopic and seismic measurements. Our trained networks show robustness over a broad range of parameter space, and critically, are entirely computationally inexpensive and fully automated. We analyze the observations of a few stars using this method and the results com- pare well to inferences obtained using other techniques. This method is both computationally cheap and inferentially accurate, paving the way for analyzing the vast quantities of stellar observations from past, current, and future missions.

Observations of the solar surface reveal the presence of flows with length scales of around $35$ Mm, commonly referred to as supergranules. Inferring the sub-surface flow profile of supergranules from measurements of the surface and photospheric wavefield is an important challenge faced by helioseismology. Traditionally, the inverse problem has been approached by studying the linear response of seismic waves in a horizontally translationally invariant background to the presence of the supergranule; following an iterative approach that does not depend on horizontal translational invariance might perform better, since the misfit can be analyzed post iterations. In this work, we construct synthetic observations using a reference supergranule, and invert for the flow profile using surface measurements of travel-times of waves belonging to modal ridges $f$ (surface-gravity) and $p_{1}$ through $p_{7}$ (acoustic). We study the extent to which individual modes and their combinations contribute to infer the flow. We show that this method of non-linear iterative inversion tends to underestimate the flow velocities as well as inferring a shallower flow profile, with significant deviations from the reference supergranule near the surface. We carry out a similar analysis for a sound-speed perturbation and find that analogous near-surface deviations persist, although the iterations converge faster and more accurately. We conclude that a better approach to inversion would be to expand the supergranule profile in an appropriate basis, thereby reducing the number of parameters being inverted for and appropriately regularizing them.

Measurements of oscillation frequencies of the Sun and stars can provide important independent constraints on their internal structure and dynamics. Seismic models of these oscillations are used to connect structure and rotation of the star to its resonant frequencies, which are then compared with observations, the goal being that of minimizing the difference between the two. Even in the case of the Sun, for which structure models are highly tuned, observed frequencies show systematic deviations from modeled frequencies, a phenomenon referred to as the "surface term." The dominant source of this systematic effect is thought to be vigorous near-surface convection, which is not well accounted for in both stellar modeling and mode-oscillation physics. Here we bring to bear the method of homogenization, applicable in the asymptotic limit of large wavelengths (in comparison to the correlation scale of convection), to characterize the effect of small-scale surface convection on resonant-mode frequencies in the Sun. We show that the full oscillation equations, in the presence of temporally stationary 3-D flows, can be reduced to an effective "quiet-Sun" wave equation with altered sound speed, Br\"{u}nt--V\"{a}is\"{a}la frequency and Lamb frequency. We derive the modified equation and relations for the appropriate averaging of three dimensional flows and thermal quantities to obtain the properties of this effective medium. Using flows obtained from three dimensional numerical simulations of near-surface convection, we quantify their effect on solar oscillation frequencies, and find that they are shifted systematically and substantially. We argue therefore that consistent interpretations of resonant frequencies must include modifications to the wave equation that effectively capture the impact of vigorous hydrodynamic convection.

Acoustic modes of oscillation are affected by stellar activity, however it is unclear how starspots contribute to these changes. Here we investigate the nonmagnetic effects of starspots on global modes with angular degree $\ell \leq 2$ in highly active stars, and characterize the spot seismic signature on synthetic light curves. We perform 3D time-domain simulations of linear acoustic waves to study their interaction with a model starspot. We model the spot as a 3D change in the sound speed stratification with respect to a convectively stable stellar background, built from solar Model S. We perform a parametric study by considering different depths and perturbation amplitudes. Exact numerical simulations allow the investigation of the wavefield-spot interaction beyond first order perturbation theory. The interaction of the axisymmetric modes with the starspot is strongly nonlinear. As mode frequency increases, the frequency shifts for radial modes exceed the value predicted by linear theory, while the shifts for the $\ell=2, m=0$ modes are smaller than predicted by linear theory, with avoided-crossing-like patterns forming between the $m=0$ and $m=1$ mode frequencies. The nonlinear behavior increases with increasing spot amplitude and/or decreasing depth. Linear theory still reproduces the correct shifts for nonaxisymmetric modes. In the nonlinear regime the mode eigenfunctions are not pure spherical harmonics, but rather a mixture of different spherical harmonics. This mode mixing, together with the frequency changes, may lead to misidentification of the modes in the observed acoustic power spectra.

We discuss recent observational, theoretical and modeling progress made in understanding the Sun's internal dynamics, including its rotation, meridional flow, convection and overshoot. Over the past few decades, substantial theoretical and observational effort has gone into appreciating these aspects of solar dynamics. A review of these observations, related helioseismic methodology and inference and computational results in relation to these problems is undertaken here.

Our Sun, primarily composed of ionized hydrogen and helium, has a surface temperature of 5777~K and a radius $R_\odot \approx 696,000$ km. In the outer $R_\odot/3$, energy transport is accomplished primarily by convection. Using typical convective velocities $u\sim100\,\rm{m\,s^{-1}}$ and kinematic viscosities of order $10^{-4}$ m$^{2}$s$^{-1}$, we obtain a Reynolds number $Re \sim 10^{14}$. Convection is thus turbulent, causing a vast range of scales to be excited. The Prandtl number, $Pr$, of the convecting fluid is very low, of order $10^{-7}$\,--\,$10^{-4}$, so that the Rayleigh number ($\sim Re^2 Pr$) is on the order of $10^{21}\,-\,10^{24}$. Solar convection thus lies in extraordinary regime of dynamical parameters, highly untypical of fluid flows on Earth. Convective processes in the Sun drive global fluid circulations and magnetic fields, which in turn affect its visible outer layers ("solar activity") and, more broadly, the heliosphere ("space weather"). The precise determination of the depth of solar convection zone, departures from adiabaticity of the temperature gradient, and the internal rotation rate as a function of latitude and depth are among the seminal contributions of helioseismology towards understanding convection in the Sun. Contemporary helioseismology, which is focused on inferring the properties of three-dimensional convective features, suggests that transport velocities are substantially smaller than theoretical predictions. Furthermore, helioseismology provides important constraints on the anisotropic Reynolds stresses that control the global dynamics of the solar convection zone. This review discusses the state of our understanding of convection in the Sun, with a focus on helioseismic diagnostics. We present our considerations with the interests of fluid dynamicists in mind.

The inference of flows of material in the interior of the Sun is a subject of major interest in helioseismology. Here we apply techniques of Full Waveform Inversion (FWI) to synthetic data to test flow inversions. In this idealized setup, we do not model seismic realization noise, training the focus entirely on the problem of whether a chosen supergranulation flow model can be seismically recovered. We define the misfit functional as a sum of L_2 norm deviations in travel times between prediction and observation, as measured using short-distance f and p_1 filtered and large-distance unfiltered $p$ modes. FWI allows for the introduction of measurements of choice and iteratively improving the background model, while monitoring the evolution of the misfit in all desired categories. Although the misfit is seen to uniformly reduce in all categories, convergence to the true model is very slow, possibly because it is trapped in a local minimum. The primary source of error is inaccurate depth localization, which, owing to density stratification, leads to wrong ratios of horizontal and vertical flow velocities (`cross talk'). In the present formulation, the lack of sufficient temporal frequency and spatial resolution makes it difficult to accurately localise flow profiles at depth. We therefore suggest that the most efficient way to discover the global minimum is to perform a probabilistic forward search, involving calculating the misfit associated with a broad range of models (generated, for instance, by a Monte-Carlo algorithm) and locating the deepest minimum. Such techniques possess the added advantage of being able to quantify model uncertainty as well as realization noise (data uncertainty).

Recently, Duvall and Hanasoge ({\it Solar Phys.} {\bf 287}, 71-83, 2013) found that large distance $[\Delta]$ separation travel-time differences from a center to an annulus $[\delta t_{\rm{oi}}]$ implied a model of the average supergranular cell that has a peak upflow of $240\rm{ms^{-1}}$ at a depth of $2.3\rm{Mm}$ and a corresponding peak outward horizontal flow of $700\rm{ms^{-1}}$ at a depth of $1.6\rm{Mm}$. In the present work, this effect is further studied by measuring and modeling center-to-quadrant travel-time differences $[\delta t_{\rm{qu}}]$, which roughly agree with this model. Simulations are analyzed that show that such a model flow would lead to the expected travel-time differences. As a check for possible systematic errors, the center-to-annulus travel-time differences $[\delta t_{\rm{oi}}]$ are found not to vary with heliocentric angle. A consistency check finds an increase of $\delta t_{\rm{oi}}$ with the temporal frequency $[\nu]$ by a factor of two, which is not predicted by the ray theory.

Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the solar interior are linearly related to measured travel-time deviations. It is widely known, however, that this assumption is invalid for sunspots and active regions, and likely for supergranular flows as well. Here, we introduce nonlinear optimization, executed iteratively, as a means of inverting for the sub-surface structure of large-amplitude perturbations. Defining the penalty functional as the $L_2$ norm of wave travel-time deviations, we compute the the total misfit gradient of this functional with respect to the relevant model parameters %(only sound speed in this case) at each iteration around the corresponding model. The model is successively improved using either steepest descent, conjugate gradient, or quasi-Newton limited-memory BFGS. Performing nonlinear iterations requires privileging pixels (such as those in the near-field of the scatterer), a practice not compliant with the standard assumption of translational invariance. Measurements for these inversions, although similar in principle to those used in time-distance helioseismology, require some retooling. For the sake of simplicity in illustrating the method, we consider a 2-D inverse problem with only a sound-speed perturbation.

Surface granulation of the Sun is primarily a consequence of thermal transport in the outer 1 % of the radius. Its typical scale of about 1 - 2 Mm is set by the balance between convection, free-streaming radiation, and the strong density stratification in the surface layers. The physics of granulation is well understood, as demonstrated by the close agreement between numerical simulation, theory, and observation. Superimposed on the energetic granular structure comprising high-speed flows, are larger scale long-lived flow systems (~ 300 m/s) called supergranules. Supergranulation has a typical scale of 24 - 36 Mm. It is not clear if supergranulation results from the interaction of granules or is causally linked to deep convection or a consequence of magneto-convection. Other outstanding questions remain: how deep are supergranules? How do they participate in global dynamics of the Sun? Further challenges are posed by our lack of insight into the dynamics of larger scales in the deep convection region. Recent helioseismic constraints have suggested that convective velocity amplitudes on large scales may be overestimated by an order of magnitude or more, implying that Reynolds stresses associated with large-scale convection, thought to play a significant role in the sustenance of differential rotation and meridional circulation, might be two orders of magnitude weaker than theory and computation predict. While basic understanding on the nature of convection on global scales and the maintenance of global circulations is incomplete, progress is imminent, given substantial improvements in computation, theory and helioseismic inferences.

Seismic noise cross correlations are used to image crustal structure and heterogeneity. Typically, seismic networks are only anisotropically illuminated by seismic noise, a consequence of the non-uniform distribution of sources. Here, we study the sensitivity of such a seismic network to structural heterogeneity in a 2-D setting. We compute finite-frequency cross-correlation sensitivity kernels for travel-time, waveform-energy and waveform-difference measurements. In line with expectation, wavespeed anomalies are best imaged using travel times and the source distribution using cross-correlation energies. Perturbations in attenuation and impedance are very difficult to image and reliable inferences require a high degree of certainty in the knowledge of the source distribution and wavespeed model (at least in the case of transmission tomography studied here). We perform single-step Gauss-Newton inversions for the source distribution and the wavespeed, in that order, and quantify the associated Cram\'{e}r-Rao lower bound. The inversion and uncertainty estimate are robust to errors in the source model but are sensitive to the theory used to interpret of measurements. We find that when classical source-receiver kernels are used instead of cross-correlation kernels, errors appear in the both the inversion and uncertainty estimate, systematically biasing the results. We outline a computationally tractable algorithm to account for distant sources when performing inversions.

Rotation is thought to drive cyclic magnetic activity in the Sun and Sun-like stars. Stellar dynamos, however, are poorly understood owing to the scarcity of observations of rotation and magnetic fields in stars. Here, inferences are drawn on the internal rotation of a distant Sun-like star by studying its global modes of oscillation. We report asteroseismic constraints imposed on the rotation rate and the inclination of the spin axis of the Sun-like star HD 52265, a principal target observed by the CoRoT satellite that is known to host a planetary companion. These seismic inferences are remarkably consistent with an independent spectroscopic observation (rotational line broadening) and with the observed rotation period of star spots. Furthermore, asteroseismology constrains the mass of exoplanet HD 52265b. Under the standard assumption that the stellar spin axis and the axis of the planetary orbit coincide, the minimum spectroscopic mass of the planet can be converted into a true mass of 1.85 (+0.52,-0.42) M_Jupiter, which implies that it is a planet, not a brown dwarf.

Seismic noise, generated by oceanic microseisms and other sources, illuminates the crust in a manner different from tectonic sources, and therefore provides independent information. The primary measurable is the two-point cross-correlation, evaluated using traces recorded at a pair of seismometers over a finite-time interval. However, raw seismic traces contain intermittent large-amplitude perturbations arising from tectonic activity and instrumental errors, which may corrupt the estimated cross-correlations of microseismic fluctuations. In order to diminish the impact of these perturbations, the recorded traces are filtered using the nonlinear one-bit digitizer, which replaces the measurement by its sign. Previous theory shows that for stationary Gaussian-distributed seismic noise fluctuations one-bit and raw correlation functions are related by a simple invertible transformation. Here we extend this to show that the simple correspondence between these two correlation techniques remains valid for {\it non-stationary} Gaussian and a very broad range of {\it non-Gaussian} processes as well. For a limited range of stationary and non-stationary Gaussian fluctuations, we find that one-bit filtering performs at least as well as spectral whitening. We therefore recommend using one-bit filtering when processing terrestrial seismic noise, with the substantial benefit that the measurements are fully compatible with current theoretical interpretation (e.g., adjoint theory). Given that seismic records are non-stationary and comprise small-amplitude fluctuations and intermittent, large-amplitude tectonic/other perturbations, we outline an algorithm to accurately retrieve the correlation function of the small-amplitude signals.

Measurements of seismic wave travel times at the photosphere of the Sun have enabled inferences of its interior structure and dynamics. In interpreting these measurements, the simplifying assumption that waves propagate through a temporally stationary medium is almost universally invoked. However, the Sun is in a constant state of evolution, on a broad range of spatio-temporal scales. At the zero wavelength limit, i.e., when the wavelength is much shorter than the scale over which the medium varies, the WKBJ (ray) approximation may be applied. Here, we address the other asymptotic end of the spectrum, the infinite wavelength limit, using the technique of homogenization. We apply homogenization to scenarios where waves are propagating through rapidly varying media (spatially and temporally), and derive effective models for the media. One consequence is that a scalar sound speed becomes a tensorial wavespeed in the effective model and anisotropies can be induced depending on the nature of the perturbation. The second term in this asymptotic two-scale expansion, the so-called corrector, contains contributions due to higher-order scattering, leading to the decoherence of the wavefield. This decoherence may be causally linked to the observed wave attenuation in the Sun. Although the examples we consider here consist of periodic arrays of perturbations to the background, homogenization may be extended to ergodic and stationary random media. This method may have broad implications for the manner in which we interpret seismic measurements in the Sun and for modeling the effects of granulation on the scattering of waves and distortion of normal-mode eigenfunctions.

Supergranulation is one of the most visible length scales of solar convection and has been studied extensively by local helioseismology. We use synthetic data computed with the Seismic Propagation through Active Regions and Convection (SPARC) code to test regularized-least squares (RLS) inversions of helioseismic holography measurements for a supergranulation-like flow. The code simulates the acoustic wavefield by solving the linearized three-dimensional Euler equations in Cartesian geometry. We model a single supergranulation cell with a simple, axisymmetric, mass-conserving flow. The use of simulated data provides an opportunity for direct evaluation of the accuracy of measurement and inversion techniques. The RLS technique applied to helioseismic-holography measurements is generally successful in reproducing the structure of the horizontal flow field of the model supergranule cell. The errors are significant in horizontal-flow inversions near the top and bottom of the computational domain as well as in vertical-flow inversions throughout the domain. We show that the errors in the vertical velocity are due largely to cross talk from the horizontal velocity.

We use analytical examples and asymptotic forms to examine the mathematical structure and physical meaning of the seismic cross correlation measurement. We show that in general, cross correlations are not Green's functions of medium, and may be very different depending on the source distribution. The modeling of noise sources using spatial distributions as opposed to discrete collections of sources is emphasized. When stations are illuminated by spatially complex source distributions, cross correlations show arrivals at a variety of time lags, from zero to the maximum surface-wave arrival time. Here, we demonstrate the possibility of inverting for the source distribution using the energy of the full cross-correlation waveform. The interplay between the source distribution and wave attenuation in determining the functional dependence of cross correlation energies on station-pair distance is quantified. Without question, energies contain information about wave attenuation. However, the accurate interpretation of such measurements is tightly connected to the knowledge of the source distribution.

As large--distance rays (say, 10\,-\,$24 ^\circ$) approach the solar surface approximately vertically, travel times measured from surface pairs for these large separations are mostly sensitive to vertical flows, at least for shallow flows within a few Mm of the solar surface. All previous analyses of supergranulation have used smaller separations and have been hampered by the difficulty of separating the horizontal and vertical flow components. We find that the large separation travel times associated with supergranulation cannot be studied using the standard phase-speed filters of time-distance helioseismology. These filters, whose use is based upon a refractive model of the perturbations, reduce the resultant travel time signal by at least an order of magnitude at some distances. More effective filters are derived. Modeling suggests that the center--annulus travel time difference $[\delta t_{\rm{oi}}]$ in the separation range $\Delta=10$\,-\,$24 ^\circ$ is insensitive to the horizontally diverging flow from the centers of the supergranules and should lead to a constant signal from the vertical flow. Our measurement of this quantity, $5.1 \pm 0.1\secs$, is constant over the distance range. This magnitude of signal cannot be caused by the level of upflow at cell centers seen at the photosphere of $10\ms$ extended in depth. It requires the vertical flow to increase with depth. A simple Gaussian model of the increase with depth implies a peak upward flow of $240\ms$ at a depth of $2.3\Mm$ and a peak horizontal flow of $700\ms$ at a depth of $1.6\Mm$.

Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of changes to isotropic wavespeeds and flow-advection related Doppler shifts. Here, we employ PDE-constrained optimization to determine the appropriate parameterization of the structural properties of the magnetic interior. Seven different wavespeeds fully characterize helioseismic wave propagation: the isotropic sound speed, a Doppler-shifting flow-advection velocity and an anisotropic magnetic velocity. The structure of magnetic media is sensed by magnetoacoustic slow and fast modes and Alfv\'{e}n waves, each of which propagates at a different wavespeed. We show that even in the case of weak magnetic fields, significant errors may be incurred if these anisotropies are not accounted for in inversions. Translation invariance is demonstrably lost. These developments render plausible the accurate seismic imaging of magnetoconvection in the Sun.

Convection in the solar interior is thought to comprise structures on a spectrum of scales. This conclusion emerges from phenomenological studies and numerical simulations, though neither covers the proper range of dynamical parameters of solar convection. Here, we analyze observations of the wavefield in the solar photosphere using techniques of time-distance helioseismology to image flows in the solar interior. We downsample and synthesize 900 billion wavefield observations to produce 3 billion cross-correlations, which we average and fit, measuring 5 million wave travel times. Using these travel times, we deduce the underlying flow systems and study their statistics to bound convective velocity magnitudes in the solar interior, as a function of depth and spherical-harmonic degree $\ell$. Within the wavenumber band $\ell<60$, Convective velocities are 20-100 times weaker than current theoretical estimates. This suggests the prevalence of a different paradigm of turbulence from that predicted by existing models, prompting the question: what mechanism transports the heat flux of a solar luminosity outwards? Advection is dominated by Coriolis forces for wavenumbers $\ell<60$, with Rossby numbers smaller than $\sim10^{-2}$ at $r/R_\odot=0.96$, suggesting that the Sun may be a much faster rotator than previously thought, and that large-scale convection may be quasi-geostrophic. The fact that iso-rotation contours in the Sun are not co-aligned with the axis of rotation suggests the presence of a latitudinal entropy gradient.

Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here we propose a mechanism to search for transits of PBHs through/nearby Earth by studying the associated seismic waves. Using a spectral-element method, we simulate and visualize this seismic wave field in Earth's interior. We predict the emergence of two unique signatures, namely, a wave that would arrive almost simultaneously everywhere on Earth's free surface and the excitation of unusual spheroidal modes with a characteristic frequency-spacing in free oscillation spectra. These qualitative characteristics are unaffected by the speed or proximity of the PBH trajectory. The seismic energy deposited by a proximal ${M^{PBH} = 10^{15}}$ g PBH is comparable to a magnitude $M_w=4$ earthquake. The non-seismic collateral damage due to the actual impact of such small PBHs with Earth would be negligible. Unfortunately, the expected collision rate is very low even if PBHs constituted all of dark matter, at ${\sim 10^{-7} {yr}^{-1}}$, and since the rate scales as ${1/M^{PBH}}$, fortunately encounters with larger, Earth-threatening PBHs are exceedingly unlikely. However, the rate at which non-colliding close encounters of PBHs could be detected by seismic activity alone is roughly two orders of magnitude larger --- that is once every hundred thousand years --- than the direct collision rate.

Stars are transparent to the passage of primordial black holes (PBHs) and serve as seismic detectors for such objects. The gravitational field of a PBH squeezes a star and causes it to ring acoustically. We calculate the seismic signature of a PBH passing through the Sun. The background for this signal is the observed spectrum of solar oscillations excited by supersonic turbulence. We predict that PBHs more massive than 10^21 g (comparable in mass to an asteroid) are detectable by existing solar observatories. The oscillations excited by PBHs peak at large scales and high frequencies, making them potentially detectable in other stars. The discovery of PBHs would have profound implications for cosmology and high-energy physics.

Inversions for local helioseismology are an important and necessary step for obtaining three-dimensional maps of various physical quantities in the solar interior. Frequently, the full inverse problems that one would like to solve prove intractable because of computational constraints. Due to the enormous seismic data sets that already exist and those forthcoming, this is a problem that needs to be addressed. To this end, we present a very efficient linear inversion algorithm for local helioseismology. It is based on a subtractive optimally localized averaging (SOLA) scheme in the Fourier domain, utilizing the horizontal-translation invariance of the sensitivity kernels. In Fourier space the problem decouples into many small problems, one for each horizontal wave vector. This multi-channel SOLA method is demonstrated for an example problem in time-distance helioseismology that is small enough to be solved both in real and Fourier space. We find that both approaches are successful in solving the inverse problem. However, the multi-channel SOLA algorithm is much faster and can easily be parallelized.

We present a new Lattice Boltzmann (LB) formulation to solve the Maxwell equations for electromagnetic (EM) waves propagating in a heterogeneous medium. By using a pseudo-vector discrete Boltzmann distribution, the scheme is shown to reproduce the continuum Maxwell equations. The technique compares well with a pseudo-spectral method at solving for two-dimensional wave propagation in a heterogeneous medium, which by design contains substantial contrasts in the refractive index. The extension to three dimensions follows naturally and, owing to the recognized efficiency of LB schemes for parallel computation in irregular geometries, it gives a powerful method to numerically simulate a wide range of problems involving EM wave propagation in complex media.

For a given {\it misfit function}, a specified optimality measure of a model, its gradient describes the manner in which one may alter properties of the system to march towards a stationary point. The adjoint method, arising from partial-differential-equation-constrained optimization, describes a means of extracting derivatives of a misfit function with respect to model parameters through finite computation. It relies on the accurate calculation of wavefields that are driven by two types of sources, namely the average wave-excitation spectrum, resulting in the {\it forward wavefield}, and differences between predictions and observations, resulting in an {\it adjoint wavefield}. All sensitivity kernels relevant to a given measurement emerge directly from the evaluation of an interaction integral involving these wavefields. The technique facilitates computation of sensitivity kernels (Fr\'{e}chet derivatives) relative to three-dimensional heterogeneous background models, thereby paving the way for non-linear iterative inversions. An algorithm to perform such inversions using as many observations as desired is discussed.

According to time-distance helioseismology, information about internal fluid motions is encoded in the travel times of solar waves. The inverse problem consists of inferring 3-D vector flows from a set of travel-time measurements. Here we investigate the potential of time-distance helioseismology to infer 3-D convective velocities in the near-surface layers of the Sun. We developed a new Subtractive Optimally Localised Averaging (SOLA) code suitable for pipeline pseudo-automatic processing. Compared to its predecessor, the code was improved by accounting for additional constraints in order to get the right answer within a given noise level. The main aim of this study is to validate results obtained by our inversion code. We simulate travel-time maps using a snapshot from a numerical simulation of solar convective flows, realistic Born travel-time sensitivity kernels, and a realistic model of travel-time noise. These synthetic travel times are inverted for flows and the results compared with the known input flow field. Additional constraints are implemented in the inversion: cross-talk minimization between flow components and spatial localization of inversion coefficients. Using modes f, p1 through p4, we show that horizontal convective flow velocities can be inferred without bias, at a signal-to-noise ratio greater than one in the top 3.5 Mm, provided that observations span at least four days. The vertical component of velocity (v_z), if it were to be weak, is more difficult to infer and is seriously affected by cross-talk from horizontal velocity components. We emphasise that this cross-talk must be explicitly minimised in order to retrieve v_z in the top 1 Mm. We also show that statistical averaging over many different areas of the Sun allows for reliably measuring of average properties of all three flow components in the top 5.5 Mm of the convection zone.

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.

Perfectly matched layers are a very efficient and accurate way to absorb waves in media. We present a stable convolutional unsplit perfectly matched formulation designed for the linearized stratified Euler equations. However, the technique as applied to the Magneto-hydrodynamic (MHD) equations requires the use of a sponge, which, despite placing the perfectly matched status in question, is still highly efficient at absorbing outgoing waves. We study solutions of the equations in the backdrop of models of linearized wave propagation in the Sun. We test the numerical stability of the schemes by integrating the equations over a large number of wave periods.

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.

We constrain the velocity spectral distribution of global-scale solar convective cells at depth using techniques of local helioseismology. We calibrate the sensitivity of helioseismic waves to large-scale convective cells in the interior by analyzing simulations of waves propagating through a velocity snapshot of global solar convection via methods of time-distance helioseismology. Applying identical analysis techniques to observations of the Sun, we are able to bound from above the magnitudes of solar convective cells as a function of spatial convective scale. We find that convection at a depth of $r/R_\odot = 0.95$ with spatial extent $\ell <20$, where $\ell$ is the spherical harmonic degree, comprise weak flow systems, on the order of 15 m/s or less. Convective features deeper than $r/R_\odot = 0.95$ are more difficult to image due to the rapidly decreasing sensitivity of helioseismic waves.

Spatial maps of the high-pass frequency filtered time-averaged root-mean-squared (RMS) Doppler velocities tend to show substantial decrements within regions of strong field and curiously, randomly distributed patches of enhancement in the vicinity. We propose that these haloes or enhancements are a consequence of magnetic-field-induced mode mixing (scattering), resulting in the preferential powering of waves that possess strong surface velocity signatures (i.e. scattering from low to high wavenumbers). Evidently, this process can occur in the reverse, and therefore in order to determine if the haloes are indeed caused by mode mixing, we must answer the question: {\it how are acoustic waves scattered by magnetic fields?} Through simulations of the interactions between waves and sunspots and models of plage, we demonstrate that the high to low modal order scattering channels are favoured. With increasing frequency and consequently, decreasing wavelength, a growing number of modes are scattered by the sunspot, thereby rendering the enhancements most visible around the high-frequency parts of the spectrum. The haloes obtained from the simulations are on the same order of magnitude but weaker than those observed. We also present observational evidence to support this theory: observations of active region AR9787 are firstly frequency filtered to isolate the 5-6 mHz signals and secondly, decomposed into three wavenumber bandpasses, $l - [0,400], [400,800], [800,2222]$. With increasing wavenumber, the extent of the halo effect is seen increase dramatically, in line with theoretical expectation.

Measuring the depth variation of the meridional flows is important for understanding the solar cycle, at least according to a number of dynamo models. While attempting to extend the early observations of \citet{giles_thesis} of time-distance measurements of flow, we have stumbled upon some systematic errors that can affect these measurements: 1) the additional distance traveled by radiation coming from points away from disk center causes an apparent `shrinking' Sun, that is an apparent flow towards the disk center and 2) in measurements away from the central longitude, the rotation signal can leak into meridional flow signals. Attempts to understand and overcome these systematic problems will be presented. Forward models based on ray theory have been applied in order to test the sensitivity of travel times to various models.

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 derive expectations for signatures in the measured travel times of waves that interact with thermal anomalies and jets. A series of numerical experiments that involve the dynamic linear evolution of an acoustic wave field in a solar-like stratified spherical shell in the presence of fully 3D time-stationary perturbations are performed. The imprints of these interactions are observed as shifts in wave travel times, which are extracted from these data through methods of time-distance helioseismology \citep{duvall}. In situations where at least one of the spatial dimensions of the scatterer was smaller than a wavelength, oscillatory time shift signals were recovered from the analyses, pointing directly to a means of resolving sub-wavelength features. As evidence for this claim, we present analyses of simulations with spatially localized jets and sound-speed perturbations. We analyze 1 years' worth solar observations to estimate the noise level associated with the time differences. Based on theoretical estimates, Fresnel zone time shifts associated with the (possible) sharp rotation gradient at the base of the convection zone are of the order 0.01 - 0.1 s, well below the noise level that could be reached with the currently available amount of data ($\sim 0.15-0.2$ s with 10 yrs of data).

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.

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.

Comprehending the manner in which magnetic fields affect propagating waves is a first step toward constructing accurate helioseismic models of active region sub-surface structure and dynamics. Here, we present a numerical method to compute the linear interaction of waves with magnetic fields embedded in a solar-like stratified background. The ideal Magneto-Hydrodynamic (MHD) equations are solved in a 3-dimensional box that straddles the solar photosphere, extending from 35 Mm within to 1.2 Mm into the atmosphere. One of the challenges in performing these simulations involves generating a Magneto-Hydro-Static (MHS) state wherein the stratification assumes horizontal inhomogeneity in addition to the strong vertical stratification associated with the near-surface layers. Keeping in mind that the aim of this effort is to understand and characterize linear MHD interactions, we discuss a means of computing statically consistent background states. Power maps computed from simulations of waves interacting with thick flux tubes of peak photospheric field strengths 600 G and 3000 G are presented. Strong modal power reduction in the `umbral' regions of the flux tube enveloped by a halo of increased wave power are seen in the simulations with the thick flux tubes. These enhancements are also seen in Doppler velocity power maps of active regions observed in the Sun, leading us to propose that the halo has MHD underpinnings.

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 study the effect of localized sound-speed perturbations on global mode frequencies by applying techniques of global helioseismology on numerical simulations of the solar acoustic wave field. Extending the method of realization noise subtraction (e.g. Hanasoge et al. 2007) to global modes and exploiting the luxury of full spherical coverage, we are able to achieve very highly resolved frequency differences that are used to study sensitivities and the signatures of the thermal asphericities. We find that (1) global modes are almost twice as sensitive to sound-speed perturbations at the bottom of the convection zone as in comparison to anomalies well in the radiative interior ($r\lesssim0.55 R_\odot$), (2) the $m$-degeneracy is lifted ever so slightly, as seen in the $a$ coefficients, and (3) modes that propagate in the vicinity of the perturbations show small amplitude shifts ($\sim 0.5%$).