Understanding Space Weather and the Sun
Understanding Space Weather and the Sun
ASTROPHYSICS AND
SPACE SCIENCE LIBRARY
VOLUME 347
EDITORIAL BOARD
Chairman
W.B. BURTON, National Radio Astronomy Observatory, Charlottesville, Virginia, U.S.A.
(bburton@[Link]); University of Leiden, The Netherlands (burton@[Link])
Executive Committee
J. M. E. KUIJPERS, University of Nijmegen, The Netherlands
E. P. J. VAN DEN HEUVEL, University of Amsterdam, The Netherlands
H. VAN DER LAAN, University of Utrecht, The Netherlands
MEMBERS
F. BERTOLA, University of Padua, Italy
J. P. CASSINELLI, University of Wisconsin, Madison, U.S.A.
C. J. CESARSKY, European Southern Observatory, Garching bei München, Germany
O. ENGVOLD, University of Oslo, Norway
A. Heck, Strasbourg Astronomical Observatory, France
R. McCRAY, University of Colorado, 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, National Tsing Hua University, Taiwan.
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, Princeton University, U.S.A.
N. O. WEISS, University of Cambridge, U.K.
THE SUN AND SPACE
WEATHER
Second Edition
by
ARNOLD HANSLMEIER
University of Graz, Institute of Physics/ IGAM, Austria
A C.I.P. Catalogue record for this book is available from the Library of Congress.
Published by Springer,
P.O. Box 17, 3300 AA Dordrecht, The Netherlands.
[Link]
Preface xiii
vii
viii CONTENTS
Bibliography 275
Internet 301
Index 305
Preface
The field of solar physics and solar--terrestrial relation, now called space
weather, is evolving rapidly. As in the first edition, it is assumed that it is inevitable for
the reader to get some basic knowledge in solar physics since the Sun is the main driver
for space weather The term space weather itself has been gaining more and more
attention during the past years as our society becomes more and more dependent on
satellites, which are vulnerable to varying conditions in space. Space weather efforts
and investigations are being made all over the world and more and more is known about
the complex relations of processes on the Sun and the Earth and its space environment.
The term space climate nowadays includes the long-term variations caused mainly by
the Sun on the Earth and the interplanetary space.
As in the first edition of the book, this edition also covers these topics but new
chapters have been introduced, e.g., a chapter on real-time space weather forecasts and
some main space weather data sources. All the chapters have updated information,
taking into account the results of new satellite missions and telescopes. The book also
includes a great amount of new literature (more than 340 original citations) so that the
reader is able to go into more details, if required in the respective chapters. Thus, the
book should be helpful to scientists as well as to students interested in overview or
finding a compendium with references to go deeper into special fields.
Furthermore, at the beginning of all the chapters, introductory books are cited,
which could be recommended for the special topics addressed there. The number of
keywords in the index has also been strongly enhanced so that the reader can find
information easily. Besides all this, suggestions from readers of the first edition have
been taken into account and are greatly acknowledged.
I want to thank all my colleagues who provided me with advice and figures and
the students who attended my lectures at Graz and Innsbruck for their help. Last but not
least I thank my family – Karoline, Roland, Christina and Alina – for the patience and
understanding when I spent lots of nights at the computer.
xiii
Chapter 1
In this introduction we briefly describe the term space weather and give motivation
why that interdisciplinary field gained high interest. Examples will demonstrate
the high relevance of space weather not only from the scientific point of view but
from the social and economic aspect of our modern civilization.
Since this is a very modern topic, there appeared several monographs about
that subject, e.g. a collection of space weather related topics1 .
1
2 CHAPTER 1. INTRODUCTION, WHAT IS SPACE WEATHER?
• The Australian IPS Radio and Space Services , the Australian Space Weather
Agency,
and many others (such as the Group in Oulu, Finland). Today, space weather is
monitored from a worldwide net of ground stations and from space. Both types of
observations are complementary. From space the whole electromagnetic spectrum
of the Sun can be observed including UV and X-rays.
An overview about space weather, environment and societies can be found in
the monograph by Lilenstein and Bornarel, 2005 [196].
1.2.1 Examples
Let us give some examples of space weather influences on satellites7
How can we study the propagation of solar disturbances through the interplan-
etary medium from Earth? A common technique is to measure scintillations in the
radio wavelength. Let us consider very distant radio sources like quasars. If the
interplanetary medium is not disturbed, the signal from this object is constant in
amplitude. But similar to the twinkling of starlight in the visible, the radio signal
becomes absorbed and refracted when passing through a plasma cloud emitted
from the Sun. By measuring many point sources distributed all over the sky, one
gets a map of areas of high scintillation which shows where the plasma wave is
propagating.
There are similarities with atmospheric weather, however the most important
differences between atmospheric and space weather systems are:
• 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.
influences such as the eruption of big volcanoes played a role, it is assumed that
during these phases the global climate on Earth was cooler than on the average.
Summarizing, the following branches strongly depend on space weather:
• Communication Systems,
8 see: I.A. Daglis, Effects of Space Weather on Technology Infrastructure, 2004, Kluwer
Chapter 2
Our Sun is the only star which is close enough to observe details on its surface such
as sunspots, faculae, prominences, coronal holes, flares etc., which are all summa-
rized as solar activity phenomena. Therefore, the study of the Sun is important
for astrophysics in general. Theories about stellar structure and evolution can be
studied in detail on the Sun1 .
On the other hand, the Sun is the driving factor for the climate on the Earth
and the structure and shape of the Earth’s magnetosphere thus determining and
influencing the near Earth space environment. Therefore, the study of solar terres-
trial relations is of great importance for our modern telecommunication systems
both based on Earth and in space.
The Cambridge Encyclopedia of the Sun, K.R. Lang, 2001, Cambridge Univ. Press; A Guide to
the Sun, K.H. Phillips, 1995, Cambridge Univ. Press; The Sun, M. Stix, 2002, Springer Verlag
2 1 Ly =1013 km, the distance light travels within one year propagating through space at a
7
8 CHAPTER 2. THE SUN A TYPICAL STAR
Figure 2.1: A typical spiral galaxy. From a distant galaxy, the Sun would be located in
one of the spiral arms. Image: A.H., private observatory.
extrasolar planets. The diameter of our galaxy is about 100 000 Ly. Galaxies are
grouped into clusters- our galaxy belongs to the so called local group of galaxies.
The small and large Magellanic cloud are two small dwarf galaxies which are
satellites of our system. The nearest large galaxy is the Andromeda galaxy which
is at a distance of more than 2 Million Ly.
Many galaxies appear as spiral galaxies. Young bright stars are found in the
spiral arms, older stars in the center and in the halo of the galaxy. An example is
given in Fig.2.1.
• Stellar radii: once the apparent diameter of a star is known than its real
diameter follows from its distance d. The problem is to measure apparent
stellar diameters since they are extremely small. One method is to use
interferometers, one other method is to use occultation of stars by the moon
2.1. THE SUN AND STARS 9
• Stellar masses: can be determined by using Kepler’s third law in the case we
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 rotation causes a
broadening of spectral lines;
-10
He Flash
Giants
MV
M
ai
n
se
qu
en
Present
ce
Sun
White
Dwarfs
+15
50 000 K 5 000 K T
O5 G
thus it depends only on the temperature TS of the object. Here, Iν is the intensity
of radiation at frequency ν; h, k, c are Planck’s constant, Boltzmann’s constant and
the speed of light. h = 6.62 × 10−34 Js−1 , k = 1.38 × 10−23 JK−1 . If that equation
is integrated over all frequencies (wavelengths), we obtain a formula for the total
power emitted by a black body, Boltzmann law:
∞
Bλ dλ = σT 4 , (2.3)
0
L = 4πr2 σTeff
4
(2.4)
For the Sun Teff = 5 785 K. This formula defines the effective temperature of a
star. σ = 5.67 × 10−8 W/m2 K4 is the Stefan Boltzmann constant.
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 70 MW per m2 of its surface3
3 The worldwide nuclear energy generation is about 350 GW. Thus an area of 5000 m2 on the
Table 2.1: Central wavelength and bandwidth of the UBVRI filter set
Name Meaning Central λ Bandwidth [nm]
U Ultraviolet 360 66
B Blue 440 98
V Visual (green) 550 87
R Red 700 207
I Infrared 900 231
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:
T λmax = 2.9 × 10−3 m K (2.5)
This is also called Wien’s law.
At about which wavelength can planets be expected to radiate most of their energy?
Answer: Let us assume the temperature of the Earth = 300 K. Then
λmax = 2.9 × 10−3 /300 ∼ 10 µ (2.6)
The Sun has a surface temperature of about 6 000 K. At what wavelength does
the Sun’s spectrum peak? Answer:
λmax = 2.9 × 10−3 /6000 ∼ 0.5 µ = 500 nm (2.7)
• negative for the hotter star. The hotter star is brighter in B than in V,
therefore for the magnitudes in these two bands: mB < mV and B-V<0.
• Main sequence stars: most stars are found along a diagonal from the upper
left (hot) to the lower right (cool).
• 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).
This leads to the question why most of the stars we observe lie on the Main
sequence. The answer is quite easy: because this denotes the longest phase in
stellar evolution. Let us discuss this briefly for the Sun:
The main steps in the evolution of the Sun are (compare with Fig. 2.2):
• 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 has passed a violent youth, the T Tauri phase.
• At the main sequence the Sun changes extremely slowly remaining there
about 1010 years. 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 the center He burning sets in. The Sun will evolve to a
red giant for some 108 ys. It will extend beyond the Earth’s orbit.
During its evolution, the Sun dramatically changes its radius (the subscript
denotes the present day value):
1 R (present Sun) → ∼ 104 R (red giant), → 0.01R (white dwarf).
For space weather long term evolutionary effects are negligible. But it is in-
teresting to investigate them especially for the early Sun (see the chapter on the
faint young Sun problem).
For the main sequence stars there exists a relation between their mass and
luminosity:
L ∼ M 3.5 (2.9)
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 history4 .
a3 G
2
= (M1 + M2 ) (2.11)
P 4π 2
In our case a denotes the distance Earth-Sun (150 × 106 km), P the revolution
period of the Earth around the Sun (1 year), M1 the mass of the Earth and M2
the mass of the Sun. One can make the assumption that M1 << M2 and therefore
M1 + M 2 ∼ M 2 .
If we know the distance of the Sun and its angular diameter the solar radius is
obtained:
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 W m−2 . UV and IR radiation from the Sun is strongly absorbed by the
Earth’s atmosphere. Therefore, accurate measurements of the solar constant have
to be done with satellites. ACRIM on SMM and ERB on Nimbus 7 showed clearly
that the presence of several large sunspots which are cooler than their surround-
ings depress the solar luminosity by ∼ 0.1%. The Variability IRradiance Gravity
Oscillation (VIRGO) experiment on the SOHO satellite is observing total solar
and spectral irradiances at 402 nm (blue channel), 500 nm (green channel), and
2.2. THE SUN 15
862 nm (red channel) since January 1996 (for a review see e.g. Pap et al. (1999)
[243]). The solar luminosity is:
By dividing the net force Fnet by −∆m = −ρ(r)∆r∆A, we find the equation of
motion of the shell:
If the acceleration is set to zero (when there is a balance), then the hydrostatic
equilibrium becomes:
dP GM (r)ρ(r)
=− (2.20)
dr r2
Therefore, the pressure at depth h 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:
P = gρh (2.21)
GM ρ M
Pc = ρ= (2.22)
R 4πR3 /3
16 CHAPTER 2. THE SUN A TYPICAL STAR
• hydrostatic equilibrium
dP GM ρ
=− 2 (2.24)
dr r
• mass continuity
dM
= 4πr2 ρ (2.25)
dr
• gradient of luminosity
dL
= 4πr2 ρ (2.26)
dr
• temperature gradient
dT 3κLρ
=− (2.27)
dr 16πacr2 T 3
In these equations r is the distance from the stellar center, P, ρ, T are the
pressure, density and temperature at radius r, M is the mass contained within
r, L the energy carried by radiation across r, 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.
E = mc2 (2.31)
If one assumes that the Sun consists of pure hydrogen which is converted into He,
then the total energy (E = 0.007mc2 ) would be 1.27 × 1045 J. The luminosity of the
Sun is L = 3.8 × 1026 J/s thus there would be energy supply for 1011 years.
The so called pp chain (Table 2.5) dominates in stars with relatively low cen-
tral temperatures (between 5 and 15×107 K, like the Sun) and the CN cycle7 is
dominant in stars with higher central temperatures.
The energy production rate, , for the pp cycle depends highly on the temper-
ature:
∼ ρT 5 (2.32)
61 eV=1.6 × 10−19 J
7 In the CN cycle C acts as a catalyst to convert H into He
18 CHAPTER 2. THE SUN A TYPICAL STAR
P dρ 1
< (2.36)
ρ dP γ
The theory which is usually used contains a free parameter, the so called mixing
length l:
dP
l = αHp = α P (2.37)
dr
where HP is the pressure scale height. It is supposed that α is of order unity. As we
will discuss later, information about the depth of the convection zone comes from
a detailed study of solar oscillations. Apart from energy transport one has also to
consider that in convection zones there is a uniform chemical composition. This
prevents any attempt of heavy chemical elements to settle in the Sun’s gravitational
field.
Table 2.6: Solar model: variation of temperature, luminosity and fusion rate throughout
the Sun
Radius fraction Radius Temperature % Luminosity Fusion rate
in R [109 ] m [106 ] K [J/kg s]
0 0.00 15.7 0 0.0175
0.09 0.06 13.8 33 0.010
0.12 0.08 12.8 55 .0068
0.14 0.10 11.3 79 .0033
0.19 0.13 10.1 91 .0016
0.22 0.15 9.0 97 0.0007
0.24 0.17 8.1 99 0.0003
0.29 0.20 7.1 100 0.00006
0.46 0.32 3.9 100 0
0.69 0.48 1.73 100 0
0.89 0.62 0.66 100 0
As it will be discussed later, the interior of the Sun can be investigated by the
propagation of waves. Solar models computed with mass loss, microscopic diffusion
of helium and heavy elements, and with updated physics have been evolved from
the pre-main sequence to present day (Morel et al., 1997 [225]); they are compared
to the observational constraints including lithium depletion and to the seismic
reference model derived by inversion. Microscopic diffusion significantly improves
the agreement with the observed solar frequencies and agree with the seismic
reference model within ±0.2% for the sound velocity and ±1% for the density,
but slightly worsens the neutrino problem. A review on the current state of solar
modeling was given by Christensen-Dalsgaard et al. (1996) [67] .
Several factors influence the image quality of solar telescopes. Sunlight can heat
up the telescope structure and the main optics causing the so called internal seeing.
Considering reflecting telescopes, in particular the main mirror absorbs up to 10%
of the collected light and its surface may heat up considerably leading to mirror
seeing. The most effective measure to prevent internal seeing is to remove the air
entirely, the telescope is evacuated. A window at the entrance and exit preserves
the vacuum. The main problem here is to have a window with high optical quality
which is thick enough to resist air pressure. Helium filling is an alternative to
evacuation. The viscosity of He and the dependence of the index of refraction from
temperature are lower than for air whereas temperature conductivity is higher. A
forced flow of He inside the telescope tube cancels inhomogeneities. The THEMIS
telescope (Télescope Héliographique pour l’étude du Magnetisme et des Instabilités
Solaires, see Arnaud et al. 1998 [12], Mein, 1997 [222]) has a He filled tube.
Other possibilities are to construct open telescopes, such as the DOT (Dutch open
telescope at La Palma, see Rutten et al., 2000 [266]). For THEMIS a joint between
the telescope tube and the dome which has an entrance window of 1 m prevents
air exchanges between outside and inside the dome.
To reduce atmospheric turbulence, systems which dynamically control the
wavefront deformations effected by the atmosphere are used (adaptive optics)8 .
Under ideal circumstances, the resolution of an optical system is limited by the
diffraction of light waves. This so-called diffraction limit is generally described by
the following angle (in radians) calculated using the light’s wavelength and optical
system’s pupil diameter d:
1.22λ
α= (2.38)
d
The turbulent atmosphere blurs images to resolution of 0.5 to 1 arcsec even at
the best sites. Adaptive optics (AO) provides a means of compensating for these
effects, leading to appreciably sharper images sometimes approaching the theo-
retical diffraction limit. One technique that has been developed for overcoming
atmospheric blurring is speckle interferometry, in which hundreds of very short ex-
posures (”specklegrams”) are processed after the observations to reconstruct the
unblurred image. However, because the specklegrams must be short exposures
and at the same time should have good signal-to-noise, speckle interferometry is
limited to imaging very bright objects. All AO systems work by determining the
shape of the distorted wavefront, and using an “adaptive” optical element – usually
a deformable mirror – to restore the uniform wavefront by applying an opposite
cancelling distortion.
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 and is
currently being upgraded.
In Fig. 2.4 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, the Swedish La Palma
Solar Telescope (SST), THEMIS in Tenerife, the NSO R.B. Dunn Telescope at
Sacramento Peak (DST, 30 inch entrance window, 20 inch solar image) and the
McMath-Pierce Facility at Kitt Peak..
GREGOR is a cooperation between three German institutes (Kiepenheuer
Institut Freiburg, Astrophysikal. Institut Potsdam and Astrophysikal. Institut
Göttingen). An existing telescope will be upgraded to a 1.5 m telescope and will
permit high resolution solar observations. The 1.5 m mirror is a lightweight mirror
consisting of Silicon Carbide (weight only 180 kg) 9 .
The SST 10 is the Swedish Vacuum solar telescope located at the Observatorio
Roque de los Muchachos in La Palma. It is a 1 m vacuum telescope. A single lens
is used thus permitting observations of high image quality through very narrow
filters. The Advanced Technology Solar Telescope (ATST) will be a collaboration
between 22 institutes11 .
9 [Link]
10 //[Link]/NatureNov2002/telescope [Link]
11 [Link]
22 CHAPTER 2. THE SUN A TYPICAL STAR
Flat
Coelostat
Concave
Exit window
Solar image
Ground level
Vertical spectrograph
ing
Grat
Figure 2.4: Optical path scheme of a vacuum telescope (e.g. Kitt Peak or VTT, Tener-
ife). Below the ground level, a vertical spectrograph is located. The solar image can be
observed at the top of the optical bank that is shown as a black box in the sketch.
2.3. OBSERVING THE SUN 23
Figure 2.5: Ground-based solar observatories that operate on a routine basis. After K.
Reardorn.
Besides these instruments for high resolution studies a number of solar moni-
toring instruments are distributed over the world (see Fig. 2.5).
A review about solar instrumentation was given by v.d. Lühe (2001) [202]. In
the next section we discuss satellite based instruments.
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
Earth12 which permits an uninterrupted view of the Sun. All previous space
solar observatories have orbited the Earth, from where their observations were
periodically interrupted as our planet ‘eclipsed’ the Sun.
The scientific objectives of SOHO are the interior of the Sun, the solar at-
mosphere and the solar wind.
The main instruments are listed in Table 2.7 together with the acronyms.
Further details about SOHO and the instruments can be found e.g. in the
review given by Fleck (2001) [98].
A drawing of the spacecraft is given in Fig. 2.6.
12 This is a stable point in the Sun-Earth system
24 CHAPTER 2. THE SUN A TYPICAL STAR
Figure 2.6: Drawing of the SOHO solar observatory (ESA & NASA)
2.3. OBSERVING THE SUN 25
RHESSI
The acronym RHESSI stands for Reuven Ramaty High Energy Solar
Spectroscopic Imager and it was launched on Feb. 5, 2002. The primary mis-
sion objective is to explore the basic physics of particle acceleration and energy
release in solar flares. The satellite is located at a circular orbit at 600 km alti-
tude. Observations can be made in the 3 keV to 17 MeV range by using cooled
hyperpure germanium crystals. The angular resolution depends on the energy
range of the observations: 2 arcseconds to 100 keV 7 arcseconds to 400 keV 36
arcseconds above 1 MeV. During solar flares large amount of energy is released,
gas being heated up to several million K radiating X-rays. High energy electrons
primarily emit hard X-rays, whereas high energy protons and ions emit primarily
gamma rays. RHESSI permits to observe both types of radiation. It seems that
the mission will be operational until 2008.
STEREO
This acronym stands for Solar Terrestrial Relations Observatory. The mission
will consist on two nearly identical space-based observatories - one ahead of Earth
in its orbit, the other trailing behind - to provide the first-ever stereoscopic mea-
surements to study the Sun and the nature of its coronal mass ejections, or CMEs.
Launch of the STEREO spacecraft is planned for no earlier than September 18,
2006. The main aim will be to study the 3-D structure and extension of coronal
mass ejections (CMEs).
26 CHAPTER 2. THE SUN A TYPICAL STAR
Workshop Held in Bangalore, India, 12-16 October, 1998, Kluwer, ASSL, 243
2.3. OBSERVING THE SUN 27
Figure 2.8: Zeeman effect: Without an external magnetic field the energy levels of an
electron in an atom a,b,c as well as d,e,f coincide. The presence of a magnetic field splits
the energy levels, however not all transitions are possible there are specific transition
rules.
When doing solar polarimetry one has to take into account for different effects such
as the Earth’s atmosphere where time fluctuations produce wavefront distortions,
instrumental polarization, spectrograph and detector polarization which all enter
as factors in the Mueller matrix.
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.
28 CHAPTER 2. THE SUN A TYPICAL STAR
Figure 2.9: Propagation of a wave throughout the outer solar atmosphere. On the
abscissa is the time, on the ordinate the frequency. Within 5 min the frequency drifts
from 80 MHz to 40 MHz indicating the propagation to the higher corona. Courtesy: H.
Aurass, Th. Mann, AIP.
trino νe like the decay of tritium or 3 H (halflife 12.33 years) into 3 He:
3
H →3 He + e− + νe (2.42)
σ ∼ 10−50 m2 (2.43)
When we compare this value to the cross sections in atomic and nuclear physics
which are about 10−20 ...10−30 m2 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 l, if the target particles have
a number density n m−3 , is given by:
1
l= (2.44)
nσ
For a solid target one has n ∼ 1029 and therefore l ∼ 1021 m for neutrinos. So neu-
trinos have an extremely large mean free path. How can we detect them? There
are many neutrinos coming from the Sun passing the Earth: about 1015 m−2 s−1 .
With the cross section and the number density given above, the number of detec-
tions N would be:
N ∼ 10−6 m−3 s−1 (2.45)
That means about one neutrino per month per cubic meter of the detector.
2.4. NEUTRINOS-TESTING THE SOLAR INTERIOR 31
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 37 Cl is in a tank containing 105 gallons of C2 Cl4 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 230 37 Cl atoms one has to expect one detection
every 5 × 105 s.
Theoretical models of the Sun predict the following count rates:
Bahcall, Pinsonneault (1992) [22] 8.0 ± 3.0 SNU,
is ∼ 1030 yr. The neutrino detector picked up a number of neutrinos from the
explosion of the supernova SN 1987A in the Large Magellanic Cloud, which is
a neighbor of our galaxy. In an updated version (Kamiokande II) 0.68 kilotons
of water were used and neutrinos from 8 B above 7.5 MeV can be detected. The
water experiment Kamiokande detects higher energy neutrinos (above 7 MeV)
by neutrino-electron scattering (ν + e → ν‘ + e ) and according to the standard
solar model the 8 B decay is the only important source of these higher-energy
neutrinos. The experiment clearly showed that the observed neutrinos come from
the sun because the electrons that are scattered by the incoming neutrinos recoil
predominantly on the direction of the sun-earth vector. The relativistic electrons
are observed by the Cerenkov radiation they produce in the water detector.
The results of the gallium experiments, GALLEX and SAGE gave an average
observed rate of 70.5±7 SNU. This is in agreement with the standard model by the
theoretical rate of 73 SNU that is calculated from the basic pp and pep neutrinos.
The 8 B neutrinos which are observed above 7.5 MeV in the Kamiokande exper-
iment, must also contribute to the gallium event rate. This contributes another
7 SNU, unless something happens to the lower energy neutrinos after they are
created in the Sun. Thus the Ga experiments are in accordance with predictions
if we exclude everything but the pp neutrinos. This is sometimes called the third
neutrino problem.
The calculated pp neutrino flux is approximately independent of solar models;
it is closely related to the total luminosity of the sun.
Summarizing the the neutrino problem we can state:
• smaller than predicted absolute event rates in the chlorine and Kamiokande
experiments.
• very low rate in the Ga experiment which implies the absence of 7 Be neutri-
nos although 8 B neutrinos are present.
Solar neutrino experiments are currently being carried out in Japan (Super
Kamiokande,Takata (1993) [309], Totsuka (1996) [313]), Canada (SNO, Sudbury,
using 1 kiloton of heavy water; Hargrove and Paterson, 1991 [127]) and in Italy
(BOREXINO, ICARUS, GNO (Gallium Neutrino Observatorium), each sensitive
to a different energy all working in Gran Sasso, Arpesella et al. (1992) [13]), in
Russia (SAGE, Caucasus) and in the United States (Homestake). The SAGE,
chlorine and GNO work radiochemical, the others electronic (recoil electrons pro-
duced by the neutrino interactions using Cherenkov effect).
The Sudbury Neutrino Observatory (SNO) (see also Krastev (2002) [173]) is
located at 6800 feet under ground in a mine in Sudbury, Ontario. The neutri-
nos react with heavy water producing flashes of light (Cerenkov radiation). The
principle is as follows:
νe + d → p + p + e− (2.50)
As the neutrino approaches the deuterium nucleus d a heavy charged particle of
the weak force (called the W boson) is exchanged. This changes the neutron in
2.4. NEUTRINOS-TESTING THE SOLAR INTERIOR 33
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) [331] and Roxburgh (1998) [263].
Is there a correlation between neutrino fluxes and solar activity? On the basis
of an analysis of the 37 Ar production rate at the Homestake station for the period
1970-1990, Basu (1992) [28] found that the solar neutrino flux varies with time
in proportion to the solar wind flux. However, Walther (1999) [330] 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 the
reported highly significant results for segments of the data are due to a statistical
fallacy.
How to explain this discrepancy between observations and theory? One explana-
tion comes from particle theory itself. There are three conserved quantities, called
electron, muon and tauon lepton numbers and correspondingly three types of neu-
trinos. There are however indications that some modifications to this standard
model are required. These involve that the neutrinos have small masses and that
the neutrinos can transform from one type to another. The Mikheyev-Smirnov-
Wolfenstein effect (MSW) explains such neutrino oscillations and by the above
mentioned experiments we can only detect electron neutrinos. Another explana-
tion of the discrepancy is that the flux is variable during the solar cycle. This
might be explained if the neutrinos possess a magnetic moment and if interaction
with the solar magnetic field is possible.
Let us give a very simplified explanation of neutrino oscillations. An indispens-
able, but counterintuitive, concept in quantum mechanics is that of superposition.
Suppose a certain particle has a property that can have several different values;
the classic example is that of Schrödinger’s cat, but let us consider a more prac-
tical one: ordinary playing cards have the property ‘suit’, with the four possible
values ‘spades’, ‘hearts’, ‘diamonds’, and ‘clubs’. In ordinary non-quantum life,
each individual card has a well-defined suit. However a quantum card may be in
a mixed state, a superposition of, say 30% spades, 60% hearts, and 10% clubs.
When you check which suit that card belongs to, you have a 30% chance of finding
that it’s a spade, 60% chance of finding it’s a heart, and so on. Note that this
is not just a matter of your ignorance of the card’s “true” suit – the point is, it
doesn’t have a single well-defined suit until you check it.
2.4. NEUTRINOS-TESTING THE SOLAR INTERIOR 35
In particle physics, the equivalent of the suits are the three families, discussed
above. A neutrino may belong to any one of the three families, making it νe , or
a νµ , or a ντ . Or, it may be a superposition of the three family flavors, mixed in
some proportions. Now, the standard model assumes that the neutrinos emitted
from the sun are in a pure νe state, without mixing. This can be understood
with the quantum mechanical concept of eigenstates. This is well known for the
K meson. An eigenstate is a state that is recognized as pure, non-mixed, without
superposition, in a certain context. In quantum mechanics different interactions
recognize and interact each with a different set of eigenstates for the particles. Try
to apply this to a card game. In different games a heart would become a spade
etc. For most particles and interactions the different eigenstates are identical.
This is not the case for the weak interaction. The weak eigenstates of quarks
are different from their strong/electromagnetic eigenstates. The K0 mesons are
produced in strong interactions of quarks, but decay through weak interactions of
their constituent quarks. Thus, the production eigenstates are different from the
travel/decay eigenstates of the K0. As far as the weak interaction is concerned,
leptons are expected to behave in the same manner as quarks. If neutrinos do have
a tiny mass, and different neutrinos have different masses, they will behave in the
same way as K0 mesons. They will be produced in a weak-interaction eigenstate,
but travel in a mass eigenstate. The mass eigenstate may be different from the
weak eigenstate. The weak-interaction eigenstates are the three neutrino flavors
discussed earlier: νe , νµ , ντ . When they arrive and interact with our detectors,
they do not arrive as the original weak eigenstate in which they were produced,
but as a mixture of two or more flavors. This is a potential solution to the solar
neutrino problem, since the experiments measure an apparent disappearance of
electron-neutrinos, without measuring the other flavors. If the neutrinos oscillate
from the 100% νe that they are produced as in the sun, to a mixture with around
40% νe electron-neutrino and 60% some other neutrinos, we get an agreement with
experimental data.
Neutrino oscillations and the solar neutrino problem are discussed by Haxton
(2001) [129] .
The search for neutrino decays during the 1999 solar eclipse is discussed in
Cecchini et al. (2000) [60] involving the emitted visible photons, while neutrinos
travel from the Moon to the Earth.
Today the most likely solution of the neutrino problem are the neutrino oscillations
suggesting that our solar model is quite correct. Bahcall and Davis (2000) [21] gave
a recent review about the solar neutrino problem and suggest further experiments.
Figure 2.10: SOHO-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 granu-
lation and supergranulation are seen.
High l
c1
Low l
c
2
Figure 2.11: 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 (c2 > c1 ).
38 CHAPTER 2. THE SUN A TYPICAL STAR
one takes a sequence of images of the oscillation pattern at fixed time intervals.
The shorter these time intervals between the images, the easier it is to identify
the oscillations. A big problem in such a project is the enormous amount of data.
Each station in the network produces more than 200 megabytes of data every day.
Details about the instrument used (a Fourier Tachometer) can be found in Beckers
et al. (1978) [30].
The BiSON (Birmingham Solar Oscillations Network) project also has six ob-
servatories, most of which are automated. As it is explained above, the GONG
observatories measure the motions on the solar surface caused by the oscillations.
The BiSON observatories do so as well, but unlike the GONG network they mea-
sure an average velocity over the solar surface (the Sun is observed as a point
source, if it were a star). The measurements therefore are sensitive only to oscilla-
tion patterns with very big wavelengths: all smaller-scale patterns are suppressed
by being averaged. The two techniques for GONG and BiSON are therefore com-
plementary.
There exist also surface waves which are called f-modes. The p-modes have fre-
quencies between 1 hour and two minutes and include the five minutes oscillations
discussed above. The g-modes have much longer periods than the p-modes. It can
be shown that they are trapped in the solar interior beneath the convection zone.
The energy generated in the sun is first transported by radiation and then at a
depth of about 200 000 km by convection. In this convection zone the amplitudes
of the g-modes are damped exponentially and thus it is extremely difficult to ob-
serve them at the solar surface. Amplitudes would be expected of a few cm/s to
mm/s but the frequency of these modes would contain valuable information about
the solar core (Turck-Chièze et al., 2004 [318]. Duvall, 2004 [83], suggested a new
method (time distance helioseismology) to detect g modes.
How can we describe the solar oscillations? First we must make some sim-
plifications. We assume that the sun is strictly spherical. This will provide a
spectrum of oscillation frequencies which will be modified by a) rotation and b)
magnetic fields. A second approximation is that the oscillations are adiabatic.
This approximation is valid since the oscillation period is in general much smaller
than the relevant thermal timescale. A third approximation is that we neglect a
change of the gravitational field of the Sun during the oscillation. This is not true
for radial oscillations: in radial oscillations all matter at any solar radius moves
inward or outward in phase. However if we consider nonspherical modes at short
wavelengths in the horizontal direction this is again a good approximation.
Any oscillation can be described by introducing three quantum numbers n, l, m.
The meaning of these numbers is as follows:
2.5. HELIOSEISMOLOGY-SOLAR OSCILLATIONS 39
• n denotes the number of points in the radial direction at which the amplitude
of the oscillation vanishes.
• l, m determine the angular behavior of the oscillation over the surface of the
Sun. In addition we have the relation −l ≤ m ≤ +l.
If Plm denotes the associated Legendre function which can be given in an analytical
form, the inward or outward motion of points on the surface is related to the value
of the real part of the function
where Θ, φ are spherical polar coordinates. If l, m are low, there is a relatively small
number of patches on the solar surface (which oscillate with different directions
of radial velocity). If l, m are large, there is a very large number of such patches.
We speak of a high degree model if l is large and conversely if l 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 l 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.
The l −ν diagram (Fig. 2.13) 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
40 CHAPTER 2. THE SUN A TYPICAL STAR
Figure 2.13: l-ν diagram from MDI high-cadence full disk data shows mode frequencies
up to 10 mHz and l=1000.
to sound waves we are used to hearing. Most of the power is concentrated in a band
near 3 mHz, that’s one oscillation every 5 minutes 16 . Higher frequencies aren’t
trapped inside the Sun, so they don’t resonate. Modes with lower ν disappear in
the background noise. The spatial scale of the modes is indicated by the angular
degree l telling how many node lines there are in the pattern at the surface of
the Sun. The l=0 modes are ‘breathing’ modes where the whole surface of the
Sun moves in and out at the same time. Higher order modes divide the surface
into a pattern like a checker board, where adjacent squares move in different
directions at any given time. A mode of a particular degree, l, at the surface can
be associated with resonances having any number of nodes in the radial direction
inside the Sun. The number of radial nodes is called the order. The curved lines in
the figure are associated with different radial orders. For a given order (line) the
16 Sound waves we can hear vibrate from tens to thousands of times per second
2.5. HELIOSEISMOLOGY-SOLAR OSCILLATIONS 41
frequency decreases with increasing spatial degree. For a give degree, the frequency
increases with order. In the Fig. 2.13, the lower left corner is most closely related
to what is happening in the core of the Sun. Moving up in frequency or degree
tells more about what is happening near the surface. Because sound waves of a
particular degree can travel in different directions the lines appear relatively broad.
If the material through which any of these modes is travelling is moving, then the
measured frequency of the mode is affected. The rotation of the Sun causes the
biggest frequency shift and makes the lines shown in the figure broad (frequency
shifting). Other motions within the Sun along the path taken by the waves cause
different types of frequency changes. Analysis of these frequency changes reveals
the internal motions of the Sun.
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 v is the fluid velocity, Γ is an effective ratio of specific
heats (ρdP/P dρ) 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
d/dt = ∂/∂t + vgrad. In an equilibrium situation:
Now we consider small disturbances about this equilibrium in which the perturbed
quantities are functions of all the spatial coordinates and the time. In the equi-
librium there is no dependence on spherical polar coordinates. For any variable f
we can write:
If the star is spherical the oscillation frequency does not depend on m. For the
Sun, the departure from sphericity is small and the real oscillation modes have a
42 CHAPTER 2. THE SUN A TYPICAL STAR
behavior close to that shown above but with different m modes having different
frequencies. The oscillation frequency ω depends on n and l. The three numbers
n, l, m are related to the numbers of times f1 vanishes in the radial-, Θ- and
φ-directions and m ≤ l.
The functions f1 (r) 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
v = dξ/dt (2.61)
• Lamb frequency Sl
Sl = cs [l(l + 1)]1/2 /r (2.66)
• Brunt-Väissälä frequency N
1 dP 1 dρ
2
N =g − (2.67)
ΓP dr ρ dr
where g = GM/r2 .
Sl is always real but ωc and N can be imaginary. It can be shown that con-
vection occurs when N 2 is negative. We can write our differential equation for ψ
as:
d2 ψ
+ Kr2 ψ = 0 (2.68)
dr2
2.5. HELIOSEISMOLOGY-SOLAR OSCILLATIONS 43
Figure 2.14: 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.
• The transition layer between the radiative and convection zone which is
called the tachocline is mostly located in the radiative zone and thin at the
equator but maybe wider at high latitudes.
• There is a sharp radial gradient of the angular velocity beneath the con-
vection zone and the narrow peak of the sound speed at 0.67 R is due to
rotationally turbulent mixing in the tachocline.
More details about these results can be found in Kosovichev et al. (1998) [171].
Helioseismology can be used also to give arguments in the question of solar
neutrinos. Turck - Chièze et al. (2001) [317] used sound-speed and density pro-
files inferred from SOHO/GOLF and SOHO/ MDI data including these modes,
together with recent improvements to stellar model computations, to build a spher-
ically symmetric seismically adjusted model in agreement with the observations.
Their model is in hydrostatic and thermal balance and produces the present ob-
served luminosity. Some fundamental ingredients were adjusted, well within the
commonly estimated errors, such as the p-p reaction rate (±1%) and the heavy-
element abundance (±3.5%); the sensitivity of the density profile to the nuclear
reaction rates was examined. The corresponding emitted neutrino fluxes demon-
strate that it is unlikely that the deficit of the neutrino fluxes measured on Earth
can be explained by a spherically symmetric classical model without neutrino fla-
vor transitions.
2.5. HELIOSEISMOLOGY-SOLAR OSCILLATIONS 45
New insight into the internal structure of the Sun can be obtained by using
time-distance helioseismology. Let us explain this technique by considering seis-
mology on earth. Here, the arrival time of the initial onset of a disturbance is
measured. If we know the variation of seismic velocity with depth within the
earth, then we can calculate the travel time of rays between an earthquake and
a receiver using geometrical approximations. So in principle, we can locate any
earthquake in both time and space by recording the arrival times of waves at
stations worldwide.
In time-distance helioseismology, the travel time of acoustic waves is measured
between various points on the solar surface. To some approximation the waves
can be considered to follow ray paths; these depend on a mean solar model. The
curvature of the ray paths is caused by increasing sound speed with depth below
the surface (see Fig. 2.11). The travel time is affected by various inhomogeneities
along the ray path, including flows, temperature inhomogeneities and magnetic
fields. The technique consists of a measurement of a large number of times between
different locations. Then an inversion method is used to construct 3-D maps of
the subsurface inhomogeneities. A review article on that technique was given by
Duvall et al. (1997) [84].
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 ρ, T
with r. 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 T, ρ... 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) [119].
Solar like oscillations found on other stars are discussed recently e.g. in Bedding
and Kjeldsen, 2006 [31] and Kjeldsen et al., 2005 [165] where 37 oscillation modes
on α Cen B were found with l=0-3.
Chapter 3
3.1 Introduction
The different layers of the Sun and its atmosphere can be defined 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 2 × 105 km
into the interior.
1
1. Photosphere: starts at the surface and extends up to 500 km.
2. Chromosphere: above the photosphere; extends to about 2 Mm.
3. Transition Region: strong increase of temperature up to 106 K over a
very small spatial range (some 104 km).
4. Corona: starts above 2 Mm, high temperature > 106 K.
47
48 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
-10
-15
Figure 3.1: Variation of electron temperature and electron density in the solar at-
mosphere
∂Iν (τν , µ)
µ = Bν (T ) − Iν (τν , µ) (3.3)
∂τν
2 Very often the solar surface is defined as the layer where τ500 nm = 1
3.2. PHENOMENA IN THE SOLAR PHOTOSPHERE 49
2hν 3 hν/kT −1
Bν (T ) = e −1 (3.4)
c2
An elementary solution yields for the intensity of radiation emerging in direction
µ: ∞
τν dτ
Bν (T )e− µ
ν
Iν (µ) = (3.5)
0 µ
The Planck function must increase with depth, since the temperature increases
with depth (see Fig. 3.1). Eddington made the following Ansatz assuming a linear
increase of the function Bν with depth:
Bν = C + Dτν (3.6)
Iν = C + Dµ (3.7)
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.
3.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 corre-
sponds roughly to 1 arcsec (as seen from the Earth). In the bright granules matter
is streaming upwards, in the darker intergranular lanes streaming downwards. Up
50 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
to now the best granulation images have been taken from the ground since no large
solar telescopes have been launched. In 2006 SOLAR B will be launched. This
will be the first large optical telescope flown in space. Its aperture is 50 cm and
angular resolution achieved will be 0.25 arcsec.
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 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 review about solar granulation was given by Muller (1999) [227] where further
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). Finally, 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. The granule lifetime
can be determined by their visual identification on successive images or by cross
correlating these images. There is a large discrepancy of the results: granular
lifetimes range from 6 to 16 min.
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-
3.2. PHENOMENA IN THE SOLAR PHOTOSPHERE 51
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) [262] 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 for 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. 3.2).
Concerning the structural properties of granules, we have to mention that their
number N increases monotonically with decreasing size. Granules of size 1.4 arcsec
are the main contributors to the total granule area. When the area A is plotted
versus their perimeter in a log-log scale, the dispersion of points (each of them
marks a granule) is small and their shape can be characterized by the relation:
P ∼ AD/2 (3.10)
where D is the fractal dimension. It seems that there are two ranges with different
fractal dimensions:
52 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
Theoretical Approaches
The simplest model of convection is the classical Rayleigh problem: suppose a
fluid (either gaseous or liquid), confined between two horizontal plates separated
by a distance h and maintained at temperature T1 (upper) and T2 (lower) with
T2 > T1 . 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 (T2 − T1 )/h is sufficiently large or, for a given gradient, when the
coefficients of the kinematic viscosity ν 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:
gαβh4
R > Rcrit R= (3.11)
κν
where β is the temperature gradient. For Rcrit 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):
4h4 Ω2
C= (3.12)
ν2
here, Ω is the vertical component of the angular velocity vector. For details see
e.g. Chandrasekhar (1961) [62].
For the solar convection zone R is extremely high, R ∼ 1010...11 .
3.2. PHENOMENA IN THE SOLAR PHOTOSPHERE 53
Important information about the origin of the solar granulation can be in-
ferred from power spectra. From 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 k −5/3 down to the scale of molecular diffusion. The
temperature power spectrum however decreases as k −5/3 only to a scale kc . At
smaller scales the spectrum decreases as k −17/3 (Espagnet et al., 1995 [90]). Thus
kc 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 de-
scribed in detail in Nordlund, 1982 [234].
Granulation-Mesogranulation
Idealized numerical experiments on turbulent convection were made by Catta-
neo et al. (2001) [58]. 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 macro-
cells (corresponding to mesogranules) are several times larger. The motion acts
as a small-scale turbulent dynamo, generating a disordered magnetic field that is
concentrated at macrocellular corners and, to a lesser extent, in the lanes that
54 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
join them. These results imply that mesogranules owe their origin to collective
interactions between the granules.
Other authors like e.g. Hoekzema et al., 1998 [135], who analyzed G band
images found that photospheric 5 min oscillations are global and rather insensitive
to local fine structure.
Using a 30-min time series of CCD spectrographs, Khomenko, Kostik and
Shchukina, 2001 [160], found different amplitudes, phases and periods of the 5-min
oscillations above granules and intergranular lanes. The most energetic intensity
oscillations occurred above intergranular lanes, the most energetic velocity oscil-
lations above granules and lanes with maximum contrast that are cospatial with
regions with maximum convective velocities.
3.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 by four men: 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.
After the initial interest and the publication of Scheiner’s major work the
interest in sunspots vanished. 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 known 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 100 as longi-
tude zero and a rotation completed by this point is known as Carrington rotation
(CR)3 . The sideral Carrington rotation is 25.38 days, the synodic value varies a
3 For example on March 14 2006 Carrington Rotation 2041 started at 14.43 UT and ended on
little during the year because of the eccentricity of the Earth’s orbit (its mean
value is 27.2753 days).
Carrington was also the first to see a white light flare on the Sun in the morning
of Sep. 1, 1859, during sketching sunspot projections with a friend. Suddenly
two crescent-shaped patches broke out, brightened, moved a distance twice their
length, then faded away as two dots within five minutes. Carrington 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 (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 RZ as a measure for solar activity. Sunspot often appear as groups. If g
denotes the number of sunspot groups and f the number of individual spots, then
RZ = k(10g + f ) (3.13)
Figure 3.4: Large sunspot showing the dark central umbra and the filamentary penum-
bra. Outside the penumbra the granulation pattern is clearly seen. Courtesy: M.
Sobotka, A.H., SST, La Palma, 2003
the energy levels. The displacement of the lines due to the Zeeman effect is given
by:
∆λ = 4.7 × 10−8 g ∗ λ2 B (3.14)
The wavelength λ is given in nm, the Landéfactor g ∗ depends on the spin and
orbital momentum of the levels and B denotes the magnetic induction given in
Tesla.
bridges, dark nuclei in the umbra etc. One problem in the study of sunspots and
their fine structure is observational stray light.
An important photometric parameter of umbral cores is the minimum intensity
(intensity of the darkest point) Imin which is usually in the range of 0.05-0.3 of the
mean photospheric intensity at λ ∼ 540 nm. There seems to be a relation between
the size of the umbrae and the temperature. Umbrae with a diameter DU < 7
have higher temperatures than the large ones. Moreover, regions with higher
magnetic field strength are darker and cooler than those with lower strength. The
darkest regions in umbral cores are dark nuclei. These are the areas with the
strongest magnetic fields and the orientation of the field is perpendicular to the
surface of the Sun. They are not necessarily centered in the umbral cores, some of
them are observed close to the edge of the penumbra. They cover 10-20% of the
total umbral core area and their size is about 1.5 arcsec4 . The penumbra shows
elongated structures which is a consequence of the strongly inclined magnetic
field. Bright penumbral filaments consist of penumbral grains. They seem to have
cometary like shapes with “heads” pointing towards the umbra and have a mean
width of only 0.36 and a length of 0.5...2 . The observed brightness approaches
the photospheric one and the lifetimes are between 40 minutes and 4 hours. They
are separated by narrow dark fibrils. The magnetic field seems to be stronger and
more horizontal in dark fibrils and weaker and more vertical in penumbral grains.
It is also interesting to note that nearly all penumbral fine structures are in
motion. The penumbral grains move towards the umbra with an average speed of
0.3-0.5 km/s. On the other hand, dark cloud like features which arise from the
dark fibrils move rapidly outwards (up to 3.5 km/s) towards the outer penumbral
border.
The last fine structure which is important to study are the light bridges . They
cross the umbra or penetrate deeply into it and can be observed for several days
although they change their shape substantially on the scale of hours. They can be
classified into faint (located inside umbral cores) and strong (separating umbral
cores). Strong light bridges separate umbral cores of equal magnetic polarities
and a subclass of them opposite polarities. The analysis of 2-D power spectra of
intensity fluctuations inside strong light bridges showed that the “granules” that
can be seen there are smaller (1.2 arcsec, normal granulation: 1.5 arcsec) and the
slopes of power spectra indicated the presence of a Kolmogorov turbulent cascade.
The magnetic field strength in strong light bridges is substantially lower than in
adjacent umbra.
A recent review about the fine structure of sunspots was given by Sobotka
(1999) [290] where other references can be found. A review on empirical modelling
and thermal structure of sunspots was given by Solanki (1997) [292].
22 year magnetic cycle. Spots appear as a magnetic flux tube rises (see magnetic
buoyancy) and intersects with the photosphere. The magnitude of the magnetic
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 (10−4 T). 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:
1 dψ ∂ψ
B= − , 0, (3.16)
ω dz ∂ω
Thus curlB = 0. Since divB = 0,
∂2ψ 1 ∂ψ ∂ 2 ψ
− + =0 (3.17)
∂ω ω ∂ω ∂z 2
The neighboring photosphere, in which the flux tube is embedded has a known
pressure variation with height Pe (z). The boundary of the flux tube is at ω =
ω0 (z), where
B 2 /2µ0 = Pe (z) (3.18)
We see that as z → ∞ the field becomes nearly horizontal and Bω ∼ F/2πω02 and
as z → −∞, the field becomes vertical and Bω ∼ F/πω02 .
There is one problem with this monolithic model: the difference in the energy
radiated by the spot and by an equivalent area of the normal photosphere is only
about a factor of 4. This is less than would be expected if convection in the spot
were completely suppressed. Therefore, it is believed that some form of convective
energy transport must occur and the field must be more complex e.g. coherent
flux tubes or a tight cluster. Reviews about these topics were given by Bogdan
(2000) [40] and Hurlburt (1999) [144].
Using the 1 m Swedish Solar Telescope a high resolution study of the inclina-
tion of magnetic fields within sunspots was performed by Langhans et al., 2005
[187]. Within sunspots, dark penumbral cores, and their extensions into the outer
penumbra, are prominent features associated with the more horizontal component
of the magnetic field from about 400 in the inner penumbra to nearly horizontal
in the middle penumbra. Bright flux component is associated with a more vertical
field component.
bra of the principal spot, and ‘c’ describes the distribution of spots in the interior
of the group. This classification scheme substituted the older scheme that was
introduced by Waldmeier (1938).
2. p-values:
x - no penumbra (group class is A or B)
r - rudimentary penumbra partially surrounds the largest spot. This penum-
bra is incomplete, granular rather than filamentary, brighter than mature
penumbra, and extends < 3 arcsec from the spot umbra. Rudimentary
penumbra may be either in a stage of formation or dissolution.
s - small, symmetric (like Zurich class J). Largest spot has mature, dark,
filamentary penumbra of circular or elliptical shape with little irregularity
to the border. The north-south diameter across the penumbra is ≤ 2.5
degrees.
a - small, asymmetric. Penumbra of the largest spot is irregular in outline
and the multiple umbra within it are separated. The north-south diameter
across the penumbra is ≤ than 2.5 degrees.
h - large, symmetric (like Zurich class H). Same structure as type ‘s’, but
north-south diameter of penumbra is more than 2.5 degrees. Area, therefore,
must be larger or equal than 250 millionths solar hemisphere.
k - large, asymmetric. Same structure as type ‘a’, but north-south diameter
of penumbra is more than 2.5 degrees. Area, therefore, must be larger or
equal than 250 millionths solar hemisphere.
3. c-values:
x - undefined for unipolar groups (class A and H)
o - open. Few, if any, spots between leader and follower. Interior spots of
very small size. Class E and F groups of ‘open’ category are equivalent to
Zurich class G.
i - intermediate. Numerous spots lie between the leading and following por-
tions of the group, but none of them possesses mature penumbra.
c - compact. The area between the leading and the following ends of the
3.2. PHENOMENA IN THE SOLAR PHOTOSPHERE 61
spot group is populated with many strong spots; at least one interior spot
shows a mature penumbra. The extreme case of compact distribution has
the entire spot group enveloped in one continuous penumbral area.
There exists also the Mount Wilson classification scheme:
α: Denotes a unipolar sunspot group.
β: A sunspot group having both positive and negative magnetic polarities, with a
simple and distinct division between the polarities.
β − γ: A sunspot group that is bipolar but in which no continuous line can be
drawn separating spots of opposite polarities.
δ: A complex magnetic configuration of a solar sunspot group consisting of oppo-
site polarity umbrae within the same penumbra.
γ: A complex active region in which the positive and negative polarities are so
irregularly distributed as to prevent classification as a bipolar group.
300
SUNSPOT NUMBER
200
100
0
1750 1760 1770 1780 1790 1800 1810 1820 1830 1840 1850
DATE
300
SUNSPOT NUMBER
200
100
0
1850 1860 1870 1880 1890 1900 1910 1920 1930 1940 1950
DATE
300
SUNSPOT NUMBER
200
100
0
1950 1960 1970 1980 1990 2000 2010 2020 2030 2040 2050
DATE
Figure 3.5: Relative Sunspot number. Two cycles can bee seen, the normal cycle with
about 11 years and the so called Gleissberg cycle with a period of about 80 years.
• 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. 3.6).
In Table 3.1 some parameters for the energetics of large sunspots are given, i.e.
spots with a diameter ≥ 3.5 × 104 km. Penumbral waves are horizontal outwards
waves (in Hα ) with velocities between 10 and 20 km/s.
30N
90N
90S
30S
EQ
0.0
0.1
0.2
0.3
0.4
0.5
[Link]
1880
1890
13
1900
1900
14
1910
1910
15
1920
1920
16
1930
1930
DATE
DATE
17
1940
1940
18
1950
1950
1960
1960
19
1970
1970
20
> 0.0%
1980
1980
21
NASA/MSFC/HATHAWAY 11/2001
> 0.1%
1990
1990
22
> 1.0%
2000
2000
Figure 3.6: Butterflydiagram illustrating the equatorward motion of spots during the
activity cycle.
sunspots they tend to overlap and can be identified further from the limb. They
appear in increased numbers in a region prior to the emergence of sunspots and
remain for a rotation or more after the spots have decayed. As it will be shown later
they are important for the energy balance between sunspots and the photosphere.
Faculae can be observed on the whole disk using filtergrams . In that case they
are often called plage and attributed to the chromosphere. Photospheric faculae
are manifestations of concentrated azimuthal magnetic fields. One possibility to
study sunspots and faculae at photospheric levels is to use the Ca II K line 0.05
nm off the center with a 0.015 nm passband.
Polar faculae appear as pointlike, bright photospheric spots near the solar limb
at latitudes of 55 degrees or more (average of 65 degrees). Polar faculae tend to
occur at lower latitudes (as low as 45 degrees) during the years in which there
are only few observable. They can be distinguished from main zone faculae by
64 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
their essentially pointlike and solitary appearance, in contrast to the more area-
and grouplike appearance of the main zone faculae (55 degrees or lower). Their
lifetime is shorter (minutes to hours) than that for ordinary faculae. The brightest
can last for a couple of days, and can be traced farther from the solar limb too.
In connection with the activity cycle it is interesting to note that polar faculae
are most numerous at times of minimum solar activity, which in turn might be an
additional hint for their relation with the upcoming new solar cycle.
• Temperature minimum: near 500 km; here the UV continuum near 160 nm,
the far IR continuum and the minima in the wings of Ca II and Mg II lines
are formed,
• κν absorption coefficient
• ην emission coefficient
66 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
Sν = ην /κν (3.20)
Sν = Bν (T ) (3.21)
which is called Kirchhoff’s law, Bν (T ) being the Planck function. We can progress
to solve the transport equation:
∞
Iν (τν = 0, µ) = Sν (τν ‘)e−τν ‘/µ dτν ‘/µ (3.22)
0
In this equation µ = cos θ, θ being the angle between the normal to the disk center
and the point where observations are done.
From a Taylor series expansion of Sν about a not specified τν∗ one gets
Iν ∼ Sν (τν ) = µ (3.23)
where τν∗ was specified to µ. That means, one observes under the angle θ to z
approximately the source function at optical depth τν = µ.
Let us consider two energy levels in an atom which have the quantum numbers
l (lower level) and u (upper level). The number of atoms per cm3 in the lower level
is Nl and in the upper level Nu . Of course a transition from l to u corresponds
to an absorption process, where a photon of energy hνl,u = χu − χl is absorbed.
Thus the number of transitions per cm3 is given by:
B(l, u) is the transition probability for the transition l → u. On the other hand let
us consider the number of spontaneous transitions from u → l which is independent
on the intensity J:
n(u → l) = Nu A(u, l) (3.25)
A(u, l) is the transition probability for spontaneous transitions. Generally, we do
not know the average intensity Jν(l,u) . However, in thermodynamic equilibrium it
is equal to the Planck function. In thermodynamic equilibrium there is a direct
balancing between the number of transitions u → l and l → u and the ratio of the
occupation numbers is governed by the Boltzmann formula:
Nu gu −(χu −χl )/kT
= e (3.26)
Nl gl
3.3. THE CHROMOSPHERE 67
and
n(l → u) = n(u → l) (3.27)
3
2hν 1
Nl B(l, u) = Nu A(u, l) (3.28)
c2 ehν/kT − 1
where we have put the Planck function. Let us also substitute the Boltzmann
formula:
2hν 3 1 gu −(χu −χl )/kT A(u, l)
= e (3.29)
2
c e hν/kT −1 gl B(l, u)
gu −hνu,l /kT A(u, l)
= e (3.30)
gl B(l, u)
where gu , gl are the statistical weights of the states u, l. This was first found by
Einstein. Besides absorption and spontaneous emission also the induced emission,
transitions from u → l depending on the intensity J, has to be considered. The
number of induced emissions is written as:
n (u → l) = Nu B(u, l)Jν(u,l) (3.31)
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:
Nl Jν(u,l) B(l, u) − Nu Jν(u,l) B(u, l) = Nu A(u, l) (3.32)
and using Jν(u,l) = Bν and the Boltzmann formula:
2hν 3 1 gl hνu,l /kT
B(l, u) e − B(u, l) = A(u, l) (3.33)
c2 ehν/kT − 1 gu
B(u, l)gu = B(l, u)gl (3.34)
3 3
gl 2hνu,l 2hνu,l
A(u, l) = B(l, u) = B(u, l) (3.35)
gu c2 c2
These relations are called Einstein transition probabilities.
B(u, l), B(l, u), A(u, l) are atomic constants. Though these relations were de-
rived from thermodynamic equilibrium, they must always hold. Therefore, they
can be used to get information for excitation conditions and the source function
in case we do not have thermodynamic equilibrium.
By these calculations one can understand the typical profile of the Ca II H and
K lines (see Fig. 3.7). There are two intensity minima on the blue and red side
of the line center (called K1v , K1r ), towards the line center two maxima (called
K2v , K2r ) and then at the line center there is a minimum (K3 ). This indicates
that the temperature increases in the chromosphere. While the source function
decouples from the Planck function it reaches a minimum K1 , exhibits a small
maximum K2 and finally drops towards the line center. The profile of the well
known Hα 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) [168].
68 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
The magnetic elements on the cell boundary are associated with intense mag-
netic fields and are heated by waves with very long periods, ranging from six
to twelve minutes; the bright points are associated with magnetic elements of low
field strength and are heated by (long-period) waves with periods near the acoustic
cutoff period of three minutes; and the quiet cell interior, which is free of mag-
netic field, may be heated by short-period acoustic waves, with periods below one
minute. This paper reviews mainly the heating of the bright points and concludes
that the large-amplitude, long-period waves heating the bright points dissipate
enough energy to account for their chromospheric temperature structure.
Skartlien et al. (2000) [287] studied the excitation of acoustic waves using
three dimensional numerical simulations of the nonmagnetic solar atmosphere and
3.3. THE CHROMOSPHERE 69
the upper convection zone. They found that transient acoustic waves in the at-
mosphere are excited at the top of the convective zone (the cooling layer) and
immediately above in the convective overshoot zone, by small granules that un-
dergo a rapid collapse, in the sense that upflow reverses to downflow, on a timescale
shorter than the atmospheric acoustic cutoff period (3 minutes). The location of
these collapsing granules is above downflows at the boundaries of mesogranules
where the upward enthalpy flux is smaller than average. An extended downdraft
between larger cells is formed at the site of the collapse. The waves produced
are long wavelength, gravity modified acoustic waves with periods close to the 3
minute cutoff period of the solar atmosphere. The oscillation is initially horizon-
tally localized with a size of about 1 Mm. The wave amplitude decays in time as
energy is transported horizontally and vertically away from the site of the event.
They also made a prediction of how to observe these “acoustic events”: a darken-
ing of intergranular lanes, which could be explained by this purely hydrodynamical
process. Furthermore, the observed “internetwork bright grains” in the Ca II H
and K line cores and associated shock waves in the chromosphere may also be
linked to such wave transients.
The coronal heating problem can be also studied by an energy release that is
associated with chromospheric magnetic reconnection. A one-dimensional circu-
larly symmetric supergranulation reconnection model was investigated by Roald et
al. (2000) [261] with typical quiet-Sun values. In this model, the assumed source
rate of elements determines heating, because all emerged elements eventually an-
nihilate.
As an example for observational evidence we cite the paper of Ryutova and
Tarbell(2000) [267]. They analyzed spectra of CII and OVI lines corresponding to
chromosphere and transition region temperatures; these showed significant broad-
ening and complex line profiles in regions overlying the sites of small scale magnetic
elements in the photospheric network. Doppler shifted multiple peaks in CII line
were always seen soon after the reconnection of magnetic flux tubes occurs and
usually consist of supersonic and subsonic components caused by shocks propagat-
ing upward. Multiple peaks in OVI line have more diverse features: they are not
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 H2V 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 H2V grains in the overlying
chromosphere were investigated by Lites et al. (1999) [197].
Cauzzi et al. (2000) [59] analyzed the temporal behavior of Network Bright
Points (NBPs) using a set of data acquired during coordinated observations be-
tween ground-based observatories (mainly at the NSO/Sacramento Peak) and the
Michelson Doppler Interferometer onboard SOHO. The NBP’s were observed in
the NaD2 line and were found to be cospatial with the locations of enhanced mag-
netic field. The “excess” of NaD2 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 NaD2 line center emission can be used
as a proxy for magnetic structures.
70 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
In a paper by Fossum and Carlsson it was shown that acoustic heating of the
chromosphere is a factor of 10 too low to balance radiative losses (Fossum and
Carlsson, 2005 [99]).
Simultaneous CaII K-line spectroheliograms and magnetic area scans were used
to search for spatial correlation between the CaII K2V 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 > 4 Mxcm−2 (Sivaraman
et al. 2000 [286].
• Granulation: size about 1 000 km, lifetime 0.2 hr, vertical flow ∼ 1 kms−1 .
Figure 3.8: Summary of basic processes of most solar flares and their emission regions.
Courtesy: B. Vrsnjak.
bold arrow-lines (grey). Electron beams are given by thin arrows (and marked
by e− ). Chromospheric evaporation from the flare kernels is indicated by thick
dotted arrows. As it is seen in the sketch, the primary energy release takes place
in the corona at heights between 104 and 105 km. DCIM indicates fast drifting
bursts in the 200-2000 MHz range. HXR (Hard X-ray) emission is related to radio
features. HXR emission at successively lower energies indicates delays of slower
electrons relative to faster ones. The power of flares is related to the height of
the energy release site. Flares are more powerful and impulsive when the energy
release site is located at low heights. This can be explained by the weakening of
the magnetic field with height. As is also seen in the Fig. 3.8, electron beams that
are produced at the primary energy release site can escape outwards exciting type
III bursts. Electrons attached to the closed field lines become trapped between the
two magnetic mirrors located near the footpoints and they excite type IV bursts.
In very strong fields also µwave emission occurs (mm-cm range). Electrons
with small pitch angles penetrate through the magnetic mirrors. They hit the
dense transition region and chromosphere and excite line emission of atoms and
ions and hard X ray emission (HXR). This process is also called thick target
Bremsstrahlung. The chromospheric plasma is heated and starts to expand. This
is the evaporation process which continues until a new hydrostatic equilibrium is
reached. This is a source of soft X-ray emission (SXR), the plasma has a density of
∼ 10−3 cm−3 and a temperature of ∼ 107 K. The evaporation and SXR emission
is a cumulative effect of precipitating electrons- the cooling is relatively slow. The
SXR curve behaves as a time integral of the HXR curve. Or it can also be stated
that the HXR curve looks like the time derivative of the SXR curve. This is called
3.4. SOLAR FLARES 73
the Neupert effect (see also Neupert, 1968 [231] and Veronig et al., 2002 [326]).
• Soft x-ray classification: since 1970 flares are also classified based on soft
x-ray observations of the Sun in the 0.1-0.8 nm band by Earth orbiting
satellites. The size of the flare is given by the peak intensity (on a logarithmic
scale) of the emission (see Table 3.4).
According to Table 3.4 a B5 flare has a peak flux of 5 × 10−7 Wm−2 . Flares
smaller than C1 can only be detected during a solar minimum phase when
the general x-ray background is low. Occasionally, flares exceed class X9 in
intensity and are referred simply to as X10, X11...
Long duration flares are linked to coronal mass ejections (CMEs) but recent
observations also showed that some short duration flares may have ejecta. Coronal
mass ejections (CMEs) leave the Sun at speeds up to 2000 km/s and can have
angular spans over several active regions whereas flares imply events that are
localized within a single active region. In CMEs the magnetic field lines are opened
in eruptive events. There occurs a closing down or reconnection within several
hours providing a prolonged energy release that 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 Hα. 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 association
with geomagnetic storms.
Confined or impulsive events may also result from loop top magnetic recon-
nection. An impulsive flare of say 1024 J is typically spread over an area of several
1014 m2 in Hα. 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.
The Wind spacecraft7 observes radio bursts in the frequency range 1-14 MHz.
Standard patrols of bursts are made above 25 MHz. With the Bruny Island Radio
Spectrometer, this gap is filled and it is studied whether radio bursts can be used in
diagnosing energetic particle generation and propagation in the inner heliosphere
(Cane, Erickson, 2006 [54]).
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 IV emission is related to magnetic
reconnection in CME.
7 was launched in 1994, part of the ISTP project
3.4. SOLAR FLARES 75
Type II radio bursts result from plasma radiation associated with a MHD
shock propagating through the corona (∼500 km/s). This can be observed by a
slow drift emission. More than 90% of type II bursts have an associated flare.
They accompany 30% of flares with an Hα importance class 2 and 3. 70% of all
type II bursts are associated with a CME.
• Homologous flares: these are earlier flares in the same location with similar
emission patterns. They occur most often in periods of frequent flare activity.
The rate of repetition ranges from a few per hour to several days.
• Soft x-ray precursors: these are transient enhancements in soft x-rays lasting
for several minutes; they occur in loops or unresolved kernels or close to flare
sites. Weak soft x-ray bursts are often observed at the time of the onset of a
76 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
CME. Sometimes several tens of minutes prior to the impulsive phase. The
location is at one foot of a large coronal arch which already exists. The
process can be interpreted by a small magnetic structure which interacts
with the large coronal arch at one of its footpoints. The whole structure
becomes then destabilized.
• Radio precursors: often tens of minutes before the onset of a flare, changes in
intensity and polarity in microwaves are observed. However the correlation
with flares is not very strict.
• Prominence eruptions: very often they precede two ribbon flares. The time
delay between the onset of the prominence eruption and the impulsive phase
is of the order of minutes. Enhanced mass motion, a slow rise of the promi-
nence and untwisting can precede the main flare by hours.
Of course in all the cases joint observations covering the whole electromagnetic
spectrum are important. In the review given by Aschwanden et al. (2001) [14] the
authors focussed on new observational capabilities (Yohkoh, SoHO, TRACE).
The formation of a radio-emitting shock wave and its precursor above a flaring
active region was investigated in Klassen et al. (2003)[166]. They used imaging
and spectral observations of radio bursts with Yohkoh soft hard and X-ray imaging
observations and identified type II precursor as a signature of the reconnection
process above the expanding soft X-ray loops.
Characteristics of flare producing sunspot groups were discussed by Ishii et al.
(2000) [146]. A review about reconnection theory and MHD of solar flares is given
by Priest (2000) [249].
3.4.4 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 Hα. Over the photosphere they appear as dark filaments, at the
limb as bright structures.
The prominence plasma contains 90% of hydrogen which is partially ionized in
the central coolest parts of prominences where the temperatures are between 6000
and 8500 K or maybe even lower. At the boundary of prominences the temperature
rapidly increases to coronal values (more than 1 million K). The plasma density
in the central cool parts is about two orders of magnitude larger than that in the
corona.
3.4. SOLAR FLARES 77
These facts imply that the magnetic field is crucial for the prominence support
and stability. The intensity of the field ranges up to a few tens of Γ.
Some prominences are short lived eruptive events (variations within minutes
to hours), others can be quiescent and survive many rotational periods of the
Sun. The upper parts are often located in the hot corona. Quiescent prominences
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 105 km.
A quiescent prominence may change into an eruptive prominence. The typical
thickness of the loop is 104 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 1016...17 m−3 which is
a hundred times greater than coronal values.
Prominences are mostly located along the so-called neutral lines where the
vertical photospheric magnetic field changes its sign. Along the neutral line, the
vertical component of that field is zero.
A possible mechanism to understand cool prominence material (temperature
about 104 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,
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, g,
in the z-direction is:
P = ρT /µ (3.40)
where is the gas constant and µ the molecular weight. In a simple model
Kippenhahn and Schlüter (1957) [163] assumed that the temperature T and the
horizontal magnetic field components Bx , By were constant and that P, ρ and Bz
were functions of x alone. The prominence is represented as a plane sheet.
Tripathi et al., 2004 [316] studied an erupting prominence with EIT and then
with LASCO when it developed into a CME.
A recent review about prominences was given by Heinzel and Anzer, 2005 [131].
78 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
Table 3.6: Tomography of the solar corona by observations at different radio frequencies
ν MHz λ cm F T
10−22 Wm−2 Hz−1
30 1000 0.17 5.1×105
300 100 14.9 7.0×105
3000 10 69 31 000
30 000 1.0 1862 10000
300 000 0.1 113 200 5900
where the emission occurs, thus also the deeper the temperature. The values are
given for the quiet Sun (Landolt, 1981 [17]).
• Coronal loops are found around sunspots and in active regions in the corona.
These structures are associated with the closed magnetic field lines that
connect magnetic regions on the solar surface. As it is shown in the chapter
on MHD, in the corona the magnetic field dominates the motion of the
plasma, and therefore the plasma is aligned in magnetic loops. These loops
last for days or weeks. Some loops, however, are associated with solar flares
and are visible for much shorter periods. These loops contain denser material
than their surroundings. The three-dimensional structure and the dynamics
of these loops is investigated for that reason.
• 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 1083.0 nm 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) [175], and accordingly they are considered the main reason for the
80 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
Figure 3.10: Progress of a Coronal Mass Ejection (CME) observed over an eight hour
period on 5-6 August 1999 by LASCO C3. The dark disk blocks the Sun so that the
LASCO instrument can observe the structures of the corona in visible light. The white
circle represents the size and position of the Sun. Courtesy: SOHO/LASCO. SOHO is
an ESA/NASA mission.
Such a structuring is seen best when CMEs erupt close to the solar limb- then they
are seen from the side. Earth- (or oppositely) directed CMEs show an outflow and
expanding brightness around the Sun- these are called halo CMEs (see also Jackson
et al., 2002 [147], where it is discussed whether Halo CMEs will hit or miss the
Earth). Therefore, Halo CMEs are of special interest for space weather.
CMEs can be observed in white light10 . In white light we see photospheric
light scattered on coronal free electrons (Thomson scattering). The brighter the
structure, the more massive. Brightness does not mean temperature. They can
also be observed in other wavelengths, where near surface structures are observed
(Hα, He 1083 nm, EUV, X-rays, microwaves to radio).
How often do CMEs occur? SOHO observations11 yield the following frequen-
cies of CME occurrences:
The CME mass shows no cycle dependence, whereas the cycle influences their lati-
tudinal distribution: during minimum CMEs are concentrated around the equator,
during maximum they originate from a wide range of latitudes.
There exist two types of CMEs:
Flare associated CMEs are, on average, faster (median speed 760 km/s) than the
ones associated with filament eruption without flare (median speed 510 km/s).
The temperature is about 8000 K in the core and more than 2 million K in the
front part and in the cavity.
SOHO/LASCO data from 1996 to 2001 were collected by Yeh et al., 2005 [343]
and they showed that the observed CMEs reveal a similar power-law behavior as
flares, and the power-law indices for both phenomena are almost identical. This
finding strongly supports the viewpoint that solar flares and CMEs are different
manifestations of the same physical process.
CMEs are an important factor in coronal and interplanetary dynamics by in-
jecting large amounts of mass and magnetic fields into the heliosphere causing
geomagnetic and interplanetary shocks which is a source of solar energetic parti-
cles.
The geoeffectiveness of CMEs is reviewed in the paper by Webb, 2002 [332]
and Kim et al., 2005 [162]. They considered more than 7000 CMEs observed by
SOHO/LASCA and also 300 frontside halo CMEs between 1997 and 2003. The
geomagnetic storm that is associated with the CME was measured by the Dst
index (see next chapter). They found that the probability of front side CME
geoeffectiveness is 40%. For speeds >400 km/s and L < 500 the probability
of detection is high (80%) but also the false alarm rate is high (60%). The most
probable areas (or coverage combinations) whose geoeffectiveness fraction is larger
than the mean probability (about 40%), are 00 < L < +300 for slower speed CMEs
(≤ 800 km/s), and −300 < L < +600 for faster CMEs (>800 km/s). Manchester
et al., 2004 [205] gave a study of a numerical simulation of a CME propagating
form the Sun to 1 AU. They found that CME is very effective in generating strong
geomagnetic activity on Earth through a strong sustained southward Bz and by
a pressure increase associated with the CME driven shock that compresses the
magnetosphere.
A recent review about CMEs can be found in van Driel-Gesztelyi, 2005 [324].
• hot temperature
• low density (only about 10−12 that of the photosphere).
The original idea for the heating of the corona was entirely non-magnetic. From
laboratory experiments we know that if a fluid is set into violent motion, it emits
sound with the amount of sound rising as a high power of the average velocity of
the fluids. As we have seen, in the photosphere convective motions occur. If these
convective motions produce sound waves, they must propagate outwards from the
surface of the Sun. The wave motion has an energy density of
1 2
Ewave = ρv (3.41)
2
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 increases.
However, it turned out that a purely acoustic heating of the corona is not suffi-
cient to explain the high temperatures there. Acoustic heating may be important
in the outer layers of some stars.
Today12 we assume that the following two processes are the main reason for
the hot corona:
• MHD waves: as it has been outlined, when a magnetic field is present there
are two characteristic speeds of wave propagation, the sound speed cs and
the Alfvén speed cH . If cs >> cH 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 stressed that MHD waves have an anisotropic propagation.
• Magnetic reconnection: The footpoints of magnetic fields often are seen to
be anchored in the photosphere. In this region they are being continually
moved around by convective motions. Thus magnetic reconnection occurs
and electric currents flow which are dissipated.
There seems to be two problems with that interpretation. MHD waves cannot
carry enough energy through the chromosphere to the corona and Alfvén waves
dissipate their energy very fast when entering the corona. Bogdan et al., 2003 [41]
have shown that Alfvén waves can transmute into other wave modes at the base
of the corona.
The first observational evidence of the presence of waves in the corona was
made by SOHO EUV observations. Waves with a frequency of 1 mHz were found
but they could only contribute to about 10% of the needed energy. The photo-
sphere is covered by small magnetic elements (size below 1000 km). These small
elements are constantly perturbed by granulation motions. The magnetic field in
the corona that is anchored at these elements therefore constantly is perturbed
12 see also the book: Mechanisms of Chromospheric and Cornal Heating, P. Ulmschneider, E.
and reconnection occurs due to the motion of the magnetic carpet. Maybe a series
of mircoflares occurs (see e.g. Benz, 2003 [33]).
RHESSI studied gamma and X-ray emission from flares and mircoflares. Mi-
croflares emit hard X-ray and it turned out that mircoflares are quite similar to
normal flares. The type III radio bursts seem to be in relation with series of
microflares- the radio signals decrease in frequency like the whistle from a depart-
ing train. In type III bursts electrons are accelerated in open magnetic field lines
and the particles escape from the Sun. RHESSI observations of microflares were
found to be coincident with TRACE observations ( showing jets in the EUV) (see
Liu et al. 2004 [198]).
13 see e.g. the classical textbook: A. J. Hundhausen, Coronal Expansion and Solar Wind, 1972,
Springer
14 He reported this observation to the Royal Astronomical Society
3.6. SOLAR WIND AND INTERPLANETARY MAGNETIC FIELD 85
Ion tail
Dust tail
Hale Bopp
Figure 3.11: Comet Hale Bopp (1997); the fainter ion tail is clearly seen.
Satellite Measurements
The first in situ measurement of the solar wind was made in 1962 by Mariner 215
First we want to mention that besides SOHO two satellite missions measure
the solar wind: Ulysses and ACE. Ulysses was launched from the space shuttle
Discovery in 1990. The spacecraft made a journey to Jupiter where the giant
planet’s gravity pulled the spacecraft into a trajectory that carried it over the
Sun’s south pole in the fall of 1994 and its north pole in the summer of 1995. The
next passes over the Sun’s south pole occurred during 2000 and over the north
pole during 2001. These two orbital passes provide views of the solar wind at
times near the minimum of solar activity and the maximum of solar activity. It
was found that in 2000 the south magnetic pole almost completely vanished at
the time of solar maximum. In November 2006 Ulysses will continue with a third
south polar pass and beginning of December 2007 with its third north polar pass.
The solar wind speed, magnetic field strength and direction and composition
were measured.
The Advanced Composition Explorer (ACE) satellite was launched in August
of 1997 and placed into an orbit about the Lagrangian L1 point between the Earth
and the Sun16 . ACE has a number of instruments that monitor the solar wind.
The SOHO/SWAN experiment (Solar Wind ANisotropies) measures the Lα
radiation that is scattered by hydrogen atoms, which flow into the solar system.
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
15 Mariner 2 also detected the slow retrograde rotation rate of Venus, its surface temperatures
Planetary Magnetospheres
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 non magnetized planets
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 discon-
nections. 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) [265].
Krymskii et al. (2000) [176] investigate the interaction of the interplanetary
magnetic field and the solar wind with Mars. Data from the Mars Global Surveyor
mission have shown that localized crustal paleomagnetic anomalies are a common
feature of the Southern Hemisphere of Mars. The magnetometer measured small-
scale magnetic fields associated with many individual magnetic anomalies (mag-
nitudes ranging from hundreds to thousands nT at altitude above 120 km). Thus
Mars is globally different from both Venus and Earth. The data collected by Lunar
Prospector near the Moon were interpreted as evidence that above regions of in-
ferred strong surface magnetic fields on the Moon the solar wind flow is deflected,
and a small-scale mini-magnetosphere exists under some circumstances. With a
factor of 100 stronger magnetic fields at Mars and a lower solar wind dynamic pres-
sure (because of the greater distance), those conditions offer the opportunity for a
larger size of small ‘magnetospheres’ which can be formed by the crustal magnetic
fields. The Martian ionosphere is controlled both by solar wind interaction and
by the crustal magnetic field. Therefore, the nature of the Martian ionosphere is
probably different from any other planetary ionospheres, and is likely to be most
complicated among the planetary ionospheres (Shinagawa, 2000 [283]).
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 1.5 × 108 km from the
Sun, the interplanetary magnetic field makes an angle of about 45 degrees to the
radial direction. Further out19 the field is nearly transverse (i.e. about 90 degrees)
to the radial direction.
The Sun’s magnetic field, that is carried out into interplanetary space is called
the interplanetary magnetic field, IMF. The interplanetary field lines are frozen in
the plasma. Because of the Sun’s rotation, the IMF like the solar wind travels out
in a spiral pattern. This can be compared to the pattern of water sprayed from
a rotating lawn sprinkler. The winding up of the magnetic field is named Parker
spiral after the scientist who first described this (see Fig.3.12). Furthermore,
sectors (typically four) with alternating inward and outward directed magnetic
fields can be identified.
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.
19 at the orbit of Saturn
3.6. SOLAR WIND AND INTERPLANETARY MAGNETIC FIELD 89
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.
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
of the solar atmosphere. In such a case the pressure gradient in the corona must
be balanced by the gravitational attraction of the Sun:
dP GM
=− 2 (3.42)
dr r
In this equation we have replaced the variable M by M since the mass of the
corona is negligible to the total mass of the Sun. We also can write:
P = nkTkin ρ = nm (3.43)
n is the number of particles per unit volume and m is the average particle mass.
Please also note that Tkin 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
5/2
κ = κ0 Tkin (3.44)
where rc , Tc are radius and temperature at some point in the corona. Combining
all four above equations one gets P and n as a function of r. When expanding this
to the Earth one gets a kinetic temperature of 5 × 105 K and a particle density of
4 × 108 m−3 . Parker pointed out that a solution of such a system to the edge of the
solar system gives nonsense. At large values of r the value of P becomes constant,
so that ρ ∼ r2/7 . 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 vr in the radial
direction, then equation 3.42 together with 3.43 becomes
dvr d GnmM
nmvr = (nkTkin ) − (3.47)
dr dr r2
Mass conservation requires:
nr2 vr = const (3.48)
90 CHAPTER 3. THE SOLAR ATMOSPHERE AND ACTIVE REGIONS
F = ρ∆v 2 vA A (3.49)
√
ρ... mass density, ∆v 2 wave velocity amplitude, vA Alfvén speed and A is the
cross section of the flow tube. The magnetic flux Φ = BA is constant, so that the
wave velocity amplitude changes with density as
∆v 2 = ∆v02 ρ/ρ0 (3.50)
Figure 3.13: The interplanetary current sheet (above) and the heliospheric current sheet
(below). Courtesy: [Link]
emitted. These energetic photons are absorbed and re-emitted in the solar interior
(mean free path between the absorption processes is only 1 cm) and therefore for
a photon generated by such nuclear reaction, it takes more than 1 million years to
diffuse out of the core region. Thus one can argue that the emergent luminosity at
the core outer boundary is effectively constant on solar cycle timescales. If there
is any luminosity variability at the surface there must be an intermediate energy
reservoir between the core and the photosphere. There are several mechanisms for
storing energy during a solar activity cycle, such as magnetic fields or gravitational
energy. Each of them leads to distinct perturbations in the equilibrium structure
of the sun. Therefore, one can argue that a sensitive determination of the solar
radius fluctuations can help to understand the solar cycle and it is clear that the
magnitude of the radius fluctuations compared to the luminosity change contains
information on where and how energy is stored.
Sofia and Endal (1979) [291] introduced the parameter W by:
δr δL
W = / (3.51)
r L
The models predict a wide range for W :
In the following we discuss briefly some measurement methods and give the
results.
3.7. VARIATIONS OF THE SOLAR DIAMETER 93
• Correlations with solar activity: nearly opposing trend; high activity means
smaller diameter.
It is important to notice that all ground based observations must take into
account the quasi biennial oscillation in the Earth’s atmosphere.
Other astrolabe measurements were done by Sanchez et al. (1995) [268] at the
San Fernando Observatory (Cadiz), Noel (1995) [233] at Santiago and Leister et
al. at Sao Paulo Observatory (1990) [193]. Ribes et al. (1991) [256] report on
photoelectric measurements made at HAO in Boulder from 1986-1990; Wittmann
et al. (1993 [340]) report on measurements using drift timing in Izana and Locarno.
Other authors measured the solar diameter from eclipse data (e.g. Kubo, 1993)
[177].
The French CNES plans to operate the satellite mission PICARD which is
supposed to operate in 2008 and the satellite will be launched with a Russian
Dnepr rocket. One of the experiments of this low mass mission is SODISM (Solar
Diameter Imager and Surface Mapper).
Chapter 4
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 cycle1 .
We sometimes mix the unit system. This was done by intention because in
some cases it is easier and more intuitive to deviate from the SI system.
Let us give a simple example. The Coulomb law is defined as2 :
q1 q 2
F =k (4.1)
r2
In the SI system k = 1/4πe ; 0 = 8.854 × 10−12 C2 N−1 m−2 . In the cgs system
k=1.
Cambridge, 1991
2 it defines the force acting between two charges q , q
1 2
97
98 CHAPTER 4. MHD AND THE SOLAR DYNAMO
Figure 4.1: Looped magnetic field lines in the solar chromosphere and corona. Photo:
TRACE mission, NASA
Equation 4.2 is the Ampere law which states that a spatially varying magnetic
field (given by ∇ × H) produces currents- in MHD often the variation of the E
fields is neglected, thus the term ∂D∂t → 0. We can also state that a current j
induces a magnetic field that is in a direction opposite to it.
Equation 4.3 is the Faraday law and states that a time varying magnetic field
produces an electric field.
If µ0 , 0 are the permeability and permittivity of free space, then for most
gaseous media in the universe:
B = µ0 H D = 0 E (4.6)
The following equation relates the electric current density to the fields producing
it (generalized Ohm’s law):
j = σ(E + u × B) (4.7)
σ is the electrical conductivity and u 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
ion:
j = ni Zi eui − ne eue ρE = ni Zi e − ne e (4.8)
ni , ui , ne , ue are the number density and velocity of the ions and electrons respec-
tively and Zi e, −e are the charges on the ion and the electron.
In astrophysics two simplifications are applied:
4.1. SOLAR MAGNETOHYDRODYNAMICS 99
The third Maxwell equation (4.4) states that there are no magnetic monopoles3 .
Electric fields can be produced by separating positive and negative charges through
the fourth Maxwell equation (4.5) 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 ∂D/∂t are
usually insignificant in astrophysical problems because electric fields are unimpor-
tant; however they can be produced by a conduction current j, if the electrical
conductivity is high enough. Such magnetic fields may be slowly variable in time
and space.
We therefore neglect ∂D/∂t, combine the equations
obtaining
∂B 1
+ ∇×∇×B=0 (4.10)
∂t µo σ
and using ∇ × ∇ × B = grad divB − ∇2 B and divB = 0:
∂B 1
= ∇2 B (4.11)
∂t µ0 σ
This is also called the induction equation- for the static case. In cartesian
coordinates this equation for the x coordinate is:
2
∂Bx 1 ∂ Bx ∂ 2 By ∂ 2 Bz
= + + (4.12)
∂t µ0 σ ∂x2 ∂y 2 ∂z 2
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 that currents vary significantly in distance L, then from (4.11) the decay
time becomes
τD = µ0 σL2 (4.13)
If at time t = 0 there exists a sinusoidal field
Bx = B0 exp(iky) (4.14)
The wavelength λ of the spatial variation of the field is 2πk, the original field
decays by a factor e in the time µ0 σλ2 /4π 2 .
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 ∼ T 3/2 . 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
1.5 ×107 K). 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 ρT /µ, where is the gas constant and µ the
mean molecular weight. With Tin = Tout we must have:
A tube of magnetic flux is lighter than its surroundings and will start to rise which
is called magnetic buoyancy.
To a good approximation, the fluid is tied to the magnetic field. For the Sun
two extremes occur:
• photosphere: the fluid motions drag the magnetic field lines around (the
magnetic field is frozen in);
4.1. SOLAR MAGNETOHYDRODYNAMICS 101
• corona: the magnetic force is so strong that it constrains the motion of the
fluid.
The tying of the fluid to the magnetic field lines also permits the propagation
of MHD waves which have some similarity to sound waves but a characteristic
speed (Alfvén speed):
cH = B 2 /µ0 ρ (4.19)
The sound speed is given by
cs = γP/ρ (4.20)
This can also be seen from the induction equation. Let us consider again the
Maxwell equations. From j = σE + u × B we can extract E:
j
E= −u×B (4.21)
σ
This is substituted into the Maxwell equation (4.3) yielding:
j ∂B
∇× −u×B =− (4.22)
σ ∂t
We have already argued that the displacement current can be neglected in the first
Maxwell equation and therefore ∇ × B = µj, from which j = 1/µ∇ × B and
1 ∂B
∇× ∇×B−u×B =− (4.23)
µσ ∂t
∂B
= ∇ × (u × B) + η∇2 B (4.25)
∂t
Here η = 1/µσ is the magnetic diffusivity. The case where the plasma is stationary
was already discussed above (u = 0), the field decays in the ohmic decay time
τ = L2 /η. Let us discuss the case when η = 0. Then, the field B is completely
determined by the plasma motions u 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:
Φ= [Link] (4.26)
S
102 CHAPTER 4. MHD AND THE SOLAR DYNAMO
If G is the material closed curve bounding S, the total rate of change of Φ is (see
also eq. 4.34):
DΦ ∂B
= .dS + B.(u × dl) (4.27)
Dt S ∂t
G
∂B
= .dS + (B × u).dl (4.28)
S ∂t
G
∂B
= − ∇ × (u × B) .dS (4.29)
S ∂t
If η = 0 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 v0 , l0 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
Rm = l0 v0 /η0 (4.31)
In an active solar surface region one has η0 = 1 m−2 s−1 , l0 = 700 km ∼ 1 arcsec
and v0 = 104 m/s we find Rm = 7 × 109 >> 1. Thus the field is frozen to the
plasma and the electric field does not drive the plasma but is simply E = −u × B.
However, if the length-scales of the system are reduced the diffusion term η∇2 B
becomes important. Then the field lines are allowed to diffuse through the plasma
and this yields to magnetic braking and changing the global topology of the field
(magnetic reconnection).
structure is buoyant because the particle density is lower there since it contains
larger magnetic energy density (see magnetic buoyancy). Thus the external pres-
sure 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 reconnec-
tion 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 [37]).
A recent review on solar MHD was given by Walsh (2001) [329].
Dρ
+ ρ∇.u = 0 (4.33)
Dt
and the total derivative means here:
D ∂
= + u.∇ (4.34)
Dt ∂t
(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 u: the momentum
equation includes the Lorentz force term j × B and other forces F, such as gravity
and viscous forces:
Du
ρ = −∇p + j × B + F (4.35)
Dt
104 CHAPTER 4. MHD AND THE SOLAR DYNAMO
Here p is the plasma pressure. Let us assume a Newtonian fluid with isotropic
viscosity, then F may be written as:
r
F = −ρg(r) + ρν∇2 u (4.36)
r
g(r) 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 r from the rotation axis:
Du dΩ 1 2
ρ = −∇p + j × B + F + ρ 2u × Ω + r × + ∇|Ω × r| (4.37)
Dt dt 2
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 u × Ω.
Ds
ρT = −L (4.39)
Dt
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:
ργ D p j2
= −∇.q + κr ∇2 T + +H (4.40)
γ − 1 Dt ργ σ
• T temperature
0 = −∇p + j × B + F (4.41)
This equation must then be solved with ∇ × B = ..., ∇.B = 0 and the ideal gas
law as well as a simplified form of the energy equation.
Let us introduce the concept of scale height. Let
dp
0=− − ρg (4.42)
dz
Substitute in the above equation ρ = pm/kT (ideal gas) and integrate:
z
dz
p = p0 exp − (4.43)
0 Hp (z)
(p0 is the pressure at z = 0). This defines the local pressure scale height Hp :
dp
0=− − ρg cos θ dz = ds cos θ (4.45)
ds
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 p0 to magnetic pressure
B02 /2µ. If β << 1, any pressure gradient is dominated by the Lorentz force and
(4.41) reduces to:
j×B=0 (4.46)
In this case the magnetic field is said to be force free. In order to satisfy (4.46)
either the current must be parallel to B (Beltrami fields) or j = ∇ × B = 0. 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 B = B(x)k,
then (4.41) becomes:
∂ B2
0= p+ (4.47)
∂x 2µ
106 CHAPTER 4. MHD AND THE SOLAR DYNAMO
∇2 A = 0 (4.48)
1 d2 X 1 d2 Z
2
=− = −n2 (4.49)
X dx Z z2
where n = const. A solution to (4.49) would be:
B0
A= sin(nx)e−nz (4.50)
n
this gives for the field components:
∂A
Bx = = B0 cos(nx)e−nz (4.51)
∂x
∂A
Bz = = −B0 sin(nx)e−nz (4.52)
∂z
The result is a two dimensional model of a potential arcade.
B0 = B0 b
E0 = E b + E⊥
v = v b + v ⊥
Here v0 is the velocity component at t = 0 in the direction of the magnetic field
line. We see that particles of opposite sign of charge q move in opposite directions,
they move along an electric field parallel to a magnetic field which destroys E .
108 CHAPTER 4. MHD AND THE SOLAR DYNAMO
The difference in mass between electrons and protons is about 1800. What follows
for the radius and frequency of gyration?
The gyration follows also from the simple statement, that in the absence of
other forces, the Lorentz force balances the centripetal force of the particle’s motion
around the field line. Let α be the pitch angle, which is the angle between the
direction of motion and the local field line. Then for the gyro radius (Larmor
radius :
cp⊥ mcv × B
RL = = (4.68)
qB qB 2
The Larmor radius for a 100 keV electron (which is typical for electrons in the
inner radiation belt of the Earth) is about 100 m.
Let us assume, that the magnetic flux through a particle’s orbit is constant-
this is certainly the case when changes of the magnetic field are small over the gyro
radius and one gyro period. From the condition that dΦB /dt = 0 the so called
first adiabatic invariant follows:
p2⊥
µB = (4.69)
2mB
or in terms of the particle’s energy:
E sin2 α
µB = (4.70)
B
From the conservation of the first adiabatic invariant it follows, that the pitch
angle increases, when the particle moves to larger field strength, until α = 900 at
the mirror point.
Summarizing the motion of a particle:
4.1. SOLAR MAGNETOHYDRODYNAMICS 109
Figure 4.3: Left: Gyration of a charged particle around magnetic field lines, the sense
depends on the + or - charge; right: charged particles are deflected by magnetic fields,
depending on their charge. Here the magnetic fieldlines point vertically inwards, denoted
by
F = −µB ∇B (4.72)
110 CHAPTER 4. MHD AND THE SOLAR DYNAMO
F = qE (4.76)
We stress here, that the ∇B and E × B drift and curvature forces dominate
the drift motions of particles in a magnetosphere.
Let us consider a large assembly of particles; these particles interact which
is called collision. If τc is the characteristic time between collisions the collision
frequency is νc = 1/τc . If νc is large, the particle motions will be disordered
and decoupled from the magnetic field, the fluid will not be tied to the field. If
4.1. SOLAR MAGNETOHYDRODYNAMICS 111
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:
σ = ne e2 τc /me (4.78)
and τc ∼ T 3/2 .
• fast and slow magnetosonic waves; their dispersion relation is given by:
Let us consider two special cases: if the waves propagate along the field k.b = k,
there are two waves with ω 2 = k 2 c2H and the sound wave ω 2 = k 2 c2s unaffected by
the field. For wave propagation perpendicular to the field only one wave survives
with ω 2 = k 2 (c2s + c2H ). When waves propagate anisotropically, it is necessary
to introduce another wave velocity in addition to the phase velocity, the group
velocity, given by ∂ω/∂k with which the wave carries energy or information. The
group velocity for Alfvén waves is always cH b. What does that mean? Regardless
of the direction in which it propagates, energy always travels along the field lines
with speed cH .
• 11 year period of the sunspot cycle; not only the number of sunspots varies
over that period but also other phenomena such as the occurrence of flares,
prominences,.... etc.
• the equator-ward drift of active latitudes which is known as Spörers law and
can be best seen in the butterfly diagram. At the beginning of a cycle active
regions appear at high latitudes and toward the end they occur near the
equator.
• 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.
Figure 4.4: Illustration of the ω effect. The field lines are wrapped around because of
the differential rotation of the Sun
The α Effect
However, the field lines are not only wrapped 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
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.
Figure 4.5: The MHD relation between flows and magnetic fields
the upper radiative zone resulting in the observed deficit of Li and Be. Gilman,
2005 [112] wrote a summary about the tachoclyne stressing its importance for in
situ generation of poloidal fields as well as creating magnetic patterns that are
seen on the surface.
The solar meridional flow is a flow of material along meridional lines from the
equator toward the poles at the surface and from the poles to the equator deep
inside. At the surface this flow is in the order of 20 m/s, but the return flow
toward the equator deep inside the Sun must be much slower since the density is
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).
∂B
= ∇ × (u × B) − ∇ × (ηe ∇ × B) (4.89)
∂t
116 CHAPTER 4. MHD AND THE SOLAR DYNAMO
The flow u 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)
[174]. The velocity and the field are expressed as:
< u >, < B > represent slowly varying mean components and u , B non axisym-
metric fluctuating components. The turbulent motion u is assumed to have a
correlation time τ and a correlation length λ which are small compared to the
scale time t0 and scale length l0 of the variations of < u > and < B >. In
other words, τ is a mean time after which the correlation between u (t = τ ) and
u (t = 0) is zero and λ is comparable to the mean eddy size. We assume that
< u >, < B >= 0.
This is substituted into the induction equation and subtracted from the com-
plete equation:
∂B
= ∇ × (< u > ×B + u × < B > +G) − ∇ × (η∇ × B ) (4.91)
∂t
where
E =< u × B > G = u × B − < u × B > (4.92)
E 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
B and here several assumptions are made. First of all we stressed that < v >= 0.
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: G ∼ 0. That is valid only if B << < B >. Then our
equation reduces to:
∂B
+ ∇ × (η∇ × B ) = ∇ × (u × < B >) (4.93)
∂t
We want to determine E. Thus only B the component of B which is correlated
with u must be considered. By definition τ, B(t + τ ) is not correlated with B(t)
for any t. B (t) may be determined by integration of the above equation from
t − τ to t. Note also, that the order of the convective turn over time τ ∼ λ/v
and thus both u and < B > may be regarded as independent of t. Thus the
integration yields:
∂ < Bj >
Ei = αij < Bij > +βijk (4.94)
∂xk
4.2. THE SOLAR DYNAMO 117
where αij , βijk depend on the local structure of the velocity field and on τ . If the
turbulent field is isotropic, then αij = αδij , βij = βijk , and
E = α < B > −β∇× < B > (4.95)
If τ is small compared to the decay time λ2 /η, the diffusive term may be neglected
and from 4.93 we get
1 1 2
α = − τ < u .∇ × u >, β= τv (4.96)
3 3
And finally:
∂B
= ∇ × (αB + u × B) − ∇ × [(η + β)∇ × B] (4.97)
∂t
Compared to the normal induction equation, this contains the term αB and the
eddy-diffusivity coefficient β. In the mean field dynamo, the magnetic diffusivity
η is replaced by a total diffusivity η‘ = η + β and the equation becomes:
∂B
= ∇ × (αB + u × B) + η‘∇2 B (4.98)
∂t
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 B is axisymmetric.
Then it can be represented by its poloidal and toroidal components A(x, z, t) and
B(x, z, t) and u = u(x, z, t)j. Neglecting the advection terms:
∂
− η∇2 B = [∇u × ∇A].j − α∇2 A (4.99)
∂t
∂
− η∇2 A = αB (4.100)
∂t
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. ∇u 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 α − ω-dynamo. 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 vc .∇ × vc of the non
axisymmetric cyclonic convection generates an azimuthal electromotive force
E which is proportional to the helicity and to Bφ .
Let us define a characteristic length scale l0 , a decay time t0 = l02 /η and u =
s0 ω, where s0 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 t = t/t0 and r = r/l0 . By an elimination of B and neglecting the α2
terms we arrive at
2
∂ 2 αl2 s0
− ∇ A = 02 [∇‘ω × ∇‘A].j (4.101)
∂t η
118 CHAPTER 4. MHD AND THE SOLAR DYNAMO
If α0 and ω0 are scale factors giving the orders of magnitudes of α and |∇‘ω| then
2
∂ α ∇ ω
− ∇2 A=D × ∇ A .j (4.102)
∂t α0 ω0
αω0 l02 s0
D= (4.103)
2η 2
Figure 4.6: 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 R located above the tachocline and
that at 0.63 R 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)
Calvo et al. (1995) [53] used the neural network technique to analyze the time
series of solar activity (given by the Wolf number).
Hernandez (1993) [132] also used neural nets to construct nonlinear models to
forecast the AL index (auroral electrojet index) given solar wind and interplanetary
magnetic field (IMF) data.
Gleisner and Lundstedt (2001) [113] used a neural network-based model for
prediction of local geomagnetic disturbances. Boberg et al. (2000) [38] made real
time Kp predictions from solar wind data using neural networks.
physics, 2000
6 See also the new edition of Solar and Stellar Activity Cycles, P. R. Wilson, A. King, Cam-
bridge 2004
4.3. STELLAR ACTIVITY 121
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.
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.
How can we measure stellar parameters like differential rotation that play a
key role in the onset of stellar dynamos? Let us assume we have a rapidly ro-
tating spotted cool star and that it is observed one week apart. By comparing
brightness/magnetic images of that star over such time intervals one can measure
the rotation rates of starspots at different latitudes over several rotation cycles
(Barnes et al. 2001) [26].
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) [124]. So far we have considered
only stars which have an activity level by orders of magnitude larger than the Sun.
123
124 CHAPTER 5. THE EARTH’S ATMOSPHERE AND CLIMATE
the Sun (covering all wavelengths from far UV to IR), and b) from IR radiation
reflected at the surface of the Earth.
In Table 5.2 we give the energy received from the Sun at a distance of 1 AU
in different wavelength regions. It is clearly seen, that 98% of the radiation from
the Sun is in visible to the near IR.
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 H2 O. The surface cools when water evaporates and heat is transmitted
to the air as vapor which recondenses to form clouds. Such phase transitions of
H2 O play a major role in the energy budget of the lower atmosphere.
5.1.2 Composition
The composition of the Earth’s atmosphere is given in Table 5.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 O2 , CO2 , CH4 , H2 . Oxygen
is produced by photosynthesis:
In this formula CH2 O denotes any variety of organic compounds. Aerobic respi-
ration and decay occur in the reverse reaction:
Volcanic gases
In the absence of this reaction, carbon would accumulate in organic form and the
fuel for photosynthesis (atmospheric CO2 ) would be depleted. If the supply of O2
is limited such as in the sediments of organic rich swamps and in the stomachs of
ruminants, we get as a product methane CH4 .
In Table 5.3 the change of the greenhouse gas and other gas concentrations of
the Earth’s atmosphere is given.
In Table 5.3 the present tropospheric concentration estimates are calculated
as annual arithmetic averages; ppm = parts per million (106 ), ppb = parts per
billion (109 ), ppt = parts per trillion (1012 ).
The Global Warming Potential (GWP) is generally used to contrast different
greenhouse gases relative to CO2 . 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: n
ai ci dt
GW P = n 0 (5.3)
a
0 CO2 CO2
c dt
where ai is the instantaneous radiative forcing due to a unit increase in the con-
centration of trace gas, i, ci is concentration of the trace gas, i, remaining at
time, t, after its release and n is the number of years over which the calculation
is performed. This formula is taken from the Intergovernmental Panel on Climate
Change (IPCC) (Houghton et al., 1990 [138]).
How was the concentration of greenhouse gases in the past? Ice cores, cylinders
of ice drilled out of glaciers and polar ice sheets, have played an important role
in to answer such questions. The drilling for The Greenland Ice Sheet Project
5.1. THE EARTH’S ATMOSPHERE 127
The measurements are from: 1 in situ air samples collected at Mauna Loa Observatory,
Hawaii (Bacastow et al.1985 [19]). 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., 1988 [92]
Two began in 1989, more than 3000 m deep. In 1992 there were data available to
reconstruct the climate over the past 200 000 years.
The CO2 data are from an ice core analyzed by Neftel et al. (1985) [230]. An
example of their measurements is given in Table 5.4.
These measurements of the CO2 gas concentration enclosed in an ice core from
Siple Station, Antarctica, indicate that atmospheric CO2 concentration around
1750 was 280 ± 5 ppmv (parts per million per volume) and has increased since,
essentially because of human factors, by 22.5 percent to 360 ppmv around 2000.
The anthropogenic emission of CO2 is about 7 Gt/yr. The natural and anthro-
pogenic changes in atmospheric CO2 over the last 1000 years from air in Antarctic
ice and firn was described in Etheridge et al. (1996) [93].
5.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.
A better knowledge about the climate and its variations in the past will enable
us to better understand what forces climate and its variations in the future. Since
there exists only a 140 years instrumental record, we have to use proxies to re-
construct climate in the past. Some widely used proxy climate data are (see 5.2):
128 CHAPTER 5. THE EARTH’S ATMOSPHERE AND CLIMATE
Table 5.4: Historical CO2 record from the Siple Station Ice Core
Depth [m] Samples Date of ice Date of Air CO2 conc. (ppmv)
measured enclosed in extracted air
187.0-187.3 10 1663 1734-1756 279
177.0-177.3 10 1683 1754-1776 279
162.0-162.3 9 1723 1794-1819 280
147.0-147.2 10 1743 1814-1836 284
128.0-129.0 47 1782 1842-1864 288
111.0-112.0 26 1812 1883-1905 297
102.0-103.0 26 1832 1903-1925 300
92.0-93.0 25 1850 1921-1943 306
82.0-83.0 28 1867 1938-1960 311
76.2-76.6 11 1876 1947-1969 312
72.4-72.7 11 1883 1954-1976 318
Paleoclimate reconstruction
Ocean sediments,
cave deposits
Lake levels,
mountain glaciers
Varved sediments,
lake sediments
bubbles, or isotopes of oxygen, that can be used to interpret the past climate
of that area2 . Let us briefly discuss one example of isotope measurements: Of
the temperature dependent markers the most important is the ratio of 18 O
to 16 O. This can be explained by the fact that water molecules composed of
H18
2 O evaporate less rapidly and condense more readily than water molecules
composed of H16 2 O. Thus, in the ice cores one obtains annual layers starting
with 18 O rich, becoming 18 O poor, and ending up 18 O rich.
• 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
versa. Accumulation rates inferred in this way are supported by measurements of
the cosmogenic isotope Beryllium 10 (10 Be), an isotope produced by the interaction
2 see also the textbook of R.B. Alley, The Two-Mile Time Machine, Princeton Univ. Press,
2002
130 CHAPTER 5. THE EARTH’S ATMOSPHERE AND CLIMATE
of cosmic rays and the upper atmosphere, can be used to determine past snow
accumulation in Vostok ice. Deposition of 10 Be is assumed to be constant. The
other two elements which are important are 18 O and deuterium. In Antarctica,
a cooling of 10 C results in a decrease of 9 per mil deuterium. The last ice age
is characterized by three minima separated by slightly warmer episodes called
interstadials.
Air initially enclosed in Vostok ice provides our only record of variations in the
atmospheric concentrations of CO2 and CH4 over a complete glacial-interglacial
cycle. For both greenhouse gases, concentrations were higher during interglacial
periods than during full glacial periods.
Crowley (2000) [70] discussed the causes of climate change over the past 1000
years. His main conclusion is that as 41-64% of pre-anthropogenic (pre-1850)
decadal-scale temperature variations were due to changes in solar irradiance and
volcanism.
Several periods of warmth (listed below) have been hypothesized to have oc-
curred in the past. However, upon close examination of these warm periods, it
becomes apparent that these periods of warmth are not similar to 20th century
warming for two specific reasons: a) the periods of hypothesized past warming do
not appear to be global in extent, or b) the periods of warmth can be explained by
known natural climatic forcing conditions that are uniquely different than those
of the last 100 years.
Examples of periods of warmth:
• Medieval: ∼ 9th to 14th centuries; this seems to be in doubt now because
the temperature anomaly at that time was very small; however the Little Ice
Age for the northern hemisphere from 15th to 19th centuries is clearly seen
(Fig. 5.3).
• mid-Holocene warm Period (approx. 6 000 years ago); this seems to be in
connection with changes of the Earth’s orbit (Theory of Milankovich).
• Penultimate interglacial period (approx. 125 000 years ago). It appears
that temperatures (at least summer temperatures) were slightly warmer than
today (by about 1 to 20 C), caused again by the changes in the Earth’s orbit
(Hughes and Diaz, 1994 [142]).
• Mid-Cretaceous Period (approx. 120-90 million years ago): Breadfruit trees
apparently grew as far north as Greenland (55◦ N), and in the oceans, warm
water corals grew farther away from the equator in both hemispheres. The
mid-Cretaceous was characterized by geography and an ocean circulation
that was vastly different from today, as well as higher carbon dioxide levels
(at least 2 to 4 times higher than today).
Figure 5.3: Temperature anomaly clearly showing the Little Ice Age. The Maunder
minimum is marked. (redrawn from Lamb, 1977 [185])
toward the Sun, at the same time the opposite hemisphere is in winter. If the
Earth’s axis were not inclined every point on the earth would receive the same
amount of sunlight each day of the year. Changes in this tilt can change the
severity of the seasons. More tilt means more severe seasons, i.e. warmer summers
and colder winters. The tilt of the Earth’s axis changes between 22 and 25 degrees
on a cycle of about 41 000 years. If the summers are cool snow and ice last
from year to year in high latitudes building up massive ice sheets. Now positive
feedbacks in the climate system start to work. Snow reflects more of the sun’s
energy into space causing additional cooling. Also the amount of the greenhouse
gas CO2 in the atmosphere falls as ice sheets grow and thus adding to the cooling.
Another astronomical effect on climate is that the orbit of the earth is not
circular. Presently, perihelion (closest approach to the Sun) occurs in January, thus
on the northern hemisphere winters are slightly milder. The perihelion changes in
a cycle of 22 000 years. Therefore, 11 000 years ago perihelion occurred in July
making seasons more severe than today. The eccentricity of the earth’s orbit varies
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
650 north (65N) (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.
132 CHAPTER 5. THE EARTH’S ATMOSPHERE AND CLIMATE
Figure 5.4: Cretaceous climate and land/sea distribution. During the Late Creta-
ceous the global climate was warmer than today’s climate. No ice existed at the
Poles. Dinosaurs migrated between the Warm Temperate and Cool Temperate Zones
as the seasons changed. The sea level was about 100 m higher than today. Courtesy:
http : //[Link]/timespace/
• 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.
5.1. THE EARTH’S ATMOSPHERE 133
Figure 5.5: Upper curve: average insolation of 65 degrees northern latitude (Watts per
one square meter of a horizontal atmosphere) in mid-July. As seen, it varies from some
390 to 490 W/m2 . Middle curve: Global temperature (Vostok ice core). Lower Curve:
Greenland, GRIP core. Image courtesy: Jan Hollan
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) [299]. Ice core results also enable the reaction of the environment
to past global changes to be investigated. The deglaciation of the northern hemi-
sphere is described in Alley and Clark (1999 [5]). 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) [308].
5.1.5 Greenhouseffect
Trace constituents of the atmosphere such as H2 O, CO2 , O3 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
134 CHAPTER 5. THE EARTH’S ATMOSPHERE AND CLIMATE
Figure 5.6: Vostok Ice core. Different depth can be attributed to different age
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.
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. 5.8). Similar greenhouse effects arise from the
gases methane, nitrous oxide and chlorofluorocarbons (CFCs). All these gases are
referred to as greenhouse gases due to their ability to trap heat.
The variation of greenhouse gases was described before.
5.1.6 Ozone
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) as the physical
thickness would be if compressed in the Earth’s atmosphere. The ozone layer is
very thin a normal range is 300 to 500 Dobson units.
We can make the simplification that throughout the stratosphere all of the
radiative energy from the sun that is absorbed by O3 is converted locally to heat.
The heating rate depends on the distribution of Ozone with height and on the
incoming solar energy.
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 (200 < λ < 300 nm). The absorption process results in
the dissociation of O3 .
5.1. THE EARTH’S ATMOSPHERE 135
400
CO2 [ppmv]
300
200
100
0
0 50 100 150 200
Years BP
Methane
800
CH4 [ppbv]
600
400
200
0
0 50 100 150 200
Years BP
Temperature
Years BP
Figure 5.7: Variation CO2 and CH4 in parallel with temperature from Vostok climate
records. Years BP are given in units of 1000. Credits: National Ice Core Laboratory.
The Chapman reactions describe the formation of ozone. First the following
reaction leads to a photodissociation of oxygen:
O2 + hν → O + O λ < 175 nm (5.4)
Then a recombination of oxygen occurs:
1. by direct two body reaction (reaction very rare!)
O + O → O2 + hν (5.5)
40000
Emission in Mill. tons
35000
30000
Industr. countries
25000
Eastern Europe, FSU
20000
Devel. countries
15000
Total
10000
5000
0
1980 1990 2000 2010 2020 2030
Year
Figure 5.8: World Carbon Dioxide Emissions. It is seen that the CO2 emission from the
Eastern European Countries and the former Soviet Union (FSU) declines, e.g. further
restructuring of the coal mining industries in Poland and the Czech Republic (US Dept.
of Energy)
O + O2 + M → O3 + M (5.7)
H + O3 → OH + O2 (5.10)
NO + O3 → NO2 + O2 (5.11)
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
5.1. THE EARTH’S ATMOSPHERE 137
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 radiation4 . 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.
So far we have discussed ozone in the stratosphere. In the Earth’s lower at-
mosphere, near ground level, ozone is formed when pollutants emitted by cars,
power plants, industrial boilers, refineries, chemical plants, and other sources re-
act 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 temperatures, normally oc-
cur. 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.
CO + OH → CO2 + H
H + O2 + M → HO2 + M
HO2 + M → OH + NO2
hν + NO2 → NO + O
O + O2 + M → O3 + M
CO + 2O2 → CO2 + O3
• 100 nm < λ < 200 nm: absorption of solar radiation leads to a dissociation
of O2 :
hν + O2 → O + O (5.12)
hν + O → O+ + e
hν + O2 → O+
2 +e
hν + N2 → +
N2 + e
4 UV radiation is divided by wavelength into UVA (320-400nm), UVB (290-320nm) and UVC
(100-290nm).
5 Chemistry of the Upper and Lower Atmosphere, Barbara J. Finlayson-Pitts, James N. Pitts,
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:
e + O+
2 → O+O
e + NO+ → N+O
O + O2 → O3 + M
H + O3 → OH + O2
O + OH → O2 + H
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 ∆z then the work mg∆z is done. m
denotes the mass of the particle, g the gravitational acceleration =9.81 m/s2 . The
work that must be done to escape the gravitational field is mgR, where R 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:
3
Ekin = kT (5.13)
2
where k is the Boltzmann constant (1.38 × 10−16 erg K−1 ). Thus an atom can
escape the gravitational field if its thermal kinetic energy ∼ kT is much larger
than mgR. Of course we must also consider collisions (except at the highest level
in the atmosphere). At the high temperatures in the thermosphere (700...2 000 K),
significant numbers of hydrogen atoms have velocities above the escape velocity
vesc ∼ 11.2 km/s. Therefore, hydrogen is lost at a rate of 108 atoms per cm2 per
second. These escaping hydrogen atoms are derived mainly from the oceans and
over the past 4.5×109 years of the Earth’s history, the sea level has declined by two
meters globally. Of course during this reaction also O2 is set into the atmosphere
which was crucial for the evolution of life.
There is also a significant loss of helium.
5.2. EARTH’S HISTORY AND ORIGIN OF THE ATMOSPHERE 139
a very active planet from the geologic point of view. Volcanism was widespread
and if we assume that volcanoes five billion years ago emitted the same gases as
they do today, the earth’s second atmosphere probably consisted of water vapor,
carbon dioxide, and nitrogen. These gases were expelled from the earth’s interior
by a process known as outgassing.
It is also possible that the impact of comets brought significant amounts of
water and other volatile gases to the Earth. The vast amounts of water vapor
expelled by the volcanic earth resulted in the formation of clouds which, in turn,
produced rain. Over a period of thousands of years, the rain accumulated as
rivers, lakes and ocean basins. This process was extremely important for the
carbon dioxide CO2 . The water reservoirs acted as sinks for that gas and through
chemical and later biological processes it became locked up in sedimentary rocks
as limestone .
On the other hand nitrogen, which is not very chemically active, continued to
accumulate in the atmosphere.
What about the most important gas oxygen we need to live? The first oxidized
rocks found in geological strata date back only 1.2 billion years. 600 million years
ago oxygen constituted only 1% of the atmosphere (currently 21%). Therefore,
Oxygen was only a trace gas in the air when life first appeared on the planet. That
was one of the reasons that life first evolved in the oceans. Single-celled bacterium
dwelling in the oceans did not need oxygen to live. Oxygen first appeared in the
environment when early bacteria developed the ability to split water molecules
apart using the energy of sunlight - a key part of photosynthesis. Photosynthesizing
organisms produced the oxygen that accumulated over geologic time.
These processes acting sequentially and simultaneously appear to have pro-
duced the delicate balance of 78% nitrogen (N2 ) and 21% oxygen (O2 ) we observe
today. By the way, oxygen is the third most abundant element in the universe and
makes up nearly half of the mass of the Earth’s crust, two thirds of the mass of
the human body and nine tenths of the mass of water.
The Earth cannot sustain more than ∼ 20% O2 in the atmosphere. Otherwise
spontaneous fires would occur that would deplete the oxygen.
The enrichment of oxygen in the atmosphere might be seen in context with
the methane content. Microbes who utilize photosynthesis existed on Earth half
a billion years or more before oxygen became prevalent, without substantially
affecting the composition of the atmosphere. The transition to an atmosphere
with noticeable oxygen content occurred about 2.4 billion years ago. According
to Catling et al. (2001) [55] 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 carbon 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
142 CHAPTER 5. THE EARTH’S ATMOSPHERE AND CLIMATE
environment we know today. This can also solve the faint young Sun problem (see
next chapter).
This problem was also recently discussed by Kasting, 2004 [154], where further
references can be found. He argues that the climate on Earth prior to 2.5 Gigayears
seems to have been even warmer than today, despite the fact, that the luminosity of
the Sun was 25-30% less luminous than today (see also next chapter on the faint
young Sun problem). Thus a warming of the atmosphere additional to present
day greenhouse gases was required and it is argued that this could have been
done by methanogens since Ammonia is unstable in low O2 atmospheres. CH4
photolyzes only at wavelengths shorter than 145 nm and it is relatively longlived
in the absence of O2 , O3 . It is produced by anaerobic bacteria that have evolved
early in the Earth’s history, the required flux was 500 Tg CH4 /yr. Now there is
a positive feedback: this flux should have increased once oxygenic photosynthesis
evolved because of increased production and recycling of organic matter. Even if
the CH4 flux would have been the same as at present and the CO2 at the same
level, this would have led to a warming of 30 degrees. However, siderite-coated
stream pebbles imply that also the CO2 concentration was 7 times the present
value. A rise in either atmospheric O2 or oceanic sulfate near the end of the
Proterozoic could have caused CH4 concentrations to decrease a second time and
may have triggered the “Snowball Earth” glaciations.
Omori et al., 2004 [238] discuss the role of plate tectonics and the amount of
carbon had carried into the mantle via the Archean subduction zone. They found
out that plate tectonics can be dated back as early as 3.8 Gigayears.
The role of the changing Sun and its evolution on the Earth’s atmosphere was
also discussed by Guinan and Ribas, 2002 [121].
The early evolution of the Earth and its atmosphere at the time when the
Earth was still growing by planetesimal impacts was studied by Abe and Matsui
(1986 [1]). They considered a magma ocean covering the Earth when the accretion
time was less than 5 × 107 y. Zahnle et al. (1988 [344]) discussed the evolution
of an impact generated steam atmosphere. Abe and Matsui (1988 [2]) reported
on the evolution of an impact-generated H2 O–CO2 atmosphere and formation of
a hot proto-ocean on Earth.
Chapter 6
The term Space Weather denotes variations of the Earth’s environment on short
terms. In analogy to meteorology, where the distinction between weather and
climate is made, space climate denotes long term variations of the Earth’s climate
mainly caused by solar variations.
143
144 CHAPTER 6. SPACE WEATHER AND CLIMATE
Figure 6.1: Penetration of different solar light waves resp. their induced particles in the
atmosphere
Table 6.1: Effects of Solar Radiation at different wavelengths on the Middle and Upper
Atmosphere.
Wavelength Variab. Variab. Effect Height
[nm] middle Atm. upper Atm. [km]
1-10, SXR sporadic Ion. all 70-100
10-100, XUV 2ppm 2x Ion. N2 , O, O2
100-120, EUV 6 ppm 30% Ion. NO 80-100
120-200, VUV 150 ppm 10% Diss. O2 40-130
200-240, UV 0.12% 5% Diss. O2 , O3 20-40
240-300, UV 1.0% <1% Diss. O3 20-40
• magnetospheric processes.
The principal part of the ionosphere is produced by XUV which is strongly ab-
sorbed there. Ionization and recombination occurs and this contributes to the
heating of the thermosphere. By comparing with table 6.1 we see that the energy
involved is small; dissociation of O2 is strong, above 120 km oxygen occurs as
atoms. Through vertical mixing the ratio O2 : N2 is constant near 0.1 throughout
the lower thermosphere. Oxygen atoms are produced down to 30 km and most of
them combine with O2 to form ozone. This attains a peak ratio of 10−5 near 30
km. Photons around 300 nm can reach the surface. They produce electronically
excited oxygen and surface ozone which drives a large fraction of urban pollution
chemistry.
6.1. THE ATMOSPHERE’S RESPONSE TO SOLAR IRRADIATION 145
Figure 6.2: Thermospheric temperature changes, a) low solar activity, F10.7 = 80, Ap =
0, b) high solar activity, F10.7 = 200, Ap = 80 (dashed line). F10.7 denotes the 2800 MHz
solar flux which is a measure for solar activity, the Ap index is a measure of the general
level of geomagnetic activity over the globe. Up to the top of the troposphere the two
curves are nearly identical. At a height of about 100 km the difference is already of the
order of several 10 K.
6.1.2 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 O2 in the 175 to 240 nm range. The 205 to 295 nm range is predominantly
absorbed by ozone O3 . 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.
6.2. Solar activity is measured in terms of the 10.7 cm radio flux and of the plage
area Ap .
• galactic cosmic rays: they originate from outside the heliosphere but their
input on Earth is partly controlled by solar activity.
During large flares, intense fluxes of energetic protons (10...104 MeV) penetrate
the Earth’s polar cap regions. They produce ionization between 100 and 20 km.
Such an event can last for a few hours to a few days. Large numbers of NOx
molecules are produced leading to a subsequent ozone depletion.
Relativistic electron precipitation are possible sources for ionization and odd
nitrogen production at altitudes above 80 km, thus well above the ozone layer.
• 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 10−22 Wm−2 Hz−1 . 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 6.2).
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 NO2 (which is
collectively called NOx ). The following reactions define the NOx production in
the thermosphere:
O+ + N2 → NO+ + N (6.1)
N+2 +O → +
NO + N (6.2)
N2 + e−
+
→ N+N (6.3)
−
+
NO + e → N+O (6.4)
• CO2 at 15 µm,
• O3 at 6.3 µm
• absorption of solar UV by O3 .
The CO2 band is the principal radiator. The O3 /CO2 ratio decreases upwards;
thus the heating to cooling also decreases and the temperature gradient is negative.
In the stratosphere ozone begins to be more and more attenuated, a temperature
maximum at 50 km occurs, the stratopause.
The solar UV flux is not variable at large scales thus the temperature changes
to be expected from a variation of that flux should also be small. Most effects
therefore come from ionization and a changing chemistry. The changes of the
UV flux from the sun have only a modest effect on ozone amounts because both
production (by a photolysis of O3 ) and destruction are affected in the same way.
Another effect is the penetration of solar protons or relativistic electrons into
the middle atmosphere. By that penetration considerable amounts of NOx are
produced; these enhancements of NOx increase the destruction of ozone at high
altitudes. This could explain the inverse correlation of ozone amounts with solar
activity found by Ruderman and Chamberlain (1975 [264]). Chakrabarty, 1982
[61] studied how Ozone is affected by solar proton events. During such events NO
is produced that destroys ozone.
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 O2 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
148 CHAPTER 6. SPACE WEATHER AND CLIMATE
Figure 6.3: Relation between global ozone, solar activity (measured by 10.7 cm radi-
ation) and volcanic eruptions (El Chichon, 1982, Mt. Pinatubo, 1991). Large volcanic
eruptions decrease solar transmission and the particles enhance ozone depletion. Adapted
from [Link] http : //[Link]/ozone/science/[Link]
[136], 1987 [137]) and Keating et al. (1987 [157]) between 30 and 0.2 mbar (24 to
60 km) and later by Clancy and Rusch (1989 [68]) 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.
The response of the stratospheric ozone concentration to solar activity was
examined in different studies. A study by Angell and Korshover (1976 [11]) es-
tablished a correlation of the ozone column with solar activity with a peak to
peak variation up to 10 % at 700 latitude and only 4 % at 470 . In the upper
stratosphere there is an increasing trend of ozone with increasing solar activity.
Enhanced photolysis of NO however reduces the ozone destruction during solar
activity maximum. An ozone variation of 3% seems to correspond to a 20% solar
UV variability at 180 nm (Keating, 1981 [156]). In general the stratospheric ozone
content varies with the solar activity cycle whereas the tropospheric ozone does
not. The stratospheric ozone decreases 2.72% - 3.79%, and total ozone 2,71 % -
4.36% when solar activity decreases; when solar activity increases the stratospheric
ozone increases 2.41% - 3.06% and total ozone increases 2.1% - 5.56% (see Asiati
et al., 2004 [15]).
There seems to be no correlation of polar stratospheric temperatures and solar
activity (Labitzke (1987 [179]), Labitzke, Van Loon (1988 [180]), Kerr (1988 [159]).
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.
6.2. THE FAINT YOUNG SUN 149
6.1.6 Troposphere
As we have seen above, only wavelengths > 300 nm penetrate to the troposphere
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:
• The effect found by Labitzke (1987 [179]): 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
[180]).
• Eddy (1976 [87], 1988 [88]) discussed the absence of sunspot activity during
the 17th century which is known as the Maunder minimum and an earlier
event, called the Spörer minimum. Both periods seem to coincide with pe-
riods of reduced global temperatures the more recent is called the Little Ice
Age. Eddy (1988 [88]) 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.
P ≈ kT /µ (6.5)
1 The amount of energy from the Sun received per unit area at the Earth
150 CHAPTER 6. SPACE WEATHER AND CLIMATE
In this formula L0 is the present solar luminosity and t0 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.
• Strong mass loss during the early phase (Willson et al. 1987 [335]).
Sagan and Mullen (1972) first pointed out the implications of this change of
solar luminosity for the Earth’s climate2 . Using a very simple model of the green-
house effect they showed that lower solar luminosity would have resulted in Ts
below the freezing point of water for roughly the first 2 Gyr of the Earth’s evo-
lution. However this cannot be correct. Already Sagan and Mullen pointed out
the presence of pillow lavas, mud cracks and ripple marks in 3.2 Gyr old rocks
suggesting strongly the presence of liquid water on the Earth’s surface at that
time. We also know that sedimentary rocks have been deposited about 3.8 Gyrs
ago and these must have formed in liquid water.
We should also stress here, that the faint young Sun was more active and
variable than today, especially in the short wavelengths (X, UV).
2 see also e.g. The Role of the Sun in Climate Change by Douglas V. Hoyt, Kenneth H.
The relevant albedo to use here is the Bond Albedo, which is the percentage
of the total incident solar radiation over radiation reflected back into space. The
present effective radiating temperature of the Earth is ∼ 255 K. If we combine 6.6
and 6.8 then the increase of Te 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 Ts and
T s > Te (6.9)
One possible explanation might be that the Earth’s albedo was significantly
lower in the past or that the greenhouse effect of its atmosphere was larger. How-
ever, as Sagan and Mullen pointed out a large change in the Earth’s albedo was
unlikely; any decrease in cloudiness that might result from lower surface tempera-
tures would likely be compensated by an increase in snow and ice cover. However
if the Earth’s surface was mostly water covered this argument does not work.
From climate research we know that there was no glaciation on Earth prior
to about 2.7 Gyr ago (e.g. oxygen isotopes imply warm surface temperatures
throughout the Precambrian (Kasting and Toon, 1989 [155])).
Thus:
CaSiO3 + CO2 → CaCO3 + SiO2 (6.13)
6.2. THE FAINT YOUNG SUN 153
When old sea floor is subducted and carbonate sediments are subjected to higher
temperatures and pressures CO2 is returned to the atmosphere/ocean. Then re-
action 6.13 goes in the opposite direction, calcium silicate is reformed and gaseous
CO2 is released. Much of this CO2 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 CO2
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 CO2 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,
pCO2 would decrease and the greenhouse effect would become smaller (Walker et
al. 1981 [328]). 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 CO2 release from volcanoes would create a 1-bar CO2
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.
varies approximately as
pCO0.3
2 (6.16)
for CO2 partial pressures of 2 to 20 bar and temperatures of 100 to 2000 C. These
data were derived by Lagache (1976 [183]) and Walker et al. (1981 [328]). Let
us assume that we can apply this relation to the Earth’s surface conditions. To
study the maximum effect let us further assume that removing land plants from
the system would reduce surface soil pCO2 by a factor of 40. Then the equation
for the silicate weathering process can be written as:
fw = [1 + 0.087(T − T0 ) + 0.0019(T − T0 )2 ]
= [PS /40P0 )]0.3 (6.17)
Here PS is the partial CO2 pressure in the soil and today we have PS = 40P0 and
obtain fw =1. On a vegetation free Earth PS ∼ P0 , and fw would be reduced by a
factor of 40−0.3 ∼ 1/3. The carbon cycle is only balanced when fw = 1. Therefore,
without vegetation, the atmospheric pCO2 and surface temperature would have
to increase to bring back the silicate weathering rate to its present value. We
substitute equation 6.15 into equation 6.17 and solve for P/P0 and obtain:
This shows that under the assumption that land plants pump up soil CO2 by
a factor of 40, the effect of eliminating them would be to increase the Earth’s
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 (1981) [275] showed that biota may
accelerate chemical weathering by stabilizing soil (silicate minerals stay in contact
with carbonated water), generating organic acids. This could lead to enhanced
weathering rates of up to 1000 instead of 3. Therefore, the CO2 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!
which is about 1 000 times that of the present Sun and radio flux is not constant.
Some of them also show molecular outflow and strong stellar winds. By their IR
and sub mm excess radiation it was found that about half of them are surrounded
by circumstellar disks.
Contrary to normal main sequence stars like the Sun their energy is not pro-
duced by nuclear fusion near the core but by a slow gravitational contraction. T
Tauri stars belong to the group of so called YSO (young stellar objects) of type II.
Type I YSO are very young protostellar objects at the age of just a few 100 000
years. An example is HR 4796. At a wavelength of 12 µm the object appears as a
point source, at a wavelength of 21 µm it is much larger and diffuse indicating a
circumstellar disk of dust. Such objects can be observed preferentially in the IR.
At an age of about 40 Million years our Sun became a zero age main sequence
star. That means that it reached the main sequence and nuclear fusion of H to
He started. At that time the Sun had about 70% of the total flux that is emitted
presently. But in the UV and X-rays the flux was higher by a factor of about
100 than now. This of course has important consequences for the formation of the
planets, their atmospheres etc. In its T Tauri and post T Tauri evolution the Sun’s
short wavelength emission was considerably higher than it is now. At that time
the terrestrial planets were formed already, the protoplanetary disk evaporated,
comets ejected out into the Oort cloud and the big bombardment period from
the remaining rocky planetesimals and comets began; this caused probably several
evaporations of the Earth’s oceans.
How can we find indications for this T Tauri and Post T Tauri Phase of the
Sun? Measurements of the 15 N to 14 N ratio in the atmosphere of the satellite Titan
(which is Saturn’s largest satellite) have shown that the bulk N is enhanced in the
heavier 15 N isotope by about 4.5 times relative to the Earth’s value. A 15 N/14 N
anomaly on Mars of about 1.6 times the terrestrial value has also been found.
These measurements can only be explained by the above mentioned T Tauri and
post T Tauri phase of the early Sun (Lammer et al. 2000 [186]). Atmospheric
sputtering and pick-up caused by a high solar wind particle outflow during a Post
T-Tauri phase could be responsible for the observed nitrogen anomaly.
So far we have only 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 [34]).
156 CHAPTER 6. SPACE WEATHER AND CLIMATE
The solar input can be represented by the already mentioned solar constant
S. In this context we give two interesting facts: First suppose that the Sun
has no nuclear sources then the present luminosity could be maintained by a
gravitational contraction of only 10−4 arcsec/year. Such a change would have
been imperceptible over historical times. The other fact is that, as we have seen
convection is a stochastic process and it is transport energy to the surface. The
efficiency of convection is given by
l/H (6.19)
where l is the mixing length and H the pressure scale height. Dearborn and
Newman (1978) [73] show that a variation of l/H by 0.02 causes a variation of
∆S/S by 1%. Such a variation is assumed to change the global temperature on
Earth by 2 K.
We now address to the question whether there exists also a variability of the
solar input on shorter timescales.
1369
ACRIM II
ACRIM I
ACRIM I
VIRGO
HF
HF
HF
1368
Solar Irradiance (Wm−2)
1367
1366
1365
0.1%
1364
1363
19781979 1980 198119821983 1984 198519861987 1988 1989 19901991 1992 1993 199419951996 1997 199819992000
Year
from: C. Fröhlich, Space Science Reviews, in preparation, and the VIRGO Team (Dec 03, 2000)
157
158 CHAPTER 6. SPACE WEATHER AND CLIMATE
A comparison of the power generated inside the cavities 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 Willson (1981 [336], 1984
[337]).
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. 6.6.
• 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 1 W/m2 (out of about 1367) between solar max-
imum and minimum was reported by Fröhlich (1987) [104], Willson and Hudson
(1988) [338] and Foukal and Lean (1988) [100]. 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 cor-
6.3. SOLAR VARIABILITY 159
Figure 6.6: Three-dimensional rendering of the angular distribution of the excess irra-
diance emitted at 500 nm by the active region studied at two stages of its development,
together with the magnetogram. A more uniform brightening of the facular region at the
later stage is apparent (after Vicente Domingo).
• Joule power
From such a model both influences a) radiation, b) particles from the sun to the
Earth can be taken into account. During 1975 and 2003 the contributions were
found to be: particles 36 GW, Joule 95 GW and 464 GW for the rest. Solar wind-
driven geomagnetic power provided 22% of the total global upper atmospheric
energy. An interesting trend is that with increasing activity (here by the term
activity short time scaled events are meant) the Joule power becomes more and
more important. In the top 15 power events, geomagnetic activity contributed to
2/3 of the total power budget.
160 CHAPTER 6. SPACE WEATHER AND CLIMATE
• auroral activity,
• sea surface temperatures (Gerard, 1990 [107]). The mechanism how this
could be related to solar irradiance variations works as follows:
Lewis et al. (1990 [195]) 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.
As it has been stated above the transparency of the tropical oceans is dependent
upon the amount of biogenic material, phytoplankton pigments, and degradation
162 CHAPTER 6. SPACE WEATHER AND CLIMATE
Figure 6.8: Temperatures derived from tree rings. Here the Maunder minimum is not
very prominent whereas a cold period between 1830 to 1870)
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 [96]).
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 [257]) 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 [141]) where other references can be found.
Indicators of solar activity such as 14 C concentration measurements and of
climate (e.g. glaciers) show a clear correlation over the last 7 000 years. This
was shown by Eddy in 1977. Considering a time series since 1860 the solar cycle
length shows an excellent correlation with northern hemisphere land temperatures
(Friis-Christensen and Lassen, 1991 [103]). For recent data however, these two
parameters diverge. If there exists a global climate contribution of solar irradiance
variations then there are three possible interactions or couplings between these
variations and the Earth’s climate:
6.3. SOLAR VARIABILITY 163
• Variations of the Sun’s spectral irradiance; this denotes changes in the lumi-
nosity of the Sun in a given wavelength range. As we have discussed above,
UV radiation influences atmospheric chemistry (production or destruction
of ozone, see also Haigh 1994 [122], 1996 [123]).
As we have mentioned above, the total solar irradiance varies by 0.1% and
these measurements have been made very accurately since 1978 (a review about
that was given by Fröhlich, 2000 [102]) The cycles covered by these measurements
are 21, 22 and 23. Of course from these time series it is impossible to extrapolate
to earlier time series when the Sun was more active (e.g. cycle 19) or less active.
One further problem of the time series available is that with the exception of
SOHO/VIRGO they are restricted to the UV.
Irradiance variations of the past solar cycles can be determined from the surface
distribution of the magnetic field if records of the field distribution or of proxies
are available. The following proxies can be used:
Using these data, one can reconstruct the cyclic component of the irradiance back
to the Maunder minimum.
As a large sunspot group passes across the solar surface, there is a dip in the
total solar irradiance. The variation is in the range of 0.02%.
Lockwood et al. (1999 [200]) reconstructed the aa-index of geomagnetic activity
and found that the interplanetary magnetic flux at minimum of solar activity (that
can be be reconstructed using the aa-index) has roughly doubled since 1900. This
is in good agreement with 10 Be concentration in Greenland ice (Beer, 2000 [32]).
10
Be 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 [191]) assumed that the background irradiance is proportional
to the amplitude of the solar cycle; Hoyt and Schatten (1993 [140]) propose a trend
corresponding to cycle length and Baliunas and Soon (1995 [25]) demonstrated that
the amplitudes of stellar cycles (observed in Ca II H and K) scale with the length
of the stellar cycle.
164 CHAPTER 6. SPACE WEATHER AND CLIMATE
A short overview of long term chances in solar irradiance was given by Solanki
and Fligge (2000 [293]).
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 [189])).
over a wide range and more than 80% are protons more than 12% He nuclei (α
particles), the rest electrons, gamma rays and neutrinos. Because the particles
are charged we have to take into account the interaction with magnetic fields,
mainly with the heliosphere and the Earth’s magnetic field. The particles are high
energetic, thus they produce showers of particles when they collide with particles
in our atmosphere (pions, kaons, mesons, muons).
There exist three components of cosmic rays:
• Galactic cosmic rays, GCR
• Anomalous cosmic rays, ACR
• Solar Energetic particles, SEP
Flashes of Cherenkov light from air showers have been studied for many
years.
• they are related to exotic particles which are predicted by several physical
theories concerning the origin of the universe; superstrings, exotic matter,
strongly interacting neutrinos,
As we have stated above, cosmic rays include a number of radioactive nuclei whose
numbers decrease through the radioactive decay. Measurements of these nuclei can
be used therefore (as in the C14 method) to determine how long it has been since
cosmic ray material was synthesized in the galactic magnetic field.
The origin of GCRs was discussed by Axford, 1994 [18].
One example of modern cosmic ray detectors is the Pierre Auger observatory
installed in Argentina in 2005. 1600 water tanks, each containing 3000 gallons
of water and separated by 1,5 km from each of its neighbors. The tanks are
completely dark inside. Therefore we can measure the Cherenkov radiation of
particles travelling faster than light. Slight differences in the detection times at
different tank positions help scientists determine the trajectory of the incoming
cosmic ray. The charged particles in an air shower also interact with atmospheric
nitrogen, causing it to emit ultraviolet light via a process called fluorescence. The
observatory’s second detection method uses these detectors to observe the trail of
nitrogen fluorescence.
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.
6.4. COSMIC RAYS 167
Figure 6.9: A collision between a high-energy cosmic ray particle and an atom in a
photographic emulsion (viewed through a microscope). Below: filling of a water tank of
the Auger observatory.
Because of the magnetic fields, only some of the GCR particles penetrate to
the inner part of the solar system. Thus the magnetic field of the heliosphere
works as a shield. The magnetic activity of the Sun changes however with the
solar cycle (every 11 year the Sun’s magnetic field reverses, the true cycle is thus
22 years). This causes a variation of the GCR flux. When the Sun is more active,
the magnetic field is stronger, and as a result, fewer GCR arrive in the vicinity of
the Earth. We can also say that the higher the energy the particles have, the less
they are modulated by the solar cycle.
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 14 C
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 14 C 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, during which 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:
low solar activity → weak magnetic field → more GCRs → more 14 C (6.21)
If that assumption is true, there is a mechanism, how the Earth’s climate can
be influenced by the Sun.
The correlation between cosmic rays and solar activity and temperatures on
Earth was studied by Usoskin et al., 2005 [322]. Comparison of the Sun-related
data sets with various reconstructions of terrestrial Northern Hemisphere mean
surface temperatures reveals consistently positive correlation coefficients for the
sunspot numbers and consistently negative correlation coefficients for the cosmic
rays. The significance levels reach up to 99% but vary strongly for the different
data sets.
Predictions of Galactic Cosmic Ray Intensity Deduced from that of Sunspot
Number were made by Lantos, 2005 [188].
The relationship between cosmic ray variability and enhanced geomagnetic
activity was summarized by Kudela and Storini, 2006 [178].
Orographic precipitation: this occurs when air is forced to rise because of the
physical presence of elevation (land). As such a parcel of air rises it cools due to
the adiabatic expansion at a rate of approximately 10 degrees per 1 000 m. The
rise of the parcel is stopped if saturation is reached. An example of this mechanism
is the west coast of Canada with large precipitation.
Convectional precipitation: this is associated with heating of the air at the
Earth’s surface. When there is enough heating, the air becomes lighter than the
surrounding masses, begins to rise (cf. a hot air balloon begins to rise), expands
and cools as above. When sufficient cooling takes place, saturation is reached
again forming precipitation. This mechanism is active in the interior of continents
and near the equator forming cumulus clouds and thunderstorms.
Convergence or frontal precipitation: this mechanism takes place when two
masses of air come together. One is usually moist and warm and the other is
cold and dry. The leading edge of the cold front acts as an inclined wall or front
causing the moist warm air to be lifted. Then the above described processes start
again: rise, cooling and saturation. This type of precipitation is common in the
mid latitudes.
Finally we have to mention the radiative cooling: this occurs when the Sun is
no longer supplying the ground and overlying air with energy due to insolation
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
17 − 35 W/m2 , 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 6.5).
It has been found that the Earth’s cloud cover follows the variation in GCR.
170 CHAPTER 6. SPACE WEATHER AND CLIMATE
Figure 6.10: Global radiative cooling and warming of the climate system (from IPCC
report).
• CRF/climate link therefore implies that the increased solar luminosity and
reduced CRF over the previous century should have contributed a warming
of 0.47 ± 0.19 K,
• Without any effect of cosmic rays, the increase in solar luminosity would
correspond to an increased temperature of 0.16 ± 0.04 K
Meteorite data on the galactic cosmic rays, the solar activity, and temperature
variations in the earth’s atmosphere lead to the conclusion that the solar activity
may be important factor exerting the influence upon the climate of the Earth
(see e.g. Alexeev, Ustinova, 2005[4]). Estimations on the long term cosmic ray
variation and possible climate on planets were made by Dorman, 2005 [78].
A summary of the effects is illustrated in Fig. 6.11. The data shown in that
Fig are i) Reconstructed NH temperature series from 1610-1980, updated with raw
data from 1981-1995 ii) Greenhouse gases (GHG) represented by atmospheric CO2
measurements (iii) Reconstructed solar irradiance (see Lean et al., 1995 [191]) (iv)
Weighted volcanic dust veil index (DVI) (v) Evolving multivariate correlation of
6.5. WHAT CAUSES THE GLOBAL WARMING? 173
NH series with the 3 forcings (i) (ii) and (iii). The data are from Mann et al.
(1999 [206], and further references therein). These authors conclude that while
the natural (solar and volcanic) forcings appear to be important factors governing
the natural variations of temperatures in past centuries, only human greenhouse
gas forcing alone, can statistically explain the unusual warmth of the past few
decades.
Chapter 7
175
176 CHAPTER 7. SPACE WEATHER AND RADIATION DAMAGE
• 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.
7.1. RADIATION DAMAGE ON LIVING ORGANISMS 177
DNA
damage
Mutations
Cancer
Replication errors
Replication
Persistent DNA
DNA damage
Repair Aging
Genomic instability
• Teratogenic effects: Teratogenic effects are effects from some agent, that are
seen in the offspring of the individual who received the agent. The agent
must be encountered during the gestation period.
Figure 7.2: Passage of ionizing radiation can result in direct effect on DNA leading to
single strand breaks (SSB), double strand breaks (DSB), associated base damage (BD),
or clusters of these damage types. Source: NASA
by heavy ion radiation. The most important late effect of energetic heavy 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.
• Blood forming organ (bone marrow) syndrome: damage to the cells which
divide at the most rapid pace; bone marrow, spleen and limphathic tissue.
Symptoms include internal bleeding, fatigue, bacterial infections and fever.
180 CHAPTER 7. SPACE WEATHER AND RADIATION DAMAGE
• Central nervous system syndrome (> 5000 rad): damage to cells which do
not reproduce such as nerve cells. Symptoms include loss of coordination,
confusion, coma, shock.
It seems now that death is not caused by radiation damage on the nervous system
but by internal bleeding and fluid and pressure build-up on the brain.
The genetic or heritable effect appears in the future generation of the exposed
person as a result of radiation damage to the reproductive cells.
We have seen that satellite systems are vulnerable to Space Weather through
its influence on energetic charged particle and plasma populations and that aircraft
electronics and air crew are subjected to atmospheric secondary radiation produced
by cosmic rays and solar particle events. This is discussed by Dyer (2001 [85]).
The Advanced Composition Explorer (ACE) continuously monitors the solar wind
and produces warnings by monitoring the high-energy particles that can produce
radiation damage in satellite systems (Zwickl et al., 1998 [347]).
3.9-2.5 Billion years ago the Earth was dominated by an oceanic lithosphere.
Cockell, 2000, [69] calculated that the DNA damage rates might have been ap-
proximately 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 inhabited by a low diversity UV-resistant biota.
Repair 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 inhabited the mixed layer.
In STS Shuttle/Mir mission experiment, recovery of bacterial cells from ra-
diation damage and the effects of microgravity were examined for Deinococcus
radiodurans (Kobayashi et al., 2000 [170]).
Lean (2000 [190]) discusses societal impacts of solar electromagnetic radiation.
Climate change and ozone depletion has significant economic and political impacts
on an international level.
The Yohkoh satellite was launched in 1991. Song and Cao (1999 [296]) discuss
CCD radiation damage. Evans et al. (1999 [94]) discuss charged-particle induced
radiation damage of a HPGe gamma-ray detector during spaceflight.
Figure 7.3: Correlation of the occurrence of solar proton events with solar activity cycle
(indicated by the sunspot number)
of radiation within 10 hrs after the start of the event. At 6:20 UT an optical flare
was observed on the Sun. At 13:00 UT the astronauts’ allowable 30- day radia-
tion exposure to skin and eyes was exceeded. At 14:00 the astronauts’ allowable
30-day radiation exposure for blood forming organs and yearly limit for eyes was
exceeded. The yearly limit for skin was exceeded at 15:00 UT. At 16:00 UT the
yearly limit for blood forming organs and the career limit for eyes was exceeded.
At 17:00 UT the career limit for skin was exceeded.
This event dramatically shows the need for space weather forecasting. The
correlation of solar proton events with activity cycle is evident (Fig. 7.1.6).
Heckman (1988 [130]) discussed proton event predictions.
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.
• 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;
Figure 7.4: Typical clear-sky UV indices over New Zealand and its surrounding region.
Seasonal variations are larger at low latitudes (denoted by numbers).
(e.g. wood, PVC) and loss of mechanical integrity (see also Andrady et al. 1998
[10]).
Caucasians have a higher risk of skin cancer because of the relative lack of
skin pigmentation. The worldwide incidence of malignant melanoma continues
to increase, and is strongly related to frequency of recreational exposure to the
sun and to history of sunburn. There is evidence that risk of melanoma is also
related to intermittent exposure to UV, especially in childhood, and to exposure
to sunlamps. However, the latter results are still preliminary.
7.2.6 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.
186 CHAPTER 7. SPACE WEATHER AND RADIATION DAMAGE
• blood and other body fluids start to boil and then freeze because there is no
air pressure,
• Lower Torso Assembly (LTA) - lower half of the EMU including pants, knee
and ankle joints, boots and lower waist Hard Upper Torso (HUT) - hard
fiberglass shell that supports several structures including the arms, torso,
helmet, life-support backpack and control module Arms Gloves - outer and
inner gloves Helmet.
• In-suit Drink Bag (IDB) - provides drinking water for the astronaut during
the spacewalk.
188 CHAPTER 7. SPACE WEATHER AND RADIATION DAMAGE
Radiation with energy less than 1 MeV/nm can not penetrate a space suit
of 1 mm thickness. Al shielding reduces the low boundary to 40 MeV. When
a high-energy ion strikes an atom in metal shielding it can produce secondary
radiation and there are cases where a small amount of shielding is worse than
none at all. Bremsstrahlung can be created (X-rays) by electrons as they interact
with spacecraft shielding.
Radiation damage of electronic components in space environment was studied
by Boscherine et al. 2003[44].
The radiation-induced degradation of polymeric spacecraft materials under
protective oxide coatings was studied by Lachance et al. 2001[181].
Polyethylene (Cn Hn ) is a relatively inexpensive, stable, and, with a low atomic
number, an effective shielding material that has been certified for use aboard the
ISS. Several designs for placement of slabs or walls of polyethylene have been eval-
uated for radiation exposure reduction- and it is shown that 20% or mor reduction
in dose in the crew quarters is achievable (see Shavers et al. 2004, [280]).
Magnetosphere, Ionosphere,
Space Weather
1, 2
A good introduction to the space environment was given in the books .
191
192 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
µ0 M 2 cos θ
BD,r = − (8.2)
4π r3
µ0 M sin θ
BD,θ = − (8.3)
4π r3
BD,Φ = 0 (8.4)
Here, M denotes the magnetic moment of the dipole, for the Earth ME ∼ 8 ×
1022 Am−2 . θ is the angle between the dipole moment and the radiusvector r.
From the above equations we can derive the magnitude of the magnetic induc-
tion at r:
B(r, θ) = Br2 + Bθ2 (8.5)
µo M
B(r, θ) = 1 + 3 cos2 θ (8.6)
4πr3
Note that the dipole strength decreases by r−3 . In the plane of the equator θ = 0
and:
µ0 M rE 3
Beq = 3
= B0 (8.7)
4πr r
The induction on the geomagnetic equator is B0 ∼ 31000 nT and rE = 6378 km
the radius of the Earth.
The space enclosed by the magnetosphere is not empty but filled with trapped
particles, namely ions and electrons. The magnetic forces are much stronger than
gravity.
The real shape of the boundary of the magnetosphere, the magnetopause, is
strongly modified by the solar wind. The distance of the magnetopause is
There exists also a neutral gas envelope of the Earth, the Geocorona that extends
to 4-5 rE .
At the side facing the Sun there must be an equilibrium between two pressures:
B2
f nmv 2 cos2 α = (8.8)
2µ0
A further condition is that the normal component of the magnetic field is zero:
Bn = 0 (8.9)
From 8.8 and the typical values for the solar wind particles: n = 10−7 m−3 ,
v = vSW = 300 kms−1 , and r = 10rE (subsolar point),we get:
1/6
B2
r = rE 2 (8.10)
2µ0 f nmvrmSW
1. bow shock: in this front region solar wind particles hit the magnetosphere.
The solar wind particles have Mach numbers > 1 that means they are su-
`
personic. This is valid for both the Alfven and the sound velocity:
B cp p
vA = √ vS = (8.11)
µ0 ρ cv ρ
and the corresponding Mach numbers are MS = vSW /vs ; MA = vSW /vA ∼
10.
2. The region between the bow shock and the magnetopause is called magne-
tosheath. Here the particles become thermalized- kinetic energy is converted
to thermal energy and the plasma is highly turbulent there.
3. The solar wind stretches the dipole field, compressing it on the side towards
the sun and stretching it into a long tail region. The field lines close at very
large distances (∼ 3000 RE ).
4. plasmasheet: this is a sheet of plasma in the tail region dividing the two lobes
of the Earth’s magnetic field. For both electrons and protons the particle
density is 0.5 cm−3 .
194 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
5. lobes: they are in the magnetotail have opposite direction and are separated
by the plasmasheet -otherwise they would cancel.
6. plasmasphere: a torus shaped region, surrounding the Earth. It was detected
in 1963 and has a very sharp edge at the plasmapause extending to 4-6 Earth
radii. It can be also regarded as an extension of the ionosphere. Inside the
plasmapause geomagnetic field lines rotate with the Earth. Outside the
plasmasphere, magnetic field lines are unable to corotate, the solar wind
influence is too large. The plasmasphere is mainly composed of hydrogen.
7. Van Allen radiation belts: in 1958 Van Allen discovered the radiation belts;
like the plasmasphere they are toroidally shaped. The inner radiation belt
extends from 400 to 12000 km above the Earth, the outer belt from 12000
to 60000 km.
In order to understand the dynamics of the current system, we recapitulate the
motions of charged particles in a magnetic field:
1. spiral motion: circling about magnetic field lines; Charged particles cannot
easily move across magnetic field lines but are forced to spiral around them.
Electrons encircle the field line in one direction, ions in the other direction.
2. Bounce motion: the particles move along the field lines from pole to pole.
Near the poles they become reflected (since the magnetic field line density
is large).
3. drift motion: Curvature of the magnetic field lines and the non-uniform
strength of the magnetic field force particles to drift around the earth, ions
in one direction, electrons in the other. For the Earth as seen from Europe:
Ions go west, electrons east.
In a magnetic field particles are being transported and this causes currents. Due
to the currents magnetic fields are generated. In a magnetosphere there are three
distinct current systems:
1. Chapman Ferraro currents: they enclose and confine the magnetosphere and
are found in the vicinity of the magnetopause.
2. cross tail currents: pass through the center of the magnetotail causing the
current sheet.
3. Field aligned currents: transient currents, short circuit through a planet’s
ionosphere and cause aurorae.
How is the magnetosphere influenced by the solar wind?
• The interaction of enhanced solar wind pressure on the dayside cause a strong
reduction of the magnetopause even below the geostationary orbit (6.6 rE ).
The observed variations of the distance of the dayside magnetopause are in
the range 4.5 to 20 rE .
• The magnetic moment of the interplanetary magnetic field (magnitude and
orientation) determines the size and extension of the magnetosphere.
8.1. GENERAL PROPERTIES 195
Figure 8.1: The Earth’s magnetosphere (above) and the plasmasphere (below). NASA
196 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
Figure 8.2: The inner and outer Van Allen radiation belt. NASA
Lobe
s .
B toward Earth
Earth
Ring current B away from Earth
x
Lobe
s
Figure 8.3: The global structure of the Earth’s magnetic field. In the two lobes the
field is opposite and the lobes are separated by a plasmasheet. For the southern half of
the magnetosphere the current is clockwise, for the northern part it is counterclockwise.
In the middle both systems add to form a neutral sheet. The right drawing is a cross
section of the left at a distance of 20 rE .
8.1. GENERAL PROPERTIES 197
Ionograms are recorded tracings of reflected high frequency radio pulses generated
by an ionosonde. There exist relationships between the sounding frequency and
the ionization densities which can reflect it. As the sounder sweeps from lower
to higher frequencies, the signal rises above the noise of commercial radio sources
and records the return signal reflected from the different layers of the ionosphere.
The top of the ionosphere is at about 1000 km, however there exists no definite
boundary between plasma in the ionosphere and the outer reaches of the Earth’s
magnetic field. In the E region the most important ions are O+ +
2 , NO , in the F
+
region it is O . In the F2 layer (at about 400 km) the electron concentration
reaches its highest values which is important for the telecommunication systems.
At high latitudes there is another source of ionization of the ionosphere– the aurora
(see next chapter).
The so called transition height starts at the height of the maximum density of
the F2 layer of the Ionosphere and extends upward with decreasing density to a
transition height where O+ ions become less numerous than H+ and He+ . The
transition height depends on day and night:
• daytime: ∼ 800 km
Above the transition height, the weak ionization has little influence on radio sig-
nals.
Some ionospheric parameters are listed in Table 8.1 where the values Ne and
Te denote electron density and electron temperature. For comparison, the values
of the solar corona are also given.
4 e.g. see Ionospheres : Physics, Plasma Physics, and Chemistry by Robert W. Schunk,
Andrew F. Nagy, Alexander J. Dessler (Editor), John T. Houghton (Editor), Michael J. Rycroft
(Editor), Cambridge Univ. Press, 2004
198 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
At low latitudes the largest electron densities are found in peaks on either side
of the magnetic equator, which is called the equatorial anomaly. Normally one
would expect that the peak concentration will occur at the equator because of the
maximum in solar ionizing radiation. This peculiarity can be explained by the
special geometry of the magnetic field and the presence of electric fields. The elec-
tric fields transport plasma and are caused by a polarizing effect of thermospheric
winds.
The ionosphere varies because of two reasons:
receivers must correct the delay of the GPS signal as it propagates through the
ionosphere to the GPS satellite (at 22 000km altitude).
The maximum usable frequency depends on the angle of the wave relative to
the horizon.
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.
1000 km
600 km
Altitude
F2
F Region
F1
150 km
E Region
90 km
D Region
10 4 10 5 10 6
Electron density
Figure 8.5: Original drawing of Chapman and Ferraro showing the interaction of plasma
from the Sun and the Earth’s magnetic field.
Generally the solar wind arriving at the Earth’s magnetopause has the following
pressure components:
Dynamic → ρv 2 (8.12)
Static → nkT (8.13)
Magnetic → B 2 /2µ0 (8.14)
8.2. SOLAR ACTIVITY AND MAGNETOSPHERE 201
Figure 8.6: An interplanetary coronal mass ejections interacts with the Earth magnetic
field (from [Link]
The pressure applied by the solar wind to the magnetopause varies with the angle
of the normal to the solar wind flow. The pressure is dominated by the dynamic
pressure. At the magnetopause the dynamic pressure is zero and the static pressure
dominates. Inside the magnetosphere the pressure is dominated by the magnetic
pressure.
The so called standoff distance i.e. the distance of the magnetopause is given
by:
2
Lmp = 107.4(nSW vSW )−1/6 (8.15)
The interaction between the interplanetary magnetic field and the Earth’s mag-
netic field depends on the orientation of the former with respect to the Earth’s
field. This was studied first by Dungey (1961 [82]) and is called Dungey’s model.
The pressure of the solar wind rises and falls. The reacting of the magnetopause
is a shrinking or expansion. When the boundary is hit by a fast flow from a CME,
the shrinking can go beyond the geosynchronous orbit of satellites (at 6.6 RE ). As
it is seen from the drawing (Fig. 8.7), a southward oriented IMF is recognized as
the most important factor promoting storms and substorms in the magnetosphere
(Fairfield and Cahill, 1966 [95]).
When the interplanetary magnetic field is oriented southward, then a flow of
plasma is predicted to the dayside of the magnetosphere after reconnection in the
tail.
202 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
Figure 8.7: Interaction of a southward oriented IMF with the magnetosphere. Possible
reconnection points are denoted by N1 , N2 , N3 . Also the formation of a disconnected
plasmoid is indicated. From “A Brief History of Magnetospheric Physics during the
Space Age” by D.P. Stern
The Sun heats the Earth’s atmosphere. Also the degree of ionization in the
ionosphere increases at the dayside and this causes convection in the ionosphere.
By this convection charged particles are transported into the magnetosphere and
by dynamo action ionospheric electric currents above the equator up to mid lat-
itudes are generated. These currents produce a magnetic field which moves with
the subsolar point. So there is a 12 h variation for a given observing site in the
measurements of the field strength.
The Sun emits particles and the solar wind compresses the magnetosphere as
it has been mentioned before. High speed particles further compress the magne-
tosphere, and a magnetic storm begins with a SSC (storm sudden commencement).
The number of charged particles trapped within regions of the magnetosphere (ra-
diation belts) is increased. These particles drift around the Earth creating a ring
current that produces a depression of the horizontal magnetic field, seen at lower
latitudes around the world as a magnetic storm. This is followed by the recovery
phase, lasting one day or more, during which the ring current subsides and the
magnetic field returns to normal. Charged particles are guided down the field
lines into the upper atmosphere. This produces auroral electrojets (large horizon-
tal currents that flow in the D and E regions of the auroral ionosphere) which are
intense east-west currents. Associated with these currents are intense magnetic
fields causing magnetic disturbances observed there.
8.2. SOLAR ACTIVITY AND MAGNETOSPHERE 203
The interaction of plasma of the solar wind with the Earth’s magnetosphere
causes currents as shown in Fig. 8.8:
• flow eastwards down the morning side around the polar regions
• flow spacewards in the evening side.
Interplanetary field lines are swept back around the Earth’s magnetic field by the
solar wind. There is an electric field according to
This equation follows from Ohm’s law (j ∼ σ(E + v × B) and in the case of a
large conductivity σ >> 1 the term j/σ → 0. This field is from dawn to dusk,
there is a field aligned current and particles move from dawn to dusk. Because
this circulation is analogous to thermal convection cells, this phenomenon is also
called convection electric field.
Finally, one also has to take into account the corotational electric field. This
is cause by the rotation of a planet’s magnetic field which induces an electric field
in the radial direction. The magnetic field moves at
v = ωrot × r (8.17)
and an electric field is induced by Ohm’s law (again we consider large conductivity):
v×B (ωrot × r) × B0
Ecor = − =− (8.18)
c cr3
The motion of charged particles in the magnetosphere is thus caused by a
• drift: due to gradient in the field strength, field curvature
• acceleration due to electric fields along the field lines, field aligned currents
Where do the particles come from? Interplanetary particles can enter the
magnetosphere via different processes:
1. spiral down into the polar cusp- there are open magnetospheric field lines
there. Atmospheric ionization is enhanced there during enhanced solar ac-
tivity → aurora.
2. reconnection is an important process. It occurs when the interplanetary field
has a component antiparallel to the planetary field. Reconnection leads to
neutral points and solar wind particles can enter there. The locations of
reconnection are the day side magnetopause and the magnetotail.
3. Kelvin-Helmholtz instability . The fast solar wind flows past the magn-
teosphere inducing ripples in the magnetospheric boundary. These ripples
induce a filed perpendicular to the solar wind flow and thus particles diffuse
into the planetary magnetosphere.
Particles are lost because of different processes, such as losses due to the mirror
points. For particles with a certain pitch angle, the mirror points lie within the
atmosphere and the particle gets lost. Another process is charge exchange of
magnetospheric ions.
8.2.3 Aurora
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
8.2. SOLAR ACTIVITY AND MAGNETOSPHERE 205
Solar
wind
Sun
Earth
Figure 8.8: Birkeland currents. The currents flow downwards on the morning side and
spacewards on the evening side.
patches blink on and off every 10 s or so. Most of aurorae are greenish yellow and
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 (see 8.7). When the solar wind
adds sufficient magnetic energy to the magnetosphere, the field lines there over-
stretch 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 substorm.
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
206 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
Table 8.5: Extension of the auroral zone. The first values given is the magnetic latitude
(Lat), the second the Kp index.
Lat Kp Lat Kp Lat Kp Lat Kp Lat Kp
66.5 0 64.5 1 62.4 2 60.4 3 58.3 4
56.3 5 54.2 6 42.2 7 50.1 8 48.1 9
station. The globally averaged Kp 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
8.4) around the earth that would otherwise rarely experience the polar lights. In
Table 8.5 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
deviation from quiet day values of the horizontal magnetic field (h) around the
auroral oval. Defined and developed by Davis and Sugiura in 1966, 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 mag-
8.2. SOLAR ACTIVITY AND MAGNETOSPHERE 207
netic field and the Earth’s magnetosphere. For these varied topics, AE possesses
advantages over other geomagnetic indices or at least shares their advantageous
properties.
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.
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 8.7
Another index devised to express geomagnetic activity on the basis of the Cp
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 8.8
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.
210 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
Figure 8.9: Different layers in the ionosphere. Reflection of radio waves occur in the
E-layer at 110 km and also the F layers (170 km, 250 km) reflect waves. In the D-layer
an absorption occurs. Within the auroral oval the nighttime E layer plasma densities can
be much higher.
• The cutoff frequency is the frequency below which a radio waves fails to
penetrate the ionosphere.
TV and FM radio stations (on VHF) are affected little by solar activity. HF
ground to air, ship to shore, amateur radio etc. are affected strongly. Also the
Faraday rotation of the plane of polarization has to be taken into account (for
satellite which employ linear polarization up to 1 GHz).
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.
8.2. SOLAR ACTIVITY AND MAGNETOSPHERE 213
• Increase of harmonics,
• voltage drops,
When flowing from the pipeline into the soil, GIC may increase corrosion of
the pipeline, and the voltages associated with GIC disturb the cathodic protection
system and standard control surveys of the pipeline.
214 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
On March 13, 1989, the most famous GIC failure occurred in the Canadian
Hydro-Quebec system during a great magnetic storm. The system suffered from
a nine-hour black-out.
A theoretical calculation of GIC in a given network (power grid, pipeline etc.)
can be divided into two steps:
• 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 [45]). A prediction of Geomag-
netically Induced Currents in Power Transmission Systems was given by Pirjola et
al. (2000 [246]).
A short description of the vulnerable Swedish power system (because being
close to the auroral oval) and pipeline system together with a historical description
of the effects that occurred at times of geomagnetically induced currents (GICs),
up to the Halloween events in 2003 and event in November 2004 was done by
Lundstedt, 2006 [203].
On 30 October 2003 50 000 customers in Southern Sweden had no electricity
due to a power failure caused by a GIC (see Pulkkinen et al., 2005 [250]). Research
on historical geomagnetic storms can help to create a good data base for intense
and super-intense magnetic storms. For the event on March 13, 1989 the Dst=-640
nT. Lakhina et al. 2005 [184] claimed to have found evidence for a superstorm
that occurred on Sep 1-2 1859 with a Dst=-1760.
• HF Communications
– Increased absorption
– Depressed MUF
– Increases LUF
– Increases fading and flutter
• Surveillance Systems
Penetration
Absorption Scattering
Reflection
Earth
• Satellite Systems
– Faraday rotation
– Scintillation
– Loss of phase lock
– Radio Frequency Interferences (RFI)
• Navigation Systems
– Position errors
solar EUV heating the GITM calculates an altitude dependent heating efficiency.
The heating efficiency starts with a value of 0.25 at a height of 100 km, reaches
a maximum of 0.6 at a height of 150 km and declines down to 0.25 at a height of
250 km.
8.3 Satellites
For a general introduction to space technology several textbooks are available6 .
• 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.
6 Gatland K, “Space Technology”, Salamander Books (London, 1981)
8.3. SATELLITES 217
Figure 8.12: Solar power systems installed in the areas defined by the dark disks could
provide a little more than the world’s current total primary energy demand (assum-
ing a conversion efficiency of 8%). That is, all energy currently consumed, including
heat, electricity, fossil fuels, etc., would be produced in the form of electricity by so-
lar cells. The colors in the map show the local solar irradiance averaged over three
years from 1991 to 1993 (24 hours a day) taking into account the cloud coverage avail-
able from weather satellites. The average electric output would be 18 TW. After:
http : //[Link]/ml/solar land area/
by Markvart et al. (1987 [210]). Defect interactions in silicon solar cells were
analyzed by Markvart et al. (1989 [209]). A study of radiation-induced defects
in silicon solar cells showing improved radiation resistance was made by Peters et
al. (1992 [245]). General information about solar cells can be found in Tada et al.
(1982 [307]).
A review on radiation damage in solar cells was given by Markvart (1990 [207]).
An analytical study has been carried out on an impact feature within a solar
cell from the Hubble Space telescope Solar array. The feature was investigated
optically, and the damage was seen as the result of a partially penetrating impact
and therefore some impact particles must have been responsible for that. The
residue in the impact was found to contain elements such as Fe, Ti, K, Ca, Si,
Mg and Na. The elements Mg, Fe and Ti are usually foreign to a solar cell and
this suggests that the impact residue may be of natural or man made origin.
Subsequent detailed analysis showed Fe and Mg in concentrations of about 10%
and Ti in only limited amounts. That implies that the residue is of natural origin.
A more detailed description can be found in Graham et al. (1997 [120])
• 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;
8.3. SATELLITES 219
– they always need to be pointed at the Sun (not being 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.
• Radioisotope thermoelectric generators: They are also reliable but tend to
be expensive to build and of course there is a risk that radioactive material
is set into the environment during a launch failure.
A radioisotope thermoelectric generator, or RTG, uses the fact that radioac-
tive materials (such as plutonium) generate heat as they decay. The heat
is converted into electricity by an array of thermocouples which then power
the spacecraft.
A thermocouple is a device which converts thermal energy directly into elec-
trical energy. Basically, it is made of two kinds of metal that can both
conduct electricity. They are connected to each other in a closed loop. If
the two metals are at different temperatures, an electric potential will exist
between them. When an electric potential occurs, electrons will start to flow,
making electric current.
Another process which belongs to this group of energy generation is nuclear
fission where unstable radioactive materials are split into smaller parts. Very
large amounts of heat are generated but the whole process is more complex
and not as reliable as using the heat produced by radioactive decay. An
RTG is steadier.
Plutonium is a very toxic heavy metal. If it is powdered and inhaled, it is a
cancer causing agent. It is sealed inside a hard, radiation proof shell. The
shell is designed to survive all conceivable accidents, so even in the unlikely
event of a launch failure, none of the radioactive particles should escape.
• Fuel cells: they are similar like batteries but they have a longer lifespan and
can be refuelled. They are already in use in the Space Shuttle. However
they run hot (400-8000 C) and the waste heat is often hard to manage.
When atoms of the two gases oxygen and hydrogen are put next to another,
they spontaneously combine to form water. This results in the release of a
220 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
lot of energy. In a fuel cell the H and O are separated by a membrane. The
refuelling means just to provide more H and O and the waste is pure water.
With an external source such as a solar panel, one can split the waste water
back into its component parts and use it again as fuel. Fuel cells were first
used by the Apollo missions since they last longer than traditional batteries
and didn’t have expensive radioactive parts.
It is extremely important to control the heat on and around a space ship. The
operating temperature is usually given between two numbers like −10◦ C to 60◦ C.
The parts of the spacecraft have been tested and will work if the temperature in the
spacecraft is between these two numbers. Why does a spacecraft have an operating
temperature? For example the rocket thruster can use hydrazine as rocket fuel.
Therefore the tanks, plumbing and pumps must be kept at a certain temperature:
Hydrazine freezes at 2◦ C and boils at 113◦ C. Most electronic components will
work only within a narrow range of temperatures, usually −50◦ C to +150◦ C,
components will stop working and make the spacecraft useless if the spacecraft
temperatures become too extreme.
Heat tends to expand material parts and the opposite happens when a part
is cooled. This problem occurs when one part of the spacecraft is pointed at the
Sun and the other one is pointed at empty space. The Sun then heats up only one
part and this uneven heating causes the spacecraft to be warped or even break or
instruments can be distorted. Another source of heating is caused by electronic
components. Heat also makes the electrical system less efficient. Electricity is
caused by the flow of electrons and the resistance grows with temperature.
Heat sources can be external (from outside the spacecraft) or internal (from
inside the spacecraft). External heat sources include:
• the Sun,
USAF uses empirically defined values to issue warnings for satellite operators.
Damaging conditions are assumed when the daily electron flux (which is given by
the number of high energy electrons (> 2MeV) per cm2 per sterad per day meets
either of the following conditions7
Such conditions often occur about 2 days after the onset of a large geomagnetic
storm.
How can we determine the probability that surface charging may occur. This
can be done by the K-index which, as we have shown in the previous chapter, is a
measure for geomagnetic storms. The values of K (3 hourly measure) range from
0-9.
• K=0: quiet
Whereas surface charging usually does not cause big problems, particles with
≥1 MeV cause Deep Dielectric Charging. When there occurs a high-speed solar
wind stream these particles are concentrated in the Van Allen belts. High energy
electrons penetrate the spacecraft’s outer surface; they penetrate the dielectric
materials such as circuit boards and the insulation in coaxial cables. This gives
rise to intense electric fields; as soon as they exceed the breakdown potential of
the material they produce sudden discharges (similar to a stroke). This discharge
damages the system: components may start to burn, semiconductors may be de-
stroyed. These dielectric charging can be avoided by a special construction of
the relevant parts however this leads to additional weight and complexity of the
system.
Again, high fluxes of these electrons vary with the 11 year solar cycle and are
most prevalent late in the cycle and at solar minimum. The GOES GEO spacecraft
measures electron fluxes in the range of 0.6 - 2 MeV (see Fig. 8.13).
Figure 8.13: Flux measurements by the GOES satellite in different energy channels.
http : //[Link]/rt plots/elec [Link]
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 Earth’s surface
magnetic field is weakest there. Particles drifting around the Earth travel much
closer to the Earth than at other latitudes and longitudes.
This higher particle concentration causes a maximum of the distribution of
errors 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. 8.14). 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
8.3. SATELLITES 223
Figure 8.14: Single event upsets; spatial distribution of errors from the UoSAT-3 space-
craft in polar orbit; please note the South Atlantic Anomaly. Adapted from C. Dyer and
D. Rodgers, 1998, Space Dep. DERA
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 (RJ , measured
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
224 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
• geomagnetic index Ap .
It is extremely difficult to predict exactly when a satellite will re-enter the at-
mosphere8 . The reason for that is that the space environment is not exactly
predictable 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. 8.15 a rough estimate of the lifetime of a satellite with effective mass
8 see also: Tobiska W K, R D Culp and C A Barth, “Predicted Solar Cycle Twenty-Two 10.7
cm Flux and Satellite Orbit Decay”, Journal of the Astronautical Sciences, 1987, pp419-433, vol
35
8.3. SATELLITES 225
to cross section ratio 100 kg /m2 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%.
• particles: 36 GW
• Joule heating: 95 GW
These values were obtained for the period 1975 to 2003 (see Knipp et al, 2004
[169]). Thus solar wind driven geomagnetic power is about 1/5 of the total budget
but becomes dominant when solar activity rises.
The Korean satellite KOMPSAT1 has an altitude of 685 km and was observed
for its daily drag acceleration which was found to correlated strongly with geo-
magnetic disturbances9 . The results are shown in Fig. 8.14.
Figure 8.16: Solar wind speed and density from the ACE spacecraft
(a), Ap and F10.7 indices (b), provisional Dst index (c), daily aver-
aged Dst (d) and daily KOMPSAT1 acceleration (e). Adapted from
[Link]
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:
1
D = ρv 2 ACd (8.26)
2
D... drag force, ρ atmospheric density, v... speed of the satellite, A... cross
sectional area perpendicular to the direction of motion, Cd ... drag coefficients.
The latter can vary; at altitudes at which satellites orbit Cd ∼ 2. We introduce
the effective cross sectional area Ae = ACd .
For a circular orbit we have the following relation:
P 2 GMe = 4π 2 a3 (8.27)
G... gravitational constant, Me ... mass of the Earth. The reduction in the period
due to atmospheric drag is given by:
dP Ae
= −3πaρ (8.28)
dt m
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 100 kg /cm2 . 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 = 400 km 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.
Satellites are perturbed to the first order by the Earth’s oblateness and at-
mospheric drag is important since it acts as energy dissipation. Space weather
influences on
8.3. SATELLITES 229
Figure 8.17: Forces acting on a satellite in a low circular orbit. From IPS Radio and
Space Services, Satellite Orbital Decay Calculations.
• atmospheric composition
These effects result to varying densities in the atmosphere causing drag effects
discussed above. The drivers are the varying solar UV radiation and coronal mass
ejections that cause geomagnetic storms. A recent review was given by Doornbos
and Klinkrad (2006, [77]).
space station had to move to the aft end of the service module where a higher
protection from enhanced radiation was expected.
Glass spherules : These may be formed in two ways. A meteorite impact melts
material, the melt is ejected from the crater and small globs of the melt
solidify. They are spherical ranging from mm to cm. The other origin is
12 There were even concerns before the landing of Apollo 11 that the regolith would not be
supportive enough to cope with the weight of the lunar module and that the module would begin
to sink beneath the surface.
8.4. SPACE WEATHER ON MOON AND MARS 231
pyroclastic glass, molten rock cools and solidifies above the Moon’s surface,
leading to glassy spherules.
The regolith contains much information on present and past solar corpuscular
radiation. Especially solar wind Ar seems to be well retained in the minerals that
were investigated by Wieler et al., 1980 [334] and Wieler et al. 1995 [333].
Therefore, the effects of “space weathering” (see also Keller et al., 1999 [158]) on
the formation of lunar soils provides the ground-truth requisite to understanding
regolith development on all atmosphereless bodies in the solar system, e.g. like
asteroids or satellites like Phobos or Deimos. Exposure to the solar wind and
solar cosmic rays for long periods of time may cause substantial alterations in the
geochemistry of the target material. Experimental results already made in 1967
(see Zeller and Ronca, 1967 [345]) indicate that hydration of oxygen-rich materials
can be expected in any surfaces which are exposed to solar protons.
8.4.3 Mars
Because Mars will be a target of future manned space missions we briefly discuss
space weather influences there.
Mars is a dry like a desert, cold as the Earth’s Antarctic and possibly lifeless
and future human colonists will be exposed to extreme sets of weather conditions.
The Earth is protected by the magnetosphere. Mars does not possess a global
magnetic field to shield the surface from SEPs and cosmic rays. The solar wind
gradually eroded the martian atmosphere and at present the surface pressure is
1-9 millibars, depending on altitude; the average is 7 mb.
In astrobiology it is believed that protection by an atmosphere and magnetic
field are essential factors for life on a [Link] that mean that because Mars
has no global magnetic field and a very thin atmosphere the planet is lifeless?
232 CHAPTER 8. MAGNETOSPHERE, IONOSPHERE, SPACE WEATHER
It can be 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: redundancy. 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 nat-
ural protection against radiation. In March 2004, small quantities of methane
(about 10 parts in a thousand million) were detected in the Martian atmosphere
by researchers operating the Planetary Fourier Spectrometer (FPS) experiment
on the orbiting Mars Express spacecraft and also by astronomers using the Keck
Telescopes and the Canada-France-Hawaii Telescope on Earth (see e.g. Krasnopol-
sky, 2006 [172]). Normally, under conditions in the Martian atmosphere Methane
should be destroyed by UV solar radiation within a few hundred years. Its presence
can be explained either by volcanic outgassing or by the action of methanogens
(bacteria that produce methane).
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
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.
Another problem on Mars which is not directly related to space weather is
martian dust. Mars is red because its surface is largely composed of iron oxide
(rust) and oxides of other minerals. Some scientists suspect that the dusty soil
on Mars may be such a strong oxidizer that it burns any organic compound such
as plastics, rubber or human skin as viciously as undiluted lye or laundry bleach.
Data from the Pathfinder mission showed that Martian dust may also contain
trace amounts of toxic metals, including arsenic and hexavalent chromium.
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 (see Fig 8.18). 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 distrib-
ution 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.
8.4. SPACE WEATHER ON MOON AND MARS 233
Figure 8.18: Artistic illustration of Earth magnetic field and Mars magnetic field (not
influenced by the solar wind). Credit: NASA
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. The strongest magnetic anomalies were found in Terra Sirenum which
might therefore be the most suitable landing site for manned missions (Alves and
Baptista, 2004[7]).
The Mars Odysssey spacecraft (2001) has onboard the Martian radiation en-
vironment experiment (MARIE) which measures the background radiation due
to galactic cosmic rays and solar protons. Good agreement between models and
measurements were found (Atwell et al. 2004,[16]).
Space Radiation Hazards on Human Missions to the Moon and Mars are de-
scribed in Townsend, 2004 [314].
Crewmembers of the ISS will be exposed to ionizing radiation (the inclination of
ISS to the Earth orbit is 51 degrees, it is in low Earth orbit, LEO). The concept is as
low as reasonably achievable (ALARA) radiation exposure and cheap Poylethylene,
Cn Hn with low atomic number can be used for shielding ( Shavers et al., 2004
[280]). Models like HZETRN allow computer simulations of radiation hazards (
Wilson et al., 2004 [339]).
The complex reactions that are caused by solar radiation and particles in dif-
ferent planetary atmospheres are discussed (mainly for the UV radiation) in the
book of Vazquez and Hanslmeier, 2006 [325].
Chapter 9
In this chapter we will give an overview over existing real-time space weather
centers and the data available there. First we will discuss a very useful classification
scheme.
• Geomagnetic storms
• Radio blackouts
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
235
236 CHAPTER 9. REAL-TIME SPACE WEATHER AND FORECASTS
Figure 9.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.
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 Kp value is about 5.
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 Kp value is about 6 and on the average one can
expect 600 events per cycle (corresponding to about 360 days per cycle).
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 Kp value
9.1. NOAA SPACE WEATHER SCALES 237
is about 7 and on the average one can expect 200 events per cycle (corresponding
to 130 days per cycle).
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 low frequency navigation can be disrupted.
Aurora has been seen down to 450 geomagnetic latitude. The Kp index is at 8
and on the average one has to count with 100 events per cycle (corresponding to
60 days per cycle).
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 amperes, HF radio
propagation may be impossible in many areas for one to two days, satellite navi-
gation may be degraded for days, low-frequency radio navigation can be blocked
for hours. The aurora has been seen down to 400 geomagnetic latitude (Italy,
southern Texas). At this level Kp = 9 and on the average one has to expect 4
events per cycle (corresponding to about 4 days per cycle).
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.
S2
moderate; there are no biological influences; for satellite operations it is important
to know that infrequent single-event upsets are possible.
238 CHAPTER 9. REAL-TIME SPACE WEATHER AND FORECASTS
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 can expect about 25 events per cycle.
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 103 and about 10 such
events per cycle may occur.
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 104 . There are about 3 such events per cycle.
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 105 ; fortunately, these events occur on a rate fewer than 1
per cycle.
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 flare classification is M1 and (10−5 ).
On the average 2000 such perturbances per cycle occur (on 950 days per cycle).
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
minutes is likely. The classification of the relevant solar event goes M5 and the
flux to 5 × 10−5 .
On the average one has 350 events per cycle (300 days per cycle).
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 classification is X1, the flux 10−4 and 175 events per cycle (140 days per
cycle) are probable.
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 classification is X10, the flux 10−3 and 8 events per cycle (8 days per cycle)
can be expected.
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.
The classification is X20, the flux 2 × 10−3 and there are less than 1 events per
cycle.
240 CHAPTER 9. REAL-TIME SPACE WEATHER AND FORECASTS
9.1.4 Summary
The classification scheme given above enables very easily to estimate the effect of
geomagnetic storms and solar radiation storms on satellites and telecommunica-
tion systems. This is also extremely important for manned space mission (ISS,
international space station). On the other hand, the solar activity has declined
after heaving reached its maximum in 2000, the new cycle will start in 2006. 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.
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).
:Product: Geophysical Alert Message [Link] :Issued: 2006 Aug 23 0911 UTC
Prepared by the US Dept. of Commerce, NOAA, Space Environment Center
Geophysical Alert Message
Solar-terrestrial indices for 22 August follow. Solar flux 81 and mid-latitude A-index
21. The mid-latitude K-index at 0900 UTC on 23 August was 3 (31 nT).
No space weather storms were observed for the past 24 hours.
No space weather storms are expected for the next 24 hours.
Furthermore, you can get a Report of Solar and Geomagnetic Activity –issued
daily at 2200Z. This summarizes solar and geomagnetic conditions for the past 24
hours and provides forecasts for the next three days and moreover contains that
day’s Solar Flux value and A-index, plus forecasts for the next three.
Another valuable information is the D-region Absorption Prediction –updated
once a minute; this color-coded world map graphic depicts the current position of
the sun and how a solar flare might have caused an HF fade.
In Today’s Space Weather more information including full-disk H-Alpha solar
images etc can be found.
The radio auroral zones are typically displaced equatorward from the optical
auroral zones (or the regions where visible auroral activity can be seen with the
eye).
A world map showing the highest vertically directed frequency being returned
from the ionosphere are the critical F2-layer maps.
There are also online prediction tools to predict HF propagation paths. Also
frequencies and times for HF communication paths between two user-entered
points anywhere in the world is provided by GRAFEX. On a world map, you
click your location, enter up to ten operating frequencies, the desired UTC hour,
and then use the mouse to drag out a target region. Now you choose either real-
time or forecast ionospheric conditions and receive a color-coded coverage map for
all frequencies plus individual maps for each chosen frequency.
During the declining phase of solar cycle 23 four epochs of extremely high flare
activity occurred. This was not observed in the previous cycles 21 and 22 (only 1
event) but similar to cycle 20 (Bai, 2006 [24]).
The Solar Mass Ejection Imager (SMEI) has been tracking coronal mass ejec-
tions (CMEs) from the Sun to the Earth and beyond since February 2003. A
prediction of the arrival time at ACE within 2 hours of its actual arrival for three
events, and within 10 hours for eight events was possible. Of these eight events,
seven were detected by SMEI more than 1 day before the transient’s arrival at the
Earth (Howard et al. 2006 [139]).
The solar cycle 24, that started in 2006 could have a 30-50% higher peak than
cycle 23. This prediction is based on a flux transport dynamo model that gave
good predictions for the peaks of cycles 16-23 (Dikpati et al., 2006 [76]). This is in
contradiction with the prediction given by Schatten, 2005 [272] using the polar field
precursor method of solar activity forecasting. In that paper, the peak amplitude
of the next solar cycle 24 is estimated at 124 ± 30 in terms of smoothed F10.7
Radio Flux and 80 ± 30 in terms of smoothed international or Zurich Sunspot
number. That would mean ‘fair space weather’.
Burov, 2006 [51] wrote a paper on the possibility of getting economically sound
forecasts of rare space weather events where he mentioned the need for precise
forecasting of major events.
Space weather prediction by cosmic rays was investigated by Mavromichalaki et
al., 2006 [213] using real-time data from a neutron monitor network. The system
collects data in real-time mode from about 15 real-time cosmic ray stations by
using the internet. The main server in Athens station collects 5-min and hourly
cosmic ray data and produces forecasting.
Medium range (1– 3 days) geomagnetic forecasting using data of slow/high
speed solar wind stream interfaces were described by Gleisner and Watermann,
2006 [114].
An approach to space weather via analysis, prediction, modeling, and classifica-
tion supported by Genetic Programming was made by Jorgensen and Karimabadi
(2005 [150]). Genetic programming can search spaces of algorithms or mathemat-
ical functions. A population of candidate solutions is allowed to evolve. Then a
natural selection and survival of the fittest starts.
Chapter 10
Asteroids, Comets,
Meteroites
10.1 Asteroids
1 According to the new definition of the international Astr. Union adopted in 2006, Ceres and
245
246 CHAPTER 10. ASTEROIDS, COMETS, METEROITES
object. However, this also requires that the asteroid’s albedo be known as well.
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 10.2 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 colliders3 .
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.
0 The likelihood of a collision is zero, or well below the chance that a random
object of the same size will strike the Earth within the next few decades. This
designation also applies to any small object that, in the event of a collision, is
unlikely to reach the Earth’s surface intact.
1 The chance of collision is extremely unlikely, about the same as a random object
of the same size striking the Earth within the next few decades.
10.2.3 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 (Q) and their semi-major axes (a).
Possible NEO missions that require spacecraft with the capability to ren-
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 [279]) and
10.2. IMPACTS BY ASTEROIDS 249
Powell et al. (1997 [248]). NEOs as near Earth resources for mining are discussed
by Gertsch et al. (1997 [108]).
Asteroid diameter
3m 10 30 100 300 1km 10km
Month
6
10 yrs
10
8
yrs K/T Impact
4
0.01 1 100 10 106 10 8
Figure 10.1: Estimated frequency of impacts on the Earth from the present population
of comets and asteroids and impact craters.
• smoke clouds,
Some more details about such scenarios can be found e.g. in Chapman and Mor-
rison (1994 [63]) or Melosh et al. (1990 [223]).
During the last 500 Million years there occurred several extinctions of marine
species: in Table 10.3 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 (∼ 1023 J energy release) and perhaps one ≤10 km
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 (≤ 10 times background values) at or near a
252 CHAPTER 10. ASTEROIDS, COMETS, METEROITES
Table 10.3: Extinction of marine species. The end of the stage is given in Myr.
Name End Name End
Pliocene 2.3 Mid Miocene 11.2
Upper Eocene 34 Maastrichtian 65
Cenomanian 93.5 Aptian 112
Tithonian 144 Callovian 159
Pliensbachian 190 Norian 206
Carnian 221 Tatarian 248
Guadelupian 250 Stephanian 290
Serpukhovian 322 Famennian 363
Frasnian 367 Eifelian 380
Ludlovian 411 Ashgillian 439
Llanvirnian 469 Tremadocian 493
Trempeleauean 505 Franconian 508
Botomian 520 Proterozoic/Cambrian 540
D...crater diameter, cf ...crater collapse factor (1.3 for craters larger than 4 km on
Earth), ge ...gravitational acceleration at the surface of the Earth, g... acceleration
at the surface of the body on which the crater is formed (in this case Earth),
W... the kinetic energy of the impacting body in kilotons TNT, ρa ...density of the
impactor (1.8 g/cm3 for a comet...7.3 g/cm3 for an iron body) and ρt ...density of
the target rock (e.g. 2.7 g/cm3 ). From this equation we see, that a carbonaceous
chondrite would produce a crater ∼ 94 km in diameter whereas an impactor of
cometary composition ∼ 150 km.
10.3 Meteorites
10.3.1 General Properties
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.
10.3. METEORITES 253
10.3.2 Classification
Meteorites can be classified into stony, stony iron and iron. The most common
meteorites are chondrites which are stony. Radiometric dating indicate an age of
about 4.5×109 years. Achondrites are also stony but they are considered differen-
tiated or reprocessed matter. They are formed by melting and recristallization on
or within meteorite parent bodies. Pallasites are stony iron meteorites composed
of olivine enclosed in metal.
The motion of meteoroids can be severely perturbed by the gravitational fields
of major planets. Jupiter’s gravitational influence is capable of reshaping an as-
teroid’s orbit from the main belt so that it dives into the inner solar system and
crosses the orbit of Earth. This is apparently the case of the Apollo and Vesta
asteroid fragments.
Particles found in highly correlated orbits are called stream components and
those found in random orbits are called sporadic components. It is thought that
most meteor streams are formed by the decay of a comet nucleus and consequently
are spread around the original orbit of the comet. When Earth’s orbit intersects
a meteor stream, the meteor rate is increased and a meteor shower results. Typ-
ically, a meteor shower will be active for several days. A particularly intense
meteor shower is called a meteor storm. Sporadic meteors are believed to have
had a gradual loss of orbital coherence with a meteor shower due to collisions and
radiative effects, further enhanced by gravitational influences. There is still some
debate concerning sporadic meteors and their relationship with showers. A well
known meteor shower are the Perseids (named after the radiant that is the name
of the constellation where the meteorites are seem to be coming from) which has
its maximum on August 11.
254 CHAPTER 10. ASTEROIDS, COMETS, METEROITES
40
30
10 v=11km/s
0
0.1. 1 10 100
Initial meteor radius [m]
Figure 10.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.
EMPs are created from the direct vaporization of impacting particles into
plasma and ESDs by a buildup of charging over the satellite surface. Both can
cause electrical and communication problems, erroneous signals in telemetry and
short circuits.
How can such problems be minimized? One simple manoeuvre is to minimize
the cross-sectional area of the satellite that is exposed to the meteor shower; e.g.
the solar panels point edge on into the meteor stream by reorienting the spacecraft.
Another technique is to turn off equipment that is particularly sensitive to ESDs.
Meteoroid impacts on spacecrafts and penetration damage are studied by
McBride and McDonnell (1999 [214]). Meteoroid morphology and density (e.g.
using NASA’s LDEF satellite results) were investigated by McDonnell and Gard-
ner (1998 [216]). The Leonid Meteor Shower and the Lunar Sodium Atmosphere
are treated in Hunten al., 1998 [143]. Charge production, impact probability, the
impact penetration depth, the amount of ionization generated upon impact on
spacecraft due to the strong Leonid shower in 1999 and 2000 are calculated by
McNeil et al., 1998 [221]
10.4 Comets
10.4.1 General Properties
Bright comets may be easily seen with the naked eye4 . Comets consist of different
parts:
• Hydrogen coma: not seen with naked eye, surrounds the visible gas/dust
coma, several million km.
Comets are named after their discoverers and are also given a number, the
year of discovery (or recovery) and a letter. In this designation 2006c denotes
the third comet discovered in 2006. There is another designation where comets
are distinguished by their perihelion passage: 2006X is the 10th comet in 2006
according to perihelion passage. Periodic comets are by a P/ and comets which
have collided with other bodies by D/.
4 see also: Comets and the Origin and Evolution of Life, by P. J. Thomas, Springer, 2006
256 CHAPTER 10. ASTEROIDS, COMETS, METEROITES
CO + hν → CO+ (10.3)
Since magnetic fields are carried by the solar wind, these ions are susceptible to
the magnetic force. Because the most common ion CO+ scatters blue light better
than red, the ion tail appears bluish. The solar wind sweeps past comets at 500
km/s
Comets from the Oort cloud may pass up to 400 times the inner solar system
before the object becomes perturbed by the inner planets and is changed into a
short period comet.
From the observed short period comets the total mass of bodies in the Kuiper
belt is estimated to be 0.0026 Earth masses which corresponds to 109 to 1010
objects. The Kuiper Belt is a disk-shaped region past the orbit of Neptune roughly
30 to 100 AU from the Sun containing many small icy bodies.
Figure 10.4: Comet observed by SOHO colliding with the Sun. SOHO/ESA, NASA
Sazonov and Yakovlev (2006 [270]) discuss comets on dangerous orbits and
how to change their orbits by explosive and sublimation methods. It is clear that
for large bodies the required charge power to destroy the dangerous bodies is too
large, and therefore, by means of an explosive impulse, such bodies can be moved
to a safe trajectory. The required destruction energy depends on
The energies are plotted in Fig. 10.5. For comparison: the energy of the Hiroshima
bombe was 15 kt. A cometary composition body of 100 m diameter requires 0.075
Mt for destruction. 1 kt TNT is the equivalent to 4,184×1012 J. The sum of
all bombs during the second world war can be estimated to 2 Mt. The biggest
H-bomb had 50 Mt.
Chapter 11
Space Debris
On Dec 3, 2001 BBC reports, that space debris lit up the sky. The spectacular
nighttime light show 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.
Orbital debris is defined as any man-made object in orbit around the Earth
which no longer serves as a useful purpose.
We will discuss estimations on the number of debris elements as well as models
that calculate their orbits and shielding mechanisms for spacecraft1 .
261
262 CHAPTER 11. SPACE DEBRIS
Figure 11.1: Orbits of GPS satellites. The satellites orbit the earth with a speed of
3.9 km/s. One revolution takes 12 h sidereal time, corresponding to 11 h 58 min earth
time. This means that the same satellite reaches a certain position about 4 minutes
earlier each day. The mean distance from the middle of the earth is 26560 km. The
system consists of at least 24 satellites, the first one started in 1978. The European
GALILEO will be operational in 2010, consisting of 30 satellites (27 working, 3 in reserve).
http : //[Link]/en/gps/[Link]
are used to study large scale phenomena such as hurricanes, or cyclones and for
communication satellites. The disadvantage of this type of orbit is that the Earth
can be observed from there with low resolution.
On the other hand, satellites in Low Earth Orbit, LEO, can only cover a small
area with high resolution and therefore often constellations are used, such as the
GPS, IRIDIUM (66 active communication satellites, first launched in 1998). The
current positions over 900 satellites can be followed using an online tool of NASA2 .
Figure 11.2: GEO and LEO objects as a source of space debris. GEO denotes geosta-
tionary orbit.
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.
From the above considerations it is clear that spacecrafts have to be protected
from collisions with space debris. Let us mention two examples: the US space com-
mand examines the trajectories of the Space Shuttle in order to identify possible
close encounters with space debris. If a dangerous object is believed to approach
a few tens of kilometers to the Space Shuttle, it will be maneuvered away from
the object (although in such a case the chances of a collision are only approxi-
mately 1:100 000). Such an operation is necessary about once every year or two
(at present).
Space debris is an inherently international problem and its solution requires
international co-operation. The Inter-Agency Space Debris Coordination Commit-
tee (IADC) whose members are ESA, NASDA (Japan), NASA, and the Russian
Space Agency RKA and the Canadian Space Agency (CSA) provides a forum for
discussion and coordination of technical space debris issues.
• Radar measurements: these have been used for space debris in low Earth
orbit (LEO).
• Optical measurements: these have been used for high Earth orbit (HEO). For
passive optical measurements the intensity of the signal from space debris is
inversely proportional to the square of its distance or altitude:
Figure 11.3: Spatial density of space debris by altitude according to ESA MASTER-
2001.
Therefore, optical telescopes 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.
• all weather,
• all day-and-night
2. Radars with electronically controlled beam direction using phased array an-
tennas. In that case multiple objects at different directions can be detected
and measured simultaneously; used for tracking and search tasks.
11.2. DETECTION OF SPACE DEBRIS 265
In the tracking mode the radar follows an object for a few minutes gaining data
on angular direction, range, range rate, amplitude and phase of the radar echoes.
From these parameters the orbital elements can be derived In the beam-park mode,
the antenna is kept fixed in a given direction and echoes are received from objects
passing within its field of view. This yields statistical information on the number
and size of detected objects; the determination of the orbit is less precise. There
is also a mixed mode. From the radar measurements the following parameters can
be derived:
1. orbital elements; thus the motion of the object’s center of mass around Earth
is defined.
2. Attitude; describes the motion of the object around its center of mass.
4. Ballistic coefficient; this describes the rate at which the orbital semi-major
axis decays.
5. Object mass,
6. material properties.
The main source of data for space debris in the size range of 1-30 cm is the
NASA Haystake radar facility operated by MIT Lincoln Laboratory. Under an
agreement with the US Air Force since 1990 data are collected. The data indicate
that there are about 100 000 fragments in orbits with sizes down to 1 cm.
When space debris (man made or meteoroids) enter the atmosphere an ion-
ization trail is created. Molecules in the upper atmosphere are ionized by the
passage of the meteor. Such ionization trails can last up to 45 minutes. Such an
ion trail will act as a mirror for radio waves. Radar measurements operating at
50 MHz permit to estimate the reentry of space debris, the ionization trail behind
reentering bodies can be detected. With such facilities one can detect meteors as
small as 100 microns.
Radar measurements of space debris have been done at Haystack (US) and
Goldstone radars (US), Russia and by Germany using the Research Establishment
for Applied Science (FGAN) radar and the Effelsberg radio telescope. Haystack
and Goldstone radars have provided a statistical picture of LEO debris at sizes
down to 0.5 cm which was confirmed by FGAN. These measurements have proven
that the debris population exceeds the natural meteoroid population for all sizes
(except between 30 and 500 m).
Radar measurements and the usage of the Sardinian Radio Telescope for space
debris detection are described e.g. in the article of Di Martino et al., 2006 [74].
The international radar space debris research was reviewed by Molotov et al.
(2005 [224]) where further literature is cited.
266 CHAPTER 11. SPACE DEBRIS
11.2.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 debris3 : 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 14 l 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.
11.2.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.
Table 11.1: Mean time between impacts on a satellite with a cross-section area of 10 m2
Height of Objects Objects Objects
circular orbit 0.1-1.0 cm 1-10 cm >10 cm
500 km 10-100 yrs 3 500-7 000 yrs 150 000 yrs
1 000 km 3-30 yrs 700 - 1 400 yrs 20 000 yrs
1 500 km 7-70 yrs 1 000-2 000 yrs 30 000 yrs
For GEO the situation is more complicated. The number of space debris of less
than 1 m in diameter is not well known. Moreover, there is no natural removal
mechanism for satellites in GEO. One can estimate an annual collision probability
for an average operational satellite with other catalogued objects at 10−5 .
Another problem concerns the re-entry. Since the last 40 years 16 000 re-
entries of catalogued space objects are recorded. 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 1 m2 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.
A family of space debris objects was found at a height of 900 km. The density
peak found there is caused by a large number of sodium-potassium liquid metals
droplets- they have been used as a coolant for the on board nuclear reactor, leaked
from the Russian ocean surveillance satellites. The estimation is about 70000 drops
with diameters between 0.5 mm and 5.5 and the detection was mainly made with
the Haystack radar. Using the Goldstone radar the so called West Ford Needles at
an altitude of 2900 km were detected. They are copper dipoles, 1.77 cm long, and
remnants that were released in 1961 and 1963 by the US MIDAS 3 and MIDAS
6 satellites for telecommunication experiments. It was first expected that they
should reenter the Earth’s atmosphere within 5 years but now a population of
40000 objects were found between 2400 and 3100 km.
In total it is estimated that more than 350000 objects larger than 1 mm crowd
the space around Earth and a particle around 5 mm is able to directly penetrate
the shuttle cabin. For more details see Valsecchi et al., 2006 [323].
268 CHAPTER 11. SPACE DEBRIS
Collision risk assessment for a spacecraft in space debris environment was stud-
ied e.g. by Tang et al., 2005 [311]
• Debris left in orbits below 600 km: normally falls back to Earth within a few
years.
• Debris left in orbits at altitudes of 800 km: the time for orbital decay is
several decades.
• Debris left in orbits at altitudes above 1 000 km: will normally continue
circling the Earth for a century or more.
• Shield development.
• 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 µm 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. 11.4 where the coordinates are
logarithmic.
270 CHAPTER 11. SPACE DEBRIS
10 4
SMM
EURECA
-4
10
-8
10
-2
10 -4 10 1
Diameter [cm]
Figure 11.4: Approximate measured debris flux in low Earth orbit by object size (sketch)
• Optical measurements; these are suited for objects in high altitude orbits.
The detectors use CCD and a 1 m telescope will be operated by ESA at the
Teide observatory in Tenerife.
The main aim of the mathematical model is a description of the debris and me-
teoroid environment at altitudes between low Earth orbit (LEO) and the geosta-
tionary orbit (GEO). The minimum size of an object is 0.1 mm. The model is
based on the catalogued population and on known break-ups of spacecraft and
rocket upper-stages in orbit. The initial distribution of fragments is described in
terms of their position, velocity, mass. The objects are then propagated forward in
time taking into account the relevant perturbations. The description of MASTER
(Meteoroid and Space Debris Terrestrial Reference Model) is given in Sdunnus et
al. 2001 [276].
272 CHAPTER 11. SPACE DEBRIS
Anomalous debris
1%
Fragmentation
31% Payloads
debris 38%
12% 18%
Operational debris Rocket bodies
11.3.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,
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
11.3. SHIELDING AND RISK ASSESSMENTS 273
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.
The necessity for collision avoidance manoeuvres was already pointed out by
Rex et al., 1991 [255].
An overview of fragmentation of LEO Upper Stages was given by Chernyavskiy
et al. (1994 [64]).
A technical report on space debris was given from the Scientific and Tech-
nical Subcommittee of the United Nations Committee on the Peaceful uses of
Outer Space (1999). The European initiatives and space mitigation standard are
reviewed by Alby et al. 2004 [3].
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.
The necessity of debris mitigation is illustrated in the ESA Space Debris Miti-
gation Handbook 2002 ( Klinkrad et al., 2004 [167]). Algorithms for the derivation
of appropriate collision avoidance strategies are presented by Sánchez-Ortiz et al.
(2004 [269]).
Bibliography
[1] Y. Abe and T. Matsui. Early evolution of the earth: Accretion, atmosphere
formation, and thermal history. Journal of Geophys. Research, 91:291, Sep-
tember 1986.
[2] Y. Abe and T. Matsui. Evolution of an Impact-Generated H2 O-CO2 At-
mosphere and Formation of a Hot Proto-Ocean on Earth. Journal of At-
mospheric Sciences, 45:3081–3101, November 1988.
[3] F. Alby, D. Alwes, L. Anselmo, H. Baccini, C. Bonnal, R. Crowther,
W. Flury, R. Jehn, H. Klinkrad, C. Portelli, and R. Tremayne-Smith. The
European Space Debris Safety and Mitigation Standard. Advances in Space
Research, 34:1260–1263, 2004.
[4] V. A. Alexeev and G. K. Ustinova. Meteorite Data on the Solar Modulation
of Galactic Cosmic Rays and an Inference on the Solar Activity Influence
on Climate of the Earth. In S. Mackwell and E. Stansbery, editors, 36th
Annual Lunar and Planetary Science Conference, page 1012, March 2005.
[5] R. B. Alley and P. U. Clark. The Deglaciation of the Northern Hemisphere:
A Global Perspective. Annual Review of Earth and Planetary Sciences,
27:149–182, 1999.
[6] L. W. Alvarez, W. Alvarez, F. Asaro, and H. V. Michel. Extraterrestrial
Cause for the Cretaceous Tertiary Extinction. Science, 208:1095, 1980.
[7] E. I. Alves and A. R. Baptista. Rock Magnetic Fields Shield the Surface of
Mars from Harmful Radiation. In S. Mackwell and E. Stansbery, editors,
Lunar and Planetary Institute Conference Abstracts, page 1540, March 2004.
[8] D. M. Amron and R. L. Moy. Stratospheric ozone depletion and its re-
lationship to skin cancer. Journal of Dermatoligic Surgery and Oncology,
17:370–372, 1991.
[9] L. Andersson, L. Eliassson, and P. Wintoft. In Proccedings of the Workshop
on Space Weather, 11-13 November 1998, ESTEC, Noordwijk WP-155, 1999.
[10] A. L. Andrady, S. H. Hamid, X. Hu, and A. Torikai. Effects of increased
solar ultraviolet radiation on materials. Journal of Photochemistry and Pho-
tobiology B, 1998.
275
276 BIBLIOGRAPHY
[27] R. E. Basher. Basic science of solar radiation and its ultraviolet components.
In Proceedings of the seminar on solar ultraviolet radiation, R. Basher (Ed),
NZMS,, 1981.
[28] D. Basu. Solar neutrino in relation to solar activity. Solar Physics, 142:205–
208, November 1992.
[29] D. Basu. Radius of the Sun in relation to solar activity. Solar Physics,
183:291–294, December 1998.
[32] J. Beer. Long-term indirect indices of solar variability. Space Science Re-
views, 94:53–66, November 2000.
[33] A. O. Benz. Coronal Heating and Microflares in Solar Active and Quiet
Regions. Astronomische Nachrichten Supplement, 324:8, July 2003.
[48] P. D. Bradley and E. Normand. IEEE Trans. Nucl. Sci., 45:6, 1998.
[50] L. F. Burlaga. Magnetic fields and plasmas in the inner heliosphere: Helios
results. Planetary and Space Science, 49:1619–1627, December 2001.
[95] D. H. Fairfield and L. J. Cahill, Jr. Transition Region Magnetic Field and
Polar Magnetic Disturbances. Journal of Geophys. Research, 71:155–+, Jan-
uary 1966.
[98] B. Fleck. Highlights from SOHO and Future Space Missions. In The Dy-
namic Sun, Proceedings of the Summerschool and Workshop held at the Solar
Observatory, page 1, May 2001.
[99] A. Fossum and M. Carlsson. High-frequency acoustic waves are not sufficient
to heat the solar chromosphere. Nature, 435:919–921, June 2005.
[101] P. Foukal and J. Lean. An empirical model of total solar irradiance variation
between 1874 and 1988. Science, 247:556–558, February 1990.
[107] J.-C. Gerard. Modelling the Climatic Response to Solar Variability. Royal
Society of London Philosophical Transactions Series A, 330:561–573, April
1990.
BIBLIOGRAPHY 283
[109] S. E. Gibson. Global Solar Wind Structure from Solar Minimum to Solar
Maximum: Sources and Evolution. Space Science Reviews, 97:69–79, 2001.
[111] R. L. Gilliland. Solar radius variations over the past 265 years. Astrophysical
Journal, 248:1144–1155, September 1981.
[121] E. F. Guinan and I. Ribas. Our Changing Sun: The Role of Solar Nu-
clear Evolution and Magnetic Activity on Earth’s Atmosphere and Climate.
In ASP Conf. Ser. 269: The Evolving Sun and its Influence on Planetary
Environments, page 85, 2002.
284 BIBLIOGRAPHY
[141] D. V. Hoyt and K. H. Schatten. The role of the sun in climate change. The
role of the sun in climate change / Douglas V. Hoyt, Kenneth H. Schat-
ten. New York : Oxford University Press, 1997. QC883.2.S6 H69 1997,
1997.
[142] M. K. Hughes and H. F. Diaz, editors. The medieval warm period, 1994.
[145] J. Hyde and E. Christiansen. Space Shuttle Meteoroid and Orbital Debris
Threat Assessment Procedure. In COSPAR, Plenary Meeting, 2002.
[147] B. V. Jackson, P. P. Hick, and A. Buffington. Halo CME’s - Will They Hit
or Miss Earth? AGU Fall Meeting Abstracts, pages A474+, December 2002.
[151] W. Kalkofen. The heating of the quiet solar chromosphere. In IAU Symp.
142: Basic Plasma Processes on the Sun, pages 197–204, 1990.
[153] J. F. Kasting. Habitable Zones around Low Mass Stars and the Search
For Extraterrestrial Life. Origins of Life and Evolution of the Biosphere,
27:291–3107, 1997.
[159] R. A. Kerr. Sunspot Weather Link Holding UP. Science, 242:1124, November
1988.
[162] R.-S. Kim, K.-S. Cho, Y.-J. Moon, Y.-H. Kim, Y. Yi, M. Dryer, S.-C. Bong,
and Y.-D. Park. Forecast evaluation of the coronal mass ejection (CME)
geoeffectiveness using halo CMEs from 1997 to 2003. Journal of Geophysical
Research (Space Physics), 110:11104, November 2005.
[163] R. Kippenhahn and A. Schlüter. Eine Theorie der solaren Filamente. Mit 7
Textabbildungen. Zeitschrift fur Astrophysics, 43:36, 1957.
[178] K. Kudela and M. Storini. Possible tools for space weather issues from cosmic
ray continuous records. Advances in Space Research, 37:1443–1449, 2006.
[179] K. Labitzke. Sunspots, the QBO, and the stratospheric temperature in the
north polar region. Geophysical Research Letter, 14:535–537, 1987.
[183] M. Lagache. New data on the kinetics of the dissolution of alkali feldspars at
200 C in CO 2 charged water. Geochimica et Cosmochimica Acta, 40:157–
161, February 1976.
[185] H. H. Lamb. Climate Present, Past and Future, Vol. 2 Climatic history and
future. Methuen, London, 1977.
[188] P. Lantos. Predictions of Galactic Cosmic Ray Intensity Deduced from that
of Sunspot Number. Solar Physics, 229:373–386, July 2005.
[189] J. Lean. The Sun’s Variable Radiation and Its Relevance For Earth. Annual
Review of Astronomy and Astrophysics, 35:33–67, 1997.
[202] O.v.d. Lühe. In The Dynamic Sun, Proceedings of the Summerschool and
Workshop held at the Solar Observatory, pages 43–68, May 2001.
[203] H. Lundstedt. The sun, space weather and GIC effects in Sweden. Advances
in Space Research, 37:1182–1191, 2006.
[204] T. J. Lydon and S. Sofia. A Method for Incorporating the Effects of Large-
Scale Magnetic Fields in the Study of Stellar Structure and Variability. As-
trophysical Journal Supplement, 101:357, December 1995.
290 BIBLIOGRAPHY
[239] M. Ossendrijver. The solar dynamo. The Astronomy and Astrophysics Re-
view, 11:287–367, 2003.
[249] E. R. Priest. Solar Flare Theory and the Status of Flare Understanding. In
ASP Conf. Ser. 206: High Energy Solar Physics Workshop - Anticipating
Hess!, page 13, 2000.
[250] A. Pulkkinen, S. Lindahl, A. Viljanen, and R. Pirjola. Geomagnetic storm
of 29-31 October 2003: Geomagnetically induced currents and their relation
to problems in the Swedish high-voltage power transmission system. Space
Weather, 3:8, August 2005.
[251] M. P. Rast. The Scales of Granulation, Mesogranulation, and Supergranu-
lation. Astrophysical Journal, 597:1200–1210, November 2003.
[252] D. M. Raup and J. J. Sepkoski. Testing for Periodicity of Extinction. Science,
241:94–96, July 1988.
[253] R. C. Reedy. Gamma Rays from Mercury Made by Solar Energetic Particles.
Bulletin of the American Astronomical Society, 29:987, July 1997.
[254] G. C. Reid. Solar total irradiance variations and the global sea surface
temperature record. Journal of Geophys. Research, 96:2835–2844, 1991.
[255] D. Rex, J. Bendisch, P. Eichler, and J. Zhang. Protecting and manoeuvring
of spacecraft in space debris environment. Advances in Space Research,
11:53–62, 1991.
[256] E. Ribes, B. Beardsley, T. M. Brown, P. Delache, F. Laclare, J. R. Kuhn,
and N. V. Leister. The variability of the solar diameter, pages 59–97. The
Sun in Time, 1991.
[257] E. Ribes, J. C. Ribes, and R. Barthalot. Evidence for a larger sun with a
slower rotation during the seventeenth century. Nature, 326:52–55, March
1987.
[258] A. J. Ridley, Y. Deng, and G. Tóth. The global ionosphere–thermosphere
model. Journal of Atmospheric and Terrestrial Physics, 68:839–864, May
2006.
[259] M. Rieutord, T. Roudier, J. M. Malherbe, and F. Rincon. On mesogranula-
tion, network formation and supergranulation. Astronomy and Astrophysics,
357:1063–1072, May 2000.
[260] E. Riklis, I. Emerit, and R. B. Setlow. New approaches to biochemical
radioprotection: antioxidants and DNA repair enhancement. Advances in
Space Research, 18:51–54, 1996.
[261] C. B. Roald, P. A. Sturrock, and R. Wolfson. Coronal Heating: Energy Re-
lease Associated with Chromospheric Magnetic Reconnection. Astrophysical
Journal, 538:960–967, August 2000.
[262] T. Roudier, J. M. Malherbe, L. November, J. Vigneau, G. Coupinot, M. La-
fon, and R. Muller. Intergranular plumes and formation of network bright
points. Astronomy and Astrophysics, 320:605–611, April 1997.
294 BIBLIOGRAPHY
[271] N. Scafetta and B. J. West. Climate sensitivity to solar activity: The con-
tribution of solar cycles 21-23 to global mean surface warming. AGU Fall
Meeting Abstracts, December 2004.
[272] K. Schatten. Fair space weather for solar cycle 24. Geophysical Research
Letter, 32:21106, November 2005.
[274] C. J. Schrijver and C. Zwaan. Solar and stellar magnetic activity. Solar and
stellar magnetic activity / Carolus J. Schrijver, Cornelius Zwaan. New York:
Cambridge University Press, 2000. (Cambridge astrophysics series ; 34),
2000.
[279] P. M. Sforza and J. L. Remo. NEO Mission Dynamics and Advanced Space
Propulsion. In Near-Earth Objects, page 432, January 1997.
[281] N. J. Shaviv. On climate response to changes in the cosmic ray flux and
radiative budget. Journal of Geophysical Research (Space Physics), 110:8105,
August 2005.
[282] H. Shibahashi. Solar structure and the neutrino problem, pages 231–246.
Stellar astrophysical fluid dynamics. Edited by Michael J. Thompson, Jørgen
Christensen-Dalsgaard. Cambridge (UK): Cambridge University Press, ISBN
0-521-81809-5, 2003, p. 231 - 246, 2003.
[290] M. Sobotka. Fine Structures in Sunspots. In ASSL Vol. 239: Motions in the
Solar Atmosphere, pages 71–97, 1999.
[294] C. P. Sonett. Sunspot time series - Spectrum from square law modulation of
the Hale cycle. Geophysical Research Letter, 9:1313–1316, December 1982.
[296] Q. Song and W. Cao. In Astrophys. in Asia and its Future, Ed. P.S. Chen,
Yunan Obs. Chin. Acadm. Sc., page 139, 1999.
[297] E. A. Spiegel and J.-P. Zahn. The solar tachocline. Astronomy and Astro-
physics, 265:106–114, November 1992.
[299] B. Stauffer. Long Term Climate Records from Polar Ice. Space Science
Reviews, 94:321–336, November 2000.
[303] Y. Suzuki. Solar Neutrinos. Space Science Reviews, 85:91–104, August 1998.
[308] E. Tajika. Physical and Geochemical Conditions Required for the Initiation
of Snowball Earth. In Lunar and Planetary Institute Conference Abstracts,
page 1573, March 2001.
[309] M. Takata. In Frontier of Neutrino Astrophysics, 1993.
[310] M. Takata and H. Shibahashi. Solar Models Based on Helioseismology and
the Solar Neutrino Problem. Astrophysical Journal, 504:1035, September
1998.
[311] Q. Tang, B. J. Pang, and W. Zhang. Collision Risk Assessment for a Space-
craft in the Space Debris Environment. In D. Danesy, editor, ESA SP-589:
4th European Conference on Space Debris, page 721, August 2005.
[312] I. Thomson. Mutat. Res., 430, 1999.
[313] Y. Totsuka. KAMIOKANDE and Super-Kainiokande. In The Sun and Be-
yond, page 221, 1996.
[314] L. Townsend. Space Radiation Hazards on Human Missions to the Moon
and Mars. AGU Fall Meeting Abstracts, pages C1+, December 2004.
[315] L. W. Townsend and J. W. Wilson. The Interplanetary Radiation Environ-
ment and Methods to Shield from it. In C. R. Stoker and C. Emmart, editors,
Strategies for Mars: A Guide to Human Exploration, page 283, 1996.
[316] D. Tripathi, V. Bothmer, S. K. Solanki, R. Schwenn, M. Mierla, and G. Sten-
borg. Plasma dynamics of a prominence associated coronal mass ejection.
In IAU Symposium, pages 401–402, November 2004.
[317] S. Turck-Chièze, S. Couvidat, A. G. Kosovichev, A. H. Gabriel,
G. Berthomieu, A. S. Brun, J. Christensen-Dalsgaard, R. A. Garcı́a, D. O.
Gough, J. Provost, T. Roca-Cortes, I. W. Roxburgh, and R. K. Ulrich. Solar
Neutrino Emission Deduced from a Seismic Model. Astrophysical Journal
Letters, 555:L69–L73, July 2001.
[318] S. Turck-Chièze, R. A. Garcı́a, S. Couvidat, R. K. Ulrich, L. Bertello,
F. Varadi, A. G. Kosovichev, A. H. Gabriel, G. Berthomieu, A. S. Brun,
I. Lopes, P. Pallé, J. Provost, J. M. Robillot, and T. Roca Cortés. Look-
ing for Gravity-Mode Multiplets with the GOLF Experiment aboard SOHO.
Astrophysical Journal, 604:455–468, March 2004.
[319] S. Turck-Chieze and I. Lopes. Toward a unified classical model of the sun
- On the sensitivity of neutrinos and helioseismology to the microscopic
physics. Astrophysical Journal, 408:347–367, May 1993.
[320] P. Ulmschneider, E. R. Priest, and R. Rosner, editors. Mechanisms of Chro-
mospheric and Coronal Heating, 1991.
[321] R. K. Ulrich and L. Bertello. Solar-Cycle Dependence of the Sun’s Appar-
ent Radius in the Neutral Iron Spectral Line at 525-NM. Nature, 377:214,
September 1995.
298 BIBLIOGRAPHY
[333] R. Wieler and H. Baur. Fractionation of Xe, Kr, and AR in the Solar
Corpuscular Radiation Deduced by Closed System Etching of Lunar Soils.
Astrophysical Journal, 453:987, November 1995.
[334] R. Wieler, P. Etique, P. Signer, and G. Poupeau. Record of the solar cor-
puscular radiation in minerals from lunar soils - A comparative study of
noble gases and tracks. In Lunar and Planetary Science Conference, pages
1369–1393, 1980.
[343] C.-T. Yeh, M. D. Ding, and P. F. Chen. Waiting time distribution of CMEs.
In Solar Magnetic Phenomena, Proceedings of the 3rd Summerschool and
Workshop held at the Solar Observatory Kanzelhöhe, Kärnten, Austria, Au-
gust 25 – September 5, 2003. Edited by A. Hanslmeier, A. Veronig, and M.
Messerotti. Astronomy and Astrophysics Space Science Library, vol. 320,
Springer, Dordrecht, The Netherlands, 2005., p.171-174, pages 171–174,
2005.
[Link]
[Link]
[Link]
[Link]
Lund Space Weather Center
[Link]
[Link]
301
List of Tables
303
304 LIST OF TABLES
11.1 Mean time between impacts on a satellite with a cross-section area of 10 m2 267
11.2 Spectacular satellite crashes on Earth. . . . . . . . . . . . . . . . . . . 269
11.3 Some examples of retrieved spacecraft and surfaces . . . . . . . . . . . 270
Index
305
306 INDEX
Volume 338: Solar Journey: The Significance of Our Galactic Environment for the
Heliospere and Earth, edited by P. Frisch.
Hardbound ISBN 978-1-4020-4397-0, September 2006
Volume 334: The New Astronomy: Opening the Electromagnetic Window and
Expanding our View of Planet Earth, edited by W. Orchiston.
Hardbound ISBN 1-4020-3723-6, October 2005
Volume 327: The Initial Mass Function 50 Years Later, edited by E. Corbelli,
F. Palla, H. Zinnecker.
Hardbound ISBN 1-4020-3406-7, June 2005
Volume 317: The Sun and the Heliosphere as an Integrated System, edited by
G. Poletto, S.T. Suess.
Hardbound ISBN 1-4020-2830-X, November 2004
Volume 315: How does the Galaxy Work - A Galactic Tertulia with Don Cox and
Ron Reynolds, edited by E. J. Alfaro, E. Pérez, J. Franco
Hardbound ISBN 1-4020-2619-6, September 2004
Volume 314: Solar and Space Weather Radiophysics- Current Status and Future
Developments, edited by D.E. Gary and C.U. Keller
Hardbound ISBN 1-4020-2813-X, August 2004
Volume 313: Adventures in Order and Chaos, by G. Contopoulos.
Hardbound ISBN 1-4020-3039-8, January 2005
Volume 309: Soft X-ray Emission from Clusters of Galaxies and Related
Phenomena, edited by R. Lieu and J. Mittaz
Hardbound ISBN 1-4020-2563-7, September 2004
Volume 304: Cosmic Gammy-ray Sources, edited by K.S. Cheng and G.E. Romero
Hardbound ISBN 1-4020-2255-7, September 2004
Volume 299: Open Issues in Local Star Fomation, edited by Jacques Lépine,
Jane Gregorio-Hetem.
Hardbound, ISBN 1-4020-1755-3, December 2003
Volume 293: Physics of the Solar System, by Bruno Bertotti, Paolo Farinella,
David Vokrouhlický
Hardbound, ISBN 1-4020-1428-7, August 2003
Paperback, ISBN 1-4020-1509-7, August 2003
Volume 287/8/9: The Future of Small Telescopes in the New Millennium, edited
by Terry D. Oswalt.
Hardbound Set only of 3 volumes, ISBN 1-4020-0951-8, July 2003
Volume 286: Searching the Heavens and the Earth: The History of Jesuit
Observatories, by Agustín Udías
Hardbound, ISBN 1-4020-1189-X, October 2003
Volume 284: Light Pollution: The Global View, edited by Hugo E. Schwarz
Hardbound, ISBN 1-4020-1174-1, April 2003
Volume 274: New Quests in Stellar Astrophysics: The Link Between Stars and
Cosmology, edited by Miguel Chávez, Alessandro Bressan, Alberto
Buzzoni,Divakara Mayya
Hardbound, ISBN 1-4020-0644-6, June 2002
Volume 270: Dayside and Polar Cap Aurora, by Per Even Sandholt,
Herbert C. Carlson, Alv Egeland
Hardbound, ISBN 1-4020-0447-8, July 2002
For further information about this book series we refer you to the following web site:
[Link]