Hannah Schunker

Physical Properties of Wave Motion in Inclined Magnetic Fields within Sunspot Penumbrae

At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically-propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic waves and the local surface motion within two sunspots. These correlations are modeled by assuming the surface motion to be elliptical, and we explore the properties of the elliptical motion on the magnetic-field inclination. We also demonstrate that the phase shift of the outward-propagating waves is opposite to the phase shift of the inward-propagating waves in stronger, more vertical fields, but similar to the inward phase shifts in weaker, more-inclined fields.

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

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

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

Comparison of observations and simulations provides a strong upper limit on the subsurface rise speed of solar active regions. Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube; it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around emerging active regions, in combination with numerical simulations of solar magnetoconvection, can constrain the subsurface rise speed of emerging magnetic flux. The observed flows imply that the rise speed of the magnetic field is no larger than 150 m/s at a depth of 20 Mm, that is, well below the prediction of the (standard) thin flux tube model but in the range expected for convective velocities at this depth. We conclude that convective flows control the dynamics of rising flux tubes in the upper layers of the Sun and cannot be neglected in models of flux emergence.

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

The content of this chapter has been published in Nielsen et al. (2015), A & A vol. 582, A10. The work was carried out and written by myself, under the supervision of L. Gizon, H. Schunker and J. Schou from the Max Planck Institute for Solar System Research.

Acoustic power absorption and enhancement generated by slow and fast MHD waves - Evidence of solar cycle velocity/intensity amplitude changes consistent with the mode conversion theory

We used long duration, high quality, unresolved (Sun-as-a star) observations collected by the ground based network BiSON and by the instruments GOLF and VIRGO on board the ESA/NASA SOHO satellite to search for solar-cycle-related changes in mode characteristics in velocity and continuum intensity for the frequency range between 2.5mHz < nu < 6.8mHz. Over the ascending phase of solar cycle 23 we found a suppression in the p-mode amplitudes both in the velocity and intensity data between 2.5mHz <nu< 4.5mHz with a maximum suppression for frequencies in the range between 2.5mHz <nu< 3.5mHz. The size of the amplitude suppression is 13+-2 per cent for the velocity and 9+-2 per cent for the intensity observations. Over the range 4.5mHz <nu< 5.5mHz the findings hint within the errors to a null change both in the velocity and intensity amplitudes. At still higher frequencies, in the so called High-frequency Interference Peaks (HIPs) between 5.8mHz <nu < 6.8mHz, we found an enhancement in the velocity amplitudes with the maximum 36+-7 per cent occurring for 6.3mHz <nu< 6.8mHz. However, in intensity observations we found a rather smaller enhancement of about 5+-2 per cent in the same interval. There is evidence that the frequency dependence of solar-cycle velocity amplitude changes is consistent with the theory behind the mode conversion of acoustic waves in a non-vertical magnetic field, but there are some problems with the intensity data, which may be due to the height in the solar atmosphere at which the VIRGO data are taken.

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.

Helioseismology of sunspots: defocusing, folding, and healing of wavefronts

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 travel times and amplitudes are measured from the cross-covariance function of the random wave field observed by SOHO/MDI around the sunspot in active region NOAO 9787. We consider separately incoming plane wave packets consisting of f modes and p modes with radial orders up to four. Observations show that the travel-time perturbations diminish with distance far away from the sunspot – a finite-wavelength phenomenon known as wavefront healing in scattering theory. Observations also show a reduction of the amplitude of the waves after their passage through the sunspot. We suggest that 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 by sunspots. We also observe an enhancement of wave amplitude away from the central path: diffracted rays intersect with unperturbed rays (caustics) and wavefronts fold and triplicate. Wave amplitude measurements in time-distance helioseismology provide independent information that can be used in concert with travel-time measurements.

Surface magnetic field effects in local helioseismology

Using helioseismic holography strong evidence is presented that the phase (or equivalent travel-time) of helioseismic signatures in Dopplergrams within sunspots depend upon the line-of-sight angle in the plane containing the magnetic eld and vertical directions. This is shown for the velocity signal in the penumbrae of two sunspots at 3, 4 and 5 mHz. Phasesensitive holography demonstrates that they are signicantly affected in a strong, moderately inclined magnetic eld. This research indicates that the effects of the surface magnetic eld are potentially very signicant for local helioseismic analysis of active regions.

Limits on radial differential rotation in Sun-like stars from parametric fits to oscillation power spectra

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.

Asteroseismic inversions for radial differential rotation of Sun-like stars: Sensitivity to uncertainties

We quantify the effect of observational spectroscopic and asteroseismic uncertainties on regularised least squares (RLS) inversions for the radial differential rotation of Sun-like and subgiant stars. We first solved the forward problem to model rotational splittings plus the observed uncertainties for models of a Sun-like star, HD 52265, and a subgiant star, KIC 7341231. We randomly perturbed the parameters of the stellar models within the uncertainties of the spectroscopic and asteroseismic constraints and used these perturbed stellar models to compute rotational splittings. We experimented with three rotation profiles: solid body rotation, a step function, and a smooth rotation profile decreasing with radius. We then solved the inverse problem to infer the radial differential rotation profile using a RLS inversion and kernels from the best-fit stellar model.We found that the inversions for Sun-like stars with solar-like radial differential rotation profiles are insensitive to the uncertainties in the stellar models. We found that when the rotation rate below the convection zone is increased to six times that of the surface rotation rate the inferred rotation profile excluded solid body rotation. With the current observational uncertainties, we found that inversions of subgiant stars are sensitive to the uncertainties in the stellar model. Our findings suggest that inversions for the radial differential rotation of subgiant stars would benefit from more tightly constrained stellar models. In Sun-like stars, the insensitivity of the inversions to stellar model uncertainties suggests that it may be possible to perform ensemble inversions for the average radial differential rotation of many stars with a range of stellar types to better constrain the inversions.