2002 - The Sun and Space Weather - HANSLMEIER
2002 - The Sun and Space Weather - HANSLMEIER
2002 - The Sun and Space Weather - HANSLMEIER
ASTROPHYSICS AND
SPACE SCIENCE LIBRARY
VOLUME 277
EDITORIAL BOARD
Chairman
W.B. BURTON, National Radio Astronomy Observatory, Charlottesville, Virginia, U.S.A.
(burton@starband.net); University of Leiden, The Netherlands (burton@strw.leidenuniv.nl)
Executive Committee
J. M. E. KUIJPERS, Faculty of Science, Nijmegen, The Netherlands
E. P. J. VAN DEN HEUVEL, Astronomical Institute, University of Amsterdam,
The Netherlands
H. VAN DER LAAN, Astronomical Institute, University of Utrecht,
The Netherlands
MEMBERS
I. APPENZELLER, Landessternwarte Heidelberg-Königstuhl, Germany
J. N. BAHCALL, The Institute for Advanced Study, Princeton, U.S.A.
F. BERTOLA, Universitá di Padova, Italy
J. P. CASSINELLI, University of Wisconsin, Madison, U.S.A.
C. J. CESARSKY, Centre d’Etudes de Saclay, Gif-sur-Yvette Cedex, France
O. ENGVOLD, Institute of Theoretical Astrophysics, University of Oslo, Norway
R. McCRAY, University of Colorado, JILA, Boulder, U.S.A.
P. G. MURDIN, Institute of Astronomy, Cambridge, U.K.
F. PACINI, Istituto Astronomia Arcetri, Firenze, Italy
V. RADHAKRISHNAN, Raman Research Institute, Bangalore, India
K. SATO, School of Science, The University of Tokyo, Japan
F. H. SHU, University of California, Berkeley, U.S.A.
B. V. SOMOV, Astronomical Institute, Moscow State University, Russia
R. A. SUNYAEV, Space Research Institute, Moscow, Russia
Y. TANAKA, Institute of Space & Astronautical Science, Kanagawa, Japan
S. TREMAINE, CITA, Princeton University, U.S.A.
N. O. WEISS, University of Cambridge, U.K.
THE SUN AND
SPACE WEATHER
by
ARNOLD HANSLMEIER
Institute for Geophysics,
Astrophysics and Meteorology,
University of Graz, Austria
No part of this eBook may be reproduced or transmitted in any form or by any means, electronic,
mechanical, recording, or otherwise, without written consent from the Publisher
Dear Reader,
Kluwer Academic Publishers is the proud publisher of ASSL Volume 277 “The Sun
and Space Weather” written by Prof. Arnold Hanslmeier of the University of Graz.
However, we must point out that a number of important references and acknowledgements
to work done by the Australian Space Weather Agency, IPS Radio and Space Services
were not included.
In particular:
1) Figures 9.5 and 9.6 are courtesy of IPS Radio and Space Services.
2) The atmospheric model (for orbital decay) in Section 9.2.7 has been developed
by IPS Principal Physicist John Kennewell.
3) The IPS Radio and Space Services Web site www.ips.gov.au should have been
listed in Section 13.2 as an important resource for information. Sections 9.2.4 to
9.2.7 present information that was taken from the educational section of this
Web site.
V
VI CONTENTS
It is generally accepted that the term space weather refers to the time-variable
conditions in the space environment that may effect space-borne or ground based
technological systems and, in the worst case, endanger human health or life. There-
fore there are social and economic aspects of this type of research: one tries to
avoid consequences of space weather events by system design or efficient warning
and prediction. During the last few years space weather activities have expanded
world-wide.
Space weather affects spacecrafts as well as ground based systems.
The main cause for space weather effects is our Sun. It emits a continuous
stream of particles which is called the solar wind. The solar wind is variable. It is
modulated by the well known 11 year solar activity cycle. Another source of space
weather effects are micrometeorites and other space debris.
Related to the solar activity are important effects on spacecraft such as space-
craft charging (surface charging and deep discharges) and single event effects. The
effects on humans in space are also to be considered (radiation, particles). Space
weather effects also play a role on high altitude/high latitude air-flight; cosmic rays
penetrate to the lower atmosphere and pose problems to humans and electronic
components of modern aeroplanes. Other influences of space weather include radio
wave propagation, satellite-ground communications, global satellite-based naviga-
tion systems, power transmission systems etc. Changes of the solar irradiance may
be one of the causes for climatic changes on the Earth.
Space debris, such as meteoroids, or parts of old satellites must be also be
taken into account and are a permanent threat for space missions.
The book is organized as follows. First we want to give a brief review about the
main source of space weather effects, our Sun. The basic physics of the Sun will
be discussed since it is essential to understand the mechanisms that cause solar
variability. This is necessary in order to make prediction models for space weather
forecasts. Then we will speak about the influence of solar variability on the Earth’s
atmosphere. There are similarities with atmospheric weather, however the most
1
2 CHAPTER 1. INTRODUCTION, WHAT IS SPACE WEATHER?
important differences between the atmospheric and space weather systems are:
Meteorological processes are localized; it is possible to make good local
Space weather events occur over a wide range of time scales: the Earth’s
magnetosphere responds to solar-originated disturbances within only a few
minutes, global reconfiguration occurs within some 10 minutes. Enhanced
fluxes of energetic particles in radiation belts decay in time scales of days,
months or even longer.
Spaceweather predictions must rely on the input of just a few isolated mea-
surements of the solar wind and the observations (both ground based and
from space) have only a global character sometimes without details.
Therefore, successful space weather activities aiming to make prediction of dan-
gerous events need to be performed on a global scale. Space-borne and ground
based observations are complementary.
Presently, the most important users of space weather research are spacecraft
engineering, spacecraft operations, RF communications. Spacecraft launchers can
make use of exact knowledge of space weather conditions and the re-entry of space-
crafts depends on the atmospheric drag conditions. When the International Space
Station is in operation, forecasts will become even more important. Other users are
telecommunication operators, users of global positioning systems, electric power
industry etc. Commercial airlines must be careful with the radiation doses to their
crews and passengers.
Our Sun is the only star which is close enough to observe details on its surface.
We can observe sunspots, faculae, prominences, coronal holes etc., which are all
summarized as solar activity phenomena. Therefore, the study of the Sun is im-
portant for astrophysics in general. Theories about stellar structure and evolution
can be studied in detail on the Sun.
On the other hand, the Sun is the driving factor for the climate on the Earth
and the structure of the Earth’s magnetosphere thus determining and influenc-
ing the near Earth space environment. Therefore, the study of solar terrestrial
relations is of great importance for our modern telecommunication systems both
based on Earth and in space. The Sun provides the main input for the so called
space weather. Other inputs are meteorites, space debris etc.
3
4 CHAPTER 2. THE SUN AS A STAR
Many galaxies appear as spiral galaxies. Young bright stars are in the spiral
arms, older stars in the center and in the halo of the galaxy. An example is given
in Fig.2.1.
Thus, from the astrophysical point of view, our sun is situated at about 2/3
from the center of the galaxy and a normal star on a bulk of stars.
observe a binary system. Stellar masses are very critical for stellar evolution,
however we know accurate masses only for some 100 stars.
Once mass and radius are known, the density and the gravitational acceler-
ation follow. These parameters are important for the stellar structure.
Stellar rotation: For simplicity we can assume that a star consists of two
halves, one half approaches to the observer and the spectral lines from that
region are blueshifted, the other half moves away and the spectral lines from
that area are redshifted. The line profile we observe in a spectrum is a
superposition of all these blue- and redshifted profiles and and the rotation
causes a broadening of spectral lines;
Stellar magnetic fields: as it will be discussed in more detail when con-
sidering the Sun, magnetically sensitive spectral lines are split into several
components when there are strong magnetic fields.
is the distance of the object in pc. For the Sun we have from a
distance of 10 pc it would be among the fainter stars visible with the naked eye.
How can we determine stellar temperatures? Stars can be considered to a
very good approximation as black body radiators. A black body is a theoretical
idealization: an object that absorbs completely all radiation at all wavelengths.
The radiation of a black body at a given temperature is given by the Planck law:
6 CHAPTER 2. THE SUN AS A STAR
thus it depends only on the temperature of the object. Here, is the intensity
of radiation at frequency are Planck’s constant, Boltzmann’s constant and
the speed of light. If that equation
is integrated over all frequencies (wavelengths), we obtain a formula for the total
power emitted by a black body:
What is the power emitted per unit area of the Sun’s surface? Answer: Put T =
6 000 K we find that the Sun radiates of its surface.
By taking the derivative with respect to of Planck’s Law and setting it equal
to zero, one can find the peak wavelength, where the intensity is at maximum:
At about what wavelength can planets be expected to radiate most of their energy?
2.1. THE SUN AND THE GALAXY 7
The Sun has a surface temperature of about 6 000 K. At what wavelength does
the Sun’s spectrum peak? Answer:
Therefore, we can measure from the spectrum directly the temperature of stars.
The temperature derived from the peak wavelength is called Wien Temperature,
the temperature derived from the difference between two intensity levels (=color)
Color temperature etc. In order to define color, a filter system must be defined.
The most commonly used system is the UBV system which has three bands that
are located in the UV (U), blue (B) and visual (V) to measure the intensity
The luminosity of stars is given in magnitudes which are defined as follows:
The color of a star is measured by comparing its magnitude through one filter (e.g.
red) with its magnitude through another (e.g. blue).
E.g. means the magnitude measured with the V filter.
Therefore, instead of determining temperatures from the comparison of the
spectrum of a star with the Planck law, one can use e.g. color indices: Let us
compare a hot star with a cool star (Betelgeuze). If we calculate B-V,
than this value will be:
positive for the cooler star, since it is brighter in V than in B (blue). If the
cool star is brighter in V it means that its magnitude has a lower value and
therefore B-V is positive.
negative for the hotter star. The hotter star is brighter in B than in V,
therefore for the magnitudes in these two bands: and B-V<0.
The HRD consists of:
8 CHAPTER 2. THE SUN AS A STAR
Main sequence stars: they lie on a diagonal from the upper left (hot) to the
lower right (cool).
giants, supergiants: they have the same temperature as the corresponding
main sequence stars but are much brighter and must have larger diameters
(see equation 2.4).
White dwarfs are faint but very hot objects thus from their location at the
lower left in the HRD it follows that they must be very compact (about
1/100 the size of the Sun). In about years our Sun will evolve into a
red giant and finally to a white dwarf.
The evolutionary path of our Sun is also given in the 2.2. We have the following
main evolutionary steps:
Pre main sequence evolution: from a protostellar gas and dust cloud the Sun
was formed and before it reaches the main sequence where it spends most of
its life, the contracting Sun passes a violent youth, the T Tauri phase.
At the main sequence the Sun changes extremely slowly and remains about
years there. In the core H is transformed to He by nuclear fusion.
The Sun evolves to a red giant, it will expand and the Earth will become part
of the solar atmosphere. The expansion starts when all H is transformed to
He in the core. Then a H burning shell supplies the energy. The He flash
sets in as soon as in center the He burning sets in. The Sun will be a red
giant for some
Finally the Sun becomes a white dwarf which slowly cools.
In terms of radius we have:
present Sun:
red giant:
white dwarf:
For the main sequence stars there exists a relation between their mass and
luminosity:
From 2.9 we see that more massive stars are very luminous and therefore they
use up their nuclear fuel much more rapidly than low massive stars like our Sun.
Massive main sequence stars that are observed today must have been formed in
very recent astronomical history.
According to their spectra, stars can be classified in the following sequence:
O-B-A-F-G-K-M. This is a sequence of temperature (see Table 2.4): O stars are
hottest, M stars coolest; the number of absorption lines increase from O to M.
Some characteristics are given in Table 2.3.
The luminosity of a star depends on a) temperature, b) surface which is
Since e. g. a K star may be a dwarf main sequence star or a giant luminosity
classes have been introduced. Class I contains the most luminous supergiants,
class II the less luminous supergiants. Class III are the normal giants, Class IV
the sub giants and class V the main sequence.
Now we understand the spectral classification of our Sun: it is a G2V star.
An application of Kepler’s third law gives us the mass of the Sun if its distance is
known which again can be derived from Kepler’s third law:
If we know the distance of the Sun and its angular diameter the solar radius is
obtained:
In order to determine it, we must use the Sun’s distance and make a measurement
of the angular diameter of the visible solar disk. This is not quite easy; one possi-
bility is to define the angular distance between the inflection points of the intensity
profiles at two opposite limbs. Such profiles can be obtained photoelectrically and
the apparent semidiameter at mean solar distance is about 960 seconds of arc (").
Knowing the mass and radius of the Sun, the mean density can be calculated:
The solar constant is the energy crossing unit area of the Earth’s surface per-
pendicular to the direction from the Earth to the Sun in every second. In SI the
units are UV and IR radiation from the Sun is strongly absorbed by the
Earth’s atmosphere. Therefore, accurate measurements of the solar constant are
only possible with satellites. ACRIM on SMM and ERB on Nimbus 7 showed
clearly that the presence of several large sunspots which are cooler than their
surroundings depress the solar luminosity by ~ 0.1%. The solar luminosity is:
If the acceleration is set to zero (when there is a balance), then the hydrostatic
equilibrium becomes:
Therefore, the pressure at depth must be high enough to support the weight of
the fluid per unit area above that depth. Let us derive an estimate for the central
pressure of a star. The pressure is given by:
What happens if a star contracts (which will be the case when the hydrostatic
equilibrium condition is not established)? According to the Virial Theorem half
of the gravitational energy which is set free is radiated away and the other half
heats the star.
12 CHAPTER 2. THE SUN AS A STAR
In these equations is the distance from the stellar center, P, T are the pressure,
density and temperature at radius M is the mass contained within the
energy carried by radiation across the nuclear energy release. The quantities
P, depend on density, temperature and composition. is the opacity and
measures the resistance of the material to energy transport.
If one assumes that the Sun consists of pure hydrogen which is converted into
He, then the total energy would be The luminosity of the Sun is
thus there would be energy supply for years.
The pp reaction dominates for stars with central temperatures between 5 and
convection occurs. For the temperature gradient we already have seen that:
If the stellar material is an ideal classical gas with constant ratio of specific heats
then:
The theory which is usually used contains a free parameter, the so called mixing
length l:
17
18 CHAPTER 3. OBSERVING THE SUN
ground and produces convection in the air just above ground. Turbulent motions
in the air near the observatory are also reduced by the smooth flow of wind across
the lake instead of turbulent flow that occurs over mountain peaks and forests.
The main instrument is a 65 cm reflector.
In Fig. 3.2 a drawing of a solar vacuum tower telescope is given. Light enters
the vacuum tank through a coelostat system and a mirror. The vertical tank is
evacuated in order to avoid turbulence in the telescope itself. At NSO, Kitt Peak,
the telescope is a 70 cm f/52 system.
The German Vacuum Tower Telescope (VTT) at the Observatorio del Teide,
Tenerife has two coelostat mirrors (80 cm) and the entrance window to the vacuum
tank (BK7) has a diameter of 75 cm and a thickness of 7 cm. The primary mirror
has 70 cm and the focal length of the system is 45.64 m. The total field of view is
700 arcsec and the scale is 4.52 arcsec/mm.
Other famous solar instruments for observing the Sun in high spatial resolution
mode are the Coupole at the Observatoire Pic du Midi and the Swedish La Palma
Solar Telescope.
A recent review about solar instrumentation was given by v.d. Lühe (2001).
3.3 SOHO
The SOlar and Heliospheric Observatory is a common project being carried out
by the European Space Agency (ESA) and the US National Aeronautics and Space
Administration (NASA) in the framework of the Solar Terrestrial Science Program
(STSP) comprising other missions like CLUSTER and the International Solar
Terrestrial Physics Program (ISTP) with Geotail, WIND and Polar. SOHO was
launched on December 2, 1995.
SOHO is located at the Lagrangian point L1 about 1.5 Million km away from
Earth which permits an uninterupted view of the Sun. All previous space solar ob-
3.3. SOHO 19
20 CHAPTER 3. OBSERVING THE SUN
servatories have orbited the Earth, from where their observations were periodically
interrupted as our planet ’eclipsed’ the Sun.
The main scientific objectives of SOHO are:
Interior of the Sun
Solar atmosphere
Solar wind
The main instruments are:
CDS: Coronal Diagnostic Spectrometer,
CELIAS: Charge, Element, and Isotope Analysis System,
COSTEP: Comprehensive Suprathermal and Energetic Particle Analyzer,
EIT: Extreme ultraviolet Imaging Telescope,
ERNE: Energetic and Relativistic Nuclei and Electron,
GOLF: Global Oscillations at Low Frequencies,
LASCO: Large Angle and Spectrometric Coronagraph Experiment,
MDI: The Michelson Doppler Imager,
SUMER: Solar Ultraviolet Measurements of Emitted Radiation
SWAN: Solar Wind A Nisotropies,
UVCS: Ultra Violet Coronagraph Spectrometer,
VIRGO: Variability of solar IRradiance and Gravity Oscillations.
Further details about SOHO and the instruments can be found e.g. in the review
given by Fleck (2001).
A drawing of the spacecraft is given in Fig. 3.3.
I intensity,
Usually the polarization parameters are given relative to the intensity, i.e.
U/I, and V/I.
A photodetector can be used to measure the intensity of light generating an
electric current proportional to the intensity I. The polarization parameters how-
ever cannot be measured directly. They have to be measured by using optical de-
vices which are polarization sensitive (linear polarizers, retarders). In a polarime-
ter, the Stokes vector is transformed into
where the transmitted intensity I' depends on
Let us consider a simple example: The intensity I of linearly polarized light is
measured with a photodetector. After introducing a linear polarizer into the light
beam with arbitrary orientation the measured intensity becomes I' < I.
22 CHAPTER 3. OBSERVING THE SUN
The light is unpolarized if for all orientations of the polarizer the intensity is
the same. The orientation of the linear polarizer is varied and thus the intensity
of the light measured with the photodetector is modulated.
Because of seeing effects, for precise polarization measurements the modulation
frequency must lie above the frequency of the intensity fluctuations caused by
disturbances e.g. atmospheric turbulences in solar observations.
The action of linear optical systems on polarized light is described by the
Mueller matrices:
The intensity of the outcoming beam is a linear combination of all four Stokes
parameters of the incoming beam. The Stokes parameters must verify the following
conditions:
When doing solar polarimetry one has to take into account for different effects:
The two basic effects that can be used for measuring magnetic fields are:
Hanle effect: useful diagnostic where the magnetic field is relatively weak
(a few to a few tens of Gauss) and where the plasma is sufficiently tenuous
that collisional excitation can be neglected in comparison to the radiative
excitation of the upper level. It introduces both a rotation of the plane of
polarization and a reduction of the net polarization of the scattered light.
possible for
where is the number of electrons per unit volume and the plasma frequency.
Let us consider the propagation of an electromagnetic wave throughout the Sun’s
outer layers. Here, the density decreases from the chromosphere to the corona.
Therefore, electromagnetic waves with higher frequencies originate in deeper layers
and by observing the Sun in different frequency channels one can measure the
propagation of a plasma wave through the atmosphere (see Fig. 3.4).
Possible emission mechanisms at cm- and mm-wavelengths are for the active
Sun component:
The quiet Sun component of the radio emission is from thermal emission from
the hot ionized gas.
At a frequency of 100 GHz (0.3 cm), the emission originates at the same height
in the photosphere than at visible wavelengths. At 1.4 GHz (21 cm) the emission
originates from the top of the chromosphere (corresponding to a black body at
100 000 K) and at longer wavelengths (e.g. 300 cm, corresponding to 0.1 GHz)
the emission arises from the corona (1-2 Million K blackbody). Thus, the size of
the Sun varies when measured at different wavelengths.
The slowly varying component also has thermal origin and arises from regions
above sunspots where the electron density is higher
Chapter 4
Phenomenology of Solar
Active Regions
4.1 Introduction
The Sun and its atmosphere can be divided as follows:
1. Core: about 1/3 of the solar radius; here energy production occurs.
2. Radiation zone: about 1/3 of the solar radius; the energy is transported
outward by innumerable emission and absorption processes transferring
the high energy photons that are produced by nuclear fusion into
longer wave photons.
3. Convection zone: starts below the surface extending about km
into the interior.
25
26 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
and hence
The Planck function must increase with depth, since the temperature increases
with depth (se Fig. 4.1). Eddington made the following ansatz:
If we look at the solar disk we immediately see that the central regions are
brighter than the limb. The function is called the limb darkening
(center to limb variation). This may be written as:
If one does an inversion of this equation information about the physical struc-
ture (temperature distribution) of the solar atmosphere is obtained. Stellar limb
functions can not be measured accurately so this method is only applicable to the
Sun.
4.2.2 Granulation
Under very good seeing conditions the Sun shows a cellular like pattern which
is called granulation. The mean diameter of the cells is about 1000 km which
corresponds roughly to 1 arcsec (as seen from the Earth). In the bright granules
matter is streaming upwards, in the darker intergranular lanes streaming down-
wards. Up to now the best granulation images have been taken from the ground
since no large solar telescopes have been launched. In order to minimize the effect
of the turbulence of the Earth’s atmosphere (seeing), the exposure times must
be shorter than 1/10 s. Usually, one makes a burst of several images and then
selects the best image for further analysis. Spectrograms show a high degree of
correlation between intensities and velocities proving the convective character of
the phenomenon. Under a spatial resolution better than 0.5 arcsec, the situation
becomes more complex. Regular granules seem to have a maximum for the upflow
near their center, so called exploding granules have a maximum upflow between
the center and the edge. Measuring the width of spectral lines one gets a hint
28 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
for turbulence. Enhanced line widths indicate enhanced turbulence. It was found
that turbulence is located in the downdrafts which is also predicted by 3 D models.
The turbulence may be generated by the shear between upflows and downflows at
granular borders and on transonic flows.
A recent review about solar granulation was given by Muller (1999) where fur-
ther references can be found. A problem to investigate the granulation is how can
we identify a granulum? One possibility is to identify them by an isophote contour
at a level close to the average intensity of the photosphere. The images must be
filtered in order to remove the intensity fluctuations at low frequency, originating
in instrumental brightness inhomogeneities and in solar large scale fluctuations
(which arise from the supergranulation, mesogranulation and oscillations). Fi-
nally, high frequency noise must be eliminated. In the Fourier domain such a filter
has the form:
The parameters are chosen, so that the maximum filter transmission stays in
between spatial scales 0.5 and 1 arcsec. Such a filter is partially restoring as it
enhances the contrast of the smallest granules which can then be identified more
clearly. Another method is to find the inflection points of the intensity distribution
in the image using a Laplacian operator.
How do granules evolve? The most common process is that of fragmentation: a
granule grows and then splits into several fragments (3...4). About 60% of granules
appear or die by this process. Some granules appear spontaneously in intergranular
spaces and grow, others result from merging of two adjacent granules. The most
spectacular evolution is observed for exploding granules.
From the physical point of view, there exists a limitation for the horizontal
expansion because of mass conservation and radiative loss. Matter is streaming
upward in a granulum, expands and horizontal flows are driven by pressure gradi-
ents; thus the central upflow is decelerated which then cannot supply the horizontal
expansion and the radiative loss. The central part cools and the granule splits into
several fragments, after a downdraft developed. On the other hand, intergranu-
lar lanes are interconnected without interruption. They contain some dark holes
which exist over 45 min and may correspond to the fingers of downflowing material
predicted by 3 D models.
Using time series with the 50 cm refractor at the turret dome of the Pic du
Midi observatory Roudier et al. (1997) showed the existence of singularities in
the intergranular lanes what they called intergranular holes which have diameters
between 0.24 arcsec and 0.45 arcsec and are visible fpr more than 45 min. These
holes appear to be systematically distributed at the periphery of mesogranular
and supergranular cells. Spectroscopic observations of the solar granulation with
high resolution yield information about velocities e.g. when observed near solar
disk center, granular profiles are blueshifted because matter rises and moves in
direction to the observer (see Fig. 4.2.2).
The granule lifetime can be determined by their visual identification on suc-
cessive images or by cross correlating these images. There is a large discrepancy
of the results: granular lifetimes range from 6 to 16 min.
4.2. PHENOMENA IN THE SOLAR PHOTOSPHERE 29
where D is the fractal dimension. It seems that there are two ranges with different
fractal dimensions:
continuum is at 1.0). The rms intensity fluctuations of the best image is 10-11%
at (50 cm refractor at La Palma) and 8-9% at (50 cm refractor
at the Pic du Midi). Restored values lie between 10 and 22%. Prom the granular
contrast we can infer the temperature variations (assuming Planck’s law) which
correspond to ~ 200 K.
Theoretical approaches
The simplest model of convection is the classical Rayleigh problem: suppose that
we have a fluid (either gaseous or liquid), confined between two horizontal plates
separated by a distance and maintained at temperature (upper) and (lower)
with If the fluid has a positive coefficient of thermal expansion as it
will be the case for a gas and for a normal fluid, the fluid near the lower plate will
tend to rise. However, this will be opposed by two effects: a) viscous dissipation,
b) thermal diffusion in the fluid. Convection will occur when the imposed temper-
ature gradient is sufficiently large or, for a given gradient, when the
coefficients of the kinematic viscosity v and of thermal diffusion are sufficiently
small. Rayleigh’s theoretical analysis of the problem in 1916 inspired Bénard to
investigate this 40 years later. It was found that convective instability occurs when
the Rayleigh number R exceeds a critical value:
where is the temperature gradient. For Rayleigh found the value 657.5.
This value depends on the boundary conditions. Later Chandrasekhar has shown
that e.g. a Coriolis force (as an effect of rotation) inhibits the onset of instability
to an extend which depends on the value of a non dimensional parameter (called
Taylor number):
here, is the vertical component of the angular velocity vector. For details see
e.g. Chandrasekhar (1961).
Let us discuss the Rayleigh number for the solar convection zone. The value is
extremely high,
Important information about the origin of the solar granulation can be in-
ferred from power spectra. Prom spectrograms we can obtain 1-D power spectra
of intensity and velocity fluctuations, from white light images, one gets 2-D power
spectra for the intensity fluctuations. The theoretical power spectrum of the ve-
locity fluctuations decreases as down to the scale of molecular diffusion. The
temperature power spectrum however decreases as only to a scale At
smaller scales the spectrum decreases as Thus separates the inertial
convective range, where heat advection dominates from the inertial conductive
range, where diffusion dominates. The former is the range of large granules, the
latter the range of small granules.
The basic set of hydrodynamic equations to describe solar convection is as
follows:
4.2. PHENOMENA IN THE SOLAR PHOTOSPHERE 31
Conservation of mass
Conservation of momentum:
Conservation of energy:
Granulation-Mesogranulation
Idealized numerical experiments on turbulent convection were made by Cattaneo et
al. (2001). The authors found two distinct cellular patterns at the surface. Energy-
transporting convection cells (corresponding to granules in the solar photosphere)
have diameters comparable to the layer depth, while macrocells (corresponding to
mesogranules) are several times larger. The motion acts as a small-scale turbulent
32 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
Other authors like e.g. Hoekzema which analyzed G band images found that
photospheric 5 min oscillations are global and rather insensitive to local fine struc-
ture (1998).
A recent paper on that topic was given by Khomenko et al.(2001) where further
references can be found.
4.2.4 Sunspots
Discovery of sunspots
When the Sun is very low just above the horizon one can make a short glimpse on it
with the unprotected naked eye. Chinese astronomers were the first who reported
on dark spots visible on the Sun. In the year 1611 sunspots were observed for the
first time through a telescope and four men can be named as the discoverers: J.
Goldsmid (Holland), G. Galilei (Italy), Ch. Scheiner (Germany) and Th. Harriot
(England). The first publication on that topic appeared from Goldsmid (he is
better known by his Latin name Fabricius). He even argued that the Sun must
rotate since the sunspots move across the disk. Since he was a Jesuit he first
suspected some defect in his telescope when he observed the spots. Then he
failed to persuade his ecclesiastical superiors who refused to allow him to publish
his discovery. However, Scheiner announced his discovery in three anonymous
letters to a friend of Galileo and Galileo responded in three letters in 1612 (the
sunspot letters) that he had discovered the sunspots. Of course Scheiner and
Galileo became enemies. Scheiner later reported his discoveries in his work Rosa
Ursinae sive Sol in 1630. Both scientists noted that the spots appear only within
zones of low latitudes at either side of the equator. There are never spots near the
poles.
Of course one never should risk a blinding of the eyes and the safest way to
observe the Sun with a telescope is when the solar image at the ocular (exit pupil)
is projected onto a screen.
After the initial interest and the publication of Scheiner’s major work there
was no big interest in sunspots. In 1977 Eddy showed that this must be seen in
connection with the fact that during 1640-1705 there was a great reduction in the
number of sunspots seen on the Sun which is now know as the Maunder Minimum.
The next significant discovery was made by Schwabe who was a German
apothecary and bought a telescope in 1826 in order to search for a planet in-
side the orbit of Mercury. He recorded the occurrence of sunspots over 43 years
and reported on a periodicity of their occurrence of about 10 years. In 1851 ap-
peared his publication on the 11 year periodicity of the annually averaged sunspot
numbers. Several years later Carrington showed from his observations that the
Sun rotates differentially; a point at the equator rotates more rapidly than one at
higher latitudes. He defined an arbitrary reference point on latitude 10° as longi-
tude zero and a rotation completed by this point is known as Carrington rotation
(CR). He was also the first to see a white light flare on the Sun in Sep. 1859
during sketching sunspot projections with a friend. Suddenly two crescent-shaped
patches broke out, brightened, moved a distance twice their length, then faded
34 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
away as two dots. All that happened within a period over five minutes. Carring-
ton reported to the Royal Astronomical Society that at 4 hours after midnight
the magnetic instruments indicated a great magnetic storm. So he was in fact
the first who noticed that there exists a connection between solar phenomena and
disturbances on Earth.
R. Wolf of Bern (1816-1893) studied all available records and derived a more
accurate estimate for the sunspot cycle. In 1848 he introduced the relative (Zurich)
sunspot number as a measure for solar activity. Sunspot often appear as
groups. If denotes the number of sunspot groups and the number of individual
spots, then
energy but the magnetic fields destroys this symmetry and one has a splitting of
the energy levels. The displacement of the lines due to the Zeeman effect is given
by:
The wavelength is given in nm, the Landéfactor depends on the spin and
orbital momentum of the levels and B denotes the magnetic induction given in
Tesla.
The strength of the magnetic field is in the order of 3 000 Gauss.
Small dark spots with diameters < 2500 km lacking penumbrae are called
pores. They exist within groups or appear also as isolated structures. Their
lifetimes are in the range of a few hours to several days.
Sunspot groups tend to emerge either sequentially at the same or similar Car-
rington longitudes, which are designated as active longitudes, or to overlap in
clusters.
induction is 0.3 T in the umbra and 0.15 T in the penumbra. In the umbra the field
is approximately vertical, and the inclination increases through the penumbra.
Hale’s observations also suggested that the Sun has an overall dipolar mag-
netic field This very weak dipolar field is reversed over the magnetic
cycle. Almost all of the photospheric field outside sunspots is concentrated in
small magnetic elements with a magnetic induction between 0.1 and 0.15 T.
Only the surface properties of the flux tube that defines a spot can be observed.
The question is, how the field structure changes with depth. The simplest model is
a monolithic column of flux. Let us assume that the pressure inside the flux tube
is negligible compared to the magnetic pressure. We also assume that the grav-
itational force is unimportant in obtaining an approximate idea of the magnetic
field structure, the magnetic field in cylindrical polar coordinates can be taken to
be current free:
Thus Since
The neighboring photosphere, in which the flux tube is embedded has a known
pressure variation with height The boundary of the flux tube is at
where
Here denotes the number of sunspot groups and the number of spots. The
factor is a correction which takes into account for the different instruments used
for the determination of R. In order to smear out effects of solar rotation, R
is given as a monthly averaged number and called the sunspot relative number.
Today there exist better methods to quantify the solar activity however sunspot
numbers are available for nearly 400 years and thus this number is still used.
The Royal Greenwich Observatory (RGO) compiled sunspot observations from
a small network of observatories to produce a data set of daily observations start-
ing in May of 1874. The observatory concluded this data set in 1976 after the
US Air Force (USAF) started compiling data from its own Solar Optical Observ-
ing Network (SOON). This work was continued with the help of the US National
Oceanic and Atmospheric Administration (NOAA) with much of the same infor-
mation being compiled through to the present.
Since 1981, the Royal Observatory of Belgium harbors the Sunspot Index Data
center (SIDC), the World data center for the Sunspot Index. Recently, the Space
Weather forecast center of Paris-Meudon was transferred and added to the ac-
tivities of the SIDC. Moreover, a complete archive of all images of the SOHO
instrument EIT has become available at the SIDC.
Let us briefly summarize the behavior of sunspots during the activity cycle:
The leader spots (i.e. by convention it is defined that the Sun rotates from
east to west; the largest spot of a group tends to be found on the western
side and is called the leader, while the second largest in a group is called
the follower) in each hemisphere are generally all of one polarity, while the
follower spots are of the opposite polarity.
40 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
4.2. PHENOMENA IN THE SOLAR PHOTOSPHERE 41
If the leaders and followers are regarded as magnetic bipoles, the orientation
of these bipoles is opposite on opposite hemispheres.
The magnetic axes of the bipoles are inclined slightly towards the equator,
the leader spot being closer. This inclination is about
Towards the end of a cycle spot groups appear at high latitudes with reversed
polarity, they belong to the new cycle whereas those with normal polarity for
the old cycle occur close to the equator. This is illustrated in the so called
butterfly diagram (see Fig. 4.6).
absorption coefficient
emission coefficient
For simplicity we consider a homogeneous, plane-parallel atmosphere stratified by
gravity. Then, the properties depend only on the height The surface of the
atmosphere in a strict mathematical sense is where no interactions take place, i.e.
the particle densities are extremely low. The optical depth is defined by:
44 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
which is called Kirchhoff’s law, being the Planck function. We can progress
to solve the transport equation:
In this equation being the angle between the normal to the disk center
and the point where observations are done.
From a Taylor series expansion of about a not specified one gets
where was specified to That means, one observes under the angle to
approximately the source function at optical depth
Let us consider two energy levels in an atom which have the quantum numbers
(lower level) and (upper level). The number of atoms per in the lower level
is und in the upper level Of course a transition from to corresponds
to an absorption process, where a photon of energy is absorbed.
Thus the number of transition per is given by:
is the transition probability for the transition On the other hand let
us consider the number of spontaneous transitions from which is independent
on the intensity J:
and
4.3. THE CHROMOSPHERE 45
where we have put the Planck function. Let us also substitute the Boltzmann
formula:
where are the statistical weights of the states This was first found by
Einstein. We must also consider the induced emission which are transitions from
depending on the intensity J. The number of induced emissions is written
as:
In an induced emission process, the photons emitted have the same directions and
phases as the inducing photons. Thus a detailed balancing in thermodynamic
equilibrium reads as:
Let us assume that the source function within each transition between two energy
levels is independent of frequency, i.e. the source function in a given line is assumed
to be frequency independent which means and have the same frequency
dependence, writing
46 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
and
this integral has to be extended over all frequencies within the line. Then we find:
Please note, that in deriving that equation we have not made any assumption of
an equilibrium condition. This expression is quite general.
How to determine the occupation numbers The number of collisional
excitation per and s is described by:
4.3. THE CHROMOSPHERE 47
thus
which yields
is the rate coefficient for collisional de-excitation and the Einstein co-
efficient for spontaneous emission. If the collisions dominate radiative
transitions and one obtains S = B(T), which is local thermodynamic equilibrium
LTE. Let us assume << 1. For a resonance line The
number of collisional transitions from level u to level 1 is
48 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
therefore, The source function is dominated by the radiation field and only
loosely coupled to the electron temperature via rare collisions.
By these calculation one can understand the typical profile of the Ca II H and
K lines (see Fig. 4.7). There are two intensity minima on the blue and red side
of the line center (called towards the line center two maxima (called
and then at the line center there is a minimum This indicates
that the temperature increases in the corona. While the source function decouples
from the Planck function it reaches a minimum exhibits a small maximum
and finally drops towards the line center. The profile of the well known line is
simpler, there is just a pure absorption. That can be explained with the structure
of the H atom.
A review about the diagnostics and dynamics of the solar chromosphere can
be found in Kneer and Uexküll (1999).
as persistent as those seen in CII line, and may have the configuration of maximum
intensity peaks corresponding either to forward or reflected shocks.
Ca II grains can also be used as indicators for shocks. Therefore spatio-
temporal correlations between enhanced magnetic fields in the quiet solar inter-
network photosphere and the occurrence of Ca II grains in the overlying
chromosphere were investigated by Lites et al. (1999).
Cauzzi et al. (2000) analyzed the temporal behavior of Network Bright Points
(NBPs) using a set of data acquired during coordinated observations between
ground-based observatories (mainly at the NSO/Sacramento Peak) and the Michel-
son Doppler Interferometer onboard SOHO. The NBP’s were observed in the
line and were found to be cospatial with the locations of enhanced magnetic field.
The “excess”’of intensity in NBPs, i.e. the emission over the average value
of quiet regions, is directly related to the magnetic flux density. Thus in analogy
with the Ca II K line, the line center emission can be used as a proxy for
magnetic structures.
Simultaneous CaII K-line spectroheliograms and magnetic area scans were used
to search for spatial correlation between the CaII bright points in the interior
of the network and corresponding magnetic elements and 60% of the bright points
spatially coincided with magnetic elements of flux density > (Sivaraman
et al. 2000).
supergranular structures.
In 1969 Parker and Jokipii hypothesized that the random fluid motions associ-
ated with solar supergranulation may influence the interplanetary magnetic field.
Magnetic footpoints anchored in the photosphere execute a random walk and the
resulting magnetic variations are carried away by the expanding solar wind. The
solar satellite mission Ulysses has observed the resulting large-scale magnetic-field
fluctuations in the solar wind.
By spatio-temporal averaging of two-dimensional velocity measurements ob-
tained in the MgI 5173 line November et al. (1981) found the“mesogranulation”,
in order to indicate the supposed convective character of the phenomenon with a
typical scale of 5 - 10 Mm and a lifetime of approximately 2 h.
The convective nature of the mesogranulation as well as the supergranulation
is not sure. E.g. Rieutord et al. (2000) assign mesogranular flows with both
highly energetic granules, which give birth to strong positive divergences (SPDs)
among which we find exploders, and averaging effects of data processing. A similar
explanation is suggested for the supergranulation.
Hathaway et al. (2000) analyzed power spectra from MDI observations. The
spectra show distinct peaks representing granules and supergranules but no dis-
tinct features at wavenumbers representative of mesogranules or giant cells. The
observed cellular patterns and spectra are well represented by a model that in-
cludes two distinct modes - granules and supergranules.
is typical for gradual or eruptive flares. The intersection of the newly formed
flare loops with the solar surface can be observed: two parallel ribbons in
Therefore, in the older literature we find the designation double ribbon
flares for eruptive flares.
Eruptive flares are very important because of their complexity and associa-
tion with geomagnetic storms; during the event the cosmic ray intensity is
also lowered.
Confined or impulsive events may also result from loop top magnetic recon-
nection. An impulsive flare of say is typically spread over an area of
several in Therefore, the main difference between eruptive and
impulsive flares may be the order of intensity.
Radio bursts and flares: solar flares are associated with radio bursts which
are observed at wavelengths ranging from mm to km. The radio classification
scheme was developed during the 1950s by Australian and French solar radio
astronomers. The different types can be easily recognized in the so called
dynamic spectrum: in such a diagram on the x-axis the time is plotted and on
the vertical axis the frequency. Since the frequency varies with height, one
can easily study the evolution with height of this phenomenon that means
the propagation throughout the solar corona.
Bursts of type III and type V are characteristic phenomena of impulsive
flares (or the impulsive or initial phase of fully developed eruptive flares).
Type III bursts and their associated type V continua are attributed to flare-
accelerated electrons moving along open field lines into the corona. Type II
and type IV bursts are most commonly identified with eruptive flares. Type
II bursts show a slow drift emission which can be interpreted by a shock
wave moving out through the corona with a speed of ~ 500km/s. Type IV
emission is related to magnetic reconnection in CME.
Type II radio bursts result from plasma radiation associated with a MHD
shock propagating through the corona. More than 90% of type II bursts have
an associated flare. They accompany 30% of flares with an importance
class 2 and 3. 70% of all type II bursts are associated with a CME.
Type III bursts occur generally in active regions and thus also without
flares. They are related to suprathermal electrons propagating upward in
54 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
the corona. Dynamic spectra of type III bursts are characterized by the
above mentioned frequency drifts from high to low frequencies as the beam
excites plasma at lower densities with increasing height in the corona. The
type III dm bursts are related to flares and described above.
4.3.7 Prominences
Prominences are great areas of luminous material extending outwards from the so-
lar atmosphere and were first observed during eclipses. They can also be observed
in the light of . Over the photosphere they appear as dark filaments, at the
limb as bright structures. Some prominences are short lived eruptive events, others
can be quiescent and survive many rotational periods of the Sun. The upper parts
are often in the hot corona. Quiescent prominences are made up of material that
is cooler than the photosphere. They often appear as huge arches of dense cool
material embedded in the hot corona. The length of the arch is typically several
100 000 km and the height up to km. A quiescent prominence may change
into an eruptive prominence. The typical thickness of the loop is km. At the
end of its life, a prominence disperses and breaks up quietly or it becomes eruptive
or matter falls back down the field lines to the photosphere. The particle densities
range from which is a hundred times greater than coronal values.
A possible mechanism to understand cool prominence material (temperature
about K) is thermal instability. The equilibrium of the corona requires:
Suppose now that this equilibrium is disturbed locally. The density of the corona
increases in such a disturbed region and it will become cooler than its surroundings.
If we assume that thermal conduction from the hotter surroundings cannot restore
equality of temperature, the dense region will continue to cool until it reaches a new
equilibrium in which its heat input balances its heat output. When a magnetic field
is present, particles can only move along the field lines, this means that thermal
conductivity parallel to the field lines is very much greater than As a result,
56 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
the longest dimension of any cool material is likely to be along the field. The
equation of equilibrium of a magnetized fluid acted on by a gravitational field,
in the is:
where is the gas constant and the molecular weight. In a simple model Kip-
penhahn and Schlüter (1957) assumed that the temperature T and the horizontal
magnetic field components were constant and that P, and were func-
tions of alone. The prominence is represented as a plane sheet.
Helmet streamers are large cap-like coronal structures with long pointed
peaks. They are found usually over sunspots and active regions. Often a
prominence or filament lying at the base of these structures can be seen.
Helmet streamers are formed by a network of magnetic loops that connect
the sunspots in active regions and help suspend the prominence material
above the solar surface. The closed magnetic field lines trap the electrically
charged coronal gases to form these relatively dense structures. The pointed
peaks are formed by the action of the solar wind blowing away from the Sun
in the spaces between the streamers.
Polar plumes are long thin streamers that project outward from the Sun’s
north and south poles. At the footpoints of these features there are bright
areas that are associated with small magnetic regions on the solar surface.
These structures are associated with the ”open” magnetic field lines at the
Sun’s poles. The plumes are formed by the action of the solar wind in much
the same way as the peaks on the helmet streamers.
Coronal Holes: From X-ray observations it was seen that the temperature of
the corona is not uniform. The lower temperature regions are called coronal
holes. They are particularly prominent near sunspot minimum and near the
solar poles. Coronal holes tend to form near the centers of large unipolar
magnetic regions; a comparison of the X-ray images with those of magnetic
field lines calculated on the assumption that the observed photospheric field
line structures extend into the corona as potential fields indicates that they
are regions of open (diverging) magnetic fields. Coronal holes can also be
observed in spectroheliograms taken in the 10 830 Å line of Helium. They
tend to rotate more slowly than sunspots or supergranular patterns and not
differentially.
The fast-speed solar wind originates form the coronal holes (e.g., Krieger
et al., 1973), and accordingly they are considered the main reason for the
”recurrent” type of geomagnetic activity. They may form at any latitude.
For the solar cycle of greatest importance are the unipolar coronal fields.
When the polar fields are strongest during sunspot minimum polar coronal
holes are well defined. They disappear during the polar field reversals near
sunspot maximum.
58 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
isotropic into all directions and this will carry off no momentum- the matter will be
pushed away from the Sun and thus the dust tails are produced. But observations
showed that there is also a plasma tail consisting of ionized gas. If the Sun emits
a continuous stream of plasma, the ionized solar gas would collide with atoms -
momentum is transferred and charge exchange reaction occur: an electron will be
exchanged between an incoming charged particles and a neutral cometary particles
which produced the plasma tail. Since the charged particles move around magnetic
field lines, the plasma tail is aligned with the local interplanetary field.
E.N. Parker predicted the existence of a solar wind from theoretical arguments
showing that a hot corona would imply a continuous stream of plasma.
The solar wind varies in strength through the solar activity cycle. It has an
average speed at the Earth of about 400 km/s. The total mass loss is a few
This is about 1 million tons of solar material flung out into space
every second. If the solar wind was the same in the past then today the total mass
loss of the Sun over that period would be in the order of Also, the mass
loss rate is comparable with that due to nuclear reactions.
The solar wind flows along the open magnetic field lines which pass through
coronal holes. Additionally to the solar wind, the Sun also looses mass by coronal
mass ejections (CME’s). Some of them but not all are accompanied by solar flares.
Low speed winds come from the regions above helmet streamers we have discussed
above while high speed winds come from coronal holes. However, if a slow moving
stream is followed by a fast moving stream the faster moving material will interact
with it. This interaction produces shock waves that can accelerate particles to
very high speeds.
As the Sun rotates these various streams rotate as well (co-rotation) and pro-
duce a pattern in the solar wind much like that of a rotating lawn sprinkler. At
the orbit of the Earth, one astronomical unit (AU) or about km from the
Sun, the interplanetary magnetic field makes an angle of about 45 degrees to the
radial direction. Further out the field is nearly transverse (i.e. about 90 degrees)
to the radial direction.
60 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
As the solar wind expands, its density decreases as the inverse of the square of
its distance from the Sun. At some large enough distance from the Sun (in a region
known as the heliopause), the solar wind can no longer ”push back” the fields and
particles of the local interstellar medium and the solar wind slows down from
400 km/s to perhaps 20 km/s. The location of this transition region (called the
heliospheric termination shock) is unknown at the present time, but from direct
spacecraft measurements must be at more than 50 AU. In 1993 observations of 3
kHz radiation from Voyagers 1 and 2 have been interpreted as coming from a radio
burst at the termination shock. This burst is thought to have been triggered by
an event in the solar wind observed by Voyager 2. From the time delay between
this triggering event and the observation of the 3 kHz radiation, the distance of
the termination shock has been put between 130 and 170 AU.
As it has been stated already, the particle density of the solar wind varies.
From May 10-12, 1999, the solar wind dropped to 2% of its normal density and
to half of its normal speed. This severe change in the solar wind also changed the
shape of Earth’s magnetic field and produced an unusual auroral display at the
North Pole.
The chemical composition of the solar wind is interesting to investigate since
it gives us hints about its origin, i.e. the sources. The most important fact is that
the solar wind composition is different from the composition of the solar surface
and shows variations that are associated with solar activity and solar features
(Bochsler, 2001).
Also magnetic clouds have been observed in the solar wind. These are produced
when solar eruptions (flares and coronal mass ejections) carry material off of the
Sun along with embedded magnetic fields. These magnetic clouds can be detected
in the solar wind through observations of the solar wind characteristics - wind
speed, density, and magnetic field strength and direction.
References on magnetic clouds can be found in Burlaga et al. (1981). About
one half of all magnetic clouds have (and usually drive) upstream interplanetary
shocks, or steep pressure pulses, that in most cases possess large energy- and dy-
namic pressure-increases across their ramps in a stationary frame of reference.
When such a sharp upstream pressure increase encounters the Earth’s magne-
tosphere it pushes it in causing a major reconfiguration of its boundary current
system measured on the ground usually some (5-10) hours before the start of the
main phase of a magnetic storm (Lepping, 2001).
Coronal mass ejections (CMEs) and solar eruptions in general are assumed to
result from quasi-static changes in the photospheric magnetic field which increase
the magnetic energy in the corona and cause sudden release of the stored energy.
This hypothesis is also called storage-release hypothesis. Chen (2001) discusses a
new theory to explain the physics of CMEs. This theory claims that the initial
structure is a magnetic flux rope that is ultimately connected to the solar dynamo
in the convection zone and that magnetic energy propagating from the source along
the submerged magnetic structure enters the corona and drives the eruption. It
predicts that CMEs evolve into interplanetary magnetic clouds (MCs).
Let us give some theoretical arguments of the solar wind and describe its prop-
erties in more detail. Suppose the hot corona sits in static equilibrium on the top
4.5. THE SOLAR WIND 61
of the solar atmosphere. In such a case the pressure gradient in the corona must
be balanced by the gravitational attraction of the Sun:
In this equation we have replaced the variable M by since the mass of the
corona is negligible to the total mass of the Sun. We also can write:
is the number of particles per unit volume and is the average particle mass.
Please also note that is the kinetic temperature of the corona which is far
from thermodynamic equilibrium.
In the corona, conduction is important for energy transport and if is the
coefficient of heat conduction, then
where are radius and temperature at some point in the corona. Combining
all four above equations one gets P and as a function of When expanding this
to the Earth one gets a kinetic temperature of K and a particle density of
Parker pointed out that a solution of such a system to the edge of the
solar system gives nonsense. At large values of the value of P becomes constant,
so that This is higher than the pressure of the interstellar medium and
thus a static model of the corona does not make sense.
If the material of the corona moves outward with a velocity in the radial
direction, then equation 4.59 together with 4.60 becomes
The thermal conduction equation must also be modified to allow for the outward
flow of kinetic energy. The resulting set of equations can only be solved numeri-
cally. Let us do the following substitutions:
62 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
Where is the radius at the base of the corona and the value of there.
We obtain then:
... mass density, wave velocity amplitude, Alfvén speed and A is the
cross section of the flow tube. The magnetic flux is constant, so that the
wave velocity amplitude changes with density as
This scattered radiation is called interplanetary Lyman alpha radiation and SWAN
observes interplanetary Lyman alpha radiation from all directions of the sky. These
Hydrogen atoms collide with solar wind protons and get ionized. This yields to
an ionization cavity around the Sun. But the form and shape of this cavity is
dependent on the solar wind. Therefore the measurement of the interplanetary
UV glow permits to determine the solar wind latitudinal distribution. If the
solar wind were isotropic, the hydrogen distribution and the Lyman alpha emission
pattern would be axisymmetric around the direction where the interplanetary
hydrogen flows into the solar system. However, this is not true.
Planetary Magnetospheres
Here we briefly describe measurements of the magnetic fields of other planets which
are useful as diagnostics of the solar wind. The magnetic field of Mercury and the
structure and dynamics of Mercury’s magnetosphere are strongly influenced by the
interaction of the solar wind with Mercury. In order to understand the internal
magnetic field, it will be necessary to correct the observations of the external field
for the distortions produced by the solar wind. The satellites Helios 1 and 2 made
a number of passes in the region traversed by the orbit of Mercury; thus it is
possible to investigate the solar wind environment of Mercury. The variables that
govern the structure and dynamics of the magnetospheres of Mercury and Earth
are approximately 5-10 times larger at Mercury than at Earth. Thus, the solar
wind interaction with Mercury will be much stronger than the interaction with
Earth (Burlaga, 2001). The solar wind is not constant and since Mercury is closer
to the origin of it, the solar wind at Mercury is probably more variable than that
at Earth.
Mercury, Earth, Jupiter, Saturn, Uranus, Neptune, and Ganymede (satellite
of Jupiter), have presently-active internal dynamos while Venus, Mars, at least
two of the Galilean moons, the Earth’s moon, comets and asteroids do not. These
active dynamos produce magnetic fields that have sufficient strength to stand off
the pressure of the exterior plasma environment and on the other hand interest-
ing interactions with the solar wind can be studied. Moreover, e.g. the jovian
magnetosphere includes a strong time-varying energy source that adds to the dy-
namics of its magnetosphere and produces a quite different circulation pattern
than that found at Earth and, presumably, Mercury. Also the unmagnetized plan-
ets Venus, Mars and even comets have induced magnetospheres associated with
the solar wind interaction with their atmospheres. Cometary magnetospheres,
parts of which can be remotely sensed, exhibit spectacular disruptions called tail
disconnections. Even the atmosphereless bodies with weak magnetic fields can
interact with the solar wind. Small magnetic anomalies on the moon and possibly
asteroids cause weak deflections of the solar wind. This is discussed in the paper
of Russell (2001).
Krymskii et al. (2000) investigate the interaction of the interplanetary mag-
netic field and the solar wind with Mars. Data from the Mars Global Surveyor
mission have shown that localized crustal paleomagnetic anomalies are a com-
mon feature of the Southern Hemisphere of Mars. The magnetometer measured
64 CHAPTER 4. PHENOMENOLOGY OF SOLAR ACTIVE REGIONS
This energy is conserved. If the wave moves into a region of lower density,
then the wave amplitude must increase. The wave turns into a shock wave and
there is a strong dissipation of energy. This is converted into heat and the local
temperature raises. However it turned out that a purely acoustic heating of the
corona is not sufficient to explain the high temperatures there. Acoustic heating
may be important in the outer layers of some stars.
Today we assume that the following two processes are the main reason for the
hot corona:
4.7. VARIATIONS OF THE SOLAR DIAMETER 65
MHD waves: as it has been outlined, when a magnetic field is present there
are two characteristic speeds of wave propagation, the sound speed and
the Alfvén speed If magnetic effects are negligible but this is
not the case for the outer solar atmosphere. The heating process by MHD
waves is analogous to the above mentioned acoustic heating. But it has to
be noted again that MHD waves have an anisotropic propagation.
yields small but measurable changes in the telescope focal length. A temperature
gradient of a few degrees from the center of the window to the aluminium cell at
the filter edge can produce a weak lens effect; that corresponds to a focal length of
a few km and changes the telescope focal length by a few parts in The secular
change is also influenced by the degradation of the front window and increased
absorptivity.
Thus the MDI data yield lower values of opposite sign. Since Sofia et al.
(1979) claimed that W ~ 0.075, solar cycle changes which affect the convective
efficiency near the photosphere will have a large effect on the solar radius; the
MDI measurements rule out this high value of W and suggest that solar cycle
luminosity changes are not caused by superficial fluctuations in the outer layers of
the Sun.
5.1 Neutrinos
5.1.1 General Properties
Let us consider the well known Beta decay. If the nucleus of an atom has too many
neutrons the most likely course is that the nucleus emits an electron. This has
the same effect as turning one of the neutrons into a proton. Such electrons are
historically referred to as beta rays having been named before they were identified
as electrons. An example of beta decay is the decay of tritium or into
69
70 CHAPTER 5. TESTING THE SOLAR INTERIOR
is the collective term for electrons and neutrinos and their relatives in the other
families. The electron and the (electron-)neutrino make up the lepton pair
of the first family. In the other two families, the electron-equivalents are called
muon and tau each with their neutrino partner, called mu-neutrino and
tau-neutrino So we have three different charged leptons: electron, muon, and
tau; and three neutrinos, one associated with each of the three charged leptons.
The was discovered in 1988.
When we compare this value to the cross sections in atomic and nuclear physics
which are about we see that neutrinos can penetrate the whole
Sun without being absorbed and therefore they can be used to test our models.
The distance between collisions, the mean free path if the target particles have
a number density is given by:
That means about one neutrino per month per cubic meter of the detector.
Only neutrinos with energies > 0.8 MeV can be detected by this reaction. This
rules out the most numerous low energy neutrinos (first reaction in the pp chain).
The is in a tank containing gallons of perchlorethylene in the
Homestake Gold Mine in Lead, South Dakota. The experiment has to be placed
deep below the surface to avoid contaminating reactions produced by cosmic rays.
Ar is an inert gas, one can extract it from the tank and observe its decay elsewhere.
Neutrino detections are measured by the solar neutrino flux unit defined by:
Since the experiment contains about atoms one has to expect one detection
every
Theoretical models of the Sun predict the following count rates:
It is important to remark that, due to the low threshold for this reaction, most of
the captured solar neutrinos are the low energy ’pp’ neutrinos.
The neutrino interaction rate is very low, of the order of one interaction per
day in the whole detector mass.
produced by neutrinos is radioactive (halflife about 16 days), and decays
by electron capture into (the reverse process of the solar neutrino capture).
5.1. NEUTRINOS 73
where T is the temperature and the mean molecular weight. Sound speeds can
be determined with the aid of helioseismology to a very high accuracy (better
than 0.2% rms throughout nearly the whole sun). Thus one can estimate tiny
errors in the model values of T and as measurable discrepancies in the precisely
determined helioseismological sound speed:
weight profiles that differ significantly from the standard values. This observational
agreement rules out in particular solar models in which deep mixing has occurred.
The best agreement is obtained when including the effect of particle diffusion-
selective sinking of heavier species in the sun’s gravitational field. Models with-
out taking into account of this effect have rms discrepancies between predicted
and measured sound speeds as large as 1% (e.g. Turck-Chièze and Lopez (1993)
whereas models including this effect have rms discrepancies of 0.1% (Bahcall et
al., 1997).
The sound-speed profile in the Sun was determined by carrying out an asymp-
totic inversion of the helioseismic data from the Low-Degree (l) Oscillation Ex-
periment (LOWL), the Global Oscillation Network Group (GONG), VIRGO on
SOHO, the High-l Helioseismometer (HLH), and observations made at the South
Pole (Takata and Shibahashi, 1998). Then the density, pressure, temperature,
and elemental composition profiles in the solar radiative interior were deduced by
solving the basic equations governing the stellar structure, with the imposition
of the determined sound-speed profile and with a constraint on the depth of the
convection zone obtained from helioseismic analysis and the ratio of the metal
abundance to the hydrogen abundance at the photosphere. Using the resulting
seismic model, neutrino fluxes were estimated and the neutrino capture rates for
the chlorine, gallium, and water Cerenkov experiments. The estimated capture
rates are still significantly larger than the observation.
Solar models with helioseismic constraints and the solar neutrino problem are
discussed in Watanabe and Shibahashi (2001) and Roxburgh (1998).
Is there a correlation between neutrino fluxes and solar activity? On the basis
of an analysis of the production rate at the Homestake station for the period
1970-1990, Basu (1992) found that the solar neutrino flux varies with time in
proportion to the solar wind flux. However, Walther (1999) found that there
exists no significant correlation between the Homestake neutrino data up to run
133 and the monthly sunspot number, according to a test that is based on certain
optimality properties for this type of problem. It is argued that priorly reported
highly significant results for segments of the data are due to a statistical fallacy.
where are spherical polar coordinates. If are low, there is a relatively small
number of patches on the solar surface (which oscillate with different directions of
radial velocity). If are large, there is a very large number of such patches. We
speak of a high degree model if is large and conversely if is small. Most of the
observable p-modes have periods between 2 and 10 minutes with 5 minutes as a
characteristic value. These p-modes are trapped near to the solar surface and in
the solar interior. For high values of the modes are trapped close to the surface.
In general the oscillation frequency of any mode depends on the internal properties
of the Sun in the region which the mode can propagate. Thus the observation of
the oscillations can provide information about the manner in which quantities vary
in the solar interior.
If we want to observe solar oscillations we must measure Doppler shifts of the
wavelength of solar spectral lines. These shifts are produced by the motion of
matter towards or away from the observer.
The diagram (Fig. 5.3) is fundamental for helioseismology. This diagram
shows how much acoustic energy there is at each frequency for every one of the
spatial modes of oscillation. A musical instrument should be tuned to a single
frequency and a few harmonious overtones, the Sun resonates in tens of millions of
ways all at the same time. The frequency of each mode reveals a slightly different
part of the Sun’s interior. The spatial modes are identified from patterns on the
dopplergrams that are made each minute. The frequencies are very low compared
5.2. HELIOSEIMOLOGY-SOLAR OSCILLATIONS 79
The first equation is the equation of motion, the second the equation of continuity,
the third the adiabatic equation and the last is the Poisson equation, denotes the
gravitational potential and is the fluid velocity, is an effective ratio of specific
heats which reduces to when is constant. The time derivative
follows the motion of the fluid. It is related to the derivative at a fixed point by
If we have an equilibrium situation:
80 CHAPTER 5. TESTING THE SOLAR INTERIOR
5.2. HELIOSEIMOLOGY-SOLAR OSCILLATIONS 81
Now we consider small disturbances about this equilibrium in which the per-
turbed quantities are functions of all the spatial coordinates and the time. In
the equilibrium there is no dependence on spherical polar coordinates. For any
variable we can write:
If the star is spherical the oscillation frequency does not depend on For the
Sun, the departure from sphericity is small and the real oscillation modes have a
behavior close to that shown above but with different modes having different
frequencies. The oscillation frequency depends on and The three numbers
are related to the numbers of times vanishes in the radial-, and
and
The functions satisfy a system of differential equations and the boundary
conditions have to be defined. Since stars do not have sharp surfaces we may
assume to a first approximation that all waves are totally reflected at the surface
which is defined as the level where density and pressure vanish. A further sim-
plification arises when the change in the gravitational potential produced by the
oscillations is unimportant; for most perturbations this is a good approximation
because some parts of the star are moving inwards and others moving outwards.
We define a perturbation vector by
Brunt-Väissälä frequency N
where
is always real but and N can be imaginary. It can be shown that con-
vection occurs when is negative. We can write our differential equation for
as:
have been gradually spinning down. This also had led to the suggestion that the
Sun might still posses a rapidly rotating core, perhaps highly magnetized which
also could explain the neutrino problem.
It is extremely important to know the internal rotation of the Sun because the
interplay between turbulent motions and rotation with magnetic fields is essen-
tial for the solar dynamo which leads to the observed 22 year cycles of magnetic
activity.
Here, helioseismology can help to understand the internal rotation rate. In a
spherically symmetric star the frequencies depend upon and but not on For
each pair, there is a fold degeneracy. Rotation breaks the spherical
symmetry and lifts the degeneracy. Advection causes a wave propagation with the
Sun’s rotation to have a higher measured frequency than a similar wave propa-
gating against rotation. Thus the difference in frequency of a pair of oppositely
propagating modes is proportional to times a weighted average of the rotation
rate where the modes have appreciable amplitude. Here, denotes
rotation at radius and latitude The resulting frequency splitting is half
the value of this difference.
Results on the study of the internal solar rotation rate from the SOHO/MIDI
instrument are given in Fig. 5.4.
Measurement of oscillations
How can we measure these oscillations? Let us briefly describe the main principles
of a Dopplerimager. Consider the intensity profile of an absorption line. If the
material from which this absorption line is emitted moves away from the observer,
the line will be redshifted according to the Doppler effect. We can use this effect
84 CHAPTER 5. TESTING THE SOLAR INTERIOR
to make velocity images of the solar surface. The light from the Sun is sent
through a filter that alternates between letting through light from a narrow range
of wavelengths on either side of the center of the line. The two light intensities
are measured at every point on the solar surface using an imaging camera. The
difference between the two intensities changes when the spectral line shifts, and
therefore that difference is a measure of the velocity.
Example of Dopplerimages are given in Fig 5.5. In the second example, the
average over 45 min was subtracted thus the Doppler effect due to solar rotation
is eliminated.
MDI is one of the scientific experiments on SOHO. The medium-l data are
spatial averages of the full disk Doppler velocity out to 90% of the solar disk’s
radius; the measurements are taken every minute and one obtains 23 000 bins of
approximately 10 arcsec resolution where solar p modes up to l=300 can be de-
termined. Low-l observables are velocity and continuum intensity images summed
into 180 bins where oscillations up to l=20 can be detected. Analyzing such data
first asymmetries are detected which arise from the interference between an out-
ward directed wave from the source and a corresponding inward wave that passes
through the region of wave propagation. The degree of asymmetry depends on
the relative locations of the acoustic sources and the upper reflection layer of the
modes. Observations of line profiles of solar modes are therefore suitable to test
theories of excitation of solar and stellar oscillations and their interaction with
turbulent convection. Concerning the rotation rate, the main result of that inves-
5.2. HELIOSEIMOLOGY-SOLAR OSCILLATIONS 85
Inversion techniques
As we have explained above, the observed oscillation frequencies depend on the
physical structure of the solar interior, e.g the variation of quantities such as
with If we assume a spherical symmetric sun and ignore rotational splitting, then
we can deduce from our model of the solar interior the corresponding oscillations.
Alternatively one can regard as unknowns and use the observed frequencies
in order to obtain them. This is called the inversion method. The total number of
quantities that can be determined in such a way is equal to the number of observed
oscillations. If more frequencies can be identified, a better model of the internal
structure can be obtained.
The Seismic Structure of the Sun from GONG data is described in Gough et
al. (1996).
Chapter 6
In this chapter we will explain the basic MHD equations which are needed to
understand solar active phenomena such as spots, prominences, flares etc. The
solar dynamo is needed to maintain the solar activity cycle.
The following equation relates the electric current density to the fields producing
it (generalized Ohm’s law):
is the electrical conductivity and is the bulk velocity of the matter. The final
equations depend on the state of matter; if it consists of electrons and one type of
87
88 CHAPTER 6. MHD AND THE SOLAR DYNAMO
ion:
are the number density and velocity of the ions and electrons respec-
tively and are the charges on the ion and the electron.
In astrophysics two simplifications are applied:
The third Maxwell equation (6.3) states that there are no magnetic monopoles.
This is a common experience: a division of a permanent magnet into two does
not separate north and south poles. Electric fields can be produced by separating
positive and negative charges through the fourth Maxwell equation (6.4) however
the attraction between these charges is so strong that charge separation is usually
cancelled out very quickly. Through the second Maxwell equation electric fields
can be produced by time varying magnetic fields. Such fields are only significant, if
there are rapid changes by time varying magnetic fields. Magnetic fields produced
by the displacement current are usually insignificant in astrophysical prob-
lems because electric fields are unimportant; however they can be produced by a
conduction current if the electrical conductivity is high enough. Such magnetic
fields may be slowly variable in time and space.
We therefore neglect and combine the equations:
and obtain:
6.1. SOLAR MAGNETOHYDRODYNAMICS 89
The solution of these equations shows that magnetic fields decay together with
the current producing them. We can derive an approximate decay time: let us
assume the currents vary significantly in distance L, then from (6.10) the decay
time becomes
The wavelength of the spatial variation of the field is the original field
decays by a factor in the time
Let us consider typical fields of stars: the dimension of the star L and the
electrical conductivity are both high (if the gas is fully ionized). Therefore, the
lifetime of a magnetic field could exceed the main sequence lifetime, such a field
is called a fossil field.
The same is not true for the Earth. Its field is produced by currents in a liquid
conducting core and continuously regenerated by a dynamo mechanism.
The electrical conductivity of an ionized gas is That means that the
characteristic time for decay of currents in the outer layers of the Sun is much
less than the solar lifetime, whereas the decay near the center exceeds the lifetime
(since the temperature near the surface is about 6 000 K and near the center about
If the field in the solar interior were a fossil field extending throughout
the Sun, the field in the outer layers would now be current free - similar to the
field of a dipole. However we don’t observe this. The surface field is very complex
and therefore it must be also regenerated by a dynamo. It is conceivable that a
fossil field of the Sun was destroyed at the very early evolution of the Sun, when
it was fully convective before reaching the main sequence. Also helioseismology
argues against a strong field.
The gas pressure can be written as where is the gas constant and the
mean molecular weight. With we must have:
A tube of magnetic flux is lighter than its surroundings and will start to rise which
is called magnetic buoyancy.
We have already argued that the displacement current can be neglected in the first
Maxwell equation and therefore from which and
6.1. SOLAR MAGNETOHYDRODYNAMICS 91
Here is the magnetic diffusivity. The case where the plasma is stationary
was already discussed above the field decays in the ohmic decay time
Let us discuss the case when Then, the field is completely
determined by the plasma motions and the induction equation is the equivalent
to the vorticity equation for an inviscid fluid. The magnetic flux through a
material surface S which is a surface that moves with the field, is:
If the total flux across any arbitrary surface moving with the fluid remains
constant, the magnetic field lines are said to be frozen in to the flow.
If are typical velocity and length-scale values for our system, then the
ratio of the two terms on the right hand side of the induction equation gives the
Magnetic Reynolds Number
density inside the current sheet, M the plasma average molecular weight.
A similar process occurs in coronal loops that were observed in hard and soft x-
rays by Yohkoh and SOHO instruments. Such a coronal loop (see right drawing in
Fig. 6.2) is stretched out by pressure which is provided by buoyancy. A magnetic
structure is buoyant because the particle density is lower there since it contains
larger magnetic energy density. Thus the external pressure is balanced by a lower
gas pressure in conjunction with a magnetic pressure. The top of the loop distends
and reconnection occurs. Particles in the reconnection region accelerate towards
the surface of the sun and out away. Those particles that are accelerated towards
the sun are confined within the loop’s magnetic field lines and follow these lines
to the footpoint of the loop where they collide with other particles and lose their
energy through x-ray emissions. Such processes are the cause of solar flares and
will be discussed in the next chapter.
Magnetic reconnection also provides a mechanism for energy to be transported
into the solar corona.
A similar process occurs in the earth’s magnetotail. The solar wind distends
the Earth’s dipole field so that the field extends far behind the Earth. Earthward
flowing plasma streams with flow velocities up to 1000 km/s (which is close to the
local Alfvén speed) have been observed (Birn et al. 1981).
A recent review on solar MHD was given by Walsh (2001)
6.1. SOLAR MAGNETOHYDRODYNAMICS 93
(See any textbook on fluid dynamics for a derivation of this formula). Now let
us consider the equation of motion in a plasma with velocity the momentum
equation includes the Lorentz force term and other forces such as gravity
and viscous forces:
Here is the plasma pressure. Let us assume a Newtonian fluid with isotropic
viscosity, then may be written as:
is the local gravity acting in the radial direction and the kinematic viscosity.
Let us make thinks more complicated: Consider a frame of reference with angular
velocity at a displacement from the rotation axis:
The three terms in [ ] denote: Coriolis force, change of rotation and centrifugal
force. Stars rotate more rapidly when they are young. Under most circumstances
the latter two terms are small compared with the Coriolis term
determines the constitution of stars. Here is the mean particle mass and the
number of particles per unit volume. If denotes the entropy per unit mass of the
plasma, then the flux of energy (heat) through a star becomes:
L is the energy loss function. This function describes the net effect of all the sinks
and sources of energy. For MHD applications this becomes:
This equation must then be solved with and the ideal gas
law as well as a simplified form of the energy equation.
Let us introduce the concept of scale height. Let
Therefore, the pressure along a given field line decreases with height, the rate of
decrease depends on the temperature structure (given by the energy equation).
If the height of a structure is much less than the pressure scale height, gravity
may be neglected. The ratio is given by gas pressure to magnetic pressure
If any pressure gradient is dominated by the Lorentz force and
(6.40) reduces to:
6.1. SOLAR MAGNETOHYDRODYNAMICS 95
In this case the magnetic field is said to be force free. In order to satisfy (6.45)
either the current must be parallel to (Beltrami fields) or In
the latter case the field is a current free or potential field.
If is not negligible and the field is strictly vertical of the form
then (6.40) becomes:
limb. Let us consider some simple model to reconstruct these features. Let us
assume:
We see immediately that Let us assume that the footpoints of the field
are anchored down into the photosphere (z=0). Projecting the resulting field onto
the plane gives:
Let us write:
This is a motion in a circle around the direction of and the motion has frequency
the magnitude of the velocity is the radius of the orbit, the gyration
radius is
the sense of the accelerated the circular motions depends on the sign of the
electric charge.
Please note that again is charge dependent. Let us assume that is the
gravitational force then
Thus the drift velocity depends on the mass/charge ratio, the ion drift is much
larger than the electron drift; the particles drift in opposite directions, a current
is produced.
Let us consider a large assembly of particles; these particles interact which
is called collision. If is the characteristic time between collisions the collision
frequency is If is large, the particle motions will be disordered
and decoupled from the magnetic field, the fluid will not be tied to the field. If
collisions are relatively rare, not only individual particles but the whole fluid will
be tied to the field. The collisions provide the electrical resistivity of matter; in a
fully ionized gas a good approximation to the value of the electrical conductivity
is:
and
The forces on the gas are the gas pressure P, the gravitational potential and
the magnetic force For a full description of the system we write down two
additional equations:
a) equation of continuity (conservation of mass):
Note that d/dt is the rate of change with time following a fluid element moving
with velocity
6.1. SOLAR MAGNETOHYDRODYNAMICS 99
Consider the simplest case: a medium with uniform density pressure con-
taining a uniform magnetic field We ignore the influence of the gravitational
field and assume that is so large that Now let us assume a
perturbation for any variable in the form of:
is the wave vector, the wave frequency. The dispersion relation between and
when there is no magnetic field is:
Therefore, in that case only one type of waves can propagate – sound waves. The
wave propagates through the fluid at the wave speed which is
called the phase velocity of the wave.
If there is a magnetic field, the force couples to the equation and also the
Maxwell equations must be taken into account. It is very important to note that
the magnetic field introduces a preferred direction into the system. In a uniform
medium, sound waves travel equally strongly in all directions from its source, this
is not true for MHD waves. If we write the magnetic field again in the form
then we find three types of MHD waves:
Alfvén waves : the dispersion relation is given by
fast and slow magnetosonic waves; their dispersion relation is given by:
Let us consider two special cases: if the waves propagate along the field
there are two waves with and the sound wave unaffected by
the field. For wave propagation perpendicular to the field only one wave survives
with When waves propagate anisotropically, it is necessary
to introduce another wave velocity in addition to the phase velocity, the group
velocity, given by with which the wave carries energy or information. The
group velocity for Alfvén waves is always What does that mean? Regardless
of the direction in which it propagates, energy always travels along the field lines
with speed
If the magnetic Reynolds number is very low, the field is unaffected by the
motions, if it is high, it is wound up many times before dissipation occurs. For
an intermediate value of the magnetic field is carried from the center of the
eddy becoming concentrated in flux ropes at the edge. This buoyant flux ropes
rise towards the surface and this leads to the appearance of sunspots. However
we must also take into account that convection involves different length scales.
Large eddies affect the overall structure of the magnetic field as it has been just
described. Others may be influenced e.g. granulation. Granulation is suppressed
in a sunspot. As it was shown earlier, in the absence of magnetic field convection
occurs in a gas, if the ration of the temperature gradient to the pressure gradient
satisfies the relation:
If a vertical magnetic field of strength B threads the fluid, then this has to be
modified to:
Thus a strong magnetic field can prevent convection and a weaker field can interfere
with convection. Note also that the magnetic field cannot prevent motions which
are oscillatory up and down the field lines hut these are likely to be less efficient
at carrying energy.
Hale’s law: as we have mentioned the leader and the follower spot have
opposite polarities. This reverses after 11 years for each hemisphere so that
the magnetic cycle is in fact 22 years.
Sunspot groups have a tilt towards the equator (this is sometimes also called
Joy’s law).
Reversal of the polar magnetic fields near the time of the cycle maximum.
As we know from fundamental physics, magnetic fields are produced by electric
currents. How are these currents generated in the Sun? The solar plasma is ionized
and it is not at rest. There are flows on the solar surface as well as in the solar
interior producing magnetic fields which contribute to the solar dynamo.
The Effect
Let us consider magnetic fields inside the Sun. There the conditions require that
the field lines are driven by the motion of the plasma. Therefore, magnetic fields
within the Sun are stretched out and wound around the Sun by differential rotation
(the Sun rotates faster at the equator than near the poles). Let us consider a north-
south orientated magnetic field line. Such a field line will be wraped once around
the Sun in about 8 months because of the Sun’s differential rotation (Fig. 6.3).
The Effect
However, the field lines are not only wraped around the Sun but also twisted by
the Sun’s rotation. This effect is caused by the coriolis force. Because the field
lines become twisted loops, this effect was called effect. Early models of the
dynamo assumed that the twisting is produced by the effects of the Sun’s rotation
on very large convective flows that transport heat to the Sun’s surface. The main
problem of that assumption was, that the expected twisting is too much and would
102 CHAPTER 6. MHD AND THE SOLAR DYNAMO
produce magnetic cycles of only a couple of years. More recent dynamo models
assume that the twisting is due to the effect of the Sun’s rotation on rising flux
tubes. These flux tubes are produced deep within the Sun.
much higher there- maybe between 1 and 2 m/s. This slow plasma flow carries
material from the polar region to the equator in about 20 years.
Thus the energy that drives the solar dynamo comes from a) rotational ki-
netic energy, b) another part in the form of small-scale, turbulent fluid motions,
pervading the outer 30% in radius of the solar interior (the convection zone).
Mathematical description
Let us discuss some basic mathematics. In the magnetohydrodynamic limit the
dynamo process is described by the induction equation:
The flow is a turbulent flow. In the mean-field electrodynamics one makes the
following assumptions: magnetic and flow fields are expressed in terms of a large-
scale mean component and a small scale fluctuating (turbulent) component. If we
average over a suitably chosen scale we obtain an equation that governs the evolu-
tion of the mean field. This is identical to the original induction equation but there
appears a mean electromotive force term associated with the (averaged) correla-
tion between the fluctuation velocity and magnetic field components. The basic
principles of mean field electrodynamics were given by Krause and Rädler(1980).
The velocity and the field are expressed as:
where
is a mean electric field that arises from the interaction of the turbulent motion
and the magnetic field. This field must be determined by solving the equation for
and here several assumptions are made. First of all we stressed that
104 CHAPTER 6. MHD AND THE SOLAR DYNAMO
This may be a good assumption when considering a fully turbulent velocity field.
However in the Sun we are dealing with a sufficiently ordered convective field
where the Coriolis force plays an important role. The other approximation is a
first order smoothing: That is valid only if Then our
equation reduces to:
where depend on the local structure of the velocity field and on If the
turbulent field is isotropic, then and
If is small compared to the decay time the diffusive term may be neglected
and from 6.86 we get
And finally:
Compared to the normal induction equation, this contains the term and the
eddy-diffusivity coefficient In the mean field dynamo, the magnetic diffusivity
is replaced by a total diffusivity and the equation becomes:
Please note that most often the prime is dropped on however, in the presence of
it is implied to use the turbulent diffusivity. It is assumed that is axisymmetric.
Then it can be represented by its poloidal and toroidal components and
and Neglecting the advection terms:
6.2. THE SOLAR DYNAMO 105
Note that the dynamo action is possible because we have a regeneration of both
toroidal and poloidal fields. Let us consider the source term in the first of the two
above equations. describes a non uniform rotation. It can be argued that this
term is larger than the next term involving This set of equations then describes
the so called The equations describe:
effect: the poloidal field is sheared by non uniform rotation to generate
the toroidal field.
effect: this is the essential feedback. The helicity of the non
axisymmetric cyclonic convection generates an azimuthal electromotive force
which is proportional to the helicity and to
Let us define a characteristic length scale a decay time and
where is of the order of the local radius of rotation and the local angular
velocity. We may rewrite the above equations in terms of the non dimensional
variables and By an elimination of B and neglection of the
terms we arrive at
If and are scale factors giving the orders of magnitudes of and then
This is called quenching. The equation tells us, that once the mean magnetic field
reaches (which is the value for equipartition), the alpha effect is suppressed.
Recently it was also proposed to take into account the extremely high Reynolds
numbers and to describe the alpha quenching by:
This takes into account that the small-scale component of the dynamo generated
magnetic field reaches equipartition with small scale turbulent fluid motions long
before the mean field does. Reviews on the solar dynamo and the emergence of
magnetic flux at the surface can be found in Fisher et al. (2000) and Moreno-
Insertis (1994).
So far we have discussed large dynamos which are invoked to explain the origin
of the solar cycle and of the large scale component of the solar magnetic field. We
should add here that the origin of small scale magnetic fields can also be understood
in terms of dynamo processes. Recent advances in the theory of dynamo operating
in fluids with high electrical conductivity – fast dynamos, indicate that most suf-
ficiently complicated chaotic flows should act as dynamos (Cattaneo, 1999). The
existence of a large scale dynamo is related to the breaking of symmetries in the
underlying field of turbulence (Cattaneo, 1997).
6.3. SOLAR ACTIVITY PREDICTION 107
explain luminosity variations of e.g. the RSCVn stars or BY Draconis stars (having
luminosities The observed lightcurves required circular spots. The
RSCVn stars occur in binary stars were tidal interactions play an important role,
therefore their starspots are quite different from the sunspots. BY Dra stars are
rapidly rotating young low massive stars characterized by intense chromospheric
emission. Large spots on the Sun cause a variation of the integrated flux < 1%,
whereas up to 30 % for RSCVn and BY Dra stars.
Also flares were detected on stars. Here it is extremely important to have
observations in the EUV/X ray window. Generally pre main sequence stars show
high levels of magnetic activity and strong flares. FU Orionis stars may be in a
phase between T Tauri and post T Tauri stars. More details about that topic can
be found in the review of Haisch et al. (1991). So far we have considered only
stars which have an activity level by orders of magnitude larger than the Sun.
Other indicators for stellar activity are:
EUV lines,
He
H and K lines of Ca II,
Mg II.
The size and extent of chromospheric active regions varies dramatically over
the course of the activity cycle. Thus by measuring the H and K lines of other stars
we can infer on stellar activity cycles. One of the programs that is being carried
out since a long time is the HK project. Almost 100 stars have been observed
continuously since 1966; at present the project is monitoring long-term changes in
chromospheric activity for approximately 400 dwarf and giant stars. In order to
compare the data with the Sun observations of reflected sunlight from the Moon
are done at Mt. Wilson and at Sac Peak and Kitt Peak National Observatory.
The sampling of the stars occurs rapidly: usually less than 10 min per star. The
accuracy of the instrument is between 1% and 2%. When plotting the HK index
against the B – V color index (which is a measure for temperature as explained in
chapter 1) then a clear trend can be seen. The HK index increases as the stellar
temperature decreases. At this point one must be careful with the interpretation.
It is not meant an absolute increase but a relative increase because in cooler stars
also the continuum decreases.
In 1972 Skumanich stated the law for the time of stellar rotation and
stellar chromospheric decay; the rotational velocity and the strength of the CaII
emission of a late type star vary inversely with the square root of the star’s age.
However later it was found that except massive T Tauri stars the majority of low
mass stars rotates slowly.
It was also found that there exists a granulation boundary in the HRD at F5
III. Stars of later spectral type begin to develop a convective envelope that grows
for the rest of their evolution. At the boundary these envelopes are extremely thin
(only 3% of the star’s radius). Stars on the right hand side in the HRD of the
granulation border have smaller rotation rates.
110 CHAPTER 6. MHD AND THE SOLAR DYNAMO
In this chapter we discuss the effects of solar variations on the Earth’s atmosphere.
The variations are small at and above 30 km height but increase to a factor of
2 above 200 km. In general, all correlations of the solar cycle with weather and
climate must be treated with great caution, because besides the Sun there are
many other influences.
111
112 CHAPTER 7. THE SUN AND CLIMATE
the Sun (covering all wavelengths from far UV to IR), and b) from IR radiation
reflected at the surface of the Earth.
The overall heat budget of the atmosphere is as follows: the surface receives
17% of its heat directly from the Sun, 15% from solar radiation scattered by clouds
and 68 % from absorption of infrared radiation emitted by the atmosphere. What
happens to the energy that is absorbed by the surface? The greater part (79%) is
returned to the atmosphere in the form of radiation. The remainder part (21%) is
transmitted to the atmosphere by conduction and as by product by the exchange
of water The surface cools when water evaporates and heat is transmitted
to the air as vapor which recondenses to form clouds. Such phase transitions of
play a major role in the energy budget of the lower atmosphere.
The region below 100 km is called homosphere, the region above 100 km the
heterosphere.
7.1.2 Composition
The composition of the Earth’s atmosphere is given in Table 7.1.
Of course there are gases that can vary considerably both in space and time like
nitric oxide, carbon monoxide and ozone. We can also consider the atmosphere as
an extension of the biosphere, especially for gases like Oxygen
is produced by photosynthesis:
114 CHAPTER 7. THE SUN AND CLIMATE
In the absence of this reaction, carbon would accumulate in organic form and the
fuel for photosynthesis (atmospheric would be depleted. If the supply of
is limited such as in the sediments of organic rich swamps and in the stomachs of
ruminants, we get as a product methane
In Table 7.2 the change of the greenhouse gas and other gas concentrations of
the Earth’s atmosphere is given.
The measurements are from:
1
in situ air samples collected at Mauna Loa Observatory, Hawaii (Keeling,
C.D., Whorf, T.P.). 2 : Etheridge, D. M.; Pearman, G. I.; Fraser, P. J. , Tellus,
Series B - Chemical and Physical Meteorology, 44B, no. 4, 282. These authors used
an ice core from the antarctic called DE08.The extracted ice-core air is analyzed
for methane using gas chromatography with flame-ionization detection. The mean
air-age was 35 yr younger than the host ice. 3 Values from Macehead, Ireland. 4
Cape Grim, Tasmania 4 data from Law Dome BHD ice core, Etheridge et al.
In Table 7.2 the present tropospheric concentration estimates are calculated
as annual arithmetic averages; ppm = parts per million ppb = parts per
billion ppt = parts per trillion
The Global Warming Potential (GWP) is generally used to contrast different
greenhouse gases relative to The GWP provides a simple measure of the rel-
ative radiative effects of the emissions of various greenhouse gases and is calculated
using the formula:
where is the instantaneous radiative forcing due to a unit increase in the con-
centration of trace gas, is concentration of the trace gas, remaining at
time, after its release and is the number of years over which the calculation
7.1. THE EARTH’S ATMOSPHERE 115
7.1.3 Paleoclimatology
First of all let us give a definition of climate: Climate is the weather we expect
over the period of a month, a season, a decade, or a century. More technically,
climate is defined as the weather conditions resulting from the mean state of the
atmosphere-ocean-land system, often described in terms of ”climate normals” or
average weather conditions. Climate Change is a departure from the expected
average weather or climate normals.
The reconstruction of past climate changes is one of the great tasks in climate
research. Since there exits only a 140 years instrumental record, we have to use
proxies to reconstruct climate in the past. Some widely used proxy climate data
types are:
Historical data: Historical documents contain a wealth of information about
past climates (diaries, records...)
Corals: Corals build their hard skeletons from calcium carbonate, a mineral
extracted from sea water. The carbonate contains oxygen and the isotopes
116 CHAPTER 7. THE SUN AND CLIMATE
of oxygen, as well as trace metals, that can be used to determine the tem-
perature of the water in which the coral grew. These temperature recordings
can then be used to reconstruct climate during that period of time that the
coral lived. Increased sea surface temperature has negative effects on the
health of coral. The most visible symptom of declining coral health is coral
bleaching.
Fossil pollen: Each species and genus of plants produces pollen grains which
have a distinct shape. These shapes can be used to identify the type of plant
from which they came. Pollen grains are well preserved in the sediment
layers that form in the bottom of a pond, lake or ocean; an analysis of the
pollen grains in each layer tells us what kinds of plants were growing at the
time the sediment was deposited. Inferences can then be made about the
climate based on the types of plants found in each layer.
Tree rings: Since tree growth is influenced by climatic conditions, patterns
in tree-ring widths, density, and isotopic composition reflect variations in
climate. In temperate regions where there is a distinct growing season, trees
generally produce one ring a year, and thus record the climatic conditions of
each year. Trees can grow to be hundreds to thousands of years old and can
contain annually-resolved records of climate for centuries to millennia.
Ice cores: Located high in mountains and deep in polar ice caps, ice has
accumulated from snowfall over many centuries. Scientists drill through the
deep ice to collect ice cores. These cores contain dust, air bubbles, or isotopes
of oxygen, that can be used to interpret the past climate of that area. Let
us briefly discuss one example of isotope measurements: Of the temperature
dependent markers the most important is the ratio of This can
be explained by the fact that water molecules composed of evaporate
less rapidly and condense more readily than water molecules composed of
Thus, in the ice cores one obtains annual layers starting with
rich, becoming poor, and ending up rich.
This process also depends on the relative temperatures of different years,
which allows comparison with paleoclimatic data.
Volcanic eruption: After the eruption of volcanoes, the volcanic ash and
chemicals are washed out of the atmosphere by precipitation and these erup-
tions leave a distinct marker within the snow which washed the atmosphere.
We can then use recorded volcanic eruptions to calibrate the age of the ice-
core (here the deuterium to hydrogen ratio is an important proxy).
Ice cores from Vostok, Antarctica, were the first to cover a full glacial-
interglacial cycle.
Ocean and lake sediments: Between 6 and 11 billion tons of sediment (tiny
fossils and chemicals) accumulate in the ocean and lake basins each year.
How can we infer e.g. from ice cores past climate? The accumulation which is
governed by saturation water pressure was lower during colder periods and vice
7.1. THE EARTH’S ATMOSPHERE 117
on cycles of 100 000 and 400 000 years. It is the combined effect of the 41 000
year tilt cycle and the 22 000 year perihelion cycle plus the small effect from the
eccentricity that influences the climate. These variations of the Earth’s orbit were
first investigated by Milankovich.
To study the effect of these astronomical variations on climate one must take
into account, that orbital changes occur over thousands of years and the climate
system also takes thousand of years to respond. The primary driver of ice ages
seems to be the total summer radiation received in northern latitude zones near
65° north (this is where the major ice sheets formed in the past) and past ice ages
correlate with the 65N summer insolation. Astronomical calculations show that
the 65N summer insulation should increase gradually over the next 25 000 years.
No decline of the 65N summer insolation that is sufficient to cause an ice age is
expected within the next 100 000 years.
The most important sources of information about such changes and the asso-
ciated composition of the atmosphere are the two large ice caps of Greenland and
Antarctica. Analysis of ice cores is the most powerful means we have to determine
how climate has changed over the last few climatic cycles. The concentrations of
the principal greenhouse gases played an important role in the
transitions from cold ice age climates to warmer interstadials. Warm interstadials
have always been accompanied by an increase of the atmospheric concentration
of the three principal greenhouse gases. This increase has been, at least for
vital for the ending of glacial epochs. A highly simplified course of events for the
past four transitions would then be as follows:
changing orbital parameters initiated the end of the glacial epoch
7.1. THE EARTH’S ATMOSPHERE 121
in the second half of the transition, warming was further amplified by de-
creasing albedo, caused by melting of the large ice sheets in the Northern
Hemisphere going parallel with a change of the ocean circulation.
The isotopic records of Greenland ice cores show evidence for fast and drastic
climatic changes during the last glacial epoch. Possible causes and mechanisms
of such changes and their significance as global climatic events are discussed by
Stauffer (2000). Ice core results also enable the reaction of the environment to past
global changes to be investigated. The deglaciation of the northern hemisphere is
described in Alley and Clark (1999). A carbon cycle model was used to reconstruct
the global mean surface temperature during the last 150 Million years showing that
during this period the tectonic forcing such as decrease in volcanic activity and
the formation and uplift of the Himalayas and the Tibetan Plateau dominated the
control of the climate (Tajika, 2001).
7.1.5 Greenhouseffect
Trace constituents of the atmosphere such as absorb energy at longer
wavelengths and thus trap heat radiated by the surface. The effect is very similar
to that of a glass pane in a greenhouse. The atmosphere is transparent to solar
radiation but it is opaque to longer wavelengths. The infrared absorbing gases
return heat to the ground and account for about 70% of the net input of energy to
the surface. If our atmosphere would contain no water vapor and carbondioxide,
the surface temperature would be about 40 K colder than today. This would imply
that large portions of the planet would be covered with ice.
122 CHAPTER 7. THE SUN AND CLIMATE
Since the 1980s there is a growing concern that the increase in the abundance of
carbondioxide caused by combustion of fossil fuels could lead to a general warming
of the global climate (see Fig. 7.7). Similar greenhouse effects arise from the gases
methane, nitrous oxide and chlorofluorcarbons (CFCs). All these gases are referred
to as greenhouse gases due to their ability to trap heat.
7.1.6 Ozone
The absorption of shortwave solar radiation in altitudes above the troposphere is
responsible for the temperature increase in these layers. Ozone absorbs most of
the UV portion of sunlight The absorption process results in
the dissociation of The recombination process involves the reactions:
M
7.1. THE EARTH’S ATMOSPHERE 123
We can make the simplification that throughout the stratosphere all of the radia-
tive energy from the sun that is absorbed by is converted locally to heat. The
heating rate depends on the distribution of Ozone with height and on the incoming
solar energy.
Tropospheric ozone is either produced by oxidation of hydrocarbons and CO or
by downward transportation of stratospheric ozone. Some examples of reactions
are given below.
Natural events such as Volcanic Eruptions can strongly influence the amount
of Ozone in the atmosphere. However, man-made chemicals such as CFCs or chlo-
rofluorocarbons are now known to have a very dramatic influence on Ozone levels
too. CFCs were once widely used in aerosol propellants, refrigerants, foams, and
industrial processes. Changes in the ozone layer caused by release of CFC’s in the
atmosphere have the potential of producing biological damage through increased
UVB radiation. While cloud cover provides protection on the ground against solar
radiation in the visible and near UV wavelengths, biologically damaging radiation
near 300 nm is controlled primarily by the total ozone content.
The ozone is measured in Dobson units. 1 Dobson Unit (DU) is defined to
be 0.01 mm thickness at STP (standard temperature and pressure). Ozone layer
thickness is expressed in terms of Dobson units, which measure what its physical
thickness would be if compressed in the Earth’s atmosphere. In those terms, it’s
very thin indeed. A normal range is 300 to 500 Dobson units.
In the Earths lower atmosphere, near ground level, ozone is formed when pollu-
tants emitted by cars, power plants, industrial boilers, refineries, chemical plants,
and other sources react chemically in the presence of sunlight. Ozone at ground
level is a harmful pollutant. Ozone pollution is a concern during the summer
months, when the weather conditions needed to form it, lots of sun, hot tempera-
tures, normally occur.
The electrons that are emitted by these reactions loose energy by collision,
elastic and inelastic. This can cause further ionization and contribute to the
production of excited states and the associated emission of airglow.
Electrons can be removed by dissociative recombination:
Hydrogen loss
Any particle in the atmosphere is bound to the Earth by the force of gravity. If
we move such a particle a vertical distance then the work is done, m
denotes the mass of the particle, the gravitational acceleration The
work that must be done to escape the gravitational field is where is the
radius of the Earth ~ 6 400 km. All atoms or molecules have a range of speeds that
is described by the Maxwell–Boltzmann distribution. The average kinetic energy
is given by:
than Of course we must also consider collisions (except at the highest level
in the atmosphere). At the high temperatures in the thermosphere (700...2000K),
significant numbers of hydrogen atoms have velocities above the escape velocity
Therefore, hydrogen is lost at a rate of atoms per per
second. These escaping hydrogen atoms are derived mainly from the oceans and
over the past years of the Earth’s history, the sea level has declined by two
meters globally. Of course during this reaction also is set into the atmosphere
which was crucial for the evolution of life.
There is also a significant of loss of helium.
Precambrian Time
Between 4 and 2.5 billion years ago the continents, atmosphere and oceans formed,
as well as one celled organisms known as prokaryotes. These are the ancestors of
present day bacteria and cyanobacteria. The atmospheric oxygen increases and
later multicellular organisms appear.
Paleozoic Era
This era is subdivided into:
Cambrian Period: 570-500 million years ago; multicellular life.
Ordovician Period: 500-435 million years ago; primitive life on land; verte-
brates in the ocean.
Silurian Period: 435-410 million years ago; first plants and insects appear on
land.
Devonian Period: 410-360 million years ago; spiders, mites and amphibians.
Carboniferous Period: 360-290 million years ago; first true reptiles appear;
coal begins to form.
126 CHAPTER 7. THE SUN AND CLIMATE
Mesozoic Era
This era is divided into:
Jurassic Period: 205-138 million years ago; huge dinosaurs, flying pterosaurs,
oldest known birds.
Cretaceous Period: 138 -65 million years ago: global warming, spread of
dinosaurs. Sudden mass extinction (probably due to asteroid impact), 70 %
of all organisms die out at the end of this period.
Cenozoic Era
This is divided into a Tertiary and a Quaternary Period. The Tertiary Period is
further divided into:
Eocene Period: 55-38 million years ago; ancestral forms of the horse, rhinoceros,
camel and other modern groups such as bats, primates etc. Mammals adapt
to marine life.
Oligocene Epoch: 38–24 million years ago; elephants, cats and dogs, mon-
keys.
Miocene Epoch: 24-5 million years ago; global climate cools; establishment
of the Antarctic ice sheet; large apes in Africa and southern Europe.
Pliocene Epoch: 5 -1.6 million years ago; climate becomes cooler and drier.
Mammals are well established as the dominant terrestrial life form; ancestors
of modern humans.
Pleistocene Epoch: 1.6 million to 10 000 years ago; most recent global ice
age; glacier ice spreads out over more than one-fourth of Earth’s land surface;
modern humans arise and begin their migrations.
Holocene Epoch 10 000 years ago to present; global climate moderates; ice
sheets retreat from Europe and North America; sea levels rise.
7.2. EARTH’S HISTORY AND ORIGIN OF THE ATMOSPHERE 127
Catling et al (2001) after photosynthesis separated the oxygen from the hydrogen,
the authors argue, the two components followed separate paths. The free oxygen
remained in the Earth’s crust, while the hydrogen went on to combine with car-
bon in a process known as ”methanogenesis,” producing methane. When methane
travelled to the upper atmosphere, ultraviolet radiation from the Sun dissolved
it into its components. The light hydrogen drifted away into space and was lost
forever to the Earth’s atmosphere.
Because the hydrogen was lost while the oxygen stayed on Earth, an excess
of oxygen gradually accumulated. When the Earth’s crust was saturated, the
oxygen spilled out and flooded the ancient atmosphere, creating the oxygen rich
environment we know today. This can also solve the faint young Sun problem (see
next chapter).
Of course these facts are extremely important to answer questions like:
How did life begin and first evolve on Earth?
How have conditions on Earth over the past four billion years changed and
affected life?
What are the most extreme conditions under which life can be found on
Earth?
Where else in our own and other planetary systems have conditions possibly
been suitable for life?
How should one search for evidence of fossil or living microbes at extrater-
restrial sites such as Mars or Jupiter’s moon Europa?
The field of Astrobiology tries to find answers to these questions.
In this formula is the present solar luminosity and the present age of the
Sun (4.6 Gyr). Other explanations of a possible different solar luminosity at the
early evolution of the Sun are:
Revisions in the standard solar model in order to solve the neutrino problem.
7.3. THE FAINT YOUNG SUN 129
Strong mass loss during the early phase (Willson et al. 1987).
For our investigations it is clear that a change in solar luminosity over time
would have affected the Earth’s radiation balance and, thus, its climate. If the
Earth is considered to radiate like a blackbody, S is the solar constant (at present
the Stefan-Boltzmann constant , A the planetary Albedo (~ 0.3),
the effective radiating temperature can be obtained by:
The relevant albedo to use here is the Bond Albedo, which is the percentage of
the total incident solar radiation/reflected back into space. The present effective
radiating temperature of the Earth is ~ 255 K. If we combine 7.8 and 7.9 then
the increase of was about 20 deg over geologic time if the albedo of the Earth
is assumed to remain constant. We must also take into account the Earth’s mean
surface temperature and
Thus:
When old sea floor is subducted and carbonate sediments are subjected to higher
temperatures and pressures is returned to the atmosphere/ocean. Then re-
action 7.14 goes in the opposite direction, calcium silicate is reformed and gaseous
is released. Much of this escapes through volcanoes. That process is
termed carbonate metamorphism and on the young Earth the rate of carbonate
metamorphism could have been augmented by faster rates of tectonic cycling and
by impact processing of carbonate rich sediments.
It is important to note that the rates of the weathering reactions are strongly
dependent on temperature. The reaction rates increase with temperature and
weathering requires liquid water. The temperature dependence of the silicate
weathering process rate leads to a negative feedback between atmospheric
and surface temperature: if the surface temperature were to decrease (because of
a faint young Sun), the weathering rate would also decrease and carbon dioxide
would begin to accumulate in the atmosphere. This increase of causes an
increase in the greenhouse effect and thus the temperature increases. The reverse
would happen if the climate became warmer: the weathering rate would increase,
would decrease and the greenhouse effect would become smaller (Walker et
al. 1981). This mechanism can explain why the temperature on Earth was high
enough for liquid water even when the solar luminosity was smaller.
The modern rate of release from volcanoes would create a 1-bar
atmosphere in only 20 Myr if carbonates were not forming. This shows that the
response time of the system is quite fast in geologic terms.
132 CHAPTER 7. THE SUN AND CLIMATE
Here is the partial pressure in the soil and today we have and
obtain On a vegetation free Earth and would be reduced by a
factor of The carbon cycle is only balanced when Therefore,
without vegetation, the atmospheric and surface temperature would have
to increase to bring back the silicate weathering rate to its present value. We
substitute equation 7.16 into equation 7.18 and solve for and obtain:
This shows that under the assumption that land plants pump up soil by
a factor of 40, the effect of eliminating them would be to increase the Earth’s
7.4. THE ATMOSPHERE’S RESPONSE TO SOLAR IRRADIATION 133
temperature by only 6 deg. The net cooling effect of the biota should be somewhat
larger because of the influence of the organic carbon cycle; today 20% of the carbon
is organic carbon rather than carbonate. One can estimate that if life suddenly
were eliminated in total the temperature would increase by 8 deg. Thus even a
lifeless Earth would apparently be no warmer than the real Earth was during the
Cretaceous, when the dinosaurs flourished.
The studies of Schwartzmann and Volk (1989) showed that biota may acceler-
ate chemical weathering by stabilizing soil (silicate minerals stay in contact with
carbonated water), generating organic acids. This could lead to enhanced weath-
ering rates of up to 1000 instead of 3. Therefore, the partial pressure on a
lifeless Earth might be as high as a few tenths of a bar and the surface temperature
may be up to 60 K warmer!
Radiative component
So far we have only the discussed the long term solar variability- summarized as
the faint young Sun problem and the influence of the changing parameters of the
Earth’s orbit on climate (Berger, 1980).
We now address to the question whether there exists also a variability of the
solar input on shorter timescales. The total solar irradiance describes the radiant
energy emitted by the sun over all wavelengths that falls vertically each second on
1 square meter outside the earth’s atmosphere. This is the definition of the solar
constant. Because of the influences of the Earth’s atmosphere this constant is
extremely difficult to measure on the surface and the most reliable measurements
can only be done from space. In Table 7.6 the satellite measurements and the
respective time spans of the measurements are summarized.
The VIRGO Experiment on the ESA/NASA SOHO Mission has two types
of radiometers to measure total solar irradiance (TSI): DIARAD and PMO6V. A
description of the instrument can be found in Fröhlich et al. (1995). Let us shortly
describe the DIARAD measurement facility which is a part of SOHO/VIRGO:
DIARAD is a Differential Absolute Radiometer. It is composed of two cylindrical
7.4. THE ATMOSPHERE’S RESPONSE TO SOLAR IRRADIATION 135
136 CHAPTER 7. THE SUN AND CLIMATE
cavities coated inside with diffuse black and mounted next to each other on the
same heat sink. The flat bottom of the cavities are in fact heat flux transducers on
which heating elements have been mounted. Both cavities see the same thermal
environment through accurately know circular apertures. A comparison of the
power generated inside the cavities (very similarly as an househould weight can
be used) is done. For instance a constant electrical power is generated in one of
the channels and the difference between the two heatflux sensors is automatically
brought back to zero by an ad hoc accurate servosystem that provides electrical
power to the other channel called ”active channel”. This one is regularly irradiated
by the Sun or closed. The difference of the electrical power fed to the active channel
when its shutter is open (exposed to the Sun) and when it is closed is proportional
to the incident solar irradiance. From time to time, the roles of the left channel
and the right channel are reversed for half an hour with the purpose of monitoring
the aging of the continuously exposed left channel. The sampling rate of the
PMO6 instrument is 1 solar total irradiance / 2 minutes, for DIARAD 1 solar
total irradiance / 3 minutes.
The ACRIM contains four cylindrical bays. Three of the bays house inde-
pendent heat detectors, called pyrheliometers, which are independently shuttered,
self calibrating, automatically controlled, and which are uniformly sensitive from
the extreme UV to the far infrared. Each pyrheliometer consists of two cavities,
and temperature differences between the two are used to determine the total so-
lar flux. One cavity is maintained at a constant reference temperature, while the
other is heated 0.5 K higher than the reference cavity and is exposed to the Sun
periodically. When the shutter covering the second cavity is open, sunlight enters,
creating an even greater difference in cavity temperatures. The power supplied to
the second cavity by the ACRIM electronics decreases automatically to maintain
the 0.5 K temperature difference between the two cavities. This decrease in the
amount of electricity is proportional to the solar irradiance entering the cavity.
Exposing the sensors to the space environment and the Solar UV radiation causes
some small changes on the surface of the cavities which may affect the measure-
ments. The ACRIM instrument monitors this type of problem by carrying three
similar sensors, two of which are normally covered. At times these are opened for
comparison purposes. Further details can be found in Wilson (1981, 1984).
Measuring the solar constant one finds:
Part of the energy is blocked by dark sunspots and subsequently released in
faculae. The screening effect by sunspots is overcompensated by the energy
storage and release. This is demonstrated in Fig. 7.10.
There are variations of the solar constant with the solar cycle.
First measurements with the ACRIM 1 (Active Cavity Radiometer Irradiance
Monitor) experiment on board the Solar Maximum Mission and the ERB exper-
iment on the Nimbus-7 satellite showed a positive correlation between the solar
cycle activity, measured by the sunspot index, and the total solar irradiance. The
peak to peak variation of about (out of about 1367) between solar maxi-
mum and minimum was reported by Fröhlich (1987), Willson and Hudson (1988)
7.4. THE ATMOSPHERE’S RESPONSE TO SOLAR IRRADIATION 137
and Foukal and Lean (1988). Somewhat larger fluctuation up to 0.2% occur over
timescales of days and weeks.
Given that the total variation between the peaks of solar cycles 21 and 22 was
about 0.1%, how much is the effect to be expected for a change of the correspond-
ing global temperature on Earth? It is expected that this change of the solar
irradiation produces a corresponding variation of about 0.2° C in globally aver-
aged equilibrium surface temperature (Hansen and Lacis, 1990). But there is some
considerable delay in the response. Because of the thermal inertia of the oceans,
the time needed to approach equilibrium is much longer than 11 years (e.g., Reid,
1991), so that the actual temperature response to the observed variation during a
solar cycle is likely to be considerably smaller, and probably insignificant from a
climatic point of view.
In order to study a long term variation of the solar output, there is no direct
observational support. It is therefore necessary to use proxy data or solar activity
indicators. Sunspot index measurements exist over a time span of roughly 350
years and they suggest the presence of a 76-80 yr cycle, the Gleissberg cycle,
modulating the 11 yr cycle (Sonett, 1982, Berry, 1987). Foukal and Lean (1990)
gave an empirical model of total solar irradiance variation between 1874 and 1988.
The presence of the solar cycle has been claimed in various sets of proxies:
auroral activity,
is ”the present is the key to the past,” originally stated by James Hutton in
1785. However, dendrochronology adds a new ”twist” to this principle: ”the
past is the key to the future.” In other words, by knowing environmental
conditions that operated in the past (by analyzing such conditions in tree
rings), we can better predict and/or manage such environmental conditions
in the future. Hence, by knowing what the climate-tree growth relationship
is in the 20th century, we can reconstruct climate from tree rings well before
weather records were ever kept! Let us give one example from Grissino-
Mayer: Fig 7.11 shows a long-term precipitation reconstruction for northern
New Mexico based on tree rings. How this reconstruction was made? The
reconstruction was developed by calibrating the widths of tree rings from
the 1900s with rainfall records from the 1900s. Because we assume that
conditions must have been similar in the past, we can then use the widths
of tree rings as a proxy (or substitute) for actual rainfall amounts prior to
the historical record.
Individual tree-growth series can be ”decomposed” into an aggregate of en-
vironmental factors:
ingly similar to records from tree rings in North America, Europe and western
Russia. The general trends reflected in the tree-ring record include cooler
conditions in the early 1700s, followed by warming that started mid-century.
An abrupt cooling occurred in the late 1700s and continued for much of
the 1800s. The coldest period was between 1830s and 1870s, after which
a steadily increasing warming trend began. An example of this analysis is
given in Fig 7.12.
Solar radius variations (Gilliland, 1981),
sedimentary rocks (Sonett and Williams, 1985)
sea surface temperatures (Gerard, 1990). The mechanism how this could be
related to solar irradiance variations works as follows:
1. absorption of solar energy by the tropical oceans in a deep surface layer,
2. transport of that energy by ocean currents,
3. transfer of that energy by evaporation into atmospheric moisture and
pressure systems leading to more precipitation (Perry, 1994).
Lewis et al. (1990) showed that solar radiation in visible frequencies, usually
assumed to be absorbed at the sea surface, penetrates to a significant depth below
the upper mixed layer of the ocean that interacts directly with the atmosphere. In
clear water, the blue wavelengths, where the greatest amount of energy is available,
penetrate the deepest, to nearly 100 m. Energy injected into the ocean at this
depth can be stored for a substantial period of time.
140 CHAPTER 7. THE SUN AND CLIMATE
As it has been stated above the transparency of the tropical oceans is dependent
upon the amount of biogenic material, phytoplankton pigments, and degradation
products that are present. In the Pacific Ocean, transparency increases from east
to west, with the greatest penetration of solar energy occurring in the western
tropical Pacific. Due to ocean currents, the North Pacific Ocean takes approxi-
mately 4 years to move temperature anomalies from the western tropical Pacific
to near North America (Favorite and McLain, 1973).
During the prolonged period between 1500 and 1850, average temperatures in
Northern Europe were much colder than they are today, this is known as the little
Ice Age. The coldest part of this period coincides with a conspicuous absence of
sunspots and other signs of solar activity, called the Maunder Minimum.
For example Gilliland reported a 76 year cycle in the solar radius, inferred
from a 258 year record of transits of the planet Mercury, solar eclipse records and
meridian transit measurements. Ribes et al. (1987) also reported as Gilliland
that the solar radius is slightly increased in times of low solar activity during the
Maunder minimum.
A review book on the role of the Sun in climate change was written by Hoyt
and Schatten (1997) where other references can be found.
UV Radiation
Solar radiation shortward 320 nm represents only 2% of the total solar irradiance;
0.01% of the incident flux is absorbed in the thermosphere at about 80 km and
0.2% in the stratosphere above 50 km. This radiation is extremely important
since the thermal structure and photochemical processes above the troposphere
are controlled by it. The stratosphere is controlled by absorption and dissociation
of in the 175 to 240 nm range. The 205 to 295 nm range is predominantly
absorbed by ozone If there is a stratosphere- troposphere coupling, this could
affect also the climate. The short term variation of UV radiation is ascribed to
the evolution and rotation of plage regions on the solar disk. The XUV induced
thermospheric temperature changes is shown for low and high solar activity in Fig.
7.13. Solar activity is measured in terms of the 10.7 cm radio flux and of the plage
area
Energetic particles
There are three main contributions:
electrons: they reach the high latitude thermosphere after interaction with
the geomagnetic field and acceleration;
high energy solar protons: their flux is enhanced during periods of large
flares;
galactic cosmic rays: they originate from outside the heliosphere but their
input on Earth is partly controlled by solar activity.
7.4. THE ATMOSPHERE’S RESPONSE TO SOLAR IRRADIATION 141
The diurnal variations show a day/night ratio of 1.28 over the equator with
the peak occurring about 2 p.m. During the night, heat is conducted down
from the top of the thermosphere to its base, the mesopause, where it is
radiated.
The effect of the sun can be expressed and measured by the 10.7 cm radio flux
which is given in the units If this quantity is multiplied
by a factor of 1.8 deg per unit of flux, we obtain the temperature.
The other variation comes from the solar activity cycle (see Table 7.7).
contribution from auroral heating- which is triggered by the solar wind- must be
taken into account.
Finally we have to stress that both UV and particle precipitation have chemical
effects and the most important is the production of N, NO and (which is
collectively called The following reactions define the production in
the thermosphere:
To test these predictions it is important to have data at the time scale of the
solar cycle; however we must also take into account the instrumental drifts as well
as the typical lifetime of the instruments which normally are below 5 yr. In the
stratosphere, the ozone response is caused primarily by changes in production from
and has a maximum value of 0.5 % for a 1 % change in the UV at 205 nm. The
study of the response of the temperature has been made by Hood (1986, 1987b)
and Keating et al. (1987) between 30 and 0.2 mbar (24 to 60 km) and later by
Clancy and Rusch (1989) up to 90 km. They establish the already mentioned 0.5
% response. The very small temperature response lags the UV by 4 to 14 days.
A study by Angell and Korshover (1975, 1976) established a correlation of the
ozone column with solar activity with a peak to peak variation up to 10 % at 70°
latitude and only 4 % at 47°.
There seems to be no correlation of polar stratospheric temperatures and solar
activity (Labitzke (1987), Labitzke, Van Loon (1988), Kerr (1988). There exists
a stratospheric biennial oscillation which is more or less periodic and reversal of
winds in the lower equatorial stratosphere with an average period of 27 months.
7.4.4 Troposphere
As we have seen above, only wavelengths below 300 nm penetrate to the tropo-
sphere and surface. We have already stressed that this part of the solar spectrum
is only slightly variable with a peak to peak variation of about 1 part in 1400. Thus
the troposphere which contains 90 % of the total mass of the Earth’s atmosphere
is subject to a nearly constant driving solar energy.
However, there have been innumerable attempts to find correlations between
solar activity and various meteorological phenomena and other variables. If the
troposphere is to be significantly influenced by the tiny changes of solar irradiation,
there should exist a very strong mechanism of amplification (trigger mechanism).
Such mechanisms were discussed:
magnetospheric effects by electric field - including also effects of thunder-
storms (Mc Cormac and Seliga, 1979)
Hines (1974) suggested a change of the transmissivity of the stratosphere to
upwardly propagating atmospheric waves (Callis et al. 1985 showed from
models that this is possibly not the case)
The effect found by Labitzke (1987): temperatures in the polar winter are
jointly influenced by the solar cycle and the quasi biennial oscillation and
the effect on the troposphere is discussed in Van Loon and Labitzke (1988).
Eddy (1976, 1988) discussed the absence of sunspot activity during the 17
th century which is known as the Maunder minimum and an earlier event,
called the Spörer minimum. Both periods seem to coincide with periods of
reduced global temperatures the more recent is called the Little Ice Age.
Eddy (1988) showed that the required solar input reduction would have to
be much greater than the tiny amplitudes detected on the time scale of a
solar cycle. Maybe also amplifying factors have to be considered.
7.4. THE ATMOSPHERE’S RESPONSE TO SOLAR IRRADIATION 145
be be reconstructed using the aa-index) has roughly doubled since 1900. This is
in good agreement with concentration in Greenland ice (Beer, 2000).
is produced by the interaction of cosmic rays with constituents of the Earth’s
atmosphere. The cosmic ray flux is modulated by the heliospheric magnetic field.
Lean et al. (1995) assumed that the background irradiance is proportional
to the amplitude of the solar cycle; Hoyt and Schatten (1993) propose a trend
corresponding to cycle length and Baliunas and Soon (1995) demonstrated that
the amplitudes of stellar cycles (observed in Ca II H and K) scale with the length
of the stellar cycle.
A short overview of long term chances in solar irradiance was given by Solanki
and Fligge (2000).
The question whether the Earth’s climate is influenced by solar activity has
a central position in the present debate about the global warming. Greenhouse
gas concentrations have a continuous increase and do not follow the observed
decrease in the 1900’s and in 1940-1970 example. These variations might be better
explained when solar activity is taken into account. During a normal sunspot cycle
the irradiance changes by 0.1% but could be greater (e.g. during the Maunder
Minimum 0.3%, Lean (1997)).
Let us now describe the interaction of cosmic ray particles with the heliosphere.
The heliosphere is defined by the interplanetary magnetic field. This shields the
interstellar plasma (charged particles). Interstellar neutral gas flows through the
solar system however, since uncharged particles are not influenced by magnetic
fields. The speed is approximately 25 km/s. When approaching the Sun, these
neutral atoms become ionized by two processes:
photo-ionization: an electron of the neutral atom is knocked off by a solar
high energy photon (e.g. a UV photon);
charge exchange: an electron is exchanged to an ionized atom of solar wind
particle.
As soon as these particles are charged the Sun’s magnetic field carries them out-
ward to the solar wind termination shock region. The ions repeatedly collide with
the termination shock, gaining energy during each collision. This continues until
they escape from the shock region and diffuse back toward the inner heliosphere.
Such particles are called anomalous cosmic rays (ACRs). ACRs are thought to
originate from the very local interstellar medium and are not related to the above
mentioned violent processes as the GCRs. They can easily be discerned from
GCRs because they have lower speed and energy. They include large quantities of
He, O, Ne and other elements which have in common high ionization potentials.
The third component of cosmic ray particles are Solar energetic particles (SEP).
They move away from the Sun due to plasma heating, acceleration and other
processes. Flares e.g. inject large amounts of energetic nuclei into space, the
composition varies from flare to flare. On the scale of cosmic radiation, SEPs have
relatively low energies.
the troposphere. With a neutron monitor one can measure the low energy part in
the GCR spectrum.
For our study here, it is important to note that by measuring cosmic rays one
can derive a proxy for solar activity very long back in time. This is possible since
isotopes in the atmosphere are produced by cosmic rays. From such recordings a
good qualitative agreement between cold and warm climatic periods and low and
high solar activity during the last 10 000 years was found. When we consider
variations during the last millennium, one can deduce, that from 1000-1300 AC
solar activity was very high which coincided with the warm medieval period. We
know from history that e.g. during that period the Vikings settled in Greenland.
The solar activity - if it is well represented by the variation- decreased and a
long period followed which is now called the little ice age (in this period falls also
the so called Maunder Minimum, 1645-1715, where practically no sunspots were
observed). This period lasted until the middle of the 19th century. From then on,
solar activity has increased and is the highest in the last 600 years.
Thus we may assume the following connections:
If that assumption is true, there is a mechanism, how the Earth’s climate can
be influenced by the Sun.
during nighttime. The surface of the Earth begins to lose energy in the form of
longwave radiation. This causes the ground and the air above it to cool. The
precipitation that results from this kind of mechanism takes the form of dew, frost
or fog.
Of course these mechanisms may act as a combination: convection and oro-
graphic uplift can cause summer afternoon showers in the mountains.
Let us compare the levels of cloud cover for summer and winter (northern
hemisphere). For summer in the northern hemisphere, highest levels of cloud
cover occur over the mid-latitude cyclone storm tracks of both hemispheres, In-
tertropical Convergence Zone over land surfaces, and the Indian Monsoon region
(orographic lifting). Lowest values occur over the subtropical deserts, the sub-
sidence regions of the subtropical oceans, and the polar regions. For winter in
the northern hemisphere highest levels of cloud cover occur over the mid-latitude
cyclone storm tracks of both hemispheres and the Intertropical Convergence Zone
over land surfaces. Lowest values occur over the subtropical deserts, the subsidence
regions of the subtropical oceans, and over the South Pole.
Clouds influence vertically integrated radiative properties of the atmosphere.
They cause a cooling through reflection of incoming shortwave radiation (sun -
light) and heating by absorption and trapping of outgoing long wave radiation
(thermal radiation). Let us consider the net radiative impact of a cloud: this
mainly depends on two parameters, on its height above the surface and its optical
thickness. High optically thin clouds tend to heat while low optically thick clouds
tend to cool. The net forcing of the global cloud cover is in the range between
as it is derived from climate models. Thus a significant influence
on the global cloud cover can be potentially very important for Earth’s climate
(see also Table 7.8).
It has been found that the Earth’s cloud cover follows the variation in GCR.
It seems to be that the ionization in the atmosphere produced by GCR is the
essential link. One can estimate that a variation in cloud cover of 3 % during an
average 11-year solar cycle could have an effect of This is a very
significant amount.
The idea that cosmic rays can influence cloud formation was first pointed out
by Svensmark (1997). He showed that there was a significant correlation between
total cloud cover over the Earth and the influx of cosmic rays. The rays ionize
7.5. COSMIC RAYS 151
particles in the low troposphere which then seed the growth of cloud water droplets.
During the past century the shielding from cosmic rays has increased since solar
activity has increased. This decreases the formation and cooling influence of low
clouds and may thus provide a possible contribution to the global warming of the
past 100 years (Marsh, Svensmark, 2000; Svensmark, 1999).
Let us consider and summarize the changes in the magnetic field in the solar
atmosphere. Shorter solar cycles facilitate a rise in the coronal source flux, longer
cycles allow it to decay. The accumulation of the coronal source flux strongly
depends on the rate of flux emergence in active regions. In general the peak and
cycle averaged sunspot numbers are larger when cycles are shorter. Therefore,
shorter cycles are associated with larger flux emergence rates, there is less time
for the open flux to decay. We can state:
shorter activity cycle increased coronal flux
The coronal source surface is where the magnetic field becomes approximately
radial. This occurs at This surface can also be regarded as the
boundary that separates the solar corona from the heliosphere. The magnetic flux
threading the corona source surface is called or open solar flux. If there is a
rise of the flux than the cosmic ray flux incident on the Earth will decrease.
Lockwood and Foster (2000) estimated that the cosmic ray flux > 3 GeV was 15%
larger around 1900 than it is now. As it was shown above, cosmic rays generate air
ions in the sub ionospheric gap which allows current to flow in the global electric
current. This connects thunderclouds with the ground via lightening.
152 CHAPTER 7. THE SUN AND CLIMATE
155
156 CHAPTER 8. SPACE WEATHER AND RADIATION DAMAGE
Gamma rays originate from the nucleus of an atom. They are capable of travelling
long distances through air and most other materials. Gamma rays require more
”shielding” material, such as lead or steel, to reduce their numbers than is required
for alpha and beta particles.
In Table 8.1 we give some definitions used in radiation physics.
The effect of radiation on any material is determined by the dose of radiation
that material receives. Radiation dose is simply the quantity of radiation energy
deposited in a material. There are several terms used in radiation protection to
precisely describe the various aspects associated with the concept of dose and how
radiation energy deposited in tissue affects humans.
Some terms related to radiation dose:
Chronic dose: A chronic dose means a person received a radiation dose over
a long period of time.
Acute dose: An acute dose means a person received a radiation dose over a
short period of time.
Somatic effects are effects from some agent, like radiation that are seen in
the individual who receives the agent.
Genetic effects: Genetic effects are effects from some agent, that are seen in
the offspring of the individual who received the agent. The agent must be
encountered pre-conception.
Teratogenic effects: Teratogenic effects are effects from some agent, that are
158 CHAPTER 8. SPACE WEATHER AND RADIATION DAMAGE
seen in the offspring of the individual who received the agent. The agent
must be encountered during the gestation period.
ions in cosmic rays and solar particle events is risk assessment in carcinogenesis.
single strand break in the DNA: this can be usually repaired and normal cell
function is restored.
breaks in both DNA strands: usually the damage is too severe to repair and
the cell dies.
3.9-2.5 Billion years ago the Earth was dominated by an oceanic lithosphere.
Cockell (2000) calculated that the DNA damage rates might have been approxi-
mately three orders of magnitude higher in the surface layer of the Archean oceans
than on present-day oceans. However, at 30 m depth, damage might have been
similar to the surface of present-day oceans. On the other hand, risk of being
transported to the surface water in the mixed layer was quite high. Thus the
mixed layer may have been inhabitated by a low diversity UV-resistant biota. Re-
pair capabilities similar to Deinococcus radiodurans would have been sufficient
to survive in the mixed layer. During the early Proterozoic ozone concentrations
increased and the UV stress would have been reduced and a greater diversity of
organisms could have inhabitated the mixed layer.
Lean (2000) discusses societal impacts of solar electromagnetic radiation.
The Yohkoh satellite was launched in 1991. Song and Cao (1999) discuss CCD
radiation damage. Evans et al. (1999) discuss charged-particle induced radiation
damage of a HPGe gamma-ray detector during spaceflight.
less than 10% of ultraviolet light, and sunblock creams work by absorbing or
reflecting UV rays. The SPF rating of sunscreens gives an indication of their
effectiveness as UV blockers. For example, an SPF of 15 means that it should take
15 times as long to before skin damage occurs (i.e., the cream should block about
93% of the radiation that causes skin damage).
UV radiation is subdivided into three wavelength bands:
UVC (220-290 nm); totally blocked by ozone and other gases, does not reach
the Earth’s surface.
A person’s potential to develop skin cancer is related to their exposure UVB
radiation (sunburn). In New Zealand, about one person in three will develop a skin
cancer during their lifetime. About half the number killed on the roads die of skin
cancer in New Zealand. New Zealand and Australia have a very high melanoma
incidence compared with other countries.
How can this be explained?
New Zealanders have an outdoor lifestyle,
wear fewer clothes now than in the past,
the ancestors of most white-skinned New Zealanders migrated from the UK,
which is at much higher latitude, and has much lower levels of UV radiation.
These people are therefore poorly adapted to the relatively high levels of UV
naturally present in New Zealand;
calculations suggest that locations in the Southern Hemisphere should receive
approximately 15% more UV than locations at a similar latitude north of the
Equator (Basher, 1981; McKenzie, 1991). This is caused by differences in
ozone between the Northern and Southern Hemispheres, and also because the
Earth is slightly closer to the Sun during the Southern Hemisphere summer
(McKenzie and Elwood, 1990);
measurements show much larger differences, with biologically-damaging UV
being 50-80% more in the Southern Hemisphere than at comparable North-
ern latitudes in Europe. The differences are caused by the buildup of tropo-
spheric pollution (tropospheric ozone and aerosols) in the North (Seckmeyer
and McKenzie, 1992);
the largest levels of UV. Much higher levels of UV are experienced in coun-
tries, such as Australia, which are closer to the equator.
The amount of UVB light at ground level is determined by three factors: a) solar
elevation, b) the amount of ozone in the atmosphere and c) the cloudiness of the
sky. Please note that during local noon the amount of background radiation is
the same as direct radiation three hours before and afterwards. At NZ’s latitude,
approximately 40% of the daily sunburn radiation occurs during the two hour
period centered on solar noon.
Since the late 1970s an ozone hole has formed over Antarctica during early
spring. The amount of ozone over New Zealand varies seasonally with a maxi-
mum in spring and a minimum in early autumn. Evidence of ozone destruction
has also been observed over the Arctic. The ozone hole is caused by the special
meteorological conditions of the cold atmosphere above polar regions which am-
plify the destructive ability of CFCs. The Antarctic ozone hole cannot shift over
New Zealand. However, ozone losses over Antarctica may contribute to changes
in ozone over the whole globe.
164 CHAPTER 8. SPACE WEATHER AND RADIATION DAMAGE
In Fig. 8.3 clear sky UV indices for different stations in NZ are given. Seasonal
variations are higher at low latitudes as well as the absolute values.
8.3.5 UV Index
The Global Solar UV Index was developed through the WHO. It provides an
estimate of the maximum solar UV exposure at the Earth’s surface. The intensity
of UV reaches a maximum around mid-day (when there is no cloud cover) at solar
noon.
It is generally presented as a forecast of the maximum amount of skin-damaging
UV expected to reach the Earth’s surface at solar noon. The values of the Index
range from zero upward; the higher the Index number, the greater the likelihood
of skin and eye damaging exposure to UV, and the less time it takes for damage
to occur.
166 CHAPTER 8. SPACE WEATHER AND RADIATION DAMAGE
Close to the equator, summer-time values reach 20. During a European summer
a value of 8 can be reached. We speak of:
low UV exposure: Index 1...2
moderate UV exposure: Index 3...4
high UV exposure: Index 5...6
very high UV exposure: Index 7...8
extreme UV exposure: Index > 9
8.4 Spacesuits
Outer space is extremely hostile and without a spacesuit the following things would
happen:
you would become unconscious within 15 s because there is no O,
blood and other body fluids start to boil and then freeze because there is no
air pressure,
tissues (skin, heart...) expand because of the boiling fluids,
extreme temperature changes: sunlight 120°C, shade -100° C.
exposure to radiation and micrometeoroids.
169
170 CHAPTER 9. THE IONOSPHERE AND SPACE WEATHER
receivers must correct the delay of the GPS signal as it propagates through the
ionosphere to the GPS satellite (at 22 000km altitude).
The ionosphere may become highly turbulent, mainly in the high latitude and
low latitude F region and at special times (often after sunset). In this context
turbulence is defined as small scaled structures (scale length cm to m) which are
irregular and embedded in the large scale ambient ionosphere (tens of kilometers).
In the equatorial region plasma irregularities are generated just after sunset and
may last for several hours. At high latitudes these irregularities may be generated
during day and night. Both effects occur most frequently during the solar cycle
maximum. Radio signals become disrupted by these perturbations and the effect
is known as ionospheric scintillation. The bigger the amplitude of the scintillated
signal the greater the impact on communication and navigation systems.
9.1.1 Aurora
Energetic particles from the Sun are carried out into space along with the ever
present hot solar wind. This wind sweeps toward Earth at super sonic speeds
ranging from 300 to 1000 km/s. It distorts the Earth’s magnetic field which forms
out a comet shaped magnetosphere.
The magnetosphere protects us from energetic particles and most of them are
deflected around the Earth. The Earth’s Van Allen Belts consist of highly energetic
ionized particles trapped in the Earth’s geomagnetic field. On the sunward side of
172 CHAPTER 9. THE IONOSPHERE AND SPACE WEATHER
the Earth, the geomagnetic field is compressed by the solar wind, on the opposite
side of the Earth, the geomagnetic field extends. Thus the field forms an elongated
cavity which is also known as the Chapman-Ferraro Cavity around the Earth.
Within this cavity are the Van Allen Radiation Belts. The radiation belts are
composed of electrons (keV) and protons (MeV).
As it is seen in Figure 9.2 there are two belts of particle concentration: a)
small inner belt between 1 and 2 Earth radii where protons of energy 50 MeV and
electrons with energies > 30 MeV reside and b) outer larger belt from 3 to 4 Earth
radii where less energetic protons and electrons are concentrated. The inner belt
is relatively stable, the outer belt varies in its number of particles by as much as
a factor of 100.
Charged particles trapped in the belts spiral along the field lines while bouncing
between the northern and southern mirror points. Particles in the inner belt may
interact with the upper atmosphere causing the auroral oval which is an annulus
centered over the magnetic poles and around 3000 km in diameter during quiet
times. The location of the auroral oval is usually found between 60 and 70 degrees
of magnetic latitude (north and south).
There are many shapes and features of aurorae. They generally start at 100
km above the surface and extend upward along the magnetic field for hundreds
of km. Auroral arcs can nearly stand still and then suddenly move (dancing,
turning). After midnight one often sees a patchy appearance of aurorae, and the
patches blink on and off every 10 s or so. Most of aurorae are greenish yellow and
9.1. GENERAL PROPERTIES 173
sometimes the tall rays turn red at their top and along their lower edge. On rare
occasions sunlight hits on the top creating a faint blue color.
The different colors depend on the specific atmospheric gas, its electrical state
and on the energy of the particle that hits the atmospheric gas. Atomic oxygen is
responsible for the two main colors of green (557.7 nm, at a height below 400 km)
and red (630.0 nm, about 400 km or higher). Excited nitrogen also emits light
(600-700 nm; below 200 km). Auroral displays are intensified if the interplanetary
magnetic field is in the opposite direction to the Earth’s magnetic field. The geo-
magnetic storms produce brightness changes and motion in the aurorae and these
are called auroral substorms. Recent models of aurorae explain the phenomenon
by a process of release of energy from the magnetotail, called magnetic reconnec-
tion. Regions of opposite magnetic fields come together and the magnetic field
lines can break and reconnect in new combinations. The point of reconnection in
the magnetotail lies usually at 100 Earth radii. When the solar wind adds suffi-
cient magnetic energy to the magnetosphere, the field lines there overstretch and
a new reconnection takes place at 15 Earth radii, the field collapses and electrons
are injected into the atmosphere.
Reconnection stores large amounts of energy in the Earth’s magnetic field until
it is released explosively. The cycle of energy storage and release is called substrom.
Multiple substorms lead to magnetic storms and acceleration of particles to very
high energies. These particles damage satellites.
The geomagnetic field is measured by magnetometers and the data are often
given as 3-hourly indices that yield a quantitative measure of the level of geomag-
netic activity. The K-index is given from 0 to 9 and depends on the observing
station. The globally averaged index is a measure for the global auroral activ-
ity.
When geomagnetic activity is low, the aurora typically is located at about 67
degrees magnetic latitude, in the hours around midnight. As activity increases,
the region of aurora expands towards the equator. When geomagnetic activity is
very high, the aurora may be seen at mid and low latitude locations (see Table
9.2) around the earth that would otherwise rarely experience the polar lights.
In Table 9.3 auroral boundaries are given as a function of the Kp index.
The magnetic activity produced by enhanced ionospheric currents flowing be-
low and within the auroral oval is measured by the Auroral Electrojet Index AE.
The definition of this index is as follows: at a certain time the total range of devi-
ation from quiet day values of the horizontal magnetic field (h) around the auroral
oval. Defined and developed by Davis and Sugiura, AE has been usefully employed
both qualitatively and quantitatively as a correlative index in studies of substorm
morphology, the behavior of communication satellites, radio propagation, radio
scintillation, and the coupling between the interplanetary magnetic field and the
Earth’s magnetosphere. For these varied topics, AE possesses advantages over
other geomagnetic indices or at least shares their advantageous properties.
174 CHAPTER 9. THE IONOSPHERE AND SPACE WEATHER
DST Index
DST stands for Disturbance Storm Time. The DST is an index of magnetic activity
derived from a network of near-equatorial geomagnetic observatories that measures
the intensity of the globally symmetrical equatorial electrojet (the ”ring current”).
Thus DST monitors the variations of the globally symmetrical ring current, which
encircles the Earth close to the magnetic equator in the Van Allen (or radiation)
belt of the magnetosphere. During large magnetic storms the signature of the ring
current can be seen in ground magnetic field recordings worldwide as so-called
main phase depression. The ring current energization which results in typical
depression of 100 nT is related to magnetic reconnection processes at the neutral
sheet.
Kp and Ap Index
The K-Index was first introduced by J. Bartels in 1938. It is a quasi-logarithmic
local index of the 3-hourly range in magnetic activity relative to an assumed quiet-
day curve for a single geomagnetic observatory site. The values consist of a single-
digit 0...9 for each 3-hour interval of the universal time day (UT).
The planetary 3-hour-range index Kp is the mean standardized K-index from 13
geomagnetic observatories between 44 degrees and 60 degrees northern or southern
geomagnetic latitude. The scale is 0...9 expressed in thirds of a unit, e.g. 5- is
4 2/3, 5 is 5 and 5+ is 5 1/3. This planetary index is designed to measure solar
particle radiation by its magnetic effects. The 3-hourly Ap (equivalent range)
index is derived from the Kp index (see Table 9.4). This table is made in such a
way that at a station at about magnetic latitude 50 degrees, Ap may be regarded
as the range of the most disturbed of the three field components, expressed in the
unit of 2 g. A daily index Ap is obtained by averaging the eight values of Ap for
each day. The Cp index, the daily planetary character figure, is defined on the
basis of Ap according to Table 9.5
Another index devised to express geomagnetic activity on the basis of the Cp
176 CHAPTER 9. THE IONOSPHERE AND SPACE WEATHER
index is the C9 index which has the range between 0 and 9. The conversion table
from the Cp index to the C9 index is given by 9.6
Sunspot Numbers
The sunspot number index is also often called Wolf number in reference to the
Swiss astronomer J. R. Wolf who introduced this index in 1848; details about how
to obtain that number can be found in the chapter about sunspots and the solar
cycle.
9.1. GENERAL PROPERTIES 177
During a solar flare event a sudden increase of X-ray emission causes a large
increase in ionization in the lower regions of the ionosphere on the sunlit side of
the Earth. Very often one observes a sudden ionospheric disturbance (SID). This
affects very low frequencies (OMEGA) as a sudden phase anomaly (SPA) or a
sudden enhancement of the signal (SES). At HF and sometimes also at VHF an
SID may appear as a short wave fade (SWF). Depending on the magnitude of the
solar flare such a disturbance may last from minutes to hours. At VHF the radio
noise created by solar flares interferes with the signal. The occurrence of solar
flare is modulated by the solar activity.
Flares may also emit energetic particles. The PCA (polar cap absorption) is
caused by high energetic particles that ionize the polar ionosphere. A PCA may
last from days to weeks depending on the size of the flare and the interaction of
the high energetic particles emitted by the flare and the Earth’s magnetosphere.
During these events polar HF communication becomes impossible. A coronal mass
ejection may be a consequence of a large solar flare or a disappearing filament and
is an ejection of a large plasma cloud into the interplanetary space. Such a coronal
mass ejection (CME) travels through the solar wind and may also reach the Earth.
This results in a global disturbance of the Earth’s magnetic field and is known as
a geomagnetic storm. High speed solar wind streams originating in coronal holes
on the Sun’s corona hits the Earth’s magnetosphere and also causes ionospheric
disturbances.
Calculate the currents produced by the geoelectric field in the circuit system
constituted by the network and its earthings.
The first step is generally more difficult, partly because the space and geophysical
input parameters are not well known.
The effects of geomagnetic disturbances on electrical systems at the earth’s
surface were studied e.g. by Boteler et al. (1998) or Lehtinen and Pirjola (1985).
A prediction of Geomagnetically Induced Currents in Power Transmission Systems
was given by Pirjola et al. (2000).
Increased absorption
Depressed MUF
Increases LUF
Increases fading and flutter
Surveillance Systems
Faraday rotation
Scintillation
Loss of phase lock
Radio Frequency Interferences (RFI)
Navigation Systems
Position errors
9.2. SATELLITES 181
9.2 Satellites
9.2.1 Solar Panels
Solar panels are devices that convert light into electricity. Some scientists call
them photovoltaics which means, basically, ”light-electricity” since the solar light
is converted into electric energy by them.
A solar panel is a collection of solar cells. Lots of small solar cells spread over
a large area can work together to provide enough power for satellites or space
stations. The more light that hits a cell, the more electricity it produces, so
spacecrafts are usually equipped with solar panels that can always be pointed at
the Sun even as the rest of the body of the spacecraft moves around.
The most efficient solar panels are the DS1 solar panels which convert about
22 % of the available energy into electrical power whereas most solar panels on
people’s houses convert only about 14 %. It is also important to note that solar
panels lose about 1-2 % of their effectiveness per year. This means after a five year
mission, the solar panels will still be making more than 90 % of what they made
at the beginning of the mission. Of course this also depends on their distance to
the Sun.
There are two major dangers to solar panels in space besides regular wear-and-
tear:
Solar flares that can damage the electronics inside the panels.
Micrometeorites, which are tiny, gravel-sized bits of rock and other space
junk floating in space can scratch or crack solar panels.
Some protection can be made by the use of a thick layer of glass. Of course,
if a satellite’s mission path takes it away from the Sun (further out into the solar
system) solar panels will become less and less efficient.
Another kind of protection to the above mentioned damaging effects can be
made by the use of a solar concentrator. This uses Fresnel lenses which collect a
182 CHAPTER 9. THE IONOSPHERE AND SPACE WEATHER
large area of sunlight and direct it towards a specific spot by bending the rays of
light and focussing them- the same principle when people use a magnifying lens
to focus the Sunlight on a piece of paper which starts a small fire.
Fresnel lenses have been invented in 1822 by Jean Fresnel. Theaters use them
for spotlights. They are shaped like a dart board with concentric rings around a
lens that is a magnifying glass. Solar concentrators put one of these lenses on top
of every solar cell. The solar cells can then be spaced farther apart since the light
is focused on each cell. Fewer cells need to be placed and the panels cost less to
construct. Thick glass or plastic cover over the solar panel are used to protect
them from micrometeorites.
DS1’s photovoltaics are made out of gallium arsenide (GaAs). GaAs is made
into a cylinder that is then sliced into cells. These solar cells are then connected
to the rest of the power network. Solar concentrators, made of clear plastic, are
placed above them to focus the Sun’s rays.
As a summary we give some literature, further references can be found therein.
Markvart et al. (1982) studied the photon and electron degradation of boron-
doped FZ silicon solar cells. Radiation-resistant silicon solar cell were investigated
by Markvart et al. (1987). Defect interactions in silicon solar cells were analyzed
by Markvart et al. (1989). A study of radiation-induced defects in silicon solar
cells showing improved radiation resistance was made by Peters et al. (1992).
General information about solar cells can be found in Tada et al. (1982).
A review on radiation damage in solar cells was given by Markvart (1990).
Solar panels: they provide abundant power for nearly all a satellite’s needs
and are safe and clean to launch. However:
solar panels are large and fragile constructions that are vulnerable to
damage from external forces or even mechanical failures;
they are rather expensive to build and put into space;
they always need to be pointed at the Sun (think about what happens
if they are blocked by planets or other objects);
the farther the satellite gets from the Sun, the less effective solar panels
work. As a rule of thumb we can state that solar powered missions
cannot travel further than the orbit of Mars.
towards Mars. The mission was over. Only then began the desperate attempt to
crash the spacecraft safely. But within minutes it had passed beyond the range of
the Crimean station, and no other Russian site ever heard from it again.
The probe contained 200 grams of plutonium 238 and crashed into the pacific
west of Easter island. People reported on an object that was brighter than Sirius
with a luminous trail five degrees in length and fragmented.
particles can cause a large current impulse sufficient to change the state of a
bistable circuit element.
Heavy Ion SEUs occur directly when a heavy ion passes through a semicon-
ductor memory element. The standard models take into account the size, shape,
and charge sensitivity of the memory element and the energy, angle, and impact
parameter of the incident particle.
For satellites around the Earth, the offset and tilt of the geomagnetic axis with
respect to the Earth’s rotation axis produces a corresponding miss-alignment of the
radiation belts. The result is the South Atlantic Anomaly. The spatial distribution
of errors show a maximum in the Atlantic ocean east of the southern part of South
America. There occurs also a significant number of errors at high latitudes due
to cosmic rays (see Fig. 9.4). These data are from UoSAT-2 which measured
from September 1988 to May 1992; UoSAT-2 monitored almost 9000 Single Event
Upsets (SEU), and the majority of these (75%) occurred in the South Atlantic
Anomaly (SAA) region.
Single event upsets pose also problems to space missions: As a result of vol-
canic action on Io, the innermost of the large Galilean moons of Jupiter, particles
(actually heavy ions) of sulphur and oxygen are present in the space surrounding
the planet. These particles form a part of the Jovian magnetosphere. Although
the origin of these particles is the moon Io, the volcanoes provide enough velocity
for them to escape from the gravitational field of the moon and to become elements
of the magnetosphere around Jupiter.
The heavy ions diffuse both inward and outward from the planet. Many of the
particles diffuse outward to 20 to 50 times the radius of Jupiter (R J , measured
188 CHAPTER 9. THE IONOSPHERE AND SPACE WEATHER
from the planet’s center), where they are accelerated by an interaction with the
massive Jovian magnetic field.
The most critical phase of mission operations for to study the Galilean satellites
of Jupiter occurs at the time of the spacecraft’s closest approach to Jupiter (4 RJ).
Heavy ions are capable of penetrating the delicate electronics in the spacecraft and
causing a stored computer bit to change its value from a ”0” to a ”1” or vice-versa,
a Single Event Upset results (SEU). A single bit flip in one of Galileo’s computer
memories could trigger a chain reaction of erroneous commands with disastrous
results.
Modern microelectronic devices can suffer from single event effects caused by
cosmic radiation neutrons in the atmosphere. The phenomenon has been observed
both on ground and at aircraft altitudes. The neutron flux at aircraft altitudes
(<15 km) is large enough to make the neutron single event effects a problem to
aircraft electronics. The most studied device type is static random access memo-
ries (SRAM) since those devices have a very high density of transistors, making
them sensitive to particle radiation. The cosmic ray neutrons are produced by
the charged primary cosmic radiation in the earth’s atmosphere. Thereby the at-
mospheric neutron flux is certainly influenced by solar activity and space weather
(see e.g. Dyer, 2001)
Normand (1996) studied the effect of SEU in avionics. Ziegler and Lanford
(1979) studied the effect of cosmic rays on computer memories. SEU in implantable
cardioverter defibrillators were studied by Bradley and Normand (1998). They
found some correlation with the expected geographical variation of the secondary
cosmic ray flux.
Prediction of times with increased risk of internal charging on spacecraft are
given by Andersson et al. (1999) and Wu et al. (1999).
and there are also unresolved variations in atmospheric density. The accuracy of
the prediction is in the order of 10 %. That means that one day before re-entry
the uncertainty is at least 2 hours. Within that time however, the satellite will
have circled the globe and thus it is difficult to predict the location of re-entry
with a reasonable warning.
The decay of a satellite’s orbit also depends on the cross section of the object
itself. In Fig. 9.5 a rough estimate of the lifetime of a satellite with effective mass
to cross section ratio in a circular orbit below 300 km is given for two
cases: a) for solar minimum conditions, b) for solar maximum conditions. The
geomagnetic field is assumed to be quiet during this period. The lifetime values
may be varied for satellites of differing mass to area ratios.
The uncertainty in the predictions is shown by a NORAD prediction in April
1979 for the expected re-entry of the SKYLAB space station between 11 June and
1 July of that year. The actual re-entry occurred on July 11, outside the stated
interval, a prediction error from mid-interval of around 15%.
actual variation in molecular mass with height and a compensation term for the
variation in temperature over the considered range from 180 to 500 km. The
variation in density due to the space environment is introduced through T, where
. Generally, the solar X-ray output incident upon the Earth
is absorbed at the base of the thermosphere (120 km) and gives rise to a heating
which propagates upward from this level. We use the solar F10.7cm flux which
can vary from 65 to 300 SFU (Solar Flux Units, 1 as a
proxy for X-rays. The other quantity to take into account is the precipitation of
particles- most of them coming from the Sun (CMEs). These are well correlated
with the large variations in the geomagnetic field measured at the ground level
and quantified by the geomagnetic indices (we use the index here). The is
computed every 24 hours and during quiet periods just above zero but may rise
up to 400.
Then we can write the following set of equations:
The output of this simple model is the density. The intermediate values are only
used to derive this density and may not correspond to true atmospheric values at
any height within the considered range. The temperature e.g. may be regarded
as the mean asymptotic value for the exosphere at large altitudes. The mean
molecular weight might be regarded as an integrated mean value from the base of
the thermosphere up to the specified height.
The solar 10.7 cm radio flux is used in averaged form (average over the last 90
days). A small correction may be made to weight the current flux more strongly.
Now let us consider the satellite drag. When a spacecraft travels through an
atmosphere it experiences a drag force opposite to the direction of its motion. This
is given by:
D... drag force, atmospheric density, speed of the satellite, A... cross
sectional area perpendicular to the direction of motion, drag coefficients.
The latter can vary; at altitudes at which satellites orbit We introduce
the effective cross sectional area
For a circular orbit we have the following relation:
G... gravitational constant, mass of the Earth. The reduction in the period
due to atmospheric drag is given by:
9.2. SATELLITES 191
Re-entry is assumed when the satellite has descended to an altitude of 180 km.
The space environmental parameters are given by the solar 10.7 cm radio flux and
the geomagnetic activity index. Furthermore one has to provide an estimate for
the satellite mass to area ratio. In the absence of any further information this
value can be taken as . This is an average value for many satellites.
If the program underestimates the actual decay of the orbit, you must de-
crease the mass to area ration, in the case of an overestimation the ratio must be
increased.
Also the situation becomes more complicated when considering satellites with
very elliptical orbits. Here, a part of the orbit is outside the current atmospheric
model. They are also subject to other perturbations (Sun, Moon). If the eccen-
tricity is not too large, one can introduce an effective height in that model:
q... is the perigee (lowest height) of the orbit and e the eccentricity. For example,
the lifetime of a satellite in an elliptical orbit with e = 0.01 and q = 400km is the
same as the lifetime of a satellite in a circular orbit of height:
Since that formula is only a rough approximation it should only be used for orbits
with e < 0.1. The solar activity should be constant during the orbit decay. Most
satellites reaching the end of their lives will have orbits with very low eccentric-
ities (i.e. nearly circular). The reason for this is that atmospheric drag acts to
circularize orbits. The apogee height is decreased whilst the perigee height is little
affected until the orbit becomes close to circular.
As we have seen the fast and effective communication of space weather effects to
the public is very important. For that reason the US NOAA (National Oceanic
and Atmospheric Administration) has introduced the space weather scales. A
summary of the different influences triggered by the Sun is shown in Fig. 10.1.
The NOAA space weather scales can be grouped into three different parts:
Geomagnetic storms
Radio blackouts
10.1.1 G1
classified as minor; the influence on power systems is weak, some grid fluctuations
can occur. Also the influence on spacecraft is negligible. It seems however that
migratory animals are affected even at this low level; the aurora is commonly
visible at high latitudes.
As an average about 1700 events per cycle ( corresponding to about 900 days
per cycle) are to be expected. The value is about 5.
193
194 CHAPTER 10. THE NOAA SPACE WEATHER SCALES
10.1.2 G2
moderate; at this activity some damage may occur in power systems: high-latitude
power systems may experience voltage alarms, long-duration storms may cause
transformer damage.
Concerning spacecraft operations, corrective actions to orientation may be re-
quired by ground control; possible changes in drag affect orbit predictions. This
imposes problems to fully automated satellites.
Concerning terrestrial telecommunication the HF radio propagation can fade
at higher latitudes, and aurora has been seen as low as New York and Idaho (
down to 55 geomagnetic latitude). The value is about 6 and on the average
one can expect 600 events per cycle (corresponding to about 360 days per cycle).
10.1.3 G3
strong; on power systems voltage corrections may be required; furthermore false
alarms can be triggered on some protection devices.
On satellite components surface charging may occur. Due to the extension of
the terrestrial atmosphere during these events drag may increase on low-Earth-
orbit satellites, and corrections may be needed for orientation problems.
It is also very important to note that intermittent satellite navigation and low-
frequency radio navigation problems may occur, HF radio may be intermittent,
and aurora has been seen as low as down to 50 geomagnetic lat. The value is
about 7 and on the average one can expect 200 events per cycle (corresponding to
130 days per cycle).
10.2. SOLAR RADIATION STORMS 195
10.1.4 G4
severe; widespread voltage control problems may occur in power systems and some
protective systems will mistakenly trip out key assets from the grid.
The problems of surface charging and tracking of satellites increase consider-
ably.
On surface pipelines, induced currents affect preventive measures; the satellite
navigation can degrade for hours and the low frequency navigation can be dis-
rupted. Aurora has been seen down to 45° geomagnetic latitude. The index is
at 8 and on the average one has to count with 100 events per cycle (corresponding
to 60 days per cycle).
10.1.5 G5
extreme; widespread voltage control problems and protective system problems can
occur; transformers may experience damages and some grid systems may experi-
ence complete collapse or blackouts.
The spacecraft operations are affected by extensive surface charging, problems
with orientation, uplink/downlink and tracking satellites.
At this activity pipeline currents can reach hundreds of amps, HF radio prop-
agation may be impossible in many areas for one to two days, satellite naviga-
tion may be degraded for days, low-frequency radio navigation can be blocked for
hours. The aurora has been seen down to 40° geomagnetic latitude (Italy, south-
ern Texas). At this level and on the average one has to expect 4 events
per cycle (corresponding to about 4 days per cycle).
10.2.1 S1
minor; there are no effects on biological systems and satellite operations; no danger
for astronauts (especially for EVAs (extravehicular activities)). There may be some
minor impacts on HF radio in the polar regions. The ion flux is about 10 (see above
units). There are about 50 events per solar cycle.
10.2.2 S2
moderate; there are no biological influences; for satellite operations it is important
to know that infrequent single-event upsets are possible.
Small effects occur on HF propagation through the polar regions and navigation
at polar cap locations is possibly affected.
The ion flux is about 100 and we have about 25 events per cycle.
196 CHAPTER 10. THE NOAA SPACE WEATHER SCALES
10.2.3 S3
strong; at this level radiation hazard avoidance is recommended for astronauts
on EVA; passengers and crew in commercial jets at high latitudes may receive
low-level radiation exposure (equivalent to approximately 1 chest x-ray).
The effects on satellite operations become important: lots of single-event up-
sets, noise in imaging systems, and slight reduction of efficiency in solar panels are
likely.
On Earth, degraded HF radio propagation through the polar regions and nav-
igation position errors are likely. The ion flux is about and we have about 10
such events per cycle.
10.2.4 S4
severe; unavoidable radiation hazard to astronauts on EVA thus it is necessary to
alarm astronauts; moreover, elevated radiation exposure to passengers and crew
in commercial jets at high latitudes (equivalent to approximately 10 chest x-rays)
is possible.
Satellites may experience memory device problems and noise on imaging sys-
tems; star-tracker problems may cause orientation problems, and solar panel effi-
ciency can be degraded.
On the surface blackout of HF radio communications through the polar regions
and increased navigation errors over several days are likely.
The ion flux is about There are about 3 such events per cycle.
10.2.5 S5
extreme; unavoidable high radiation hazard to astronauts on EVA; high radiation
exposure to passengers and crew in commercial jets at high latitudes (equivalent
to approximately 100 chest x-rays) is possible.
Satellites may be put out of operation, memory impacts can cause loss of
control, may cause serious noise in image data, star-trackers may be unable to
locate sources; permanent damage to solar panels is possible.
At the surface complete blackout of HF communications is possible through the
polar regions, and position errors make navigation operations extremely difficult.
The ion flux is at fortunately, these events occur on a rate fewer than 1
per cycle.
10.3.1 R1
minor; we have to take into account a weak or minor degradation of HF radio
communication on the sunlit side, as well as occasional loss of radio contact.
Concerning navigation we have to consider that low-frequency navigation sig-
nals may be degraded for brief intervals. The physical measurement is M1 and
On the average 2000 such perturbances per cycle occur (on 950 days per
cycle).
10.3.2 R2
moderate; limited blackout of HF radio communication on sunlit side occur, loss
of radio contact for tens of minutes.
Navigation: a degradation of low-frequency navigation signals for tens of min-
utes is likely. The physical classification of the relevant solar event goes M5 and
the flux to
On the average one has 350 events per cycle (300 days per cycle).
10.3.3 R3
strong; a wide area blackout of HF radio communication, as well as a loss of radio
contact for about an hour on the sunlit side of Earth is likely.
Since low-frequency navigation signals are being degraded for about an hour
this also has serious consequences for navigation.
The physical classification is X1, the flux and one has 175 events per cycle
(140 days per cycle).
10.3.4 R4
severe; HF radio communication blackout occurs mostly on the sunlit side of Earth
for one to two hours and a HF radio contact loss during this time has to be
expected.
Outages of low-frequency navigation signals cause increased error in position-
ing of navigational systems for one to two hours. Minor disruptions of satellite
navigation are likely on the sunlit side of Earth.
The physical classification is X10, the flux and one has 8 events per cycle
(8 days per cycle).
10.3.5 R5
extreme; a Complete HF (high frequency) radio blackout on the entire sunlit side
of the Earth lasting for a number of hours may occur . This results in no HF radio
contact with mariners and en route aviators in this sector.
Navigation: Low-frequency navigation signals used by maritime and general
aviation systems experience outages on the sunlit side of the Earth for many hours,
causing loss in positioning. Satellite navigation errors in positioning increase for
several hours on the sunlit side of Earth, which may spread into the night side.
198 CHAPTER 10. THE NOAA SPACE WEATHER SCALES
The physical classification is X20, the flux and one has less than 1
events per cycle.
10.4 Summary
The classification scheme given above enables very easily to estimate the effect of
geomagnetic storms and solar radiation storms on satellites and telecommunication
systems. This is also extremely important for manned space mission (ISS, inter-
national space station). On the other hand, the solar activity is declining again
after heaving reached its maximum in 2000. One can estimate that there will be
about 25 EVA/year necessary for the construction of the space station. For that
reason, it is extremely important to alert astronauts for S4 and S5 storms. The
predicted sales figures for GPS systems rise from 5000 Million USD for 1998 and
more than 9000 Million USD for 2000. This means that more and more systems
are equipped with these navigation systems but on the other hand we must take
into account that small degradations may even occur at R1 levels. The frequency
of such events is however more than 2000 per cycle.
Also the number of satellites will increase.
Skylab is an example of a spacecraft re-entering Earth’s atmosphere prema-
turely as a result of higher-than-expected solar activity because of the increased
drag.
Systems such as LORAN and OMEGA are adversely affected when solar ac-
tivity disrupts their radio wavelengths. The OMEGA system consists of eight
transmitters located through out the world. Airplanes and ships use the very low
frequency signals from these transmitters to determine their positions. During
solar events and geomagnetic storms, the system can give navigators information
that is inaccurate by as much as several miles. If navigators are alerted that a
proton event or geomagnetic storm is in progress, they can switch to a backup
system. GPS signals are affected when solar activity causes sudden variations in
the density of the ionosphere.
We have seen before that some military detection or early-warning systems are
also affected by solar activity. The Over-the-Horizon Radar bounces signals off
the ionosphere in order to monitor the launch of aircraft and missiles from long
distances. During geomagnetic storms, this system can be severely hampered by
radio clutter. That can occur at even low activity (R1 perturbances).
4 billion years ago. After that, the solar wind gradually eroded the martian atmo-
sphere until, today, it has less than 1% of the thickness of the Earth’s.
No global magnetic field and a very thin atmosphere – those are the two factors
that render Mars vulnerable to space radiation. Does that mean that because of
this unprotected exposure Mars is lifeless? It is assumed that certain life forms
could be radiation resistant like the terrestrial microbe Deinococcus radiodurans.
D. radiodurans has a feature that is considered all-important in aerospace: redun-
dancy. Its genetic code repeats itself many times so that damage in one area can
be recognized and repaired quickly. It withstands attacks from acid baths, high
and low temperatures, and even radiation doses, e.g. the microbe can withstand
without loss of viability a dosage that is 3 000 times greater than what would kill
a human.
Tiny Martians might also live in rocks or soil, substances that provide natural
protection against radiation. In addition to that, magnetic fields provide some
protection, ancient global, regional and sub surface fields.
NASA sent a radiation monitor, MARIE (Mars Radiation Environment Ex-
periment) to the Red Planet on April 7th with the 2001 Mars Odyssey spacecraft.
MARIE is one of three scientific instruments on board – the other two will search
for signs of water and interesting minerals on Mars. MARIE can detect charged
particles (electrons, protons, ions) with energies between 15 MeV and 500 MeV.
As we have seen in the previous chapter, particulate radiation poses the greater
threat to humans. Particles from solar flares (protons) are of greater concern -
here particles with relatively low energies (around 70 MeV) are produced. Such
protons lose energy in tissue. Cosmic ray nuclei have energies between 300 and
500 MeV and they penetrate the human body in a short time so that there is
200 CHAPTER 10. THE NOAA SPACE WEATHER SCALES
not enough time to transfer their energy to the surrounding tissue. Solar protons
when passing through humans ionize molecules along their track- the ionization
creates free radicals causing modification or even break of the DNA strands and
if the cell survives it can become cancerous.
Therefore, human settlers on Mars must be protected from these energetic
protons. The air density at the surface of Mars is equivalent to that of the Earth’s
atmosphere at 20 km altitude. Astronauts must be protected by shelter walls.
Since Mars has no global magnetic field the surface is eroded by the solar wind
as well as by the planet’s atmosphere. Currently Mars looses approximately 2
kg/s of its atmosphere. In 1998 magnetometers discovered a network of magnetic
loops arrayed across Mars’s southern hemisphere. Locally, the magnetic fields arch
over the surface like umbrellas, hundreds of km high. In such an area you would
measure a field about as strong as the Earth’s (a few tenths of a gauss). Elsewhere
the field is extremely weak.
The martian ionosphere traces the distribution of the surface magnetic field,
and there seems to be a 1-to-1 correspondence: places where magnetic umbrellas
deflect the solar wind are also spots where the ionosphere is retained high above
the surface of the planet.
The Earth’s global magnetic field is caused and maintained from an active
dynamo – that is, circulating currents at the planet’s liquid metallic core. A
similar dynamo once churned inside Mars, but for reasons unknown it stopped
working four billion years ago. The patchwork fields mentioned above, we see now
are remnants of that original magnetic field. Hellas and Argyre, two large impact
basins on Mars are about four billion years old and are demagnetized. If the
dynamo was still operating when those impact features formed, the crust would
have re-magnetized as they cooled. Hence, the dynamo must have stopped before
then.
Chapter 11
Asteroids, Comets,
Meteroites
11.1 Asteroids
11.1.1 General Properties
On the first day of January 1801, Giuseppe Piazzi discovered an object which he
first thought was a new comet. But after its orbit was better determined it was
clear that it was not a comet but more like a small planet and it was therefore
named asteroid. The proper name of the first asteroid detected is Ceres. Three
other small bodies were discovered in the next few years (Pallas, Vesta, and Juno).
By the end of the 19th century several hundred asteroids were known.
Several thousand asteroids have been discovered and given provisional desig-
nations so far. Thousands more are discovered each year. There are undoubtedly
hundreds of thousands more that are too small to be seen from the Earth. There
are 26 known asteroids larger than 200 km in diameter. About 99 % of all objects
> 100 km are known however of the total number of asteroids with diameters
between 10 and 100 km we know only 50%. It is difficult to estimate the total
number of asteroids, perhaps as many as a million 1 km sized asteroids may exist.
Since most of the asteroids have orbits between Jupiter and Mars, it was first
assumed that they are remnants of a larger planet that broke up. However, the
total mass of all the asteroids is less than that of the Moon.
Ceres has a diameter of 933 km, the next largest are Pallas, Vesta and Hygiea
which are between 400 and 525 km in diameter. All other known asteroids are less
than 340 km.
Asteroids are classified into:
C-type: extremely dark (albedo 0.03), similar to carbonaceous chondrite
meteorites; 75% of known asteroids.
S-type: 17% of asteroids; bright (albedo 0.1-0.2); metallic Ni, Fe and Mg
silicates.
201
202 CHAPTER 11. ASTEROIDS, COMETS, METEROITES
Since the albedo for most asteroids is not known, an albedo range between 0.25 to
0.05 is usually assumed. This results in a range for the diameter of the asteroid.
The table 11.4 shows the diameter ranges for an asteroid based on its absolute
magnitude, assuming an albedo ranging from 0.25 to 0.05.
In other words, asteroids that can’t get any closer to the Earth (i.e. MOID)
than 0.05 AU (roughly 7,480,000 km) or are smaller than about 150 m in diameter
(i.e. H = 22.0 with assumed albedo of 13%) are not considered PHAs. The current
list of PHAs is obtained from the Minor Planet Center on a daily basis. Asteroids
with a small MOID to Earth should be carefully followed because they can become
Earth colliders.
Because of long-range planetary gravitational perturbations and, particularly,
close planetary approaches, asteroid orbits change with time. Consequently, MOID
also changes. As a rule of thumb, MOID can change by up to 0.02 AU per century,
except for approaches within 1 AU of massive Jupiter, where the change can be
larger. Thus, an asteroid that has a small MOID with any planet should be
monitored. Currently there are about 350 known PHA’s.
11.4 NEOs
Near-Earth Objects (NEOs) are comets and asteroids that have been nudged by
the gravitational attraction of nearby planets into orbits that allow them to enter
the Earth’s neighborhood. Composed mostly of water ice with embedded dust
particles, comets originally formed in the cold outer planetary system while most
of the rocky asteroids formed in the warmer inner solar system between the orbits
of Mars and Jupiter. The scientific interest in comets and asteroids is due largely
to their status as the relatively unchanged remnant debris from the solar system
formation process some 4.6 billion years ago. The giant outer planets (Jupiter,
Saturn, Uranus, and Neptune) formed from an agglomeration of billions of comets
and the left over bits and pieces from this formation process are the comets we see
today. Likewise, today’s asteroids are the bits and pieces left over from the initial
agglomeration of the inner planets that include Mercury, Venus, Earth, and Mars.
As the primitive, leftover building blocks of the solar system formation process,
comets and asteroids offer clues to the chemical mixture from which the planets
formed some 4.6 billion years ago. If we wish to know the composition of the
primordial mixture from which the planets formed, then we must determine the
chemical constituents of the leftover debris from this formation process - the comets
and asteroids.
In terms of orbital elements, NEOs are asteroids and comets with perihelion
distance q less than 1.3 AU. Near-Earth Comets (NECs) are further restricted to
include only short-period comets (i.e orbital period P less than 200 years). The
vast majority of NEOs are asteroids, referred to as Near-Earth Asteroids (NEAs).
NEAs are divided into groups (Aten, Apollo, Amor) according to their perihelion
distance (q), aphelion distance and their semi-major axes (a).
Possible NEO missions that require spacecraft with the capability to ren-
11.5.. THE CRETACEOUS-TERTIARY IMPACT 205
dezvous at great distances (1 AU) from the Earth within a releatively short amount
of time (on the order of a year) are discussed by Sforza and Remo (1997) and Pow-
ell et al. (1997). NEOs as near Earth resources for mining are discussed by Gertsch
et al. (1997).
craters,
earthquakes,
tsunami
from a stone asteroid is negligible if it is less than 200 m in diameter. Small and
relatively frequent impactors such as Tunguska produce only air blast damage and
leave no long term scars. Objects 2.5 times larger which hit every few thousand
years cause coherent destruction over many thousand km of coast. Let us assume
an asteroid > 200 m hits an ocean. A water wave generated by such an impactor
has a long range because it is two-dimensional; its height falls off inversely with
distance from the impact. When the wave strikes a continental shelf, its speed
decreases and its height increases to produce tsunamis. Tsunamis produce most
of the damage from asteroids between 200 m and 1 km. An impact anywhere in
the Atlantic by an asteroid 400 m in diameter would devastate the coasts on both
sides of the ocean by tsunamis over 100 m high. An asteroid 5 km in diameter
hitting the mid Atlantic would produce tsunami that would inundate the entire
East coast of the US to the Appalachian mountains (see the paper of Hill and
Mader, 1997).
In Fig.11.1 the estimated frequency of impacts as a function of asteroid diam-
eter is shown.
An asteroid or comet in diameter (which releases or Mt
TNT) would cause a global catastrophe:
smoke clouds,
Some more details about such scenarios can be found e.g. in Chapman and Mor-
rison (1994) or Melosh et al. (1990).
During the last 500 Million years there occurred several extinctions of marine
species: in Table 11.1 we give the formal end of stage in Myr.
From the data we can deduce that the Earth should be hit by several asteroids
and comets larger than a few km ( energy release) and perhaps one
in a period of ~ 100 Myr. Thus for the last 500 Myr years 5 events of extinctions
are to be expected and 20 minor events which is in agreement with astronomical
predictions.
From normal meteoroid ablation an iridium anomaly is observed and is one
of the most significant signatures for impact. The search for iridium has resulted
in reports of elevated iridium levels ( times background values) at or near a
number of extinction boundaries. The Ir levels are generally significantly weaker
than the K/T anomaly.
The conversion from projectile mass to crater size was given by Shoemaker et
al. (1990):
11.6 Meteorites
The term meteor comes from the Greek meteoron, meaning phenomenon in the
sky. Meteors are small solid particles that enter the Earth’s atmosphere from
interplanetary space. They move at high speeds and the friction they encounter
in the air vaporizes them (typically at heights between 80 and 110 km above the
surface). The light caused by the luminous vapors formed in such an encounter
appears like a star moving rapidly across the sky, fading within a few seconds. To
be visible, a meteor must be within 200 km of the observer. The total number of
meteors bright enough to be visible is estimated to be about 25 million per day.
A meteoroid is matter revolving around the sun or any object in interplanetary
space that is too small to be called an asteroid or a comet. Even smaller particles
are called micrometeoroids or cosmic dust grains, which includes any interstellar
material that should happen to enter our solar system. A meteorite is a meteoroid
that reaches the surface of the Earth without being completely vaporized.
One of the primary goals of studying meteorites is to determine the history
and origin of their parent bodies. Several achondrites sampled from Antarctica
since 1981 have conclusively been shown to have originated from the moon based
on compositional matches of lunar rocks obtained by the Apollo missions of 1969-
1972. Sources of other specific metorites remain unproven, although another set of
eight achondrites are suspected to have come from Mars. These meteorites contain
atmospheric gases trapped in shock melted minerals which match the composition
of the Martian atmosphere as measured by the Viking landers in 1976. All other
groups are presumed to have originated on asteroids or comets; the majority of
meteorites are believed to be fragments of asteroids.
A typical bright meteor is produced by a particle with a mass less than 1 g.
A particle the size of a golf ball produces a bright fireball. The total mass of
meteoritic material entering the Earth’s atmosphere is estimated to be about 100
11.6. METEORITES 209
Space Debris
Orbital debris is defined as any man-made object in orbit around the Earth which
no longer serves as a useful purpose.
In 1957 Sputnik 1 was launched as the first man made spacecraft. In the years
of space activities some 3 750 launches led to more than 23 000 observable space
objects (larger than 10 cm) of which currently 7 500 are still in orbit. Only 6% of
the catalogued orbit population comprise operational spacecraft, while 50% can be
attributed to decommissioned satellites, spent upper stages, and mission related
objects (launch adapters, lens covers, etc.). The remainder of 44% is originating
from 129 on-orbit fragmentations which have been recorded since 1961. These
events, all but 1 or 2 of them explosions of spacecraft and upper stages, are assumed
to have generated a population of objects larger than 1 cm on the order of 70 000 to
120 000. Only at sizes of in the range of 0.1 mm the sporadic flux from meteoroids
prevails over man-made debris. From a statistical point of view we have to note
that most orbital debris reside within 2 000 km of the Earth’s surface. Within
this volume, the amount of debris varies significantly with altitude and regions of
debris concentration are found near 800 km, 1 000 km and 1 500 km.
211
212 CHAPTER 12. SPACE DEBRIS
On Dec 3, 2001 BBC reports, that space debris lit up the sky. The spectacular
nighttime light show was seen over parts of southern England is now believed to
have been caused by burning Russian space debris. Observers said the fragments,
which could be seen over parts of Essex and Sussex, were very bright and traced
across the sky for up to four minutes.
Of course the International Space Station (ISS) will be the most heavily shielded
spacecraft ever flown. Critical components (e.g., habitable compartments and high
press tanks) will normally be able to withstand the impact of debris as large as
1 cm in diameter. ISS will also have manoeuvering capability to avoid hazardous
objects.
12.1. REENTRY OF ORBITAL DEBRIS 213
Impacts on spacecraft: all spacecraft collide with very small orbital debris
particles and meteoroids. The Long Duration Exposure Facility (LDEF) was
a bus sized spacecraft. It was returned after 5.7 years in low Earth orbit.
The LDEF was placed in low Earth orbit (LEO) by the space shuttle Chal-
lenger in April 1984 and retrieved by the space shuttle Columbia in January
1990. On the LDEF over 30 000 impacts were found (these craters were
visible to the naked eye and larger than 0.5 mm). Form that sample about
1000 were chemically analyzed in order to investigate the origin of the pro-
jectiles. The largest crater found on LDEF had a diameter of 5 mm and
was probably caused by a particle of 1 mm. Some impacts were clustered in
time. On the European Retrievable Carrier (EURECA), the largest impact
crater diameter was 6.4 mm. The returned solar array of the HST (Hubble
Space Telescope, NASA/ESA) had been the one with the highest orbit al-
titude. It was found that the impact flux for HST was considerably higher
(factor 2-8) than for EURECA. The infra-red astronomical satellite (IRAS),
launched in 1983 to perform a sky survey at wavelengths ranging from 8 to
120 was operational during 10 months near altitude of 900 km. 200 000
potential debris sightings are stored in a database. About 10 000 sightings
are attributed to real objects. A plot of debris flux in low Earth orbit as a
function of object size (cm) is given in Fig. 12.2 where the coordinates are
logarithmic.
in orbit. They measured about 100 possible impact sites and rated them by
size on an arbitrary scale of 1 to 10 (10 being the largest). They found 14
impact craters with a diameter greater than 450 microns.
At NASA a hypervelocity impact technology facility is under operation (HITF).
since radars must provide their own illumination. Therefore, optical tele-
scopes of modest size are more suitable than most radars for detection of
debris at high altitudes. On the other hand, radars are better suited to
detect objects in LEO.
12.4.2 Telescopes
Space debris can be categorized into objects that reflect radar well but sunlight
poorly. The other group reflects sunlight well but radar poorly. Thus, radar
and optical telescopes see somewhat different debris populations. With the use
of optical telescopes, debris at very high altitudes (e.g. in geosynchronous orbits,
GEO) can be detected.
The US Space Command employs aperture telescopes of 1 m to track HEO
objects. With these telescopes objects of 1 m at geosynchronous altitudes, corre-
sponding to a limiting stellar magnitude of 16 can be detected. A limiting stellar
magnitude of 17 or greater is needed to detect debris smaller than 1 m near GEO.
Most objects in GEO are intact; in 1978 a Russian Ekran satellite in GEO was
observed to explode.
NASA is using two optical telescopes for measuring orbital debris: a 3 m
diameter liquid mirror telescope which is referred to as the LMT, and a charged
coupled device-equipped 0.3 m Schmidt camera, which is commonly referred to as
the CCD Debris Telescope or CDT. The LMT consists of a 3 m diameter parabolic
dish that holds four gallons of liquid mercury. The dish is spun up to a rate of 10
revolutions per minute. Centrifugal force and gravity cause the mercury to spread
out in a thin layer over the dish creating a reflective parabolic surface that is as
good as many polished glass mirrors.
12.4.3 Catalogues
There are two catalogues of space objects that are frequently updated:
United States Space Command catalogue,
Space Object catalogue of the Russian Federation.
Based on those two catalogues data are also archived in the Database and Infor-
mation System Characterizing Objects in Space (DISCOS) of ESA. The National
Space Development Agency (NASDA) of Japan is studying a debris database.
Current catalogues contain information on satellites and debris as small as 10-30
cm in diameter. Some recent activities are aimed to provide detection of 5 cm
objects at altitudes below 600 km. For smaller sizes modelers must use statistical
measurements.
mechanism for satellites in GEO. One can estimate an annual collision probability
for an average operational satellite with other catalogued objects at
Another problem concerns the re-entry. Since the last 40 years 16 000 known
re-entries of catalogued space objects are known. No significant damage or injury
occurred which can be attributed to the large expanse of ocean surface and sparse
population density in many land regions. During the past years, approximately
once each week an object with a cross section of or more entered the Earth’s
atmosphere. The risk of re-entry comes from:
Mechanical impact,
chemical contamination,
radiological contamination.
Since about 12% of the present catalogued space debris population consists
of objects discarded during normal satellite deployment (fasteners, yaw, weights,
nozzle covers, lens caps, tethers,...) one should take mitigation measures against
these objects. 85% of all space debris larger than 5 cm result from fragmentation
of upper stages. In 1996 the French CERISE spacecraft was struck and partially
disabled by an impact fragment which most probably came from an exploded
Ariane upper stage.
12.4.5 Shielding
Protection against particles 0.1-1 cm size can be achieved by shielding spacecraft
structures. Objects 1-10 cm in size cannot be shielded nor can they be rou-
tinely tracked by surveillance networks. Protection against these particles can be
achieved through special features in the design (e.g. redundant systems, frangible
structures...). Physical protection against particles larger 10 cm is not technically
feasible. In front of the spacecraft wall single sheet Whipple bumbers or complex
layers of metal and ceramic/polymer fabrics can be used for shieldings. They break
up the impacting particle and absorb the energy of the resulting ejecta. Bumper
shields should be positioned at a sufficient distance from the shielded object.
The penetration depth (damage potential) of an impacting object depends on:
mass,
velocity,
220 CHAPTER 12. SPACE DEBRIS
There are several models (NASA: BUMPER, ESA: ESABASE, Russia: BUFFER
and COLLO).
For manned spacecraft shield designs offer protection against objects smaller
than 1 cm. The PNP (probability of penetration) is an important criterion for
shield design. One can also install automatic detection systems to locate damage.
For EVA (extravehicular activities) current spacesuits have many features with
inherent shielding qualities to offer protection from objects of sizes up to 0.1 mm.
By properly orientating their spacecraft, astronauts may also be able to use their
vehicles against the majority of space debris or direct meteoroid streams. The
United States Space Surveillance Network (SSN) and the Russian Space Surveil-
lance System(SSS) monitor the LEO environment to warn crewed spacecraft if an
object is projected to approach within a few km. If an object is predicted to pass
through a box of 5 × 225 × 5 km oriented along the flight path of the United States
Space Shuttle, the SSN sensor intensifies its tracking of the potential risk object. If
the improved fly-by prediction indicates a conjunction within a box of 2 × 5 × 2 km
an avoidance manoeuvre is performed. During 1986-1997 4 such evasive manoeu-
vres were executed. Collision avoidance manoeuvres were performed by the ESA
satellite ERS-1 in June 1997 and March 1998 and by the CNES satellite SPOT-2
in July 1997.
Calculations made prior to the launch of spacecrafts permit the establishment
of safe launch windows.
For unmanned spacecraft, lower PNPs are tolerable.
12.5. LITERATURE 221
12.5 Literature
An overview of fragmentations of LEO Upper Stages was given by Chernyavskiy
et al. (1994).
A technical report on space debris was given from the Scientific and Technical
Subcommittee of the United Nations Committee on the Peaceful uses of Outer
Space (1999).
Chapter 13
Appendix
13.1 Bibliography
Adams L., Harboe-Sorensen R., Daly E. and Ward J., 1989, IEEE Trans.
Nucl. Sci. N-S36, 6, 2339 (1989).
Adams J.H., Jr., 1986, NRL Memorandum Report 5901 (1986).
Alley, R.B., Clark, P. U., 1999, Annual Review Of Earth And Planetary
Sciences, Volume 27, 149
Alvarez, L.W., Alvaraez, W., Asaro, F., Michel, H.V., 1980, Science, 208,
1095
Andersson, L., Eliasson, L., Wintoft, P., 1999, in Proceedings of the Work-
shop on Space Weather, 11-13 November 1998, ESTEC, Noordwijk, The
Netherlands, WPP-155
Angell, J.K., Korshover, J., 1976, Month Weather Rev., 101, 426-443
Arnaud, J., Mein, P., Rayrole, J., 1998, in Crossroads for European Solar
and Heliospheric Physics. Recent Achievements and Future Mission Possibil-
ities, Puerto de la Cruz, Tenerife, Canary Islands, Spain. March 23-27,1998.
Organised by ESA Solar Physics Planning Group (SPPG) and Instituto de
Astrofisica de Canarias (IAC). European Space Agency, 1998., 213
Arpesella, C., et al. BOREXINO Prop. Vols. 1,2 ed. R. Raghavan et al.,
Univ. of Milano, Milano, 1992
Aschwanden, Markus J.; Poland, Arthur I.; Rabin, Douglas M., 2001 Annual
Review of Astronomy and Astrophysics, 39, 175
Badhwar, G.D., 1997, Radiat. Res., 148, S3
Bahcall, J. N.; Pinsonneault, M. H., 1992, Rev. Mod. Phys., Volume 64,
885
223
224 CHAPTER 13. APPENDIX
Bahcall, J.N., Davis, R., Jr. , 2000 PASP, 112, Issue 770, 429
Bahcall, J.N., in Proc. of the 2nd int. workshop on Low Energy Solar
Neutrinos, Univ. of Tokyo, Tokyo, Japan, World Scientific
Bahcall, J.N., Pinonncault, M.H., Basu, S., Christensen-Dalsgaard, J., 1997,
Phys. Rev. Lett. 78, 171
Baliunas, S., Soon, W.H., 1995, Astrophys. J., 450, 896
Barnes J.R., Collier Cameron A., James D.J., Donati J.F. , 2001 MNRAS
324, 231
Basher, R. E. Basic science of solar radiation and its ultraviolet components,
Proceedings of the seminar on solar ultraviolet radiation, R. Basher (Ed),
NZMS, 1981.
Basu, D., 1992, Solar Physics, 142, 205
Basu, D., 1998, Sol. Phys. 183, 291
Basu, Sarbani, Christensen-Dalsgaard, J., Schou, J., Thompson, M. J., Tom-
czyk, S., 1996, Astrophysical Journal, 460, 1064
Beckers, J.M., and Brown, T.M. 1978, Oss. Mem. d. Oss. Astrofis. d.
Arcetri, 106, 189
Beer, J., 2000, Space Science Rev., 94, 53
Benz A. O., 1993 Plasma Astrophysics, Kinetic Processes in Solar and Stellar
Bradley, P.D., Normand, E. , 1998, , IEEE Trans. Nucl. Sci., 45, no. 6
Coronae (Dordrecht: Kluwer)
Berger, A., 1980, Vistas in Astron., 24(2),103-122
Berner, R.A., Lasaga, A.C., Garrels, R.M., 1983, Amer. J. Sci. 283,641-683
Berry, P.A.M., Vistas in Astron., 30, 97-108
Birn, J., Hones, E. W., Jr., 1981, Journal of Geophysical Research, 86, 6802
Boberg, F., Wintoft, P., Lundstedt, H., 2000, Physics and Chemistry of
Earth, Vol. 25, 4
Bogdan, T.J.,2000, Solar Physics, 192, 373
Bohor, B.F., 1990, in Global catastrophes in Earth history, V.L. Sharpton,
P.D. Ward Eds., 335, Geolo. Society of America, Special Paper 247
Bochsler, P., 2001, Space Science Review, 97, 113
Boteler, D.H., R.J. Pirjola and H. Nevanlinna, 1998 Advances in Space Re-
search, 22, 17
13.1 BIBLIOGRAPHY 225
Bourgeois, J., Hansen, T.A., Wiberg, P.L., Kauffman, E.G., 1988, Science
241,567
Burlaga, L. F., E. Sittler, F. Mariani, and R. Schwenn,1981, J. Geophys.
Res., 86, 6673, 1981.
Burlaga, L.F., 2001, Planetary and Space Science, 49, 1619
Calvo, R. A., Ceccato, H. A., Piacentini, R. D., 1995, Astrophysical Journal,
444, 916 Publication Date:
Canuto, V.M., Levine, J.S., Augustsson, T.R., Imhoff, C.L., Nature, 296,
816-820
Catling, D.C., Zahnle, K.J., McKay, C.P., 2001, Science, 293, 839
Cattaneo, F., Lenz, D., Weiss, N., 2001, The Astrophysical Journal, 563,
L91-L94
Cattaneo, F., in SCORe’96 : Solar Convection and Oscillations and their
Relationship, Proceedings of a workshop, held in Aarhus, Denmark, May 27
- 31, 1996, Eds.: F.P. Pijpers, J. Christensen-Dalsgaard, and C.S. Rosenthal,
Kluwer Academic Publishers (Astrophysics and Space Science Library Vol.
225),201
Cattaneo, F., 1999, in Motions in the Solar Atmosphere, Proceedings of the
Summerschool and Workshop held at the Solar Observatory Kanzelhöhe,
Kärnten, Austria, September 1-12, 1997, edited by Arnold Hanslmeier and
Mauro Messerotti, Kluwer Academic Publishers, Dordrecht/Boston/London,
Astrophysics and Space Science Library, v. 239, 1999, ISBN 0-7923-5507-5,
119
Cauzzi, G.; Falchi, A.; Falciani, R., 2000, Astronomy and Astrophysics, 357,
1093
Cecchini, S., Giacomelli, G., Hasegan, D., Mandrioli, G., Maris, O., Patrizii,
L., Plaian, A., Popa, V., Stefanov, L., Valeanu, V., 2000, Astrophysics and
Space Science, 273, 35
Chandrasekhar, S., 1961, Hydrodynamik and Hydromagnetic Stability, Ox-
ford Univ. Press
Chapman, C.R., Morrison, D., 1994, Nature, 367, 33
Charbonneau, P., and MacGregor, K.B. 1996, The Astrophysical Journal,
473, L59
Chen, J., Space Science Reviews, 95, 165
Chernyavskiy, G., Morozov, N., Johnson, N., McKnight, D., Maclay, T.,1994,
IAA-94-IAA6.5.696
226 CHAPTER 13. APPENDIX
Fisher, G. H.; Fan, Y.; Longcope, D. W.; Linton, M. G.; Pevtsov, A. A.,
2000, Solar Physics, 192, 119
Fleck, B., 2001 The Dynamic Sun, Proceedings of the Summerschool and
Workshop held at the Solar Observatory Kanzelhöhe, Kärnten, Austria, Au-
gust 30 - September 10, 1999, edited by Arnold Hanslmeier, Mauro Messerotti
and Astrid Veronig, Kluwer Academic Publishers, Dordrecht/Boston/London,
Astrophysics and Space Science Library, v. 259, ISBN 0-7923-6915-7, 336 ,
1
Foukal, P., and Lean, J., 1990, Science, 247, 556
Foukal, P.V., Lean, J., 1988, Astrophysical Journal, 328, 347-357
Friis-Christensen, E., Lassen, K., 1991, Science, 254, 698
Fröhlich, C., 2000, Space Science Review, 94,15
Emilio, M., Kuhn, J.R., Bush, R.I., Scherrer, P., 2000, Astrophysical Journal,
543, 1007-1010
Fröhlich, C., 1987, J. Geophys. Res., 92, 796-800
Fröhlich, C., Romero, J., Roth, H., Wehrli, C., Andersen, B. N., Appour-
chaux, T., Domingo, V., Telljohann, U., Berthomieu, B., Delache, P., Provost,
J., Toutain, T., Crommelynck, D., Chevalier, A., Fichot, A., Däppen, W.,
Gough, D. O., Hoeksema, T., Jiménez, Gómez, M., Herreros, J., Roca-
Cortés, T., Jones, A. R., Pap, J.and Willson, R. C.: 1995, Solar Physics,
162, 101
Gaizauskas V 1989 Preflare activity Solar Physics, 121, 135-52
Gerard, J.C., 1990, Phil. Trans. Roy. Soc., London, in Press
Gertsch, R., Remo, J.L., Sour-Gertsch, L., 1997, in Near Earth Objects,
Ann. of the New York Academy of Sciences 822, 468
Gibson, S.E., Space Science Review, 97, 69
Gilliland, R.L., 1981, Astrophysical Journal, 248, 1144-1155
Gilliland, R.L., 1980, Astrophysical Journal, 248, 1144
Gleisner, H., Lundstedt, H.,2001, Journal of Geophysical Research, Volume
106, 8425
Goldreich, P., Murray, N., Kumar, P. 1994, Astrophysical Journal, 424, 466
Goode, P. R., Strous, L. H. Rimmele, Th. R., Stebbins, R. T. , 1998,
Astrophysical Journal Letters 495, L27
Gough, D.O., 1981, Solar Physics, 74, 21-34
228 CHAPTER 13. APPENDIX
Gough, D. O., Kosovichev, A. G., Toomre, J., Anderson, E., Antia, H. M.,
Basu, S., Chaboyer, B., Chitre, S. M., Christensen-Dalsgaard, J., Dziem-
bowski, W. A., Eff-Darwich, A., Elliott, J. R., Giles, P. M., Goode, P. R.,
Guzik, J. A., Harvey, J. W., Hill, F., Leibacher, J. W., Monteiro, M. J. P.
F. G., Richard, O., Sekii, T., Shibahashi, H., Takata, M., Thompson, M. J.,
Vauclair, S., Vorontsov, S. V., 1996, Science, 272, 1296
Graham, G.A., Sexton, A., Grady, M.M., Wright, I.P.,1997, Advances of
Space Research, 20, 1461-1465
Haigh, J.D., 1994, Nature, 370, 544-546
Haigh, 1996, Science, 272, 981
Haisch, B., Strong, K.T., Rodono, M., 1991, Ann. Rev. of Astron. and
Astrophys., 29, 275
Hansen, J.E., Lacis, A.A., 1990, Nature, 346, 713
Hathaway, D. H.; Beck, J. G.; Bogart, R. S.; Bachmann, K. T.; Khatri, G.;
Petitto, J. M.; Han, S.; Raymond, J. ,2000, Solar Physics, 193, 299
Haxton, W.C., 2001, in: Current aspects of neutrino physics. David O.
Caldwell (ed.). Physics and astronomy online library. Berlin: Springer,
ISBN 3-540-41002-3, 2001, p. 65 - 88
Heckman, G.R., 1988, in JPL, Cal. Inst. of Techn., Interplanetary Particle
Environment, proc., N8928454-22-90, 91
Hernandez, J. V., Tajima, T. Horton, W., 1993, Geophysical Research Let-
ters, 20, 23, 2707
Hills, J.G., Mader, Ch.,L., 1997, Near Earth Objects, Ann. New York Acad.
of Sciences, Vol. 822, 381
Hines, C.O., 1974, Atmos. Sci., 31, 589-591
Hoekzema, N. M., Brandt, P. N., Rutten, R. J., 1998, Astronomy and As-
trophysics, 333, 322
Hoeksema, J. T. ,Schou, J. , Saba, J. L. R. , Tarbell, T. D. , Title, A. M. ,
Wolfson, C. J. Milford, P. N. 1997, Solar Physics, 170, 63.
Holland, H.D., 1978, The Chemistry of the Atmosphere and the Oceans,
New York, Wiley
Hood, L.L., J. Geophys. Res., 91; 5264-5276
Hood, L.L., J. Geophys. Res., 92, 876-888
Houghton, G.J. Jenkins, and J.J. Ephraums (eds.),Climate Change: The
IPCC Scientific Assessment. Cambridge University Press, Cambridge, UK.
13.1. BIBLIOGRAPHY 229
Hoyt, D.V., and K. H. Schatten, 1997. The Role of the Sun in Climate
Change, Oxford University Press, 279 pp.
Hoyt, D.V., Schatten, K.H., 1993, J. Geophys. Res., 98, 18895
Hughes, M.K. and H.F. Diaz, 1994, Climatic Change, Volume 26, 109
Hurlburt, N.E., 1999, in SOHO-9 Workshop “Helioseismic Diagnostics of
Solar Convection and Activity”, Stanford, California, July 12-15
Ishii, T. T.; Kurokawa, H.; Takeuchi, T. T. , 2000, in The Sun and Space
Weather, 24th meeting of the IAU, Manchester, England
Johnson, S.J., Dansgaard, W., Clausen, H.B., Langway, C.C., 1970, Nature
227, 482-483
Kalkofen, W., 1997, Astrophysical Journal Lett. 436,L145
Kalkofen, W., 1990, in Basic Plasma Processes on the Sun (A92-30901 12-
92). Dordrecht, Netherlands, Kluwer Academic Publishers, 1990, 197
Kalkofen, W.; Keller, C. U.; Smartt, R.; Hasan, S. S., 2000, Astronomy and
Astrophysics, 363, 279
Kasting, 1989, Paleogeog. Paleoclimat. Paleoecol., 75, 83-95
Kasting, J.F., Whittet, D.C.B., Sheldon, W.R., 1997, Origins of Life and
Evolution of the Biosphere, v. 27, 413
Keating, G.M., Pitts, M.C., Brasseur, G., De Rudder, A., 1987, J. Geophys.
Res. 92, 889-902
Kerr, R.A., 1988, Science, 242, 124-125
Kiefer, J., 2001, Phys. Med. 17 Suppl 1, 1
Khomenko, E.V., Kostik, R.I., Shchukina, N.G., 2001, Astronomy and As-
trophysics, 369,660
Kneer, F., von Uexküll, M., 1999, in Motions in the Solar Atmosphere,
Proceedings of the Summerschool and Workshop held at the Solar Ob-
servatory Kanzelhöhe, Kärnten, Austria, September 1-12, 1997, edited by
Arnold Hanslmeier and Mauro Messerotti, Kluwer Academic Publishers,
Dordrecht/Boston/London, Astrophysics and Space Science Library, v. 239,
1999, ISBN 0-7923-5507-5, 99
Kosovichev A. G. et al., 1998, in Sounding solar and stellar interiors, Pro-
ceedings of the 181st symposium of the International Astronomical Union,
held in Nice, France, September 30 - October 3, edited by Janine Provost,
Francois-Xavier Schmider. 1996 Publisher: Dordrecht Kluwer Academic
Publishers, 1998 438 p.
230 CHAPTER 13. APPENDIX
Lydon, T.J., Sofia, S., 1995, Astrophysical Journal Supplement, 101, 357
McBride, N., McDonnell, J.a.m., 1999, Planetary and Space Science, 47,
1005
Mc Cormac, B.N., Seliga, T.A., eds., 1979, Solar-Terrestrial influences on
Weather and Climate, Dordrecht,: D. Reidel.
Mc Donald, A.B.,1991, in Proc. of the 9th Lake Luise Winter Inst., A.
Astbury et al. (eds), World Scientific,1
McDonnell, J.A.M., Gardner, D.J., Icarus, 133, 25
McDougall, J.D., 1988, Science, 239, 485
McGregor, K.B., and Charbonneau, 1997, Astrophysical Journal, 486, 484
McKay, G.F., Dubeau, J., Thomson, I., 1993 in Proceedings of Spacebound
93 Conference, Canadian Space Agency, Ottawa, Canada, 11
McKenzie, R. L., and J. M. Elwood, Intensity of solar ultraviolet radiation
and its implications for skin cancer,1990, NZ Medical Journal, 103, 152 1990.
McLean, D. J. and Labrum, N. R. ed 1985 Solar Radiophysics (Cambridge
Univ. Press)
McIntosh,1990, Sol. Phys. 125, 251
Mann, M.E., Bradley, R.S., and Hughes, M.K., Geophys. Res. Lett., 1999,
26, 759
Markvart, T., Cumberbatch, T.J., Dollery, A.A., Walkden, M., Photon and
electron degradation of boron-doped FZ silicon solar cells, Proc. Third Eu-
ropean Symposium on Photovoltaic Generators in Space (1982), p.109
Markvart, T., Willoughby, A.F.W., Dollery, A.A., Radiation-resistant silicon
solar cell, Proc. 19th IEEE Photovoltaic Specialists Conference (1987), p.
709 J.W. Peters, T.
Markvart, T., Willoughby, A.F.W., Goodbody, G.C., Dollery, A.A., Defect
interactions in silicon solar cells, Second Workshop on Radiation-Induced
and/or Process-Related Electrically Active Defects in Semiconductors-Insulator
Systems, North Carolina, U.S.A. (1989).
Markvart, T., Radiation Damage in Solar Cells: a review, Journal of Material
Science: Materials in Electronics, Vol. 1, No. 1 (1990)
Marsh, N. D., Svensmark, H., 2000, Physical Review Letters 85, 5004
Mein, P., Rayrole, J., 1997, in 1st Advances in Solar Physics Euroconference.
Advances in Physics of Sunspots, ASP Conf. Ser. Vol. 118., Eds.: B.
Schmieder, J.C. del Toro Iniesta and M. Vazquez, p. 334.
232 CHAPTER 13. APPENDIX
Melosh, H.J., Schneider, N.M., Zahnle, K.J., Latham, D., 1990, Nature, 343,
251
Morel, P., Provost, J., Berthomieu, G., 1997, Astronomy and Astrophysics,
327, 349
Moreno-Insertis, F. ,1994, in Solar Magnetic Fields, eds. M. Schüssler and
W. Schmidt, (Cambridge: Cambridge Univ. Press), 117
Muller, R., 1999, in Motions in the Solar Atmosphere, Proceedings of the
Summerschool and Workshop held at the Solar Observatory Kanzelhöhe,
Kärnten, Austria, September 1-12, 1997, edited by Arnold Hanslmeier and
Mauro Messerotti, Kluwer Academic Publishers, Dordrecht/Boston/London,
Astrophysics and Space Science Library, v. 239, 1999, ISBN 0-7923-5507-5,
p. 35-70
Narain, U.; Ulmschneider, P. , 1990, Space Science Reviews, vol. 54, 377
National Research Council, Space Studies Board, 2000, Radiation and the
International Space Station: Recommendations to reduce risk. Washington,
DC, National Academy Press
Neckel, H., 1995, Sol. Phy. 156,7
Neftel, A.; Moor, E.; Oeschger, H.; Stauffer, B. , 1985, Nature, 315, 45
Noel, F., 1997, Astronomy and Astrophysics, 325, 825
Normand, E., IEEE Trans. Nucl. Sci.,1996 43, 461
November, L.J., Toomre, J., Gebbie, K.B., Simon, G., 1981, Astrophysical
Journal, 245, L123
Ribes, E., Ribes, J.C., Barthalot, R., Nature, 326, 52-55
Ribes, E., Beardsley, B., Brown, T.M.,Delach, Ph., Laclare, F., Leister, V.N.,
1991, in The Sun in Time, ed. C.P. Sonett, M.S. Giampapa, M.S. Matthews
(Tucson Ariz. Univ. Press), 59
Otto, A., 1995, Rev. Geophys. Vol. 33 Issue S1, 657
Palle B., E., Butler, C.J., 2000, ESA SP-463
Perry, C.A., 1994, International Journal of Climatology, 14, 969
Peters, J.W., Markvart, T., Willoughby, A.F.W., 1992, A study of radiation-
induced defects in silicon solar cells showing improved radiation resistance,
Materials Science Forum 83-87 1539
Pirjola, R., Boteler, D., Viljanen, A. and Amm, O., 2000 , Advances in Space
Research, 26, 514.
Potgieter, M.S., 1998, Space Science Rev. 83, 147
13.1. BIBLIOGRAPHY 233
Powell, J., Maise, G., Ludewig, H., Todosow, M., 1997, in Near Earth Ob-
jects, Ann. of the New York Academy of Sciences 822, 447
Priest, E., 2000, in The Sun and Space Weather, 24th meeting of the IAU,
Manchester, England
Raup, D., 1989, Phil. Transactions of the Royal Society of London..., 325,
421
Reedy, R.C., Second Conference on the High Energy Radiation Background
in Space, Snowmass, Co, 22-23 July 1997, 41
Reid, G.C., 1991, J. Geophys. Res. 96, 2835
Riklis, E., Emerit, I., Setlow, R.B.,1996, Advances in Space Research, Vol.
18, 51
Rieutord, M.; Roudier, T.; Malherbe, J. M.; Rincon, F., 2000, Astronomy
and Astrophysics, 357, 1063
Roald, Colin B.; Sturrock, P. A.; Wolfson, Richard, 2000, Astrophysical
Journal, 538, 960
Roudier, T., Malherbe, J. M., November, L., Vigneau, J., Coupinot, G.,
Lafon, M., Muller, R., 1997, Astronomy and Astrophysics, 320, 605.
Roxburgh, I.W., 1998, Astrophysics and Space Science, 261, 57
Ruderman, M.A., Chamberlain, J.W., Planet. Space Sci., 1975, 23, 247-268
Russell, C.T., 2001, Planetary and Space Science, 49, 1005
Rutten, R. J., Hammerschlag, R. H., Bettonvil, F. M., Suetterlin, P., 2000,
American Astronomical Society, SPD meeting 32, 02.107
Ryutova, M. P.; Tarbell, T. D., 2000, American Astronomical Society, SPD
meeting 32, 01.41
Sanchez, M., Parra, F., Soler, M., Soto, R., 1995, Astronomy and Astro-
physics Suppl. 110, 351
Sello, S., 2001, Astronomy and Astrophysics, 377, 312
Schrijver, C.J., Zwaan, C. , 2000, Solar and Stellar Magnetic Activity Cam-
bridge U. Press, New York, 2000.
Schwartzmann, D.W., Volk, T., 1989, Nature, 340, 457-460
Seckmeyer, G., McKenzie, R.L., 1992, Nature, 359, 135
Shibata, K., 1996, Adv. Space Res 17, 9-18
Settecerri 2, M.J. Matney 2, J.C. Zhang 2 , Haystack Radar Measurements
of the Orbital Debris Environment, 1990-1994, JSC - 27436
234 CHAPTER 13. APPENDIX
Sforza, P.M., Remo, J.L., 1997, in Near Earth Objects, Ann. of the New
York Academy of Sciences 822, 432
Simpson, J.A., 1998, Space Science Rev., 83, 7
Shinagawa, H., 2000, Advances in Space Research, 26, 1599
Shoemaker, E.M., Wolfe, R.F., Shoemaker, C.S., 1990, Geol. Soc. of Amer-
ica, Spec. Paper 247, 155
Snodgrass, H. B., 2001, American Astronomical Society Meeting 198, 71.02
Sobotka, M., 1999, in Motions in the Solar Atmosphere, Proceedings of the
Summerschool and Workshop held at the Solar Observatory Kanzelhöhe,
Kärnten, Austria, September 1-12, 1997, edited by Arnold Hanslmeier and
Mauro Messerotti, Kluwer Academic Publishers, Dordrecht/Boston/London,
Astrophysics and Space Science Library, v. 239, 1999, ISBN 0-7923-5507-5,
, 71
Sofia, S., Okeefe, J., Lesh, J.R., Endal, A.S., 1979, Science, 204, L306
Sonett, C.P., Williams, G.E., 1985, J. Geophys. Res., 90, 12019-12026
Sonett, C.P., 1982, Geophys. Res. Let., 9, 1313-1316
Sofia, S., Heaps, W., Twigg, L.W., 1994, ApJ, 427, 1048
Song, Q., Cao, W., 1999, in Obs. Astrophys. in Asia and its Future, Ed.
P.S. Chen, Yunan Obs. Chin. Acadm. Sc., 1999, 139
Solanki, S.K., 1997, 1st Advances in Solar Physics Euroconference. Advances
in Physics of Sunspots, ASP Conf. Ser. Vol. 118., Eds.: B. Schmieder, J.C.
del Toro Iniesta, M. Vazquez, p. 178
Solanki, S.K., Fligge, M., 2000, ESA SP-463, 51
Spruit, H.C., 1982, Astronomy and Astrophysics, 108, 348
Stauffer, B., 2000, Space Science Reviews, 94, 321
Stein, Robert F.; Nordlund, Å , 2000, Solar Physics, 192, 91
Strassmeier K.G. and Linsky J.L., 1996, Stellar surface structure, conference
proceedings of IAU symposium 176, Kluwer Academic Publishers, Dordrecht,
The Netherlands
Strous, L. H., Goode, P. R., Rimmele, Th. R., 2000, Astrophysical Journal,
535, 1000
Suzuki, Y., 1998, Space Science Reviews, 85, 91
Svensmark, H., Friis-Christensen, E., 1997, Journal of Atmospheric Solar-
Terrestrial Physics 59, 1225
13.1. BIBLIOGRAPHY 235
Walker, J.C.G., Hays, P.B., Kasting, J.F., 1981, J. Geophys. Res. 86, 9776-
9782
Walsh, R.W., 2001 The Dynamic Sun, Proceedings of the Summerschool
and Workshop held at the Solar Observatory Kanzelhöhe, Kärnten, Aus-
tria, August 30 - September 10, 1999, edited by Arnold Hanslmeier, Mauro
Messerotti and Astrid Veronig,
Kluwer Academic Publishers, Dordrecht/Boston/London,
Astrophysics and Space Science Library, v. 259, ISBN 0-7923-6915-7, 336
Watanabe, S., Shibahashi, H., 2001, Publications of the Astronomical Society
of Japan, vol.53, no. 3, p. 565
Walther, G., 1999, Astrophysical Journal, 513, 990
Willson, R. C., 1981, Solar Physics, 74, 217-229.
Willson, R. C., 1984, Space Science Reviews, 38, 203-242.
Willson, R.C., Hudson, H., Nature, 332, 810-812
Willson, R.C., Hudson, H.S., Fröhlich, C., Brusa, R.W., Science, 1987, 234,
1114-1117
Wittmann, A.D., Alge, E., Bianda, M., 1993, Solar Physics, 145, 205
Wu, J.-G., Lundstedt, H., Eliasson, L., Hilgers, A., 1999, in Proceedings of
the Workshop on Space Weather, 11-13 November 1998, ESTEC, Noordwijk,
The Netherlands, WPP-155
Yang T.C., Mei, M., George, K.A., Craise, L.M., 1996, Advances of Space
Research, 18, 149
Zahnle, K.J., Walker, J.C. G., 1992, Rev. Geophys. Space Phys., 20, 280-292
Ziegler, J.F., Lanford, W.A., 1979, Science, 20, 776
Zwickl, R. D.; Doggett, K. A.; Sahm, S.; Barrett, W. P.; Grubb, R. N.;
Detman, T. R.; Raben, V. J.; Smith, C. W.; Riley, P.; Gold, R. E.; Mewaldt,
R. A.; Maruyama, T., 1998, Space Science Rev., 86, 633
13.2 Internet
Today’s space weather can be found under:
http://www.sel.noaa.gov/today.html
The web site of the National Oceanic and Atmospheric Administration:
http://www.sel.noaa.gov/
ESA Space Weather Site:
13.2. INTERNET 237
http://www.estec.esa.nl/wmwww/spweather/
NASA Space weather resources:
http://spdf.gsfc.nasa.gov/space_weather/Space_Weather_at_SSDOO.html
Space Science Institute/ NASA and NSF site:
http://www.spacescience.org/SWOP/
Lund Space Weather Center
http://www.irfl.lu.se/
239
240 LIST OF TABLES
2.1 A typical spiral galaxy. From a distant galaxy our own system would
appear similar, the Sun would be located in one of the spiral arms 4
2.2 Sketch of the Hertzsprung-Russell-diagram with evolutionary path of the
Sun. 6
4.1 Variation of electron temperature and electron density in the solar atmo-
sphere 26
4.2 Spectroscopic observation of solar granulation. The entrance of a spec-
trograph slit covers different granular/intergranular areas. Line profiles
emanating from granules are blueshifted because matter moves upwards
and profiles from intergranular areas are redshifted because matter moves
away from the observer. This is valid for solar granulation observed near
the disk center 29
4.3 Solar granulation and small network bright points 32
4.4 Large sunspot showing the dark central umbra and the filamentary penum-
bra. Outside the penumbra the granulation pattern is clearly seen 35
4.5 Relative Sunspot number 40
4.6 Butterflydiagram illustrating the equatorward motion of spots during the
activity cycle 41
4.7 Profile of the CaII line 48
4.8 Coronal hole seen by the solar satellite YOHKOH 58
4.9 Comet Hale Bopp (1997); the fainter ion tail is clearly seen 59
241
242 LIST OF FIGURES
5.1 Left: waves with low and high l; the low l modes are reflected deeper than
the high l modes. Right: Explanation how the waves are reflected in the
solar interior. The wavefront (normal to the propagation) is deflected
since the sound velocity is higher in deeper layers 77
5.2 Examples of several modes 78
5.3 l-v diagram from MDI high-cadence full disk data shows mode frequencies
up to 10 mHz and l=1000 80
5.4 This diagram shows the solar rotation rate inferred from two months
of MDI Medium-l data as a function of radius at three latitudes, 0
degrees, 30 degrees, and 60 degrees 83
5.5 MDI Dopplerimage; left: the rotation of the Sun is clearly seen; right:
the rotation of the Sun was eliminated and therefore only velocities due
to granulation and supergranulation are seen 84
6.1 Looped magnetic field lines in the solar chromosphere and corona. Photo:
NASA 88
6.2 Principle of magnetic reconnection 92
6.3 Illustration of the effect. The field lines are wraped around because of
the differential rotation of the Sun 101
6.4 The MHD relation between flows and magnetic fields 102
6.5 a) Cutaway images of solar rotation showing a peak and a trough of the
0.72R variation, with black indicating slow rotation, grey intermediate,
and white fast, b) Variations with time of the difference of the rotation
rate from the temporal mean at two radii deep within the Sun, with the
site at 0.72 located above the tachocline and that at 0.63 below
it, both sampling speeding up and slowing down in the equatorial region.
Results obtained from GONG data for two different inversions are shown
with black symbols, those from MDI with red symbols. (Image courtesy
NSF's National Solar Observatory) 107
10.1 After a flare or coronal mass ejection erupts from the Sun's surface, major
disturbances arrive with a range of time delays and a storm begins to
build in the space surrounding the Earth. 194
10.2 Summary of space weather effect in the Earth's environment 199
11.1 Estimated frequency of impacts on the Earth from the present population
of comets and asteroids and impact craters. 207
11.2 Height in the atmosphere at which half the kinetic energy of a stony
meteoroid is dissipated. Note that asteroids with > 100 m hit the ground
with most of their original kinetic energy 209