Showing 1–50 of 78 results for author: Birch, A C

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

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

Submitted 8 June, 2024; originally announced June 2024.

Impact of spatially correlated fluctuations in sunspots on metrics related to magnetic twist

Authors: C. Baumgartner, A. C. Birch, H. Schunker, R. H. Cameron, L. Gizon

Abstract: The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxi… ▽ More The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxies is based on $αあるふぁ_z$, the ratio of the vertical current to the vertical magnetic field. Another class of proxies is based on the so-called twist density, $q$, which depends on the ratio of the azimuthal field to the vertical field. However, the sensitivity of these proxies to temporal fluctuations of the magnetic field has not yet been well characterized. We aim to determine the sensitivity of twist proxies to temporal fluctuations in the magnetic field as estimated from time-series of SDO/HMI vector magnetic field maps. To this end, we introduce a model of a sunspot with a peak vertical field of 2370 Gauss at the photosphere and a uniform twist density $q= -0.024$ Mm$^{-1}$. We add realizations of the temporal fluctuations of the magnetic field that are consistent with SDO/HMI observations, including the spatial correlations. Using a Monte-Carlo approach, we determine the robustness of the different proxies to the temporal fluctuations. The temporal fluctuations of the three components of the magnetic field are correlated for spatial separations up to 1.4 Mm (more than expected from the point spread function alone). The Monte-Carlo approach enables us to demonstrate that several proxies for the twist of the magnetic field are not biased in each of the individual magnetograms. The associated random errors on the proxies have standard deviations in the range between $0.002$ and $0.006$ Mm$^{-1}$, which is smaller by approximately one order of magnitude than the mean value of $q$. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: accepted in A&A

Journal ref: A&A 664, A183 (2022)

Testing solar surface flux transport models in the first days after active region emergence

Authors: Nils Gottschling, Hannah Schunker, Aaron C. Birch, Robert Cameron, Laurent Gizon

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

Submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Astronomy & Astrophysics. 13 pages, 11 figures

Journal ref: A&A 660, A6 (2022)

Solar inertial modes: Observations, identification, and diagnostic promise

Authors: Laurent Gizon, Robert H. Cameron, Yuto Bekki, Aaron C. Birch, Richard S. Bogart, Allan Sacha Brun, Cilia Damiani, Damien Fournier, Laura Hyest, Kiran Jain, B. Lekshmi, Zhi-Chao Liang, Bastian Proxauf

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

Submitted 20 July, 2021; originally announced July 2021.

Comments: Forthcoming article, Astronomy and Astrophysics

Evolution of solar surface inflows around emerging active regions

Authors: N. Gottschling, H. Schunker, A. C. Birch, B. Löptien, L. Gizon

Abstract: Solar active regions are associated with Evershed outflows in sunspot penumbrae, moat outflows surrounding sunspots, and extended inflows surrounding active regions. The latter have been identified on established active regions by various methods. The evolution of these inflows and their dependence on active region properties as well as their impact on the global magnetic field are not yet underst… ▽ More Solar active regions are associated with Evershed outflows in sunspot penumbrae, moat outflows surrounding sunspots, and extended inflows surrounding active regions. The latter have been identified on established active regions by various methods. The evolution of these inflows and their dependence on active region properties as well as their impact on the global magnetic field are not yet understood. We aim to understand the evolution of the average inflows around emerging active regions and to derive an empirical model for these inflows. We analyze horizontal flows at the surface of the Sun using local correlation tracking of solar granules observed in continuum images of SDO/HMI. We measure average flows of a sample of 182 isolated active regions up to seven days before and after their emergence onto the solar surface with a cadence of 12 hours. We investigate the average inflow properties with respect to active region characteristics of total flux and latitude. We fit a model to these observed inflows for a quantitative analysis. We find that converging flows of around $20$ to $30$ m/s are first visible one day prior to emergence, in agreement with recent results. These converging flows are present independently of active region properties of latitude or flux. We confirm a recently found prograde flow of about $40$ m/s at the leading polarity during emergence. We find that the time after emergence when the latitudinal inflows increase in amplitude depends on the flux of the active region, ranging from one to four days after emergence and increasing with flux. The largest extent of the inflows is up to about $7 \pm 1^\circ$ away from the center of the active region within the first six days after emergence. The inflow velocities have amplitudes of about $50$ m/s. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: Accepted for publication in Astronomy & Astrophysics. 16 pages, 17 figures

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

A journey of exploration to the polar regions of a star: probing the solar poles and the heliosphere from high helio-latitude

Authors: Louise Harra, Vincenzo Andretta, Thierry Appourchaux, Frédéric Baudin, Luis Bellot-Rubio, Aaron C. Birch, Patrick Boumier, Robert H. Cameron, Matts Carlsson, Thierry Corbard, Jackie Davies, Andrew Fazakerley, Silvano Fineschi, Wolfgang Finsterle, Laurent Gizon, Richard Harrison, Donald M. Hassler, John Leibacher, Paulett Liewer, Malcolm MacDonald, Milan Maksimovic, Neil Murphy, Giampiero Naletto, Giuseppina Nigro, Christopher Owen , et al. (7 additional authors not shown)

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

Submitted 22 April, 2021; originally announced April 2021.

Comments: 29 pages, 9 figures, White paper for ESA Voyage 2050

Probing the Solar Meridional Circulation using Fourier Legendre Decomposition

Authors: D. C. Braun, A. C. Birch, Y. Fan

Abstract: We apply the helioseismic methodology of Legendre Function Decomposition to 88 months of Dopplergrams obtained by the Helioseismic and Magnetic Imager (HMI) as the basis of inferring the depth variation of the mean meridional flow, as averaged between 20 and 60 degrees latitude and in time, in both the northern and southern hemispheres. We develop and apply control procedures designed to assess an… ▽ More We apply the helioseismic methodology of Legendre Function Decomposition to 88 months of Dopplergrams obtained by the Helioseismic and Magnetic Imager (HMI) as the basis of inferring the depth variation of the mean meridional flow, as averaged between 20 and 60 degrees latitude and in time, in both the northern and southern hemispheres. We develop and apply control procedures designed to assess and remove center-to-limb artifacts, using measurements obtained by performing the analysis with respect to artificial poles at the east and west limbs. Forward modeling is carried out, using sensitivity functions proportional to the mode kinetic energy density, to evaluate the consistency of the corrected frequency shifts with models of the depth variation of the meridional circulation in the top half of the convection zone. The results, taken at face value, imply substantial differences between the meridional circulation in the northern and southern hemisphere. The inferred presence of a return (equator-ward propagating) flow at a depth of approximately 40 Mm below the photosphere in the northern hemisphere is surprising and appears to be inconsistent with many other helioseismic analyses. This discrepancy may be the result of an inadequacy of our methodology to remove systematic errors in HMI data. Our results appear to be at least qualitative similar to those by Gizon et al. (2020) which point to an anomaly in HMI data that is not present in MDI or GONG data. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: Accepted for publication in the Astrophysical Journal. 20 pages, 16 figures, 3 appendices

Helioseismological determination of the subsurface spatial spectrum of solar convection: Demonstration using numerical simulations

Authors: Vincent G. A. Böning, Aaron C. Birch, Laurent Gizon, Thomas L. Duvall Jr

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

Submitted 17 February, 2021; originally announced February 2021.

Comments: 16 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 649, A59 (2021)

Solar east-west flow correlations that persist for months at low latitudes are dominated by active region inflows

Authors: Chris S. Hanson, Thomas L. Duvall Jr., Aaron C. Birch, Laurent Gizon, Katepalli R. Sreenivasan

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

Submitted 25 October, 2020; originally announced October 2020.

Comments: 9 pages, 10 figures, accepted into A&A

Journal ref: A&A 644, A103 (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

Predicting frequency changes of global-scale solar Rossby modes due to solar cycle changes in internal rotation

Authors: C. R. Goddard, A. C. Birch, D. Fournier, L. Gizon

Abstract: Context. Large-scale equatorial Rossby modes have been observed on the Sun over the last two solar cycles. Aims. We investigate the impact of the time-varying zonal flows on the frequencies of Rossby modes. Methods. A first-order perturbation theory approach is used to obtain an expression for the expected shift in the mode frequencies due to perturbations in the internal rotation rate. Results. U… ▽ More Context. Large-scale equatorial Rossby modes have been observed on the Sun over the last two solar cycles. Aims. We investigate the impact of the time-varying zonal flows on the frequencies of Rossby modes. Methods. A first-order perturbation theory approach is used to obtain an expression for the expected shift in the mode frequencies due to perturbations in the internal rotation rate. Results. Using the time-varying rotation from helioseismic inversions we predict the changes in Rossby mode frequencies with azimuthal orders from m = 1 to m = 15 over the last two solar cycles. The peak-to-peak frequency change is less than 1 nHz for the m = 1 mode, grows with m, and reaches 25 nHz for m = 15. Conclusions. Given the observational uncertainties on mode frequencies due to the finite mode lifetimes, we find that the predicted frequency shifts are near the limit of detectability. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Journal ref: A&A 640, L10 (2020)

Average motion of emerging solar active region polarities II: Joy's law

Authors: Hannah Schunker, Christian Baumgartner, Aaron C. Birch, Robert H. Cameron, Douglas C. Braun, Laurent Gizon

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

Submitted 9 June, 2020; originally announced June 2020.

Comments: accepted in Astronomy and Astrophysics

Journal ref: A&A 640, A116 (2020)

Estimating the nonstructural component of the helioseismic surface term using hydrodynamic simulations

Authors: J. Schou, A. C. Birch

Abstract: As the amount of asteroseismic data available continues to grow, the inability to accurately model observed oscillation frequencies is becoming a critical problem for interpreting these frequencies. A major component of this problem is the modeling of the near-surface layers. Our aim is to develop a method to estimate the effect of the near-surface layers on oscillation frequencies. In the propose… ▽ More As the amount of asteroseismic data available continues to grow, the inability to accurately model observed oscillation frequencies is becoming a critical problem for interpreting these frequencies. A major component of this problem is the modeling of the near-surface layers. Our aim is to develop a method to estimate the effect of the near-surface layers on oscillation frequencies. In the proposed method we numerically estimate eigenfunctions in 3D hydrodynamic simulations. We match those to the eigenfunctions calculated from the classic equations applied to the horizontal averages of the structure variables. We use this procedure to calculate the frequency perturbation resulting from the dynamical part of the interaction of the oscillations with near-surface convection. As the last step we scale the numbers to the Sun. To provide a qualitative test of our method we performed a series of simulations, calculated the perturbations using our procedure, and compared them to previously reported residuals relative to solar models. We find that we can largely reproduce the observed frequency residuals without resorting to poorly justified theoretical models. We find that, while the calculations of Houdek et al. (2017, MNRAS, 464, L124) produce similar frequency perturbations, the density-pressure phase differences computed here do not match those of that work. △ Less

Submitted 28 April, 2020; originally announced April 2020.

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

Journal ref: A&A 638, A51 (2020)

Rossby modes in slowly rotating stars: depth dependence in distorted polytropes with uniform rotation

Authors: C. Damiani, R. H. Cameron, A. C. Birch, L. Gizon

Abstract: Large-scale Rossby waves have recently been discovered from measurements of horizontal surface and near-surface solar flows (Löptien at al. 2018). We are interested in understanding why only the sectoral modes are seen in the observations and also in modelling the radial structure of the observed modes. To do so, we characterise here the radial eigenfunctions of r modes for slowly-rotating polytro… ▽ More Large-scale Rossby waves have recently been discovered from measurements of horizontal surface and near-surface solar flows (Löptien at al. 2018). We are interested in understanding why only the sectoral modes are seen in the observations and also in modelling the radial structure of the observed modes. To do so, we characterise here the radial eigenfunctions of r modes for slowly-rotating polytropes in uniform rotation. We find that for free-surface boundary conditions on a spheroid of non-vanishing surface density, r modes can only exist for $\ell=m$ spherical harmonics in the inviscid case, and we compute their depth dependence and frequencies to leading order. For quasi-adiabatic stratification the sectoral modes with no radial nodes are the only modes which are almost toroidal and the depth dependence of the corresponding horizontal motion scales as $r^m$. For all r modes except the zero radial order sectoral ones, non-adiabatic stratification plays a crucial role in the radial force balance. The lack of quasi-toroidal solutions when stratification is close to neutral, except for the sectoral modes without nodes in radius, follows from the statement that the system needs to be in both horizontal and radial force balance. In the absence of super- or subadiabatic stratification and viscosity, both the horizontal and radial force balances independently determine the pressure perturbation. The only quasi-toroidal cases in which the two determinations of the pressure perturbation are consistent are the special cases where $\ell=m$, and the horizontal displacement scales with $r^m$. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: 10 pages, 8 figures, accepted in A&A

Journal ref: A&A 637, A65 (2020)

Characterizing the spatial pattern of solar supergranulation using the bispectrum

Authors: Vincent G. A. Böning, Aaron C. Birch, Laurent Gizon, Thomas L. Duvall Jr., Jesper Schou

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

Submitted 19 February, 2020; originally announced February 2020.

Comments: 16 pages, 12 figures, accepted for publication by A&A

Journal ref: A&A 635, A181 (2020)

Exploring the latitude and depth dependence of solar Rossby waves using ring-diagram analysis

Authors: B. Proxauf, L. Gizon, B. Löptien, J. Schou, A. C. Birch, R. S. Bogart

Abstract: Global-scale Rossby waves have recently been unambiguously identified on the Sun. Here we study the latitude and depth dependence of the Rossby wave eigenfunctions. By applying helioseismic ring-diagram analysis and granulation tracking to SDO/HMI observations, we compute maps of the radial vorticity of flows in the upper solar convection zone (down to depths of more than $16$ Mm). We use a Fourie… ▽ More Global-scale Rossby waves have recently been unambiguously identified on the Sun. Here we study the latitude and depth dependence of the Rossby wave eigenfunctions. By applying helioseismic ring-diagram analysis and granulation tracking to SDO/HMI observations, we compute maps of the radial vorticity of flows in the upper solar convection zone (down to depths of more than $16$ Mm). We use a Fourier transform in longitude to separate the different azimuthal orders m in the range $3 \le m \le 15$. At each $m$ we obtain the phase and amplitude of the Rossby waves as a function of depth using the helioseismic data. At each $m$ we also measure the latitude dependence of the eigenfunctions by calculating the covariance between the equator and other latitudes. We then study the horizontal and radial dependences of the radial vorticity eigenfunctions. The horizontal eigenfunctions are complex. As observed previously, the real part peaks at the equator and switches sign near $\pm 30^\circ$, thus the eigenfunctions show significant non-sectoral contributions. The imaginary part is smaller than the real part. The phase of the radial eigenfunctions varies by only roughly $\pm 5^\circ$ over the top $15$ Mm. The amplitude of the radial eigenfunctions decreases by about $10\%$ from the surface down to $8$ Mm (the region where ring-diagram analysis is most reliable, as seen by comparing with the rotation rate measured by global-mode seismology). The radial dependence of the radial vorticity eigenfunctions deduced from ring-diagram analysis is consistent with a power-law down to $8$ Mm and is unreliable at larger depths. However, the observations provide only weak constraints on the power-law exponents. For the real part, the latitude dependence of the eigenfunctions is consistent with previous work (using granulation tracking). The imaginary part is smaller than the real part but significantly nonzero. △ Less

Submitted 8 January, 2020; v1 submitted 4 December, 2019; originally announced December 2019.

Comments: Submitted to and accepted by Astronomy and Astrophysics; includes minor changes from language editing report

Journal ref: A&A 634, A44 (2020)

An improved multi-ridge fitting method for ring-diagram helioseismic analysis

Authors: Kaori Nagashima, Aaron C. Birch, Jesper Schou, Bradley W. Hindman, Laurent Gizon

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

Submitted 18 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 633, A109 (2020)

Average motion of emerging solar active region polarities I: Two phases of emergence

Authors: Hannah Schunker, Aaron C. Birch, Robert H. Cameron, Douglas C. Braun, Laurent Gizon, Raymond B. Burston

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

Submitted 30 April, 2019; v1 submitted 28 March, 2019; originally announced March 2019.

Comments: accepted A&A

Journal ref: A&A 625, A53 (2019)

Time-distance helioseismology of solar Rossby waves

Authors: Zhi-Chao Liang, Laurent Gizon, Aaron C. Birch, Thomas L. Duvall Jr

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

Submitted 29 April, 2019; v1 submitted 18 December, 2018; originally announced December 2018.

Comments: accepted for publication in A&A

Journal ref: A&A 626, A3 (2019)

Signal and noise in helioseismic holography

Authors: Laurent Gizon, Damien Fournier, Dan Yang, Aaron C. Birch, Hélène Barucq

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

Submitted 30 September, 2018; originally announced October 2018.

Comments: Submitted to Astronomy and Astrophysics

Journal ref: A&A 620, A136 (2018)

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

Authors: Zhi-Chao Liang, Laurent Gizon, Aaron C. Birch, Thomas L. Duvall Jr., S. P. Rajaguru

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

Submitted 17 October, 2018; v1 submitted 27 August, 2018; originally announced August 2018.

Comments: accepted for publication in A&A

Journal ref: A&A 619, A99 (2018)

Validating Forward Modeling and Inversions of Helioseismic Holography Measurements

Authors: K. DeGrave, D. C. Braun, A. C. Birch, A. D. Crouch, B. Javornik

Abstract: Here we use synthetic data to explore the performance of forward models and inverse methods for helioseismic holography. Specifically, this work presents the first comprehensive test of inverse modeling for flows using lateral-vantage (deep-focus) holography. We derive sensitivity functions in the Born approximation. We then use these sensitivity functions in a series of forward models and inversi… ▽ More Here we use synthetic data to explore the performance of forward models and inverse methods for helioseismic holography. Specifically, this work presents the first comprehensive test of inverse modeling for flows using lateral-vantage (deep-focus) holography. We derive sensitivity functions in the Born approximation. We then use these sensitivity functions in a series of forward models and inversions of flows from a publicly available magnetohydrodynamic quiet-Sun simulation. The forward travel times computed using the kernels generally compare favorably with measurements obtained by applying holography, in a lateral-vantage configuration, on a 15-hour time series of artificial Dopplergrams extracted from the simulation. Inversions for the horizontal flow components are able to reproduce the flows in the upper 3Mm of the domain, but are compromised by noise at greater depths. △ Less

Submitted 10 July, 2018; originally announced July 2018.

Comments: accepted for publication by the Astrophysical J

Evolution and wave-like properties of the average solar supergranule

Authors: J. Langfellner, A. C. Birch, L. Gizon

Abstract: Solar supergranulation presents us with many mysteries. For example, previous studies in spectral space found that supergranulation has wave-like properties. Here we study, in real space, the wave-like evolution of the average supergranule over a range of spatial scales (from 10 to 80 Mm). We complement this by characterizing the evolution of the associated network magnetic field. We use one year… ▽ More Solar supergranulation presents us with many mysteries. For example, previous studies in spectral space found that supergranulation has wave-like properties. Here we study, in real space, the wave-like evolution of the average supergranule over a range of spatial scales (from 10 to 80 Mm). We complement this by characterizing the evolution of the associated network magnetic field. We use one year of data from the Helioseismic and Magnetic Imager (HMI) to measure horizontal near-surface flows near the solar equator by applying time-distance helioseismology on Dopplergrams and granulation tracking on intensity images. The average supergranule outflow (or inflow) is constructed by averaging over 10000 individual outflows (or inflows). The contemporaneous evolution of the magnetic field is studied with HMI line-of-sight observations. We confirm and extend previous measurements of the supergranular wave dispersion relation to angular wavenumbers in the range 50<kR<270. We find a plateau for kR>120. In real space, larger supergranules undergo oscillations with longer periods and lifetimes than smaller cells. We find excellent agreement between TD and LCT and obtain wave properties that are independent of the tracking rate. The observed network magnetic field follows the oscillations of the supergranular flows with a six-hour time lag. This behavior can be explained by computing the motions of corks carried by the supergranular flows. Signatures of supergranular waves in surface horizontal flows near the solar equator can be observed in real space. These oscillatory flows control the evolution of the network magnetic field, in particular they explain the recently discovered east-west anisotropy of the magnetic field around the average supergranule. Background flow measurements that we obtain from Doppler frequency shifts do not favor shallow models of supergranulation. △ Less

Submitted 31 May, 2018; originally announced May 2018.

Comments: 13 pages, 14 figures (plus appendix), accepted for publication in A&A

Journal ref: A&A 617, A97 (2018)

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

Authors: Björn Löptien, Laurent Gizon, Aaron C. Birch, Jesper Schou, Bastian Proxauf, Thomas L. Duvall Jr., Richard S. Bogart, Ulrich R. Christensen

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

Submitted 18 May, 2018; originally announced May 2018.

Comments: This is the submitted version of the paper published in Nature Astronomy. 23 pages, 8 figures, 1 table

Journal ref: Nature Astronomy (2018)

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)

Helioseismic Holography of Simulated Sunspots: dependence of the travel time on magnetic field strength and Wilson depression

Authors: T. Felipe, D. C. Braun, A. C. Birch

Abstract: Improving methods for determining the subsurface structure of sunspots from their seismic signature requires a better understanding of the interaction of waves with magnetic field concentrations. We aim to quantify the impact of changes in the internal structure of sunspots on local helioseismic signals. We have numerically simulated the propagation of a stochastic wave field through sunspot model… ▽ More Improving methods for determining the subsurface structure of sunspots from their seismic signature requires a better understanding of the interaction of waves with magnetic field concentrations. We aim to quantify the impact of changes in the internal structure of sunspots on local helioseismic signals. We have numerically simulated the propagation of a stochastic wave field through sunspot models with different properties, accounting for changes in the Wilson depression between 250 and 550 km and in the photospheric umbral magnetic field between 1500 and 3500 G. The results show that travel-time shifts at frequencies above approximately 3.50 mHz (depending on the phase-speed filter) are insensitive to the magnetic field strength. The travel time of these waves is determined exclusively by the Wilson depression and sound-speed perturbation. The travel time of waves with lower frequencies is affected by the direct effect of the magnetic field, although photospheric field strengths below 1500 G do not leave a significant trace on the travel-time measurements. These results could potentially be used to develop simplified travel-time inversion methods. △ Less

Submitted 25 May, 2017; originally announced May 2017.

Comments: Accepted for publication in A&A

Journal ref: A&A 604, A126 (2017)

Measuring solar active region inflows with local correlation tracking of granulation

Authors: B. Löptien, A. C. Birch, T. L. Duvall Jr., L. Gizon, B. Proxauf, J. Schou

Abstract: Context. Local helioseismology has detected spatially extended converging surface flows into solar active regions. These play an important role in flux-transport models of the solar dynamo. Aims. We aim to validate the existence of the inflows by deriving horizontal flow velocities around active regions with local correlation tracking of granulation. Methods. We generate a six-year long-time s… ▽ More Context. Local helioseismology has detected spatially extended converging surface flows into solar active regions. These play an important role in flux-transport models of the solar dynamo. Aims. We aim to validate the existence of the inflows by deriving horizontal flow velocities around active regions with local correlation tracking of granulation. Methods. We generate a six-year long-time series of full-disk maps of the horizontal velocity at the solar surface by tracking granules in continuum intensity images provided by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Results. On average, active regions are surrounded by inflows extending up to 10 deg from the center of the active region of magnitudes of 20-30 m/s, reaching locally up to 40 m/s, which is in agreement with results from local helioseismology. By computing an ensemble average consisting of 243 individual active regions, we show that the inflows are not azimuthally symmetric but converge predominantly towards the trailing polarity of the active region with respect to the longitudinally and temporally averaged flow field. △ Less

Submitted 24 May, 2017; originally announced May 2017.

Comments: 10 pages, 9 figures. Submitted to Astronomy & Astrophysics

Journal ref: A&A 606, A28 (2017)

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

Authors: Zhi-Chao Liang, Aaron C. Birch, Thomas L. Duvall Jr., Laurent Gizon, Jesper Schou

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

Submitted 3 April, 2017; originally announced April 2017.

Comments: accepted for publication in A&A

Journal ref: A&A 601, A46 (2017)

The amplitude of the cross-covariance function of solar oscillations as a diagnostic tool for wave attenuation and geometrical spreading

Authors: Kaori Nagashima, Damien Fournier, Aaron C. Birch, Laurent Gizon

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

Submitted 28 December, 2016; originally announced December 2016.

Comments: accepted for publication in A&A

Journal ref: A&A 599, A111 (2017)

Shape of a slowly rotating star measured by asteroseismology

Authors: Laurent Gizon, Takashi Sekii, Masao Takata, Donald W. Kurtz, Hiromoto Shibahashi, Michael Bazot, Othman Benomar, Aaron C. Birch, Katepalli R. Sreenivasan

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

Submitted 19 November, 2016; originally announced November 2016.

Comments: 15 pages, 2 tables, 3 figures

Journal ref: Science Advances 16 Nov 2016, Vol. 2, no. 11, e1601777

Computational helioseismology in the frequency domain: acoustic waves in axisymmetric solar models with flows

Authors: Laurent Gizon, Hélène Barucq, Marc Duruflé, Chris S. Hanson, Michael Leguèbe, Aaron C. Birch, Juliette Chabassier, Damien Fournier, Thorsten Hohage, Emanuele Papini

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

Submitted 5 November, 2016; originally announced November 2016.

Comments: Accepted on 3 Nov 2016 for publication in Astronomy and Astrophysics

Journal ref: A&A 600, A35 (2017)

Intensity contrast of the average supergranule

Authors: J. Langfellner, A. C. Birch, L. Gizon

Abstract: While the velocity fluctuations of supergranulation dominate the spectrum of solar convection at the solar surface, very little is known about the fluctuations in other physical quantities like temperature or density at supergranulation scale. Using SDO/HMI observations, we characterize the intensity contrast of solar supergranulation at the solar surface. We identify the positions of… ▽ More While the velocity fluctuations of supergranulation dominate the spectrum of solar convection at the solar surface, very little is known about the fluctuations in other physical quantities like temperature or density at supergranulation scale. Using SDO/HMI observations, we characterize the intensity contrast of solar supergranulation at the solar surface. We identify the positions of ${\sim}10^4$ outflow and inflow regions at supergranulation scales, from which we construct average flow maps and co-aligned intensity and magnetic field maps. In the average outflow center, the maximum intensity contrast is $(7.8\pm0.6)\times10^{-4}$ (there is no corresponding feature in the line-of-sight magnetic field). This corresponds to a temperature perturbation of about $1.1\pm0.1$ K, in agreement with previous studies. We discover an east-west anisotropy, with a slightly deeper intensity minimum east of the outflow center. The evolution is asymmetric in time: the intensity excess is larger 8 hours before the reference time (the time of maximum outflow), while it has almost disappeared 8 hours after the reference time. In the average inflow region, the intensity contrast mostly follows the magnetic field distribution, except for an east-west anisotropic component that dominates 8 hours before the reference time. We suggest that the east-west anisotropy in the intensity is related to the wave-like properties of supergranulation. △ Less

Submitted 29 September, 2016; originally announced September 2016.

Comments: 5 pages, 3 figures (plus appendix), accepted for publication in A&A

Journal ref: A&A 596, A66 (2016)

SDO/HMI survey of emerging active regions for helioseismology

Authors: H. Schunker, D. C. Braun, A. C. Birch, R. B. Burston, L. Gizon

Abstract: Observations from the Solar Dynamics Observatory (SDO) have the potential for allowing the helioseismic study of the formation of hundreds of active regions, which would enable us to perform statistical analyses. Our goal is to collate a uniform data set of emerging active regions observed by the SDO/HMI instrument suitable for helioseismic analysis up to seven days before emergence. We restricted… ▽ More Observations from the Solar Dynamics Observatory (SDO) have the potential for allowing the helioseismic study of the formation of hundreds of active regions, which would enable us to perform statistical analyses. Our goal is to collate a uniform data set of emerging active regions observed by the SDO/HMI instrument suitable for helioseismic analysis up to seven days before emergence. We restricted the sample to active regions that were visible in the continuum and emerged into quiet Sun largely avoiding pre-existing magnetic regions. As a reference data set we paired a control region (CR), with the same latitude and distance from central meridian, with each emerging active region (EAR). We call this data set, which is currently comprised of 105 emerging active regions observed between May 2010 and November 2012, the SDO Helioseismic Emerging Active Region (SDO/HEAR) survey. To demonstrate the utility of a data set of a large number of emerging active regions, we measure the relative east-west velocity of the leading and trailing polarities from the line-of-sight magnetogram maps during the first day after emergence. The latitudinally averaged line-of-sight magnetic field of all the EARs shows that, on average, the leading (trailing) polarity moves in a prograde (retrograde) direction with a speed of 121 +/- 22 m/s (-70 +/- 13 m/s) relative to the Carrington rotation rate in the first day. However, relative to the differential rotation of the surface plasma, the east-west velocity is symmetric, with a mean of 95 +/- 13 m/s. The SDO/HEAR data set will not only be useful for helioseismic studies, but will also be useful to study other features such as the surface magnetic field evolution of a large sample of EARs. △ Less

Submitted 7 September, 2016; v1 submitted 29 August, 2016; originally announced August 2016.

Comments: Accepted by Astronomy and Astrophysics, 11 figures, one longtable; update corrects units in Figure 1

Helioseismic holography of simulated sunspots: magnetic and thermal contributions to travel times

Authors: T. Felipe, D. C. Braun, A. D. Crouch, A. C. Birch

Abstract: Wave propagation through sunspots involves conversion between waves of acoustic and magnetic character. In addition, the thermal structure of sunspots is very different than that of the quiet Sun. As a consequence, the interpretation of local helioseismic measurements of sunspots has long been a challenge. With the aim of understanding these measurements, we carry out numerical simulations of wave… ▽ More Wave propagation through sunspots involves conversion between waves of acoustic and magnetic character. In addition, the thermal structure of sunspots is very different than that of the quiet Sun. As a consequence, the interpretation of local helioseismic measurements of sunspots has long been a challenge. With the aim of understanding these measurements, we carry out numerical simulations of wave propagation through sunspots. Helioseismic holography measurements made from the resulting simulated wavefields show qualitative agreement with observations of real sunspots. We use additional numerical experiments to determine, separately, the influence of the thermal structure of the sunspot and the direct effect of the sunspot magnetic field. We use the ray approximation to show that the travel-time shifts in the thermal (non-magnetic) sunspot model are primarily produced by changes in the wave path due to the Wilson depression rather than variations in the wave speed. This shows that inversions for the subsurface structure of sunspots must account for local changes in the density. In some ranges of horizontal phase speed and frequency there is agreement (within the noise level in the simulations) between the travel times measured in the full magnetic sunspot model and the thermal model. If this conclusion proves to be robust for a wide range of models, it would suggest a path towards inversions for sunspot structure. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: Accepted for publication in The Astrophysical Journal

A low upper limit on the subsurface rise speed of solar active regions

Authors: Aaron C. Birch, Hannah Schunker, Douglas C. Braun, Robert Cameron, Laurent Gizon, Björn Löptien, Matthias Rempel

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

Submitted 18 July, 2016; originally announced July 2016.

Comments: Published in Science Advances. DOI: 10.1126/sciadv.1600557

Journal ref: Sci. Adv. 2, e1600557 (2016)

Statistics of the two-point cross-covariance function of solar oscillations

Authors: Kaori Nagashima, Takashi Sekii, Laurent Gizon, Aaron C. Birch

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

Submitted 21 June, 2016; originally announced June 2016.

Comments: accepted for publication in A&A

Journal ref: A&A 593, A41 (2016)

The shrinking Sun: a systematic error in local correlation tracking of solar granulation

Authors: B. Löptien, A. C. Birch, T. L. Duvall Jr., L. Gizon, J. Schou

Abstract: Context. Local correlation tracking of granulation (LCT) is an important method for measuring horizontal flows in the photosphere. This method exhibits a systematic error that looks like a flow converging towards disk center, also known as the shrinking-Sun effect. Aims. We aim at studying the nature of the shrinking-Sun effect for continuum intensity data and at deriving a simple model that can… ▽ More Context. Local correlation tracking of granulation (LCT) is an important method for measuring horizontal flows in the photosphere. This method exhibits a systematic error that looks like a flow converging towards disk center, also known as the shrinking-Sun effect. Aims. We aim at studying the nature of the shrinking-Sun effect for continuum intensity data and at deriving a simple model that can explain its origin. Methods. We derived LCT flow maps by running the local correlation tracking code FLCT on tracked and remapped continuum intensity maps provided by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory. We also computed flow maps from synthetic continuum images generated from STAGGER code simulations of solar surface convection. We investigated the origin of the shrinking-Sun effect by generating an average granule from synthetic data from the simulations. Results. The LCT flow maps derived from HMI and from the simulations exhibit a shrinking-Sun effect of comparable magnitude. The origin of this effect is related to the apparent asymmetry of granulation originating from radiative transfer effects when observing with a viewing angle inclined from vertical. This causes, in combination with the expansion of the granules, an apparent motion towards disk center. △ Less

Submitted 19 April, 2016; v1 submitted 15 April, 2016; originally announced April 2016.

Comments: 8 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 590, A130 (2016)

Sensitivity Kernels for Flows in Time-Distance Helioseismology: Extension to Spherical Geometry

Authors: Vincent G. A. Böning, Markus Roth, Wolfgang Zima, Aaron C. Birch, Laurent Gizon

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

Submitted 14 April, 2016; v1 submitted 13 April, 2016; originally announced April 2016.

Comments: 17 pages, 5 figures, 2 tables, accepted for publication in ApJ. v2: typo in arXiv author list corrected

Data compression for local correlation tracking of solar granulation

Authors: Björn Löptien, Aaron C. Birch, Tom L. Duvall Jr., Laurent Gizon, Jesper Schou

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

Submitted 10 December, 2015; originally announced December 2015.

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

Journal ref: A&A 587, A9 (2016)

Anisotropy of the solar network magnetic field around the average supergranule

Authors: J. Langfellner, L. Gizon, A. C. Birch

Abstract: Supergranules in the quiet Sun are outlined by a web-like structure of enhanced magnetic field strength, the so-called magnetic network. We aim to map the magnetic network field around the average supergranule near disk center. We use observations of the line-of-sight component of the magnetic field from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). The a… ▽ More Supergranules in the quiet Sun are outlined by a web-like structure of enhanced magnetic field strength, the so-called magnetic network. We aim to map the magnetic network field around the average supergranule near disk center. We use observations of the line-of-sight component of the magnetic field from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). The average supergranule is constructed by coaligning and averaging over 3000 individual supergranules. We determine the positions of the supergranules with an image segmentation algorithm that we apply on maps of the horizontal flow divergence measured using time-distance helioseismology. In the center of the average supergranule the magnetic (intranetwork) field is weaker by about 2.2 Gauss than the background value (3.5 Gauss), whereas it is enhanced in the surrounding ring of horizontal inflows (by about 0.6 Gauss on average). We find that this network field is significantly stronger west (prograde) of the average supergranule than in the east (by about 0.3 Gauss). With time-distance helioseismology, we find a similar anisotropy. The observed anisotropy of the magnetic field adds to the mysterious dynamical properties of solar supergranulation. △ Less

Submitted 11 June, 2015; v1 submitted 6 May, 2015; originally announced May 2015.

Comments: 4 pages, 4 figures, accepted for publication in A&A Letters

Journal ref: A&A 579, L7 (2015)

Spatially resolved vertical vorticity in solar supergranulation using helioseismology and local correlation tracking

Authors: J. Langfellner, L. Gizon, A. C. Birch

Abstract: Flow vorticity is a fundamental property of turbulent convection in rotating systems. Solar supergranules exhibit a preferred sense of rotation, which depends on the hemisphere. This is due to the Coriolis force acting on the diverging horizontal flows. We aim to spatially resolve the vertical flow vorticity of the average supergranule at different latitudes, both for outflow and inflow regions. T… ▽ More Flow vorticity is a fundamental property of turbulent convection in rotating systems. Solar supergranules exhibit a preferred sense of rotation, which depends on the hemisphere. This is due to the Coriolis force acting on the diverging horizontal flows. We aim to spatially resolve the vertical flow vorticity of the average supergranule at different latitudes, both for outflow and inflow regions. To measure the vertical vorticity, we use two independent techniques: time-distance helioseismology (TD) and local correlation tracking of granules in intensity images (LCT) using data from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Both maps are corrected for center-to-limb systematic errors. We find that 8-h TD and LCT maps of vertical vorticity are highly correlated at large spatial scales. Associated with the average supergranule outflow, we find tangential (vortical) flows that reach about 10 m/s in the clockwise direction at 40° latitude. In average inflow regions, the tangential flow reaches the same magnitude, but in the anti-clockwise direction. These tangential velocities are much smaller than the radial (diverging) flow component (300 m/s for the average outflow and 200 m/s for the average inflow). The results for TD and LCT as measured from HMI are in excellent agreement for latitudes between $-$60° and 60°. From HMI LCT, we measure the vorticity peak of the average supergranule to have a full width at half maximum of about 13 Mm for outflows and 8 Mm for inflows. This is larger than the spatial resolution of the LCT measurements (about 3 Mm). On the other hand, the vorticity peak in outflows is about half the value measured at inflows (e.g. 4/(10^6 s) clockwise compared to 8/(10^6 s) anti-clockwise at 40° latitude). Results from MDI/SOHO obtained in 2010 are biased compared to the HMI/SDO results for the same period. △ Less

Submitted 19 June, 2015; v1 submitted 1 April, 2015; originally announced April 2015.

Comments: 12 pages, 13 figures (plus appendix), accepted for publication in A&A

Journal ref: A&A 581, A67 (2015)

Interpretation of Helioseismic Travel Times - Sensitivity to Sound Speed, Pressure, Density, and Flows

Authors: Raymond Burston, Laurent Gizon, Aaron C. Birch

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

Submitted 11 August, 2015; v1 submitted 31 March, 2015; originally announced March 2015.

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)

Image compression in local helioseismology

Authors: Björn Löptien, Aaron C. Birch, Laurent Gizon, Jesper Schou

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

Submitted 1 December, 2014; v1 submitted 15 September, 2014; originally announced September 2014.

Comments: 6 pages, 6 figures. Published in Astronomy & Astrophysics

Journal ref: A&A 571, A42 (2014)

Time-distance helioseismology: A new averaging scheme for measuring flow vorticity

Authors: Jan Langfellner, Laurent Gizon, Aaron C. Birch

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

Submitted 2 January, 2015; v1 submitted 20 August, 2014; originally announced August 2014.

Comments: 7 pages, 6 figures

Journal ref: A&A 570, A90 (2014)

Helioseismology with Solar Orbiter

Authors: Björn Löptien, Aaron C. Birch, Laurent Gizon, Jesper Schou, Thierry Appourchaux, Julián Blanco Rodríguez, Paul S. Cally, Carlos Dominguez-Tagle, Achim Gandorfer, Frank Hill, Johann Hirzberger, Philip H. Scherrer, Sami K. Solanki

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

Submitted 20 June, 2014; originally announced June 2014.

Comments: 35 pages, 12 figures. Accepted for publication in Space Science Reviews

Generalization of the noise model for time-distance helioseismology

Authors: Damien Fournier, Laurent Gizon, Thorsten Hohage, Aaron C. Birch

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

Submitted 20 June, 2014; originally announced June 2014.

MSC Class: 62M10 (Primary) 60H40 (Secondary)

Journal ref: A&A 567, A137 (2014)

Evaluation of the capability of local helioseismology to discern between monolithic and spaghetti sunspot models

Authors: T. Felipe, A. D. Crouch, A. C. Birch

Abstract: The helioseismic properties of the wave scattering generated by monolithic and spaghetti sunspots are analyzed by means of numerical simulations. In these computations, an incident f or p1 mode travels through the sunspot model, which produces absorption and phase shift of the waves. The scattering is studied by inspecting the wavefield, computing travel-time shifts, and performing Fourier-Hankel… ▽ More The helioseismic properties of the wave scattering generated by monolithic and spaghetti sunspots are analyzed by means of numerical simulations. In these computations, an incident f or p1 mode travels through the sunspot model, which produces absorption and phase shift of the waves. The scattering is studied by inspecting the wavefield, computing travel-time shifts, and performing Fourier-Hankel analysis. The comparison between the results obtained for both sunspot models reveals that the differences in the absorption coefficient can be detected above noise level. The spaghetti model produces an steep increase of the phase shift with the degree of the mode at short wavelengths, while mode-mixing is more efficient for the monolithic model. These results provide a clue for what to look for in solar observations to discern the constitution of sunspots between the proposed monolithic and spaghetti models. △ Less

Submitted 30 April, 2014; originally announced May 2014.

Comments: Accepted for publication in The Astrophysical Journal

Interpreting the Helioseismic and Magnetic Imager (HMI) Multi-Height Velocity Measurements

Authors: Kaori Nagashima, Björn Löptien, Laurent Gizon, Aaron C. Birch, Robert Cameron, Sebastien Couvidat, Sanja Danilovic, Bernhard Fleck, Robert Stein

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

Submitted 14 April, 2014; originally announced April 2014.

Comments: accepted for publication in Solar Physics

Journal ref: Solar Physics 289 (2014) 3457

Propagating Linear Waves in Convectively Unstable Stellar Models: a Perturbative Approach

Authors: Emanuele Papini, Laurent Gizon, Aaron C. Birch

Abstract: Linear time-domain simulations of acoustic oscillations are unstable in the stellar convection zone. To overcome this problem it is customary to compute the oscillations of a stabilized background stellar model. The stabilization, however, affects the result. Here we propose to use a perturbative approach (running the simulation twice) to approximately recover the acoustic wave field, while preser… ▽ More Linear time-domain simulations of acoustic oscillations are unstable in the stellar convection zone. To overcome this problem it is customary to compute the oscillations of a stabilized background stellar model. The stabilization, however, affects the result. Here we propose to use a perturbative approach (running the simulation twice) to approximately recover the acoustic wave field, while preserving seismic reciprocity. To test the method we considered a 1D standard solar model. We found that the mode frequencies of the (unstable) standard solar model are well approximated by the perturbative approach within $1~μみゅー$Hz for low-degree modes with frequencies near $3~μみゅー$Hz. We also show that the perturbative approach is appropriate for correcting rotational-frequency kernels. Finally, we comment that the method can be generalized to wave propagation in 3D magnetized stellar interiors because the magnetic fields have stabilizing effects on convection. △ Less

Submitted 15 December, 2013; originally announced December 2013.

Comments: 14 pages. Published online in Solar Physics, available at http://link.springer.com/article/10.1007/s11207-013-0457-7