Time-distance Helioseismology Research Articles

VERIFICATION OF THE HELIOSEISMOLOGY TRAVEL-TIME MEASUREMENT TECHNIQUE AND THE INVERSION PROCEDURE FOR SOUND SPEED USING ARTIFICIAL DATA

We performed three-dimensional numerical simulations of the solar surface acoustic wave field for the quiet Sun and for three models with different localized sound-speed perturbations in the interior with deep, shallow, and two-layer structures. We used the simulated data generated by two solar acoustics codes that employ the same standard solar model as a background model, but utilize different integration techniques and different models of stochastic wave excitation. Acoustic travel times were measured using a time-distance helioseismology technique, and compared with predictions from ray theory frequently used for helioseismic travel-time inversions. It is found that the measured travel-time shifts agree well with the helioseismic theory for sound-speed perturbations, and for the measurement procedure with and without phase-speed filtering of the oscillation signals. This testing verifies the whole measuring-filtering-inversion procedure for static sound-speed anomalies with small amplitude inside the Sun outside regions of strong magnetic field. It is shown that the phase-speed filtering, frequently used to extract specific wave packets and improve the signal-to-noise ratio, does not introduce significant systematic errors. Results of the sound-speed inversion procedure show good agreement with the perturbation models in all cases. Due to its smoothing nature, the inversion procedure may overestimate sound-speed variations in regions with sharp gradients of the sound-speed profile.

MERIDIONAL FLOW IN THE SOLAR CONVECTION ZONE. I. MEASUREMENTS FROM GONG DATA

Large-scale plasma flows in the Sun's convection zone likely play a major role in solar dynamics on decadal timescales. In particular, quantifying meridional motions is a critical ingredient for understanding the solar cycle and the transport of magnetic flux. Because the signal of such features can be quite small in deep solar layers and be buried in systematics or noise, the true meridional velocity profile has remained elusive. We perform time-distance helioseismology measurements on several years worth of GONG Doppler data. A spherical harmonic decomposition technique is applied to a subset of acoustic modes to measure travel-time differences to try to obtain signatures of meridional flows throughout the solar convection zone. Center-to-limb systematics are taken into account in an intuitive, yet ad hoc manner. Travel-time differences near the surface that are consistent with a poleward flow in each hemisphere and are similar to previous work are measured. Additionally, measurements in deep layers near the base of the convection zone suggest a possible equatorward flow, as well as partial evidence of a sign change in the travel-time differences at mid-convection zone depths. This analysis on an independent data set using different measurement techniques strengthens recent conclusions that the convection zone may have multiple "cells" of meridional flow. The results may challenge the common understanding of one large conveyor belt operating in the solar convection zone. Further work with helioseismic inversions and a careful study of systematic effects are needed before firm conclusions of these large-scale flow structures can be made.

FULL WAVEFORM INVERSION FOR TIME-DISTANCE HELIOSEISMOLOGY

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.

Helioseismology of sunspots: defocusing, folding, and healing of wavefronts

We observe and characterize the scattering of acoustic wave packets by a sunspot in a regime where the wavelength is comparable to the size of the sunspot. Spatial maps of wave travel times and amplitudes are measured from the cross-covariance function of the random wave field observed by SOHO/MDI around the sunspot in active region NOAO 9787. We consider separately incoming plane wave packets consisting of f modes and p modes with radial orders up to four. Observations show that the travel-time perturbations diminish with distance far away from the sunspot – a finite-wavelength phenomenon known as wavefront healing in scattering theory. Observations also show a reduction of the amplitude of the waves after their passage through the sunspot. We suggest that a significant fraction of this amplitude reduction is due to the defocusing of wave energy by the fast wave-speed perturbation introduced by the sunspot. This “geometrical attenuation” will contribute to the wave amplitude reduction in addition to the physical absorption of waves by sunspots. We also observe an enhancement of wave amplitude away from the central path: diffracted rays intersect with unperturbed rays (caustics) and wavefronts fold and triplicate. Wave amplitude measurements in time-distance helioseismology provide independent information that can be used in concert with travel-time measurements.

A Study of Connections Between Solar Flares and Subsurface Flow Fields of Active Regions

We investigate the connections between the occurrence of major solar flares and subsurface dynamic properties of active regions. For this analysis, we select five active regions that produced a total of 11 flares with peak X-ray flux intensity higher than M5.0. The subsurface velocity fields are obtained from time–distance helioseismology analysis using SDO/HMI (Solar Dynamics Observatory/Helioseismic and Magnetic Imager) Doppler observations, and the X-ray flux intensity is taken from GOES (Geostationary Operational Environmental Satellites). It is found that among the eight amplitude bumps in the evolutionary curves of subsurface kinetic helicity, five (62.5%) of them had a flare stronger than M5.0 occurring within 8 hours, either before or after the bumps. Another subsurface parameter is the Normalized Helicity Gradient Variance (NHGV), reflecting kinetic helicity spread in different depth layers; it also shows bumps near the occurrence of these solar flares. Although there is no one-to-one correspondence between the flare and the subsurface properties, these observational phenomena are worth further studies to better understand the flares’ subsurface roots, and to investigate whether the subsurface properties can be used for major flare forecasts.

SOLAR WAVE-FIELD SIMULATION FOR TESTING PROSPECTS OF HELIOSEISMIC MEASUREMENTS OF DEEP MERIDIONAL FLOWS

The meridional flow in the Sun is an axisymmetric flow that is generally poleward directed at the surface, and is presumed to be of fundamental importance in the generation and transport of magnetic fields. Its true shape and strength, however, is debated. We present a numerical simulation of helioseismic wave propagation in the whole solar interior in the presence of a prescribed, stationary, single-cell, deep meridional circulation serving as a test-bed for helioseismic measurement techniques. A deep-focusing time-distance helioseismology technique is applied to the artificial data showing that it can in fact be used to measure the effects of the meridional flow very deep in the solar convection zone. It is shown that the ray-approximation which is commonly used for interpretation of helioseismology measurements remains a reasonable approximation even for the very long distances between 12 and 42 degrees corresponding to depths between 52 and 195 Mm considered here. From the measurement noise we extrapolate that on the order of a full solar cycle may be needed to probe the flow all the way to the base of the convection zone.

Horizontal Flows in the Photosphere and Subphotosphere of Two Active Regions

We compare horizontal flow fields in the photosphere and in the subphotosphere (a layer 0.5 Mm below the photosphere) in two solar active regions: AR 11084 and AR 11158. AR 11084 is a mature, simple active region without significant flaring activity, and AR 11158 is a multipolar, complex active region with magnetic flux emerging during the period studied. Flows in the photosphere are derived by applying the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) on HMI-observed vector magnetic fields, and the subphotospheric flows are inferred by time–distance helioseismology using HMI-observed Dopplergrams. Similar flow patterns are found for both layers for AR 11084: inward flows in the sunspot umbra and outward flows surrounding the sunspot. The boundary between the inward and outward flows, which is slightly different in the photosphere and the subphotosphere, is within the sunspot penumbra. The area having inward flows in the subphotosphere is larger than that in the photosphere. For AR 11158, flows in these two layers show great similarities in some areas and significant differences in other areas. Both layers exhibit consistent outward flows in the areas surrounding sunspots. On the other hand, most well-documented flux-emergence-related flow features seen in the photosphere do not have counterparts in the subphotosphere. This implies that the horizontal flows caused by flux emergence do not extend deeply into the subsurface.

Subsurface Supergranular Vertical Flows as Measured Using Large Distance Separations in Time–Distance Helioseismology

As large-distance rays (say, 10 – 24°) 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 [δt oi] in the separation range Δ=10 – 24∘ 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±0.1 seconds, is constant over the distance range. This magnitude of the signal cannot be caused by the level of upflow at cell centers seen at the photosphere of 10 m s−1 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 m s−1 at a depth of 2.3 Mm and a peak horizontal flow of 700 m s−1 at a depth of 1.6 Mm.

Anomalously weak solar convection.

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 ℓ. Within the wavenumber band ℓ < 60, convective velocities are 20-100 times weaker than current theoretical estimates. This constraint 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 ℓ < 60, with Rossby numbers smaller than approximately 10(-2) at r/R([symbol: see text]) = 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 isorotation contours in the Sun are not coaligned with the axis of rotation suggests the presence of a latitudinal entropy gradient.

SYSTEMATIC CENTER-TO-LIMB VARIATION IN MEASURED HELIOSEISMIC TRAVEL TIMES AND ITS EFFECT ON INFERENCES OF SOLAR INTERIOR MERIDIONAL FLOWS

We report on a systematic center-to-limb variation in measured helioseismic travel times, which must be taken into account for an accurate determination of solar interior meridional flows. The systematic variation, found in time-distance helioseismology analysis using SDO/HMI and SDO/AIA observations, is different in both travel-time magnitude and variation trend for different observables. It is not clear what causes this systematic effect. Subtracting the longitude-dependent east-west travel times, obtained along the equatorial area, from the latitude-dependent north-south travel times, obtained along the central meridian area, gives remarkably similar results for different observables. We suggest this as an effective procedure for removing the systematic center-to-limb variation. The subsurface meridional flows obtained from inversion of the corrected travel times are approximately 10 m/s slower than those obtained without removing the systematic effect. The detected center-to-limb variation may have important implications in the derivation of meridional flows in the deep interior, and needs a better understanding.

HELIOSEISMOLOGY OF A REALISTIC MAGNETOCONVECTIVE SUNSPOT SIMULATION

We compare helioseismic travel-time shifts measured from a realistic magnetoconvective sunspot simulation using both helioseismic holography and time-distance helioseismology, and measured from real sunspots observed with the Helioseismic and Magnetic Imager instrument on board the Solar Dynamics Observatory and the Michelson Doppler Imager instrument on board the Solar and Heliospheric Observatory. We find remarkable similarities in the travel-time shifts measured between the methodologies applied and between the simulated and real sunspots. Forward modeling of the travel-time shifts using either Born or ray approximation kernels and the sound-speed perturbations present in the simulation indicates major disagreements with the measured travel-time shifts. These findings do not substantially change with the application of a correction for the reduction of wave amplitudes in the simulated and real sunspots. Overall, our findings demonstrate the need for new methods for inferring the subsurface structure of sunspots through helioseismic inversions.

Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology

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.

Validated helioseismic inversions for 3D vector flows

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.

Time–Distance Helioseismology Data-Analysis Pipeline for Helioseismic and Magnetic Imager Onboard Solar Dynamics Observatory (SDO/HMI) and Its Initial Results

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

INVERSION OF FULL ACOUSTIC WAVEFIELD IN LOCAL HELIOSEISMOLOGY: A STUDY WITH SYNTHETIC DATA

We present the first results from the inversion of full acoustic wavefield in the helioseismic context. In contrast to time-distance helioseismology, which involves analyzing the travel times of seismic waves propagating into the solar interior, wavefield tomography models both the travel times and amplitude variations present in the entire seismic record. Unlike the use of ray-based, Fresnel-zone, Born, or Rytov approximations in previous time-distance studies, this method does not require any simplifications to be made to the sensitivity kernel in the inversion. In this study, the acoustic wavefield is simulated for all iterations in the inversion. The sensitivity kernel is therefore updated while lateral variations in sound-speed structure in the model emerge during the course of the inversion. Our results demonstrate that the amplitude-based inversion approach is capable of resolving sound-speed structures defined by relatively sharp vertical and horizontal boundaries. This study therefore provides the foundation for a new type of subsurface imaging in local helioseismology that is based on the inversion of the entire seismic wavefield.

Wave absorption and moat flow in AR9787

We present the results of a study of the wave absorption and the moat flow around the sunspot in AR9787 (the "HELAS sunspot", so-called because it has been the object of in-depth study at a HELAS workshop) using Fourier-Hankel decomposition of the oscillatory wave field. We carry out inversions both for the flow speed and for the sound speed perturbation as a function of depth. The results are compared with results from time-distance helioseismology.

HMI time-distance pipeline: An overview and data products

The Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory provides uninterrupted high-resolution observations of solar oscillations over the entire disk. This gives a unique opportunity for mapping subsurface flows and wave-speed structures and investigating their role in the Sun's dynamics and magnetic activity on various scales by methods of local helioseismology. A data analysis pipeline for the time-distance helioseismology analysis has been developed and implemented at the SDO Joint Science Operation Center (JSOC) at Stanford. It provides near-real time processing of the helioseismology data. We provide an overview of this pipeline, including the data flow procedures, measurement and inversion codes, and our data products.

Local helioseismology of sunspot regions: Comparison of ring-diagram and time-distance results

Local helioseismology provides unique information about the subsurface structure and dynamics of sunspots and active regions. However, because of complexity of sunspot regions local helioseismology diagnostics require careful analysis of systematic uncertainties and physical interpretation of the inversion results. We present new results of comparison of the ring-diagram analysis and time-distance helioseismology for active region NOAA 9787, for which a previous comparison showed significant differences in the subsurface sound-speed structure, and discuss systematic uncertainties of the measurements and inversions. Our results show that both the ring-diagram and time-distance techniques give qualitatively similar results, revealing a characteristic two-layer seismic sound-speed structure consistent with the results for other active regions. However, a quantitative comparison of the inversion results is not straightforward. It must take into account differences in the sensitivity, spatial resolution and the averaging kernels. In particular, because of the acoustic power suppression, the contribution of the sunspot seismic structure to the ring-diagram signal can be substantially reduced. We show that taking into account this effect reduces the difference in the depth of transition between the negative and positive sound-speed variations inferred by these methods. Further detailed analysis of the sensitivity, resolution and averaging properties of the local helioseismology methods is necessary for consolidation of the inversion results. It seems to be important that both methods indicate that the seismic structure of sunspots is rather deep and extends to at least 20 Mm below the surface, putting constraints on theoretical models of sunspots.

AN ESTIMATE OF THE DETECTABILITY OF RISING FLUX TUBES

The physics of the formation of magnetic active regions (ARs) is one of the most important problems in solar physics. One main class of theories suggests that ARs are the result of magnetic flux that rises from the tachocline. Time–distance helioseismology, which is based on measurements of wave propagation, promises to allow the study of the subsurface behavior of this magnetic flux. Here, we use a model for a buoyant magnetic flux concentration together with the ray approximation to show that the dominant effect on the wave propagation is expected to be from the roughly 100 m s −1 retrograde flow associated with the rising flux. Using a B-spline-based method for carrying out inversions of wave travel times for flows in spherical geometry, we show that at 3 days before emergence the detection of this retrograde flow at a depth of 30 Mm should be possible with a signal-to-noise level of about 8 with a sample of 150 emerging ARs.

Investigation of a sunspot complex by time-distance helioseismology

AbstractSunspot regions often form complexes of activity that may live for several solar rotations, and represent a major component of the Sun's magnetic activity. It had been suggested that the close appearance of active regions in space and time might be related to common subsurface roots, or “nests” of activity. EUV images show that the active regions are magnetically connected in the corona, but subsurface connections have not been established. We investigate the subsurface structure and dynamics of a large complex of activity, NOAA 10987-10989, observed during the SOHO/MDI Dynamics run in March-April 2008, which was a part of the Whole Heliospheric Interval (WHI) campaign. The active regions in this complex appeared in a narrow latitudinal range, probably representing a subsurface toroidal flux tube. We use the MDI full-disk Dopplergrams to measure perturbations of travel times of acoustic waves traveling to various depths by using time-distance helioseismology, and obtain sound-speed and flow maps by inversion of the travel times. The subsurface flow maps show an interesting dynamics of decaying active regions with persistent shearing flows, which may be important for driving the flaring and CME activity, observed during the WHI campaign. Our analyses, including the seismic sound-speed inversion results and the distribution of deep-focus travel-time anomalies, gave indications of diverging roots of the magnetic structures, as could be expected from Ω-loop structures. However, no clear connection in the depth range of 0-48 Mm among the three active regions in this complex of activity was detected.