Showing 1–50 of 77 results for author: Hanasoge, S

Detectability of axisymmetric magnetic fields from the core to the surface of oscillating post-main sequence stars

Authors: Shatanik Bhattacharya, Srijan Bharati Das, Lisa Bugnet, Subrata Panda, Shravan M. Hanasoge

Abstract: Magnetic fields in the stellar interiors are key candidates to explain observed core rotation rates inside solar-like stars along their evolution. Recently, asteroseismic estimates of radial magnetic field amplitudes near the hydrogen-burning shell (H-shell) inside about 24 red-giants (RGs) have been obtained by measuring frequency splittings from their power spectra. Using general Lorentz-stress… ▽ More Magnetic fields in the stellar interiors are key candidates to explain observed core rotation rates inside solar-like stars along their evolution. Recently, asteroseismic estimates of radial magnetic field amplitudes near the hydrogen-burning shell (H-shell) inside about 24 red-giants (RGs) have been obtained by measuring frequency splittings from their power spectra. Using general Lorentz-stress (magnetic) kernels, we investigated the potential for detectability of near-surface magnetism in a 1.3 $M_{\odot}$ star of super-solar metallicity as it evolves from a mid sub-giant to a late sub-giant into an RG. Based on these sensitivity kernels, we decompose an RG into three zones - deep core, H-shell, and near-surface. The sub-giants instead required decomposition into an inner core, an outer core, and a near-surface layer. Additionally, we find that for a low-frequency g-dominated dipolar mode in the presence of a typical stable magnetic field, ~25% of the frequency shift comes from the H-shell and the remaining from deeper layers. The ratio of the subsurface tangential field to the radial field in H-burning shell decides if subsurface fields may be potentially detectable. For p-dominated dipole modes close to $νにゅー_\rm{max}$, this ratio is around two orders of magnitude smaller in subgiant phases than the corresponding RG. Further, with the availability of magnetic kernels, we propose lower limits of field strengths in crucial layers in our stellar model during its evolutionary phases. The theoretical prescription outlined here provides the first formal way to devise inverse problems for stellar magnetism and can be seamlessly employed for slow rotators. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Accepted for publication in The Astrophysical Journal, 34 pages, 22 figures, 3 tables

Probing Depth Variations of Solar Inertial Modes through Normal Mode Coupling

Authors: Krishnendu Mandal, Shravan M. Hanasoge

Abstract: Recently discovered inertial waves, observed on the solar surface, likely extend to the deeper layers of the Sun. Utilizing helioseismic techniques, we explore these motions, allowing us to discern inertial-mode eigenfunctions in both radial and latitudinal orientations. We analyze $8$ years of space-based observations ($2010 - 2017$) taken by the Helioseismic and Magnetic Imager (HMI) onboard the… ▽ More Recently discovered inertial waves, observed on the solar surface, likely extend to the deeper layers of the Sun. Utilizing helioseismic techniques, we explore these motions, allowing us to discern inertial-mode eigenfunctions in both radial and latitudinal orientations. We analyze $8$ years of space-based observations ($2010 - 2017$) taken by the Helioseismic and Magnetic Imager (HMI) onboard the Solar dynamic observatory (SDO) using normal-mode coupling. Coupling between same and different-degree acoustic modes and different frequency bins are measured in order to capture the various length scales of inertial modes. We detect inertial modes at high latitude with azimuthal order $t=1$ and frequency $\sim -80$ nHz. This mode is present in the entire convection zone. The presence of Rossby modes may be seen down to a depth of $\sim 0.83R_\odot$ and the Rossby signal is indistinguishable from noise below that depth for high azimuthal order. We find that the amplitudes of these modes increase with depth down to around $0.92 R_\odot$ and decrease below that depth. We find that the latitudinal eigenfunctions of Rossby modes deviate from sectoral spherical harmonics if we use a similar approach as adopted in earlier studies. We found that spatial leakage and even pure noise in the measurements of non-sectoral components can also explain the above-mentioned characteristics of the latitudinal eigenfunctions. This realization underscores the necessity for careful interpretation when considering the latitudinal eigenfunctions of Rossby modes. Exploring the depth-dependent characteristics of these modes will enable us to capture interior dynamics distinctly, separate from p-mode seismology. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 17 pages, 9 figures, submitted to ApJ

Magnetic flux in the Sun emerges unaffected by supergranular-scale surface flows

Authors: Prasad Mani, Chris S. Hanson, Siddarth Dhanpal, Shravan Hanasoge, Srijan Bharati Das, Matthias Rempel

Abstract: Magnetic flux emergence from the convection zone into the photosphere and beyond is a critical component of the behaviour of large-scale solar magnetism. Flux rarely emerges amid field-free areas at the surface, but when it does, the interaction between magnetism and plasma flows can be reliably explored. Prior ensemble studies identified weak flows forming near emergence locations, but the low si… ▽ More Magnetic flux emergence from the convection zone into the photosphere and beyond is a critical component of the behaviour of large-scale solar magnetism. Flux rarely emerges amid field-free areas at the surface, but when it does, the interaction between magnetism and plasma flows can be reliably explored. Prior ensemble studies identified weak flows forming near emergence locations, but the low signal-to-noise required averaging over the entire dataset, erasing information about variation across the sample. Here, we apply deep learning to achieve improved signal-to-noise, enabling a case-by-case study. We find that these associated flows are dissimilar across instances of emergence and also occur frequently in the quiet convective background. Our analysis suggests diminished influence of supergranular-scale convective flows and magnetic buoyancy on flux rise. Consistent with numerical evidence, we speculate that small-scale surface turbulence and / or deep-convective processes play an outsize role in driving flux emergence. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 20 pages, 13 figures, Accepted in ApJ

Asteroseismology applied to constrain structure parameters of δでるた Scuti stars

Authors: Subrata Kumar Panda, Siddharth Dhanpal, Simon J. Murphy, Shravan Hanasoge, Timothy R. Bedding

Abstract: Asteroseismology is a powerful tool to probe the structure of stars. Space-borne instruments like CoRoT, Kepler and TESS have observed the oscillations of numerous stars, among which δでるた Scutis are particularly interesting owing to their fast rotation rates and complex pulsation mechanisms. In this work, we inferred model-dependent masses, metallicities and ages of 60 δでるた Scuti stars from their photom… ▽ More Asteroseismology is a powerful tool to probe the structure of stars. Space-borne instruments like CoRoT, Kepler and TESS have observed the oscillations of numerous stars, among which δでるた Scutis are particularly interesting owing to their fast rotation rates and complex pulsation mechanisms. In this work, we inferred model-dependent masses, metallicities and ages of 60 δでるた Scuti stars from their photometric, spectroscopic and asteroseismic observations using least-squares minimization. These statistics have the potential to explain why only a tiny fraction of δでるた Sct stars pulsate in a very clean manner. We find most of these stars with masses around 1.6 {M_\odot} and metallicities below Z = 0.010. We observed a bimodality in age for these stars, with more than half the sample younger than 30 Myr, while the remaining ones were inferred to be older, i.e., hundreds of Myrs. This work emphasizes the importance of the large-frequency separation ({Δでるたνにゅー}) in studies of δでるた Scuti stars. We also designed three machine learning (ML) models that hold the potential for inferring these parameters at lower computational cost and much more rapidly. These models further revealed that constraining dipole modes can help in significantly improving age estimation and that radial modes succinctly encode information pertaining to stellar luminosity and temperature. Using the ML models, we also gained qualitative insight into the importance of stellar observables in estimating mass, metallicity, and age. The effective surface temperature T eff strongly affects the inference of all structure parameters and the asteroseismic offset parameter εいぷしろん plays an essential role in the inference of age. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Inferring coupling strengths of mixed-mode oscillations in red-giant stars using deep learning

Authors: Siddharth Dhanpal, Othman Benomar, Shravan Hanasoge, Masao Takata, Subrata Panda, Abhisek Kundu

Abstract: Asteroseismology is a powerful tool that may be applied to shed light on stellar interiors and stellar evolution. Mixed modes, behaving as acoustic waves in the envelope and buoyancy modes in the core, are remarkable because they allow for probing the radiative cores and evanescent zones of red-giant stars. Here, we have developed a neural network that can accurately infer the coupling strength, a… ▽ More Asteroseismology is a powerful tool that may be applied to shed light on stellar interiors and stellar evolution. Mixed modes, behaving as acoustic waves in the envelope and buoyancy modes in the core, are remarkable because they allow for probing the radiative cores and evanescent zones of red-giant stars. Here, we have developed a neural network that can accurately infer the coupling strength, a parameter related to the size of the evanescent zone, of solar-like stars in $\sim$5 milliseconds. In comparison with existing methods, we found that only $\sim$43\% inferences were in agreement to within a difference of 0.03 on a sample of $\sim$1,700 \textit{Kepler} red giants. To understand the origin of these differences, we analyzed a few of these stars using independent techniques such as the Monte Carlo Markov Chain method and Echelle diagrams. Through our analysis, we discovered that these alternate techniques are supportive of the neural-net inferences. We also demonstrate that the network can be used to yield estimates of coupling strength and large period separation in stars with structural discontinuities. Our findings suggest that the rate of decline in the coupling strength in the red-giant branch is greater than previously believed. These results are in closer agreement with calculations of stellar-evolution models than prior estimates, further underscoring the remarkable success of stellar-evolution theory and computation. Additionally, we show that the uncertainty in measuring large-period separation increases rapidly with diminishing coupling strength. △ Less

Submitted 2 October, 2023; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 11 figures, Accepted in The Astrophysical Journal

A linear model for inertial modes in a differentially rotating Sun

Authors: Jishnu Bhattacharya, Chris S. Hanson, Shravan M. Hanasoge, Katepalli R. Sreenivasan

Abstract: Inertial wave modes in the Sun are of interest owing to their potential to reveal new insight into the solar interior. These predominantly retrograde-propagating modes in the solar subsurface appear to deviate from the thin-shell Rossby-Haurwitz model at high azimuthal orders. We present new measurements of sectoral equatorial inertial modes at $m>15$ where the modes appear to become progressively… ▽ More Inertial wave modes in the Sun are of interest owing to their potential to reveal new insight into the solar interior. These predominantly retrograde-propagating modes in the solar subsurface appear to deviate from the thin-shell Rossby-Haurwitz model at high azimuthal orders. We present new measurements of sectoral equatorial inertial modes at $m>15$ where the modes appear to become progressively less retrograde compared to the canonical Rossby-Haurwitz dispersion relation in a co-rotating frame. We use a spectral eigenvalue solver to compute the spectrum of solar inertial modes in the presence of differential rotation. Focussing specifically on equatorial Rossby modes, we find that the numerically obtained mode frequencies lie along distinct ridges, one of which lies strikingly close to the observed mode frequencies in the Sun. We also find that the $n=0$ ridge is deflected strongly in the retrograde direction. This suggests that the solar measurements may not correspond to the fundamental $n=0$ Rossby-Haurwitz solutions as was initially suspected, but to a those for a higher $n$. The numerically obtained eigenfunctions also appear to sit deep within the convection zone -- unlike those for the $n=0$ modes -- which differs substantially from solar measurements and complicates inference. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 16 pages, 12 figures

Minuscule corrections to near-surface solar internal rotation using mode-coupling

Authors: Srijan Bharati Das, Samarth G. Kashyap, Deniz Oktay, Shravan M. Hanasoge, Jeroen Tromp

Abstract: The observed solar oscillation spectrum is influenced by internal perturbations such as flows and structural asphericities. These features induce splitting of characteristic frequencies and distort the resonant-mode eigenfunctions. Global axisymmertric flow -- differential rotation -- is a very prominent perturbation. Tightly constrained rotation profiles as a function of latitude and radius are p… ▽ More The observed solar oscillation spectrum is influenced by internal perturbations such as flows and structural asphericities. These features induce splitting of characteristic frequencies and distort the resonant-mode eigenfunctions. Global axisymmertric flow -- differential rotation -- is a very prominent perturbation. Tightly constrained rotation profiles as a function of latitude and radius are products of established helioseismic pipelines that use observed Dopplergrams to generate frequency-splitting measurements at high precision. However, the inference of rotation using frequency-splittings do not consider the effect of mode-coupling. This approximation worsens for high-angular-degree modes, as they become increasingly proximal in frequency. Since modes with high angular degrees probe the near-surface layers of the Sun, inversions considering coupled modes could potentially lead to more accurate estimates of rotation very close to the surface. In order to investigate if this is indeed the case, we perform inversions for solar differential rotation, considering coupling of modes for angular degrees $160 \leq \ell \leq 300$ in the surface gravity $f$-branch and first-overtone $p$ modes. In keeping with the character of mode coupling, we carry out a non-linear inversion using an eigenvalue solver. Differences in inverted profiles for frequency splitting measurements from MDI and HMI are compared and discussed. We find that corrections to the near-surface differential rotation profile, when accounting for mode-coupling effects, are smaller than 0.003 nHz and hence are insignificant. These minuscule corrections are found to be correlated with the solar cycle. We also present corrections to even-order splitting coefficients, which could consequently impact inversions for structure and magnetic fields. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 21 pages; 15 figures; Accepted for publication in the Astrophysical Journal Supplement Series

Exploring the Solar Poles: The Last Great Frontier of the Sun

Authors: Dibyendu Nandy, Dipankar Banerjee, Prantika Bhowmik, Allan Sacha Brun, Robert H. Cameron, S. E. Gibson, Shravan Hanasoge, Louise Harra, Donald M. Hassler, Rekha Jain, Jie Jiang, Laurène Jouve, Duncan H. Mackay, Sushant S. Mahajan, Cristina H. Mandrini, Mathew Owens, Shaonwita Pal, Rui F. Pinto, Chitradeep Saha, Xudong Sun, Durgesh Tripathi, Ilya G. Usoskin

Abstract: Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodyna… ▽ More Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere. △ Less

Submitted 30 December, 2022; originally announced January 2023.

Comments: This White Paper was submitted in 2022 to the United States National Academies Solar and Space Physics (Heliophysics) Decadal Survey

A spectral solver for solar inertial waves

Authors: Jishnu Bhattacharya, Shravan M. Hanasoge

Abstract: Inertial waves, which are dominantly driven by the Coriolis force, likely play an important role in solar dynamics, and additionally, provide a window into the solar subsurface. The latter allows us to infer properties that are inaccessible to the traditional technique of acoustic-wave helioseismology. Thus, a full characterization of these normal modes holds promise in enabling the investigation… ▽ More Inertial waves, which are dominantly driven by the Coriolis force, likely play an important role in solar dynamics, and additionally, provide a window into the solar subsurface. The latter allows us to infer properties that are inaccessible to the traditional technique of acoustic-wave helioseismology. Thus, a full characterization of these normal modes holds promise in enabling the investigation of solar subsurface dynamics. In this work, we develop a spectral eigenvalue solver to model the spectrum of inertial waves in the Sun. We model the solar convection zone as an anelastic medium, and solve for the normal modes of the momentum and energy equations. We demonstrate that the solver can reproduce the observed mode frequencies and line-widths well, not only of sectoral Rossby modes, but also the recently observed high-frequency inertial modes. In addition, we believe that the spectral solver is a useful contribution to the numerical methods on modeling inertial modes on the Sun. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: 6 Figures, accepted for publication in ApJS

Deep learning reconstruction of sunspot vector magnetic fields for forecasting solar storms

Authors: Dattaraj B. Dhuri, Shamik Bhattacharjee, Shravan M. Hanasoge, Sashi Kiran Mahapatra

Abstract: Solar magnetic activity produces extreme solar flares and coronal mass ejections, which pose grave threats to electronic infrastructure and can significantly disrupt economic activity. It is therefore important to appreciate the triggers of explosive solar activity and develop reliable space-weather forecasting. Photospheric vector-magnetic-field data capture sunspot magnetic-field complexity and… ▽ More Solar magnetic activity produces extreme solar flares and coronal mass ejections, which pose grave threats to electronic infrastructure and can significantly disrupt economic activity. It is therefore important to appreciate the triggers of explosive solar activity and develop reliable space-weather forecasting. Photospheric vector-magnetic-field data capture sunspot magnetic-field complexity and can therefore improve the quality of space-weather prediction. However, state-of-the-art vector-field observations are consistently only available from Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) since 2010, with most other current and past missions and observational facilities such as Global Oscillations Network Group (GONG) only recording line-of-sight (LOS) fields. Here, using an inception-based convolutional neural network, we reconstruct HMI sunspot vector-field features from LOS magnetograms of HMI as well as GONG with high fidelity (~ 90% correlation) and sustained flare-forecasting accuracy. We rebuild vector-field features during the 2003 Halloween storms, for which only LOS-field observations are available, and the CNN-estimated electric-current-helicity accurately captures the observed rotation of the associated sunspot prior to the extreme flares, showing a striking increase. Our study thus paves the way for reconstructing three solar cycles worth of vector-field data from past LOS measurements, which are of great utility in improving space-weather forecasting models and gaining new insights about solar activity. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 19 Pages, 11 Figures, Accepted for publication in ApJ

Measuring frequency and period separations in red-giant stars using machine learning

Authors: Siddharth Dhanpal, Othman Benomar, Shravan Hanasoge, Abhisek Kundu, Dattaraj Dhuri, Dipankar Das, Bharat Kaul

Abstract: Asteroseismology is used to infer the interior physics of stars. The \textit{Kepler} and TESS space missions have provided a vast data set of red-giant light curves, which may be used for asteroseismic analysis. These data sets are expected to significantly grow with future missions such as \textit{PLATO}, and efficient methods are therefore required to analyze these data rapidly. Here, we describ… ▽ More Asteroseismology is used to infer the interior physics of stars. The \textit{Kepler} and TESS space missions have provided a vast data set of red-giant light curves, which may be used for asteroseismic analysis. These data sets are expected to significantly grow with future missions such as \textit{PLATO}, and efficient methods are therefore required to analyze these data rapidly. Here, we describe a machine learning algorithm that identifies red giants from the raw oscillation spectra and captures \textit{p} and \textit{mixed} mode parameters from the red-giant power spectra. We report algorithmic inferences for large frequency separation ($Δでるたνにゅー$), frequency at maximum amplitude ($νにゅー_{max}$), and period separation ($ΔでるたΠぱい$) for an ensemble of stars. In addition, we have discovered $\sim$25 new probable red giants among 151,000 \textit{Kepler} long-cadence stellar-oscillation spectra analyzed by the method, among which four are binary candidates which appear to possess red-giant counterparts. To validate the results of this method, we selected $\sim$ 3,000 \textit{Kepler} stars, at various evolutionary stages ranging from subgiants to red clumps, and compare inferences of $Δでるたνにゅー$, $ΔでるたΠぱい$, and $νにゅー_{max}$ with estimates obtained using other techniques. The power of the machine-learning algorithm lies in its speed: it is able to accurately extract seismic parameters from 1,000 spectra in $\sim$5 seconds on a modern computer (single core of the Intel Xeon Platinum 8280 CPU). △ Less

Submitted 20 February, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Imaging the Sun's near-surface flows using mode-coupling analysis

Authors: Prasad Mani, Chris Hanson, Shravan Hanasoge

Abstract: The technique of normal-mode coupling is a powerful tool with which to seismically image non-axisymmetric phenomena in the Sun. Here we apply mode coupling in the Cartesian approximation to probe steady, near-surface flows in the Sun. Using Doppler cubes obtained from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory, we perform inversions on mode-coupling measurements to… ▽ More The technique of normal-mode coupling is a powerful tool with which to seismically image non-axisymmetric phenomena in the Sun. Here we apply mode coupling in the Cartesian approximation to probe steady, near-surface flows in the Sun. Using Doppler cubes obtained from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory, we perform inversions on mode-coupling measurements to show that the resulting divergence and radial vorticity maps at supergranular length scales ($\sim$30 Mm) near the surface compare extremely well with those obtained using the Local Correlation Tracking method. We find that the Pearson correlation coefficient is $\geq$ 0.9 for divergence flows, while $\geq$ 0.8 is obtained for the radial vorticity. △ Less

Submitted 1 January, 2022; originally announced January 2022.

Comments: Accepted for publication in Astrophysical Journal

Investigating toroidal flows in the Sun using normal-mode coupling

Authors: Prasad Mani, Shravan Hanasoge

Abstract: Helioseismic observations have provided valuable datasets with which to pursue the detailed investigation of solar interior dynamics. Among various methods to analyse these data, normal-mode coupling has proven to be a powerful tool, used to study Rossby waves, differential rotation, meridional circulation, and non-axisymmetric multi-scale subsurface flows. Here, we invert mode-coupling measuremen… ▽ More Helioseismic observations have provided valuable datasets with which to pursue the detailed investigation of solar interior dynamics. Among various methods to analyse these data, normal-mode coupling has proven to be a powerful tool, used to study Rossby waves, differential rotation, meridional circulation, and non-axisymmetric multi-scale subsurface flows. Here, we invert mode-coupling measurements from Helioseismic Magnetic Imager (HMI) and Michelson Doppler Imager (MDI) to obtain mass-conserving toroidal convective flow as a function of depth, spatial wavenumber, and temporal frequency. To ensure that the estimates of velocity magnitudes are proper, we also evaluate correlated realization noise, caused by the limited visibility of the Sun. We benchmark the near-surface inversions against results from Local Correlation Tracking (LCT). Convective power likely assumes greater latitudinal isotropy with decrease in spatial scale of the flow. We note an absence of a peak in toroidal-flow power at supergranular scales, in line with observations that show that supergranulation is dominantly poloidal in nature. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 10 pages, 4 figures, forthcoming article in Astrophysical Journal Letters

Detection of Rossby modes with even azimuthal orders using helioseismic normal-mode coupling

Authors: K. Mandal, S. M. Hanasoge, L. Gizon

Abstract: Retrograde Rossby waves, measured to have significant amplitudes in the Sun, likely have notable implications for various solar phenomena. Rossby waves create small-amplitude, very-low frequency motions (on the order of the rotation rate and lower), which in turn shift the resonant frequencies and eigenfunctions of the acoustic modes of the Sun. The detection of even azimuthal orders Rossby modes… ▽ More Retrograde Rossby waves, measured to have significant amplitudes in the Sun, likely have notable implications for various solar phenomena. Rossby waves create small-amplitude, very-low frequency motions (on the order of the rotation rate and lower), which in turn shift the resonant frequencies and eigenfunctions of the acoustic modes of the Sun. The detection of even azimuthal orders Rossby modes using mode coupling presents additional challenges and prior work therefore only focused on odd orders. Here, we successfully extend the methodology to measure even azimuthal orders as well. We analyze 4 and 8 years of Helioseismic and Magnetic Imager (HMI) data and consider coupling between different-degree acoustic modes (of separations 1 and 3 in harmonic degree). The technique uses couplings between different frequency bins to capture the temporal variability of the Rossby modes. We observe significant power close to the theoretical dispersion relation for sectoral Rossby modes (where the azimuthal order is same as harmonic degree, s = |t|). Our results are consistent with prior measurements of Rossby modes with azimuthal orders over the range t = 4 to 16 with maximum power occurring at mode t = 8. The amplitudes of these modes vary from 1 to 2 m/s. We place an upper bound of 0.2 m/s on the sectoral t = 2 mode, which we do not detect in our measurements. This effort adds credence to the mode-coupling methodology in helioseismology △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Accepted for publication in Astronomy & Astrophysics. 9 pages, 5 figures

Journal ref: A&A 652, A96 (2021)

Characterising solar surface convection using Doppler measurements

Authors: Samarth G Kashyap, Shravan M Hanasoge

Abstract: The Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) records line-of-sight Dopplergram images of convective flows on the surface. These images are used to obtain the multi-scale convective spectrum. We design a pipeline to process the raw images to remove large-scale features like differential rotation, meridional circulation, limb shift and imaging artefacts. The Hi… ▽ More The Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) records line-of-sight Dopplergram images of convective flows on the surface. These images are used to obtain the multi-scale convective spectrum. We design a pipeline to process the raw images to remove large-scale features like differential rotation, meridional circulation, limb shift and imaging artefacts. The Hierarchical Equal Area Pixelization scheme (HEALPix) is used to perform spherical harmonic transforms on the cleaned image. Because we only have access to line-of-sight velocities on half the solar surface, we define a "mixing matrix" to relate the observed and true spectra. This enables the inference of poloidal and toroidal flow spectra in a single step through the inversion of the mixing matrix. Performing inversions on a number of flow profiles, we find that the poloidal flow recovery is most reliable among all the components. We also find that the poloidal spectrum is in qualitative agreement with inferences from Local Correlation Tracking (LCT) of granules. The fraction of power in vertical motions increases as a function of wavenumber and is at 8% level for $\ell = 1500$. In contrast to seismic results and LCT, the flows show nearly no temporal-frequency dependence. Poloidal flow power peaks in the range of $\ell - |m| \approx 150-250$ which may potentially hint at a latitudinal preference for convective flows. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 21 pages, 23 figures

Journal ref: ApJ 916 (2021), 87

Analyzing supergranular power spectra using helioseismic normal-mode coupling

Authors: Chris S. Hanson, Shravan Hanasoge, Katepalli R. Sreenivasan

Abstract: Normal-mode coupling is a technique applied to probe the solar interior using surface observations of oscillations. The technique, which is straightforward to implement, makes more use of the seismic information in the wavefield than other comparable local imaging techniques and therefore has the potential to significantly improve current capabilities. Here, we examine supergranulation power spect… ▽ More Normal-mode coupling is a technique applied to probe the solar interior using surface observations of oscillations. The technique, which is straightforward to implement, makes more use of the seismic information in the wavefield than other comparable local imaging techniques and therefore has the potential to significantly improve current capabilities. Here, we examine supergranulation power spectra using mode-coupling analyses of intermediate-to-high-degree modes by invoking a Cartesian-geometric description of wave propagation under the assumption that the localized patches are much smaller in size than the solar radius. We extract the supergranular power spectrum and compare the results with prior helioseismic studies. Measurements of the dispersion relation and life times of supergranulation, obtained using near surface modes (f and p$_1$), are in accord with the literature. We show that the cross-coupling between the p$_2$ and p$_3$ acoustic modes, which are capable of probing greater depths, are also sensitive to supergranulation. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 7 pages, 3 figures, Accepted into ApJ

Inferring solar differential rotation through normal-mode coupling using Bayesian statistics

Authors: Samarth G. Kashyap, Srijan Bharati Das, Shravan M. Hanasoge, Martin F. Woodard, Jeroen Tromp

Abstract: Normal-mode helioseismic data analysis uses observed solar oscillation spectra to infer perturbations in the solar interior due to global and local-scale flows and structural asphericity. Differential rotation, the dominant global-scale axisymmetric perturbation, has been tightly constrained primarily using measurements of frequency splittings via "$a$-coefficients". However, the frequency-splitti… ▽ More Normal-mode helioseismic data analysis uses observed solar oscillation spectra to infer perturbations in the solar interior due to global and local-scale flows and structural asphericity. Differential rotation, the dominant global-scale axisymmetric perturbation, has been tightly constrained primarily using measurements of frequency splittings via "$a$-coefficients". However, the frequency-splitting formalism invokes the approximation that multiplets are isolated. This assumption is inaccurate for modes at high angular degrees. Analysing eigenfunction corrections, which respect cross coupling of modes across multiplets, is a more accurate approach. However, applying standard inversion techniques using these cross-spectral measurements yields $a$-coefficients with a significantly wider spread than the well-constrained results from frequency splittings. In this study, we apply Bayesian statistics to infer $a$-coefficients due to differential rotation from cross spectra for both $f$-modes and $p$-modes. We demonstrate that this technique works reasonably well for modes with angular degrees $\ell=50-291$. The inferred $a_3-$coefficients are found to be within $1$ nHz of the frequency splitting values for $\ell > 200$. We also show that the technique fails at $\ell < 50$ owing to the insensitivity of the measurement to the perturbation. These results serve to further establish mode coupling as an important helioseismic technique with which to infer internal structure and dynamics, both axisymmetric (e.g., meridional circulation) and non-axisymmetric perturbations. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 20 pages, 11 figures

Journal ref: ApJS 253 (2021) 47

A general formulation for computing spherical helioseismic sensitivity kernels while incorporating systematical effects

Authors: Jishnu Bhattacharya, Shravan M. Hanasoge, Katepalli R. Sreenivasan

Abstract: As helioseismology matures and turns into a precision science, modeling finite-frequency, geometric and systematical effects is becoming increasingly important. Here we introduce a general formulation for treating perturbations of arbitrary tensor rank in spherical geometry using fundamental ideas of quantum mechanics and their extensions in geophysics. We include line-of-sight projections and cen… ▽ More As helioseismology matures and turns into a precision science, modeling finite-frequency, geometric and systematical effects is becoming increasingly important. Here we introduce a general formulation for treating perturbations of arbitrary tensor rank in spherical geometry using fundamental ideas of quantum mechanics and their extensions in geophysics. We include line-of-sight projections and center-to-limb differences in line-formation heights in our analysis. We demonstrate the technique by computing a travel-time sensitivity kernel for sound-speed perturbations. The analysis produces the spherical harmonic coefficients of the sensitivity kernels, which leads to better-posed and computationally efficient inverse problems. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Journal ref: ApJ 895 117 (2020)

Application and interpretation of deep learning for identifying pre-emergence magnetic-field patterns

Authors: Dattaraj B. Dhuri, Shravan M. Hanasoge, Aaron C. Birch, Hannah Schunker

Abstract: 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… ▽ More 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. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 24 pages, 24 figures, accepted for publication in ApJ

Validating inversions for toroidal flows using normal-mode coupling

Authors: Prasad Subramanian, Shravan Hanasoge

Abstract: Normal-mode coupling is a helioseismic technique that uses measurements of mode eigenfunctions to infer the interior structure of the Sun. This technique has led to insights into the evolution and structure of toroidal flows in the solar interior. Here, we validate an inversion algorithm for normal-mode coupling by generating synthetic seismic measurements associated with input flows and comparing… ▽ More Normal-mode coupling is a helioseismic technique that uses measurements of mode eigenfunctions to infer the interior structure of the Sun. This technique has led to insights into the evolution and structure of toroidal flows in the solar interior. Here, we validate an inversion algorithm for normal-mode coupling by generating synthetic seismic measurements associated with input flows and comparing the input and inverted velocities. We study four different cases of input toroidal flows and compute synthetics that take into account the partial visibility of the Sun. We invert the synthetics using Subtractive Optimally Localized Averages (SOLA) and also try to mitigate the systematics of mode leakage. We demonstrate that, ultimately, inversions are only as good as the model we assume for the correlation between flow velocities. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: 14 pages, 20 figures

Solar wind prediction using deep learning

Authors: Vishal Upendran, M. C. M Cheung, Shravan Hanasoge, Ganapathi Krishnamurthi

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

Submitted 10 June, 2020; originally announced June 2020.

Comments: 29 pages, 12 figures, accepted for publication in AGU: Spaceweather

Supervised convolutional neural networks for classification of flaring and nonflaring active regions using line-of-sight magnetograms

Authors: Shamik Bhattacharjee, Rasha Alshehhi, Dattaraj B. Dhuri, Shravan M. Hanasoge

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

Submitted 27 May, 2020; originally announced May 2020.

Comments: 16 pages, 13 figures, accepted for publication in ApJ

Sensitivity kernels for inferring Lorentz stresses from normal-mode frequency splittings in the Sun

Authors: Srijan Bharati Das, Tuneer Chakraborty, Shravan M. Hanasoge, Jeroen Tromp

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

Submitted 29 April, 2020; originally announced April 2020.

Comments: 32 pages, 8 figures, Astrophysical Journal

Validating time-distance helioseismic inversions for non-separable subsurface profiles of an average supergranule

Authors: Vedant Dhruv, Jishnu Bhattacharya, Shravan Hanasoge

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

Submitted 22 August, 2019; originally announced August 2019.

Comments: 20 pages, 7 figures

Properties of solar Rossby waves from normal mode coupling and characterizing its systematics

Authors: Krishnendu Mandal, Shravan Hanasoge

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

Submitted 5 February, 2020; v1 submitted 16 August, 2019; originally announced August 2019.

Comments: Added references, appendix A and errorbar in Figure 4. Accepted for publication in ApJ

Machine learning reveals systematic accumulation of electric current in lead-up to solar flares

Authors: Dattaraj B. Dhuri, Shravan M. Hanasoge, Mark C. M. Cheung

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

Submitted 24 May, 2019; originally announced May 2019.

Comments: Related code is available at https://github.com/dbdhuri/ML-SolarFlare-Analysis

Journal ref: Published in Proceedings of the National Academy of Sciences of the United States of America (PNAS) May, 2019

Latitudinal differential rotation in the solar analogues 16 Cyg A and B

Authors: M. Bazot, O. Benomar, J. Christensen-Dalsgaard, L. Gizon, S. Hanasoge, M. Nielsen, P. Petit, K. R. Sreenivasan

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

Submitted 9 February, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: Accepted in A&A

Detection of Rossby waves in the Sun using normal-mode coupling

Authors: Shravan Hanasoge, Krishnendu Mandal

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

Submitted 19 January, 2019; originally announced January 2019.

Comments: 7 pages, 3 figures, accepted for publication in ApJL

Supervised Neural Networks for Helioseismic Ring-Diagram Inversions

Authors: Rasha Alshehhi, Chris S. Hanson, Laurent Gizon, Shravan Hanasoge

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

Submitted 6 January, 2019; originally announced January 2019.

Comments: 10 pages, 10 Figures, Accepted by A&A

Journal ref: A&A 622, A124 (2019)

Asteroseismic detection of latitudinal differential rotation in 13 Sun-like stars

Authors: O. Benomar, M. Bazot, M. B. Nielsen, L. Gizon, T. Sekii, M. Takata, H. Hotta, S. Hanasoge, K. R. Sreenivasan, J. Christensen-Dalsgaard

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

Submitted 21 September, 2018; originally announced September 2018.

Comments: 45 pages, 11 figures, 4 tables, published in Science

Journal ref: Science, Vol. 361, Issue 6408, pp. 1231-1234, 21 Sep 2018

Helioseismic inversion to infer depth profile of solar meridional flow using spherical Born kernels

Authors: K. Mandal, S. M. Hanasoge, S. P. Rajaguru, H. M. Antia

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

Submitted 1 July, 2018; originally announced July 2018.

Comments: 24 pages, 8 figures, accepted for publication in ApJ

Asymmetry of line profiles of stellar oscillations measured by Kepler for ensembles of solar-like oscillators: impact on mode frequencies and dependence on effective temperature

Authors: O. Benomar, Mjo. Goupil, K. Belkacem, T. Appourchaux, M. B. Nielsen, M. Bazot, L. Gizon, S. Hanasoge, K. R. Sreenivasan, B. Marchand

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

Submitted 17 April, 2018; originally announced April 2018.

Comments: 8 pages, 7 Figures, 1 Table, Accepted to ApJ

Measurement process and inversions using helioseismic normal-mode coupling

Authors: Shravan Hanasoge

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

Submitted 8 May, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: 23 pages, 11 figures, ApJ

Iterative inversion of synthetic travel times successful at recovering sub-surface profiles of supergranular flows

Authors: Jishnu Bhattacharya, Shravan M. Hanasoge, Aaron C. Birch, Laurent Gizon

Abstract: 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… ▽ More 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. △ Less

Submitted 11 August, 2017; originally announced August 2017.

Comments: Accepted for publication in in A&A, 11 pages, 6 figures

Journal ref: A&A 607, A129 (2017)

Seismic sensitivity of Normal-mode Coupling to Lorentz stresses in the Sun

Authors: Shravan M. Hanasoge

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

Submitted 28 July, 2017; v1 submitted 26 May, 2017; originally announced May 2017.

Comments: 18 pages, 4 figures; MNRAS

Sensitivity of Helioseismic Measurements of Normal-mode Coupling to Flows and Sound-speed Perturbations

Authors: Shravan M. Hanasoge, Martin Woodard, H. M. Antia, Laurent Gizon, Katepalli R. Sreenivasan

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

Submitted 23 May, 2017; originally announced May 2017.

Comments: 24 pages, 8 Figures; MNRAS

Finite-frequency sensitivity kernels in spherical geometry for time-distance helioseismology

Authors: Krishnendu Mandal, Jishnu Bhattacharya, Samrat Halder, Shravan Hanasoge

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

Submitted 11 May, 2017; originally announced May 2017.

Comments: 32 pages, 10 figures, accepted for publication in ApJ

Strategies in seismic inference of supergranular flows on the Sun

Authors: Jishnu Bhattacharya, Shravan M. Hanasoge

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

Submitted 30 May, 2016; originally announced May 2016.

Comments: Accepted in ApJ, 12 pages, 10 figures, 1 table

Asteroseismic determination of fundamental parameters of sun-like stars using multi-layered neural networks

Authors: Kuldeep Verma, Shravan Hanasoge, Jishnu Bhattacharya, H M Antia, Ganapathy Krishnamurthi

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

Submitted 1 July, 2016; v1 submitted 2 February, 2016; originally announced February 2016.

Comments: Submitted. 9 pages, 7 Figures, 1 Table, MNRAS

Frequency shifts of resonant modes of the Sun due to near-surface convective scattering

Authors: Jishnu Bhattacharya, Shravan. M. Hanasoge, H. M. Antia

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

Submitted 15 May, 2015; originally announced May 2015.

Comments: 9 pages, 4 figures

Solar Dynamics, Rotation, Convection and Overshoot

Authors: S. Hanasoge, M. S. Miesch, M. Roth, J. Schou, M. Schuessler, M. J. Thompson

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

Submitted 30 March, 2015; originally announced March 2015.

Comments: 31 pages, 10 figures, Space Science Reviews

Seismic Sounding of Convection in the Sun

Authors: Shravan Hanasoge, Laurent Gizon, Katepalli R. Sreenivasan

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

Submitted 27 March, 2015; originally announced March 2015.

Comments: 29 pages, 12 figures, in review, Annual Reviews of Fluid Mechanics

Simulating acoustic waves in spotted stars

Authors: Emanuele Papini, Aaron C. Birch, Laurent Gizon, Shravan M. Hanasoge

Abstract: 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 stu… ▽ More 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. △ Less

Submitted 13 April, 2015; v1 submitted 20 March, 2015; originally announced March 2015.

Comments: 7 pages, 5 figures, accepted for publication in A&A on 19 March 2015

Journal ref: A&A 577, A145 (2015)

Full Waveform Inversion of Solar Interior Flows

Authors: Shravan M Hanasoge

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

Submitted 8 October, 2014; originally announced October 2014.

Comments: 21 pages, 6 figures; in press, ApJ

Additional Evidence Supporting a Model of Shallow, High-Speed Supergranulation

Authors: T. L. Duvall Jr., S. M. Hanasoge, S. Chakraborty

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

Submitted 9 April, 2014; originally announced April 2014.

Full Waveform Inversion for Time-Distance Helioseismology

Authors: Shravan Hanasoge, Jeroen Tromp

Abstract: 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 so… ▽ More 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. △ Less

Submitted 29 January, 2014; originally announced January 2014.

Comments: 24 pages, 10 figures, to appear in ApJ

The Quest to Understand Supergranulation and Large-Scale Convection in the Sun

Authors: Shravan M. Hanasoge, Katepalli R. Sreenivasan

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

Submitted 31 December, 2013; originally announced January 2014.

Comments: 19 pages, 4 figures, Invited Review, Accepted for publication in Solar Physics Topical Edition

Measurements and Kernels for Source-Structure Inversions in Noise Tomography

Authors: Shravan M. Hanasoge

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

Submitted 2 October, 2013; originally announced October 2013.

Comments: 19 pages, 12 figures, Geophysical Journal International

Seismic constraints on rotation of Sun-like star and mass of exoplanet

Authors: Laurent Gizon, Jérome Ballot, Eric Michel, Thorsten Stahn, Gérard Vauclair, Hans Bruntt, Pierre-Olivier Quirion, Othman Benomar, Sylvie Vauclair, Thierry Appourchaux, Michel Auvergne, Annie Baglin, Caroline Barban, Fréderic Baudin, Michaël Bazot, Tiago Campante, Claude Catala, William Chaplin, Orlagh Creevey, Sébastien Deheuvels, Noël Dolez, Yvonne Elsworth, Rafael García, Patrick Gaulme, Stéphane Mathis , et al. (10 additional authors not shown)

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

Submitted 20 August, 2013; originally announced August 2013.

Comments: Published in Proceedings of the National Academy of Sciences (5 pages, 5 figures, 3 tables). Available at http://www.pnas.org/cgi/doi/10.1073/pnas.1303291110

Journal ref: PNAS 2013, vol. 110, no. 33, pp. 13267-13271

Interpreting Cross-correlations of One-bit Filtered Seismic Noise

Authors: Shravan Hanasoge, Michal Branicki

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

Submitted 16 August, 2013; originally announced August 2013.

Comments: 23 pages, 13 figures, accepted for publication in Geophysical Journal International