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

L. Gizon

L. 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.IMastro-ph.EPastro-phastro-ph.GAastro-ph.HEphysics.ins-detphysics.space-ph

Catalog footprint

What is connected

54works
8topics
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

54 published item(s)

preprint2022arXiv

Faculae cancel out on the surfaces of active Suns

Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The magnetic drivers of campfires seen by the Polarimetric and Helioseismic Imager (PHI) on Solar Orbiter

The Extreme Ultraviolet Imager (EUI) on board the Solar Orbiter (SO) spacecraft observed small extreme ultraviolet (EUV) bursts, termed campfires, that have been proposed to be brightenings near the apexes of low-lying loops in the quiet-Sun atmosphere. The underlying magnetic processes driving these campfires are not understood. During the cruise phase of SO and at a distance of 0.523\,AU from the Sun, the Polarimetric and Helioseismic Imager on Solar Orbiter (SO/PHI) observed a quiet-Sun region jointly with SO/EUI, offering the possibility to investigate the surface magnetic field dynamics underlying campfires at a spatial resolution of about 380~km. In 71\% of the 38 isolated events, campfires are confined between bipolar magnetic features, which seem to exhibit signatures of magnetic flux cancellation. The flux cancellation occurs either between the two main footpoints, or between one of the footpoints of the loop housing the campfire and a nearby opposite polarity patch. In one particularly clear-cut case, we detected the emergence of a small-scale magnetic loop in the internetwork followed soon afterwards by a campfire brightening adjacent to the location of the linear polarisation signal in the photosphere, that is to say near where the apex of the emerging loop lays. The rest of the events were observed over small scattered magnetic features, which could not be identified as magnetic footpoints of the campfire hosting loops. The majority of campfires could be driven by magnetic reconnection triggered at the footpoints, similar to the physical processes occurring in the burst-like EUV events discussed in the literature. About a quarter of all analysed campfires, however, are not associated to such magnetic activity in the photosphere, which implies that other heating mechanisms are energising these small-scale EUV brightenings.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The on-ground data reduction and calibration pipeline for SO/PHI-HRT

The ESA/NASA Solar Orbiter space mission has been successfully launched in February 2020. Onboard is the Polarimetric and Helioseismic Imager (SO/PHI), which has two telescopes, a High Resolution Telescope (HRT) and the Full Disc Telescope (FDT). The instrument is designed to infer the photospheric magnetic field and line-of-sight velocity through differential imaging of the polarised light emitted by the Sun. It calculates the full Stokes vector at 6 wavelength positions at the Fe I 617.3 nm absorption line. Due to telemetry constraints, the instrument nominally processes these Stokes profiles onboard, however when telemetry is available, the raw images are downlinked and reduced on ground. Here the architecture of the on-ground pipeline for HRT is presented, which also offers additional corrections not currently available on board the instrument. The pipeline can reduce raw images to the full Stokes vector with a polarimetric sensitivity of $10^{-3}\cdot I_{c}$ or better.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The PLATO field selection process I. Identification and content of the long-pointing fields

PLATO (PLAnetary Transits and Oscillations of stars) is an ESA M-class satellite planned for launch by end 2026 and dedicated to the wide-field search of transiting planets around bright and nearby stars, with a strong focus on discovering habitable rocky planets hosted by solar-like stars. The choice of the fields to be pointed at is a crucial task since it has a direct impact on the scientific return of the mission. In this paper we describe and discuss the formal requirements and the key scientific prioritization criteria that have to be taken into account in the Long-duration Observation Phase (LOP) field selection, and apply a quantitative metric to guide us in this complex optimization process. We identify two provisional LOP fields, one for each hemisphere (LOPS1, LOPN1), and discuss their properties and stellar content. While additional fine-tuning shall be applied to LOP selection before the definitive choice (to be made two years before launch), we expect their position will not move by more than a few degrees with respect to what is proposed in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Power spectrum of turbulent convection in the solar photosphere

The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An all-sky catalog of solar-type dwarfs for exoplanetary transit surveys

Most future surveys designed to discover transiting exoplanets, including TESS and PLATO, will target bright (V<13) and nearby solar-type stars having a spectral type later than F5. In order to enhance the probability of identifying transits, these surveys must cover a very large area on the sky, because of the intrinsically low areal density of bright targets. Unfortunately, no existing catalog of stellar parameters is both deep and wide enough to provide a homogeneous input list. As the first Gaia data release exploitable for this purpose is expected to be released not earlier than late 2017, we have devised an improved reduced-proper-motion method to discriminate late field dwarfs and giants by combining UCAC4 proper motions with APASS DR6 photometry, and relying on RAVE DR4 as an external calibrator. The output, named UCAC4-RPM, is a publicly-available, complete all-sky catalog of solar-type dwarfs down to V<13.5, plus an extension to logg>3.0 subgiants. The relatively low amount of contamination (defined as the fraction of false positives; <30%) also makes UCAC4-RPM a useful tool for the past and ongoing ground-based transit surveys, which need to discard candidate signals originating from early-type or giant stars. As an application, we show how UCAC4-RPM may support the preparation of the TESS (that will map almost the entire sky) input catalog and the input catalog of PLATO, planned to survey more than half of the whole sky with exquisite photometric precision.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MESA meets MURaM: Surface effects in main-sequence solar-like oscillators computed using three-dimensional radiation hydrodynamics simulations

... [C]urrent stellar models predict oscillation frequencies that are systematically affected by simplified modelling of the near-surface layers. We use three-dimensional radiation hydrodynamics simulations to better model the near-surface equilibrium structure of dwarfs with spectral types F3, G2, K0 and K5, and examine the differences between oscillation mode frequencies. ... We precisely match stellar models to the simulations' gravities and effective temperatures at the surface, and to the temporally- and horizontally-averaged densities and pressures at their deepest points. We then replace the near-surface structure with that of the averaged simulation and compute the change in the oscillation mode frequencies. We also fit the differences using several parametric models currently available in the literature. The surface effect in the stars of solar-type and later is qualitatively similar and changes steadily with decreasing effective temperature. In particular, the point of greatest frequency difference decreases slightly as a fraction of the acoustic cut-off frequency and the overall scale of the surface effect decreases. The surface effect in the hot, F3-type star follows the same trend in scale (i.e. it is larger in magnitude) but shows a different overall variation with mode frequency. We find that the two-term fit by Ball & Gizon (2014) is best able to reproduce the surface terms across all four spectral types, although the scaled solar term and a modified Lorentzian function also match the three cooler simulations reasonably well. ... Our simplified results suggest that the surface effect is generally larger in hotter stars (and correspondingly smaller in cooler stars) and of similar shape in stars of solar type and cooler. However, we cannot presently predict whether this will remain so when other components of the surface effect are included.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the uncertain nature of the core of $\boldsymbolα$ Cen A

High-quality astrometric, spectroscopic, interferometric and, importantly, asteroseismic observations are available for $α$ Cen A, which is the closest binary star system to earth. Taking all these constraints into account, we study the internal structure of the star by means of theoretical modelling. Using the Aarhus STellar Evolution Code (ASTEC) and the tools of Computational Bayesian Statistics, in particular a Markov chain Monte Carlo algorithm, we perform statistical inferences for the physical characteristics of the star. We find that $α$ Cen A has a probability of approximately 40% of having a convective core. This probability drops to few percents if one considers reduced rates for the $^{14}$N(p,$γ$)$^{15}$O reaction. These convective cores have fractional radii less than 8% when overshoot is neglected. Including overshooting also leads to the possibility of a convective core mostly sustained by the ppII chain energy output. We finally show that roughly 30% of the stellar models describing $α$ Cen A are in the subgiant regime.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SDO/HMI survey of emerging active regions for helioseismology

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anisotropy of the solar network magnetic field around the average supergranule

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asteroseismology of solar-type stars with K2

We present the first detections by the NASA K2 Mission of oscillations in solar-type stars, using short-cadence data collected during K2 Campaign\,1 (C1). We understand the asteroseismic detection thresholds for C1-like levels of photometric performance, and we can detect oscillations in subgiants having dominant oscillation frequencies around $1000\,\rm μHz$. Changes to the operation of the fine-guidance sensors are expected to give significant improvements in the high-frequency performance from C3 onwards. A reduction in the excess high-frequency noise by a factor of two-and-a-half in amplitude would bring main-sequence stars with dominant oscillation frequencies as high as ${\simeq 2500}\,\rm μHz$ into play as potential asteroseismic targets for K2.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraining differential rotation of Sun-like stars from asteroseismic and starspot rotation periods

In previous work we identified six Sun-like stars observed by Kepler with exceptionally clear asteroseismic signatures of rotation. Here, we show that five of these stars exhibit surface variability suitable for measuring rotation. In order to further constrain differential rotation, we compare the rotation periods obtained from light-curve variability with those from asteroseismology. The two rotation measurement methods are found to agree within uncertainties, suggesting that radial differential rotation is weak, as is the case for the Sun. Furthermore, we find significant discrepancies between ages from asteroseismology and from three different gyrochronology relations, implying that stellar age estimation is problematic even for Sun-like stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measuring stellar rotation periods with Kepler

We measure rotation periods for 12151 stars in the Kepler field, based on the photometric variability caused by stellar activity. Our analysis returns stable rotation periods over at least six out of eight quarters of Kepler data. This large sample of stars enables us to study the rotation periods as a function of spectral type. We find good agreement with previous studies and vsini measurements for F, G and K stars. Combining rotation periods, B-V color, and gyrochronology relations, we find that the cool stars in our sample are predominantly younger than ~1Gyr.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new correction of stellar oscillation frequencies for near-surface effects

Space-based observations of solar-like oscillations present an opportunity to constrain stellar models using individual mode frequencies. However, current stellar models are inaccurate near the surface, which introduces a systematic difference that must be corrected. We introduce and evaluate two parametrizations of the surface corrections based on formulae given by Gough (1990). The first we call a cubic term proportional to $ν^3/\mathcal{I}$ and the second has an additional inverse term proportional to $ν^{-1}/\mathcal{I}$, where $ν$ and $\mathcal{I}$ are the frequency and inertia of an oscillation mode. We first show that these formulae accurately correct model frequencies of two different solar models ... We then incorporate the parametrizations into a modelling pipeline that simultaneously fits the surface effects and the underlying stellar model parameters. We apply this pipeline to synthetic observations of a Sun-like stellar model, solar observations degraded to typical asteroseismic uncertainties, and observations of the well-studied CoRoT target HD52265. For comparison, we also run the pipeline with the scaled power-law correction proposed by Kjeldsen et al. (2008). The fits to synthetic and degraded solar data show that the method is unbiased and produces best-fit parameters that are consistent with the input models and known parameters of the Sun. Our results for HD52265 are consistent with previous modelling efforts and the magnitude of the surface correction is similar to that of the Sun. The fit using a scaled power-law correction is significantly worse but yields consistent parameters ... We find that the cubic term alone is suitable for asteroseismic applications and it is easy to implement in an existing pipeline. ... This parametrization is thus a useful new way to correct model frequencies so that observations of individual mode frequencies can be exploited.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetic Flux Transport at the Solar Surface

After emerging to the solar surface, the Sun's magnetic field displays a complex and intricate evolution. The evolution of the surface field is important for several reasons. One is that the surface field, and its dynamics, sets the boundary condition for the coronal and heliospheric magnetic fields. Another is that the surface evolution gives us insight into the dynamo process. In particular, it plays an essential role in the Babcock-Leighton model of the solar dynamo. Describing this evolution is the aim of the surface flux transport model. The model starts from the emergence of magnetic bipoles. Thereafter, the model is based on the induction equation and the fact that after emergence the magnetic field is observed to evolve as if it were purely radial. The induction equation then describes how the surface flows -- differential rotation, meridional circulation, granular, supergranular flows, and active region inflows -- determine the evolution of the field (now taken to be purely radial). In this paper, we review the modeling of the various processes that determine the evolution of the surface field. We restrict our attention to their role in the surface flux transport model. We also discuss the success of the model and some of the results that have been obtained using this model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rotational splitting as a function of mode frequency for six Sun-like stars

Asteroseismology offers the prospect of constraining differential rotation in Sun-like stars. Here we have identified six high signal-to-noise main-sequence Sun-like stars in the Kepler field, which all have visible signs of rotational splitting of their p-mode frequencies. For each star, we extract the rotational frequency splitting and inclination angle from separate mode sets (adjacent modes with l=2, 0, and 1) spanning the p-mode envelope. We use a Markov chain Monte Carlo method to obtain the best fit and errors associated with each parameter. We are able to make independent measurements of rotational splittings of ~8 radial orders for each star. For all six stars, the measured splittings are consistent with uniform rotation, allowing us to exclude large radial differential rotation. This work opens the possibility of constraining internal rotation of Sun-like stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Seismic analysis of HD43587Aa, a solar-like oscillator in a multiple system

Context. The object HD 43587Aa is a G0V star observed during the 145-day LRa03 run of the COnvection, ROtation and planetary Transits space mission (CoRoT), for which complementary High Accuracy Radial velocity Planet Searcher (HARPS) spectra with S/N>300 were also obtained. Its visual magnitude is 5.71, and its effective temperature is close to 5950 K. It has a known companion in a highly eccentric orbit and is also coupled with two more distant companions. Aims. We undertake a preliminary investigation of the internal structure of HD 43587Aa. Methods. We carried out a seismic analysis of the star, using maximum likelihood estimators and Markov Chain Monte Carlo methods. Results. We established the first table of the eigenmode frequencies, widths, and heights for HD 43587Aa. The star appears to have a mass and a radius slightly larger than the Sun, and is slightly older (5.6 Gyr). Two scenarios are suggested for the geometry of the star: either its inclination angle is very low, or the rotation velocity of the star is very low. Conclusions. A more detailed study of the rotation and of the magnetic and chromospheric activity for this star is needed, and will be the subject of a further study. New high resolution spectrometric observations should be performed for at least several months in duration.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Seismic constraints on the radial dependence of the internal rotation profiles of six Kepler subgiants and young red giants

Context : We still do not know which mechanisms are responsible for the transport of angular momentum inside stars. The recent detection of mixed modes that contain the signature of rotation in the spectra of Kepler subgiants and red giants gives us the opportunity to make progress on this issue. Aims: Our aim is to probe the radial dependance of the rotation profiles for a sample of Kepler targets. For this purpose, subgiants and early red giants are particularly interesting targets because their rotational splittings are more sensitive to the rotation outside the deeper core than is the case for their more evolved counterparts. Methods: We first extract the rotational splittings and frequencies of the modes for six young Kepler red giants. We then perform a seismic modeling of these stars using the evolutionary codes CESAM2k and ASTEC. By using the observed splittings and the rotational kernels of the optimal models, we perform inversions of the internal rotation profiles of the six stars. Results: We obtain estimates of the mean rotation rate in the core and in the convective envelope of these stars. We show that the rotation contrast between the core and the envelope increases during the subgiant branch. Our results also suggest that the core of subgiants spins up with time, contrary to the RGB stars whose core has been shown to spin down. For two of the stars, we show that a discontinuous rotation profile with a deep discontinuity reproduces the observed splittings significantly better than a smooth rotation profile. Interestingly, the depths that are found most probable for the discontinuities roughly coincide with the location of the H-burning shell, which separates the layers that contract from those that expand. These results will bring observational constraints to the scenarios of angular momentum transport in stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Structure and evolution of solar supergranulation using SDO/HMI data

Context: Studying the motions on the solar surface is fundamental for understanding how turbulent convection transports energy and how magnetic fields are distributed across the solar surface. Aims: From horizontal velocity measurements all over the visible disc of the Sun and using data from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI), we investigate the structure and evolution of solar supergranulation. Methods: Horizontal velocity fields were measured by following the proper motions of solar granules using a newly developed version of the coherent structure tracking (CST) code. With this tool, maps of horizontal divergence were computed. We then segmented and identified supergranular cells and followed their histories by using spatio-temporal labelling. With this dataset we derived the fundamental properties of supergranulation, including their motion. Results: We find values of the fundamental parameters of supergranulation similar to previous studies: a mean lifetime of 1.5 days and a mean diameter of 25~Mm. The tracking of individual supergranular cells reveals the solar differential rotation and a poleward circulation trend of the meridional flow. The shape of the derived differential rotation and meridional flow does not depend on the cell size. If there is a background magnetic field, the diverging flows in supergranules are weaker. Conclusions: This study confirms that supergranules are suitable tracers that may be used to investigate the large-scale flows of the solar convection as long as they are detectable enough on the surface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The PLATO 2.0 Mission

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Asteroseismic fundamental properties of solar-type stars observed by the NASA Kepler Mission

We use asteroseismic data obtained by the NASA Kepler Mission to estimate the fundamental properties of more than 500 main-sequence and sub-giant stars. Data obtained during the first 10 months of Kepler science operations were used for this work, when these solar-type targets were observed for one month each in a survey mode. Stellar properties have been estimated using two global asteroseismic parameters and complementary photometric and spectroscopic data. Homogeneous sets of effective temperatures were available for the entire ensemble from complementary photometry; spectroscopic estimates of T_eff and [Fe/H] were available from a homogeneous analysis of ground-based data on a subset of 87 stars. [Abbreviated version... see paper for full abstract.]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparison of solar horizontal velocity fields from SDO/HMI and Hinode data

The measurement of the Sun's surface motions with a high spatial and temporal resolution is still a challenge. We wish to validate horizontal velocity measurements all over the visible disk of the Sun from Solar Dynamics Observatory/ Helioseismic and Magnetic Imager (SDO/HMI) data. Horizontal velocity fields are measured by following the proper motions of solar granules using a newly developed version of the Coherent Structure Tracking (CST) code. The comparison of the surface flows measured at high spatial resolution (Hinode, 0.1 arcsec) and low resolution (SDO/HMI, 0.5 arcsec) allows us to determine corrections to be applied to the horizontal velocity measured from HMI white light data. We derive horizontal velocity maps with spatial and temporal resolutions of respectively 2.5 Mm and 30 min. From the two components of the horizontal velocity Vx and Vy measured in the sky plane and the simultaneous line of sight component from SDO/HMI dopplergrams v_D, we derive the spherical velocity components (Vr, Vtheta, Vphi). The azimuthal component Vphi gives the solar differential rotation with a high precision (+-0.037km/s) from a temporal sequence of only three hours. By following the proper motions of the solar granules, we can revisit the dynamics of the solar surface at high spatial and temporal resolutions from hours to months and years with the SDO data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparison of solar surface flows inferred from time--distance helioseismology and coherent structure tracking using HMI/SDO observations

We compare measurements of horizontal flows on the surface of the Sun using helioseismic time--distance inversions and coherent structure tracking of solar granules. Tracking provides 2D horizontal flows on the solar surface, whereas the time--distance inversions estimate the full 3-D velocity flows in the shallow near-surface layers. Both techniques use HMI observations as an input. We find good correlations between the various measurements resulting from the two techniques. Further, we find a good agreement between these measurements and the time-averaged Doppler line-of-sight velocity, and also perform sanity checks on the vertical flow that resulted from the 3-D time--distance inversion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First high-resolution images of the Sun in the 2796 Å Mg II k line

We present the first high-resolution solar images in the Mg II k 2796 Å line. The images, taken through a 4.8 Å broad interference filter, were obtained during the second science flight of SUNRISE in June 2013 by the SuFI instrument. The Mg II k images display structures that look qualitatively very similar to images taken in the core of Ca II H. The Mg II images exhibit reversed granulation (or shock waves) in the internetwork regions of the quiet Sun, at intensity contrasts that are similar to those found in Ca II H. Very prominent in Mg II are bright points, both in the quiet Sun and in plage regions, particularly near disk center. These are much brighter than at other wavelengths sampled at similar resolution. Furthermore, Mg II k images also show fibril structures associated with plage regions. Again, the fibrils are similar to those seen in Ca II H images, but tend to be more pronounced, particularly in weak plage.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kepler White Paper: Asteroseismology of Solar-Like Oscillators in a 2-Wheel Mission

We comment on the potential for continuing asteroseismology of solar-type and red-giant stars in a 2-wheel Kepler Mission. Our main conclusion is that by targeting stars in the ecliptic it should be possible to perform high-quality asteroseismology, as long as favorable scenarios for 2-wheel pointing performance are met. Targeting the ecliptic would potentially facilitate unique science that was not possible in the nominal Mission, notably from the study of clusters that are significantly brighter than those in the Kepler field. Our conclusions are based on predictions of 2-wheel observations made by a space photometry simulator, with information provided by the Kepler Project used as input to describe the degraded pointing scenarios. We find that elevated levels of frequency-dependent noise, consistent with the above scenarios, would have a significant negative impact on our ability to continue asteroseismic studies of solar-like oscillators in the Kepler field. However, the situation may be much more optimistic for observations in the ecliptic, provided that pointing resets of the spacecraft during regular desaturations of the two functioning reaction wheels are accurate at the < 1 arcsec level. This would make it possible to apply a post-hoc analysis that would recover most of the lost photometric precision. Without this post-hoc correction---and the accurate re-pointing it requires---the performance would probably be as poor as in the Kepler-field case. Critical to our conclusions for both fields is the assumed level of pointing noise (in the short-term jitter and the longer-term drift). We suggest that further tests will be needed to clarify our results once more detail and data on the expected pointing performance becomes available, and we offer our assistance in this work.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pathway to the Square Kilometre Array - The German White Paper -

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rotation periods of 12 000 main-sequence Kepler stars: Dependence on stellar spectral type and comparison with v sin i observations

Aims: We aim to measure the starspot rotation periods of active stars in the Kepler field as a function of spectral type and to extend reliable rotation measurements from F-, G-, and K-type to M-type stars. Methods: Using the Lomb-Scargle periodogram we searched more than 150 000 stellar light curves for periodic brightness variations. We analyzed periods between 1 and 30 days in eight consecutive Kepler quarters, where 30 days is an estimated maximum for the validity of the PDC_MAP data correction pipeline. We selected stable rotation periods, i.e., periods that do not vary from the median by more than one day in at least six of the eight quarters. We averaged the periods for each stellar spectral class according to B - V color and compared the results to archival vsini data, using stellar radii estimates from the Kepler Input Catalog. Results: We report on the stable starspot rotation periods of 12 151 Kepler stars. We find good agreement between starspot velocities and vsini data for all F-, G- and early K-type stars. The 795 M-type stars in our sample have a median rotation period of 15.4 days. We find an excess of M-type stars with periods less than 7.5 days that are potentially fast-rotating and fully convective. Measuring photometric variability in multiple Kepler quarters appears to be a straightforward and reliable way to determine the rotation periods of a large sample of active stars, including late-type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asteroseismology of the solar analogs 16 Cyg A & B from Kepler observations

The evolved solar-type stars 16 Cyg A & B have long been studied as solar analogs, yielding a glimpse into the future of our own Sun. The orbital period of the binary system is too long to provide meaningful dynamical constraints on the stellar properties, but asteroseismology can help because the stars are among the brightest in the Kepler field. We present an analysis of three months of nearly uninterrupted photometry of 16 Cyg A & B from the Kepler space telescope. We extract a total of 46 and 41 oscillation frequencies for the two components respectively, including a clear detection of octupole (l=3) modes in both stars. We derive the properties of each star independently using the Asteroseismic Modeling Portal, fitting the individual oscillation frequencies and other observational constraints simultaneously. We evaluate the systematic uncertainties from an ensemble of results generated by a variety of stellar evolution codes and fitting methods. The optimal models derived by fitting each component individually yield a common age (t=6.8+/-0.4 Gyr) and initial composition (Z_i=0.024+/-0.002, Y_i=0.25+/-0.01) within the uncertainties, as expected for the components of a binary system, bolstering our confidence in the reliability of asteroseismic techniques. The longer data sets that will ultimately become available will allow future studies of differential rotation, convection zone depths, and long-term changes due to stellar activity cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Precise Modeling of the Exoplanet Host Star and CoRoT Main Target HD 52265

This paper presents a detailed and precise study of the characteristics of the Exoplanet Host Star and CoRoT main target HD 52265, as derived from asteroseismic studies. The results are compared with previous estimates, with a comprehensive summary and discussion. The basic method is similar to that previously used by the Toulouse group for solar-type stars. Models are computed with various initial chemical compositions and the computed p-mode frequencies are compared with the observed ones. All models include atomic diffusion and the importance of radiative accelerations is discussed. Several tests are used, including the usual frequency combinations and the fits of the échelle diagrams. The possible surface effects are introduced and discussed. Automatic codes are also used to find the best model for this star (SEEK, AMP) and their results are compared with that obtained with the detailed method. We find precise results for the mass, radius and age of this star, as well as its effective temperature and luminosity. We also give an estimate of the initial helium abundance. These results are important for the characterization of the star-planet system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasi full-disk maps of solar horizontal velocities using SDO/HMI data

For the first time, the motion of granules (solar plasma on the surface on scales larger than 2.5 Mm) has been followed over the entire visible surface of the Sun, using SDO/HMI white-light data. Horizontal velocity fields are derived from image correlation tracking using a new version of the coherent structure tracking algorithm.The spatial and temporal resolutions of the horizontal velocity map are 2.5 Mm and 30 min respectively . From this reconstruction, using the multi-resolution analysis, one can obtain to the velocity field at different scales with its derivatives such as the horizontal divergence or the vertical component of the vorticity. The intrinsic error on the velocity is ~0.25 km/s for a time sequence of 30 minutes and a mesh size of 2.5 Mm.This is acceptable compared to the granule velocities, which range between 0.3 km/s and 1.8 km/s. A high correlation between velocities computed from Hinode and SDO/HMI has been found (85%). From the data we derive the power spectrum of the supergranulation horizontal velocity field, the solar differential rotation, and the meridional velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Seismic evidence for a rapidly rotating core in a lower-giant-branch star observed with Kepler

Rotation is expected to have an important influence on the structure and the evolution of stars. However, the mechanisms of angular momentum transport in stars remain theoretically uncertain and very complex to take into account in stellar models. To achieve a better understanding of these processes, we desperately need observational constraints on the internal rotation of stars, which until very recently were restricted to the Sun. In this paper, we report the detection of mixed modes - i.e. modes that behave both as g modes in the core and as p modes in the envelope - in the spectrum of the early red giant KIC7341231, which was observed during one year with the Kepler spacecraft. By performing an analysis of the oscillation spectrum of the star, we show that its non-radial modes are clearly split by stellar rotation and we are able to determine precisely the rotational splittings of 18 modes. We then find a stellar model that reproduces very well the observed atmospheric and seismic properties of the star. We use this model to perform inversions of the internal rotation profile of the star, which enables us to show that the core of the star is rotating at least five times faster than the envelope. This will shed new light on the processes of transport of angular momentum in stars. In particular, this result can be used to place constraints on the angular momentum coupling between the core and the envelope of early red giants, which could help us discriminate between the theories that have been proposed over the last decades.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accurate p-mode measurements of the G0V metal-rich CoRoT target HD 52265

The star HD 52265 is a G0V metal-rich exoplanet-host star observed in the seismology field of the CoRoT space telescope from November 2008 to March 2009. The satellite collected 117 days of high-precision photometric data on this star, showing that it presents solar-like oscillations. HD 52265 was also observed in spectroscopy with the Narval spectrograph at the same epoch. We characterise HD 52265 using both spectroscopic and seismic data. The fundamental stellar parameters of HD 52265 were derived with the semi-automatic software VWA, and the projected rotational velocity was estimated by fitting synthetic profiles to isolated lines in the observed spectrum. The parameters of the observed p modes were determined with a maximum-likelihood estimation. We performed a global fit of the oscillation spectrum, over about ten radial orders, for degrees l=0 to 2. We also derived the properties of the granulation, and analysed a signature of the rotation induced by the photospheric magnetic activity. Precise determinations of fundamental parameters have been obtained: Teff = 6100 +- 60 K, log g = 4.35 +- 0.09, [M/H] = 0.19 +- 0.05, as well as vsini = 3.6 +0.3 -1.0 km/s. We have measured a mean rotation period P_rot = 12.3 +- 0.15 days, and find a signature of differential rotation. The frequencies of 31 modes are reported in the range 1500-2550 micro-Hz. The large separation exhibits a clear modulation around the mean value <Dnu> = 98.3 +- 0.1 micro-Hz. Mode widths vary with frequency along an S-shape with a clear local maximum around 1800 micro-Hz. We deduce lifetimes ranging between 0.5 and 3 days for these modes. Finally, we find a maximal bolometric amplitude of about 3.96 +- 0.24 ppm for radial modes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multi-Channel Three-Dimensional SOLA Inversion for Local Helioseismology

Inversions for local helioseismology are an important and necessary step for obtaining three-dimensional maps of various physical quantities in the solar interior. Frequently, the full inverse problems that one would like to solve prove intractable because of computational constraints. Due to the enormous seismic data sets that already exist and those forthcoming, this is a problem that needs to be addressed. To this end, we present a very efficient linear inversion algorithm for local helioseismology. It is based on a subtractive optimally localized averaging (SOLA) scheme in the Fourier domain, utilizing the horizontal-translation invariance of the sensitivity kernels. In Fourier space the problem decouples into many small problems, one for each horizontal wave vector. This multi-channel SOLA method is demonstrated for an example problem in time-distance helioseismology that is small enough to be solved both in real and Fourier space. We find that both approaches are successful in solving the inverse problem. However, the multi-channel SOLA algorithm is much faster and can easily be parallelized.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Validated helioseismic inversions for 3-D vector flows

According to time-distance helioseismology, information about internal fluid motions is encoded in the travel times of solar waves. The inverse problem consists of inferring 3-D vector flows from a set of travel-time measurements. Here we investigate the potential of time-distance helioseismology to infer 3-D convective velocities in the near-surface layers of the Sun. We developed a new Subtractive Optimally Localised Averaging (SOLA) code suitable for pipeline pseudo-automatic processing. Compared to its predecessor, the code was improved by accounting for additional constraints in order to get the right answer within a given noise level. The main aim of this study is to validate results obtained by our inversion code. We simulate travel-time maps using a snapshot from a numerical simulation of solar convective flows, realistic Born travel-time sensitivity kernels, and a realistic model of travel-time noise. These synthetic travel times are inverted for flows and the results compared with the known input flow field. Additional constraints are implemented in the inversion: cross-talk minimization between flow components and spatial localization of inversion coefficients. Using modes f, p1 through p4, we show that horizontal convective flow velocities can be inferred without bias, at a signal-to-noise ratio greater than one in the top 3.5 Mm, provided that observations span at least four days. The vertical component of velocity (v_z), if it were to be weak, is more difficult to infer and is seriously affected by cross-talk from horizontal velocity components. We emphasise that this cross-talk must be explicitly minimised in order to retrieve v_z in the top 1 Mm. We also show that statistical averaging over many different areas of the Sun allows for reliably measuring of average properties of all three flow components in the top 5.5 Mm of the convection zone.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

3D Numerical Simulations of f-Mode Propagation Through Magnetic Flux Tubes

Three-dimensional numerical simulations have been used to study the scattering of a surface-gravity wave packet by vertical magnetic flux tubes, with radii from 200 km to 3 Mm, embedded in stratified polytropic atmosphere. The scattered wave was found to consist primarily of m=0 (axisymmetric) and m=1 modes. It was found that the ratio of the amplitude of these two modes is strongly dependant on the radius of the flux tube: The kink mode is the dominant mode excited in tubes with a small radius while the sausage mode is dominant for large tubes. Simulations of this type provide a simple, efficient and robust way to start understanding the seismic signature of flux tubes, which have recently began to be observed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An absorbing boundary formulation for the stratified, linearized, ideal MHD equations based on an unsplit, convolutional perfectly matched layer

Perfectly matched layers are a very efficient and accurate way to absorb waves in media. We present a stable convolutional unsplit perfectly matched formulation designed for the linearized stratified Euler equations. However, the technique as applied to the Magneto-hydrodynamic (MHD) equations requires the use of a sponge, which, despite placing the perfectly matched status in question, is still highly efficient at absorbing outgoing waves. We study solutions of the equations in the backdrop of models of linearized wave propagation in the Sun. We test the numerical stability of the schemes by integrating the equations over a large number of wave periods.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constructing semi-empirical sunspot models for helioseismology

One goal of helioseismology is to determine the subsurface structure of sunspots. In order to do so, it is important to understand first the near-surface effects of sunspots on solar waves, which are dominant. Here we construct simplified, cylindrically-symmetric sunspot models, which are designed to capture the magnetic and thermodynamics effects coming from about 500 km below the quiet-Sun $τ_{5000}=1$ level to the lower chromosphere. We use a combination of existing semi-empirical models of sunspot thermodynamic structure (density, temperature, pressure): the umbral model of Maltby et al. (1986) and the penumbral model of Ding and Fang (1989). The OPAL equation of state tables are used to derive the sound speed profile. We smoothly merge the near-surface properties to the quiet-Sun values about 1mm below the surface. The umbral and penumbral radii are free parameters. The magnetic field is added to the thermodynamic structure, without requiring magnetostatic equilibrium. The vertical component of the magnetic field is assumed to have a Gaussian horizontal profile, with a maximum surface field strength fixed by surface observations. The full magnetic field vector is solenoidal and determined by the on-axis vertical field, which, at the surface, is chosen such that the field inclination is 45$^\circ$ at the umbral-penumbral boundary. We construct a particular sunspot model based on SOHO/MDI observations of the sunspot in active region NOAA 9787. The helioseismic signature of the model sunspot is studied using numerical simulations of the propagation of f, p$_1$, and p$_2$ wave packets. These simulations are compared against cross-covariances of the observed wave field. We find that the sunspot model gives a helioseismic signature that is similar to the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convectively stabilised background solar models for local helioseismology

In local helioseismology numerical simulations of wave propagation are useful to model the interaction of solar waves with perturbations to a background solar model. However, the solution to the equations of motions include convective modes that can swamp the waves we are interested in. For this reason, we choose to first stabilise the background solar model against convection by altering the vertical pressure gradient. Here we compare the eigenmodes of our convectively stabilised model with a standard solar model (Model S) and find a good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Helioseismology of Sunspots: A Case Study of NOAA Region 9787

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Linear Sensitivity of Helioseismic Travel Times to Local Flows

Time-distance helioseismology is a technique for measuring the time for waves to travel from one point on the solar surface to another. These wave travel times are affected by advection by subsurface flows. Inferences of plasma flows based on observed travel times depend critically on the ability to accurately model the effects of subsurface flows on time-distance measurements. We present a Born approximation based computation of the sensitivity of time distance travel times to weak, steady, inhomogeneous subsurface flows. Three sensitivity functions are obtained, one for each component of the 3D vector flow. We show that the depth sensitivity of travel times to horizontally uniform flows is given approximately by the kinetic energy density of the oscillation modes which contribute to the travel times. For flows with strong depth dependence, the Born approximation can give substantially different results than the ray approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Local Helioseismology: Three Dimensional Imaging of the Solar Interior

The Sun supports a rich spectrum of internal waves that are continuously excited by turbulent convection. The GONG network and the MDI/SOHO space instrument provide an exceptional data base of spatially-resolved observations of solar oscillations, covering an entire sunspot cycle (11 years). Local helioseismology is a set of tools for probing the solar interior in three dimensions using measurements of wave travel times and local mode frequencies. Local helioseismology has discovered (i) near-surface vector flows associated with convection (ii) 250 m/s subsurface horizontal outflows around sunspots (iii) ~50 m/s extended horizontal flows around active regions (converging near the surface and diverging below), (iv) the effect of the Coriolis force on convective flows and active region flows (v) the subsurface signature of the 15 m/s poleward meridional flow, (vi) a +/-5 m/s time-varying depth-dependent component of the meridional circulation around the mean latitude of activity, and (vii) magnetic activity on the far side of the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Outstanding problems in local helioseismology

Time-distance helioseismology and related techniques show great promise for probing the structure and dynamics of the subphotospheric layers of the Sun. Indeed time-distance helioseismology has already been applied to make inferences about structures and flows under sunspots and active regions, to map long-lived convective flow patterns, and so on. Yet certainly there are still many inadequacies in the current approaches and, as the data get better and the questions we seek to address get more subtle, methods that were previously regarded as adequate are no longer acceptable. Here we give a short and partial description of outstanding problems in local helioseismology, using time-distance helioseismology as a guiding example.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Scattering of helioseismic waves by a sunspot: wavefront healing and folding

We observe and characterize the scattering of acoustic wave packets by a sunspot, in a regime where the wavelength is comparable to the size of the sunspot. Spatial maps of wave traveltimes and amplitudes are measured from the cross-covariance function of the random wave field. The averaging procedure is such that incoming wave packets are plane wave packets. Observations show that the magnitude of the traveltime perturbation caused by the sunspot diminishes as waves propagate away from the sunspot -- a finite-wavelength phenomenon known as wavefront healing. Observations also show a reduction of the amplitude of the waves after their passage through the sunspot. A significant fraction of this amplitude reduction is due to the defocusing of wave energy by the fast wave-speed perturbation introduced by the sunspot. This ``geometrical attenuation'' will contribute to the wave amplitude reduction in addition to the physical absorption of waves. In addition, we observe an enhancement of wave amplitude away from the central path: diffracted rays intersect with unperturbed rays (caustics) and wavefronts fold and triplicate. Thus we find that ray tracing is useful to interpret these phenomena, although it cannot explain wavefront healing.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

SLiM: a code for the simulation of wave propagation through an inhomogeneous, magnetised solar atmosphere

In this paper we describe the semi-spectral linear MHD (SLiM) code which we have written to follow the interaction of linear waves through an inhomogeneous three-dimensional solar atmosphere. The background model allows almost arbitrary perturbations of density, temperature, sound speed as well as magnetic and velocity fields. We give details of several of the tests we have used to check the code. The code will be useful in understanding the helioseismic signatures of various solar features, including sunspots.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Tomography of the Solar Interior

Solar oscillations consist of a rich spectrum of internal acoustic waves and surface gravity waves, stochastically excited by turbulent convection. They have been monitored almost continuously over the last ten years with high-precision Doppler images of the solar surface. The purpose of helioseismology is to retrieve information about the structure and the dynamics of the solar interior from the frequencies, phases, and amplitudes of solar waves. Methods of analysis are being developed to make three-dimensional images of subsurface motions and temperature inhomogeneities in order to study convective structures and regions of magnetic activity, like sunspots.

0 citations0 reviewsNo code linkedOpen access