Helioseismology Explained

Helioseismology, a term coined by Douglas Gough, is the study of the structure and dynamics of the Sun through its oscillations. These are principally caused by sound waves that are continuously driven and damped by convection near the Sun's surface. It is similar to geoseismology, or asteroseismology (also coined by Gough), which are respectively the studies of the Earth or stars through their oscillations. While the Sun's oscillations were first detected in the early 1960s, it was only in the mid-1970s that it was realized that the oscillations propagated throughout the Sun and could allow scientists to study the Sun's deep interior. The modern field is separated into global helioseismology, which studies the Sun's resonant modes directly, and local helioseismology, which studies the propagation of the component waves near the Sun's surface.

Helioseismology has contributed to a number of scientific breakthroughs. The most notable was to show that the anomaly in the predicted neutrino flux from the Sun could not be caused by flaws in stellar models and must instead be a problem of particle physics. The so-called solar neutrino problem was ultimately resolved by neutrino oscillations. The experimental discovery of neutrino oscillations was recognized by the 2015 Nobel Prize for Physics.[1] Helioseismology also allowed accurate measurements of the quadrupole (and higher-order) moments of the Sun's gravitational potential, which are consistent with General Relativity. The first helioseismic calculations of the Sun's internal rotation profile showed a rough separation into a rigidly-rotating core and differentially-rotating envelope. The boundary layer is now known as the tachocline and is thought to be a key component for the solar dynamo. Although it roughly coincides with the base of the solar convection zone — also inferred through helioseismology — it is conceptually distinct, being a boundary layer in which there is a meridional flow connected with the convection zone and driven by the interplay between baroclinicity and Maxwell stresses.

Helioseismology benefits most from continuous monitoring of the Sun, which began first with uninterrupted observations from near the South Pole over the austral summer. In addition, observations over multiple solar cycles have allowed helioseismologists to study changes in the Sun's structure over decades. These studies are made possible by global telescope networks like the Global Oscillations Network Group (GONG) and the Birmingham Solar Oscillations Network (BiSON), which have been operating for over several decades.

Types of solar oscillation

Solar oscillation modes are interpreted as resonant vibrations of a roughly spherically symmetric self-gravitating fluid in hydrostatic equilibrium. Each mode can then be represented approximately as the product of a function of radius

r

and a spherical harmonic
m
Y
l(\theta,\phi)
, and consequently can be characterized by the three quantum numbers which label:

n

;

\ell

; and

m

.

It can be shown that the oscillations are separated into two categories: interior oscillations and a special category of surface oscillations. More specifically, there are:

Pressure modes (p modes)

Pressure modes are in essence standing sound waves. The dominant restoring force is the pressure (rather than buoyancy), hence the name. All the solar oscillations that are used for inferences about the interior are p modes, with frequencies between about 1 and 5 millihertz and angular degrees ranging from zero (purely radial motion) to order

103

. Broadly speaking, their energy densities vary with radius inversely proportional to the sound speed, so their resonant frequencies are determined predominantly by the outer regions of the Sun. Consequently it is difficult to infer from them the structure of the solar core.

Gravity modes (g modes)

Gravity modes are confined to convectively stable regions, either the radiative interior or the atmosphere. The restoring force is predominantly buoyancy, and thus indirectly gravity, from which they take their name. They are evanescent in the convection zone, and therefore interior modes have tiny amplitudes at the surface and are extremely difficult to detect and identify. It has long been recognized that measurement of even just a few g modes could substantially increase our knowledge of the deep interior of the Sun. However, no individual g mode has yet been unambiguously measured, although indirect detections have been both claimed and challenged. Additionally, there can be similar gravity modes confined to the convectively stable atmosphere.

Surface gravity modes (f modes)

Surface gravity waves are analogous to waves in deep water, having the property that the Lagrangian pressure perturbation is essentially zero. They are of high degree

\ell

, penetrating a characteristic distance

R/\ell

, where

R

is the solar radius. To good approximation, they obey the so-called deep-water-wave dispersion law:
2=gk
\omega
\rmh
, irrespective of the stratification of the Sun, where

\omega

is the angular frequency,

g

is the surface gravity and

k\rm=\ell/R

is the horizontal wavenumber, and tend asymptotically to that relation as

k\rminfty

.

What seismology can reveal

The oscillations that have been successfully utilized for seismology are essentially adiabatic. Their dynamics is therefore the action of pressure forces

p

(plus putative Maxwell stresses) against matter with inertia density

\rho

, which itself depends upon the relation between them under adiabatic change, usually quantified via the (first) adiabatic exponent

\gamma1

. The equilibrium values of the variables

p

and

\rho

(together with the dynamically small angular velocity

\Omega

and magnetic field

\rmB

) are related by the constraint of hydrostatic support, which depends upon the total mass

M

and radius

R

of the Sun. Evidently, the oscillation frequencies

\omega

depend only on the seismic variables

\rho(p,\Omega,\rmB)

,

\gamma1

,

\Omega

and

\rmB

, or any independent set of functions of them. Consequently it is only about these variables that information can be derived directly. The square of the adiabatic sound speed,

c2=\gamma1p/\rho

, is such commonly adopted function, because that is the quantity upon which acoustic propagation principally depends. Properties of other, non-seismic, quantities, such as helium abundance,

Y

, or main-sequence age

t\odot

, can be inferred only by supplementation with additional assumptions, which renders the outcome more uncertain.

Data analysis

Global helioseismology

The chief tool for analysing the raw seismic data is the Fourier transform. To good approximation, each mode is a damped harmonic oscillator, for which the power as a function of frequency is a Lorentz function. Spatially resolved data are usually projected onto desired spherical harmonics to obtain time series which are then Fourier transformed. Helioseismologists typically combine the resulting one-dimensional power spectra into a two-dimensional spectrum.

The lower frequency range of the oscillations is dominated by the variations caused by granulation. This must first be filtered out before (or at the same time that) the modes are analysed. Granular flows at the solar surface are mostly horizontal, from the centres of the rising granules to the narrow downdrafts between them. Relative to the oscillations, granulation produces a stronger signal in intensity than line-of-sight velocity, so the latter is preferred for helioseismic observatories.

Local helioseismology

Local helioseismology—a term coined by Charles Lindsey, Doug Braun and Stuart Jefferies in 1993[2] —employs several different analysis methods to make inferences from the observational data.

Inversion

Introduction

The Sun's oscillation modes represent a discrete set of observations that are sensitive to its continuous structure. This allows scientists to formulate inverse problems for the Sun's interior structure and dynamics. Given a reference model of the Sun, the differences between its mode frequencies and those of the Sun, if small, are weighted averages of the differences between the Sun's structure and that of the reference model. The frequency differences can then be used to infer those structural differences. The weighting functions of these averages are known as kernels.

Structure

The first inversions of the Sun's structure were made using Duvall's law and later using Duvall's law linearized about a reference solar model. These results were subsequently supplemented by analyses that linearize the full set of equations describing the stellar oscillations about a theoretical reference model and are now a standard way to invert frequency data. The inversions demonstrated differences in solar models that were greatly reduced by implementing gravitational settling: the gradual separation of heavier elements towards the solar centre (and lighter elements to the surface to replace them).

Rotation

If the Sun were perfectly spherical, the modes with different azimuthal orders m would have the same frequencies. Rotation, however, breaks this degeneracy, and the modes frequencies differ by rotational splittings that are weighted-averages of the angular velocity through the Sun. Different modes are sensitive to different parts of the Sun and, given enough data, these differences can be used to infer the rotation rate throughout the Sun. For example, if the Sun were rotating uniformly throughout, all the p modes would be split by approximately the same amount. Actually, the angular velocity is not uniform, as can be seen at the surface, where the equator rotates faster than the poles. The Sun rotates slowly enough that a spherical, non-rotating model is close enough to reality for deriving the rotational kernels.

Helioseismology has shown that the Sun has a rotation profile with several features:

Relationship to other fields

Geoseismology

See main article: Seismology.

Helioseismology was born from analogy with geoseismology but several important differences remain. First, the Sun lacks a solid surface and therefore cannot support shear waves. From the data analysis perspective, global helioseismology differs from geoseismology by studying only normal modes. Local helioseismology is thus somewhat closer in spirit to geoseismology in the sense that it studies the complete wavefield.

Asteroseismology

See main article: Asteroseismology.

Because the Sun is a star, helioseismology is closely related to the study of oscillations in other stars, known as asteroseismology. Helioseismology is most closely related to the study of stars whose oscillations are also driven and damped by their outer convection zones, known as solar-like oscillators, but the underlying theory is broadly the same for other classes of variable star.

The principal difference is that oscillations in distant stars cannot be resolved. Because the brighter and darker sectors of the spherical harmonic cancel out, this restricts asteroseismology almost entirely to the study of low degree modes (angular degree

\ell\leq3

). This makes inversion much more difficult but upper limits can still be achieved by making more restrictive assumptions.

History

Solar oscillations were first observed in the early 1960s as a quasi-periodic intensity and line-of-sight velocity variation with a period of about 5 minutes. Scientists gradually realized that the oscillations might be global modes of the Sun and predicted that the modes would form clear ridges in two-dimensional power spectra. The ridges were subsequently confirmed in observations of high-degree modes in the mid 1970s, and mode multiplets of different radial orders were distinguished in whole-disc observations. At a similar time, Jørgen Christensen-Dalsgaard and Douglas Gough suggested the potential of using individual mode frequencies to infer the interior structure of the Sun. They calibrated solar models against the low-degree data finding two similarly good fits, one with low

Y

and a corresponding low neutrino production rate

L\nu

, the other with higher

Y

and

L\nu

; earlier envelope calibrations against high-degree frequencies preferred the latter, but the results were not wholly convincing. It was not until Tom Duvall and Jack Harvey connected the two extreme data sets by measuring modes of intermediate degree to establish the quantum numbers associated with the earlier observations that the higher-

Y

model was established, thereby suggesting at that early stage that the resolution of the neutrino problem must lie in nuclear or particle physics.

New methods of inversion developed in the 1980s, allowing researchers to infer the profiles sound speed and, less accurately, density throughout most of the Sun, corroborating the conclusion that residual errors in the inference of the solar structure is not the cause of the neutrino problem. Towards the end of the decade, observations also began to show that the oscillation mode frequencies vary with the Sun's magnetic activity cycle.

To overcome the problem of not being able to observe the Sun at night, several groups had begun to assemble networks of telescopes (e.g. the Birmingham Solar Oscillations Network, or BiSON, and the Global Oscillation Network Group[9]) from which the Sun would always be visible to at least one node. Long, uninterrupted observations brought the field to maturity, and the state of the field was summarized in a 1996 special issue of Science magazine. This coincided with the start of normal operations of the Solar and Heliospheric Observatory (SoHO), which began producing high-quality data for helioseismology.

The subsequent years saw the resolution of the solar neutrino problem, and the long seismic observations began to allow analysis of multiple solar activity cycles. The agreement between standard solar models and helioseismic inversions was disrupted by new measurements of the heavy element content of the solar photosphere based on detailed three-dimensional models. Though the results later shifted back towards the traditional values used in the 1990s, the new abundances significantly worsened the agreement between the models and helioseismic inversions. The cause of the discrepancy remains unsolved and is known as the solar abundance problem.

Space-based observations by SoHO have continued and SoHO was joined in 2010 by the Solar Dynamics Observatory (SDO), which has also been monitoring the Sun continuously since its operations began. In addition, ground-based networks (notably BiSON and GONG) continue to operate, providing nearly continuous data from the ground too.

See also

External links

Satellite instruments

Ground-based instruments

Further reading

Notes and References

  1. News: Neutrino 'flip' wins physics Nobel Prize . BBC News . 6 October 2015 . Webb, Jonathan.
  2. Book: Lindsey, C. . Braun, D.C. . Jefferies, S.M. . "Local Helioseismology of Subsurface Structure" in "GONG 1992. Seismic Investigation of the Sun and Stars" . GONG 1992. Seismic Investigation of the Sun and Stars. Proceedings of a Conference Held in Boulder . January 1993 . Astronomical Society of the Pacific Conference Series . 42 . T.M. Brown . 81–84 . 978-0-937707-61-6 . 1993ASPC...42...81L.
  3. Braun, D.C. . Duvall Jr., T.L. . Labonte, B.J. . Acoustic absorption by sunspots . The Astrophysical Journal . August 1987 . 319 . L27–L31 . 10.1086/184949 . 1987ApJ...319L..27B.
  4. Hill, F. . Rings and trumpets - Three-dimensional power spectra of solar oscillations . Astrophysical Journal . October 1988 . 333 . 996–1013 . 10.1086/166807 . 1988ApJ...333..996H.
  5. Basu, S. . Antia, H.M. . Bogart, R.S. . Ring-Diagram Analysis of the Structure of Solar Active Regions . The Astrophysical Journal . August 2004 . 610 . 2 . 1157–1168 . 10.1086/421843 . 2004ApJ...610.1157B. free .
  6. Duvall Jr., T.L. . Jefferies, S.M. . Harvey, J.W. . Pomerantz, M.A. . 4244835 . Time-distance helioseismology . Nature . April 1993 . 362 . 6419 . 430–432 . 10.1038/362430a0 . 1993Natur.362..430D. 2060/20110005678 . free .
  7. Donea, A.-C. . Braun, D.C. . Lindsey, C. . Seismic Images of a Solar Flare . The Astrophysical Journal . March 1999 . 513 . 2 . L143–L146 . 10.1086/311915 . 1999ApJ...513L.143D. free .
  8. Woodard, M. F. . Solar Subsurface Flow Inferred Directly from Frequency-Wavenumber Correlations in the Seismic Velocity Field . The Astrophysical Journal . January 2002 . 565 . 1 . 634–639 . 10.1086/324546 . 2002ApJ...565..634W. 10.1.1.513.1704 . 122970114 .
  9. )