Search

The Sun's global inertial modes are very sensitive to the solar differential rotation and to properties of the deep solar convection zone which are currently poorly constrained. These properties include the superadiabatic temperature gradient, the latitudinal entropy gradient, and the turbulent viscosity. The inertial modes also play a key role in controlling the Sun's large-scale structure and dynamics, in particular the solar differential rotation. This paper summarizes recent observations and advances in the (linear and nonlinear) modeling of the solar inertial modes., Comment: Submitted on 05-Mar-2024 for publication in procs. IAU Symp 365, Eds. A. Getling and L. Kitchatinov

Emerging active regions are associated with convective flows on the spatial scale and lifetimes of supergranules. To understand how these flows are involved in the formation of active regions, we aim to identify where active regions emerge in the supergranulation flow pattern. We computed supergranulation scale flow maps at the surface for all active regions in the Solar Dynamics Observatory Helioseismic Emerging Active Region Survey. We classified each of the active regions into four bins, based on the amplitude of their average surface flow divergence at emergence. We then averaged the flow divergence over the active regions in each bin as a function of time. We also considered a corresponding set of control regions. We found that, on average, the flow divergence increases during the day prior to emergence at a rate independent of the amount of flux that emerges. By subtracting the averaged flow divergence of the control regions, we found that active region emergence is associated with a remaining converging flow at 0.5-1 days prior to emergence. This remnant flow, $\Delta \, \mathrm{div} \, \mathbf{v_h} = (-4.9 \pm 1.7) \times 10^{-6}$ 1/s, corresponds to a flow speed of 10-20 m/s (an order of magnitude less than supergranulation flows) out to a radius of about 10 Mm. We show that these observational results are qualitatively supported by simulations of a small bipole emerging through the near-surface convective layers of the Sun. The question remains whether these flows are driving the emergence, or are caused by the emergence., Comment: 22 pages, 27 figures, accepted to be published in MNRAS

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

Solaris is a transformative Solar Polar Discovery-class mission concept to address crucial outstanding questions that can only be answered from a polar vantage. Solaris will image the Sun's poles from ~75 degree latitude, providing new insight into the workings of the solar dynamo and the solar cycle, which are at the foundation of our understanding of space weather and space climate. Solaris will also provide enabling observations for improved space weather research, modeling and prediction, revealing a unique, new view of the corona, coronal dynamics and CME eruptions from above., Comment: This White Paper was submitted in 2022 to the United States National Academies Solar and Space Physics (Heliophysics) Decadal Survey

Active regions (ARs) play an important role in the magnetic dynamics of the Sun. Solar surface flux transport models (SFTMs) are used to describe the evolution of the radial magnetic field at the solar surface. There is however uncertainty about using these models in the early stage of AR evolution. We aim to test the applicability of SFTMs in the first days after the emergence of ARs by comparing them with observations. The models we employ range from passive evolution to models where the inflows around ARs are included. We simulate the evolution of the surface magnetic field of 17 emerging active regions using a local surface flux transport simulation. We selected regions that do not form fully-fledged sunspots that exhibit moat flows. The simulation includes diffusion and advection. We use observed flows from local correlation tracking of solar granulation, as well as parametrizations of the inflows around ARs. To evaluate our simulations, we measure the cross correlation between the observed and the simulated magnetic field, as well as the total unsigned flux of the ARs, over time. We also test the validity of our simulations by varying the starting time relative to the emergence of flux. We find that the simulations using observed surface flows can reproduce the evolution of the observed magnetic flux. The effect of buffeting of the field by supergranulation can be described as a diffusion process. The SFTM is applicable after 90% of the peak total unsigned flux of the AR has emerged. Diffusivities in the range between $D=250$ to $720$ km$^2$/s are consistent with the evolution of the AR flux in the first five days after this time. We find that the converging flows around emerging ARs are not important for the evolution of the total flux of the AR in these first five days; their effect of increasing flux cancellation is balanced by the decrease of flux transport away from the AR., Comment: Accepted for publication in Astronomy & Astrophysics. 13 pages, 11 figures

The oscillations of a slowly rotating star have long been classified into spheroidal and toroidal modes. The spheroidal modes include the well-known 5-min acoustic modes used in helioseismology. Here we report observations of the Sun's toroidal modes, for which the restoring force is the Coriolis force and whose periods are on the order of the solar rotation period. By comparing the observations with the normal modes of a differentially rotating spherical shell, we are able to identify many of the observed modes. These are the high-latitude inertial modes, the critical-latitude inertial modes, and the equatorial Rossby modes. In the model, the high-latitude and critical-latitude modes have maximum kinetic energy density at the base of the convection zone, and the high-latitude modes are baroclinically unstable due to the latitudinal entropy gradient. As a first application of inertial-mode helioseismology, we constrain the superadiabaticity and the turbulent viscosity in the deep convection zone., Comment: Forthcoming article, Astronomy and Astrophysics

A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by any other mission. One of the major mysteries of the Sun is the solar cycle. The activity cycle of the Sun drives the structure and behaviour of the heliosphere and is, of course, the driver of space weather. In addition, solar activity and variability provides fluctuating input into the Earth climate models, and these same physical processes are applicable to stellar systems hosting exoplanets. One of the main obstructions to understanding the solar cycle, and hence all solar activity, is our current lack of understanding of the polar regions. In this White Paper, submitted to the European Space Agency in response to the Voyage 2050 call, we describe a mission concept that aims to address this fundamental issue. In parallel, we recognise that viewing the Sun from above the polar regions enables further scientific advantages, beyond those related to the solar cycle, such as unique and powerful studies of coronal mass ejection processes, from a global perspective, and studies of coronal structure and activity in polar regions. Not only will these provide important scientific advances for fundamental stellar physics research, they will feed into our understanding of impacts on the Earth and other planets' space environment., Comment: 29 pages, 9 figures, White paper for ESA Voyage 2050

Understanding convection is important in stellar physics, for example as an input in stellar evolution models. Helioseismic estimates of convective flow amplitudes in deeper regions of the solar interior disagree by orders of magnitude among themselves and with simulations. We aim to assess the validity of an existing upper limit of solar convective flow amplitudes at a depth of 0.96 solar radii obtained using time-distance helioseismology and several simplifying assumptions. We generated synthetic observations for convective flow fields from a magnetohydrodynamic simulation (MURaM) using travel-time sensitivity functions and a noise model. We compared the estimates of the flow with the actual values. For the scales of interest ($\ell<100$), we find that the current procedure for obtaining an upper limit gives the correct order of magnitude of the flow for the given flow fields. We also show that this estimate is not an upper limit in a strict sense because it underestimates the flow amplitude at the largest scales by a factor of about two because the scale dependence of the signal-to-noise ratio has to be taken into account. After correcting for this and after taking the dependence of the measurements on direction in Fourier space into account, we show that the obtained estimate is indeed an upper limit. We conclude that time-distance helioseismology is able to correctly estimate the order of magnitude (or an upper limit) of solar convective flows in the deeper interior when the vertical correlation function of the different flow components is known and the scale dependence of the signal-to-noise ratio is taken into account. We suggest that future work should include information from different target depths to better separate the effect of near-surface flows from those at greater depths. The measurements are sensitive to all three flow directions, which should be taken into account., Comment: 16 pages, 13 figures, accepted for publication in A&A

Context: Giant-cell convection is believed to be an important component of solar dynamics. For example, it is expected to play a crucial role in maintaining the Sun's differential rotation. Aims: We re-examine early reports of giant convective cells detected using correlation analysis of Dopplergrams. We extend this analysis using 19 years of space and ground-based observations of near-surface horizontal flows. Method: Flow maps are derived through local correlation tracking of granules and helioseismic ring-diagram analysis. We compute temporal auto-correlation functions of the east-west flows at fixed latitude. Results: Correlations in the east-west velocity can be clearly seen up to five rotation periods. The signal consists of features with longitudinal wavenumbers up to m=9 at low latitudes. Comparison with magnetic images indicates that these flow features are associated with magnetic activity. The signal is not seen above the noise level during solar minimum. Conclusions: Our results show that the long-term correlations in east-west flows at low latitudes are dominantly due to inflows into active regions and not to giant convective cells., Comment: 9 pages, 10 figures, accepted into A&A

Magnetic flux generated within the solar interior emerges to the surface, forming active regions (ARs) and sunspots. Flux emergence may trigger explosive events, such as flares and coronal mass ejections and therefore understanding emergence is useful for space-weather forecasting. Evidence of any pre-emergence signatures will also shed light on sub-surface processes responsible for emergence. In this paper, we present a first analysis of emerging ARs from the Solar Dynamics Observatory/Helioseismic Emerging Active Regions (SDO/HEAR) dataset (Schunker et al. 2016) using deep convolutional neural networks (CNN) to characterize pre-emergence surface magnetic-field properties. The trained CNN classifies between pre-emergence (PE) line-of-sight magnetograms and a control set of non-emergence (NE) magnetograms with a True Skill Statistic (TSS) score of ~85%, 3h prior to emergence and ~40\%, 24h prior to emergence. Our results are better than a baseline classification TSS obtained using discriminant analysis of only the unsigned magnetic flux. We develop a network pruning algorithm to interpret the trained CNN and show that the CNN incorporates filters that respond positively as well as negatively to the unsigned magnetic flux of the magnetograms. Using synthetic magnetograms, we demonstrate that the CNN output is sensitive to the length-scale of the magnetic regions with small-scale and intense fields producing maximum CNN output and possibly a characteristic pre-emergence pattern. Given increasing popularity of deep learning, techniques developed here for interpretation of the trained CNN -- using network pruning and synthetic data -- are relevant for future applications in solar and astrophysical data analysis., Comment: 24 pages, 24 figures, accepted for publication in ApJ

The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes Joy's law, acting on east-west flows as they rise towards the surface. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we can constrain the origins of Joy's law. We measured the tilt angle of the primary bipoles in 153 emerging active regions in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey. We used line-of-sight magnetic field measurements averaged over 6 hours to define the polarities and measure the tilt angle up to four days after emergence. We find that at the time of emergence the polarities are on average aligned east-west, and that neither the separation nor the tilt depends on latitude. We do find, however, that ARs at higher latitudes have a faster north-south separation speed than those closer to the equator at the emergence time. After emergence, the tilt angle increases and Joy's law is evident about two days later. The scatter in the tilt angle is independent of flux until about one day after emergence, when higher-flux regions have a smaller scatter in tilt angle than lower-flux regions. Our finding that active regions emerge with an east-west alignment is surprising since thin flux tube models predict that tilt angles of rising flux tubes are generated below the surface. Previously reported tilt angle relaxation of deeply anchored flux tubes can be largely explained by the change in east-west separation. We conclude that Joy's law is caused by an inherent north-south separation speed present when the flux first reaches the surface, and that the scatter in the tilt angle is consistent with buffeting of the polarities by supergranulation., Comment: accepted in Astronomy and Astrophysics

Context. The spatial power spectrum of supergranulation does not fully characterize the underlying physics of turbulent convection. For example, it does not describe the non-Gaussianity in the horizontal flow divergence. Aims. Our aim is to statistically characterize the spatial pattern of solar supergranulation beyond the power spectrum. The next-order statistic is the bispectrum. It measures correlations of three Fourier components and is related to the nonlinearities in the underlying physics. Methods. We estimated the bispectrum of supergranular horizontal surface divergence maps that were obtained using local correlation tracking (LCT) and time-distance helioseismology (TD) from one year of data from the Helioseismic and Magnetic Imager on-board the Solar Dynamics Observatory starting in May 2010. Results. We find significantly nonzero and consistent estimates for the bispectrum. The strongest nonlinearity is present when the three coupling wave vectors are at the supergranular scale. These are the same wave vectors that are present in regular hexagons, which were used in analytical studies of solar convection. At these Fourier components, the bispectrum is positive, consistent with the positive skewness in the data and with supergranules preferentially consisting of outflows surrounded by a network of inflows. We use the bispectrum to generate synthetic divergence maps that are very similar to the data by a model that consists of a Gaussian term and a weaker quadratic nonlinear component. Thereby, we estimate the fraction of the variance in the divergence maps from the nonlinear component to be of the order of 4-6%. Conclusions. We propose that bispectral analysis is useful for understanding solar turbulent convection, for example for comparing observations and numerical models of supergranular flows. This analysis may also be useful to generate synthetic flow fields., Comment: 16 pages, 12 figures, accepted for publication by A&A

Context: There is a wide discrepancy in current estimates of the strength of convection flows in the solar interior obtained using different helioseismic methods applied to observations from SDO/HMI. The cause for these disparities is not known. Aims: As one step in the effort to resolve this discrepancy, we aim to characterize the multi-ridge fitting code for ring-diagram helioseismic analysis that is used to obtain flow estimates from local power spectra of solar oscillations. Methods: We updated the multi-ridge fitting code developed by Greer et al.(2014) to solve several problems we identified through our inspection of the code. In particular, we changed the merit function to account for the smoothing of the power spectra, model for the power spectrum, and noise estimates. We used Monte Carlo simulations to generate synthetic data and to characterize the noise and bias of the updated code by fitting these synthetic data. Results: The bias in the output fit parameters, apart from the parameter describing the amplitude of the p-mode resonances in the power spectrum, is below what can be measured from the Monte-Carlo simulations. The amplitude parameters are underestimated; this is a consequence of choosing to fit the logarithm of the averaged power. We defer fixing this problem as it is well understood and not significant for measuring flows in the solar interior. The scatter in the fit parameters from the Monte-Carlo simulations is well-modeled by the formal error estimates from the code. Conclusions: We document and demonstrate a reliable multi-ridge fitting method for ring-diagram analysis. The differences between the updated fitting results and the original results are less than one order of magnitude and therefore we suspect that the changes will not eliminate the aforementioned orders-of-magnitude discrepancy in the amplitude of convective flows in the solar interior., Comment: Accepted for publication in A&A

Our goal is to constrain models of active region formation by tracking the average motion of active region polarity pairs as they emerge onto the surface. We measured the motion of the two main opposite polarities in 153 emerging active regions (EARs) using line-of-sight magnetic field observations from the Solar Dynamics Observatory Helioseismic Emerging Active Region (SDO/HEAR) survey (Schunker et al. 2016). We first measured the position of each of the polarities eight hours after emergence and tracked their location forwards and backwards in time. We find that, on average, the polarities emerge with an east-west orientation and the separation speed between the polarities increases. At about 0.1 days after emergence, the average separation speed reaches a peak value of 229 +/- 11 m/s, and then starts to decrease, and about 2.5 days after emergence the polarities stop separating. We also find that the separation and the separation speed in the east-west direction are systematically larger for active regions with higher flux. Our results reveal two phases of the emergence process defined by the rate of change of the separation speed as the polarities move apart. Phase 1 begins when the opposite polarity pairs first appear at the surface, with an east-west alignment and an increasing separation speed. We define Phase 2 to begin when the separation speed starts to decrease, and ends when the polarities have stopped separating. This is consistent with the picture of Chen, Rempel, & Fan (2017): the peak of a flux tube breaks through the surface during Phase 1. During Phase 2 the magnetic field lines are straightened by magnetic tension, so that the polarities continue to move apart, until they eventually lie directly above their anchored subsurface footpoints., Comment: accepted A&A

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&A

Helioseismic holography is an imaging technique used to study heterogeneities and flows in the solar interior from observations of solar oscillations at the surface. Holograms contain noise due to the stochastic nature of solar oscillations. We provide a theoretical framework for modeling signal and noise in Porter-Bojarski helioseismic holography. The wave equation may be recast into a Helmholtz-like equation, so as to connect with the acoustics literature and define the holography Green's function in a meaningful way. Sources of wave excitation are assumed to be stationary, horizontally homogeneous, and spatially uncorrelated. Using the first Born approximation we calculate holograms in the presence of perturbations in sound-speed, density, flows, and source covariance, as well as the noise level as a function of position. This work is a direct extension of the methods used in time-distance helioseismology to model signal and noise. To illustrate the theory, we compute the hologram intensity numerically for a buried sound-speed perturbation at different depths in the solar interior. The reference Green's function is obtained for a spherically-symmetric solar model using a finite-element solver in the frequency domain. Below the pupil area on the surface, we find that the spatial resolution of the hologram intensity is very close to half the local wavelength. For a sound-speed perturbation of size comparable to the local spatial resolution, the signal-to-noise ratio is approximately constant with depth. Averaging the hologram intensity over a number $N$ of frequencies above 3 mHz increases the signal-to-noise ratio by a factor nearly equal to the square root of $N$. This may not be the case at lower frequencies, where large variations in the holographic signal are due to the individual contributions of the long-lived modes of oscillation., Comment: Submitted to Astronomy and Astrophysics

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&A

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 table

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., Comment: Accepted for publication in in A&A, 11 pages, 6 figures

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&A

Context. In time-distance helioseismology, wave travel times are measured from the two-point cross-covariance function of solar oscillations and are used to image the solar convection zone in three dimensions. There is, however, also information in the amplitude of the cross-covariance function, for example about seismic wave attenuation. Aims. Here we develop a convenient procedure to measure the amplitude of the cross-covariance function of solar oscillations. Methods. In this procedure, the amplitude of the cross-covariance function is linearly related to the cross-covariance function and can be measured even for high levels of noise. Results. As an example application, we measure the amplitude perturbations of the seismic waves that propagate through the sunspot in active region NOAA 9787. We can recover the amplitude variations due to the scattering and attenuation of the waves by the sunspot and associated finite-wavelength effects. Conclusions. The proposed definition of cross-covariance amplitude is robust to noise, can be used to relate measured amplitudes to 3D perturbations in the solar interior under the Born approximation, and will provide independent information from the travel times., Comment: accepted for publication in A&A

Stars are not perfectly spherically symmetric. They are deformed by rotation and magnetic fields. Until now, the study of stellar shapes has only been possible with optical interferometry for a few of the fastest-rotating nearby stars. We report an asteroseismic measurement, with much better precision than interferometry, of the asphericity of an A-type star with a rotation period of 100 days. Using the fact that different modes of oscillation probe different stellar latitudes, we infer a tiny but significant flattening of the star's shape of $\Delta R/R = (1.8 \pm 0.6) \times 10^{-6}$. For a stellar radius $R$ that is $2.24$ times the solar radius, the difference in radius between the equator and the poles is $\Delta R = 3 \pm 1$ km. Because the observed $\Delta R/R$ is only one-third of the expected rotational oblateness, we conjecture the presence of a weak magnetic field on a star that does not have an extended convective envelope. This calls to question the origin of the magnetic field., Comment: 15 pages, 2 tables, 3 figures

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

Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube, it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around emerging active regions, in combination with numerical simulations of solar magnetoconvection, can constrain the subsurface rise speed of emerging magnetic flux. The observed flows imply that the rise speed of the magnetic field is no larger than 150 m/s at a depth of 20 Mm, that is, well below the prediction of the (standard) thin flux tube model but in the range expected for convective velocities at this depth. We conclude that convective flows control the dynamics of rising flux tubes in the upper layers of the Sun and cannot be neglected in models of flux emergence., Comment: Published in Science Advances. DOI: 10.1126/sciadv.1600557

Context: The cross-covariance of solar oscillations observed at pairs of points on the solar surface is a fundamental ingredient in time-distance helioseismology. Wave travel times are extracted from the cross-covariance function and are used to infer the physical conditions in the solar interior. Aims: Understanding the statistics of the two-point cross-covariance function is a necessary step towards optimizing the measurement of travel times. Methods: By modeling stochastic solar oscillations, we evaluate the variance of the cross-covariance function as function of time-lag and distance between the two points. Results: We show that the variance of the cross-covariance is independent of both time-lag and distance in the far field, i.e., when they are large compared to the coherence scales of the solar oscillations. Conclusions: The constant noise level for the cross-covariance means that the signal-to-noise ratio for the cross-covariance is proportional to the amplitude of the expectation value of the cross-covariance. This observation is important for planning data analysis efforts., Comment: accepted for publication in A&A

We extend an existing Born approximation method for calculating the linear sensitivity of helioseismic travel times to flows from Cartesian to spherical geometry. This development is necessary for using the Born approximation for inferring large-scale flows in the deep solar interior. In a first sanity check, we compare two $f-$mode kernels from our spherical method and from an existing Cartesian method. The horizontal and total integrals agree to within 0.3 %. As a second consistency test, we consider a uniformly rotating Sun and a travel distance of 42 degrees. The analytical travel-time difference agrees with the forward-modelled travel-time difference to within 2 %. In addition, we evaluate the impact of different choices of filter functions on the kernels for a meridional travel distance of 42 degrees. For all filters, the sensitivity is found to be distributed over a large fraction of the convection zone. We show that the kernels depend on the filter function employed in the data analysis process. If modes of higher harmonic degree ($90\lesssim l \lesssim 170$) are permitted, a noisy pattern of a spatial scale corresponding to $l\approx 260$ appears near the surface. When mainly low-degree modes are used ($l\lesssim70$), the sensitivity is concentrated in the deepest regions and it visually resembles a ray-path-like structure. Among the different low-degree filters used, we find the kernel for phase-speed filtered measurements to be best localized in depth., Comment: 17 pages, 5 figures, 2 tables, accepted for publication in ApJ. v2: typo in arXiv author list corrected

Context. Several upcoming and proposed space missions, such as Solar Orbiter, will be limited in telemetry and thus require data compression. Aims. We test the impact of data compression on local correlation tracking (LCT) of time-series of continuum intensity images. We evaluate the effect of several lossy compression methods (quantization, JPEG compression, and a reduced number of continuum images) on measurements of solar differential rotation with LCT. Methods. We apply the different compression methods to tracked and remapped continuum intensity maps obtained by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory. We derive 2D vector velocities using the local correlation tracking code FLCT and determine the additional bias and noise introduced by compression to differential rotation. Results. We find that probing differential rotation with LCT is very robust to lossy data compression when using quantization. Our results are severely affected by systematic errors of the LCT method and the HMI instrument. The sensitivity of LCT to systematic errors is a concern for Solar Orbiter., Comment: 9 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

Time-distance helioseismology uses cross-covariances of wave motions on the solar surface to determine the travel times of wave packets moving from one surface location to another. We review the methodology to interpret travel-time measurements in terms of small, localized perturbations to a horizontally homogeneous reference solar model. Using the first Born approximation, we derive and compute 3D travel-time sensitivity (Fr\'echet) kernels for perturbations in sound-speed, density, pressure, and vector flows. While kernels for sound speed and flows had been computed previously, here we extend the calculation to kernels for density and pressure, hence providing a complete description of the effects of solar dynamics and structure on travel times. We treat three thermodynamic quantities as independent and do not assume hydrostatic equilibrium. We present a convenient approach to computing damped Green's functions using a normal-mode summation. The Green's function must be computed on a wavenumber grid that has sufficient resolution to resolve the longest lived modes. The typical kernel calculations used in this paper are computer intensive and require on the order of 600 CPU hours per kernel. Kernels are validated by computing the travel-time perturbation that results from horizontally-invariant perturbations using two independent approaches. At fixed sound-speed, the density and pressure kernels are approximately related through a negative multiplicative factor, therefore implying that perturbations in density and pressure are difficult to disentangle. Mean travel-times are not only sensitive to sound-speed, density and pressure perturbations, but also to flows, especially vertical flows. Accurate sensitivity kernels are needed to interpret complex flow patterns such as 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., Comment: 7 pages, 5 figures, accepted for publication in A&A on 19 March 2015

Context. Several upcoming helioseismology space missions are very limited in telemetry and will have to perform extensive data compression. This requires the development of new methods of data compression. Aims. We give an overview of the influence of lossy data compression on local helioseismology. We investigate the effects of several lossy compression methods (quantization, JPEG compression, and smoothing and subsampling) on power spectra and time-distance measurements of supergranulation flows at disk center. Methods. We applied different compression methods to tracked and remapped Dopplergrams obtained by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. We determined the signal-to-noise ratio of the travel times computed from the compressed data as a function of the compression efficiency. Results. The basic helioseismic measurements that we consider are very robust to lossy data compression. Even if only the sign of the velocity is used, time-distance helioseismology is still possible. We achieve the best results by applying JPEG compression on spatially subsampled data. However, our conclusions are only valid for supergranulation flows at disk center and may not be valid for all helioseismology applications., Comment: 6 pages, 6 figures. Published in Astronomy & Astrophysics

Time-distance helioseismology provides information about vector flows in the near-surface layers of the Sun by measuring wave travel times between points on the solar surface. Specific spatial averages of travel times have been proposed for distinguishing between flows in the east-west and north-south directions and measuring the horizontal divergence of the flows. No specific measurement technique has, however, been developed to measure flow vorticity. Here we propose a new measurement technique tailored to measuring the vertical component of vorticity. Fluid vorticity is a fundamental property of solar convection zone dynamics and of rotating turbulent convection in particular. The method consists of measuring the travel time of waves along a closed contour on the solar surface in order to approximate the circulation of the flow along this contour. Vertical vorticity is related to the difference between clockwise and counter-clockwise travel times. We applied the method to characterize the vortical motions of solar convection using helioseismic data from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) and from the Michelson Doppler Imager onboard the Solar and Heliospheric Observatory (SOHO/MDI). Away from the equator, a clear correlation between vertical vorticity and horizontal divergence is detected. Horizontal outflows are associated with negative vorticity in the northern hemisphere and positive vorticity in the southern hemisphere. The signal is much stronger for HMI than for MDI observations. We characterize the spatial power spectrum of the signal by comparison with a noise model. Vertical vorticity at horizontal wavenumbers below 250/R_Sun can be probed with this helioseismic technique., Comment: 7 pages, 6 figures

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

In time-distance helioseismology, information about the solar interior is encoded in measurements of travel times between pairs of points on the solar surface. Travel times are deduced from the cross-covariance of the random wave field. Here we consider travel times and also products of travel times as observables. They contain information about e.g. the statistical properties of convection in the Sun. The basic assumption of the model is that noise is the result of the stochastic excitation of solar waves, a random process which is stationary and Gaussian. We generalize the existing noise model (Gizon and Birch 2004) by dropping the assumption of horizontal spatial homogeneity. Using a recurrence relation, we calculate the noise covariance matrices for the moments of order 4, 6, and 8 of the observed wave field, for the moments of order 2, 3 and 4 of the cross-covariance, and for the moments of order 2, 3 and 4 of the travel times. All noise covariance matrices depend only on the expectation value of the cross-covariance of the observed wave field. For products of travel times, the noise covariance matrix consists of three terms proportional to $1/T$, $1/T^2$, and $1/T^3$, where $T$ is the duration of the observations. For typical observation times of a few hours, the term proportional to $1/T^2$ dominates and $Cov[\tau_1 \tau_2, \tau_3 \tau_4] \approx Cov[\tau_1, \tau_3] Cov[\tau_2, \tau_4] + Cov[\tau_1, \tau_4] Cov[\tau_2, \tau_3]$, where the $\tau_i$ are arbitrary travel times. This result is confirmed for $p_1$ travel times by Monte Carlo simulations and comparisons with SDO/HMI observations. General and accurate formulae have been derived to model the noise covariance matrix of helioseismic travel times and products of travel times. These results could easily be generalized to other methods of local helioseismology, such as helioseismic holography and ring diagram analysis.

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe I 6173.3 {\AA} line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multi-height velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams. We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30 - 40 km above (line-center) and 30 - 40 km below (average-wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic cutoff frequency., Comment: accepted for publication in Solar Physics

We study Doppler velocity measurements at multiple heights in the solar atmosphere using a set of six filtergrams obtained by the Helioseismic magnetic Imager on board the Solar Dynamics Observatory. There are clear and significant phase differences between core and wing Dopplergrams in the frequency range above the photospheric acoustic cutoff frequency, which indicates that these are really "multi-height" datasets., Comment: 4 pages, 3 figures, in the proceedings "Progress in Physics of the Sun and Stars: A New Era in Helio- and Asteroseismology" (ASP Conference Series)

Linear time-domain simulations of acoustic oscillations are unstable in the stellar convection zone. To overcome this problem it is customary to compute the oscillations of a stabilized background stellar model. The stabilization, however, affects the result. Here we propose to use a perturbative approach (running the simulation twice) to approximately recover the acoustic wave field, while preserving seismic reciprocity. To test the method we considered a 1D standard solar model. We found that the mode frequencies of the (unstable) standard solar model are well approximated by the perturbative approach within $1~\mu$Hz for low-degree modes with frequencies near $3~\mu$Hz. We also show that the perturbative approach is appropriate for correcting rotational-frequency kernels. Finally, we comment that the method can be generalized to wave propagation in 3D magnetized stellar interiors because the magnetic fields have stabilizing effects on convection., Comment: 14 pages. Published online in Solar Physics, available at http://link.springer.com/article/10.1007/s11207-013-0457-7

In order to assess the ability of helioseismology to probe the subsurface structure and magnetic field of sunspots, we need to determine how helioseismic travel times depend on perturbations to sunspot models. Here we numerically simulate the propagation of f, p1, and p2 wave packets through magnetic sunspot models. Among the models we considered, a ~50 km change in the height of the Wilson depression and a change in the subsurface magnetic field geometry can both be detected above the observational noise level. We also find that the travel-time shifts due to changes in a sunspot model must be modeled by computing the effects of changing the reference sunspot model, and not by computing the effects of changing the subsurface structure in the quiet-Sun model. For p1 modes the latter is wrong by a factor of four. In conclusion, numerical modeling of MHD wave propagation is an essential tool for the interpretation of the effects of sunspots on seismic waveforms., Comment: 5 pages, 3 figures: submitted to A&A

Convection is the mechanism by which energy is transported through the outermost 30% of the Sun. Solar turbulent convection is notoriously difficult to model across the entire convection zone where the density spans many orders of magnitude. In this issue of PNAS, Hanasoge et al. (2012) employ recent helioseismic observations to derive stringent empirical constraints on the amplitude of large-scale convective velocities in the solar interior. They report an upper limit that is far smaller than predicted by a popular hydrodynamic numerical simulation., Comment: Printed in the Proceedings of the National Academy of Sciences (2 pages, 1 figure). Available at http://www.pnas.org/cgi/doi/10.1073/pnas.1208875109

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., Comment: 4 pages, 4 figures, accepted Physical Review Letters

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., Comment: 23 pages, 9 Figures, accepted by ApJ