Helioseismology is the study of the solar interior using observations of oscillations at the surface. It suffers from systematic errors, such as a center-to-limb error in travel-time measurements. Understanding these errors requires a good understand
ing of the nontrivial relationship between wave displacement and helioseismic observables. The wave displacement causes perturbations in the atmospheric thermodynamical quantities which perturb the opacity, the optical depth, the source function, and the local ray geometry, thus affecting the emergent intensity. We aim to establish the most complete relationship up to now between the displacement and the intensity perturbation by solving the radiative transfer problem in the atmosphere. We derive an expression for the intensity perturbation caused by acoustic oscillations at any point on the solar disk by applying the first-order perturbation theory. As input, we consider adiabatic modes of oscillation of different degrees. The background and the perturbed intensities are computed considering the main sources of opacity in the continuum. We find that, for all modes, the perturbations to the thermodynamical quantities are not sufficient to model the intensity. In addition, the geometrical effects due to the displacement must be taken into account as they lead to a difference in amplitude and a phase shift between the temperature at the surface and intensity perturbations. The closer to the limb, the larger the differences. This work presents improvements for the computation of the intensity perturbations, in particular for high-degree modes, and explains differences in intensity computations in earlier works. The phase shifts and amplitude differences between the temperature and intensity perturbations increase towards the limb. This should help to interpret some of the systematic center-to-limb effects observed in local helioseismology.
Retrograde Rossby waves, measured to have significant amplitudes in the Sun, likely have notable implications for various solar phenomena. Rossby waves create small-amplitude, very-low frequency motions (on the order of the rotation rate and lower),
which in turn shift the resonant frequencies and eigenfunctions of the acoustic modes of the Sun. The detection of even azimuthal orders Rossby modes using mode coupling presents additional challenges and prior work therefore only focused on odd orders. Here, we successfully extend the methodology to measure even azimuthal orders as well. We analyze 4 and 8 years of Helioseismic and Magnetic Imager (HMI) data and consider coupling between different-degree acoustic modes (of separations 1 and 3 in harmonic degree). The technique uses couplings between different frequency bins to capture the temporal variability of the Rossby modes. We observe significant power close to the theoretical dispersion relation for sectoral Rossby modes (where the azimuthal order is same as harmonic degree, s = |t|). Our results are consistent with prior measurements of Rossby modes with azimuthal orders over the range t = 4 to 16 with maximum power occurring at mode t = 8. The amplitudes of these modes vary from 1 to 2 m/s. We place an upper bound of 0.2 m/s on the sectoral t = 2 mode, which we do not detect in our measurements. This effort adds credence to the mode-coupling methodology in helioseismology
Context. The frequencies, lifetimes, and eigenfunctions of solar acoustic waves are affected by turbulent convection, which is random in space and in time. Since the correlation time of solar granulation and the periods of acoustic waves ($sim$5 min)
are similar, the medium in which the waves propagate cannot a priori be assumed to be time independent. Aims. We compare various effective-medium solutions with numerical solutions in order to identify the approximations that can be used in helioseismology. For the sake of simplicity, the medium is one dimensional. Methods. We consider the Keller approximation, the second-order Born approximation, and spatial homogenization to obtain theoretical values for the effective wave speed and attenuation (averaged over the realizations of the medium). Numerically, we computed the first and second statistical moments of the wave field over many thousands of realizations of the medium (finite-amplitude sound-speed perturbations are limited to a 30 Mm band and have a zero mean). Results. The effective wave speed is reduced for both the theories and the simulations. The attenuation of the coherent wave field and the wave speed are best described by the Keller theory. The numerical simulations reveal the presence of coda waves, trailing the coherent wave packet. These late arrival waves are due to multiple scattering and are easily seen in the second moment of the wave field. Conclusions. We find that the effective wave speed can be calculated, numerically and theoretically, using a single snapshot of the random medium (frozen medium); however, the attenuation is underestimated in the frozen medium compared to the time-dependent medium. Multiple scattering cannot be ignored when modeling acoustic wave propagation through solar granulation.
Retrograde-propagating waves of vertical vorticity with longitudinal wavenumbers between 3 and 15 have been observed on the Sun with a dispersion relation close to that of classical sectoral Rossby waves. The observed vorticity eigenfunctions are sym
metric in latitude, peak at the equator, switch sign near $20^circ$-$30^circ$, and decrease at higher latitudes. We search for an explanation that takes into account solar latitudinal differential rotation. In the equatorial $beta$ plane, we study the propagation of linear Rossby waves (phase speed $c <0$) in a parabolic zonal shear flow, $U = - overline{U} xi^2<0$, where $overline{U} = 244$ m/s and $xi$ is the sine of latitude. In the inviscid case, the eigenvalue spectrum is real and continuous and the velocity stream functions are singular at the critical latitudes where $U = c$. We add eddy viscosity in the problem to account for wave attenuation. In the viscous case, the stream functions are solution of a fourth-order modified Orr-Sommerfeld equation. Eigenvalues are complex and discrete. For reasonable values of the eddy viscosity corresponding to supergranular scales and above (Reynolds number $100 le Re le 700$), all modes are stable. At fixed longitudinal wavenumber, the least damped mode is a symmetric mode with a real frequency close to that of the classical Rossby mode, which we call the R mode. For $Re approx 300$, the attenuation and the real part of the eigenfunction is in qualitative agreement with the observations (unlike the imaginary part of the eigenfunction, which has a larger amplitude in the model. Conclusion: Each longitudinal wavenumber is associated with a latitudinally symmetric R mode trapped at low latitudes by solar differential rotation. In the viscous model, R modes transport significant angular momentum from the dissipation layers towards the equator.
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 ob
servations, 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.
Radial velocity (RV) measurements are used to search for planets orbiting late-type main-sequence stars and confirm the transiting planets. The most advanced spectrometers are approaching a precision of $sim 10$ cm/s that implies the need to identify
and correct for all possible sources of RV oscillations intrinsic to the star down to this level and possibly beyond. The recent discovery of global-scale equatorial Rossby waves in the Sun, also called r modes, prompted us to investigate their possible signature in stellar RV measurements. R modes are toroidal modes of oscillation whose restoring force is the Coriolis force and propagate in the retrograde direction in a frame that corotates with the star. The solar r modes with azimuthal orders $3 leq m lesssim 15$ were identified unambiguously because of their dispersion relation and their long e-folding lifetimes of hundreds of days. Here we simulate the RV oscillations produced by sectoral r modes with $2 leq m leq 5$ assuming a stellar rotation period of 25.54 days and a maximum amplitude of the surface velocity of each mode of 2 m/s. This amplitude is representative of the solar measurements, except for the $m=2$ mode which has not yet been observed. Sectoral r modes with azimuthal orders $m=2$ and $3$ would produce RV oscillations with amplitudes of 76.4 and 19.6 cm/s and periods of 19.16 and 10.22 days, respectively, for a star with an inclination of the rotation axis $i=60^{circ}$. Therefore, they may produce rather sharp peaks in the Fourier spectrum of the radial velocity time series that could lead to spurious planetary detections. Sectoral r~modes may represent a source of confusion in the case of slowly rotating inactive stars that are preferential targets for RV planet search. The main limitation of the present investigation is the lack of observational constraint on the amplitude of the $m=2$ mode on the Sun.
During the solar magnetic activity cycle the emergence latitudes of sunspots change, leading to the well-known butterfly diagram. This phenomenon is poorly understood for other stars as starspot latitudes are generally unknown. The related changes in
starspot rotation rates caused by latitudinal differential rotation can however be measured. Using the set of 3093 Kepler stars with activity cycles identified by Reinhold et al. (2017), we aim to study the temporal change in starspot rotation rates over magnetic activity cycles, and how this relates to the activity level, mean rotation rate, and effective temperature of the star. We measure the photometric variability as a proxy for the magnetic activity and the spot rotation rate in each quarter over the duration of the Kepler mission. We phase-fold these measurements with the cycle period. We perform averages over stars with comparable mean rotation rates and effective temperature at fixed activity-cycle phases. We detect a clear correlation between the variation of activity level and the variation of the starspot rotation rate. The sign and amplitude of this correlation depends on the mean stellar rotation and, to a lesser extent, on the effective temperature. For slowly rotating stars (with periods between 15-28 days) the starspot rotation rates are clearly anti-correlated with the level of activity during the activity cycles. A transition is seen at periods of 10-15 days, where stars with effective temperature above 4200K instead show positive correlation. Our results can be interpreted in terms of a stellar butterfly diagram, but these appear different from the Suns as the starspot rotation rates are either in phase or anti-phase with the activity level. Alternatively, the activity cycles seen by Kepler are short (around 2.5 years) and may therefore be secondary cycles, perhaps analogous to the solar quasi-biennial oscillations.
Rotational shear in Sun-like stars is thought to be an important ingredient in models of stellar dynamos. Thanks to helioseismology, rotation in the Sun is characterized well, but the interior rotation profiles of other Sun-like stars are not so well
constrained. Until recently, measurements of rotation in Sun-like stars have focused on the mean rotation, but little progress has been made on measuring or even placing limits on differential rotation. Using asteroseismic measurements of rotation we aim to constrain the radial shear in five Sun-like stars observed by the NASA Kepler mission: KIC004914923, KIC005184732, KIC006116048, KIC006933899, and KIC010963065. We used stellar structure models for these five stars from previous works. These models provide the mass density, mode eigenfunctions, and the convection zone depth, which we used to compute the sensitivity kernels for the rotational frequency splitting of the modes. We used these kernels as weights in a parametric model of the stellar rotation profile of each star, where we allowed different rotation rates for the radiative interior and the convective envelope. This parametric model was incorporated into a fit to the oscillation power spectrum of each of the five Kepler stars. This fit included a prior on the rotation of the envelope, estimated from the rotation of surface magnetic activity measured from the photometric variability. The asteroseismic measurements without the application of priors are unable to place meaningful limits on the radial shear. Using a prior on the envelope rotation enables us to constrain the interior rotation rate and thus the radial shear. In the five cases that we studied, the interior rotation rate does not differ from the envelope by more than approximately +/-30%. Uncertainties in the rotational splittings are too large to unambiguously determine the sign of the radial shear.
Accurate modelling of solar-like oscillators requires that modelled mode frequencies are corrected for the systematic shift caused by improper modelling of the near-surface layers, known as the surface effect. ... We investigate how much additional u
ncertainty is introduced to stellar model parameters by our uncertainty about the functional form of the surface effect. At the same time, we test whether any of the parametrizations is significantly better or worse at modelling observed subgiants and low-luminosity red giants. We model six stars observed by Kepler that show clear mixed modes. We fix the input physics of the stellar models and vary the choice of surface correction ... Models using a solar-calibrated power law correction consistently fit the observations more poorly than the other four corrections. Models with the remaining four corrections generally fit ... about equally well, with the combined surface correction by Ball & Gizon perhaps being marginally superior. The fits broadly agree on the model parameters within about the $2sigma$ uncertainties, with discrepancies between the modified Lorentzian and free power law corrections occasionally exceeding the $3sigma$ level. Relative to the best-fitting values, the total uncertainties on the masses, radii and ages of the stars are all less than 2, 1 and 6 per cent, respectively. A solar-calibrated power law ... appears unsuitable for use with more evolved solar-like oscillators. Among the remaining surface corrections, the uncertainty in the model parameters introduced by the surface effects is about twice as large as the uncertainty in the individual fits for these six stars. Though the fits are thus somewhat less certain because of our uncertainty of how to manage the surface effect, these results also demonstrate that it is feasible to model the individual mode frequencies of subgiants and low-luminosity red giants. ...
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 furthe
r 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.
