Search

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe I 6173.3 {\AA} 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., Comment: accepted for publication in Solar Physics

Quantitative helio- and asteroseismology require very precise measurements of the frequencies, amplitudes, and lifetimes of the global modes of stellar oscillation. It is common knowledge that the precision of these measurements depends on the total length (T), quality, and completeness of the observations. Except in a few simple cases, the effect of gaps in the data on measurement precision is poorly understood, in particular in Fourier space where the convolution of the observable with the observation window introduces correlations between different frequencies. Here we describe and implement a rather general method to retrieve maximum likelihood estimates of the oscillation parameters, taking into account the proper statistics of the observations. Our fitting method applies in complex Fourier space and exploits the phase information. We consider both solar-like stochastic oscillations and long-lived harmonic oscillations, plus random noise. Using numerical simulations, we demonstrate the existence of cases for which our improved fitting method is less biased and has a greater precision than when the frequency correlations are ignored. This is especially true of low signal-to-noise solar-like oscillations. For example, we discuss a case where the precision on the mode frequency estimate is increased by a factor of five, for a duty cycle of 15%. In the case of long-lived sinusoidal oscillations, a proper treatment of the frequency correlations does not provide any significant improvement; nevertheless we confirm that the mode frequency can be measured from gapped data at a much better precision than the 1/T Rayleigh resolution., Comment: Accepted for publication in Solar Physics Topical Issue "Helioseismology, Asteroseismology, and MHD Connections"

Flows in the upper convection zone are measured by helioseismology on a wide variety of scales. These include differential rotation and meridional circulation, local flows around complexes of magnetic activity and sunspots, and convective flows. The temporal evolution of flows through cycle 23 reveals connections between mass motions in the solar interior and the large-scale characteristics of the magnetic cycle. Here I summarize the latest observations and their implications. Observations from local helioseismology suggest that subsurface flows around active regions introduce a solar-cycle variation in the meridional circulation.

The variation of rotation with latitude is poorly known on stars other than the Sun. Several indirect techniques, photometric and spectroscopic, have been used to search for departure from rigid rotation for sufficiently fast rotators. Here we investigate the possibility of measuring stellar differential rotation for solar-type stars through asteroseismology. Rotationally split frequencies of global oscillation provide information about rotation at different latitudes depending on the azimuthal order, m, of the mode of pulsation. We present a method to estimate differential rotation based on the realization that the m =± 1a ndm =± 2 components of quadrupole oscillations can be observed simultaneously in asteroseismology. Rotational frequency splittings can be inverted to provide an estimate of the difference in stellar angular velocity between the equator and 45 ◦ latitude. The precision of the method, assessed through Monte Carlo simulations, depends on the value of the mean rotation and on the inclination angle between the rotation axis and the line of sight.

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 $\tau_{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., Comment: 18 pages, 9 Figures submitted to Solar Physics

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., Comment: 73 pages, 19 figures, accepted by Solar Physics

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

The propagation of solar waves through the sunspot of AR 9787 is observed using temporal cross-correlations of SOHO/MDI Dopplergrams. We then use three-dimensional MHD numerical simulations to compute the propagation of wave packets through self-similar magneto-hydrostatic sunspot models. The simulations are set up in such a way as to allow a comparison with observed cross-covariances (except in the immediate vicinity of the sunspot). We find that the simulation and the f-mode observations are in good agreement when the model sunspot has a peak field strength of 3 kG at the photosphere, less so for lower field strengths. Constraining the sunspot model with helioseismology is only possible because the direct effect of the magnetic field on the waves has been fully taken into account. Our work shows that the full-waveform modeling of sunspots is feasible., Comment: 21 pages, Accepted in Solar Physics