BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Laurent Gizon

Laurent Gizon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPphysics.flu-dynastro-ph.IMApplicationsastro-ph.GAmath-phmath.MPnlin.PSphysics.ao-phphysics.space-ph

Catalog footprint

What is connected

45works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2023arXiv

Solaris: A Focused Solar Polar Discovery-class Mission to achieve the Highest Priority Heliophysics Science Now

Solaris is a transformative Solar Polar Discovery-class mission concept to address crucial outstanding questions that can only be answered from a polar vantage. Solaris will image the Sun's poles from ~75 degree latitude, providing new insight into the workings of the solar dynamo and the solar cycle, which are at the foundation of our understanding of space weather and space climate. Solaris will also provide enabling observations for improved space weather research, modeling and prediction, revealing a unique, new view of the corona, coronal dynamics and CME eruptions from above.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Contribution of flows around active regions to the north-south helioseismic travel-time measurements

Context. In local helioseismology, the travel times of acoustic waves propagating in opposite directions along the same meridian inform us about horizontal flows in the north-south direction. The longitudinal averages of the north-south helioseismic travel-time shifts vary with the sunspot cycle. Aims. We aim to study the contribution of inflows into solar active regions to this solar-cycle variation. Methods. To do so, we identify the local flows around active regions in the horizontal flow maps obtained from correlation tracking of granulation in SDO/HMI continuum images. We compute the forward-modeled travel-time perturbations caused by these inflows using 3D sensitivity kernels. In order to compare with the observations, we average these forward-modeled travel-time perturbations over longitude and time in the same way as the measured travel times. Results. The forward-modeling approach shows that the inflows associated with active regions may account for only a fraction of the solar-cycle variations in the north-south travel-time measurements. Conclusions. The travel-time perturbations caused by the large-scale inflows surrounding the active regions do not explain in full the solar-cycle variations seen in the helioseismic measurements of the meridional circulation. Keywords: Sun: activity -- Sun: helioseismology

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evolution of Dipolar Mixed-mode Coupling Factor in Red Giant Stars: Impact of Buoyancy Spike

Mixed modes observed in red giants allow for investigation of the stellar interior structures. One important feature in these structures is the buoyancy spike caused by the discontinuity of the chemical gradient left behind during the first dredge-up. The buoyancy spike emerges at the base of the convective zone in low-luminosity red giants and later becomes a glitch when the g-mode cavity expands to encompass the spike. Here, we study the impact of the buoyancy spike on the dipolar mixed modes using stellar models with different properties. We find that the applicability of the asymptotic formalisms for the coupling factor, q, varies depending on the location of the evanescent zone, relative to the position of the spike. Significant deviations between the value of q inferred from fitting the oscillation frequencies and either of the formalisms proposed in the literature are found in models with a large frequency separation in the interval 5 to 15 microHz, with evanescent zones located in a transition region which may be thin or thick. However, it is still possible to reconcile q with the predictions from the asymptotic formalisms, by choosing which formalism to use according to the value of q. For stars approaching the luminosity bump, the buoyancy spike becomes a glitch and strongly affects the mode frequencies. Fitting the frequencies without accounting for the glitch leads to unphysical variations in the inferred q, but we show that this is corrected when properly accounting for the glitch in the fitting.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Imaging individual active regions on the Sun's far side with improved helioseismic holography

Helioseismic holography is a useful method to detect active regions on the Sun's far side and improve space weather forecasts. We aim to improve helioseismic holography by using a clear formulation of the problem, an accurate forward solver in the frequency domain, and a better understanding of the noise properties. Building on the work of Lindsey et al., we define the forward- and backward-propagated wave fields (ingression and egression) in terms of a Green's function. This Green's function is computed using an accurate forward solver in the frequency domain. We analyse overlapping segments of 31 hr of SDO/HMI dopplergrams, with a cadence of 24 hr. Phase shifts between the ingression and the egression are measured and averaged to detect active regions on the far side. The phase maps are compared with direct EUV intensity maps from STEREO/EUVI. We confirm that medium-size active regions can be detected on the far side with high confidence. Their evolution (and possible emergence) can be monitored on a daily time scale. Seismic maps averaged over 3 days provide an active region detection rate as high as 75% and a false discovery rate only as low as 7%, for active regions with areas above one thousandth of an hemisphere. For a large part, these improvements can be attributed to the use of a complete Green's function (all skips) and to the use of all observations on the front side (full pupil). Improved helioseismic holography enables the study of the evolution of medium-size active regions on the Sun's far side.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Theory of solar oscillations in the inertial frequency range: Linear modes of the convection zone

On the one hand, several types of global-scale inertial modes of oscillation have been observed on the Sun. They include the equatorial Rossby modes, critical-latitude modes, and high-latitude modes. On the other hand, the columnar convective modes (predicted by simulations; also known as banana cells or thermal Rossby waves) remain elusive. We aim to investigate the influence of turbulent diffusivities, non-adiabatic stratification, differential rotation, and a latitudinal entropy gradient on the linear global modes of the rotating solar convection zone. We solve numerically for the eigenmodes of a rotating compressible fluid inside a spherical shell. We identify modes in the inertial frequency range including the columnar convective modes, as well as modes of mixed character. The corresponding mode dispersion relations and eigenfunctions are computed for azimuthal orders $m \leq 16$. The three main results are as follows. Firstly, we find that, for $m \gtrsim 5$, the radial dependence of the equatorial Rossby modes with no radial node ($n=0$) is radically changed from the traditional expectation ($r^m$) for turbulent diffusivities $\gtrsim 10^{12}$ cm$^2$ s$^{-1}$. Secondly, we find mixed modes, i.e. modes that share properties of the equatorial Rossby modes with one radial node ($n=1$) and the columnar convective modes. Thirdly, we show that the $m=1$ high-latitude mode in the model is consistent with the solar observations when the latitudinal entropy gradient corresponding to a thermal wind balance is included (baroclinally unstable mode). To our knowledge, this work is the first realistic eigenvalue calculation of the global modes of the rotating solar convection zone. This calculation reveals a rich spectrum of modes in the inertial frequency range, which can be directly compared to the observations. In turn, the observed modes can inform us about the solar convection zone.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Viscous inertial modes on a differentially rotating sphere: Comparison with solar observations

In a previous paper we studied the effect of latitudinal rotation on solar equatorial Rossby modes in the beta-plane approximation. Since then, a rich spectrum of inertial modes has been observed on the Sun, which is not limited to the equatorial Rossby modes and includes high-latitude modes. Here we extend the computation of toroidal modes in 2D to spherical geometry, using realistic solar differential rotation and including viscous damping. The aim is to compare the computed mode spectra with the observations and to study mode stability. At fixed radius, we solve the eigenvalue problem numerically using a spherical harmonics decomposition of the velocity stream function. Due to the presence of viscous critical layers, the spectrum consists of four different families: Rossby modes, high-latitude modes, critical-latitude modes, and strongly damped modes. For each longitudinal wavenumber m<4, up to three Rossby-like modes are present on the sphere, in contrast to the equatorial beta plane where only the equatorial Rossby mode is present. The least damped modes in the model have eigenfrequencies and eigenfunctions that resemble the observed modes; the comparison improves when the radius is taken in the lower half of the convection zone. For radii above 0.75R and Ekman numbers E<10^{-4}, at least one mode is unstable. For either m=1 or m=2, up to two Rossby modes are unstable when the radial dependence of the Ekman number follows a quenched diffusivity model (E=2. 10^{-5} at the base of the convection zone). For m=3, up to two Rossby modes can be unstable, including the equatorial Rossby mode. Although the 2D model discussed here is highly simplified, the spectrum of toroidal modes appears to include many of the observed solar inertial modes. The self-excited modes in the model have frequencies close to those of the observed modes with the largest amplitudes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asteroseismology of luminous red giants with Kepler I: Long Period Variables with radial and non-radial modes

While long period variables (LPVs) have been extensively investigated, especially with MACHO and OGLE data for the Magellanic Clouds, there still exist open questions in their pulsations regarding the excitation mechanisms, radial order and angular degree assignment. Here, we perform asteroseismic analyses on LPVs observed by the 4-year Kepler mission. Using a cross-correlation method, we detect unambiguous pulsation ridges associated with radial fundamental modes ($n=1$) and overtones ($n\geqslant2$), where the radial order assignment is made by using theoretical frequencies and observed frequencies. Our results confirm that the amplitude variability seen in semiregulars is consistent with oscillations being solar-like. We identify that the dipole modes, $l=1$, are dominant in the radial orders of $3\leq n \leq6$, and that quadrupole modes, $l=2$, are dominant in the first overtone $n=2$. A test of seismic scaling relations using Gaia DR2 parallaxes reveals the possibility that the relations break down when $ν_{\rm max}$ $\lesssim$ 3 $μ$Hz (R $\gtrsim$ 40 R$_{\odot}$, or log $\rm L/L_{\odot}$ $\gtrsim$ 2.6). Our homogeneous measurements of pulsation amplitude and period for 3213 LPVs will be very valuable for probing effects of pulsation on mass loss, in particular in those stars with periods around 60 days, which has been argued as a threshold of substantial pulsation-triggered mass loss.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asteroseismology of luminous red giants with Kepler. II. Dependence of mass loss on pulsations and radiation

Mass loss by red giants is an important process to understand the final stages of stellar evolution and the chemical enrichment of the interstellar medium. Mass-loss rates are thought to be controlled by pulsation-enhanced dust-driven outflows. Here we investigate the relationships between mass loss, pulsations, and radiation, using 3213 luminous Kepler red giants and 135000 ASAS-SN semiregulars and Miras. Mass-loss rates are traced by infrared colours using 2MASS and WISE and by observed-to-model WISE fluxes, and are also estimated using dust mass-loss rates from literature assuming a typical gas-to-dust mass ratio of 400. To specify the pulsations, we extract the period and height of the highest peak in the power spectrum of oscillation. Absolute magnitudes are obtained from the 2MASS Ks band and the Gaia DR2 parallaxes. Our results follow. (i) Substantial mass loss sets in at pulsation periods above ~60 and ~100 days, corresponding to Asymptotic-Giant-Branch stars at the base of the period-luminosity sequences C' and C. (ii) The mass-loss rate starts to rapidly increase in semiregulars for which the luminosity is just above the Red-Giant-Branch tip and gradually plateaus to a level similar to that of Miras. (iii) The mass-loss rates in Miras do not depend on luminosity, consistent with pulsation-enhanced dust-driven winds. (iv) The accumulated mass loss on the Red Giant Branch consistent with asteroseismic predictions reduces the masses of red-clump stars by 6.3%, less than the typical uncertainty on their asteroseismic masses. Thus mass loss is currently not a limitation of stellar age estimates for galactic archaeology studies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterizing the spatial pattern of solar supergranulation using the bispectrum

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detection of exomoons in simulated light curves with a regularized convolutional neural network

Many moons have been detected around planets in our Solar System, but none has been detected unambiguously around any of the confirmed extrasolar planets. We test the feasibility of a supervised convolutional neural network to classify photometric transit light curves of planet-host stars and identify exomoon transits, while avoiding false positives caused by stellar variability or instrumental noise. Convolutional neural networks are known to have contributed to improving the accuracy of classification tasks. The network optimization is typically performed without studying the effect of noise on the training process. Here we design and optimize a 1D convolutional neural network to classify photometric transit light curves. We regularize the network by the total variation loss in order to remove unwanted variations in the data features. Using numerical experiments, we demonstrate the benefits of our network, which produces results comparable to or better than the standard network solutions. Most importantly, our network clearly outperforms a classical method used in exoplanet science to identify moon-like signals. Thus the proposed network is a promising approach for analyzing real transit light curves in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Exomoon indicators in high-precision transit light curves

While the solar system contains about 20 times more moons than planets, no moon has been confirmed around any of the thousands of extrasolar planets known so far. Tools for an uncomplicated identification of the most promising exomoon candidates could be beneficial to streamline follow-up studies.} Here we study three exomoon indicators that emerge if well-established planet-only models are fitted to a planet-moon transit light curve: transit timing variations (TTVs), transit duration variations (TDVs), and apparent planetary transit radius variations (TRVs). We re-evaluate under realistic conditions the previously proposed exomoon signatures in the TTV and TDV series. We simulate light curves of a transiting exoplanet with a single moon. These model light curves are then fitted with a planet-only transit model, pretending there were no moon, and we explore the resulting TTV, TDV, and TRV series for evidence of the moon. The previously described ellipse in the TTV-TDV diagram of an exoplanet with a moon emerges only for high-density moons. Low-density moons distort the sinusoidal shapes of the TTV and the TDV series due to their photometric contribution to the combined planet-moon transit. Sufficiently large moons can produce periodic apparent TRVs of their host planets that could be observable. We find that Kepler and PLATO have similar performances in detecting the exomoon-induced TRV effect around simulated bright ($m_V=8$) stars. These stars, however, are rare in the Kepler sample but will be abundant in the PLATO sample. Moreover, PLATO's higher cadence yields a stronger TTV signal. The periodogram of the sequence of transit radius measurements can indicate the presence of a moon. The TTV and TDV series of exoplanets with moons can be more complex than previously assumed. We propose that TRVs could be a more promising means to identify exomoons in large exoplanet surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solar Rossby waves observed in GONG++ ring-diagram flow maps

Context: Solar sectoral Rossby waves have only recently been unambiguously identified in Helioseimsic and Magnetic Imager (HMI) and Michelson Doppler Imager (MDI) maps of flows near the solar surface. So far this has not been done with the Global Oscillation Network Group (GONG) ground-based observations, which have different noise properties. Aims: We utilize 17 years of GONG++ data, to identify and characterize solar Rossby waves using ring-diagram helioseismology. We compare directly with HMI ring-diagram analysis. Methods: Maps of the radial vorticity are obtained for flows within the top 2 Mm of the surface for 17 years of GONG++. The data is corrected for systematic effects including the annual periodicity related to the $B_0$ angle. We then compute the Fourier components of the radial vorticity of the flows in the co-rotating frame. We perform the same analysis on the HMI data that overlap in time. Results: We find that the solar Rossby waves have measurable amplitudes in the GONG++ sectoral power spectra for azimuthal orders between $m=3$ and $m=15$. The measured mode characteristics (frequencies, lifetimes and amplitudes) from GONG++ are consistent with the HMI measurements in the overlap period from 2010 to 2018 for $m\le9$. For higher-$m$ modes the amplitudes and frequencies agree within two sigmas. The signal-to-noise ratio of modes in GONG++ power spectra is comparable to HMI for $8\le m\le11$, but is lower by a factor of two for other modes. Conclusions: The GONG++ data provide a long and uniform data set to study solar global-scale Rossby waves from 2001.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solar-cycle irradiance variations over the last four billion years

The variability of the spectral solar irradiance (SSI) over the course of the 11-year solar cycle is one of the manifestations of solar magnetic activity. There is a strong evidence that the SSI variability has an effect on the Earth's atmosphere. The faster rotation of the Sun in the past lead to a more vigorous action of solar dynamo and thus potentially to larger amplitude of the SSI variability on the timescale of the solar activity cycle. This could led to a stronger response of the Earth's atmosphere as well as other solar system planets' atmospheres to the solar activity cycle. We calculate the amplitude of the SSI and TSI variability over the course of the solar activity cycle as a function of solar age. We employ the relationship between the stellar magnetic activity and the age based on observations of solar twins. Using this relation we reconstruct solar magnetic activity and the corresponding solar disk area coverages by magnetic features (i.e. spots and faculae) over the last four billion years. These disk coverages are then used to calculate the amplitude of the solar-cycle SSI variability as a function of wavelength and solar age. Our calculations show that the young Sun was significantly more variable than the present Sun. The amplitude of the solar-cycle Total Solar Irradiance (TSI) variability of the 600 Myr old Sun was about 10 times larger than that of the present Sun. Furthermore, the variability of the young Sun was spot-dominated (the Sun being brighter at the activity minimum than in the maximum), i.e. the Sun was overall brighter at activity minima than at maxima. The amplitude of the TSI variability decreased with solar age until it reached a minimum value at 2.8 Gyr. After this point, the TSI variability is faculae-dominated (the Sun is brighter at the activity maximum) and its amplitude increases with age.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shape of a slowly rotating star measured by asteroseismology

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Solar-cycle variation of the rotational shear near the solar surface

Helioseismology has revealed that the angular velocity of the Sun increases with depth in the outermost 35 Mm of the Sun. Recently, we have shown that the logarithmic radial gradient ($\rm d\lnΩ/\rm d\ln r $) in the upper 10~Mm is close to $-1$ from the equator to $60^\circ$ latitude.We aim to measure the temporal variation of the rotational shear over solar cycle 23 and the rising phase of cycle 24 (1996-2015). We used f mode frequency splitting data spanning 1996 to 2011 from the Michelson Doppler Imager (MDI) and 2010 to 2015 from the Helioseismic Magnetic Imager (HMI). In a first for such studies, the f mode frequency splitting data were obtained from 360-day time series. We used the same method as in our previous work for measuring $\rm d\lnΩ/d\ln r $ from the equator to $80^\circ$ latitude in the outer 13~Mm of the Sun. Then, we calculated the variation of the gradient at annual cadence relative to the average over 1996 to 2015. We found the rotational shear at low latitudes ($0^\circ$ to $30^\circ$) to vary in-phase with the solar activity, varying by $\sim \pm 10$\% over the period 1996 to 2015. At high latitudes ($60^\circ$ to $80^\circ$), we found rotational shear to vary in anti-phase with the solar activity. By comparing the radial gradient obtained from the splittings of the 360-day and the corresponding 72-day time series of HMI and MDI data, we suggest that the splittings obtained from the 72-day HMI time series suffer from systematic errors. We provide a quantitative measurement of the temporal variation of the outer part of the near surface shear layer which may provide useful constraints on dynamo models and differential rotation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asteroseismic inversions for radial differential rotation of Sun-like stars: ensemble fits

Radial differential rotation is an important parameter for stellar dynamo theory and for understanding angular momentum transport. We investigate the potential of using a large number of similar stars simultaneously to constrain their average radial differential rotation gradient: we call this 'ensemble fitting'. We use a range of stellar models along the main sequence, each with a synthetic rotation profile. The rotation profiles are step functions with a step of -350 nHz, which is located at the base of the convection zone. These models are used to compute the rotational splittings of the p modes and to model their uncertainties. We then fit an ensemble of stars to infer the average step size. All the uncertainties on the inferred step size for individual stars are of the order 1 micro Hz. Using 15 stellar models in an ensemble fit, we show that the uncertainty on the average step size is reduced to less than the input step size, which allows us to constrain the sign of the radial differential rotation. We show that a solar-like step size (approximately 30 nHz) can be constrained by an ensemble fit of thousands of main-sequence stars. Observing the number of stars required to successfully exploit the ensemble fitting method will be possible with future asteroseismology missions, such as PLATO. We demonstrate the potential of ensemble fitting by showing that any systematic differences in the average step size between F, G, and K-type stars larger than 100 nHz can be detected.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Data compression for local correlation tracking of solar granulation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rotation, differential rotation, and gyrochronology of active Kepler stars

The high-precision photometry from the CoRoT and Kepler satellites has led to measurements of surface rotation periods for tens of thousands of stars. Our main goal is to derive ages of thousands of field stars using consistent rotation period measurements in different gyrochronology relations. Multiple rotation periods are interpreted as surface differential rotation (DR). We re-analyze the sample of 24,124 Kepler stars from Reinhold et al. (2013) using different approaches based on the Lomb-Scargle periodogram. Each quarter (Q1-Q14) is treated individually using a prewhitening approach. Additionally, the full time series, and different segments thereof are analyzed. For more than 18,500 stars our results are consistent with the rotation periods from McQuillan et al. (2014). Thereof, more than 12,300 stars show multiple significant peaks, which we interpret as DR. Gyrochronology ages between 100 Myr and 10 Gyr were derived for more than 17,000 stars using different gyrochronology relations. We find a bimodal age distribution for Teff between 3200-4700 K. The derived ages reveal an empirical activity-age relation using photometric variability as stellar activity proxy. Additionally, we found 1079 stars with extremely stable (mostly short) periods. Half of these periods may be associated with rotation stabilized by non-eclipsing companions, the other half might be due to pulsations. The derived gyrochronology ages are well constrained since more than 93.0 % of the stars seem to be younger than the Sun where calibration is most reliable. Explaining the bimodality in the age distribution is challenging, and limits accurate stellar age predictions. The existence of cool stars with almost constant rotation period over more than three years of observation might be explained by synchronization with stellar companions, or a dynamo mechanism keeping the spot configurations extremely stable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Seismic Sounding of Convection in the Sun

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulating acoustic waves in spotted stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The EChO science case

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalization of the noise model for time-distance helioseismology

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Helioseismic and Magnetic Imager Multi-height Dopplergrams

We study Doppler velocity measurements at multiple heights in the solar atmosphere using a set of six filtergrams obtained by the Helioseismic magnetic Imager on board the Solar Dynamics Observatory. There are clear and significant phase differences between core and wing Dopplergrams in the frequency range above the photospheric acoustic cutoff frequency, which indicates that these are really "multi-height" datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Helioseismology with Solar Orbiter

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Image compression in local helioseismology

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The radial gradient of the near-surface shear layer of the Sun

Helioseismology has provided unprecedented information about the internal rotation of the Sun. One of the important achievements was the discovery of two radial shear layers: one near the bottom of the convection zone (the tachocline) and one near the surface. These shear layers may be important ingredients for explaining the magnetic cycle of the Sun. We measure the logarithmic radial gradient of the rotation rate ($\rm d\lnΩ/\rm d\ln r$) near the surface of the Sun using 15 years of f mode rotational frequency splittings from the Michelson Doppler Imager (MDI) and four years of data from the Helioseismic and Magnetic Imager (HMI). We model the angular velocity of the Sun in the upper $\sim 10$ Mm as changing linearly with depth and use a multiplicative optimally localized averaging inversion to infer the gradient of the rotation rate as a function of latitude. Both the MDI and HMI data show that $\rm d\lnΩ/\rm d\ln r$ is close to $-1$ from the equator to 60$^{\circ}$ latitude and stays negative up to 75$^{\circ}$ latitude. However, the value of the gradient is different for MDI and HMI for latitudes above $60^{\circ}$. Additionally, there is a significant difference between the value of $d\lnΩ/d\ln r$ using an older and recently reprocessed MDI data for latitudes above $30^\circ$. We could reliably infer the value of $\rm d\lnΩ/\rm d\ln r$ up to 60$^{\circ}$, but not above this latitude, which will hopefully constrain theories of the near-surface shear layer and dynamo. Furthermore, the recently reprocessed MDI splitting data are more reliable than the older versions which contained clear systematic errors in the high degree f modes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Helioseismology of sunspots: how sensitive are travel times to the Wilson depression and to the subsurface magnetic field?

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Precise and accurate interpolated stellar oscillation frequencies on the main sequence

High-quality data from space-based observatories present an opportunity to fit stellar models to observations of individually-identified oscillation frequencies, not just the large and small frequency separations. But such fits require the evaluation of a large number of accurate stellar models, which remains expensive. Here, we show that global-mode oscillation frequencies interpolated in a grid of stellar models are precise and accurate, at least in the neighbourhood of a solar model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Propagation of seismic waves through a spatio-temporally fluctuating medium: Homogenization

Measurements of seismic wave travel times at the photosphere of the Sun have enabled inferences of its interior structure and dynamics. In interpreting these measurements, the simplifying assumption that waves propagate through a temporally stationary medium is almost universally invoked. However, the Sun is in a constant state of evolution, on a broad range of spatio-temporal scales. At the zero wavelength limit, i.e., when the wavelength is much shorter than the scale over which the medium varies, the WKBJ (ray) approximation may be applied. Here, we address the other asymptotic end of the spectrum, the infinite wavelength limit, using the technique of homogenization. We apply homogenization to scenarios where waves are propagating through rapidly varying media (spatially and temporally), and derive effective models for the media. One consequence is that a scalar sound speed becomes a tensorial wavespeed in the effective model and anisotropies can be induced depending on the nature of the perturbation. The second term in this asymptotic two-scale expansion, the so-called corrector, contains contributions due to higher-order scattering, leading to the decoherence of the wavefield. This decoherence may be causally linked to the observed wave attenuation in the Sun. Although the examples we consider here consist of periodic arrays of perturbations to the background, homogenization may be extended to ergodic and stationary random media. This method may have broad implications for the manner in which we interpret seismic measurements in the Sun and for modeling the effects of granulation on the scattering of waves and distortion of normal-mode eigenfunctions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Helioseismology challenges models of solar convection

Convection is the mechanism by which energy is transported through the outermost 30% of the Sun. Solar turbulent convection is notoriously difficult to model across the entire convection zone where the density spans many orders of magnitude. In this issue of PNAS, Hanasoge et al. (2012) employ recent helioseismic observations to derive stringent empirical constraints on the amplitude of large-scale convective velocities in the solar interior. They report an upper limit that is far smaller than predicted by a popular hydrodynamic numerical simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Seismic probes of solar interior magnetic structure

Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of changes to isotropic wavespeeds and flow-advection related Doppler shifts. Here, we employ PDE-constrained optimization to determine the appropriate parameterization of the structural properties of the magnetic interior. Seven different wavespeeds fully characterize helioseismic wave propagation: the isotropic sound speed, a Doppler-shifting flow-advection velocity and an anisotropic magnetic velocity. The structure of magnetic media is sensed by magnetoacoustic slow and fast modes and Alfvén waves, each of which propagates at a different wavespeed. We show that even in the case of weak magnetic fields, significant errors may be incurred if these anisotropies are not accounted for in inversions. Translation invariance is demonstrably lost. These developments render plausible the accurate seismic imaging of magnetoconvection in the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constructing and Characterising Solar Structure Models for Computational Helioseismology

In this paper, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science, 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose to make additional corrections to the sound speed of Model S before stabilisation. No stabilised model can be perfectly solar-like, so we present three stabilised models with slightly different eigenmodes. The models are appropriate to study the f and p1 to p4 modes with spherical harmonic degrees in the range from 400 to 900. Background model CSM has a modified pressure gradient for stabilisation and has eigenfrequencies within 2% of Model S. Model CSM_A has an additional 10% increase in sound speed in the top 1 Mm resulting in eigenfrequencies within 2% of Model S and eigenfunctions that are, in comparison with CSM, closest to those of Model S. Model CSM_B has a 3% decrease in sound speed in the top 5 Mm resulting in eigenfrequencies within 1% of Model S and eigenfunctions that are only marginally adversely affected. These models are useful to study the interaction of solar waves with embedded three-dimensional heterogeneities, such as convective flows and model sunspots. We have also calculated the response of the stabilised models to excitation by random near-surface sources, using simulations of the propagation of linear waves. We find that the simulated power spectra of wave motion are in good agreement with an observed SOHO/MDI power spectrum. Overall, our convectively stabilised background models provide a good basis for quantitative numerical local helioseismology. The models are available for download from http://www.mps.mpg.de/projects/seismo/NA4/.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Adjoint Method Applied to Time-Distance Helioseismology

For a given {\it misfit function}, a specified optimality measure of a model, its gradient describes the manner in which one may alter properties of the system to march towards a stationary point. The adjoint method, arising from partial-differential-equation-constrained optimization, describes a means of extracting derivatives of a misfit function with respect to model parameters through finite computation. It relies on the accurate calculation of wavefields that are driven by two types of sources, namely the average wave-excitation spectrum, resulting in the {\it forward wavefield}, and differences between predictions and observations, resulting in an {\it adjoint wavefield}. All sensitivity kernels relevant to a given measurement emerge directly from the evaluation of an interaction integral involving these wavefields. The technique facilitates computation of sensitivity kernels (Fréchet derivatives) relative to three-dimensional heterogeneous background models, thereby paving the way for non-linear iterative inversions. An algorithm to perform such inversions using as many observations as desired is discussed.

0 citations0 reviewsNo code linkedOpen access