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 $\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.Comment: PNAS; 5 figures, 5 page

Time-Distance Helioseismology

Time-distance helioseismology is a method of ambient noise imaging using the solar oscillations. The basic realization that led to time-distance helioseismology was that the temporal cross correlation of the signals at two 'surface' (or photospheric) locations should show a feature at the time lag corresponding to the subsurface travel time between the locations. The temporal cross correlation, as a function of the location separation, is the Fourier transform of the spatio-temporal power spectrum of the solar oscillations, a commonly used function in helioseismology. It is therefore likely the characteristic ridge structure of the correlation function had been seen before without appreciation of its significance. Travel times are measured from the cross correlations. The times are sensitive to a number of important subsurface solar phenomena. These include sound speed variations, flows, and magnetic fields. There has been much interesting progress in the 17 years since the first paper on this subject (Duvall et al., Nature, 1993, 362, 430-432). This progress will be reviewed in this paper

Sub-Wavelength Resolution Imaging of the Solar Deep Interior

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).Comment: Accepted, ApJ; 17 pages, 12 figure

Simulations of Convection Zone Flows and Measurements from Multiple Viewing Angles

A deep-focusing time-distance measurement technique has been applied to linear acoustic simulations of a solar interior perturbed by convective flows. The simulations are for the full sphere for r/R greater than 0.2. From these it is straightforward to simulate the observations from different viewing angles and to test how multiple viewing angles enhance detectibility. Some initial results will be presented

Comparison of acoustic travel-time measurement of solar meridional circulation from SDO/HMI and SOHO/MDI

Time-distance helioseismology is one of the primary tools for studying the solar meridional circulation. However, travel-time measurements of the subsurface meridional flow suffer from a variety of systematic errors, such as a center-to-limb variation and an offset due to the P-angle uncertainty of solar images. Here we apply the time-distance technique to contemporaneous medium-degree Dopplergrams produced by SOHO/MDI and SDO/HMI to obtain the travel-time difference caused by meridional circulation throughout the solar convection zone. The P-angle offset in MDI images is measured by cross-correlating MDI and HMI images. The travel-time measurements in the south-north and east-west directions are averaged over the same observation period for the two data sets and then compared to examine the consistency of MDI and HMI travel times after correcting the systematic errors. The offsets in the south-north travel-time difference from MDI data induced by the P-angle error gradually diminish with increasing travel distance. However, these offsets become noisy for travel distances corresponding to waves that reach the base of the convection zone. This suggests that a careful treatment of the P-angle problem is required when studying a deep meridional flow. After correcting the P-angle and the removal of the center-to-limb effect, the travel-time measurements from MDI and HMI are consistent within the error bars for meridional circulation covering the entire convection zone. The fluctuations observed in both data sets are highly correlated and thus indicate their solar origin rather than an instrumental origin. Although our results demonstrate that the ad hoc correction is capable of reducing the wide discrepancy in the travel-time measurements from MDI and HMI, we cannot exclude the possibility that there exist other systematic effects acting on the two data sets in the same way.Comment: accepted for publication in A&

Solar meridional circulation from twenty-one years of SOHO/MDI and SDO/HMI observations: Helioseismic travel times and forward modeling in the ray approximation

The south-north travel-time differences are measured by applying time-distance helioseismology to the MDI and HMI medium-degree Dopplergrams covering May 1996-April 2017. Our data analysis corrects for several sources of systematic effects: P-angle error, surface magnetic field effects, and center-to-limb variations. An interpretation of the travel-time measurements is obtained using a forward-modeling approach in the ray approximation. The travel-time differences are similar in the southern hemisphere for cycles 23 and 24. However, they differ in the northern hemisphere between cycles 23 and 24. Except for cycle 24's northern hemisphere, the measurements favor a single-cell meridional circulation model where the poleward flows persist down to $\sim$0.8 $R_\odot$, accompanied by local inflows toward the activity belts in the near-surface layers. Cycle 24's northern hemisphere is anomalous: travel-time differences are significantly smaller when travel distances are greater than 20$^\circ$. This asymmetry between northern and southern hemispheres during cycle 24 was not present in previous measurements (e.g., Rajaguru & Antia 2015), which assumed a different P-angle error correction where south-north travel-time differences are shifted to zero at the equator for all travel distances. In our measurements, the travel-time differences at the equator are zero for travel distances less than $\sim$30$^\circ$, but they do not vanish for larger travel distances. This equatorial offset for large travel distances need not be interpreted as a deep cross-equator flow; it could be due to the presence of asymmetrical local flows at the surface near the end points of the acoustic ray paths.Comment: accepted for publication in A&

Global-scale equatorial Rossby waves as an essential component of solar internal dynamics

The Sun's complex dynamics is controlled by buoyancy and rotation in the convection zone and by magnetic forces in the atmosphere and corona. While small-scale solar convection is well understood, the dynamics of large-scale flows in the solar convection zone is not explained by theory or simulations. Waves of vorticity due to the Coriolis force, known as Rossby waves, are expected to remove energy out of convection at the largest scales. Here we unambiguously detect and characterize retrograde-propagating vorticity waves in the shallow subsurface layers of the Sun at angular wavenumbers below fifteen, with the dispersion relation of textbook sectoral Rossby waves. The waves have lifetimes of several months, well-defined mode frequencies below 200 nHz in a co-rotating frame, and eigenfunctions of vorticity that peak at the equator. Rossby waves have nearly as much vorticity as the convection at the same scales, thus they are an essential component of solar dynamics. We find a transition from turbulence-like to wave-like dynamics around the Rhines scale of angular wavenumber of twenty; this might provide an explanation for the puzzling deficit of kinetic energy at the largest spatial scales.Comment: This is the submitted version of the paper published in Nature Astronomy. 23 pages, 8 figures, 1 tabl

Time-distance helioseismology of solar Rossby waves

Context. Solar Rossby waves (r modes) have recently been discovered in the near-surface horizontal flow field using the techniques of granulation-tracking and ring-diagram analysis applied to six years of SDO/HMI data. Aims. Here we apply time-distance helioseismology to the combined SOHO/MDI and SDO/HMI data sets, which cover 21 years of observations from May 1996 to April 2017. The goal of this study is to provide an independent confirmation over two solar cycles and in deeper layers of the Sun. Methods. We have measured south-north helioseismic travel times along the equator, which are sensitive to subsurface north-south flows. To reduce noise, the travel times were averaged over travel distances from 6$^\circ$ to 30$^\circ$; the mean distance corresponds to a p-mode lower turning point of 0.91 $R_\odot$. The 21-year time series of travel-time measurements was split into three seven-year subsets and transformed to obtain power spectra in a corotating frame. Results. The power spectra all show peaks near the frequencies of the classical sectoral Rossby waves for azimuthal wavenumbers in the range $3 \leq m \leq 15$. The mode frequencies and linewidths of the modes with $m \leq 9$ are consistent with a previous study whereas modes with $m \geq 10$ are shifted toward less negative frequencies by 10--20 nHz. While most of these modes have e-folding lifetimes on the order of a few months, the longest lived mode, $m=3$, has an e-folding lifetime of more than one year. For each mode, the rms velocity at the equator is in the range of 1--3 m s$^{-1}$ , with the largest values for $m\sim10$. No evidence for the $m=2$ sectoral mode is found in the power spectrum, implying that the rms velocity of this mode is below $\sim$0.5 m s$^{-1}$.Comment: accepted for publication in A&