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

Aaron C. Birch

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

ResearcherUnclaimed source record
astro-ph.SRphysics.flu-dynApplicationsnlin.PS

Catalog footprint

What is connected

20works
4topics
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

20 published item(s)

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

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

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

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

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

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

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

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

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