COSMIC RAYS AND THE EARTH’S ATMOSPHERE
A.D. Erlykin1,2 and A.W. Wolfendale2 1 P.N. Lebedev Physics Institute, Leninsky Prospekt, Moscow, Russia. 2 Department of Physics, University of Durham, South Road, Durham, DH1 3LE, U.K.
Abstract Avery brief summary is given of aspects of cosmic ray physics which have rel- evance to the possible effects of cosmic rays on ‘climate’. It is concluded that a more detailed look at the effect of fast ionizing particles on the atmosphere from the standpoint of cloud production would be advantageous.
1. INTRODUCTION There have been many claims for a correlation between solar properties (e.g. sunspot number) and climate but all have suffered from the absence of a reasonable physical cause, the point being that the energy changes in the solar variations are deemed to be too small to account for the necessary climate forcing. This is not to say that there are no well-documented effects on very long time scale; there are. Those due to variations in the sun-earth distance and the inclination of the earth’s axis (Milankovic effects), which relate to 103 − 106 y periods, are generally agreed. What is not (yet) agreed is that the 11-year solar cycle has a significant correlation with climate. The best evidence favouring a specific cause for a sunspot (SS) — climate connection relates to the apparent role of cosmic rays which are, themselves, modulated by solar activity. The observation is that by the Danish group in which there is a correlation of cloud cover and CR intensity over the oceans. The likelihood of this being genuine comes from the fact that CR are the major source of ionization away from land and CR, of course, provide ionization. Insofar as the cloud cover/CR intensity results are considered in detail elsewhere, more discussion will not be given here; rather, we will concentrate on the CR aspects.
2. COSMIC RAY INTENSITY AS A FUNCTION OF ATMOSPHERIC DEPTH It is relevant to consider the manner in which the vertical intensity of cosmic rays varies with height in the atmosphere. Considering the three major components: protons, electrons and muons, the values of −2 −1 −1 IV , the vertical intensity in cm s sr ,atheights above sea level of 2, 5 and 10 km, respectively, are:
− − − p :2× 10 4, 1.6 × 10 3 and 2 × 10 2 − − − e :5× 10 3, 3 × 10 2 and 2 × 10 1 − − − µ :1× 10 2, 2 × 10 2 and 5.5 × 10 2.
Some comments can be made, as follows: (i) The peak of the ionization (for µ and e)isinthe 10 km region, much higher than the common cloud level. Such an observation is not ‘the kiss of death’ to the correlation idea, because of uncertainty in transport phenomena for the products of CR ionization between the 10 km level and much lower levels. (ii) Although the ionization produced by protons is lower than that produced by e it is important to point out that the rare ‘stars’, produced by proton interactions in the atmosphere, contain very highly ionizing nuclear fragments. It is not inconceivable that subtle effects leading to cloud droplets are associated with these highly ionizing fragments. 3. KNOWN METEOROLOGICAL EFFECTS Concerning the CR intensity at ground level there are three major ‘meteorological variations’: (i) The pressure coefficient, due simply to absorption of the secondaries, of −2% per cm Hg. (ii) A correlation with the height of the 100 mb level, amounting to ∼−5% per km. The reason is to do with µ − e decay. (iii) The mean air temperature between the 100 and 200 mb levels. The dependence is +0.1% per ◦K and is due to π − µ decay. All the effects are well understood by the cosmic ray fraternity. However, they should be borne in mind by the wider community when correlations are sought.
4. COSMIC RAY ORIGIN — AND EFFECTS 4.1 Galactic particles Below the ‘knee’ in the energy spectrum, at ∼ 3 × 1015 eV, it is probable that most CR come from supernova remnants (SNR) by way of shock acceleration. It is these particles, mainly — specifically, below about 1012 eV — that are modulated by the solar wind which has, itself, an 11-year cycle. It is inevitable that there should be intensity variations due to the stochastic nature of SNR but these should be rare. There has been a claim for a 2-fold increase in intensity some 35 thousand years ago but it seems likely (Beer, private communication) that the increase is due not to an SNR but to a variation of the earth’s magnetic field and/or solar variability.
4.2 Solar particles With the advent of space vehicles, the study of the (mainly) low energy solar CR has become a ‘growth industry’. The range-energy relation is such that most protons, or their progeny, do not reach ground level; even at 10 GeV the particles only reach a height of about 10 km. Nevertheless, solar CR are important, particularly in the polar regions where effects on the ozone layer have been claimed, and their presence is eminently reasonable. Finally, we can remark on the possibility of very rare solar flares having serious effects on climate and, indeed, mankind itself. An extrapolation of the log N − log S curve for energy deposited on earth above S would indicate very serious effects every million years, or so. However, this interval is surely too short (otherwise we would not have survived for so long!). What can be said is that significant effects might be expected every 1000 years, or so. Effects on climate might occur if the atmosphere happened to be in an unstable phase at the time.
5. CONCLUSIONS Cosmic ray effects offer the possibility of being relevant to climate change. Although it is premature to be dogmatic, the likelihood of significant climatic effects is high enough for a detailed analysis of the physics — and meteorology — of CR-air interactions to be not just desirable, but vital.
ACKNOWLEDGEMENTS The authors are grateful to Jasper Kirkby for re-kindling our interest in this topic. ICE CORE DATA ON CLIMATE AND COSMIC RAY CHANGES
J. Beer Federal Institute of Environmental Science and Technology, EAWAG, CH-8600 Dübendorf, Switzerland, tel: +41 1 823 51 11 / fax: +41 1 823 52 10, email: [email protected]
Abstract Ice cores represent archives which contain unique information about a large variety of environmental parameters. Climatic information is stored in the form of stable isotopes, greenhouse gases and various chemical substances. The content of cosmogenic nuclides such as 10Be and 36Cl provide long-term records of the intensity of the cosmic ray flux and its modulation by solar activity and the geomagnetic dipole field. Cosmogenic nuclides are produced by the interaction of cosmic ray particles with the atmosphere. After production, these nuclides are transported and distributed within the environment, depending on their geochemical properties. Some of them are removed from the atmosphere by snow and incorporated into ice sheets and glaciers. The analysis of the Greenland ice cores GRIP and GISP2 are discussed in terms of climate and cosmic ray changes during the past 50’000 years.
1. ARCHIVE ICE Polar ice sheets are formed from snow. The snowflakes grow together to grains which slowly increase in size. Due to the pressure of the overlying new snow layers, the grains become more and more compacted and finally turn into ice. The consequence of this formation process is that the ice not only preserves all the atmospheric constituents such as aerosols and dust, it also contains air bubbles that enable to determine the atmospheric composition and in particular the reconstruction of greenhouse gases in the past. This unique property makes ice the only archive that virtually stores all the climate forcing factors (greenhouse gases, aerosols and volcanic dust, solar irradiance) except internal variability. Ice cores also contain information on the corresponding climate response (temperature, precipitation rate, wind speed, atmospheric circulation). Another important property of ice is that it flows. This can be seen in Fig. 1, which schematically depicts an ice-sheet. The ice slowly flows towards the margin of the ice sheet, where it partly melts and partly breaks up as icebergs. Under steady-state conditions, the ice lost in the ablation area is replaced by snow falling on the accumulation area where new layers are formed continuously. As a consequence of the horizontal movement of the ice, the annual layers become thinner with increasing depth, as indicated in Fig. 1. This leads to another special property of the archive ice. The depthÐage relationship is non- linear, which has the advantage that the uppermost part of the core is well resolved and the total time period covered is long (of the order of 105 years for polar ice cores). The disadvantage of this non-linear time-scale is, however, that dating ice is difficult and relies strongly on correct modeling of the ice-flow. The main ice sheets are situated in polar regions (Greenland, with a maximum thickness of approx. 3 km and Antarctica, with a thickness of up to 4 km). Smaller ice sheets at lower latitudes can only be found at high altitudes (Andes, Himalayas, Alps) [1]. There is a steadily growing number of parameters which can be measured in ice cores. It is beyond the scope of this paper to discuss all these parameters. In Table 1, a small selection of those related to climate forcing and climate response is given. ❄❄❄❄❄❄❄❄ ❄ ❄ ❄❄❄ ❄ ❄❄❄
ACCUMULATION
ABLATION ABLATION
BEDROCK
Figure 1: Formation of an ice sheet. The snow falling in the accumulation region turns into ice that slowly flows towards the ablation area where it breaks up into ice-bergs or melts. As a consequence of the flow characteristics the thickness of annual layers decreases with increasing depth.
Table 1. Climate parameters measured in ice cores.
Parameter Proxy for
CO2 Greenhouse gases
CH4 Greenhouse gases
SO4 Volcanic eruptions Ash Volcanic eruptions 10Be, 36Cl Solar activity 18 d O Temperature Borehole temperature Temperature Annual layer thickness Precipitation rate Dust Wind speed Anions / cations Atmospheric circulation
18 18 As an example, d O of the GRIP ice core is shown in Fig. 2. d O (relative deviation of the 18O/16O ratio in ice from a standard in ‰) reflects mainly the temperature at which snow is formed.
-30
-35
18 d O [‰]
-40
-45 01020304050 Age [ky BP]
18 Figure 2: d O measured in the GRIP ice core from Greenland. Low values indicate cold climate. During the last 10'000 years the temperature was relatively stable compared to the preceding glacial period. Figure 2 shows that during glacial times the temperature in Greenland was characterized by abrupt changes (so-called Dansgaard-Oeschger events) of up to 20∞ C within a few decades. The last 10’000 years, the so-called Holocene, however, looks comparatively stable. The Dansgaard- Oeschger events were probably caused by abrupt changes in the ocean circulation, transporting heat to high latitudes. In the following, we will concentrate on what cosmogenic radionuclides in ice cores can tell us.
2. COSMOGENIC RADIONUCLIDES IN ICE The cosmic ray particles (87% protons, 12% helium nuclides, 1% heavier particles) that enter the Earth’s atmosphere react with Nitrogen, Oxygen and Argon, producing a cascade of secondary particles. These nuclear processes produce a variety of cosmogenic nuclides such as 10Be, 14C and 36Cl. These nuclides are listed in Table 2 together with their main properties.
Table 2. Some cosmogenic radionuclides and their main properties.
Nuclide Half-life Target Production rate (years) (atoms cm-2 s-1) 10Be 1.5 106 N, O 0.018 14C 5730 N, O 2.0 36Cl 3.01 105 Ar 0.0019 The physics of the production processes is well understood and therefore the production rate can be calculated for each point in the atmosphere, depending on the heliospheric modulation and the geomagnetic field intensity, provided the involved nuclear cross-sections are known [2]. As an example, Fig. 3 shows the dependence of the mean global production rate of 10Be as a function of solar modulation parameter F (F = 0: quiet sun, F = 1’000: very active sun) and the geomagnetic field intensity B in relative units (B = 1 corresponds to the present field intensity). As can be seen, the dynamic range between no magnetic field (B = 0), no solar modulation (F = 0) and doubled magnetic field (B = 2), very active sun (F = 1’000) is about one order of magnitude. Note that the dependencies are non-linear and that production changes by only a factor 3-4 were observed so far.
4
3
2
Be Production Rate 1 10
Rel. 0 0.5 Geomagnetic Field 1 0 200 1.5 400 600 800 2 1000 F) Solar activity (
Figure 3: Dependence of the relative mean global 10Be production rate on the geomagnetic field intensity and the solar activity parameter F. The production rate 1 corresponds to a geomagnetic field 1 and a F of 550 corresponding to the average solar activity. The transport of the cosmogenic nuclides produced in the atmosphere is not as well 14 understood as the production processes. C forms CO2 and exchanges between the main reservoirs of the carbon cycle (atmosphere, ocean, biosphere). 10Be and 36Cl become attached to aerosols or exist in gaseous form (H36Cl). After a mean residence time of 1 to 2 years they are removed from the atmosphere mainly by wet precipitation. In Polar Regions, the aerosols are removed by the snow that forms the ice sheet. Assuming a production rate of 0.018 10Be atoms cm-2 s-1 (Table 2) and a precipitation rate of 100 cm y-1, a simple calculation reveals an average 10Be concentration of approximately 107 atoms per kg of ice. Extremely sensitive detection techniques are necessary to measure 107 atoms. Due to the long half- life, decay counting is not feasible. However, accelerator mass spectrometry (AMS), using single atom detection is suitable to do the job [3]. A known amount of stable 9Be (typically 0.5 mg) is added to each sample. This leads to a 10Be/9Be ratio in the range of 10-13 to 10-12. After extraction of the Be from the water by ion exchange technique, a BeO sample is produced. This sample is put into the ion source of the AMS system and an ion beam is produced and accelerated to high energy (20 MeV) by means of a tandem accelerator. This high energy destroys the molecular background and enables suppression of the isobaric background (10B in the case of 10Be). In the following, some of the results obtained so far are discussed:
3. GEOMAGNETIC MODULATION To reconstruct the geomagnetic field from the 10Be and 36Cl fluxes we assume that the 10Be and 36Cl fluxes at Summit are proportional to their average global production rate.
geomagneticK field intensity (deduced from Be-10 and Cl-36 data) 1.5 geomagnetic field intensity (Mediterranean Sea)
1
0.5 Geomagnetic field intensity [relative to the present level]
0 20 30 40 50 60 Age BP [kyrs]
Figure 4: Comparison between the geomagnetic field reconstructed from a combined 10Be - 36Cl record from the GRIP ice core [4] with the paleomagntic data derived from a mediterranean sediment core [5].
Figure 4 shows the geomagnetic field intensity for the period 20-60 kyr BP, reconstructed from the combined 10Be and 36Cl flux in the GRIP ice core [4]. The shaded area indicates the uncertainty in the calculated field intensity. Also shown is the field measured on a Mediterranean sediment core [5]. The correlation between the geomagnetic field intensities obtained from these two independent reconstructions is very high (r2=70%). In our calculation of the geomagnetic field intensity, the combined 10Be and 36Cl flux was normalized in such a way that a value of 10% of its current value is assumed for the minimum of the calculated geomagnetic field intensity at about 40 kyr BP (Laschamp event). The normalization is also supported by new data from sediment cores of the Atlantic ocean [6]. The Laschamp event corresponds to a period of increased cosmic ray flux and therefore provides a test case for the proposed relationship between cosmic ray flux, cloud cover, and climate change [20]. The maximum of the combined flux of 10Be and 36Cl should be correlated 18 with the d O data (note the inverse scale) and CH4 data (Fig. 5). However, this is clearly not the case. During the Laschamp event (36-41.5 kyr B.P.) the combined flux of 10Be and 36Cl is not 18 2 2 significantly (p < 0.1) correlated with either d O (r = 0.07%) or CH4 (r = 0.09%). The same applies over the entire time interval shown in Figure 5 (r2 = 0.3% and 0.4%, respectively) [7]. A 2 (a) (10Be + 36Cl) (proxy for 1.8 cosmic ray flux)
Cl)-flux 1.6 36 1.4 (rel. units) Be + 1.2 10 ( 1
(b) B 18O -42 d (proxy for climate) -41
-40
O [‰] -39
18 corr.(A,B): d -38 r2=0.2%
(c) 350 C CH4 (proxy for climate) 400
450 [ppbv] 4 500 corr.(A,C): CH 550 r2=0.4%
20 25 30 35 40 45 50 55 60
Age BP [kyr]
10 36 18 Figure 5: Comparison of the combined Be- Cl flux with the climate parameters d O and CH4. According to the proposed relationship between cosmic ray flux and climate (Svensmark, this volume) a correlation between the three parameters is expected for the Laschamp geomagnetic minimum (shaded area) which is not present [7].
4. SOLAR MODULATION Direct observations clearly reveal that part of the solar variability is cyclic. In the following, we will concentrate only on cycles with time scales of years and longer. Cycles with periodicities from centuries to millennia are based on indirect or proxy data. Since these data (e.g. 10Be, 14C) represent a complex combination of different signals it is not always possible to unambiguously attribute a cycle to solar variability. One way of distinguishing between solar variability induced signals and others is to compare 10Be and 14C. Both radionuclides are produced by similar nuclear reactions in the atmosphere. Their respective production rate and their dependence on solar activity can be calculated [2]. However, after production their geochemical behaviour differs completely. A comparison of the two radionuclide records therefore allows us to distinguish between the production signal caused by solar and geomagnetic modulation and the system signal caused by the climate affecting the transport and the exchange processes between the different reservoirs. The results from such comparisons indicate that, for the past several millennia, the short- 14 term (decades to centuries) fluctuations in the D C record are mainly due to production variations, most probably caused by solar modulation. It is important to note that cycles associated with solar activity do not have a fixed periodicity. For example in the case of the sunspot cycle, the periodicity varies between 9 and 17 years. This raises the important question whether the periodicity averaged over longer times remains constant or not [8, 9]. To answer this question, longer and very precisely dated records of solar activity are needed than are presently available. The most prominent solar cycle is the 11-y Schwabe cycle discovered by Schwabe in 1843 when analysing his 17 year-long sunspot data. In Fig. 6, the sunspot cycle based on sunspot groups [10] is shown for the period 1600-1999 together with the inversely plotted 10Be concentration measured in the Dye 3 ice core from South Greenland [9].
0.5 10 Be [10 4 g -1 ] 1
100 Dalton Maunder 50 Sunspots
0 1600 1650 1700 1750 1800 1850 1900 1950 2000 YEAR
Figure 6: Comparison of sunspot numbers with 10Be concentration. Periods of local reduced solar activity are dashed.
In view of the fact that sunspot numbers and heliospheric modulation of the 10Be production rate are different representations of a common cause, i.e. solar activity, the agreement is good. A detailed analysis shows that the 10Be signal lags behind the sunspot signal by about 1 year, corresponding to the mean residence of 10Be in the atmosphere. It is interesting to note that the Schwabe cycle is still present in the 10Be record during the Maunder minimum [11]. A 90-year cycle was discussed by Gleissberg when analysing the auroral record [12]. The Dye 3 annual 10Be record going back to 1423 also shows the 90-year Gleissberg cycle [13].
14 The 205-year DeVries cycle is the most prominent periodicity in the D C record during the Holocene. However, as with other periodicities, its amplitude and periodicity are variable with time. Since the sunspot record is too short to detect the 205-year DeVries cycle, its attribution to solar variability is based on indirect evidence. Cycles with longer periodicities (e.g. 1000-2000 years) could not yet be attributed to solar modulation. An especially interesting feature of the sunspot record is the period from 1645 to 1715 A.D. which is characterized by an almost complete absence of sunspots (Fig. 6), the so-called Maunder minimum. Since then, solar activity has steadily grown with the exceptions of a few less pronounced minima: the Dalton minimum (1790-1830) and some weaker minima around 1890 and 1960. 40
30
20
10 C [‰]
14 0 D
-10
-20 7000 6000 5000 4000 3000 2000 1000 0 Year BP
Figure 7: 14C peaks corresponding to periods of low solar activity and possibly also reduced solar irradiance
Maunder type minima occurred earlier throughout the Holocene and are called grand minima. Since only little is known about these grand minima from direct observations, their occurrence is documented mainly by cosmogenic nuclide records. Fig. 7 shows the detrended 14 D C record [14]: The grand minima that correspond to maxima with regard to the cosmogenic 14 nuclide production are marked with arrows. How do we know that these peaks in D C are of solar origin and not caused by climatic effects or geomagnetic modulation? Firstly, the similarity in amplitude and duration of the peaks with the one corresponding to the Maunder minimum points to a common cause. Secondly, it seems rather unlikely that the geomagnetic dipole field would exhibit such strong excursions within only approximately a century. Finally, the good agreement 14 14 10 between the measured D C and the calculated D C based on Be data from ice cores convincingly shows that these peaks are due to production and not climatic effects. This brings us to the last topic, solar forcing and its detection [15].
5. SOLAR FORCING OF CLIMATE CHANGE The two main problems related to solar forcing and climate change are: 1. The lack of a quantitative solar forcing function. The physical processes responsible for changes in solar irradiance are not yet well understood, especially as far as long-term changes are concerned. All attempts so far are therefore based mainly on various assumptions leading to differences of about a factor of 2. Longer forcing records are based on simple linear regression models [16]. There may be other effects on the atmosphere caused by the interaction of the heliosphere with the magnetosphere and by cosmic rays with the atmosphere which could also contribute to climate change [17]. 2. The response function of the climate system to solar forcing is probably variable in time and not well known. There is an increasing number of experiments with global circulation models (GCM) including solar forcing. However, these model runs do not take into account the change in the spectral energy distribution and its potential effects on the atmosphere (e.g. ozone). In view of all these uncertainties, which would be the best strategy to detect solar induced climate changes? One approach is to use the Milankovic forcing that is caused by planetary gravitational effects on the orbital parameters of the earth [18]. Although only changes of the eccentricity causes changes in the total solar irradiance, the fact that the latitudinal forcing function can be calculated precisely for any time offers the unique opportunity to study the response function on longer time-scales (≥ 10 ky). Another straightforward approach is to search for fingerprints of solar forcing. All we know for sure is that solar irradiance changed in phase with solar activity over the past two Schwabe cycles. It is reasonable to assume that longer solar activity changes are associated with larger changes in solar irradiance [15]. Therefore, good candidates for solar forcing effects are solar minima, in particular grand minima. In fact, instrumental temperature records reveal cold events during the local minima around 1810, 1890 and 1960 (Fig. 5). The Maunder and Spoerer minima occurred during the so-called “little ice-age”, a period characterized by a general advance of glaciers. The more high-resolution climate records become available, the more evidence is found that abrupt climate changes indeed often coincide with solar minima (Van Geel, this volume)[19]. With regard to the question of the underlying physical mechanisms of solar forcing, a crucial test is the phase relationship. While the proposed mechanism of cosmic ray induced cloud formation tolerates no phase shift between cosmic ray flux and climate response [20], this is less the case for changes in solar irradiance that may be coupled by slow processes with the modulation of cosmic rays.
6. CONCLUSIONS Ice cores contain a large number of proxies for different climate parameters such as for 18 temperature (d O), greenhouse gases (CO2, CH4) and aerosols (chemical constituents). In the form of cosmogenic nuclides (10Be, 36Cl) they also provide unique information about the cosmic ray flux which is modulated by the geomagnetic dipole field and the solar activity that can be traced back in time over the past ca. 60'000 years. The suggested relationship between geomagnetic field, galactic cosmic rays, and climate could not be confirmed for the period of the Laschamp event (36-41.5 kyr B.P.). 10Be measurements show that solar variability has a cyclic component with periodicities of 11, 90, 205 and possibly more years. However, the relationship between solar activity and solar irradiance is not yet understood in detail.
ACKNOWLEGMENTS The author thanks W. Mende, R. Muscheler and G. Wagner for helpful discussions and C. Wedema for typing the manuscript and improving the English. This work was financially supported by the Swiss National Science Foundation.
REFERENCES [1] Thompson, L.G., et al., Tropical climate instability: the Last Glacial cycle from a Qinghai- Tibetan ice core. Science, 276: p. 1821-1825 (1997) [2] Masarik, J. and J. Beer, Simulation of particle fluxes and cosmogenic nuclide production in the Earth's atmosphere. J. Geophys. Res., 104(D10): p. 12,099 - 13,012. (1999) [3] Suter, M., et al., Advances in AMS at Zurich. Nucl. Instrum. Meth., B40/41: p. 734-740. (1989) [4] Wagner, G., et al., Reconstruction of the geomagnetic field between 20 and 60 ka BP from cosmogenic radionuclides in the GRIP ice core. Nucl. Instr. Meth., B 172: p. 597-604.( 2000) [5] Tric, E., et al., Paleointensity of the geomagnetic field during the last 80,000 years. J. Geophys. Res., 97: p. 9337-9351.( 1992) [6] Laj, C., et al., North Atlantic paleointensity stack since 75 ka (NAPIS-75) and the duration of the Laschamp event. Phil. Trans. R. Soc. Lond., 358: p. 1009-1025 (2000) [7] Wagner, G., et al., Some results relevant to the discussion of a possible link between cosmic rays and the Earth’s climate. J. Geophys. Res., 106(D4): p. 3381-3388.( 2001) [8] Dicke, R.H., Is there a chronometer hidden deep in the Sun? Nature, 276: p. 676-680.( 1978) [9] Beer, J., et al., Solar Variability Traced by Cosmogenic Isotopes, in The Sun as a Variable Star: Solar and Stellar Irradiance Variations, J.M. Pap, et al., Cambridge University Press: Cambridge. p. 291-300.Editors. (1994) [10] Hoyt, D.V. and K.H. Schatten, Group sunspot numbers: a new solar activity reconstruction. Solar Physics, 179: p. 189-219.( 1998) [11] Beer, J., S.M. Tobias, and N.O. Weiss, An active Sun throughout the Maunder minimum. Solar Physics, 181: p. 237-249.( 1998) [12] Gleissberg, W., The eighty-year solar cycle in auroral frequency numbers. J. Br. Astron. Assoc., 75: p. 227.( 1965) [13] Beer, J., et al., 10Be as an indicator of solar variability and climate, in The solar engine and its influence on terrestrial atmosphere and climate, E. Nesme-Ribes, Springer-Verlag: Berlin. p. 221-233. Editor. (1994) [14] Stuiver, M., et al., INTCAL98 Radiocarbon age calibration, 24,000-0 cal BP. Radiocarbon, 40(3): p. 1041-1083.( 1998) [15] Beer, J., W. Mende, and R. Stellmacher, The role of the Sun in climate forcing. Quat. Sci. Rev., 19(1-5): p. 403-415.( 2000) [16] Bard, E., et al., Solar irradiance during the last 1200 years based on cosmogenic nuclides. Tellus, 52B: p. 985-992.( 2000) [17] Tinsley, B.A., G.M. Brown, and P.H. Scherrer, Solar Variability Influences on Weather and Climate: Possible Connections Through Cosmic Ray Fluxes and Storm Intensification. J. Geophys. Res., 94(D12): p. 14783-14792.( 1989) [18] Berger, A., et al., eds. Milankovitch and Climate. NATO ASI Series. Vol. 126, D. Reidel Publishing Company: Dordrecht. Vol 1 und 2 zusammen: 895.( 1984) [19] Magny, M., Solar influences on Holocene climatic changes illustrated by correlations between past lake-level fluctuations and the atmospheric 14C record. Quaternary Research, 40: p. 1-9.( 1993) [20] Svensmark, H. and E. Friis-Christensen, Variation of cosmic ray flux and global cloud coverage - a missing link in solar-climate relationships. J. Atm. Terr. Phys., 59(11): p. 1225-1232.( 1997) LONG-TERM VARIATIONS IN COSMIC RAY FLUXES AND TOTAL SOLAR IRRADIANCE AND THEIR POSSIBLE INFLUENCE ON GLOBAL CLIMATE CHANGE
M. Lockwood *) Space Science and Technology Department, Rutherford Appleton Laboratory, Oxfordshire, UK
Abstract Studies of how geomagnetic activity is excited by the solar wind flow have allowed quantification of the open magnetic flux of the Sun, revealing it to have more than doubled during the 20th century. This flux fills the heliosphere out to the termination shock and shields Earth from galactic cosmic rays: thus, were air ions produced by cosmic rays to facilitate the formation of clouds in any way, this magnetic field would modulate terrestrial cloud cover. We here confirm there is a strong and statistically significant anticorrelation between the heliospheric field and the global coverage of low-altitude (<3.2 km) clouds and discuss the implications for extrapolating cloud-cover estimates back in time. We also show that the correlation between clouds and cosmic rays and the anticorrelation between clouds and total solar irradiance (TSI) are very similar in their strength and significance, making distinction between potential TSI and cosmic ray effects difficult to achieve.
1. INTRODUCTION The aa index of geomagnetic activity was devised by Mayaud in 1972 and, on annual timescales at least, successfully quantifies global geomagnetic activity from just two, antipodal observatories [1]. The importance of this index lies in the fact that it is a homogenous data series that extends back to 1868. Lockwood et al. have recently used the aa data to infer long-term changes in the open flux of the Sun that threads the coronal source surface and is dragged into the heliosphere by the solar wind flow [2]. The method they devised was an inversion of the analysis of Stamper et al. [3] who used data from solar cycles 20, 21 and 22, for which regular spacecraft measurements of the near-Earth heliosphere are available. The method was refined by Lockwood and Stamper [4] who used only cycles 21 and 22 to determine the required coefficients and held back the heliospheric field measurements from cycle 20 as an independent test of the method. The RMS differences between the inferred and observed radial components of the heliospheric field Br for the test cycle 20 were actually smaller than for the fitted cycles 21 and 22. Further confirmation of the method comes from the very high and highly significant anticorrelation between the inferred open solar flux and the counts detected by various neutron monitors due to cosmic ray bombardment of the atmosphere [5]: 80% of the variation of the cosmic ray flux could be associated with the heliospheric field strength. A number of different processes contribute to the shielding of cosmic rays [6], but scattering by irregularities in the heliospheric field is a dominant effect [7], such that the shielding is dependent on the open solar flux. Included in the remaining 20% that is not explained by the variation in open solar flux, is the known effect of solar cycle number on cosmic ray fluxes at Earth. This is expected theoretically as a consequence of the gradient and curvature drifts associated with large-scale heliospheric structure [6]. The effect reverses with the polarity of the polar solar field, roughly 1 year after the
*) Also at Department of Physics and Astronomy, Southampton University, Southampton, UK. peak of each cycle, and is apparent at the data, predominantly at solar minimum when the heliospheric field is weakest [7, 8, 9]. The method used to derive the heliospheric field is based on the theory of solar wind Ð magnetosphere coupling by Vasyluinas et al. [10], and thus by extrapolating back in time to before cycle 21 we are assuming that no there is no additional unknown factor, not included in this theory, the behaviour of which is different on decadal and century timescales. Given the correlation obtained for cycles 21 and 22 is 0.97, to be relevant this factor would need to have introduced variability before the start of cycle 21, but not subsequently. An important point in this respect is that the correlation between the open solar flux inferred from geomagnetic activity and cosmic ray neutron products was equally high and significant for solar cycles 19, 20, 21 and 22. Thus the method has been confirmed by independent data from cycles 19 and 20 Ð despite the fact that cycle 19 was the largest solar activity cycle ever observed and that cycle 20 was surprisingly weak. The key finding from the method is that the heliospheric field, averaged over full solar cycles, increased by 140% over the 20th century [2]. Extrapolating using the correlations with the cosmic ray fluxes discussed above, yields that the flux of primary cosmic rays was some 15% higher on average in 1900 than at present for energies above 3 GeV and 4% higher for >13GeV [5]. Support for this inferred drift comes from the abundance of the 10Be found, for example, in the Dye-3 Greenland ice core [11, 12]. This is formed as a spallation product when cosmic rays impact O and N in the atmosphere, and is then deposited in the ice sheet by precipitation. The dependence of precipitation on climatic conditions introduces scatter, nevertheless a clear anticorrelation with the inferred open solar flux is found [5]. In addition, the variation of 14C found, for example, in tree rings is consistent with the change seen in 10Be [13]. The complication for both these isotopes is that the abundances detected are subject to climate change. However, the effect is very different in the two cases: 10Be is precipitated into ice sheets, a process that introduces lag and a dependence on climate, whereas 14C is directly absorbed in gaseous state but has reservoirs in the biomass and oceans, exchange with which masks the true cosmogenic production rate and is expected to vary with climate. However, the similarity of the inferred century-scale changes in 10Be and 14C production rates strongly implies that the cause is a variation in cosmic rays and not climatic. Svensmark and Friis-Christensen [14], Svensmark [15] and Marsh and Svensmark [16] have discussed a correlation between cosmic ray fluxes and global cloud cover on Earth. The correlation is best with higher energy cosmic ray fluxes and low-altitude cloud cover. The present paper contains three studies. Given that the open solar flux quantifies 80% of the cosmic ray variation, section 2 looks at the direct correlation between open solar flux and cloud cover. We then use the long-term variation in the open flux derived from the aa index to look at the possible change in global cloud cover since 1900, assuming the correlation were real and not influenced by any other factors. One possibility for such a factor is the total solar irradiance (TSI) of the Sun, which is now known to show a solar cycle variation [17, 18, 19] and which also shows an upward drift over the 20th century in a variety of reconstructions that employ proxy data [4, 20, 21, 22]. The open solar flux, for the interval of the global cloud cover data at least, is well correlated with the TSI [4, 23]. This correlation was originally found in annual mean data and before observations for the rising phase of solar cycle 23 became available [4]. However, recent work [23] has shown that the correlation, although somewhat lower in monthly averages (correlation coefficient, r = 0.61), is highly significant (>99.99%), and has been maintained in solar cycle 23. Section 3 compares the correlations between cosmic rays and cloud cover and TSI and cloud cover. Section 3 also considers the effect of temporal smoothing on the significance of these correlations. 2. CLOUD COVER AND OPEN SOLAR FLUX Figure 1 shows the time series of monthly means for 1984-1994 of the strength of the heliospheric field near Earth (|B|), the aa index and the percentage change in global cloud cover from the average,
Figure 1. Monthly means of : (a) the magnitude of the heliospheric field near Earth ( |B| , blue line) ; (b) the aa index (green and black line); and (c) -
2.1 Correlation Analysis The correlogram for annual means of
It can be seen that a good match to ccf is again achieved. The peak anticorrelation is at a lag of Ð1 month, more consistent with a solar wind propagation delay. The peak is r = Ð0.53, a weaker correlation coefficient than in Figure 2, but a statistically much more significant result because there is much less “persistence” or “conservation” in the unsmoothed data. We can quantify the significance using the Student’s t-test, by making a correction to the number of degrees of freedom to allow for persistence [30]:
1/2 Ne = N (1 - a1)/(1 + a1) , t = |r| {(Ne - 2)/(1 - r)} , (1) where N is the number of samples, Ne is the effective number of samples, and a1 is the mean autocorrelation at lag 1 of the input and the output. The acf at lag 1 for
0.40 for monthly data, giving a1 = 0.505, the t statistic derived from (1) then yields a significance -4 of 99.98% (i.e. the probability that this correlation is a chance occurrence is just 2¥10 ). 2.2 Extrapolation back in Time Data on the IMF magnitude |B| only extends back to 1963 [26]. However, we can go further back in time if we use the coronal source flux, Fs, which has been estimated from the aa index [2] and is related to |B| by: 2 2 Fs = (1/2). 4pR1 |Br| = 2pR1 |B| cos( g ) (2) where R1 is 1AU (the mean Earth-Sun distance), Br is the radial component of the IMF and g is the IMF garden hose angle at Earth [2, 23]: on annual timescales g is almost constant [2, 3] and thus |B| and Fs are approximately linearly related.
Figure 3. Same as Figure 2, but for monthly mean data
The uncertainties inherent in the use of equation (2) have recently been analysed in detail and are of order ±20% in monthly data and ±10% in annual data [23]. The implications for the past behaviour of global cloud cover depend on whether the relationship implied by the above correlation is linear or not. Figure 4 shows the scatter plot for the annual means of
Thus the data strongly imply that global cloud cover was higher around 1900 than it is now. However, before we can say this with certainty and quantify the factor involved, we need to understand the physical and chemical mechanisms of the interaction and so we can understand the regression fit and know which is the most appropriate functional form to use.
3. CLOUD COVER AND TOTAL SOLAR IRRADIANCE The correlation of global, low-altitude cloud cover is significantly higher for cosmic rays than for the 10.7 cm radio flux from the Sun [14, 15, 16]. However, at this wavelength, the solar emission is far from representative of the IR, optical, and UV emissions that dominate energy input into the terrestrial climate system (rather, F10.7 is very closely related to sunspot activity and is most relevant to Earth’s upper atmosphere, the thermosphere). A sequence of total solar irradiance (TSI) values covering more than 2 solar cycles has been compiled by Fröhlich and co-workers [18, 19]. This compilation requires careful intercalibration of the various space-based radiometers used and allowance for their degradations with time and exposure. The data reveal a solar cycle variation with TSI being of order 0.1% greater at sunspot maximum than at sunspot minimum. Figures 7 and 8 show the scatter plots of the integrated, global, low-altitude cloud cover
Figure 5. Extrapolated low-altitude global cloud cover estimates for 1868-1998: (black) from a linear fit to observation; and (blue) from a square-law loss. The observed data are shown in red.
Figure 6. Detail from Figure 5 for 1980-2000. Figure 7. Scatter plot of monthly means of global cloud cover
Figure 8. Scatter plot of monthly means of
Figure 9 shows the time series of monthly means for the
Figure 10. Same as figure 7, for 12-point running means of the monthly data
Figures 10, 11 and 12 correspond to 7, 8 and 9, respectively, but are for 12-point running means of the monthly data. It can be seen that the long-term variations are very similar in the data and the scatter in the scatter plots has been greatly reduced. For
Figure 11. Same as figure 8, for 12-point running means of the monthly data
4. CONCLUSIONS AND DISCUSSION The fraction of the globe covered by low-altitude clouds has shown a solar cycle variation in recent data [14, 15, 16]. With little more than one cycle of heterogeneous data, we cannot be sure that this is truly an oscillation that will continue to match the solar cycle variations so well. The temporal correlation between the low-altitude global cloud cover and cosmic ray counts is highly significant in monthly data. Introducing smoothing increases the correlation coefficient magnitude and makes the time-series plots appear to be in greater agreement, but also removes the statistical significance from the correlation. An anti-correlation between the total solar irradiance (TSI) and global cloud cover is found to have the same strength and significance as the correlation with cosmic ray fluxes. Thus we cannot tell if it is more likely to be the cosmic rays or the TSI that are influencing cloud cover. In addition, because of the nature of all correlation studies, cannot we be sure that either TSI or cosmic rays have a causal effect on cloud cover. A parametric study using a global coupled ocean-atmosphere model is required to see if we would expect this anticorrelation between TSI and low-altitude cloud cover. Whether caused by TSI or comic ray variations, the presence of a solar cycle signal in global cloud cover would effectively be an amplification of the solar effect on Earth’s climate. Recent modelling using the UK Hadley Centre’s HAD3CM global coupled ocean-atmosphere model has pointed towards the presence of such an amplification, both when fitting the 11-year solar cycle variation in the amplitude of the global spatial pattern of average tropospheric temperatures, and when fitting the 150-year drift in the global average surface temperature [M.R. Allen et al., Private communication, 2000]. In both cases an amplification factor of about 2.5 was needed to gain the best fit, compared to the value from radiative forcing arguments. To within the 90% confidence level, this factor varied between 1 and 6.
Figure 12. Same as figure 8, for 12-point running means of the monthly data
Somewhat surprisingly, this amplification of the solar influence calls for an amplification factor for the man-made influences that go into the model [M.R. Allen, private communication, 1999]. This amplification factor is much smaller than for the solar effect (of order 1.1). The key reason for this behaviour is that the main drift in the longer-term TSI variation took place between 1900-1950. This also happened to coincide with a period of reduced volcanic activity (a global cooling phenomenon). On the other hand, anthropogenic greenhouse gasses have had their dominant effect in the past 30 years. Underestimating the solar effect early in the 20th century effectively causes the model to fit an anthropogenic effect that starts earlier but is less steep.
ACKNOWLEDGEMENTS The author thanks Nigel Marsh and Henrik Svensmark for the provision of the D2 cloud cover data and the World Data Centre system for the cosmic ray, heliospheric and geomagnetic data. He is also grateful to Myles Allen, Nigel Marsh, Henrik Svensmark, Jasper Kirkby, and many other scientists for valuable discussions. This work was funded by the UK Particle Physics and Astronomy Research Council.
REFERENCES [1] P.N. Mayaud, J. Geophys. Res., 77, 6870-6874, 1972. [2] M. Lockwood et al., Nature, 399, 437-439, 1999. [3] R. Stamper, et al., J. Geophys Res., 104, 28,325-28,342, 1999 [4] M Lockwood and R. Stamper, Geophys Res. Lett., 26, 2461-2464, 1999. [5] M. Lockwood, J. Geophys Res., 106, 16021-16038, 2001 [6] J.R. Jokipii in “The Sun in Time”, eds. C.P. Sonnet, M.S. Giampapa and M.S. Matthews, Univ. of Arizona Press, pp. 205-221, 1991 [7] H.V. Cane, Geophys. Res. Lett., 26, 565-568, 1999 [8] I.G. Usoskin, J. Geophys.Res., 103, 9567-9574, 1998. [9] H.S. Ahluwalia, J. Geophys. Res., 102, 24,229-24,236, 1997 [10] V.M. Vasyliunas, Planet Space Sci., 30, 359-365, 1982. [11] K.G. McCracken and F.B. McDonald, The long-term modulation of the galactic cosmic radiation, 1500-2000, in press, in Proc. 27th. Int. Cosmic Ray Conference, Hamburg, 2001 [12] J. Beer et al., Sol. Phys., 181, 237-249, 1998. [13] E. Bard, et al., Earth and Planet. Sci. Lett., 150, 453-462, 1997. [14] H. Svensmark, and E. Friis-Christensen, J. Atmos. Sol. Terr. Phys., 59, 1225-1232, 1997. [15] H. Svensmark, Phys. Rev. Lett., 81, 5027-5030, 1998. [16] N. Marsh, and H. Svensmark, Space Sci. Rev., 94, (1/2), 215-230, 2000. [17] R.C. Willson, Science, 277, 1963-1965, 1997. [18] C. Fröhlich, and J. Lean, Geophys. Res. Lett., 25, 4377-4380, 1998. [19] C. Fröhlich, Space Sci. Rev., 94, (1/2), 15-24, 2000. [20] J. Lean, et al., Geophys. Res. Lett., 22, 3195-3198, 1995. [21] S.K. Solanki and M. Fligge, Geophys. Res. Lett., 26, 2465-2468, 1999. [22] Hoyt, D., and K. Schatten, J. Geophys. Res., 98, 18,895-18,906, 1993. [23] M. Lockwood, An evaluation of the correlation between open solar flux and total solar irradiance, Astron and Astrophys., in press, 2001. [24] Y.-M. Wang, Geophys. Res. Lett., 27, 621-624, 2000. [25] W.B. Rossow, et al., International Satellite Cloud Climatology Project (ISCCP): Documentation of new datasets, WMO/TD 737, World Meteorol. Organ., Geneva, 1996. [26] D. A. Couzens and J. H. King, Interplanetary Medium Data Book - Supplement 3, National Space Science Data Center, Goddard Space Flight Center, Greenbelt, Maryland, USA, 1986. [27] M. Lockwood, The relationship between the near-Earth Interplanetary field and the coronal source flux: dependence on timescale, J. Geophys. Res., in press, 2001. [28] M.S. Potgieter, Adv. in Space Res, 16(9), 191-203, 1995. [29] A.C. Cummings et al., J. Geophys. Res., 99, 11,547-11,552, 1994. [30] Wilks, D.S., Statistical methods in the atmospheric sciences, Academic Press, San Diego, California, USA, 1995. EVIDENCE FROM THE PAST: SOLAR FORCING OF CLIMATE CHANGE BY WAY OF COSMIC RAYS AND/OR BY SOLAR UV?
Bas van Geel1, Hans Renssen2 and Johannes van der Plicht3 1 Institute for Biodiversity and Ecosystem Dynamics, Universiteit van Amsterdam, Kruislaan 318, 1098 SM Amsterdam, The Netherlands [email protected] 2 Institut d«Astronomie et de Géophysique G. Lemaître, Université Catholique de Louvain, 2 Chemin du Cyclotron, B-1348 Louvain-la-Neuve, Belgium [email protected] 3 Centre for Isotope Research, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands [email protected]
Abstract Major Holocene shifts to cool and wet climate types in the temperate zones correspond to suddenly increasing values of the atmospheric 14C content, suggesting a link between changing solar activity and climate change. In the temperate zones the transition from the Subboreal to the Subatlantic (ca 850 cal BC) represents a sudden, strong shift from a relatively dry and warm climate to a humid and cool episode. The moment of change occurred at, or maybe even just before the start of a sharp rise of the atmospheric 14C content. In previous studies, we postulated two amplification mechanisms: a) increased cosmic ray flux causes an increase in atmospheric 14C content, and also a climate shift, b) a decline of solar UV causes a reduced stratospheric ozone concentration, leading to climate change at the earth surface. Two phenomena indicate that mechanism a) is much less likely than mechanism b): 1) The enhancement of cosmic ray intensity to relatively high levels took place several decades after the climate shift. 2) In Central Africa and in Western India there was a shift to dryness. Chronological differentiation in solar output may play a role, but this is purely hypothetical.
1. INTRODUCTION Over the last few hundred years, changes in solar irradiance have been relatively small (less than 1 W/m2). As a consequence, solar forcing of abrupt climate change has been controversial [1]. However, there is strong evidence from the past for an important role of the sun upon climate change [2-5]. To explain this past evidence of solar forcing, we postulated two possible amplifying mechanisms that could explain how relatively small changes in solar irradiance could lead to abrupt climatic shifts [6]. a) Changes of cosmic ray intensity (modulated by fluctuating solar wind) might have an effect on cloud formation and thus on the planetary albedo and on temperature [7], and/or b) Within the small changes of solar activity, changes in UV are important [1]. Changes in solar UV have an effect on ozone formation in the lower stratosphere. Variations in the ozone concentration modulate the stratospheric temperature, leading to changes in the stratospheric circulation that could be propagated downwards to the Earth’s surface, thus influencing atmospheric circulation patterns world-wide [8, 9]. We review the evidence for solar forcing of climate change at the Subboreal-Subatlantic transition, as found in raised bogs and other paleodata, and we evaluate the possible contribution of both mechanisms mentioned above.
2. RAISED BOG AS ARCHIVES OF PAST CLIMATE Peat deposits are valuable archives for paleoclimate studies. The so-called raised bogs in NW- Europe are rainwater fed and the paleohydrological changes of such bogs mainly reflect climate shifts. A climate shift around 850 calendar years BC is visible in raised bog profiles as a transition from peat which was formed during a period of a relatively warm climate (darker, more decomposed peat) to lighter coloured upper peat, formed during a period of cooler, wetter climatic conditions. We use the Radiocarbon (14C) method for precise dating of climate-induced transitions in peat layers. Radiocarbon ages are expressed in BP, Radiocarbon "years" relative to 1950 AD. Radiocarbon years are different from calendar years because the production of 14C has not been constant in the past due to changes in both the geomagnetic field strength and in solar activity. The 14C time scale is calibrated by measuring the 14C content of tree rings, dated absolutely by means of dendrochronology [10]. The solar activity changes characterise the calibration curve by means of fluctuations (the so-called "wiggles"). Calibration of a single Radiocarbon date usually yields an irregular probability distribution in calendar age, quite often over a long time interval. This is problematic in paleoclimatological studies, especially when a precise temporal comparison between different climate proxies is required. However, a sequence of (uncalibrated) 14C dates can be matched to the wiggles in the calibration curve (wiggle-match dating [11, 12]). A high-resolution 14C sample sequence can result in a precise chronology of the peat core. This dating strategy also revealed relationships between atmospheric 14C variations and short-term climatic fluctuations (as detected in peat deposits) caused by solar variations. Data from Holocene lake deposits in the Jura Mountains also strongly point to a relationship between 14C fluctuations and paleohydrological shifts under the influence of climate change [2]. The climate shift around 850 cal BC (Subboreal-Subatlantic transition) was one of the most important climate shifts during the Holocene. We focused on this transition, which was immediately followed by a sharp rise of the atmospheric 14C content during the period between 850 and 760 cal BC. We identified the peat-forming mosses (representatives of the genus Sphagnum) in peat profiles of Northwest and Central European raised bogs. Knowing the ecological preferences of the mosses, we could interpret the recorded changes in species composition in terms of hydrological changes, related to climate change [13, 14]. Before the climate shift from the Subboreal to the Subatlantic period, Sphagna of the section Acutifolia were important peat formers in the Dutch bogs. Then Sphagnum papillosum and Sphagnum imbricatum took over. This change of the peat-forming plants indicates a shift from relatively warm, to cooler, wetter climatic conditions. The paleo-record from raised bogs shows that the abrupt climate shift happened at, or maybe even shortly before, the start of the period of the sharply rising atmospheric 14C content (Figure 1). In various lowland regions in the Netherlands where settlement sites were present, the climate shift at the Subboreal-Subatlantic transition caused a considerable rise of the ground water table so that arable land was transformed into wetland, where peat growth started. Bronze Age farmers living in such areas had to migrate because they could no longer produce enough food in their original settlement areas. Like the raised bog evidence, the archaeological evidence also points to a climate shift just preceding the enhanced cosmic ray intensity [13]. We also found strong evidence for climate change around 850 cal BC in other parts of the world [15, 16 and references therein]. In the temperate zones of Europe, North America and South America there is evidence for an equatorward shift of suddenly enhanced Westerlies, while the climate changed (cooling, higher effective precipitation). Fig. 1: The radiocarbon calibration curve (lower diagram) for the period between 1000 to 500 BC and 14 corresponding atmospheric fluctuations (D C, upper diagram). The moment of the climate shift, which precedes the rise of the atmospheric radiocarbon content, is indicated with an arrow. 3. THE CONTRIBUTION OF AMPLIFYING MECHANISMS The observed climate changes around 850 cal BC may have been caused by the lowering of solar irradiation through two amplifying factors, namely, (1) increased cosmic ray intensity stimulating
Polar Cell Polar Easterlies PFJ 60°N polar front Prevailing Ferrel Westerlies Cell 30°N STJ
N.E. Trades Hadley Cell
0° ITCZ
S.E. Trades
° 30 S STJ Prevailing Westerlies
60°S Polar PFJ Easterlies
Fig. 2A Simplified model of the tropospheric circulation (similar to the present situation) before the discussed climate change around 850 cal BC. ITCZ: Intertropical Convergence Zone; STJ: Subtropical Jet; PFT: Polar Front Jet.
Polar Polar Easterlies Cell 60°N polar front PFJ Prevailing Westerlies Ferrel Cell 30 N ° STJ N.E. Trades Hadley Cell 0 ° ITCZ
S.E. Trades STJ 30°S Prevailing Westerlies PFJ 60 S ° Polar Easterlies
Fig. 2B: As in Fig. 2A, but for the period directly after the climate change around 850 cal BC. Grey arrows denote changes which may be summarised as follows: equatorward shift of location of Jets, expansion of polar cells (i.e., cooling in mid-latitudes), relocation of mid-latitude storm tracks (regional increase in precipitation), and reduction of strength of Hadley Cells (i.e., drier conditions in the tropics). cloud formation and possibly also precipitation in certain regions, and (2) reduced solar UV intensity, causing a decline of stratospheric ozone production and cooling as a result of less absorption of sunlight. Figure 2 [after Ref. 6] shows the effect a decline of solar UV would have on the atmospheric circulation near the Earth’s surface [compare Refs. 8 and 9]: a decrease in the latitudinal extent of Hadley Cell circulation (weakening of the monsoon) may have occurred with concomitant equatorward relocation of mid-latitude storm tracks [see also Ref. 17]. This picture fits in with the paleoclimatological evidence from the northern and southern temperate zones (cooler, wetter) and the contemporaneous dryness crisis in Central Africa and Western India, which is evident from pollen records and archaeological evidence [13, 16]. The evidence during the Subboreal-Subatlantic transition strongly supports the "Haigh model" [8, 9] as an effective amplification mechanism for changes in solar activity. The combination of detailed paleoclimatological data from different parts of the world delivers circumstantial evidence for the suggestion that the UV-ozone mechanism had more effect on climate than the mechanism related to the increase of the cosmic ray intensity. In summary, we conclude that there is paleo-evidence for solar forcing of climate change around 850 cal BC. Of the two possible amplification mechanisms, the reduced UV scenario was most likely the effective one. There seem to be two arguments against an important role of cosmic rays (cloud formation) in relation to climate change: 1) Detailed series of radiocarbon dates from archaeological sites and raised bogs [11, 13, 14, 18] show that the abrupt climate change around 850 cal BC had already occurred when the cosmogenic isotope 14C only started to show an initially insignificant rise. This is supported by data for 10Be (another cosmogenic isotope, the production of which more directly reflects changing cosmic ray intensities than 14C), showing a corresponding and more or less contemporaneous rise as shown by 14C. For this event there might be a delay of approximately 10 years in the 14C rise only [J. Beer, pers. comm.; compare Ref. 19]. Consequently, the time- lag in the rise of the 14C content (compared to climate change) cannot be attributed to possible delaying processes related to the carbon cycle. In other words: the strong rise in cosmic ray intensity only followed climate change, and thus cannot have triggered the change [compare Ref. 12 for major climate shifts during the Little Ice Age in relation to similar increases of atmospheric Radiocarbon]. 2) The widespread dryness in the tropics (weaker monsoon in Central Africa and Western India) after 850 cal BC is not an effect that is expected to occur with enhanced cloud formation under the influence of increased cosmic ray intensity. However, the dryness in the tropics may not be inconsistent, as climatic teleconnections are not always straightforward (e.g., in the case of El Niño) and it could be that cooling in the mid-latitudes (where enhanced cloud formation due to increased cosmic ray intensities may be favoured) has resulted in drying in some regions in the tropics. On the other hand, it must be noted that the observed world-wide, but strongly contrasting changes in climate at the Subboreal-Subatlantic transition fit remarkably well in the model for an important role of solar UV (see Fig. 2). An important role for the reduced UV-scenario would raise one, yet unanswered, question: could a considerable decline of solar activity indeed have chronologically different phases (effective electromagnetic signal before magnetic signal; so first a UV decline with strong effects on climate, and later a more gradual decline of solar wind affecting the increased production of cosmogenic isotopes)? Solar physicists might be able to answer this question. Alternatively, detailed observations of variations in solar activity in the near future may reveal a solution to the question about which amplification mechanism plays a role in solar forcing of climate change.
ACKNOWLEDGEMENTS We thank Jürg Beer and Raimund Muscheler for critical reading of the manuscript and Dmitri Mauquoy for correction of the text. REFERENCES [1] D.V. Hoyt and K.H. Schatten, The role of the sun in climate change, Oxford, 1997 (Oxford University Press) [2] M. Magny, Quat. Res. 40 (1993) 1. [3] B. van Geel et al., Quat. Sc. Rev. 18 (1999) 331. [4] D.A. Hodell et al., Science 292 (2001) 1367. [5] U. Neff et al., Nature 411 (2001) 291. [6] B. van Geel and H. Renssen, In: Water, Environment and Society in Times of Climatic Change (Kluwer, Dordrecht, 1998) p. 21-41. [7] H. Svensmark and E. Friis-Christensen, J. Atm. Sol. Terr. Phys. 59 (1997) 1225. [8] J.D. Haigh, Nature 370 (1994) 544. [9] J.D. Haigh, Science 272 (1996) 981. [10] M. Stuiver et al., Radiocarbon 40 (1998) 1041. [11] M.R. Kilian et al., , 2000. Quat. Sc. Rev. 19 (2000) 1011. [12] D. Mauquoy et al., Evidence from North-West European bogs showing that Little Ice Age climatic changes were driven by changes in solar activity. Holocene 12: in press. [13] B. van Geel et al., Radiocarbon 40 (1998) 535. [14] A. Speranza et al., Quat. Sc. Rev. 19 (2000) 1589. [15] B. van Geel et al., 2000. Holocene 10 (2000) 659. [16] B. van Geel et al., 2001. In: Y. Yasuda and V. Shinde (Eds), Monsoon and Civilization, Extended Abstracts of the 2nd International Workshop of the Asian Lake Drilling Programme (Pune, India, 2001) p. 35. [17] D. Shindell et al., Science 284 (1999) 305. [18] B. van Geel et al., J. Quat. Sc. 11 (1996) 451. [19] R. Muscheler et al., Terra Nostra 3 (2001) 156. THE ROLE OF CLOUD COVER VARIATIONS ON THE SOLAR 13 ILLUMINATION SIGNAL RECORDED BY d C OF A SHALLOW WATER IONIAN SEA CORE (1147-1975 AD)
G.Cini Castagnoli, G. Bonino, D.Cane and C.Taricco Dipartimento di Fisica Generale dell'Università, Via P.Giuria 1, 10125 Torino, and Istituto di Cosmogeofisica del CNR, Corso Fiume 4, 10133 Torino, Italy
Abstract 13 We show the d C profile of Globigerinoides ruber, measured in the GT90/3 shallow-water Ionian sea core, dated with high accuracy (better than 1%) using radiometric and tephroanalysis methods,. It is commonly 13 accepted that d C variations in symbiotic foraminifera mainly record the effects of symbiont density and of photosynthetic activity, varying with ambient light level. The core, extracted from the Gallipoli platform, was sampled at contiguous steps of thickness 2.5 mm, corresponding to 3.87 13 years. The d C profile covers the period 1147-1975 AD. During the first seven centuries it appears fairly flat, while it shows a steep increase between 1760 and 1950 of ~0.3‰. The analysis of the time series performed using different methods shows a dominant decadal periodicity throughout the record. The 11-year component is identified at high significance level by Monte Carlo singular spectrum analysis (MC-SSA); the SSA-reconstructed-11-year component is in phase with the sunspot
solar cycle. The average amplitude of this component is A11y=0.04‰. 13 The modern d C increase (induced by a light level increase) of about 0.3 ‰ is concomitant with the decrease of the number of cloudy days per 13 year of about 11% at the site of the core deposition. If also the d C 11-y cycle has its origin in the modulation of cloudiness, the observed variation of 0.08 ‰ (peak-to-trough) requires an 11 y-cloud cover cycle (paced by the sun) of about 11%*0.08‰/0.3‰=3%. This is of the same order of the 11 y solar cloudiness cycle proposed by Svensmark and Friis-Christensen for the recent solar cycles, on a global scale (1980- 1995).
1. INTRODUCTION The carbon isotopic ratio 13C/12C in the shells of symbiont-bearing foraminifera is controlled by symbiont density and by their photosynthetic activity [1], i.e. by the primary productivity of the habitat. Provided the isotopic ratio of the bath is known at the time of the shell growth [2, 3], the isotopic ratio can be utilised for the quantitative study of paleoceanic and paleoatmospheric processes. Since Urey proposal in 1947 [4], many isotopic measurements have been performed for elucidating the climatic changes in the geological past. But only few stable-isotope time series covering the last millennium are available, which may be used to determine recent-past changes. This happens in spite of their importance for understanding the evolution of the present-day environmental conditions. A few carbon stable isotope time series were studied in different archives, mainly corals (see e.g. Ref. [5]), but they cover only the last few centuries. Sediments with a high sedimentation rate, which allows high resolution, may offer the opportunity to study in detail the millennial time scale; however it is difficult to find a suitable site for an absolute dating. We have found the right characteristics in the Gallipoli Terrace (Gulf of Taranto, Ionian sea) at a water depth of about 200 m. The carbonatic sediment is deposited at a sedimentation rate constant over the last 2 millennia [6-8].
13 In this paper we present the d C signal in specimens of the symbiotic planktonic foraminifera Globigerinoides ruber of the shallow water Ionian sea core GT90/3.
13 The d C time series, covering the last 828 years with the time resolution of 3.87 years, gives us the possibility to acquire information a) on the modulation of light at sea surface by the solar irradiance and by the cloud coverage in this region over the past millennium; b) on the presence of an 11-year signal most likely forced by the solar cycle and c) on the importance of this interdecadal variation of solar origin with respect to the variations of the trend.
2. THE IONIAN SEDIMENTS
The coast and the whole Salentina Peninsula are very flat and there is no direct river discharge on the platform where we took the cores. We have extracted many cores in different coring campaigns. We performed an accurate dating by radiometric and tephroanalysis methods [6,9]. The sedimentation in the cores shows no obvious laminae or discontinuities; dating is based upon 210 226 evaluation of Pb (T1/2 = 22.3 years) "excess" with respect to the activity supported in situ by Ra 210 222 (T1/2=1600 years). The "excess" Pb is atmospheric fallout from decay of Rn (T 1/2 =3.82 days). Core dating by this method is restricted to ages not greater than 150 years. The high correlation coefficient between the profile of "excess" 210Pb in the sediment and a decreasing exponential provides information on the constancy of the sedimentation rate over the past two centuries. Checks on the Pb age profile and on its extrapolation to the whole core are obtained from a 137Cs spike at 1963-1964 AD, due to a peak in nuclear weapons testing, on the one hand, and from tephroanalysis, on the other. The latter identifies clinopyroxene sedimentation peaks corresponding to the well-known historical volcanic eruptions at Pompei (79 AD), Pollena (472 AD), Ischia (1301 AD), Monte Nuovo (1538 AD) and starting from 1631 AD up to the present identifies the minor peaks corresponding to the detailed registration of the volcanic activity of the Vesuvius by the Vesuvian Observatory [10]. The position of the Gallipoli Terrace is particularly favourable for the collection, in the core mud, of the volcanic markers, fallout of the Campanian area activity, because the westerly winds bring the ashes towards the Gulf of Taranto. The sedimentation rate s was found to be quite constant along the cores and uniform throughout the -1 whole platform in the last millennia: we determined s = (0.0645±0.0007) cm year ; therefore the core depth scale may be transformed in a time scale, accurate better than 1%. The volcanic markers allow also to infer that bioturbation is not effective in the region at least within the adopted sampling interval of 0.25 cm, corresponding to 3.87 years. In fact the number density of pyroxenes in the volcaniclastic layers, with sharp boundaries, is identical in different cores taken in the area (see Fig.2 in Ref.[6]). The presence of 137Cs at the proper core-depth guarantees that the top of the core has not been perturbed. In Fig.1 we show the carbonate profiles of the cores GT14, GT89/3 and GT90/3 (all from the same area), sampled at contiguous steps of the same thickness 0.25 cm to determine the total carbonate content (CaCO3). We may notice the remarkable correlation between the carbonate records of different cores, demonstrating the uniformity of the deposition of the platform. In the same figure, we present (at the base) the number density pyroxene record measured in the upper 130 cm of the cores. It provides exact time "benchmarks" starting from the first historical eruption of the Vesuvius, described by Plinius, which destroyed the cities of Pompei and Ercolano in 79 AD. We notice that we have found large peaks only in correspondence to the volcanic events historically recorded. In these sediments we have studied the profiles of different bulk properties of the mud, with the primary aim of providing time series useful for investigating solar-terrestrial relationships in the past millennia (see e.g. Ref.[11]). Recently, we have chosen to measure the stable isotope composition of G. ruber planktonic foraminiferal tests [12-14]. G. ruber, a surface-warm water- dwelling foraminifer, shows maximum abundance in the top 20 m of the mixed layer in early autumn when the thermocline begins to break down [15]. This symbiont bearing foraminifer and 13 13 therefore the d C of its shell, like the d C of other species [16], is mainly controlled by the symbiont photosynthetic activity and by the ambient irradiance levels.
sample number 1400 1300 1200 1100 1000 900 800 700 600 500 400 300 200 100 0 45 40 35 GT89-3 GT14
30 45 25 40 35 45 30 40
35 GT90-3 30 Ischia 1301 A.D. 25 200 Pompei Pollena Montenuovo 79A.D. 472 A.D. 1538 A.D. 100
0 0 -3500 -3000 -2500 -2000 -1500 -1000 -500 0 500 1000 1500 year (AD)
Fig.1. Carbonate profiles (percentage of CaCO3 in the sediment’s mud) of the three shallow-water Ionian sea cores GT14, GT89/3, GT90/3. We may read the CaCO3 concentration as a function of sample number (upper scale) and as a function of time (lower scale). The reference (top) level is 1979 AD. In the lower part of the figure, the pyroxenes profile is also plotted, clearly showing that the principal peaks of the last 2 millennia are those caused by the eruptions of Pompei, Pollena, Ischia and Montenuovo.
3. EXPERIMENTAL PROCEDURE
We sampled the GT90/3 core (39∞45'53"N, 17∞53'33"E, water depth 174 m) at 0.25 cm thickness intervals, from the top down to sample 215, in a continuous sequence, covering the time interval 1147 AD-1979 AD. Samples of about 5 g of sediment were soaked in 5% calgon solution over night then treated in 10% H2O2 to remove any residual organic material, subsequently washed with distilled water jet through a sieve (150 mm). The fraction > 150 mm was kept and oven-dried at 50∞C. G-ruber were picked up under microscope. For each sample, 20-30 specimens of G.ruber of the same size were selected for the isotopic measurements, which were performed using a VG- PRISM mass spectrometer fitted with an automated ISOCARB preparation device. Analytical precision based on internal standards was better than 0.1‰ [13]. Calibration of the mass Spectrometer to VPDB scale was done using NBS19 and NBS18 carbonate standards.
13 4. THE d C TIME SERIES AND ITS ANALYSIS 13 In Fig.2 we show the d C profile (mean value xm = 0.84‰, s = 0.14‰), consisting of a continuous record of N = 215 points from 1147 to 1975; the sampling interval is ts = 3.87 years. 1.3
1.1
0.9
0.7
0.5
0.3 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 years (AD)
13 Fig.2 d C profile measured in Globigerinoides Ruber of the shallow water Ionian sea core GT90/3 (red line). This profile covers the period 1147-1975 AD, with a resolution of 3.87 years.
The 5 points running average (heavy-black line) and the SSA-reconstructed trend from PCs 1-2-5 are superposed to the data. The most evident feature of the series is the rapid enrichment in 13 d C starting from about 1760. In order to obtain reliable results, the analysis of this time series was performed using different spectral methods, like periodogram, correlogram, maximum entropy method, superposition of epochs (SE) and singular spectrum analysis (SSA; see, e.g., the review paper of Ref.[17]; here we present the results obtained by the classical Blackman-Tukey (BT) correlogram ([18]; see also Ref.[19]) and by SSA [20-22]; the results are also tested using a Monte Carlo approach (MC-SSA) [23,24]. The typical problems of the classical spectral estimates (power leakage and high variance) has therefore to be tackle for the correlogram. In order to reduce the effect of power leakage (due to the implicit window of the time series) and to give a consistent estimate of the true spectrum, a Bartlett (triangular) window was applied to the autocorrelation function. We observe that different window (like, for example, the Hamming or Hanning window) give very similar spectral results. The variance reduction is obtained using a window length M 0.1 0.09 ~400 y 0.08 0.07 0.06 0.05 11.3 y 0.04 0.03 0.02 0.01 0 0 0.02 0.04 0.06 0.08 0.1 0.12 frequency (year-1) 13 Fig.3. Blackman-Tukey power spectrum of the d C time series. A Bartlett window (width 50) was applied to the autocorrelation function. The FFT spectrum was computed using 512 frequencies. Note the prominent peak at 11.3 years and the power at the low frequencies of the trend. 10 9 8 7 6 5 11 y (EOFs3-4) 4 3 2 1 0 010203040506070 eigenvalue order k Fig.4. SSA-variance spectrum, plotted as percent of the total variance (normalized eigenvalues lk) associated to each of the 65 eigenvectors EOFs of the CD matrix, in decreasing order of variance. EOFs 3-4, associated with the 11.3 years oscillation, carries ~9% of the total variance. The trend is represented by EOFs 1-2-5 (~20% of the total variance). In order to perform a good signal-to-noise separation, we use the Monte Carlo method (MC-SSA) [23,24]. In this approach, we assume a model for our time series (null-hypothesis) and we determine the parameters using a maximum-likelihood criterion. Then a Monte Carlo ensemble of surrogate series is generated from the model and SSA is applied to data and surrogates (EOFs of the null-hypothesis basis are used), in order to test whether it is possible to distinguish the series from the ensemble. Since a large class of geophysical processes generates series with larger power at lower frequencies, we have assumed AR(1) noise in evaluating evidence for trend and oscillations. This is done to avoid overestimation of the system predictability, having underestimated the amplitude of the stochastic component of the time series [24]. In our case, we adopt the size 10000 for the Monte Carlo ensembles and we assume at first a pure noise AR(1) null-hypothesis: in this case we note anomalous power at frequencies corresponding to periods of ª11 years and to the trend; this is significant at 98% level. We confidently reject this hypothesis since the eigenvalues corresponding to the 11-year oscillation and to the trend component stand above the Monte Carlo range; moreover the 11-year range is of interest a priori and we expect a low frequency enhancement due to the presence of an evident trend in the series. Now we assume AR(1)+trend (EOFs 1-2-5) null-hypothesis. The result of the test against this hypothesis is shown in Figure 5.1, where we have plotted the eigenvalues and the surrogate bars as a function of the dominant frequency associated with the corresponding EOFs of the composite null-hypothesis basis. The vertical bars indicate the range in which lie the 98% of the eigenvalues determined from the ensemble of Monte Carlo simulations. We note that EOFs with period of ª11 years still show more variance than expected on this null-hypothesis. This is significant at 98% level. Finally, we test against AR(1)+trend+11years (EOFs 3-4) null-hypothesis; since in Figure 5.2 there are no excursions above the 98% MC-bars, we cannot reject this hypothesis. 1 5.1) 99th & 1st percentiles 0.1 11 years trend 0.01 null-hypothesis AR(1)+EOFs 1-2-5 (trend) 0.001 0 0.1 0.2 0.3 0.4 0.5 1 5.2) 99th & 1st percentiles 0.1 11 years trend 0.01 null-hypothesis AR(1)+EOFs 1-2-5 (trend) + 3-4 (11 years) 0.001 0 0.1 0.2 0.3 0.4 0.5 frq. assoc. with EOF-k (cycles in 3.87 years) 13 Fig.5 Application of Monte Carlo SSA to the d C time series (the empty squares show the eigenvalues associated with the EOFs included in the null-hypothesis). The Monte Carlo ensemble size is 10000. 5.1) Test 13 of d C series against the composite null-hypothesis of AR(1) noise plus trend. EOFs corresponding to the 11.3- year oscillation show exceeding power with respect to this hypothesis. This is significant at the 98% level. 5.2) 13 Test of d C series against the AR(1) noise plus 11.3-year oscillation and trend. No excursions occur outside the 98% limits, than we cannot reject this null-hypothesis and therefore the series can be explained by this model. 13 Therefore, our spectral analysis suggests that the d C time series is composed of a trend, on which an oscillation of about 11 years is superimposed; the background consists of AR(1) noise. Applying the method of SE to the reconstructed components (principal components PCs 3- 4), we determine the average amplitude A11y=0.04‰ of the interdecadal signal over the 828 years covered by the series. 5. DISCUSSION 13 The connection of the d C-11-year component to the solar cycle may be inferred by inspection 13 of Fig.6, where the reconstructed d C component (from PCs 3-4) is compared with the sunspot number series from 1700 to 1975. 200 ~11 y cycle (RCs 3,4) 0.1 Sunspots 150 100 0.0 50 -0.1 0 1750 1800 1850 1900 1950 2000 year (AD) Fig.6. Comparison of the SSA-reconstructed component (~ 11.3-years, from PCs 3-4), heavy line, with the 13 sunspot number series, light line, from 1750 to 1975. The variations of d C and of sunspot number are in 13 phase, in such a way high solar activity corresponds to high d C values. 13 The two signals are in phase, suggesting a solar forcing of d C, in such a way that high solar 13 output corresponds to high d C values. This favours the point of view that higher solar irradiance, as seen in space at sunspot maximum in the last 20 years [25-26], induces higher primary productivity of the symbionts (algae living on the spines of the G.ruber), which, using preferentially 12C for growing their tissue, leave carbon enriched of 13C in the chamber of G.ruber, for the construction of the shell. The primary productivity could in principle be modified also by other factors (nutrients availability, etc); however changes in illumination seem to be the most effective on carbon isotope fractionation processes: in this case, one of the important sources of 13 d C variations, beside irradiance solar variability, are the changes from year to year of the cloud coverage. The cloudiness seems to be paced by the solar cycle through galactic cosmic ray (GCR) modulation, as suggested by Svensmark and Friis-Christensen [27] and Svensmark [28]. Higher solar activity corresponds to lower GCR flux, giving lower cloud coverage, thus reinforcing the effect of the enhanced irradiation. Unfortunately a quantitative evaluation of all the above mechanisms involved in the formation of the 13C signal in G.ruber is not yet available. A decrease in cloud coverage has been observed from 1875 to 1975 in the region in which we took the core. The meteorologist De Giorgi has collected an accurate homogeneous archive of rainfall data from 1875 to 1921 in Lecce (Puglia). In the paper of Mangia et al. [29] those data have been integrated up to 1980 by using data taken in Bari. A decrease of annual rainy days (with rainfall >0.2 cm) of about 11% between 1875 and 1980 has been reported. In figure 7 we 13 compare the d C data (dashed line) with the number of rainy days per year (solid line). 1.3 40 50 1.2 60 1.1 70 80 1 90 0.9 100 13C - GT90/3 d d d d 0.8 110 120 rainy days per year 0.7 130 0.6 140 1860 1880 1900 1920 1940 1960 1980 year (AD) 13 Fig.7. d C modern increase in GT90/3 core (dashed line) and decrease in rainy days per year (inverted scale, solid line) from 1875 to 1980. We note that the rainy days decrease of 40 days per year, i.e. of ~11%, is concomitant with 13 the d C modern increase of ~0.3‰. The Apulian observations are confirmed by other evidences. Buffoni et al. [30] have analysed the precipitation in Italy from 1833 to 1996. On a yearly basis a decreasing trend in precipitation of the order of 10% has been found statistically significant. Furthermore Russo et al. [31] have observed in the long running time series of rainy days over Genova a decrease of about 10% between 1833 and 1992 with a negative tendency throughout the whole period. Therefore the experimental decrease in our area of about 11% of rainy days between 1875 13 and 1980 can be used for interpreting the carbon isotope effect. If the d C modern increase of 13 0.3‰ is due to change in illumination produced by a cloud cover decrease of 11% and if the d C 11y cycle of 0.08‰ (peak-to-trough) is attributed to the same process, we deduce a cloud cover variation over the 11-y cycle of 0.08‰*11%/0.3‰=3%. This value is in agreement with the variation of the global cloud coverage, paced by the solar cycle, between 1980 and 1995, as given by Svensmark and Friis-Chrstensen [27] and Svensmark [28]. ACKNOWLEDGEMENTS We are grateful to Prof. Carlo Castagnoli for his support and discussions and to Alberto Romero for continuous technical assistance. This work was supported by MURST-Co-fin 98, 2000 and by CNR. REFERENCES [1] Spero, H.J. and Williams, D.F., 1988, Extracting environmental information from 13 planktonic foraminiferal d C data, Nature 335, 717-719. [2] Spero, H.J., Lerche, J. & Williams, D.F., 1991, Opening the carbon isotope "vital effect" black box, 2: Quantitative model for interpreting foraminiferal carbon isotope data. Paleoceanography 6, 639-655. [3] Spero, H.J., Bijma, J., Lea, D.W. and Bemis, B.E., 1997, Effect of seawater carbonate concentration on foraminiferal carbon and oxygen isotopes Nature 390, 497-500. [4] Urey, J.C., 1947, J.Chem.Soc. 562. [5] Boiseau, M,, Ghil, M., and Juillet-Leclerc, A., 1999, Climatic trends and interdecadal variability from South-Central Pacific coral records, Geophys.Res.Lett. 26, 2281-2884. [6] Cini Castagnoli, G., Bonino, G., Caprioglio, F., Provenzale, A., Serio, M. & Zhu, G.M., 1990, The carbonate profile of two recent Ionian sea cores: evidence that the sedimentation rate is constant over the last Millennia Geophys.Res.Lett. 17, 1937-1940. [7] Cini Castagnoli, G., Bonino, Provenzale, A., Serio, M. & Callegari, E., 1992a, The CaCO3 profiles of deep and shallow Mediterranean sea cores as indicators of past solar-terrestrial relationships, Nuovo Cimento C15, 547-563. [8] Cini Castagnoli, G., Bonino, G., Serio, M. & Sonett, C.P., 1992b, Common spectral features in the 5500-year record of total carbonate in sea sediments and Radiocarbon in tree rings. Radiocarbon 34, 798-805. [9] Bonino, G., Cini Castagnoli, G., Callegari E., and Zhu, G.M., 1993, Radiometric and tephroanalysis dating of Recent Ionian Sea cores, Nuovo Cimento C 16, 155-161. [10] Arno’, V., Principe, C., Rosi, M., Santacroce, R., Sbrana, A. and Sheridan M.F., 1987, Eruptive History, in Summa-Vesuvius, Quaderni de La Ricerca Scientifica, CNR, 114, (ed. Santacroce R.) 53-103. [11] Cini Castagnoli, G., 1997, The 200y variability in terrrestrial archives: solar forcing? in Past and Present variability of the Solar-Terrestrial System: measurements, data analysis and theoretical models (eds. Cini Castagnoli, G. and Provenzale, A. 47-57, I0S press, Amsterdam). [12] Cini Castagnoli, G., Bernasconi, S.M., Bonino, G., Della Monica, P., Taricco, C., 1999a, 700 year record of the 11 year solar cycle by planktonic foraminifera of a shallow water Mediterranean core, Adv. Space Res. 24, 233-236. [13] Cini Castagnoli, G., Bernasconi, S.M., Bonino, G., Della Monica, P., and Taricco, C., 18 1999b, Solar activity in the last millennium recorded in the d O profile of planktonic foraminifera of a shallow water Ionian sea core, Sol. Phys. 188, 191-202. [14] Cini Castagnoli, G., Bonino, G., Taricco, C., 2001, Long term solar-terrestrial records from sediments: carbon isotopes in planktonic foraminifera in the last millennium, Adv. Space Res., in press. [15] Beveridge, N.A.S., and Shackleton, N.J., 1994, Carbon isotopes in recent planktonic foraminifera: a record of anthropogenic CO2 invasion of the surface ocean, Earth and Planet. Sci. Lett. 126, 259-273. [16] Spero, H.J. and Lea, D.W., 1993, Intraspecific stable isotope variability in the planktonic foraminifera Globigerinoides sacculifer: Results from laboratory experiments. Mar. Micropaleontol. 22, 221-234.0 [17] Ghil, M., and Taricco, C., 1997, Advanced spectral-analysis methods, in Past and Present Variability of the Solar-terrestrial System: Measurement, Data Analysis and Theoretical Models Proc. of the Int. School of Physics E.Fermi, Varenna, 1996, Course CXXXIII, edited by G.Cini Castagnoli and A.Provenzale (IOS Press, Amsterdam), p.137-159. [18] Blackman, R.B., and Tukey, J.W., 1958, The Measurement of Power Spectra from the Point of View of Communication Engineering (Dover, New York). [19] Kay, S.M., 1988, Modern Spectral Estimation: Theory and Applications (Prentice-Hall, Cliffs, N.J.). [20] Broomhead, D.S., and King, G.P., 1986, Extracting qualitative dynamics from experimental data, Physica D 20, 217-236. [21] Vautard, R., and Ghil, M., 1989, Singular spectrum analysis in nonlinear dynamics, with applications to paleoclimatic time series, Physica D 35, 395-424. [22] Dettinger, M.D., Ghil, M., Strong, C.M., Weibel, W., and Yiou, P., 1995, Software expedites singular-spectrum analysis of noisy time serie, EOS Trans. AGU, 76, 12. [23] Allen, M.R., 1992, Interactions between the atmosphere and oceans on timescales of weeks to years, PhD Thesis, Clarendon Laboratory, Oxford. [24] Allen, M.R., and Smith, L.A., 1996, Monte Carlo SSA: detecting irregular oscillations in the presence of coloured noise, J.Clim.9, 3373. [25] Willson, R.C., and Hudson, H.S., 1991, The Sun’s luminosity over a complete solar cycle, Nature 351, 42-44. [26] Pap, J.M., 1997, in Past and Present Variability of the Solar-terrestrial System: Measurement, Data Analysis and Theoretical Models, Proc. of the Int. School of Physics E.Fermi, Varenna, 1996, Course CXXXIII, edited by G.Cini Castagnoli and A.Provenzale (IOS Press, Amsterdam), p.1-24. [27] Svensmark, H., and Friis-Christensen, E., 1997, Variation of cosmic ray flux and global cloud coverage- a missing link in solar-climate relationships, J.Atmos.Sol.Terr.Phys. 59, 1225-1232. [28] Svensmark, H., 1998, Influence of cosmic rays on Earthís climate, Phys.Rev.Lett., 81, 5027-5030. [29] Mangia, C., De Sanctis, L.V., Ruggiero, L., Zito, G., Zuanni, F., 1991, Un secolo di precipitazioni piovose a Lecce, 4th Workshop Progetto strategico: Clima, ambiente e territorio nel mezzogiorno, ISIATA-CNR, Lecce, 153-161. [30] Buffoni, L., Maugeri, M., Nanni, T., 1999, Precipitation in Italy from 1833 to 1996, Teor. Appl. Climatol. 63, 33-40. [31] G.Russo, C.Eva, C.Palau, A.Caneva and A.Sacchini, 2000, Nuovo Cimento, C23, 39-52. [32] Vautard, R., Yiou, P., and Ghil, M., 1992, Singular Spectrum Analysis: a toolkit for short, noisy, chaotic time series, Physica D58, 95-126. APPENDIX 1 The SSA approach involves 3 basic steps: 1) embedding the time series in a vector space of dimension M (for the choice of M, see Ref.[32]); 2) computing the MxM lag- covariance matrix CD of the data (see the two different approach of Broomhead and King [20] and of Vautard and Ghil, [21]); 3) diagonalising CD : T LD = EDCDED, where LD = diag(l1,l2,...,lM ), with l1 ≥ l2 ≥...≥ lM ≥ 0 and ED is the MxM matrix having the corresponding eigenvectors Ek , k=1,...,M as its columns. For each Ek we construct the time series (of length N-M+1), called k-th principal component (PC), representing the projection of the original time series in the direction determined by the eigenvector Ek (also called empirical orthogonal function, EOF). Each eigenvalue l k gives the variance of the corresponding PC; its square root is called singular value (SV). Having chosen a subset of eigenvalues, it is possible to extract time series of length N, combining PCs; these time series are called reconstructed components (RCs) and they can be superimposed to the original signal. COSMIC RAY MEASUREMENTS IN THE ATMOSPHERE Y.I. Stozhkov, N.S. Svirzhevsky, and V.S. Makhmutov Lebedev Physical Institute, Russian Academy of Sciences, Moscow, Russia Abstract We present the main characteristics of cosmic ray population in the atmosphere and its variability (11-and 22-year solar cycle variations, solar protons originating from powerful solar flares, energetic electron precipitation during geomagnetic disturbances and Forbush decreases of cosmic rays). The experimental data were obtained from the long-term cosmic ray monitoring in the atmosphere from 1957 to now. The relationship between cosmic ray fluxes, and atmospheric processes are also discussed. 1. METHOD OF REGULAR MEASUREMENTS OF COSMIC RAYS IN THE ATMOSPHERE Regular cosmic ray measurements of cosmic rays are carried out with different instruments: ionization chambers, meson telescopes, and neutron monitors on the ground. The idea of cosmic ray monitoring in the atmosphere with radiosondes was suggested by Prof. S.N. Vernov in the middle of fifties and was realized by him and Prof. A.N. Charakhchyan in 1957. Now the cosmic ray monitoring cover a wide range of cosmic ray energy spectrum. It is schematically shown in Fig. 1. 5 SCR 4 GCR 3 atm. cosmic ray flux 2 nm 1 7891011 12 lg E, eV Fig. 1 Schematic view of galactic and solar cosmic ray spectra (GCR, SCR, thick and thin curves accordingly). The dotted vertical lines show the minimal energy of primary particles, which are detected by radiosondes in the atmosphere (E>0.1 GeV, upper arrow labeled atm) and by neutron monitors on the ground level (E>1.5 GeV, arrow with nm). The ground-based ionization chambers and meson telescopes record the primary particles with E>9 GeV. The cosmic ray measurements in the atmosphere are made with standard radiosondes in which the charged particle detectors are Geiger counter (hereafter counter) and telescope consisting of two counters and with 7 mm Al plate between them. Single counter records charged particles (electrons with energy E>0.2 MeV, protons with E>5 MeV) and g-rays with E>0.02 MeV (efficiency <1 %). Telescope records electrons with E>5 MeV, protons with E>30 MeV and is not sensitive to g-rays. For the isotropic angular distribution of particles in the upper hemisphere the geometrical factors of these detectors are 15.1 cm2 and 17.8 cm2 sr. The long-term cosmic ray measurements in the atmosphere have been started at the several latitudes with the different geomagnetic cutoff rigidities Rc in the middle of the last century and they are continued till now [1]. Every day or several times per week balloon flights have been made. Also several sea expeditions had been organized where the measurements of cosmic ray fluxes in the atmosphere in a wide range of Rc had been made. Till now more than 70.000 balloon flights have been performed. In Table 1 the sites and periods of observations are given. The cosmic ray fluxes are measured from the ground level up to 30-35 km. Table 1. The sites and periods of cosmic ray measurements in the atmosphere. Site of observations Geographic Rc, GV Period of observations coordinates Mirny, Antarctica 66∞33¢ S; 93∞00¢ 0.04 03.1963 - present time Tixie 71∞33¢ N; 128∞54¢ 0.4 02.1978 Ð 09.1987 Murmansk region 68∞59¢ N; 33∞05¢ 0.6 07.1957 - present time Norilsk 69∞00¢ N; 88∞00¢ 0.6 11.1974 Ð 06.1982 Moscow region 55∞28¢ N; 37∞19¢ 2.4 07.1957 - present time Alma-Ata 43∞12¢ N; 76∞56¢ 6.7 03.1962 Ð 02.1992 Erevan 40∞10¢ N; 44∞30¢ 7.6 01.1976 Ð 06.1989 Sea expeditions 60∞ N - 60∞ S 0.1-17 1963-65; 1968-72; 1975-76; 1986-87 In the atmosphere the main part of charged particles is secondary ones except of altitudes h≥20 km in polar regions where there are low energy primary protons. Below h<20 km cosmic rays mainly consist of secondary electrons and muons. 2. GALACTIC COSMIC RAYS To study galactic cosmic ray flux variations in different energy intervals the magnetic field of the Earth and the atmosphere are used as separators of charged particles according to their rigidity and energy. As an example in Fig. 2 the data obtained during the flights of radiosondes at Murmansk, Mirny, and Moscow on September 1997 are presented. During 1997 solar activity level was low and cosmic ray fluxes in the atmosphere were maximal ones. In Fig. 3 and 4 the samples of data obtained at the northern and southern latitudes with the different values of Rc during the Antarctic sea expedition of 1986-87 are shown [2]. The several radiosonde launchings were made at the each latitude and averaged data are presented in these figures. One can see a noticeable dependence of cosmic ray fluxes on Rc. Also the atmospheric depth (or pressure) X where maximum fluxes of charged particles, Nm, are observed increases with the growth of Rc. The examples of the time dependencies of charged particle fluxes (monthly averaged values) measured at the polar (northern and southern) and middle latitudes in the stratosphere and troposphere are given in Fig. 5 and 6. The period of observations covers ~19-23 solar activity cycles. 40 35 30 25 20 h, km 15 10 5 0 0 1000 2000 3000 4000 N, min-1 Fig. 2. The count rate of single counter vs. altitude in the atmosphere: at the northern polar latitude, Murmansk region, Rc=0.6 GV (the radiosonde flights on 2 and 4 September 1997 - open circles and black points, accordingly); at Mirny in the Antarctic, Rc=0.04 GV (the flights on 3 and 8 September - open and black triangles, accordingly); at the middle latitude, Moscow region, Rc=2.4 GV (the flight on 3 September - open squares). 3.5 3.0 0.6 north 2.5 -1 2.4 s 2.0 -2 1.5 5.3 N, cm 1.0 6.7 0.5 10.7 13.7 0.0 110100 1000 X, g cm-2 Fig. 3. The cosmic ray fluxes vs. atmospheric pressure X measured at different latitudes of the northern hemisphere during the sea expedition of 1986-1987. The values of Rc in GV are shown for each curve. The vertical bars show standard errors. 3.5 3.0 0.04 south 2.5 1.7 -1 s 2.0 2.4 -2 1.5 3.4 N, cm 5.4 1.0 7.3 0.5 10 13.7 0.0 110100 1000 X, g cm-2 Fig. 4. The cosmic ray fluxes vs. atmospheric pressure X measured at different latitudes of the southern hemisphere during the sea expedition of 1986-1987. The values of Rc in GV are shown for each curve. The vertical bars show standard errors. From Fig. 6 it is seen that the cosmic ray latitude effect between polar and middle latitudes disappears in the troposphere that is the cosmic ray fluxes observed at these latitudes are equal. 3.5 3.0 -1 2.5 s -2 2.0 N, cm 1.5 1.0 55 65 75 85 95 Year (after 1900) Fig. 5. Time dependence of monthly averaged cosmic ray fluxes in the stratosphere at h=31 km (X=10 g/cm2) measured at the northern and southern polar latitudes (Rc=0.6 and 0.04 GV, upper solid and dotted curves, accordingly) and at the middle latitude (bottom gray curve, Rc=2.4 GV). 1.3 -1 1.1 s -2 N, cm 0.9 0.7 55 65 75 85 95 Year (after 1900) Fig. 6. Time dependence of cosmic ray fluxes averaged per month in the troposphere at h=10.5 km (X=250 g/cm2) measured at the northern polar latitudes (Rc=0.6 GV, solid curve) and at the middle latitude (dotted curve, Rc=2.4 GV). 2 At each atmospheric pressure level, X, g/cm , only particles with the energy E>Ea (or rigidity R>Ra) where Ea is the atmospheric cutoff energy can contribute to the count rate of our detectors. The atmospheric cutoff Ea or Ra is defined by the characteristics of nuclear interactions of primary cosmic rays with air atoms. From the latitude measurements (at the different Rc) one can get the values of atmospheric cutoff as a function of X. In Fig. 7 the relationship of Ra. and X is presented. 10 8 6 , GV a R 4 2 0 10 100 1000 X, g/cm2 Fig. 7. The atmospheric cutoff Ra. vs. atmospheric pressure X. Open points were obtained from the sea expedition data and black points Ð from the long-term data obtained at the stationary sites (see Table 1 and Fig, 3, 4). Solid curve is the approximation: 0.8 2 Ra.=0.04X where R is in GV and X is in g/cm . Thus, the measurements of cosmic ray fluxes at the different atmospheric depths give the information on the variations of primary cosmic ray integral fluxes from R>0.5 GV at Xª10 g/cm2 up to R>(9-10) GV at sea level. R =0.04 and 0.6 G 3.5 c 3.0 -1 s -2 , cm 2.5 max N Rc=2.4 GV 2.0 _ _ + + + 1.5 55 60 65 70 75 80 85 90 95 100 Year (after 190 Fig.8. Monthly values of Nmax (maximum cosmic ray fluxes in the atmosphere) recorded at the Murmansk region, Rc=0.6 GV (upper solid curve), at Mirny station, Antarctic, Rc=0.04 GV (upper gray curve) and at Moscow region, Rc=2.4 GV (bottom curve). The narrow vertical stripes between two dotted lines show the periods of solar polar magnetic field inversions and +/- signs denote the magnetic field polarity of the solar northern polar region. In Fig. 8 the long-term experimental data on maximum fluxes of cosmic rays in the atmosphere obtained at the latitudes with Rc=0.04, 0.6, and 2.4 GV are depicted. In this figure and Figs. 5 and 6 the 11-year changes of cosmic ray fluxes are seen: in 1965, 1977, 1987, and 1997 the Nmax values are maximum ones and in 1957, 1970, 1982, 1991 they are minimal. The peaks observed at the northern polar and middle latitudes in 1962-63 were due to the radioactivity from the nuclear explosions in the atmosphere. In Fig. 8 (also in Figs. 5 and 6) the 22-year solar magnetic cycle is seen in the time dependence of Nmax: during the negative phases of solar magnetic cycles (~1960-70 and ~1980-90) cosmic ray time dependence shows a peaked form and it has a plateau during the positive phases (~1970-80 and ~1990- 2000). The difference in cosmic ray drift current directions in the heliosphere during positive and negative phases of the solar magnetic cycle explains the peaked and smoothed curves observed [3]. The amplitude of the 11-year solar cycle variations decreases with the growth of atmospheric pressure X. In Fig. 9 the changes of yearly averaged cosmic ray minimum fluxes, Nmin, observed in the 11- year solar cycles relative to 1965 are shown as a function of X: A(x)=[100(N65(x) - Nmin (x)]/N65(x), % . The N(x) value had a maximum in 1965 (see Fig. 8). The value A decreases with the increase of X and at X>600 g/cm2 becomes about 3 %. In June- August, 1991 the absolute minimum cosmic ray fluxes were recorded for the whole period of the observation from 1957 till present time (see squares in Fig.9).. 60 (7-9).59 9 (7-9).70 8 50 (10-12).82 (6-8).91 7 40 6 , 5 a 30 R 4 A, % 20 3 2 10 1 0 0 0 200 400 600 800 X, g/cm2 Fig. 9. The amplitude of 11-year cosmic ray changes relative to 1965 vs. atmospheric pressure X. The periods considered (months and year) correspond to the minimum cosmic ray fluxes and are given in the insert of this Figure. Cosmic ray fluxes were averaged for these periods. The atmospheric cutoff Ra is shown by solid line. -0.1 3.5 0 0.1 0.2 -1 3 s 0.3 -2 0.4 , cm 2.5 0.5 max solar activity N 0.6 2 0.7 0.8 0.9 1.5 55 65 75 85 95 Year (after 190 Fig. 10. Time dependence of maximum cosmic ray flux in the polar atmosphere, Nmax (black curve) and solar activity level (gray curve) defined as h/j where h is sunspot number and j is sunspot average helio-altitude. The monthly averaged data smoothed with the period of T=3 months were used. The values of cosmic ray fluxes N(x) in the heliosphere and, in turn, in the atmosphere are defined by solar activity level. The close relationship is observed between N(x) and solar activity parameter (h/j) where h is a sunspot group number, j is sunspot averaged helio-altitude [4]. This relationship is demonstrated in Fig. 10. The correlation coefficient is R(Nmax, h/j) = Ð0.83±0.03. The galactic cosmic ray modulation in the heliosphere is produced by magnetic irregularities of interplanetary magnetic field (IMF). In turn, the value of IMF and its irregularity density are defined by solar activity. The density of these irregularities increases with the growth of IMF strength. So, we can expect a relationship between IMF strength and cosmic ray fluxes in the heliosphere as well as in the atmosphere. This relationship is given in Fig. 11. The correlation coefficient R(Nmax, IMF)= -0.71±0.04. 4 3.5 5 -1 3 6 s -2 7 IMF, , cm 2.5 8 max N 9 2 10 1.5 11 55 65 75 85 95 Year (after 190 Fig.11. Time dependence of maximum cosmic ray flux in the polar atmosphere, Nmax (black curve) and IMF (gray curve). The monthly averaged data smoothed with the period of T=3 months were used. The data on IMF were taken from INTERNET: http://nssdc.gsfc.nasa.gov/omniweb/. 30 25 20 15 Start time, UT, date h, km GCR 10 0822 11.09.0 1053 11.09.0 5 1942 11.09.0 0326 11.10.0 0 0.1 1.0 10.0 100.0 N, cm-2 s-1 Fig. 12. Cosmic ray fluxes recorded in the atmosphere at Mirny in the Antarctic during the solar proton event on 9 November 2000. Solid curve is a charged particle flux background produced by galactic cosmic rays in the atmosphere (GCR). Different symbols show the data obtained during the different flights of radiosondes (the date and start time of launchings are given in the insert).. 3. SOLAR COSMIC RAYS, PRECIPITATION, AND RADIOACTIVITY Since the beginning of cosmic ray measurements in the atmosphere in 1957 several tens of solar proton events were recorded (e.g. [5]). As a rule such events are observed in the polar atmosphere where Rc are rather low (see Table 1). During solar proton events total cosmic ray flux at high altitudes in the atmosphere increases in several (sometimes in tens) times. As example, the solar proton fluxes generated by solar flare on November 9, 2000 and recorded in the Earth’s atmosphere are given in Fig. 12. Solar proton fluxes were observed in the atmosphere at h>17 km and their values increase with altitude. 100 1 3 -1 sr 10 -1 2 s -2 1 J(>E), cm 4 0.1 100 1000 E, MeV Fig. 13. The energy spectra of solar protons on 9 November 2000. The start times of radiosonde launchings and exponents g of solar proton energy spectra I(>E)~E-g are given below: 1 ÐDate - 11.09.00, Start time (UT) - 8:22, g=7.2; 2 Ð Date - 11.09.00, Start time (UT) - 10:53, g=5.8; 3 Ð Date - 11.09.00, Start time (UT) - 13:36, g=7.8; 4 Ð Date Ð 11.10.00, Start time (UT) - 19:42, g=4.9. From these data the fluxes and energy spectra of solar flare particles in the energy range of E=100- 500 MeV were obtained. They are depicted in Fig. 13. In this solar proton event the particles with E>500 MeV were not observed and this event was not recorded by neutron monitors. We note that the observed solar proton “soft” energy spectra could be due to the additional acceleration of particles in the interplanetary space by shock waves as it was happened in the past during the solar proton events in August 1972 [6]. 200 10 8 150 6 100 4 sunspot number 50 2 solar proton event number 0 0 1955 1965 1975 1985 1995 Year Fig. 14. The time dependences of yearly number of the solar proton events recorded in the atmosphere (black points and right axes) and yearly sunspot number (open points and left axes). For the 45-year period of cosmic ray monitoring in the atmosphere 105 solar proton events were recorded. In Fig. 14 the time dependence of the yearly number of these events and solar activity level (sunspot number) are shown. It is seen that the solar proton events mainly occur during the ascending and descending phases of solar activity. 35 35 30 30 25 25 20 20 h, km 15 h, km 15 10 10 Start time UT Date Start time UT Date 5 1250 05.03.0 5 1250 05.03.0 0930 05.05.0 0930 05.05.0 0 0 0 1234 0123 N, cm-2 s-1 N, cm-2 s-1 Fig. 15. Precipitation of high energy electrons into the northern polar atmosphere recorded by single counter on 5 March 2000 (left panel, black points). The background from galactic cosmic rays is shown by open points. The telescope recorded galactic cosmic ray background only (right panel). The inserts show the dates of radiosonde flights and launching times. During the geomagnetic disturbed periods in the polar atmosphere at high altitudes high-energy electron precipitation events are detected [7, 8, 9]. Our northern polar station in Murmansk region is near the polar oval (McIlwain’s parameter L=5.6) where the precipitations are observed rather often. In Fig. 15 the example of precipitation detected on November 2000 is shown. The counter records secondary g-rays produced by precipitating electrons in the atmosphere.. But at the same time the telescope records only the background from galactic cosmic rays and it allows to separate precipitation from solar proton events. It is significant that g-rays recorded in the atmosphere at h=25-35 km are produced by the precipitating electrons with E>several MeV. The time dependence of the yearly precipitation number and the sunspot number are given in Fig. 16. Into the precipitation data the corrections for patrol efficiency were introduced. The data obtained show that the precipitation take place most often during the descending phase of solar activity (remind that solar proton events are observed most often during the ascending and descending solar activity phases, see Fig. 14). This fact was established earlier in other papers [10]. The total number of precipitation recorded at the station in Murmansk region during 1957-2000 equals 549 events. For almost the same period (1963-2000) at the Antarctic station Mirny (66∞33¢ S; 93∞00¢; Rc=0.04 GV; Lª11) 10 precipitation events were recorded only. In Table 2 the yearly precipitation numbers for 20 - 23rd solar cycles are given. The last line in this Table includes the ascending phase of 23rd solar activity cycle only. The regular monitoring of charged particle fluxes in the atmosphere provides the prompt control of radiation conditions and allows to detect radiation clouds from nuclear explosions or nuclear plant failures. In Fig. 17 and 18 the observations of radioactive clouds in the polar northern atmosphere and over Moscow are shown as an examples. The excesses of particles over the cosmic ray background recorded by single counter were due to the radioactivity particles. The telescope data showed normal count rate from cosmic rays only. A powerful surface nuclear explosion was produced near lake Lobnor in China on the 14th of October 1970. In the atmosphere radioactive clouds were observed near Murmansk on October 25-26 (see Fig. 17) and on November 11-12, near Moscow and Alma-Ata on November 5 and 6. In the polar region it was observed in the altitude range hª15-25 km and charged particle flux was increased in 7-8 times in comparison with the cosmic ray background. At the first registration the cloud had the shape of a disk with vertical size of nearly 4 km and horizontal size along the wind direction of ~1100 km. The maximal activity was equal ~10-3 Bq/cm3 as measured at hª18-22 km. After 26 October the radioactive cloud passed away from the observation site. A radioactive cloud in the atmosphere over Moscow was recorded on 12-14 April 1993 (see Fig. 18). The cloud was seen at hª10-30 km and had the horizontal extension of nearly 1000 km and maximum activity ~10-4 Bq/cm3. We do not know this cloud origin but on 6 April the failure at the chemical plant in Siberia (Seversk town near Tomsk) occurred [11]. A more dense cloud was observed over Moscow on the 27th of October 1999, the maximum activity being 7◊10-3 Bq/cm3 at h>15 km. The cloud horizontal extension was about 200 km. The source of this last cloud is unknown also. 100 200 80 150 60 100 40 50 sunspot number 20 precipitation number 0 0 1955 1965 1975 1985 1995 Year Fig. 16. The time variations of the yearly number of precipitation (black points) and sunspot number (open points). The observations were made at the northern polar station in Murmansk region (68∞59¢ N; 33∞05¢; Rc=0.6 GV; L=5.6) during the time interval of (6-12) UT. Table 2. The precipitation number in different solar activity cycles. solar cycle number, sunspot number precipitation Precipitation period per year number (total) number per year 20 (1964-1975) 58.8 144 12 21 (1976-1985) 82.9 140 14 22 (1986-1995) 78.5 118 12 23 (1996-2000) 46.9 85 17 30 25 20 15 h, km background 10 1602 LT 25 Oct. 197 1802 LT 25 Oct. 197 5 2207 LT 25 Oct. 197 1157 LT 26 Oct. 197 0 051015 20 N, cm-2 s-1 Fig. 17. Charged particle fluxes in the northern polar atmosphere detected by single counter. At h>15 km the excess of flux over the cosmic ray background (solid gray line) is due to radioactive cloud particles produced by the nuclear explosion in China on 14 October 1970. In the insert the legend on the balloon start time is given. 40 30 20 h, km 10 0 012345 N, cm-2 s-1 Fig. 18. The count rate of single counter vs. altitude in the atmosphere over Moscow on 12-14 April 1993:æ background from galactic cosmic rays; charge particle flux measurements on 13 April, launching local time 0830 LT (◊), on 13 April, 1430 LT (D); on 14 April, 0830 LT ( );. 4. COSMIC RAY FLUXES AND ATMOSPHERIC PROCESSES 10 If one compares the flux of solar electromagnetic radiation falling on the top of the atmosphere (Fsunª10 -2 -1 2 -2 -1 erg m s ) with the flux of cosmic ray energy (FCRª10 erg m s for particles with energy E≥0.1 GeV) the evident conclusion could be made: the influence of charged cosmic ray particles on the processes in the atmosphere is negligible in comparison with influence of the electromagnetic radiation coming from the Sun. However, let us imagine for a moment that cosmic rays stopped to intrude into the Earth’s atmosphere. The ion production will be aborted and the global electric circuit will be destroyed. The production of thundercloud electricity and lightning will be over. The cloud area will be decreased and precipitation level will fall down. The cosmic rays with energy E=0.1-15 GeV carry about 60 % of all cosmic ray energy and these particles constitute about 95 % of all cosmic ray flux. These particles undergo the influence of the geomagnetic field in such way that the fluxes of primary cosmic rays at polar latitudes is higher than the ones at equatorial regions as much as ~30-35 times. In the atmosphere this difference is about 4 times. Below some aspects of influence of charged particle fluxes on atmospheric processes are considered (see also [12]). In our analysis we use the experimental data of the long-term measurements of cosmic ray fluxes at the different atmospheric depths (from the Earth’s surface up to 30-35 km) and at the different latitudes. 4.1 The global electric circuit and ion production It is well known that the Earth has a negative electric charge about 6¥105 C and the strength of electric field produced by this charge near the Earth’s surface measured during fair-well weather equals Eª-130 V/m (directed to the Earth’s surface). The value of average current flowing between equalizing layer to be in the ionosphere at the altitude hª55-80 km and the Earth’s surface is Iª10-12A/m2 [13, 14]. The light ions provide this current in the atmosphere. The ions are produced by cosmic particles (radioactivity of soil also gives ions but only in the lower atmosphere at h< 3 km). If cosmic ray flux changes the ion density the air conductivity changes also. The lightning in thunderstorms and precipitation form another branch of the closed electric circuit charging the Earth by negative electricity and providing electric current from the Earth to the ionosphere. The sketch of global electric circuit is given in Fig. 19 (see, e.g. [16]). h =60 km Ip Rp Im, Rm Earth Ig Fig. 19. The sketch of the global electric circuit: h=60 km Ð equalizing layer; Ip, Rp and Im, Rm Ð atmospheric electric currents and resistances in the atmosphere at polar and middle latitudes, accordingly; Ig Ð current of thunderstorms and precipitation charging the Earth by negative electricity. The equation describing the relation between ion production rate, q, and their recombination in the atmosphere under quasi-state conditions is usually taken in the form q(h)=a(h) [n(h)]2, (1) where n is ion concentration, a is recombination coefficient, h is atmospheric altitude [17]. Using the experimental data on cosmic ray fluxes and ion concentrations in the atmosphere one can test the validity of this equation. In Fig. 20 the ion concentrations, n, and the charged particle fluxes, N, measured at several latitudes vs. altitude are presented [18]. From the experimental data on ion concentration n and cosmic ray flux N one can get that the ion production rate q is proportional to charged particle flux: q(h)=msN(h), where m and s are the number of air particles per cm3 and ionization cross-section. The values of m and s are the same for different latitudes and depends on the altitude only. It isn’t true for the case of polar latitudes and h>20 km where s is increased. At h<20 km the value of s equals 2*1018 cm2 within 10-15 % for all latitudes. 17.3 5.3 3 17.3 5.6 3.4 0.04 30 0.6 30 20 20 h, km h, km 10 10 0 0 0123 012 n, 103 cm-3 N, cm-2 s-1 Fig. 20. Ion concentration n (left panel) and cosmic ray flux N (right panel) as a functions of altitude h in the atmosphere at the latitudes with the geomagnetic cutoff rigidities Rc=17.3, 5.6, 5.3, 3.0, 3.4, 0.6 and 0.04 GV. Horizontal bars show the standard errors. Let us consider the measurements of n and N performed at two different latitudes. According to the expression (1), we can construct the following ratio: 2 2 [q1(h)/q2(h)]=[a(h) n (h)]1/[a(h) n (h)]2, (2) where the subscripts 1 and 2 correspond to the latitudes with different geomagnetic cutoff rigidities Rc1 and Rc2. Substituting msN(h) instead of q and taking (ms)1= (ms)2 and a1ªa2 (these suggestions are fulfilled in the atmosphere rather well) one can get 2 [N1(h)/N2(h)]=[n1(h)/n2(h)] . (3) 0.8 1 0.6 0.4 2 Ratio 0.2 3 0.0 51015202530 h, km Fig. 21. The ratio of cosmic ray fluxes (curve 1), ion concentrations (curve 2) and squared ion concentrations (curve 3) as a function of altitude. These values were calculated from the experimental data obtained at the equatorial (Rc=17.3 GV) and middle (Rc=3.3 GV) latitudes (see Fig.20) without any normalization of the data. The standard errors are given by vertical bars. In Fig. 21 the ratios of charged particle fluxes (curve 1-open points), ion concentrations (curve 2- dark points), and squared ion concentrations (curve 3-crosses) obtained from the experimental data presented in Fig. 20 are given. The data obtained in the equatorial (Rc=17.3 GV) and middle latitudes (Rc=3.3 GV) were used. It is seen that the ratio of cosmic ray fluxes (curve 1) coincides with ion concentration one (curve 2) and differs significantly from squared ion concentration ratio (curve 3). The details of such consideration are given in [18]. Thus, from this analysis the important conclusion must be made that the ion balance in the atmosphere under quiet conditions is described by linear equation (not quadratic one) q(h)=b(h) n(h), (4) where b(h) is the linear recombination coefficient. From the available experimental data on cosmic rays in the atmosphere and light ion concentrations the value of b(h) and q(h) can be calculated for any site of the Earth and any level of solar activity. The ion production rate q can be written as q(h)=N(h) s(h) r(h)/M, (5) where N(h) is cosmic ray flux at the altitude h, s is the ionization cross-section in air, r(h) is the air density and M is the average mass of air atom. The relationship between atmospheric electric current J, electric field strength E and conductivity l is J=l(h) E(h)=n(h) k(h) E(h), (6) where k(h) is the mobility of light ions at the altitude h. Thus, using the expressions 4, 5 and 6 one can find J=N(h) s(h) r(h) k(h) E(h)/[M b(h)]. (7) 2 0.66 N -2 1.6 -1 0.62 s A m -2 -12 J N, cm J, 10 1.2 0.58 0.8 0.54 1965 1970 1975 1980 1985 Year Fig. 22. The yearly average values of atmospheric electric current J(h) (from [19]) and cosmic ray flux N(h) at hª8 km in the polar region. On the right side of this equation all values are constant except cosmic ray flux N(h) and electric field strength E(h). If one supposes that E(h) is constant or weakly changes in the periods of fair-well weather then there is the linear relationship of cosmic ray flux N(h) and atmospheric electric current J(h). Such conclusion is confirmed by the experimental data showed in Fig.22. The data on J(h) were taken from [19]. The correlation coefficient between J(h) and N(h) is positive and equals r(J, N)= +0.77±0.10. The correlation of atmospheric electric current and solar activity level (sunspot number W) is low, r(J,W)=-0.32±0.22. 4.2 Thundercloud electricity and lightning production In 1920 Wilson put forward fascinating idea suggested that thunderstorms act as a global generator of electric current maintaining the Earth’s electric charge [20]. Since the experimental evidences supporting this hypothesis were obtained (see references in [14]). However, the mechanisms of thunderstorm electricity production (separation of negative and positive charges in thundercloud) and lightning generation are not clear till present time, although there are a number of hypotheses on the thundercloud electricity origin (see, e.g., [21-23]). Cosmic rays could be responsible for the thunderstorm electrification [24]. Secondary charged particle fluxes generated in the atmosphere by primary cosmic rays are the only source of positive and negative ion production in the atmosphere at h>(2-3) km. The problem consists in the spatial separation of negative and positive ions in the process of thundercloud formation. The thunderclouds are formed from ascending wet air mass when the fronts of cold and warm air meet each other. The air masses contain heavy ions (charged aerosols) because light ions produced by cosmic rays adhere to neutral heavy particles. As it is known from the observations the concentration of aerosols has a maximum in the low atmosphere near the Earth’s surface and its value is ~2¥104 cm-3. The half of these particles carries out the positive or negative electric charges [25]. Ascending air mass picks up the aerosols. During ascending air mass is cooled and processes of condensation of water molecules on neutral and charged aerosols take place. The condensation rate depends essentially on the charge presence and its sign. Namely, negative charged aerosols grow faster than positive ones as much as ~104 times [26, 27]. The rapid growth of aerosols with negative charge makes them heavy and their lift with the rising air mass is stopped at the low altitudes. At the same time aerosols with positive charge continue to rise with ascending wet air mass and stop their rising at higher altitudes than negative charged aerosols. In this way the spatial separation of electric charges inside the cloud occurs (in detail see [24]). Inside the thundercloud the strength of electric field can grow up to Eª3 kV/cm and the distance between separated positive and negative charges is roughly estimated as Dhª3-4 km. The high value of E is observed under thundercloud also. But the observed values of E are much less than the puncture voltage at the altitudes where thunderclouds exist (hª2-7 km). At hª3 km the value of puncture voltage is 15-30 kV/cm [28]. In [29] Ermakov put forward the idea that in such electric fields the discharges (lightning) are produced by extensive air showers arising from high energy cosmic ray particles with e =1014-10 15 eV. These high-energy cosmic rays interact with nuclei of ambient air and give rise to many thousands of charged secondaries. Along ionized tracks of these secondary particles in a strong electric field the avalanches develop and propagate. The high energy cosmic ray particle flux is enough to explain the number of lightning observed. As cosmic rays hit the Earth’s atmosphere accidentally in all directions the lightning arise by chance also. There is another mechanism of lightning production suggested by Gurevich [30, 31] in which relativistic electron is accelerated in the electric field of thundercloud and produces avalanche. In the process of thundercloud formation one can recognize initial, maturity and decay phases. In Fig. 23 these phases and the processes of thundercloud electricity and lightning production are shown schematically [24]. _ J ______9 + + + + + + + + + + + + + + + + + + _ + + + + + _ _+ _+ _+ _ + _ + + + + + _ _ _ _ _ + + 5 + + + _6 ______5_ _ _ + + + + + + 1 4 1 ______2 + + + + + + + _ _ _ _ _ 2 + + + + + + 3 + + + + + + 3 7 + + + + _ + _ + + _+ _ _+ 8 + a + _ + + 10 b c Fig. 23. The phases of thundercloud formation: a Ð initial phase; b Ð maturity phase; c Ð decay phase. Notations: 1 and 2 Ð warm and cold fronts of air; 3-ascending flux of wet air with ions; 4 and 5 Ð extensive air showers produced by primaries with energies e ≥1014 eV and e ≥1015 eV accordingly; 6 Ð intracloud lightning; 7 Ð cloud-to-ground lightning; 8 Ð ground-to-cloud lightning; 9 Ð the negative screen layer; 10 Ð positive charge on the cloud base; J Ð current of negative ions from the ionosphere to the top of the thundercloud. 4.3 The relationship between cosmic rays and other atmospheric phenomena There are several publications in which the changes of cosmic ray fluxes are considered to be responsible for some processes in the atmosphere. Some of such phenomena are listed by Tinsley [16]. Below we analyze the relationships of cosmic ray fluxes with cloudiness, precipitation and lightning. The influence of charged particle fluxes on cloudiness was found by Veretenenko and Pudovkin [32]. They found that the value of cloudiness reduced when cosmic ray fluxes in the interplanetary space and in the atmosphere decreased during so-called Forbush-effects of cosmic rays. As was shown by Stozhkov et al. [33, 34] during Forbush-effects the value of precipitation decreased also. In contrast, when the ionization level is increased due to the invasion of solar flare protons into the atmosphere the precipitation level also increased. These results were obtained from the analyses of the precipitation data recorded at the numerous meteorological stations located in Brazil and the Former Soviet Union. More than two hundreds of Forbush-effects and several tens of solar flare events were analyzed. In the analyses the superposed epoch method was applied. Figures 24 and 25 demonstrate the changes of precipitation in the cases of decreases and increases of cosmic ray fluxes in the atmosphere [33, 34]. The value of precipitation level decrease obtained by superposed-epoch method for more than 70 events of Forbush effects is D0= -(17.4±2.7) %. The probability of the occasional appearance of effect is less than 10-4 if the values of D have a normal distribution. The results on relative increase of rainfall level of D during solar proton events were obtained by superposed-epoch method for more than 53 events of solar proton enhancements. The amplitude of positive increase is D0=(13.3±5.3) %. The probability of effect appearance by chance is less than 10-2. D, % 10 5 0 -35 -25 -15 -5 5 15 25 35 -5 Days -10 -15 -20 Fig. 24. The changes of the daily precipitation level, D, %, relative to mean value evaluated from precipitation data during 1 month before (-30 to –1 days) and 1 month after (1 to 30 days) Forbush decrease event. The day “0” correspond to the Forbush decrease main phase. The precipitation data used in the analyses were obtained in the Former Soviet Union and Brazil. D, % 15 10 5 0 -30 -20 -10 0 10 20 30 -5 Day -10 -15 Fig. 25. The changes of the daily precipitation level, D, %, relative to mean value evaluated from precipitation data during 1 month before and 1 month after solar proton events recorded by ground-based neutron monitors (“0”-day). The link of cosmic ray intensity and global cloud coverage was found by Svensmark and Friis- Christensen [35]. Their results demonstrate the relationship between charged particle fluxes on the Earth’s surface and cloudiness during long-term cosmic ray modulation in the 11-year solar activity cycle. When cosmic ray flux increases cloudiness increases and one can expect that the number of thundercloud (or thundercloud coverage) will increase also. The ion production rate and ion concentration in air grow; the total electric charge in the thunderclouds increases. Thus, the number of lightning has to grow and the relationship between cosmic ray flux or ion production rate and thundercloud discharge number has to take place. Now there are the long-term experimental data on lightning flashes over the United States [36] and the link of lightning number and cosmic ray fluxes can be checked. In Fig. 26 the relationship of lightning number L with ion production rate q is shown. The correlation coefficient between these values is r(L, q)=+0.85±0.09. The values of q were calculated from the data on cosmic ray flux measured in the atmosphere at the middle latitudes. 12 30 s) 2 11 25 q 10 events/year 20 6 ion pairs/(cm 15 9 6 L, 10 L q, 10 10 8 1988 1992 1996 2000 Year Fig. 26. The yearly number of lightning L detected in United States in 1989-1998 (black points, from [36]) and ion production rate q in the air column (h=2-10 km) of the middle latitudes (open points). 4.4 Artificial influence on precipitation The results obtained by Veretenenko and Pudovkin [32], Stozhkov et al. [33], Svensmark and Friis- Christensen [35], Ermakov and Stozhkov [24] show clearly the important role of charged particle fluxes on the cloud, thundercloud formation, and precipitation processes. In the lower atmosphere the changes of cosmic ray flux during Forbush effects (decrease of cosmic ray flux) or solar proton events (increase of flux) is about (2-15) %. In the first case the decrease of precipitation level is observed, in the second one the growth of precipitation takes place (see Figs. 24 and 25). It is possible to increase the flux of charged particles in the lower atmosphere using an electron accelerator onboard airplane. The accelerated electrons can irradiate the cloud increasing ionization level inside the cloud. In turn, it could increase precipitation. The modern linear machines accelerating electrons up to the energy of several tens of MeV have a suitable weight and the energy consumption to be installed onboard airplane, e.g. a machine accelerating electrons up to energy 10 MeV with the current 10 mA has the weight about 1 ton, sizes of 3¥0.5¥1 m3 and the energy consumption of ~1 kW. Let us consider a cloud of 3¥3¥2 km3 in sizes, the top of which is at the altitude ~3 km. The flux of cosmic ray secondaries (mainly relativistic electrons and muons) falling on the upper surface of such cloud is ~7¥109 particles/s. The total energy released by these particles inside the cloud to ionize air atoms equals ~4.5¥106 erg/s. In contrast, the total energy released by the accelerated electrons (~6¥1013 electrons/s) is ~109 erg/s. Thus, the accelerator with the parameters given above could increase the ionization inside the cloud in 10-100 times in comparison with the natural background produced by cosmic rays. In comparison with the methods of artificial influence on the clouds used in practice [37, 38] the irradiation of clouds by accelerated particles is safe because the accelerated electrons and their secondaries will be absorbed in ambient air because of its energy losses. The proposed method could be useful to struggle with droughts and downpours causing floods. 5. CONCLUSION We present the main characteristics of cosmic ray population in the atmosphere and its variability. The experimental data were obtained from the long-term cosmic ray monitoring in the atmosphere from 1957 to now. The main features of cosmic ray variability are the following: ¥ 11-and 22-year solar cycle variations; ¥ solar proton events originating from powerful solar flares; ¥ energetic electron precipitation during geomagnetic disturbances and Forbush decreases of cosmic rays; The cosmic ray monitoring in the atmosphere allows to detect radioactive clouds producing by nuclear explosions or failures at atomic plants. Cosmic ray particle fluxes play important role in many atmospheric processes and only now this role begins to be elucidated. The thundercloud electricity and lightning production, cloud formation, influence on the value of global cloudiness and precipitation on the short (days) and long (11-year solar activity cycle) time scales, operation of global electric circuit and long-scale global climate changes depend on the values of cosmic ray flux. ACKNOWLEDGEMENTS We are very grateful to our colleagues who made hard work to get the experimental data on the long-term cosmic ray variations in the atmosphere. This work is partly supported by Russian Foundation for Basic Research (grants No. 01-02-31005 and No. 99-02-18222). REFERENCES [1] A.N. Charakhchyan, G.A. Bazilevskaya, Y.I. Stozhkov, and T.N. Charakchyan, Cosmic rays in the atmosphere and in space near the Earth during 19 and 20 solar activity cycles, Tr. Fis. Inst. im P.N. Lebedeva, Russian Akad. Nauk, Nauka, 83 (1976) 3 (in Russian). [2] A.E. Golenkov , A.K. Svirzhevskaya, N.S. Svirzhevsky, and Y.I. Stozhkov, Cosmic ray latitude survey in the stratosphere during the 1987 solar minimum, Conf. Pap., Int. Cosmic Ray Conf., XXIst, 7 (1990) 14. [3] H. Moraal, Observations of the eleven-year cosmic-ray modulation cycle, Space Sci. Rev., 19 (1999) 845. [4] Y.I. Stozhkov and T.N. Charakhchyan, On the role of the heliolatitudes of the sunspots in the 11- year galactic cosmic ray modulation, Acta Physics Academial Sciantiarum Hungaricae, (Suppl. 2), 29 (1970). [5] G.A. Bazilevskaya, M.B. Krainev, Yu.I. Stozhkov , A.K. Svirzhevskay, and N.S. Svirzhevsky, Long-term Soviet program for the measurement of ionizing radiation in the atmosphere, Journ. Geomag. and Geoelectr., 43 (1991) 893. [6] G.A. Bazilevskaya., A.N. Charakhchyan, Y.I. Stozhkov, and T.N. Charakhchyan, The energy spectrum and the conditions of propagation in the interplanetary space for solar protons during the cosmic ray events on August 4 to 9, 1972, Conf. Pap., Int. Cosmic Ray Conf., XIIIrd, Denver, USA, 2 (1973). [7] V.S Makhmutov, G.A. Bazilevskaya, A.I. Podgorny, Y.I. Stozhkov, and N.S. Svirzhevsky, The precipitation of electrons into the Earth's atmosphere during 1994, Proc. 24 ICRC, Italy, Rome, 4 (1995) 1114. [8] G.A. Bazilevskaya and V.S. Makhmutov, The electron precipitation into the atmosphere according to cosmic ray experiment in the stratosphere, Izv. AN SSSR, Ser. Fiz., 63 (1999) 1670 (in Russian). [9] V.S. Makhmutov, G.A. Bazilevskaya, M.B. Krainev, Characteristics of energetic electron recipitation into the Earth's polar atmosphere and geomagnetic conditions, Adv. Space. Res., (2001) (in press). [10] G.D. Reeves, Relativistic electrons and magnetic storms: 1992-1995, Geoph. Res. Lett., 25 (1998) 1817. [11] G.A. Bazilevskaya, A.K. Svirzhevskay, N.S. Svirzhevsky, Y.I. Stozhkov, Radioactive cloud in the atmosphere at Moscow site on 12-14 April 1993, Kratkie soobsheniya po fizike, Moscow, Lebedev Instituite, 7-8 (1994) 36 (in Russian). [12] Y.I. Stozhkov, V.I. Ermakov, and P.E. Pokrevsky, Cosmic rays and atmospheric processes, Izv. Russian Akad. Nauk, ser. fiz., 65 (2001) 406 (in Russian). [13] J. Alan Chalmers, Atmospheric Electricity, Pergamon press (1967). [14] R. Reiter, Phenomena in Atmosphere and Environmental Electricity, Amsterdam, Elsvier (1992). [15] R. Markson, Solar modulation of atmospheric electrification and possible implications for the Sun-Weather relationship, Nature, 273 (1978) 103. [16] Brain A. Tinsley, Correlations of atmospheric dynamics with solar wind-induced changes of air- earth current density into cloud tops, Journ. Geophys. Res., 101 (1996) 29,701. [17] L.B. Loeb, Basic Processes of Gaseous electronics, New-York (1960). [18] V.I. Ermakov, G.A. Bazilevskaya, P.E. Pokrevsky, and Y.I. Stozhkov. Ion balance equation in the atmosphere, Journ. Geoph. Res., 102 (1997) 23,413. [19] R.G. Roble, On solar-terrestrial relationships in the atmospheric electricity, Journ. Geoph. Res., 90 (1985) 6000. [20] C.T. Wilson, The maintenance of the Earth's electric charge, Observatory, 45 (1922). [21] E.R. Williams, Electricity of thunderclouds, Scientific American, 1 (1989) 34. [22] M.B. Baker and J.G. Dash, Mechanism of charge transfer between colliding ice particle in thunderstorms, Journ. Geoph. Res., 99 (1994) 10,621. [23] V. Brooks, C.P.R. Saunders, An experimental investigation of the inductive mechanism of thundercloud electrification, Journ. Geoph. Res., 99 (1994) 10,627. [24] V.I. Ermakov and Y.I. Stozhkov, New mechanism of thundercloud electricity and lightning production, Proc. 11-th Intern. Conf. Atmospher. Elect., Alabama, USA (1999) 242. [25] P.N. Tverscoi, , Course of meteorology, Leningrad, Gidrometeoizdat, (1962) (in Russian). [26] A.I. Rusanov and V.L. Kusmin, On the influence of electric field on the surface tension of the polar liquid, Kolloidnyi Journal, 39 (1977) 388 (in Russian). [27] A.I. Rusanov, To thermodynamics of nucleation on charged centers, Doklady Academii Nauk, USSR, 238 (1978) 831 (in Russian). [28] J.M. Meek and J. Craggs, Electrical Breakdown of Gases, Oxford, the Claredon Press (1953). [29] Ermakov , Molnii-sledy chastiz sverchvysokich energii, Nayka I zhizn, Moskva, Prosveshenie (1993) 92 (in Russian). [30] A.V. Gurevich , G.M. Molikh, and R.A. Roussel-Dupre, Runaway mechanism of air breakdown and preconditioning during a thunderstorm, Phys. Lett., A, 165 (1992) 463. [31] A.V. Gurevich, K.P. Zubin, R.A. Roussel-Dupre, Lightning initiation by simultaneous effect of runaway breakdown and cosmic ray showers, Phys. Lett., A, 254 (1999) 79. [32] S.V. Veretenenko and M.I. Pudovkin, Effects of forbush-decreases in cloudiness variations, Geomagn. and Aeronomy, 34 (1994) 38 (in Russian). [33] Y.I. Stozhkov , J. Zullo, Jr., I.M. Martin, G.Q. Pellegrino, H.S. Pinto, P.C. Bezerra, G.A. Bazilevskaya, V.S. Machmutov, N.S. Svirzevskii, A. Turtelli, Jr., Rainfalls during great Forbush-decreases, Nuovo Cimento, 18C (1995) 335. [34] Y.I. Stozhkov , P.E. Pokrevsky, J. Zullo, Jr., I.M. Martin, V.P. Ohlopkov, G.Q. Pellegrino, H.S. Pinto, P.C. Bezerra, A. Turtelli, Jr., Influence of charged particle fluxes on precipitation, Geomagn. and Aeronomy, 36 (1996) 211 (in Russian). [35] H. Svensmark and E. Friis-Christensen, Variation of cosmic ray flux and global coverage - a missing link in solar-climate relationships, Journ. Atmospheric and Solar-Terrest. Physics, 59 (1997) 1225. [36] R.E. Orville, G.R. Huffines, Lightning ground flash measurements over contiguous United States: a ten-year summary 1989-1998, Proc. 11th Intern. Conf. Atmospher. Electr., Alabama, USA, (1999) 412. [37] L.G. Kachurin, Fisicheskie osnovy vosdeistviya na atmosfernye processy, Gidrometeoisdat, Leningrad (1990) (in Russian). [38] New Scientist, New method of production artificial precipitation, 151 (1996) 10. CLOUD PROPERTY SURVEY FROM SATELLITE OBSERVATIONS USING VERTICAL SOUNDERS (TOVS PATH-B) AND IMAGERS (ISCCP) C. J. Stubenrauch and F. Eddounia Laboratoire de Météorologie Dynamique, Ecole Polytechnique, Palaiseau, France Abstract Since about 1980 the use of satellite radiometers allows a continuous survey of cloud properties over the whole globe. We compare the evolution of cloud amount and effective cloud amount from two different cloud climatologies : ISCCP (International Satellite Cloud Climate Project), using imagers onboard geostationary and polar orbiting weather satellites, and TOVS Path-B, obtained from TIROS-N Operational Vertical Sounder (TOVS) measurements onboard polar orbiters. Over a period of about ten years (1983-1995), the average cloud amount of about 67% and the effective cloud amount of about 53% are quite stable. The most important disturbance within this period was the volcanic eruption of Mount Pinatubo in June 1991. The slight increase of the visible reflectances by the volcanic aerosols has led to a slight overestimation of cloud optical thickness by ISCCP and hence to a slight overestimation of low clouds and underestimation of the amount of high clouds during the following year. Since infrared radiation is less affected by volcanic aerosols, TOVS cloud properties should be more reliable. The stable TOVS Path-B effective high cloud amount over the whole period in the tropics indicates that the volcanic aerosols do not change the properties of high clouds on this scale. 1. CLOUD OBSERVATIONS FROM SATELLITE Only satellite observations are capable to give a continuous survey of the state of the atmosphere over the whole globe. At present, twenty years of these measured radiances are available. In order to convert these radiances into cloud properties, complex inversion algorithms are necessary. These algorithms consist of two parts : i) cloud detection and ii) determination of cloud properties using radiative transfer models. In the following we present two global cloud climatologies. Whereas ISCCP is the cloud climatology with the best diurnal sampling and spatial resolution, TOVS Path-B cirrus (semi- transparent ice cloud) properties, obtained from vertical sounders with a relatively high spectral resolution, are especially reliable, day and night. However, one has to keep in mind that both climatologies give information only on the uppermost cloud. For climate studies, using one of these datasets, it is important to understand how cloud properties are perceived by these different instruments and inversion methods. A detailed comparison has shown that both datasets agree quite well [1]. Discrepancies can be explained by differences in temperature profiles, horizontal (partial cloud cover) and vertical (multi-layer clouds) heterogeneities. For example, in the case of thin cirrus overlying low clouds, one determines with TOVS the cirrus properties, whereas ISCCP determines a mixture of both clouds, from the visible channel. 1.1 ISCCP climatology For its global cloud climatology, the International Satellite Cloud Climatology Project (ISCCP) [2,3] puts emphasis on temporal and spatial resolution, rather than on spectral resolution, by using one visible (VIS, day only) and one infrared (IR) atmospheric window radiance measurement from imagers on the suite of geostationary and polar orbiting weather satellites. Time sampling is three hourly and the initial spatial resolution of about 5 km is sampled to intervals of about 30 km, which means that about one pixel out of 36 is kept for cloud information. The first ISCCP dataset has been thoroughly studied (e. g. [4]-[9]). Some of these studies have led to a recent re-analysis [3], mostly improving the treatment of cirrus and polar clouds. ISCCP has reprocessed eleven years of data (D- series) from 1983 until 1993. The processing has just taken up again; by the time of writing data are available until 1997, and the whole data period should be available within the next months. Clouds are detected through a variable IR-VIS threshold test which compares the measured radiances to ‘clear sky composite’ radiances that have been inferred from a series of statistical tests on the space and time variations of the IR and VIS radiances [10]. These clear sky conditions are associated with low IR and VIS spatial and temporal variability. ISCCP cloud properties are determined for each pixel by comparing the observed radiances with a detailed radiative transfer model. The model includes the effects of the atmosphere, with properties specified from the operational analysis of the TOVS data (only one profile per day is available), and surface determined from the clear radiances. Cloudy pixels are assumed to be covered completely by a single homogeneous cloud layer. Cloud-top temperature, Tcld, is first retrieved assuming that all clouds are black bodies. During daytime, when VIS radiances are available to retrieve cloud optical thicknesses, t, the cloud-top temperature of ‘transmissive’ clouds (t < 11) is corrected to account for the radiation transmitted from below. This means that Tcld is decreased as a function of t. In the first version of the ISCCP data (C-series), all clouds were represented in the radiative model by a cloud composed of 10 mm radius spherical liquid water droplets; in the new ISCCP dataset (D-series), clouds with Tcld ≥ 260 K are treated with the same liquid cloud model, but all clouds with Tcld < 260 K are treated with a model cloud composed of 30 mm ice polycrystals [3]. Tcld is converted into cloud-top pressure, pcld, using the operational TOVS atmospheric profiles. Clouds can be classified according to three pcld intervals (separated at 440 and 680 hPa). During daytime, clouds are classified into nine types, by separating each of the three cloud height categories into thin, medium and thick clouds according to three t intervals (divided at 3.6 and 23). Among many other variables, the D2 dataset gives statistics on cloud amount, CA, as calculated in Equation (1), and cloud amount separately for high, midlevel and lowlevel clouds at a spatial resolution of 2.5∞. CA = Ncld / Ntot. (1) wh ere Ncld is the number of cloudy pixels within the grid, and Ntot is the total number of pixels within the grid. 1.2 TOVS Path-B climatology The Improved Initialization Inversion (3I) algorithms [11] convert infrared and microwave radiation measured from the TIROS-N Operational Vertical Sounder (TOVS) onboard the NOAA polar orbiters into atmospheric temperature and humidity profiles and into cloud and surface properties. Within the framework of the NOAA/NASA Pathfinder Program, eight years of TOVS data (1987- 1995) have already been processed at LMD. This TOVS Path-B dataset provides these atmospheric parameters at a spatial resolution of 1∞ [12]. Results for the whole TOVS observation period from 1979 until now should be available by the end of 2001, since the re-calibration of the HIRS brightness temperatures obtained by comparing airmass-averaged brightness temperatures computed from radiosonde measurements to collocated observed brightness temperatures has just taken up again for the extending period. The 3I algorithms are based on i) the Thermodynamic Initial Guess Retrieval (TIGR) dataset, describing ~2000 different atmospheric conditions extracted from a huge collection of radiosonde measurements and ii) a fast line-by-line radiative transfer model, Automatized Atmospheric Absorption Atlas (4A, [13]), simulating clear sky and cloudy radiances at 30 pressure levels. Cloud det ect ion is performed at HIRS spatial resolution (~17 km at nadir) by eight (seven) threshold tests during daytime (nighttime). An important part of the cloud detection is the use of simultaneous MSU radiance measurements. Since the latter probe through the clouds, they are used to predict clear sky IR brightness temperatures which are compared to those of the HIRS instrument for all individual pixels to decide if they are cloudy. A summary of the 3I cloud detection scheme is given in Table 1 of [14]. To insure more coherence with the MSU spatial resolution (~100 km at nadir), the HIRS radiances are averaged separately over clear pixels and over cloudy pixels within 100 km x 100 km regions. Cloud properties are determined from the averaged cloudy pixel radiances assuming that all cloudy pixels are covered by a single homogeneous cloud layer. The average cloud-top pressure, pcld, and the average effective cloud amount over cloudy pixels, Necld, are obtained from four radiances in the 14 mm CO2-absorption band (with peak responses from 400 to 900 hPa levels in the atmosphere) 2 and one in the 11 mm IR atmospheric win dow by minimizing a weighted c [15]. Empirical weights reflect the effect of the brightness temperature uncertainty within a TIGR airmass class on these radiances at the various cloud levels. The method is based on the coherence of the effective cloud emissivity, Necld, in Equation (2), obtained from the five wavelengths at the pressure level of the real cloud. IImi()ll- clr () i Npeelcld() cld@= Np cld (,) k i for i = 4,8 (2) Ipcld(,) kll i- I clr () i wh ere li is the wavelength of HIRS channel i, pk is the pressure level k out of 30 levels, Im is the measured radiance, Iclr is the retrieved clear sky radiance, and Icld is the calculated radiance emitted by a homogeneous opaque single cloud layer. Tcld is obtained from pcld using the retrieved 3I atmospheric temperature profiles. The cloud amount, CA, is again determined as in Equation (1), but this time in each 1∞ grid. The effective cloud amount over a 1∞ grid, eN, is the product of cloud amount and effective cloud emissivity, Necld: eN = CA x Necld (3) Cloud types are defined by the cloud-top pressure and effective cloud amount. High clouds (pcld < 440 hPa) are divided into three categories: opaque (Necld > 90%), cirrus (90% < Necld < 50%) or thin cirrus (Necld < 50%). Since midlevel (440 h Pa < pcld < 680 hPa) and low-level (pcld > 680 hPa) clouds have a smaller horizontal extension, only two classes in each height category can be distinguished: mostly cloudy or overcast (eN > 50%) and partly cloudy (eN < 50%) fields. Their relatively high spectral resolution make infrared sounders very useful for the determination of cloud properties (frequency, altitude, cloud top temperature and effective emissivity), day and night. Their coarse spatial resolution (20 km) has less effect on clouds with large spatial extents like cirrus clouds. In addition to cloud height and effective emissivity, we start to retrieve mean effective ice crystal sizes for cirrus clouds, taking advantage of the fact that spectral cirrus emissivity differences between 8 and 11 mm depend on this parameter [16]. An eight year survey of these cirrus properties will be available within the framework of the European project CIRAMOSA (CIrrus microphysical properties and their effect on RAdiation: survey and integration into climate MOdels using combined SAtellite observations ; web-site : http ://www.lmd.polytechnique.fr/CIRAMOSA/Welcome.html). 2. CLOUD CLIMATE STUDIES The decade of retrieved cloud data is certainly not yet enough to study climate change, but these datasets give a good starting point to study cloud properties and their variations in correlation with natural events like volcanic eruptions or the El Niño event and La Niña event. Whereas in normal conditions the trade winds blow towards the west across the Tropical Pacific and therefore pile up warm surface water in the west Pacific, during El Niño the trade winds relax in the central and western Pacific leading to a penetration of warm water towards the east. El Niño and La Niña are opposite phases of the El Niño-Southern Oscillation (ENSO) cycle, with El Niña sometimes referred to as the cold phase of ENSO and El Niño as the warm phase of ENSO. Within the period from 1983 to 1995, there were one El Niña event (1985) and a following El Niño event (1986/87) and then a rapid succession of El Niño events in the 90's (1991/92, 1993 and 1994). The eruption of Mount Pinatubo in the Philippines in June 1991 has spread a huge amount of sulfate aerosols into the stratosphere which staid in the atmosphere for more than two years [17]. Recently, a link between the variation of galactic cosmic rays intensity and cloud amount on earth has been found by Svensmark et al. [18], which could not be confirmed by Kristjansson et al. [19]. 2.1 Cloud amount variation In order to give a first impression of cloud amount and its variations with time and seasons during the period from 1983 until 1995, we present in Fig. 1 monthly averages of cloud amount from ISCCP (D2) and TOVS Path-B ('3I') compared to the TOVS Path-B effective cloud amount. The ISCCP D2 cloud amount is obtained during daytime from IR and VIS data, and the cloud amount during nighttime, using only IR data, is adjusted to the daytime results after a comparison of both methods during daytime [20]. TOVS Path-B data for this analysis are used only from the NOAA- 10 and NOAA-12 satellites with local observation times at 7h30 am and pm. The NOAA-11 satellite which was launched in 1988 for a local observation time of 1h30 am and pm, has drifted strongly during its operation with a local observation time of 5h30 in 1991. This drift has strong consequences on the diurnal sampling of the data, especially over land where the diurnal cycle of clouds can be strong [7]. In Fig. 1 we observe average global cloud amounts of 67% (ISCCP) and 77% (TOVS Path-B), whereas the effective cloud amount, taking into account also the cloud opacity, is only 53%. The 10% difference between ISCCP and TOVS Path-B cloud amount can be explained 1) by a higher sensitivity to thin cirrus clouds of TOVS due to its better spectral resolution and 2) by the larger HIRS pixel size for which it can be declared cloudy even if it is only half covered by clouds [14]. Cloud amount alone is not a sufficient variable to look for climate changes. One also should look at the thickness of clouds, which is possible only during day with ISCCP. On the other hand, the TOVS Path-B effective cloud amount which is reliable day and night, combines cloud cover and cloud thickness. Radiative effects of clouds depend on effective cloud amount and cloud height. The inversion method also takes the smaller measured radiances of partly covered pixels into consideration. We observe a seasonal cycle in cloud amount with a maximum in northern hemisphere winter and a minimum in summer. A minimum of effective cloud amount in northern hemisphere winter leads to the assumption of a larger amount of thin cirrus clouds during this season. Within a few percent these global cloud properties seem stable during this whole period. globe 0.85 3I effective 0.8 ISCCP 3I 0.75 0.7 0.65 0.6 cloud amount 0.55 0.5 0.45 84 86 88 90 92 94 year Fig. 1. Monthly mean ISCCP cloud amount, TOVS Path-B (3I) cloud amount and effective cloud amount as a function of time. 2.2 Variation of low, midlevel and high clouds In this section we explore cloud amount and effective cloud amount separately for low, midlevel and high clouds, as defined in sections 1.1 and 1.2. We compare in Figs. 2a to 2c monthly mean cloud amount over the globe for these cloud types from ISCCP using IR and VIS radiances ('ISCCP day') and using IR radiances only ('ISCCP IR', day and night). The latter data have been used by Marsh and Svensmark [21] to reveal a correlation between galactic cosmic ray intensity and low cloud amount. TOVS Path-B effective cloud amount ('3I') for these cloud types, which is defined as their frequency times effective cloud amount during appearance, is also shown in Figs. 2a to 2c. For low and midlevel clouds, we notice a slightly larger cloud amount from the 'ISCCP IR' analysis than from the ISCCP day measurements. However, for high clouds, the 'ISCCP IR' analysis yields a 8% smaller cloud amount than the 'ISCCP day' analysis. As the IR analysis converts the IR brightness temperature into cloud-top temperature under the hypothesis that clouds are black bodies, the cloud-top temperature of semi-transparent clouds is overestimated corresponding to an underestimation of their height. Therefore, the 'ISCCP IR' low and midlevel cloud amounts contain in addition also higher, semi-transparent clouds, whereas the 'ISCCP IR' high cloud amount contains only high opaque clouds. A separation into cloud height does only make sense when using the more reliable 'ISCCP day' analysis. globe 35 30 25 20 3I effective low cloud amount 15 ISCCP day ISCCP IR 10 84 86 88 90 92 94 year globe 35 3I effective ISCCP day 30 ISCCP IR 25 20 mid cloud amount 15 10 84 86 88 90 92 94 year globe 35 3I effective ISCCP day 30 ISCCP IR 25 20 high cloud amount 15 10 84 86 88 90 92 94 year Fig. 2. Monthly mean ISCCP cloud amount from VIS and IR radiances (day) and from IR radiances only (IR), as well as TOVS Path-B (3I) effective cloud amount as a function of time a) for low clouds, b) for midlevel clouds and c) for high clouds. The difference between ‘ISCCP day’ cloud amount and TOVS Path-B effective cloud amount should give an estimation of cloud thickness. This difference is negligible for low clouds in northern hemisphere winter, but 5% for low clouds in northern hemisphere summer. Also midlevel clouds and high clouds have in general an about 6% smaller effective cloud amount. This means that on average low clouds are slightly thicker in northern hemisphere winter, and there is a high occurrence of semi-transparent cirrus clouds over the globe (about 25%, not shown). The slight increase of ‘ISCCP day’ low cloud amount and TOVS Path-B effective low cloud amount and decrease of ‘ISCCP day’ high cloud amount are related to the volcanic eruption of Mount Pinatubo as we will show in the following section. 2.3 Evolution of cloud amount In order to look for systematic changes in cloud amount, we analyze running means over twelve months which take out seasonal variations and subtract them from the average over the whole data period (1983-1993 for ISCCP and 1987-1995 for TOVS Path-B). Figs. 3a to 3d show the evolution of ISCCP cloud amount and TOVS Path-B effective cloud amount ('3I') a) over the globe, b) in the tropics (15∞N to 15∞S), c) over northern hemisphere midlatitudes (30∞N to 60∞N) and d) over southern hemisphere midlatitudes (30∞S to 60∞S). globe tropics 2 3I effective 4 3I effective ISCCP 3 ISCCP 1 2 1 0 0 -1 -1 cloud amount cloud amount -2 -3 -2 -4 84 86 88 90 92 94 84 86 88 90 92 94 year year NH mid latitude SH mid latitude 2 3I effective 2 3I effective ISCCP ISCCP 1 1 0 0 -1 -1 cloud amount cloud amount -2 -2 84 86 88 90 92 94 84 86 88 90 92 94 year year Fig. 3. Difference between running mean over twelve months of ISCCP cloud amount and TOVS Path-B (3I) effective cloud amount and mean over the whole data period (1983-1993 for ISCCP and 1987-1995 for TOVS Path-B) as a function of time. Each point is plotted against the seventh month of each one year period. Results are shown a) over the globe, b) in the tropics, c) over northern hemisphere midlatitudes and d) over southern hemisphere midlatitudes. Over the globe, variations are within 1%. Differences in behavior between ISCCP cloud amount and TOVS Path-B effective cloud amount appear after the volcanic eruption of Mount Pinatubo in 1991. This appears even clearer by looking at the tropics where the TOVS Path-B effective cloud amount increased by 4% and decreases only a year later. The maxima of cloud amount in 1987 is probably related to the El Niño event [3]. This effect seems to be the highest in the northern hemisphere midlatitudes. Normally, the El Niño event shifts the tropical convection from West to East, but perhaps the convection also shifts slightly northwards. This has to be studied more in detail by studying geographical maps. 2.4 Evolution of high and low clouds If one wants to study correlations between cloud amount and the variation of galactic cosmic ray intensity, one should look for them at high clouds and at higher latitudes, since the intensity decreases by entering the atmosphere and since the Earth magnetic field is lower at these latitudes. Figs. 4a to 4d shows the evolution of ISCCP high cloud amount from ISCCP using IR and VIS radiances ('ISCCP day') and using only IR radiances ('ISCCP IR', day and night) and of TOVS Path-B effective high cloud amount ('3I') a) over the globe, b) in the tropics, c) over northern hemisphere midlatitudes and d) over southern hemisphere midlatitudes. globe tropics 3 4 3I effective 3I effective 3 2 ISCCP IR ISCCP IR ISCCP day 2 ISCCP day 1 1 0 0 -1 -1 -2 high cloud amount -2 high cloud amount -3 -4 -3 84 86 88 90 92 94 84 86 88 90 92 94 year year NH mid latitude SH mid latitude 3 3 3I effective 3I effective 2 ISCCP IR 2 ISCCP IR ISCCP day ISCCP day 1 1 0 0 -1 -1 high cloud amount -2 high cloud amount -2 -3 -3 84 86 88 90 92 94 84 86 88 90 92 94 year year Fig. 4. Difference between running mean over twelve months of ISCCP high cloud amount from VIS and IR radiances (day) and from IR radiances only (IR), as well as TOVS Path-B (3I) effective high cloud amount and mean over the whole corresponding data period as a function of time. Each point is plotted against the seventh month of each one year period. Results are shown a) over the globe, b) in the tropics, c) over northern hemisphere midlatitudes and d) over southern hemisphere midlatitudes. From these figures we conclude that there are no significant variations of the effective high cloud amount and of the 'ISCCP IR' high cloud amount, corresponding to high opaque clouds only. However, the 'ISCCP day' high cloud amount decreases by 2.5% globally and 4% in the tropics after the volcanic eruption of Mount Pinatubo in 1991. This effect can be explained by the increase of the VIS reflectance by the volcanic aerosols of this eruption which have had an optical thickness of more than 0.1. In the ISCCP cloud property retrieval this increase in VIS reflectance affects then the optical thickness of the clouds, and an overestimation of optical thickness leads then to an underestimation of cloud height [3]. This can be seen in Figs. 5a and 5b which show the evolution ISCCP low cloud amount from ISCCP using IR and VIS radiances ('ISCCP day') and using only IR radiances and of TOVS Path-B effective low cloud amount ('3I') a) over the globe and b) in the tropics. The low cloud amount increase from 'ISCCP day' by 3% in the tropics (and midlevel cloud amount increase by 1.5%, not shown) shows this effect clearly, whereas the 'ISCCP IR' low cloud amount, containing a mixture of low and semi-transparent higher clouds over the whole period, does not show such an increase. Therefore, the increase of the TOVS Path-B effective low cloud amount which is nearly identical to the increase of the 'ISCCP day' low cloud amount should also be slightly overestimated. Whereas the TOVS Path-B cloud properties should in general not be affected by the volcanic aerosols since they are retrieved from IR radiances, there is one cloud test out of eight during day which makes use of the VIS reflectance [14]. This test should be affected by the volcanic aerosols, and since low clouds show a stronger contrast with the surface in albedo than in temperature, the effective cloud amount of low clouds should be more affected. Also, the threshold is lower over ocean (15%) than over land (20%). Therefore, cloud amount over ocean should be more affected. This effect has been analyzed by separating ocean and land and NOAA- 10/NOAA-12 observations and NOAA-11 observations. globe tropics 2 3I effective 4 3I effective ISCCP IR 3 ISCCP IR ISCCP day ISCCP day 1 2 0 1 0 -1 -1 low cloud amount -2 low cloud amount -2 -3 84 86 88 90 92 94 84 86 88 90 92 94 year year NH mid latitude SH mid latitude 2 3I effective 2 3I effective ISCCP IR ISCCP IR ISCCP day ISCCP day 1 1 0 0 -1 -1 low cloud amount low cloud amount -2 -2 84 86 88 90 92 94 84 86 88 90 92 94 year year Fig. 5. Difference between running mean over twelve months of ISCCP low cloud amount from VIS and IR radiances (day) and from IR radiances only (IR), as well as TOVS Path-B (3I) effective low cloud amount and mean over the whole corresponding data period as a function of time. Each point is plotted against the seventh month of each one year period. Results are shown a) over the globe, b) in the tropics, c) over northern hemisphere midlatitudes and d) over southern hemisphere midlatitudes. Figs. 6 show the evolution of these different TOVS Path-B cloud amounts a) over ocean and b) over land. Indeed, over ocean, where there is no effect with NOAA-11 observations at 1h30 am, taken during nighttime, one observes a strong cloud amount increase for the NOAA- 10/NOAA-12 observations which are taken when the sun has a low zenith angle, already a difficult time for analyzing VIS reflectances. Over land, there is no effect, with exception at 7h30 pm where a much smaller cloud amount overestimation than over ocean can be seen. However, due to the basis of cloud emissivity coherence in our retrieval method the strong overestimation in cloud amount is nearly compensated by a smaller retrieved effective cloud emissivity, so that the final effect for low clouds is only about 3% in the tropics. As we have seen in Figs. 4, effective high cloud amount is not affected by volcanic aerosols. Nevertheless, we will improve this cloud detection test in the next TOVS Path-B re-analysis. tropical ocean tropical land 100 100 NOAA1012 am NOAA1012 am 95 NOAA11 am 95 NOAA11 am 90 NOAA1012 pm 90 NOAA1012 pm NOAA11 pm NOAA11 pm 85 85 80 80 75 75 cloud amount 70 cloud amount 70 65 65 60 60 88 89 90 91 92 93 94 95 88 89 90 91 92 93 94 95 year year Fig. 6. Running mean over TOVS Path-B (3I) cloud amount from different satellite observations as a function of time. Each point is plotted against the seventh month of each one year period. Results are shown for a) tropical ocean and b) tropical land. 3. CONCLUSION AND OUTLOOK Satellite observations provide a unique possibility to survey cloud properties over a long period of time. During the observed decade (1983-1995), the average cloud amount of about 67% and effective cloud amount of about 53% are stable within 2% over the globe. Within this period, the most important disturbance was the volcanic eruption of Mount Pinatubo in the Philippines in June 1991. The volcanic aerosols which had an optical thickness of more than 0.1 slightly increased the VIS reflectances. Whereas the ISCCP cloud amount and TOVS Path-B effective high cloud amount were not affected by this event, the 'ISCCP day' high cloud amount is slightly underestimated (4.5% in the tropics and 2.5% over the globe) and the 'ISCCP day' low amount overestimated (4% in the tropics and 1% over the globe). Nevertheless, a separation into cloud height does only make sense when using the more reliable 'ISCCP day' analysis, since the 'ISCCP IR' low and midlevel cloud amounts contain in addition also higher, semi-transparent clouds, whereas the 'ISCCP IR' high cloud amount contains only high opaque clouds. The TOVS Path-B effective low cloud amount is also slightly overestimated (4% in the tropics and 1% over the globe), because one of the cloud detection tests makes use of the VIS reflectance. As has been pointed out earlier, if there are any correlations between the galactic cosmic ray intensity and cloud properties these should occur for high clouds and at higher latitudes. The TOVS Path-B effective high cloud amount which is stable within 1% over the whole observation period does not show any correlation with the cosmic ray intensity variations. It should be noted that cloud radiative effects are determined not only by cloud amount but depend also on cloud thickness and effective cloud amount includes both variables. Both datasets will be extended. By the time of writing, ISCCP data processing has already started again, and the whole data period (1983 until now) should be available within the next months. TOVS Path-B processing just started with the re-calibration of the HIRS brightness temperatures, and the whole TOVS observation period from 1979 until now should be available by the end of 2001. ACKNOWLEDGEMENTS Our thanks to W. B. Rossow for many stimulating discussions within the last years in order to advance in understanding both cloud datasets. We also want to thank the rest of our ARA (Analyse du Rayonnement Atmosphérique) group for their support, and especially S. Serrar and R. Armante for their help in computational matters and G. Rädel for bringing our attention to CERN again. REFERENCES [1] Stubenrauch, C. J., W. B. Rossow , N. A. Scott, and A. Chédin, Clouds as Seen by Satellite Sounders (3I) and Imagers (ISCCP): III) Spatial Heterogeneity and Radiative Effects, J. Climate 12 (1999) 3419-3442. [2] W. B. Rossow and R. A. Schiffer, ISCCP Cloud Data Products. Bull. Amer. Meteor. Soc. 72 (1991) 1-20. [3] W. B. Rossow and R. A. Schiffer, Advances in understanding clouds from ISCCP. Bull. Amer. Meteor. Soc. 80 (1999) 2261-2287. [4] R. Fu, A. D. Del Genio, and W. B. Rossow, Behavior of deep convective clouds in the tropical Pacific deduced from ISCCP radiances. J. Climate 3 (1990) 1129-1152. [5] G. Tselioudis, W. B. Rossow, and D. Rind, Global patterns of cloud optical thickness variation with temperature, J. Climate 5 (1992) 1484-1495. [6] S. A. Klein and D. L. Hartmann, The seasonal cycle of low stratiform clouds. J. Climate 6 (1993) 1587-1606. [7] B. Cairns, Diurnal variations of cloud from ISCCP data. Atm. Res. 37 (1995) 133-146. [8] X. Liao, W. B. Rossow and D. Rind,: Comparison between SAGE II and ISCCP High-Level Cloulds. Part II: Locating Cloud Tops. J. Geophys. Res. 100 (1995) 1137-1147. [9] Y. Jin, W. B. Rossow and D. P. Wylie, Comparison of the Climatologies of High-Level Clouds from HIRS and ISCCP. J. Climate. 9 (1996) 2850-2879. [10] W. B. Rossow and L. C. Garder, Cloud detection using satellite measurements of infrared and visible radiances for ISCCP. J. Climate 6 (1993) 2341-2369. [11] A. Chédin, N. A. Scott, C. Wahiche and P. Moulinier, The Improved Initialized Inversion method: A high resolution physical method for temperature retrievals from the TIROS-N Series. J. Clim. Appl. Meteor. 24 (1985) 124-143. [12] N. A. Scott, A. Chédin, R. Armante, J. Francis, C. J. Stubenrauch, J.-P. Chaboureau, F. Chevallier, C. Claud and F. Chéruy, Characteristics of the TOVS Pathfinder Path-B Dataset. Bull. Amer. Meteor. Soc. 80 (1999) 2679-2701. [13] N. A. Scott, and A. Chédin, A fast line-by-line method for atmospheric absorption computations: The Automized Atmospheric Absorption Atlas. J. Appl. Meteor. 20 (1981) 802-812. [14] C. J. Stubenrauch, W. B. Rossow , F. Chéruy, N. A. Scott and A. Chédin, Clouds as Seen by Satellite Sounders (3I) and Imagers (ISCCP): I) Evaluation of Cloud Parameters, J. Climate 12 (1999) 2189-2213. [15] C. J. Stubenrauch, A. Chédin, R. Armante and N. A. Scott, Clouds as Seen by Satellite Sounders (3I) and Imagers (ISCCP): II) A New Approach for Cloud Parameter Determination in the 3I Algorithms, J. Climate 12 (1999) 2214-2223. [16] C. J. Stubenrauch, R. Holz, A. Chédin, D. Mitchell and A. J. Baran, Retrieval of Cirrus Ice Crystal Sizes from 8.3 and 11.1 mm Emissivities Determined by the Improved Initialization Inversion of TIROS-N Operational Vertical Sounder Observations, J. Geophys. Res. 104 (1999) 31793-31808. [17] M. P. McCormick, L. W. Thomason, and C. R. Trepte, Atmospheric effects of the Mt Pinatubo eruption, Nature 373 (1995) 399-404. [18] H. Svensmark and E. Friis-Christensen, Variation of cosmic ray flux and global cloud coverage: A missing link in solar climate relationships, J. Atmos. Sol. Terr. Phys. 59 (1997) 1225-1232. [19] J. E. Kristjansson, and J. Kristiansen, Is there a cosmic ray signal in recent variations in global cloudiness and cloud radiative forcing?, J. Geophys. Res. 105 (2000) 11,851- 11,863. [20] W. B. Rossow, A. W. Walker, D. Beuschel, and M. Roiter, International Satellite Cloud Climatology Project (ISCCP) description of new cloud datasets, World Climate Research Programme (ICSU and WMO) WMO/TD-No.737 (1996) 115pp. [21] N. Marsh, and H. Svensmark, Cosmic rays, clouds, and climate, Space Science Reviews 94 (2000) 215-230. ATMOSPHERIC ELECTRICITY AND CLOUD MICROPHYSICS R. G. Harrison Department of Meteorology, The University of Reading, P.O. Box 243, Reading RG6 6BB, UK Abstract The terrestrial atmospheric electrical system covers a range of dimensional scales, from charged molecular clusters to convective cloud systems. Charge-exchange associated with thunderclouds leads to positive charge in the upper conductive regions of the atmosphere and a net negative charge on the planetary surface. In non-thunderstorm regions, a vertical ionic current flows, replenishing the air with molecular ions otherwise removed by attachment, recombination or nucleation processes. Ions may have indirect effects on non-thunderstorm clouds, and therefore conceivably on climate, via cloud microphysical processes. Cloud Condensation Nuclei (CCN) and Ice Nuclei (IN) are necessary for the formation of water clouds and freezing of ice clouds respectively. In both cases, ionisation may be important: it is now known that ultrafine aerosol can be formed from ionisation, probably providing an additional source of CCN. It is also known that electrified aerosol, perhaps active as IN, can be collected by droplets more effectively than neutral particles. 1. INTRODUCTION Electrical processes in atmospheric air arise from the combined effect of natural ionisation and the natural electric fields generated indirectly by charge separation in thunderclouds. In non- thunderstorm regions, which probably constitute the majority of the global cloud area, the electrical processes will not generate the large breakdown electric fields associated with lightning, but microscopic aerosol particles acquire charges by diffusion of the molecular cluster ions formed from ionisation. In this overview, the effect of small charges on aerosol particles and droplets are considered. Since the charge arises from radiolysis of air by cosmic rays and natural radioactivity, the discussion here is structured in terms of the processes associated with charge generation and removal, including tutorial material on microphysical cloud processes. It has been observed (Marsh and Svensmark, 2000) that there is a correlation between low cloud properties and the neutrons produced by cosmic rays. 2. THUNDERSTORMS AND GLOBAL ELECTRIFICATION The atmospheric electrical system originally discussed by Wilson (1929) can be simplified into an electric circuit in which thunderstorms separate charge in convective regions. The charge separation leads to a potential difference between conductive regions of the upper atmosphere and the surface, which causes an ionic leakage current to flow vertically (figure 1). Currents of order 2000A flow in the circuit, with an upper atmosphere potential of ~300kV. The conduction current density in undisturbed regions is ~2pA.m-2. The charge-exchange processes within thunderclouds are complicated, and probably result from the interaction between rising ice crystals rising and riming soft hail (graupel). Typical microphysical collisions exchange charge with typical magnitudes of tens of femtoCoulombs, but the precise magnitude and polarity is greatly influenced by the liquid water content and temperature (MacGorman and Rust, 1998). Figure 1. The global atmospheric electrical circuit (from Harrison, 1997). 3. ATMOSPHERIC PROPERTIES AND CLOUD MICROPHYSICS 3.1 Atmospheric properties 3.1.1 Bulk properties The troposphere (lower atmosphere) shows variations in temperature and water content, and partitioning of the water concentration between liquid, solid or vapour forms is critical to the formation and distribution of clouds. Figure 2 shows a vertical sounding of temperature and humidity, which illustrates the atmospheric structure. The presence of low cloud (which was observed from the surface) is evident from the sharp increase in relative humidity, marked as A. B shows a slight temperature inversion associated with the top of the planetary boundary layer, and at C the temperature ceases to fall with height, at the tropopause. It is clear that there is considerable variability in the relative humidity during the ascent, and in the region where cloud was identified optically. Radiosonde ascent on 22nd March 1998 at 1400h, 1000 (launched from Medina Valley, Isle of Wight, UK) 80 pressure 60 800 air temperature 40 relative humidity Relative Humidity(%) 20 600 temperature (degC) and 0 C 400 A -20 -40 200 B -60 pressure /mbar altitude/m 0 -80 0 5000 10000 15000 20000 Figure 2. Vertical atmospheric sounding in non-frontal synoptic conditions, showing relative humidity, temperature and pressure variations with height. (A, B and C are discussed in the text.) 3.1.2 Microphysical properties In addition to variability in temperature and humidity, there is a considerable variety in the sizes and abundance of aerosol particles and cloud droplets present in the atmosphere. Figure 3 shows a comparison of the sizes of cloud droplets, raindrops, and a condensation nucleus. The typical molecular cluster comprising an atmospheric small ion will have a diameter less than one nanometre. Figure 3. Size spectrum of particles present in a typical cloud. (from Rogers and Yau, 1989). 3.2 Cloud Microphysics The concentration of water vapour in air can be determined by its gaseous partial pressure, and, at any given temperature there is an associated maximum value of partial pressure due to water vapour, the saturation vapour pressure. Air containing sufficient water vapour to generate the saturation vapour pressure is saturated, with a relative humidity of 100%. Slightly greater relative humidities (supersaturations) can occur in localised regions, but they are never greater than a few percent, because of the abundance of aerosol particles on which the water can condense. Many different kinds of aerosol particles are capable of acting as condensation nuclei. Below 0¡C however, liquid water droplets may persist without freezing, although 0¡C is the temperature at which ice melts. Any liquid water droplet with a temperature below 0¡C is supercooled, in a thermodynamically unstable state in which freezing may be readily initiated by heterogeneous or homogeneous nucleation. In heterogeneous nucleation, the supercooled water freezes as a result of the presence of a suitable ice nucleus. Homogeneous nucleation occurs if cooling is continued further, and all supercooled water in atmospheric clouds becomes ice at temperatures colder than - 40¡C by this process. 3.2.1 Saturation vapour pressure, temperature and relative humidity At any given temperature T, the maximum partial pressure of water vapour, the saturation vapour pressure es(T) is given by the Clausius-Clapeyron equation as 1 des l = 2 (1) es dT RTv where l is the latent heat of vaporisation of water, and Rv is the gas constant for water vapour -1 -1 (461.5 J.kg .K ). es(T) can in principle be found by integration from equation (1), but l is also a function of temperature, which leads to many empirical formulae for es(T). Common forms include the exponential (Magnus) equation e.g. es(T)=6.112 exp [17.67 T / (243.5 + T) ] (2) where es is given in millibars and T is in Celsius. The relative humidity is the actual vapour pressure expressed as a fraction of es, at the same temperature. Supersaturation is expressed either as a percentage relative humidity greater than 100%, or as a saturation ratio S. (101% RH = 1% supersaturation = saturation ratio S =1.01). 3.2.2 Activation of condensation nuclei In the troposphere, supersaturations are never greater than a few percent, and are typically rather less. Consequently direct condensation onto ions, which permits visualisation of particle tracks in a Cloud Chamber (S ~ 4), cannot occur in the lower atmosphere. Condensation on aerosol particles, which are larger, does occur, however, and the minimum size of particle necessary depends on the degree of supersaturation. All aerosol particles are therefore potentially able to act as condensation nuclei (CN), if the supersaturation is sufficiently large, but it is the subset of particles able to cause condensation at atmospheric supersaturations which is of interest in cloud physics. These condensation nuclei are known as Cloud Condensation Nuclei (CCN). The vapour pressure over the curved water surface of a particle of radius r, es(r) is greater than es over a plane surface at the same temperature. If condensation occurs on a particle, its growth rate is proportional to the difference in between the bulk vapour pressure e and es(r). For e - es(r) > 0 the cloud droplet grows. This situation is rather more complicated in a mixed-phase cloud due to the differences in vapour pressure over ice and supercooled water. Ice particles grow at the expense of supercooled water in a mixed-phase system. shrinks grows saturation ratio S È b ˘ Êaˆ S(r) = 1- exp Í 3 ˙ Á ˜ Î r ˚ Ë r¯ 0.6% activated at 0.13µm ultrafines and ions Figure 4. Activation of particles at typical atmospheric saturation ratios. The maximum in the saturation ratio curve S(r) defines the minimum radius of particle required to act as a nucleus on which a cloud droplet to grow. For a supersaturation of 0.6%, a 0.13µm radius particle is required. A droplet smaller than this will evaporate. The function S(r) principally depends on a “curvature” term a, and a “solution” (dissolved salt) term b. (after Rogers and Yau, 1989). 3.2.3 Supercooling and ice nucleation Supercooled droplets are common under atmospheric conditions, and result from water droplets cooling in the absence of suitable ice nuclei (IN) to permit heterogeneous ice nucleation. At temperatures cooler than -40¼C all supercooled droplets begin to freeze by homogeneous nucleation. Only very few atmospheric aerosol particles can act as IN, typically less than 1% although the exact fraction increases as the droplets become colder. The ability of a particle to act as an ice nucleus depends on a variety of physical properties including its shape, solubility, crystal structure and its history in cloud processing. At warmer temperatures (-6 to -10¼C) ice multiplication occurs by mechanical production of ice splinters on freezing, generating additional ice fragments which are also able to act as IN. Figure 5 shows the dramatic temperature change occurring when a supercooled water droplet freezes, releasing latent heat. Temperature within 1µL water drop during supercooling and freezing 5.00 0.00 Td (drop) Te (environment) -5.00 -10.00 mperature / degC te -15.00 -20.00 700 750 800 850 time (from beginning of cooling) /seconds Figure 5. Time series of temperature Td within a supercooled water droplet, in an environment at temperature Te, as it freezes and releases latent heat (from Harrison and Lodge, 1998). 4. ATMOSPHERIC IONISATION AND ELECTRIFICATION 4.1 Steady-state ion concentrations Ion-pairs are continually produced in the atmosphere by radiolysis of air molecules, figure 6. The ions produced are rarely single species but clusters of water molecules around a central ion. Typical atmospheric ion concentrations in unpolluted air and fine weather are about 500 ions.cm-3 (Chalmers, 1967). There are three principal sources of high-energy particles which cause radiolysis: Radon isotopes, cosmic rays and terrestrial gamma radiation. The partitioning between the sources varies vertically. Near the surface, ionisation from turbulent transport of radon and other radioactive isotopes is important, together with gamma radiation from isotopes below the surface. Ionisation from cosmic rays is always present, comprising about 20% of the ionisation at the surface. The cosmic fraction increases with increasing height in the atmosphere and dominates above the planetary boundary layer. Figure 6. Formation of small ions by radiolysis of air molecules. Small ions consist of clusters of water molecules collected around a singly charged ion. They have a lifetime of the order of a hundred seconds. Clusters such as H3O+(H2O)n, H+(H2O)n, NO+(H2O)n and NO2+(H2O)n are common for the positive ions and O2-(H2O)n, CO4-(H2O)n, NO-(H2O)n or NO2-(H2O)n for the negative ions (Volland, 1984). The chemical difference between the species in the positive and negative ions leads to some physical asymmetries in the ion properties, with the negative ions more mobile. The ratio of mobilities m-/m+ ~ 1.2. 4.2 Ion balance equation Atmospheric small ions of both signs with number concentrations n+ and n- are governed by • • dn± =-qnnnabm - () aNada () (3) dt ±±Ú Â ±1,jj a =0 j =-• where the ions are produced at a rate q per unit volume. Ions (which are assumed to carry unit charges) are removed by ion-ion recombination (with recombination coefficient a), and by attachment to aerosol particles, which causes charge transfer to the aerosol. The aerosol attachment rate b±1,j(a) depends on aerosol particle radius a and the number of elementary charges j present on the aerosol particle of radius a (Gunn, 1954). In equation (3), the size and charge distributions of atmospheric aerosol particles are accounted for by the integral of number concentration N(r) over all particle radii, and by a sum across all possible particle charges at each radius. Recombination is the principal loss mechanism of ions in clean, aerosol free air. If aerosol is present, then ions are also lost by aerosol attachment. It is instructive to simplify the ion balance equation by neglecting the ion sign (i.e. n+ ª n- = n) and replacing the aerosol particle size distribution by an equivalent monodisperse particle number concentration Z. The ion-aerosol equation can then be written as dn =-qnnZab2 - (4) dt 4.2.1 Time dependent solution Integrating this equation gives the ion concentration n as a function of time t, for a zero initial ion concentration at time zero, as ÈÊ -+()ba22Zqt4 ˆ ˘ 22 - []-+()babZqZ4 - ÍË1 e ¯ ˙ nt()= Í ˙ (5) a Ê -+()ba22Zqt4 ˆ 2 Í 1 + e ˙ ÎË ¯ ˚ which highlights two interesting points. Firstly, if the ion-pair production rate q is uniform and the removal rates are also steady, the ion concentration tends to a steady value for large values of t. Secondly the equation can be simplified according to the situations in which attachment or recombination dominates as the removal mechanisms, according to whether an2 or nbZ is the bigger term. In the atmosphere in polluted air, these terms are roughly comparable, and therefore all the terms in equation (5) have to be evaluated. 4.2.2 Recombination Limit In the case of ion loss solely by recombination, such as in relatively aerosol-free regions of the atmosphere, equation (5) reduces to È 1- e-2 aqt ˘ È q ˘ Í()˙ nt()= Í ˙ (6) Îa ˚ Í 1 + e-2 aqt ˙ Î()˚ 1/2 and the steady-state concentration after a long time has elapsed is given by n• = (q/a) . Inserting typical atmospheric values of q ª10 ion-pairs cm-3 s-1and a =1.6 x 10-6 cm3 s-1 gives -3 n• = 2500 ion-pairs cm . Typical values of small ion concentrations observed in mountain air are about 500 ions cm-3 of each sign, suggesting that attachment processes are almost always significant in modulating the ion concentrations in the lower troposphere. 4.3 Aerosol electrification Collisions between the ions and atmospheric aerosol lead to charge-exchange and electrification of the aerosol, and the ion asymmetry ensures that the collisions do not lead to an average charge of zero. Local electric fields can cause further asymmetries, by depletion of one sign of ion concentration, and consequently substantial aerosol electrification can occur in such regions. The number concentration Nj of monodisperse aerosol particles carrying j elementary charges is given by the Modified Boltzmann Distribution (Clement and Harrison, 1992), as N j È 2 ˘ È 22 ˘ j Èn++m ˘ 8pe0 akT je -je = Í ˙ 2 sinhÍ ˙ expÍ ˙ (7) N0 În--m ˚ je Î88pe0 akT ˚ Î pe0 akT ˚ where m± are the positive and negative small ion mobilities, n± their number concentrations, T the temperature, e the modulus of the electronic charge, k Boltzmann’s constant and e0 the permittivity of free space. The mean charge J (Gunn, 1955) is given by 4peakT Èn++ m ˘ J = 0 lnÍ ˙ (8) 2 m e În--˚ Charge distribution on water droplets radius (3.32 ± 0.65)µm, ion ratio 0.82 0.0800 0.0600 N ( j 0.0400 ) / Z 0.0200 0.0000 -40 -35 -30 -25 -20 -15 -10 -5 0 5 10 15 number of particle charges,j experiment MBD theory Figure 7. Charge distribution on water droplets in the presence of ion asymmetry using the Modified Boltzmann Distribution (MBD). (Experimental data from Gunn and Woessner, 1956) 5. OBSERVED TROPOSPHERIC ELECTRICAL PROPERTIES UNDER NON- THUNDERSTORM CONDITIONS Many electrical soundings of the atmosphere have been made during disturbed (thunderstorm) conditions, but few such measurements have been made under more quiescent atmospheric conditions. There are difficulties with in-cloud measurements, as a balloon or aircraft platform will be required: this may itself introduce difficulties with sampling, particularly under small electric fields. Regions of aerosol particles, some of which may acquire appreciable charges, are the principal perturbation to atmospheric electric fields under fair weather conditions. In general the upper and lower surfaces of a horizontal region of particles or droplets will charge from ions flowing vertically as a result of the fair weather conduction current, and the region within the layer will be of low conductivity compared with that of the surrounding air. Figure 8 shows electric field profiles observed in non-thunderstorm clouds, which are typically two orders of magnitude smaller than equivalent profiles determined in thunderstorms. Figure 8. Typical electric field profiles found in non-thunderstorm clouds under (a) liquid water and (b) supercooled conditions (from MacGorman and Rust, 1998). Optical link V to f PLL electrometer SONDE electrodes radio link P,T,U RECEIVER PLL f to V Logger Figure 9. Summary schematic of an atmospheric electric field sensor using vertically-spaced spherical electrodes, flown under a conventional radiosonde balloon. The electrometer circuit and batteries are mounted within the lower electrode, with the data conveyed via a voltage-to-frequency converter and an optical link operating at 100kHz. A phase-locked loop (PLL), recovers the data from an optical receiver, and the signal is injected into a standard meteorological RS80 radiosonde measuring pressure P, temperature T and relative humidity U. The uhf receiver recovers the 100kHz signal, which a further PLL converts back to a voltage. The voltage is logged at 20Hz by a computer and analogue to digital converter (from Harrison, 2001). A modern sensor suitable for use in low fields using a standard meteorological radiosonde has recently been described (Harrison, 2001), using the displacement current to detect charged aerosol particles from the changes caused in electric field. Figure 9 shows a schematic of the system, which is disposable. Regions of space charge of ~10nC.m-3 were reported in shallow layers, with electrical structures suggesting electrification by the conduction current. 6. DIRECT INFLUENCES OF ELECTRIFICATION AND RADIOACTIVITY ON CLOUDS In considering how natural atmospheric ionisation might influence cloud physics or cloud formation, Harrison (2000) identified two possible routes: (1) direct processes, such as production of new aerosol (e.g. sulphate) by gas-to-particle conversion (GPC) or homogeneous nucleation (2) indirect processes, such as the modification of existing heterogeneous nucleation processes by affecting the condensation nuclei (CN) or ice nuclei (IN). 6.1 Direct processes 6.1.1 CN production In the presence of high levels of radioactivity, the radiolytic formation of particles has been shown to occur. Bricard et al (1968) found the CN concentration could be made to cycle by the regular addition of Thoron (a short-lived source of a-particles), figure 10. Vohra et al (1984) observed particle formation in artificial air in the presence of trace concentrations of sulphur dioxide, ozone and ethene, with naturally-occurring radon concentration levels, suggesting that ultrafine aerosol production could occur in atmospheric air under natural conditions. Recent theoretical work by Yu and Turco (2001), further strengthens the expectation that radiolytic particle production will be found in the atmosphere, and that the conventional ion-aerosol balance equations are incomplete. Figure 10. Particles formed in filtered Parisian air, with the addition of regular cycles of Thoron (from Bricard et al, 1968). 6.2 Indirect processes 6.2.1 Scavenging of charged particles The removal of aerosol particles by cloud droplets, scavenging, is known to be influenced by many factors, including electrical forces (Pruppacher and Klett, 1997). Figure 11 shows the partitioning of the water drop charge Qd in response to the image charge induced by the charged aerosol particle (carrying a charge Qa) brought close to the drop. Charge conservation requires that Qd = I + D. In magnitude, I = -(A/s)Qa, where s is the separation distance between the aerosol and drop centres. The image charge is located at a distance c from the centre of the drop (Jackson, 1975). Summing the Coulomb and image forces, the net electrical force acting between the particles' centres is 1 È QI QD˘ F = aa+ (9) e pe Í 22˙ 4 0 Î()bc- b ˚ where a positive Fe is repulsive. CONDUCTING waterdrop charge Qd=Je drop radius A WATERDROP drop charge D c image charge I Charged aerosol trajectory s aerosol charge Qa = j e RADIOACTIVE radioactive decay rate h CHARGED ion pair yield per decay I AEROSOL residual aerosol charge m aerosol radius a Figure 11. Schematic of the (radioactively) charged aerosol and the image charge I construction within a water drop of radius A. The aerosol and falling drop carry charges Qa and Qd respectively, with Q d = I + D , where D is the non-image charge considered at the centre of the drop (from Tripathi and Harrison, 2001). 6.2.2 Electrofreezing Tinsley et al (2000) have shown that the electrical image force is very significant in aerosol- droplet collisions as, unlike the Coulomb force, it is always attractive between the charged aerosol and water droplet at small separations. This process, electroscavenging, is a subset of many processes described more generally as electrofreezing, in which electrical fields or charges influence the freezing of supercooled droplets. Tinsley and Dean (1991) argued that modification of the electrical properties of aerosols might change their efficacy of aerosol as contact ice nuclei, ultimately leading to storm intensification by triggering latent heat release. Direct ionisation has, however, recently been shown not to lead to freezing of supercooled water (Seeley et al, 2001). There is currently no definitive evidence that charging influences ice nuclei efficiency or that contact nucleation of ice is the dominant freezing mechanism. 7. DISCUSSION Ionisation in the atmosphere is ubiquitous and part of the atmospheric electrical system, which transports ions in low electrical field regions of the atmosphere, such as clear air and non- electrified clouds. There appear at least two ionisation-related processes of relevance to cloud formation: (1) aerosol electrification (2) ultrafine aerosol production 7.1 Aerosol electrification Stratified regions of aerosol will charge in the atmosphere under quiescent conditions as a result of ion transport by the conduction current. Although the charges carried are unlikely to be sufficiently large to initiate bulk discharge processes such as lightning, the distribution of charges expected on aerosols under natural ion asymmetry may yield a small fraction of particles with significant charge levels. Such highly-charged aerosol would normally be rapidly neutralised by atmospheric ionisation, but in ion-depleted regions arising from large aerosol concentrations, moderately-high charge levels might persist. Scavenging processes are influenced by aerosol charge; the water drop charge has negligible effect by comparison. Since heterogeneous ice nucleation requires the collection of suitable aerosol able to operate as ice nuclei, it is therefore conceivable that aerosol charging could influence ice formation. If charge were itself shown to be a enhancing effect for ice nuclei, then the synergy between the effects of increased collection and efficient nuclei could be potent. 7.2 Ultrafine particle production Radiolytic particle production has been observed in laboratory air at atmospheric levels of ionisation, but the particles formed are small and not, at their formation, able to act as cloud condensation nuclei. The recent theoretical work of Yu and Turco (2001) is, however, compelling, in that it shows that in aerosol-deficient regions, such as marine stratus cloud, cosmic ray ionisation could provide an appreciable source of particles. Further microphysical modelling is required to show that sufficient ultrafine particles can survive to become cloud condensation nuclei before the effect on suitable clouds can be assessed. 8. CONCLUSIONS The physical processes, if any, leading to the cosmic ray-low cloud correlation observed by Marsh and Svensmark (2000) remain to be established in the atmosphere. As discussed above, there are atmospheric electrical mechanisms relating ionisation to cloud which remain relatively unexplored in atmospheric physics, and in suitable cloud, could conceivably offer physical explanations for the observed correlation. However without numerical and theoretical estimates of their significance, it is currently impossible to regard ionisation effects as irrelevant to cloud processes. REFERENCES Bricard F., Billard F. and Madelaine G., (1968), Formation and evolution of nuclei of condensation that appear in air initially free of aerosols, J. Geophys. Res., 54, 39-52 Chalmers J.A. (1967), Atmospheric Electricity, 2nd edition, Pergamon Press, Oxford Clement C.F. and Harrison R.G. (1992), The charging of radioactive aerosols J. Aerosol Sci. 23, 5, 481-504 Gunn R. (1954), Diffusion charging of atmospheric droplets by ions, and the resulting combination coefficients, J. Meteorol., 11, 339 Gunn R. (1955), The statistical electrification of aerosols by ionic diffusion, J. Coll. Sci., 10, 107- 119 Gunn R. and Woessner R.H. (1956), Measurements of the systematic electrification of aerosols, J. Coll. Sci., 11, 254-259 Harrison R.G. (1997), Climate change and the global atmospheric electrical system Atmos. Environ. 31, 20, 3483-3484 Harrison R.G. (2000), Cloud formation and the possible significance of charge for atmospheric condensation and ice nuclei Space Science Reviews, 94, 381-396 Harrison R.G. (2001), A balloon-carried electrometer for high-resolution atmospheric electric field measurements in clouds Rev Sci Inst, 72, 6, 2738-2741 Harrison R.G. and Lodge B.N. (1998), A calorimeter to detect freezing in supercooled water droplets Rev Sci Inst, 69, 11, 4004-4005 Jackson J.D. (1975), Classical Electrodynamics Wiley Marsh N.D. and Svensmark H. (2000), Low cloud properties influenced by cosmic rays, Phys. Rev. Lett., 85, 23, 5004-5007 MacGorman D.R. and Rust W.D. (1998), The electrical nature of storms, OUP Pruppacher, H.R. and Klett, J.D., (1997). Microphysics of clouds and precipitation, 2nd edition, Kluwer Rogers and Yau (1989), A short course in Cloud Physics, Pergamon Press Seeley L.H., Seidler G.T. and Dash J.G., (2001), Laboratory investigation of possible ice nucleation by ionizing radiation in pure water at tropospheric temperatures, J. Geophys. Res., 106 ,D3, 3033-3036,2001 Tinsley B.A. and Dean G.W (1991) Apparent tropospheric response to MeV-GeV particle flux variations: a connection via electrofreeezing of supercooled water in high-level clouds? J Geophys Res 96, pp22283-22296 Tinsley, B.A., Rohrbaugh R.P., Hei M., and Beard K.V., (2000), Effects of image charges on the scavenging of aerosol particles by cloud droplets, and on droplet charging and possible ice nucleation processes, J. Atmos Sci., 57, 2118-2134. Tripathi S.N. and Harrison R.G. (2001), Scavenging of electrified radioactive aerosol, Atmos Environ, (in press) Vohra K.G., Subba Ramu M.C.and Muraleedharan T.S. (1984), An experimental study of the rôle of radon and its daughters in the conversion of sulphur dioxide into particles in the atmosphere, Atmos. Env., 18, 8, 1653-1656 Volland H. (1984), Atmospheric electrodynamics, Springer-Verlag, Berlin Wilson C.T.R. (1929), Some thundercloud problems, J. Franklin Institute, 208, 1-12 Yu F. and Turco R.P. (2001) From molecular clusters to nanoparticles: the rôle of ambient ionisation in tropospheric aerosol formation J. Geophys. Res. 106, D5, 4797-4814 TROPOSPHERIC ION MEASUREMENTS K L. Aplin* and R. G. Harrison Department of Meteorology, The University of Reading, PO Box 243, Earley Gate, Reading, RG6 6BB UK Abstract To investigate ion-induced nucleation in the atmosphere experimentally, it is necessary to select only the circumstances when ion-induced effects are likely to dominate, and identify and exclude times when other meteorological factors are influencing ion concentrations. To do this, reliable ion, aerosol and meteorological measurements at the same site are required, with analysis of different weather conditions. A methodology for identifying the effects of meteorological conditions on the local atmospheric electrical environment is discussed, based on ion and meteorological measurements made at Reading in the spring and summer of 2000. It is found that even under virtually identical synoptic meteorological conditions, there is significant variability in the ion concentration. 1. INTRODUCTION Atmospheric small ions are produced near the Earth's surface by natural radioactivity and cosmic rays; the electrical conductivity of the air s is proportional to the total ion concentration n. In equilibrium, the concentration of small ions is modulated by the ion production rate and the atmospheric aerosol concentration (e.g. MacGorman and Rust, 1998). Meteorological factors also have an important effect on s, for example, suspended aerosol particles are wind-borne, and may be advected to locally reduce the ion concentration by attachment. The links between meteorology and atmospheric electricity are tacitly acknowledged, but surprisingly little effort has been made to quantify the relationship between them. Chalmers' (1967) discussion was typical of much of the atmospheric electricity literature, clearly recognising meteorological effects on conductivity, but giving little explanation of the causal links. Atmospheric water drops carry an electric charge, and therefore conditions such as fog and haze can cause perturbations to atmospheric electrical measurements, both directly (electrically) and by condensation onto the insulators essential in measuring apparatus. Periods of stable atmospheric electrical conditions, however, are necessary for consistent measurements. In the classical paradigm of “fair-weather” atmospheric electricity, Ohm’s law relates the air conductivity s, vertical charge flux density J and vertical potential gradient E, as J = sE (1) (e.g. MacGorman and Rust, 1998). Conditions under which fair-weather properties can be expected were only summarised relatively recently by Reiter (1992), who excluded periods when hydrometeors were present at the surface. High cloud and fair-weather cumulus were permitted with this classification, until the cumulus started to become grey at the base indicating that it was beginning to charge. However, Barlow and Harrison (1999) showed experimentally that non- electrified clouds perturb the surface atmospheric electric field by thermal influences on the turbulent transport of charged particles and ions. * Now at Rutherford Appleton Laboratory, Space Science and Technology Department, Chilton Didcot, Berks OX11 OQX, UK. The need to investigate meteorological effects on atmospheric ion variability has recently become highly relevant for climate studies, following the published correlation between cosmic ray ionisation and clouds (Svensmark and Friis-Christensen, 1997). A theoretical mechanism linking ions and clouds has been described, involving the nucleation of atmospheric ultrafine aerosol onto atmospheric ions (Yu and Turco, 2001). Explicit observation of this effect in the atmosphere requires the distinction of ion-mediated nucleation from other factors affecting ions, and knowledge of favourable atmospheric conditions for ion-mediated processes. To do this, the effects of particular meteorological conditions on the atmospheric conductivity must be classified. 2. METEOROLOGY AND IONS 2.1 Expected effects of meteorological conditions on the ion concentration Some aspects of the diurnal variation of conductivity (s) (which is directly proportional to the ion concentration) under different meteorological conditions can be inferred. For example, the nocturnal inversion traps radioactive gases near the surface where they cause increased local ionisation and an associated higher conductivity. When the sun rises turbulence sharply increases, mixing the air and dispersing radioactive gases trapped in the surface layer. This “sunrise effect” (Chalmers, 1967) should be most pronounced on clear days, as there is a greater difference between the daytime and nocturnal atmospheric stability. Buoyant convection and the presence of clouds damping turbulence are already thought to influence the atmospheric potential gradient (Barlow and Harrison, 1999), well after sunrise. Smaller “sunrise effects”, showing similar behaviour, may occur throughout the day if the sun is in and out of cloud. Extending similar reasoning, s would be expected to vary more over a clear day than a completely cloudy day, as the former case has a greater diurnal variation in the solar radiation and associated convection. Convective turbulent mixing at the surface reduce ion concentrations to a minimum during the warmest part of the day, suggesting that on cloud-free days there should be an inverse relationship between s and solar radiation. On cloudy days mechanical generation of atmospheric turbulence dominates, with less direct dependence on the solar radiation. 2.2 Ion measurements at Reading Negative conductivity was measured using a self-calibrating instrument (Aplin and Harrison, 2001) at The University of Reading Meteorology Field Site from April-July 2000. Automatic meteorological measurements are made at this site at 1 Hz, and manual weather observations are made daily. Three days are selected for detailed study here: 13th, 17th and 18th June 2000, (year days 165, 169 and 170), for which synoptic charts are shown in Figure 1. These three days would all be traditionally classified as having fair-weather atmospheric electrical conditions. Days 169 and 170 were both characterised by high pressure and weak southerly flow. Day 170 was completely cloud-free and the hottest day of the year, with a maximum of 29.7 ¼C, whereas Day 165 was cooler, with intermittent cloud and westerly flow. a) b) c) Figure 1: Synoptic charts of the three selected days a) 13th June (day 165) b) 17th June (day 169) c) 18th June (day 170) (data from www.wetterzentrale.de). The conductivity was sampled at nominally two-minute intervals and processed as described in Aplin (2000); hourly averages are discussed here and are shown in Figure 2 below. Two cases have been selected for detailed analysis. The cloud-free and hot day 170 is compared with day 165, which had stratocumulus cloud with sunny intervals until 1600 followed by sunshine until 1930. The effect of aerosol on conductivity during cloud-free periods on the consecutive days 169 and 170 is also investigated. 40 165 169 170 30 20 10 negative conductivity (fS/m) 0 45678910111213141516171819202122 hour of day BST Figure 2: Hourly averaged diurnal variation of negative conductivity for 13th, 17th and 18th June 2000 2.3 Conductivity on a cloudy day (day 165) and cloud-free day (day 170) A sunrise effect is clearly apparent on day 170 with a peak at 0400 followed by a decrease in s as the solar radiation increases. The minimum in s occurs between 1600-1700, which coincides with the maximum temperature (1550). Over the whole day, s is negatively correlated to > 95% with the global solar radiation Sg. On day 165, the mean conductivity † (Chilbolton, 48 km SW of Reading) http://www.met.rdg.ac.uk/radar/ some cloud present at Reading from the radiometer measurements, but it is likely to have been thin or high cloud with little impact on the surface conditions. Despite the very similar weather on days 169 and 170 (see Table 1), the conductivities were different and did not show the same variation, with a correlation coefficient r = Ð0.26 between them. This suggests that other processes cause the variability, such as changes in the aerosol concentration. s could also be sensitive to small changes in meteorological variables such as the wind speed. PM10 mass concentrations M are measured at a site 2 km NW of the University field site‡, and are therefore used to infer bulk properties of the air mass. M was almost twice as high on day 169 as day 170, probably because of the increased traffic on a Saturday. The correlation coefficient between M and s was negative on day 169 and positive to > 98% on day 170. Ion-aerosol theory (e.g. Clement and Harrison, 1992) predicts a negative correlation between aerosol and s. Although the theory links aerosol number concentration and s, these observations show the expected general relationship between mass concentration and s. The lower s may therefore have been caused by higher local wind speeds associated with high aerosol concentrations. Table 1: Average conductivity, meteorological conditions, aerosol mass concentration and correlation (r) between M and s for 0530-1830 on 17th and 18th June 2000 -1 Day ‡ http://www.aeat.co.uk/netcen/airqual/ The mass concentration M (mass per unit volume) and the number concentration Z of a monodisperse aerosol population are related by Ê 4 ˆ MZ= Á pr r3 ˜ (4) Ë 3 ¯ where r is the particle density. This equation was used to estimate Z from M assuming the average radius of the population (2 mm) and that the particles are ammonium sulphate, with r = 1.77 gcm-3 (Khlystov et al, 2001). The coefficients a and q were assumed to be 1.6 x 10-6 cm3s-1 and 10 cm- 3s-1 respectively. b was calculated from Gunn’s (1954) expression, and is a function of the mean radius, temperature and mean charge j (assumed to be Ð1e). The aerosol was split into three size modes, nucleation with mean radius rn = 0.25 mm, accumulation with ra = 0.88 mm and coarse with rc = 5.6 mm (Seinfeld and Pandis, 1998). The average sensitivity of the conductivity to a 1% change in the aerosol number concentration in each mode was calculated from equation 5. For the initially assumed size distribution, (shown in Table 2) with 80% of the particles in the fine mode and only 5% in the coarsest size range, the conductivity is most sensitive to changes in the coarse mode. This is not what would intuitively be expected, because there are fewest particles in the coarse mode, and from the ion balance equation (2), the number of particles controls the conductivity, by attachment. Another factor affecting the availability of particles to attach is their surface area. This was tested by fixing the surface area-Z product, by increasing the number of particles in the nucleation and accumulation modes. The nucleation mode was most sensitive to changes in Z; therefore for typical aerosol size distributions, the surface area of the coarse aerosol dominates over the number. Table 2 Characteristics of the assumed aerosol distribution, and calculated sensitivities to a 1% change in aerosol number concentration. Mode Nucleation Accumulation Coarse Assumed fraction of particles 0.8 0.15 0.05 Mean radius (mm) 0.25 0.88 5.6 Sensitivity to 1% change in Z 0.13% 0.11% 0.22% Sensitivity to 1% change in Z 0.2% 0.0% 0.0% (surface area x Z fixed) One finding is that conductivity is relatively insensitive to changes in the aerosol population, which may be significant because of the discussion whether conductivity can be used as an indicator of urban aerosol pollution (see e.g. Aplin, 2000). In this case, changes in the aerosol concentration are not significant enough to account for the observed conductivity variation on day 169, therefore other sources of variability in the conductivity dominate. 3. CONCLUSIONS The conductivity of air in the surface layer is relatively insensitive to the aerosol concentration at the atmospheric concentrations observed; therefore most of its variability is likely to be due to meteorological factors. The different conductivities measured on two days which were a) almost identical synoptically and b) conformed to “fair weather” in the strictest classical sense, imply that the ion concentration is sensitive to small changes in meteorological variables. The traditional “fair-weather” classification, which aims to ensure consistent and comparable measurements, therefore appears inadequate. Furthermore, it does not allow for electro-meteorological interactions and excludes other sources of atmospheric electrical variability. The majority of atmospheric electricity and aerosol measurements made at the surface in England probably occur in non fair-weather conditions, for which there is little understanding of the variability in anything other than general terms. There is evidence that cloud affects surface atmospheric electrical conditions by modulating atmospheric stability. We propose an extension to the existing classification: (1) fair-weather: as Reiter’s (1992) definition, but cloud-free, with variability caused almost entirely by micrometeorological factors (2) semi-fair-weather: presenting similar atmospheric electrical conditions to (1), but identified primarily by meteorological stability criteria with no local charge generation, rather than solely the absence of electrified clouds and (3) non-fair-weather. Extension of the classical fair-weather paradigm to include data obtained on cloudy days, and more detailed investigation of micrometeorological effects on conductivity is important if the ion-induced effects hypothesised to contribute to aerosol production are to be unambiguously identified. ACKNOWLEDGEMENTS The experimental work was supported by the UK Natural Environment Research Council. KLA’s attendance at the IACI meeting was funded by the Environmental Sciences Department, University of Hertfordshire. REFERENCES Aplin K.L. (2000), Instrumentation for atmospheric ion measurements, PhD. Thesis, The University of Reading, UK Aplin K.L. and Harrison R.G. (2001), Rev. Sci. Instrum., 72, 8, in press Barlow J.F and Harrison R.G. (1999), In Christian H.J. (ed) Proceedings 11th International Conference on Atmospheric Electricity, Guntersville, Alabama 7th-11th June 1999 NASA/CP-1999-209261, 575-578 Chalmers J.A. (1967), Atmospheric Electricity, 2nd edition, Pergamon Press, Oxford Clement C.F. and Harrison R.G. (1992), J. Aerosol Sci. 23, 5, 481-504 Dhanorkar S. and Kamra A.K. (1997), J. Geophys. Res. 102, D25, 30147-30159 Dolezalek H. (ed.) (1974), Electrical processes in atmospheres, Springer Verlag, Darmstadt Gunn R. (1954), J. Meteorol., 11, 339 Khlystov A., Kos G. P. A., ten Brink H.M., Mirme A., Tuch T., Roth C. and Kreyling W.G. (2001), Atm. Env., 35, 11, 2045-2051 MacGorman D.R. and Rust W.D. (1998), The electrical nature of storms, OUP Reiter R. (1992), Phenomena in atmospheric and environmental electricity, Elsevier, Amsterdam Seinfeld J.H. and Pandis S.N. (1998), Atmospheric chemistry and physics, Wiley, New York Svensmark H. and Friis-Christensen E. (1997), J. Atmos. Solar-Terrestrial Phys, 59, 1225-1232 Yu F. and Turco R.P. (2001). J. Geophys. Res. 106, D5, 4797-4814. ATMOSPHERIC AEROSOLS: FORMATION AND GROWTH Markku Kulmala University of Helsinki, Department of Physical Sciences, Division of Atmospheric Sciences, P.O. Pox 64, FIN-00014, University of Helsinki, Finland 1. INTRODUCTION It is widely recognised that the increasing atmospheric concentrations of greenhouse gases such as carbon dioxide and methane can potentially drive a significant warming process of the earth’s climate. However, a topic of more recent attention is the possibility that increased atmospheric concentrations of aerosol particles might drive a significant radiative forcing process of the planet (see, for example, Charlson et al., 1992; and Charlson and Wigley, 1994). The increased aerosol concentrations are largely due to secondary particle production i.e. homogeneous nucleation from vapour precursors. The secondary aerosols have both natural and anthropogenic origin. Aerosol particles influence the climate by two distinct mechanisms: the direct reflection of solar radiation by aerosol particles, and the indirect increase in cloud reflectivity caused by enhanced number of cloud condensation nuclei. IPCC (1996) has reported that uncertainties in the estimation of direct and indirect aerosol effects on global climate are big (see Fig. 1.). These uncertainties arise largely from the limited information on the spatial and temporal distribution of aerosols and clouds. However, recently some progress has been made in evaluating the radiative effects of various aerosol components such as sulfate, organics, black carbon, sea-salt, and crustal species (Chuang et al., 1997; Haywood and Ramaswamy, 1998; Kaufman and Fraser, 1997; Winter and Chylek, 1997; Sokolik and Toon, 1996). Despite these efforts, substantial uncertainties still remain in quantifying the contribution from each source, particularly, for biogenic and natural emissions, including organic vapours. Without understanding the contribution of natural emissions of aerosols and particles to radiative forcing, we can never hope to accurately predict or understand the true effect of anthropogenic emissions. Among the key questions in reducing error bars are how aerosol particles are formed, how they will grow from clusters of a few molecules to CCN sizes (>100 nm) and how they will form cloud droplets. Once formed, clouds have a very extensive influence on the Earth's radiation budget through their albedo and greenhouse effects. With global warming, future cloud properties are likely to change due to the warmer and moister conditions, and possibly due to increased aerosol particle emissions from both primary (e.g. wind generated sea-spray) and secondary aerosols (from biogenically and anthropogenically influenced gas-to-particle conversion processes). Clouds are however rather crudely presented in global and regional climate models (GCM, RCM). Processes, such as nucleation, droplet activation during condensation, diffusive growth, droplet evaporation, droplet coalescense and conversion to raindrops, are very crudely taken into account in present-day atmospheric large-scale models. For example, we have recently shown the importance of aerosol formation and growth processes to CCN concentrations (Kulmala et al., 2000) as well as the effect of nitric acid and other semivolatile gases in influencing cloud formation processes, in particular, in enhancing the cloud droplet population, thereby increasing cloud reflectance (Kulmala et al., 1993; Laaksonen et al., 1997). The importance of including multi-component aerosol populations, and the dynamic feedback in the cloud forming processes, along with the importance of coupling chemical and physical processes in predicting cloud droplet populations have been illustrated by O’Dowd et al., (1999a; 1999b). Particle formation and growth in the atmosphere have recently received growing experimental and theoretical interest. Therefore, instrumental techniques for measuring concentrations of freshly formed particle have been developed, and particles with diameter of about 3 nm can be detected. These small particles have been found in large variety of environments: in the free troposphere (Clarke, 1992; Schröder and Ström, 1997; Raes et al., 1997), in the marine boundary layer (Covert et al., 1992; Hoppel et al., 1994; O’Dowd et al., 1998), in the vicinity of evaporating clouds (Hegg et al., 1991), in Arctic and Antarctic areas (Wiedensohler et al., 1996; Pirjola et al., 1998; O’Dowd et al., 1997), in urban areas and in stack plumes (Kerminen and Wexler, 1994; Kerminen and Wexler, 1996; Väkevä et al., 2000). Starting during the mid-nineties, aerosol formation and growth events have been observed also in forested areas e.g. over boreal forest in Finland (Mäkelä et al., 1997, 2000; Kulmala et al., 1998), and in other type of forests in Portugal (Kavouras et al., 1998), Greece (Kavouras et al., 1999), Canada (Leaitch et al., 1999), and in USA (Marti et al., 1997). In all these cases particle formation and growth events took place in remote forested areas, where the release of highly reactive volatile organic carbons (VOCs) from trees followed by a rapid oxidation to low volatile products, has to be considered as a potential source for nucleating vapours. 3 Tropospheric ) Halocarbons 2 aerosols - direct effect 2 N 2 O C H 4 Tropospheric Fossil 1 C O 2 Ozone fuel Solar soot Biomass Sulphate burning 0 Stratospheric -1 Ozone Tropospheric aerosols -indirect effect -2 Global mean radiative forcing (W /m Confidence level Very Very High L o w Low Low Very Very Low Low Low Low Figure 1: Estimates of globally and annually averaged anthropogenic radiative forcing (in Wm-2) due to the changes in concentrations of greenhouse gases and aerosols from pre-industrial times to present day and to natural changes in solar output from 1850 to present day (IPCC, 1996). Atmospheric aerosol particles in urban areas, on the other hand, cause the loss of visibility (e.g. Finlayson-Pitts and Pitts, 2000) and health effects (Dockery and Pope, 1994). Heavily industrialized areas suffer from pollution fogs (smogs) that are often related to coal burning and nowadays also to traffic. The most well-known example of such smogs is the London ”pea- souper” smog, which occurred every once in a while until the 50’s, when coal burning was forbidden. Besides visibility degradation, the London smog episodes caused serious health effects and ”excess deaths”. One significant part of health problems related to atmospheric aerosols and fog droplets, since particles having diameters less than 10 mm can penetrate deep into the respiratory system (Dockery and Pope, 1994). Recently, the effect of ultra-fine particles have been discussed and their local variations have been investigated (e.g. Buzorius et al., 1999). 2. AEROSOL DYNAMICS During the processes of formation and growth of atmospheric aerosols the aerosol dynamics, atmospheric chemistry and meteorology form a coupled system. The importance of atmospheric chemistry (e.g. Pirjola and Kulmala, 1998; Pirjola 1999) as well as meteorological conditions (Nilsson and Kulmala, 1998; Nilsson et al., 2000; Väkevä et al., 2000) on particle formation and growth have been demonstrated under tropospheric conditions. Although ternary nucleation of water-ammonia-sulphuric acid vapours (Korhonen et al., 1999) has shown to be able to explain atmospheric nucleation Ð i.e. formation of ~1 nm particles - in many cases (Kulmala et al., 2000), the exact routes for formation of 3 nm particles are still unclear, because besides nucleation, also the growth from 1 nm size to 3 nm size is needed. In order to be able to understand the formation and growth processes of atmospheric aerosols and cloud droplets their thermodynamic properties should be known. For example, in the condensation process, the driving force is the vapour pressure difference between gas phase and surface. However, in the atmosphere where there are multicomponent, multiphase mixtures, their thermodynamic state and phase diagrams are typically very complex. It is very important to obtain thermodynamically consistent vapour pressures, chemical activities, surface tensions and densities for organic and inorganic compounds and their water solutions (for the importance see e.g. Korhonen et al., 1999) as a function of temperature and composition. In future, development of nucleation theories, modelling and nucleation rate parameterizations are needed. So far, conclusions on whether or not certain substances cause nucleation in the atmosphere conditions are usually based on predictions given by the classical nucleation theory (CNT). CNT treats the nucleating molecular clusters as macroscopic droplets which is a questionable approach since the nucleating clusters often contain less than fifty molecules. Nucleation of various vapors using molecular dynamics (MD) and Monte Carlo (MC) simulation techniques is needed to investigate. So far, some investigations were carried out using ab initio calculations on small sulfuric acid-water clusters (Arstila et al 1998), classical MD (Laasonen et al, 2000) and MC (Vehkamäki and Ford, 1999) simulations of argon nucleation, as well as DFT calculations of nucleation in binary systems imitating water and different organic molecules (Laaksonen et al., 1995, Napari and Laaksonen 2000). Also, a new nucleation mechanism based on stable dimers (Lushnikov and Kulmala, 1998) has been proposed. In contrast to laboratory conditions, the formation of aerosol in the atmosphere can be kinetically limited by some of the intermediate steps of its formation processes. The equilibrium state is thus not necessarily the aerosol itself but can be, for example, thermodynamically stable clusters (TSC), as we have recently shown (Kulmala et al., 2000). Although there is strong indication that the water-sulphuric acid-ammonia nucleation mechanism (Korhonen et al., 1999) explains the formation of new atmospheric aerosols (diameter < 3 nm) in many circumstances, the condensation of these vapors does not explain the observed growth rates of the particles (Kulmala et al., 2000), and in atmospheric conditions nucleation and growth are decoupled (Kulmala et al., 2000). The other possible relevant nucleation mechanism is ion-induced nucleation. Aerosol dynamic modelling (nucleation, condensation, coagulation, deposition) with gas phase chemistry to obtain the atmospheric significance of nucleation and condensation of different vapours have been and will be performed. The aerosol dynamics and atmospheric chemistry model used in the present research is based on the model recently developed by our research group (Pirjola and Kulmala, 1998; Pirjola 1999). In these models aerosol formation and growth including aerosol dynamics to evaluate sink terms for condensable molecules and gas phase chemistry to include source terms for these molecules will be used. Process models will be coupled with dispersion models. In the chemistry part of the model the chemistry of O3, NOx, VOC and other relevant species will be related to aerosol formation. The effects of meteorological dynamics on aerosol processes will be studied by applying the aerosol dynamic models in a Lagrangian approach including wave motions and atmospheric mixing. The results shows that ternary water-ammonia-sulphuric acid system is proper candidate for atmospheric aerosol formation. 3. FORMATION AND GROWTH OF ATMOSPHERIC AEROSOLS, FIELD EXPERIMENTS Formation and growth of aerosol particles have been observed and will be observed at atmospheric conditions. Our research group has participated in several field campaigns. These includes continuous measurements performed at our field stations and several international intensive campaigns like Aerosol Characterisation Experiment 1 and 2 in 1994 and 1997 (ACE-1 and ACE- 2 organised by IGAC), International Arctic Ocean Expeditions 1991 and 1996, Biogenic aerosol formation in the boreal forest (BIOFOR, 1997-1999, SMEAR stations, Finland, Hyytiälä), New particle formation and fate in the coastal environment (PARFORCE, Mace Head, Ireland), and ongoing the OSOA (Origin and Formation of Secondary Organic Aerosol) experiment. As an example we consider here BIOFOR results in more detailed. All data measured during the BIOFOR campaigns are available on the Biofor web pages http://mist.helsinki.fi/Biofor/index.html (ask for usercode and password from the corresponding author). In addition to the numerical data there are also a number of plots produced as a result of the analysis of the data. The data are classified into 9 subgroups: 1) aerosol total number concentration and size distribution measurements in the size range 3-800 nm, 2) aerosol chemistry, 3) aerosol and gas fluxes by eddy covariance and gradient methods, 4) measurements of meteorological parameters and gas concentrations at six different levels from the mast, 5) meteorology of boundary layer and trajectories, 6) concentrations and emissions of BVOC (biological volatile organic compounds), 7) ground level concentrations of inorganic gases, 8) measurements of the size distribution of wet (ambient) aerosol from 0.5-32 _m at 18 m height, and 9) solar radiation measurements. The detailed descriptions of the instruments used are given on the web pages. When the particle formation event occurs, the mode of the fresh particles appears into the measurement range. In Figure 2. aerosol number size distributions measured using Differential Mobility Particle Sizer (DMPS) during a typical nucleation event day are shown. The nucleation mode practically dominates the spectrum with its high number concentration during the nucleation burst. For this event, particle growth from nucleation mode up to accumulation mode is clearly observable. The growth is frequently seen to continue during the following days up to accumulation mode (see also Kulmala et al., 2001). 105 104 103 00 03 06 09 12 15 18 21 00 Time of day (hr) Figure 2: Aerosol size distributions measured by DMPS from 2m height inside the forest (6.4. 1999). 105 TSI3010 67m TSI3025 67m 104 103 00 03 06 09 12 15 18 21 00 Time of day (hr) Figure 3: Aerosol concentrations measured by CPC’s above the forest heights 67m and 18m (6.4. 1999). During the events, aerosol fluxes determined using an eddy covariance technique are observed to be downwards. Also the measurements made by Condensation Particle Counters (CPC) and DMPS at different heights support this finding. From particle flux data, using the eddy covariance method (Buzorius et al., 1998), usually a small overall downward flux is observed. The downward flux clearly increases during nucleation events, with an exception of the cases when the surface wind was from direction of 220-250∞ (direction of the Tampere city and the Hyytiälä institute buildings). Then a strong upward particle flux is observed due to local surface-level pollution. The difficulties in absolute calibration of the DMPS set ups as well as sampling losses in the lines suggested that the gradient of particles will be best determined placing two identical CPC pairs in the mast (18 m and 67 m height). The CPC pairs consisted of the ultrafine CPC (TSI Inc 3025) for determination of the particles larger than 3 nm in diameter and conventional CPC (TSI Inc. 3010) for particles larger than 10 nm in diameter. The difference of the reading of the CPC’s gives an approximate value for the ultrafine mode particle concentrations in the beginning of the burst. The data from the CPC pairs is shown for the event day of 6 April 1999 in Figure 3. The difference between the CPC readings from the two levels shows that the ultrafine particles have higher concentrations in higher level during the nucleation burst. This result will support the particle flux data that illustrate a net loss of particles to the canopy; however, it does not necessarily indicate a particle source at the top of the boundary layer or higher altitudes, even though nucleation is more probable in these regions. 4. CONCLUSIONS According to recent results on atmospheric aerosol formation some preliminary conclusions can be made on atmosphere aerosol formation. (see Kulmala et al., 2001) The most probable formation mechanism is ternary nucleation (water Ð sulphuric acid Ð ammonia) and the growth to observable sizes takes place mainly owing to condensation of organic vapours. Nevertheless, there is no direct proof of this phenomenon because the composition of 1Ð5 nm size particles is very difficult to determine using present state-of-art instrumentation. The other possible nucleation mechanism is ion induced nucleation with sulphuric acid and water vapours. Nucleation takes place typically in very specific weather conditions: e.g. in Hyytiälä in cold air advection in Polar and Arctic air masses, at low cloudiness, and no precipitation. Furthermore, the nucleation was closely connected to the onset of strong turbulence in the morning-noon transition from stable to unstable stratification, which should also correspond to the onset of convection and entrainment from aloft. The emission rates for several gaseous compounds have been determined (Kulmala et al., 2001). Using four independent ways the amount of the condensable vapour needed for observed growth of aerosol particles was estimated to 2-10 x 107 vapour molecules cm-3. The estimations for source rate gives 7.5-11 x 104 cm-3s-1. REFERENCES Arstila H., Laasonen K. and Laaksonen A, (1998) J. Chem. Phys., 108, 1031. Buzorius, G., U. Rannik, J.M. Mäkelä, T. Vesala and M. Kulmala (1998) J. Aerosol. Sci. 29, 157. Buzorius, G., Hämeri, K., Pekkanen, J. and Kulmala M. (1999) Atmospheric Environment 33, 553. Buzorius, G., Ü. Rannik, J. M. Mäkelä, P. Keronen, T. Vesala, M. Kulmala. (2000) J. Geophys. Res. 105, 19905. Charlson, R.J., and Wigley, T.M.L. 1994. Sulphate aerosol and climatic change. Scientific American, 270, 48-57. Charlson, R.J., Schwartz, S.E., Hales, J.M., Cess, R.D., Coakley, J.A., Hansen, J.E. and Hofmann, D.J. 1992. Climate forcing by anthropogenic aerosols. Science, 255, 423-430. Clarke, A.D. 1992. Atmospheric nuclei in the remote free troposphere. J. Atmos. Chem., 14, 479- 488. Covert, D.S., Kapustin, V.N., Quinn, P.K., and Bates, T.S. 1992. New particle formation in the marine boundary layer. J. Geophys. Res., 97, 20581 Ð 20587. Chuang C. C., Penner J. E., Taylor K. E., Grossman A. S. and Walton J. J. (1997). J. Geophys. Res. 102, 3761. Dickerson, R.R., Kondragunta, S., Stenchikov, G., Civerolo, K.L., Doddridge, B.G. and Holben, B.N. (1997) Science 278, 827. Dockery, D.W. and Pope, C.A. (1994) Annu. Rev. Public. Health, 15, 107. Finlayson-Pitts B.J. and Pitts, Jr. J.N. (2000) Chemistry of the upper and lower atmosphere. Academic Press, San Diego. Hämeri, K., Väkevä, M., Hansson, H.-C. and Laaksonen A. (2000) J. Geophys. Res. 105, 22231. Haywood J. M. and Ramaswamy V. (1998) J. Geophys. Res. 103, 6043. Hegg, D.A., Radke, L.F. and Hobbs, P.V. 1991. Measurements of Aitken nuclei and cloud condensation nuclei in the marine atmosphere and their relation to the DMS-cloud-climate hypothesis. J. Geophys. Res., 96, 18727-18733. Hoppel, W.A., Frick, G.M., Fitzgerand, J.W. and Larson, R.E. 1994. Marine boundary layer measurements of new particle formation and the effects nonprecipitating clouds have on aerosol size distribution. J. Geophys. Res., 99, 14443-14459. IPCC (Intergovernmental Panel on Climate Change) 1996. Climate Change 1995 (Edited by Houghton J.T., Meira Filho L.G., Gallander B.A:, Harris N., Kattenberg A. and Maskell K.). Cambridge University Press, Cambridge Kaufman Y. J. and Fraser R. S. (1997) Science 277, 1636. Kavouras, I.G., Mihalopoulos, N. and Stephanou, E.G. 1998. Formation of atmospheric particles from organic acids produced by forests. Nature, 395, 683-686 Korhonen, P., Kulmala, M., Laaksonen, A., Viisanen, Y., McGraw, R. and Seinfeld, J.H. 1999. Ternary nucleation of H2SO4, NH3 and H2O in the atmosphere. J. Geophys. Res., 104, 26,349-26,353. Kulmala, M., A. Laaksonen, P. Korhonen, T. Vesala, T. Ahonen, J.C. Barrett, (1993) J. Geophys. Res. 98, 22949. Kulmala, M., Laaksonen A., Charlson R.J. and Korhonen P. (1997) Nature, 388, 336. Kulmala, M., A. Toivonen, J. M. Mäkelä and A. Laaksonen, 1998. Analysis of the growth of nucleation mode particles observed in Boreal forest. Tellus, 50B, 449-462. Kulmala, M., Pirjola, L. and Mäkelä, J.M. 2000. Stable sulphate clusters as a source of new atmospheric particles. Nature, 404, 66-69. Kulmala, M., K. Hämeri, P.P. Aalto, J.M. Mäkelä, L. Pirjola, E.D. Nilsson, G. Buzorius, U. Rannik, M. dal Maso, W. Seidl, T. Hoffman, R. Janson, H.-C. Hansson, Y. Viisanen, A. Laaksonen, and C.D. O’Dowd. 2001. Overviwe of the internal project on biogenic aerosol formation in the boreal forest (BIOFOR). Tellus 53B, 324-343. Laaksonen, A. Talanquer, V., and Oxtoby, D.W (1995) Annu. Rev. Phys. Chem., 46, 489. Laaksonen, J. Hienola, M. Kulmala, F. Arnold, (1997) Geophys. Res. Lett. 24, 3009. Laaksonen, P. Korhonen, M. Kulmala, R.J. Charlson, (1998) J. Atmos. Sci. 55, 853. Laasonen, K., , S. Wonzack, R. Strey and A. Laaksonen, J. Chem. Phys. (in press) Leaitch, W.R., Bottenheim, J.W., Biesenthal, T.A., Li, S.M., Liu, P.S.K., Asalian, K., Dryfhoutclark, H., Hopper, F. and Brechtel, F. 1999. A case study of gas-to-particle conversion in an eastern Canadian forest. J. Geophys. Res. 104, 8095-8111. Lushnikov, A.A., Kulmala, M. (1998) Phys. Rev. Lett. 81, 5165. Marti, J.J., Weber, R.J., McMyrry, P.H., Eisele, F., Tanner, D., and Jefferson A. 1997. New particle formation at a remote continental site: Assessing the contributions of SO2 and organic precursors. J. Geophys. Res. 102, 6331-6339. Mäkelä, J.M., Aalto, P., Jokinen, V., Pohja, T., Nissinen, A., Palmroth, S., Markkanen, T., Seitsonen, K., Lihavainen, H., and Kulmala, M. 1997. Observations of ultrafine aerosol particle formation and growth in boreal forest, Geophys. Res. Lett., 24, 1219-1222. Mäkelä, J.M., Koponen, I.K., Aalto, P. and Kulmala, M. 2000. One-year data of submicron size modes of tropospheric background aerosol in Southern Finland. J. Aerosol Sci., 31, 595- 611. Napari, A. Laaksonen, L. (2000) Phys. Rev. Lett. 84, 2184. Nilsson, E.D. and Kulmala, M. 1998. The potential for atmospheric mixing processes to enhance the binary nucleation rate. J. Geophys. Res., 103, 1381-1389. Nilsson, E.D., Pirjola, L. and Kulmala, M. 2000. The effect of atmospheric waves on aerosol nucleation and size distribution. J. Geophys. Res., 105, 19,917-19,926. O’Dowd, C.D., Lowe, J.A., Smith, M.H., Davison, B., Hewitt, C.N. and Harrison, R.M. 1997. Biogenic sulphur emissions and inferred non-sea-salt-sulphate cloud condensation nuclei in and around Antarctica, J. Geophys. Res., 102, 12,839-12,854. O’Dowd, C.D., Geever, M., Hill, M.K., Jennings, S.K., and Smith, M.K. 1998. New particle formation and spatial scales in the clean marine coastal environment. Geophys. Res. Lett., 25, 1661-1664 O’Dowd, C. D., J. A. Lowe, M. H. Smith and A. D. Kaye (1999a) Q. J. R.. Meteorol. Soc., 125, 1295. O'Dowd, CD, J.A. Lowe, and MH Smith (1999b) Geophys. Res. Letts, 26, 1311. Seckmeyer, G., Mayer, B., Bernhard, G., Erb, R., Albold, A., Jäger, H. and Stockwell, W.R. (1997) Atmospheric Environment 34, 2971. Pirjola, L. 1999. Effects of the increased UV radiation and biogenic VOC emissions on ultrafine sulphate aerosol formation. J. Aerosol Sci., 30, 355-367 Pirjola, L. and Kulmala, M. 1998. Modelling the formation of H2SO4-H2O particles in rural, urban and marine conditions. Atmos. Res., 46, 321-347. Pirjola, L., Laaksonen, A., Aalto, P., and Kulmala, M. 1998. Sulfate aerosol formation in the Arctic boundary layer. J. Geophys. Res., 103, 8309-8322. Raes, F., Van Dingenen, R., Cuevas, E., Van Velthoven, P.F.J., and Prospero, J.M. 1997. Observations of aerosols in the free troposphere and marine boundary layer of the subtropical Northeast Atlantic: Discussion of processes determining their size distribution, J. Geophys. Res., 102, 21315-21328. Schröder, F. and Ström, J. 1997. Aircraft measurements of submicrometer aerosol particles (> 7 nm) in the midlatitude free troposphere and tropopause region. Atmos. Res., 44, 333-356. Seinfeld, J.H. and Pandis (1998) Atmospheric Chemistry and Physics. (Wiley) Sokolik I. N. and Toon O. B. (1996) Nature 381, 681. Vehkamäki, H., and I. J. Ford, (1999) J. Chem. Phys., 112, 4193. Viisanen, Y. (1991) Ph.D.Thesis, University of Helsinki, Department of Physics. Viisanen, Y., Strey, R., Reiss, H., (1993) J. Chem. Phys., 99, 4680. Viisanen, Y., Strey, R., Reiss, H. (2000) J. Chem. Phys., 112, 8205 Viisanen, Y., Strey, R. (1996) J. Chem. Phys., 105, 8293 Viisanen, Y., Kulmala, M., Laaksonen, A., (1997) J. Chem. Phys., 107, 920 Väkevä, M., Hämeri, K., Puhakka, T., Nilsson, E.D., Hohti, H. and Mäkelä, J.M. 2000. Effects of meteorological processes on aerosol particle size distribution in an urban background area, J. Geophys. Res., 105, 9807-9821. Winter B. and Chylek P. (1997)Tellus 49B, 72. COSMIC RAYS, PARTICLE FORMATION, NATURAL VARIABILITY OF GLOBAL CLOUDINESS, AND CLIMATE IMPLICATIONS Fangqun Yu Atmospheric Sciences Research Center, State University of New York, Albany, New York, USA Abstract Via its role in aerosol formation, cosmic ray may affect the global cloudiness and hence climate. Here we show that an increase in cosmic ray fluxes may lead to an increase in particle production in the lower troposphere but a decrease in particle production in the upper troposphere. In addition to the reported positive correlation between cosmic ray variations and low cloudiness, our analysis of satellite-based cloud cover data reveals that high cloudiness may be anti-correlated with cosmic ray variations if volcano and El Niño impacts are excluded. The observed different correlations between cosmic ray variations and low, middle and high cloud anomalies are consistent with the predicted different sensitivities of particle production to cosmic ray changes at various altitudes. The influence of the solar-modulated cosmic ray fluxes on global cloudiness, if confirmed, may provide the external forcing needed to reconcile the apparent differences between observed surface and troposphere temperature trends. 1. INTRODUCTION Clouds play a key role in the energy budget of Earth’s surface and lower atmosphere, and are probably the largest contributor to the uncertainty concerning the global climate change1. Small modifications of the amount, distribution, or radiative properties of clouds can have significant impacts on the predicted climate2. To detect and attribute anthropogenic influences on climate, it is crucial to quantify the natural fluctuations of cloudiness and the associated radiative forcing. In 1997, Svensmark and Friis-Christensen3 reported a surprising discovery that total cloud cover over midlatitude ocean correlates closely with the galactic cosmic ray (GCR) intensity. The cloud data analyzed include the C2 data sets from the International Satellite Cloud Climatology Project (ISCCP)4. Recently it has become possible to infer global cloud properties at different altitudes from the ISCCP-D2 data, which come from an improvement of procedures leading to the C2 data5. Analyses of the ISCCP-D2 data indicate that a clear correlation can only be seen between GCR fluxes and the global average of low cloud cover6,7. Due to its potential importance and implication, the GCR-cloud-climate hypothesis3 has been under close scrutiny8-10. Two of the main questions and doubts raised against the hypothesis are: (1) no convincing physical mechanism is available to explain the correlation, (2) there is no obvious correlation between solar activity and high cloudiness (where, it is argued, if GCR ionization has any impact on cloud microphysics, it would most likely be found in the upper troposphere where GCR incidence is greatest). Here we first try to address these two raised issues by investigating the role of GCR ionization in particle formation and the potential altitude-dependent influence of GCR variations on particle production and global cloudiness. We then explore the possibility of GCR-induced global cloud changes as an external forcing that may reconcile the apparent differences in global mean temperature trends between ground and atmosphere measurements. Over the last two decades, the temperature records taken at the Earth's surface show rapid warming (globally 0.15 ± 0.05 oC per decade), however the data produced by satellite and balloon studies indicate little if o any warming (globally 0.05 ± 0.10 C per decade) of the low to mid-troposphere - the atmospheric layer extending up to about 8 km from the Earth's surface11,12. Climate models generally predict that this atmospheric layer should warm faster than the surface if increased concentrations of greenhouse gases are causing the warming. Model simulations taking into account the effects of sulfate aerosols, stratospheric ozone depletion, and volcano eruptions were not able to reconcile these inconsistencies13-15. Gaffen et al16 suggested that these inconsistencies may be associated with external forcings of climate system that result in different surface and low tropospheric temperature changes. It is of interest and importance to investigate if the GCR- induced global cloud changes can provide the kind of external forcing needed to reconcile the inconsistencies. 2. GCR-CN-CCN-CLOUD HYPOTHESIS GCR-CN-CCN-Cloud Hypothesis GCR CN CCN Cloud GCR Variations Production rate Abundance Properties Ions Clusters CN CCN Cloud Droplets Ion Condensation Condensation Activation Nucleation Coagulation Coagulation H2SO4 H2SO4, Organics Ions, H2SO4 (?? NH3, Organics) H2O SO2 (?? NH3, HNO3, HCl) DMS Figure 1. Schematic illustrating of GCR-CN-CCN-Cloud Hypothesis that, if confirmed, might explain the correlation between variations of GCR flux and low cloud cover. The possible dominating species involved in the different phases of CN formation and growth processes are also indicated. The organics species may play an important role in growing the CN into the size of CCN. Figure 1 shows the GCR-CN-CCN-Cloud hypothesis that, if confirmed, might offer a physically- based link between GCR fluxes and global low-level cloud properties. Several steps are involved in this hypothesis. First, the modulation of galactic cosmic radiation by the solar cycle will cause a notable variation in aerosol production and condensation nuclei (CN) population in the lower atmosphere. Second, a systematic change in the ultrafine particle production rate will affect the population of cloud condensation nuclei (CCN). Third, a change in CCN abundance will affect the cloud properties. Clouds that form in air containing high CCN concentrations tend to have high droplet concentrations, which leads to both an increase in albedo and an increase in absorption. Increase in the CCN concentration also inhibits rainfall and therefore increases cloud lifetimes (cloud coverage). These effects Ð which are due to more, smaller droplets at fixed liquid water content Ð are particularly significant in marine air, where the CCN concentrations are generally quite low. In our proposed hypothesis, the first key process connecting GCR flux and low cloud cover is that ions generated by GCR ionization play an important role in new particle formation in the lower atmosphere which is the focus of discussion of the next section. 3. COSMIC RAY IONIZATION AND PARTICLE FORMATION Ambient ions are continuously generated by galactic cosmic rays at the rate of ~2 ion-pairs cm-3s-1 at ground level and up to ~20-30 ion-pairs cm-3s-1 in the upper troposphere17,18. Due to enhanced growth and stability of charged clusters (as a consequence of electrostatic interactions), air ions produced from GCR ionization may play an important role in the production of new particles under typical tropospheric conditions19,20. The proposed ion-mediated nucleation (IMN) theory can physically explain the enhanced growth rate (a factor of ~ 10) of sub-nanometer clusters and 2 the square of sulfuric vapor concentration ([H2SO4] ) dependence of nucleation rate as observed by Weber et al21, and seems to account consistently for ultrafine aerosol formation in jet plumes, in clean continental air and in marine boundary layer, as well as for the diurnal variation in the atmospheric mobility spectrum, as demonstrated by Yu and his colleagues 19,20,22-24. It is generally known that sulfuric acid vapor concentration ([H2SO4]), temperature (T), relative humidity (RH), pressure (P), and the surface area of preexisting particles are among the list of parameters controlling the particle formation in the troposphere. The IMN theory adds another important parameter-ion concentration ([ion], or ionization rate Q)-to this list. Here we focus on investigating the influence of GCR variations on particle formation and CN abundance at different altitudes. We employ an advanced particle microphysics (APM) model that simulates a size-resolved multicomponent aerosol system via a unified collisional mechanism involving both neutral and charged particles down to molecular sizes20. The size-resolved ion-ion recombination coefficients, ion-neutral collision kernels, and neutral-neutral interaction coefficients calculated in the model are physically consistent and naturally altitude (temperature, pressure, and relative humidity) dependent20. For the simulations presented below, the ion concentration is initialized as Q /a where a is ion-ion recombination coefficient. The pre-existing particles are initialized as two log-normal modes with total number densities of 19.5/cm3 and 0.6/cm3, median dry diameters of 0.09 mm and 0.3 mm, and standard deviations of 1.6 and 1.5, respectively. This gives an initial 2 3 wet surface area of ~ 4.5 mm /cm at 90% relative humidity, corresponding to a cloud-processed clean air mass where typical significant aerosol nucleation has been observed. Figure 2 shows the total condensation nuclei bigger than 3nm (Nd>3 nm) after three hours of simulations as a function of ionization rates at three different altitudes (0, 5, 8 km). The values of [H2SO4], T, P, and RH for each altitude (as specified in the figure legend) are fixed during the three-hour simulations. The shaded areas in Figure 1 are the ranges of Q values corresponding to low (>680 mb), middle (440-680 mb), and high (<440 mb) cloud regions as defined in ISCCP cloud data according to the cloud top pressures. It is clear from Figure 2 that significant number of ultrafine particles have formed under all the considered conditions. Most of these newly formed particles began as electrically charged clusters that have the advantage of enhanced growth and stability due to electrostatic effects. The neutral sub-critical clusters, on the other hand, grow too slowly to exceed the critical size under the prevailing conditions. The production rate of ultrafine particles is most sensitive to [H2SO4] and [ion] (or ionization rate). [H2SO4] controls the growth rate of ion clusters, while [ion] determines the lifetime of charged clusters and the availability of ions. The neutralization by ion-ion recombination will make the growing charged clusters lose their growth advantage and the resulting neutral clusters may dissociate if smaller than the critical size. At typical [H2SO4] where nucleation has been observed, for very low Q most of the ion clusters have sufficient time to reach the larger stable sizes prior to recombination and the nucleation rate is limited by Q. As Q increases, ion concentration increases but the lifetime of ions decreases and hence the fraction of ions having sufficient time to grow to the larger stable sizes decreases. As a result, the total number of particles nucleated first increases but later on decreases as Q increases. Figure 2 demonstrates that, as Q increases, Nd>3 nm increases rapidly in the low cloud region but decreases in the high cloud region. The Q value at turning point (i.e., dN/dQ=0) is sensitive to [H2SO4] and is most likely located in middle cloud region. 7 3 0 km (T=288 K, P=1013 mb, RH=90%) with [H2 SO 4 ]= 2.0x10 / cm 7 3 5 km (T=256 K, P= 541 mb, RH=77%) with [H2 SO 4 ]= 1.5x10 / cm 7 3 8 km (T=236 K, P= 357 mb, RH=68%) with [H2 SO 4 ]= 1.1x10 / cm 20000 18000 16000 5 km 0 km 14000 ) 3 12000 (#/cm 10000 d>3 nm N 8000 6000 8 km 4000 Low Middle High 2000 051015 20 25 30 Ionization rate Q (ion-pairs/cm3) Figure 2. Simulated concentrations of total condensation nuclei larger than 3nm (Nd>3 nm) after three hour of simulations for various ionization rates (Q) at three altitudes (0, 5, and 8 km). The shaded areas are the ranges of Q corresponding to low (>680 mb), middle (440-680 mb), and high (<440 mb) cloud regions as defined in ISCCP cloud data. Nd>3 nm increases rapidly in the low cloud region but decreases rapidly in the high cloud region as Q increases. During a solar cycle, the values of Q vary by ~20-25% in the upper troposphere and ~5- 10% in the lower troposphere. To study the effect of such systematic change of ionization rates on particle production at different altitudes, we increase the baseline ionization rate at each chosen altitude by 20% and compare the CN abundance after three hours of simulations. The altitude- dependent values of [H2SO4], Q, T, RH, P, and the surface area of preexisting particles are specified and some of them are shown in Figure 3. The baseline values of Q at different altitudes are from observations17,18, and the temperature and pressure are according to the US standard atmosphere. The [H2SO4] and RH are parameterized in a way so that they are constant in the lowest 2 km of atmosphere (2x107/cm3 and 90%, respectively) and gradually decrease with altitude above 2 km. These parameterizations are reasonable and are within the range of the observed values in various field campaigns25,26. Figure 4 shows the total condensation nuclei bigger than 3nm (Nd>3 nm) after three hours of simulations at different altitudes. The black line (with opaque circles) is for the baseline Q values while the green line (with filled circles) is for Q values 20% over the corresponding baseline values. The shaded areas in Figure 4 are low, middle, and high cloud regions as defined in ISCCP cloud data. [H2SO4], Q, T, and RH at each altitude (see Figure 3) are fixed during the three-hour simulations. It is clear from Figure 3 that an increase in GCR ionization rate associated with solar activity leads to an increase in the ultrafine production rate (i.e., dN/dQ>0) in the lower troposphere (as indicated by the red arrows) but a decrease in the ultrafine production rate (i.e., dN/dQ<0) in the upper troposphere (as indicated by the blue arrows). In the middle troposphere, dN/dQ changes sign and the average value of dN/dQ is small compared to that of lower and upper troposphere. It is interesting to note that the optimum particle formation layer is located in the middle troposphere (3-5 km altitude, likely in cloud outflows or top of low clouds), which is consistent with the measurements obtained in recent field campaigns such as ACE-126. 3 T (K) H2SO4 (#/cm ) 190 210 230 250 270 290 5x106 1077 1.5x10 2x107 2.5x107 12 12 11 11 H2SO4 10 10 9 9 8 T 8 Altitude (km) 7 7 6 6 5 RH 5 Altitude (km) 4 4 3 3 2 Surface Area/5 Q 2 1 1 0 0 0.5 0.6 0.7 0.8 0.9 1.0 051015 20 25 30 35 40 2 3 3 Relative Humidity, Surface Area/5 ( µm /cm ) Q (ion-pairs/cm s) Figure 3. The vertical profiles of [H2SO4], baseline Q, T, and RH used in the model to study the effect of a systematic change in ionization rates on particle production at different altitudes. 12 11 10 9 High 8 7 440 mb 6 5 Middle Altitude (km) 4 3 680 mb 2 Low 1 0 0 5000 10000 15000 3 Nd>3 nm (#/cm ) Figure 4. Simulated concentrations of total condensation nuclei larger than 3nm after three hours of simulations at different altitudes. The black line is for the baseline Q values while the green line is for Q values 20% over corresponding baseline values. The arrows indicate the changes in Nd>3 nm as ionization rates increase by 20%. The shaded areas are ISCCP low, middle, and high cloud regions. 4. NATURAL VARIABILITY OF GLOBAL CLOUDINESS It is well known that the abundance of cloud condensation nuclei (CCN) affects cloud properties27- 29. Clouds that form in air containing high CCN concentrations tend to have high droplet concentrations, which lead to an increase in both cloud albedo and absorption. Increases in the CCN concentration also inhibit rainfall and therefore increase cloud lifetimes (cloud coverage). Since the dominating number of CCN is evolved from newly formed ultrafine particles, a systematic change in the ultrafine particle production rate will affect the population of CCN. It is physically plausible that an increase in ultrafine production rate will increase the CCN abundance and cloudiness. During a solar cycle, the values of Q vary by ~20-25% in the upper troposphere and ~5-10% in the lower troposphere. Therefore, based on the influence of GCR ionization change on particle formation rate at different altitudes as shown in Figure 4, we can expect that if GCR variations have any impact on cloudiness, they should correlate positively with low cloud amount and negatively with high cloud amount. For middle clouds, such a correlation (if any) is likely to be weak. With these insights, we analyze the ISCCPÐD2 cloud data sets to study if the expected different correlations between GCR fluxes and low, middle, and high cloudiness exist. ISCCP-D2 data sets are considered as the most reliable measure of global cloud cover30 and are widely used for diagnostic studies of the climate system31 as well as verification of climate model simulations32. We look into the infra-red (IR) cloud data because they provide spatially and temporal unrestricted measurements that include clouds over the entire globe during both day and night6,7. Figure 5 shows the global average monthly mean anomalies of (a) high, (b) middle, and (c) low IR cloud cover during last solar cycle. To smooth out the seasonal variations, the monthly anomalies are calculated by subtracting the climatic monthly average from each month on an equal area grid before averaging over the globe6,7. The variations of GCR fluxes as measured from CLIMAX (normalized to May, 1965) are also indicated in each panel (dot-dashed lines). It is very clear that the low cloud anomalies highly co-vary with the change of GCR fluxes as has been reported by Marsh and Svensmark6,7. During one solar cycle, the absolute amount of low cloudiness changes by ~1.5-2%. The fluctuation of middle cloud anomalies is small compared to that of low cloud, and no obvious correlation exists between middle cloudiness and GCR variations. For the high cloud anomalies, there is no obvious correlation for the whole solar cycle. There may have several explanations for this. First, it takes much longer time for new particles to grow to the size of CCN or ice nuclei (IN) in high altitude than in low altitude due to much lower precursor vapor concentrations. As a result, the initial difference in CN production rate may not lead to obvious difference in CCN/IN abundance as a result of coagulation, scavenging, and mixing. Second, the properties of high cloud are determined by ice nuclei abundance which may be insensitive to CN production rate. The processes controlling IN abundance in high altitude are currently not well known. Third, there may exit a negative correlation but it does not appear in the ISCCP-D2 data of last solar cycle because of the influence of other processes such as volcano eruptions and El Niño events. We note that there were two major volcano eruptions during the period (El Chichón in April 1982 and Mt Pinatubo in June 1991). Volcano eruptions can inject large amount of SO2 into the stratosphere which leads to the formation of sulfate aerosols. On one hand, the cooling of upper troposphere as a result of volcano eruption may enhance the high cloud formation. On the other hand, the volcano aerosols descending from the stratosphere to the upper troposphere are likely to increase the frequency and lifetime of cirrus clouds33-35 and hence high cloudiness. The timescale to disperse the volcanic stratospheric aerosols around the whole globe through meridional circulations is 1-2 years36-38. Therefore, the effect of volcanic eruptions on global high cloudiness may become obvious 1-2 years after the eruptions. This is consistent with the observed increase in high cloudiness in 1984 and 1993 (i.e., 1-2 years after the El Chichón and Mt Pinatubo eruptions). A detailed analysis of Stratospheric Aerosol and Gas Experiment (SAGE) I and II aerosol extinction data for the upper troposphere39,40 indicates that a substantial enhancement of aerosols down to 2-3 km below the tropopause persisted until 1986 for the El Chichón eruption (i.e., ~4 years after the eruptions). The high cloudiness in 1987 may have been affected by the El Niño event during that year41. El Chichón Mt Pinatubo 1.5 High cloud GCR flux variations (%) 1 0 0.5 -10 0 -0.5 -20 -1 Change in cloudiness (%) -1.5 (a) -30 1.5 Middle cloud 1 0 GCR flux variations (%) 0.5 -10 0 -0.5 -20 -1 Change in cloudiness (%) -1.5 (b) -30 1.5 GCR flux variations (%) 1 Low cloud 0 0.5 -10 0 -0.5 -20 -1 Change in cloudiness (%) -1.5 (c) -30 1982 1984 1986 1988 1990 1992 1994 1996 Year Figure 5. The global average monthly mean anomalies of (a) high, (b) middle, and (c) low IR cloud cover during last solar cycle. The variations of galactic cosmic ray (GCR) fluxes as measured from CLIMAX (normalized to May, 1965) are also indicated in each panel (dot-dashed lines). The shaded areas in Figure 2(a) corresponding to the years that global high cloudiness might have been affected by volcano eruptions and El Niño event. The shaded areas in Figure 5(a) corresponding to the years that global high cloudiness might have been affected by volcano eruptions and El Niño event. From 1988 to 1993, the impact of volcano eruptions and El Niño on global high clouds is likely negligible, and it is during this period that we find a significant anti-correlation between GCR fluxes and high cloud anomalies. The increase of high cloudiness during 1988-1989 and the decrease of high cloudiness during 1991-1992 can be readily explained by the potential role of GCR in aerosol formation and CCN abundance. Furthermore, if we look at the whole period from 1984 to 1994, we do not see obvious enhancements of high cloudiness over average values due to volcano eruptions. If we consider what the volcano eruptions might have superimposed on the natural variable high cloudiness, it is clear that volcano eruptions might have enhanced the global high cloudiness by up to ~ 1.5% (without volcano effect, the high cloudiness during the solar minimum 1986-1987 is expected to be 1.5-2% less the value during solar maximum 1990-1992). In summary, the predicted different sensitivities of the particle production to cosmic ray changes at different altitudes seem to be consistent with the observed different correlations between cosmic ray variations and low, middle and high cloud anomalies. However, due to the limit of cloud cover data available and uncertainties in the volcano and El Niño impacts, our conclusions, especially with regard to the existence of anti-correlation between high cloudiness and cosmic ray variations, are not definitive. More research is obviously needed. 5. CLIMATE IMPLICATIONS While the first key step of GCR-CN-CCN-Cloud hypothesis seems to be consistent with the observed spatially-dependent correlations of GCR variations and global cloudiness, much more work is needed to clearly establish the GCR-Cloud connection. Nevertheless, it is meaningful to discuss the climate implications associated with the possible GCR-induced cloud changes. We assume that the anomalies of high cloud cover correlate negatively while that of low cloud cover correlate positively with GCR variations, and the magnitudes of the fluctuations are similar (1.5- 2% absolute change). As a result of opposite systematic variations of low and high clouds associated with solar activity, the total global cloud cover may show no obvious correlation with GCR variations. However, the radiative effects are unlikely to cancel each other. First, the net radiative forcing of clouds depends on their altitude and optical thickness. High optically thin clouds tend to warm while optically thick high and low clouds tend to cool2. Since cloud plays an import role in the Earth radiation budget, a systematic absolute increase of high cloud amount by ~1.5-2% and a decrease of low cloud amount by ~1.5-2% from solar minimum to solar maximum, if confirmed, may represent an important mechanism to amplify the effect of solar variability on Earth’s climate. Second, the opposite change in high and low clouds may change the atmosphere heating profile and the distribution of energy between the atmosphere and the surface, and hence may have far-reaching dynamical and climatic consequences. A systematic increase in high cloud may either warm or cool the atmosphere and the earth’s surface below, depending on the types of high clouds and the underlying atmospheric properties. For example, it has been shown that the presence of a cirrostratus (with a base height of 16 km and thickness of 1.5 km) in an otherwise clear tropical atmosphere has a net cooling effect for the atmosphere below ~6 km but has a net heating effect for the atmosphere above ~ 6 km when solar zenith angles are small (< ~ 60o)42. A systematic decrease in low cloud is likely to warm the surface by allowing more sunlight to reach the earth surface, while the same decrease will cool the lower troposphere by reducing the visible absorption in the cloud layer and infra-red absorption in the cloud layer and the atmosphere below42. The long-term trend of global low and high cloud cover as a result of GCR variations may become an important external forcing of Earth climate system. Based on observations, Lockwood et al43 have shown that from 1964 to 1996 the strength of the solar magnetic flux, shielding the Earth from GCR, has increased by ~ 41% while GCR has decreased by ~3.7%43. The ion chamber measurements44 also indicate that the sea level GCR intensity has decreased by ~2% from 1979 to 1994. The GCR intensity decrease is expected to be larger at higher altitudes in the troposphere. From the data available, we estimate that the decrease in GCR fluxes during the past two decades (1979-1999) is about 1/3 to 1/2 of the maximum variations during the last solar cycle. Thus, if the connection between low and high cloudiness exists, the global mean low cloud amount might have been decreasing (0.25-0.5% per decade) and high cloud amount increasing (0.25-0.5% per decade) during the past two decades. The net impact of GCR variations during the past two decades are likely to have warmed the earth surface but cooled the lower troposphere. Note that the potential GCR-induced change in cloud albedo and absorption may enhance such an impact (a decrease in cloud droplet concentration due to fewer CCN as a result of reduced GCR fluxes may also imply a low cloud albedo and absorption). While the exact amount of net radiative forcing associated with GCR-induced low and high cloud changes remains to be investigated, it is physically plausible that the decrease in GCR fluxes during the past two decades has led to a net warming of ~0.05 oC per decade at the surface while a net cooling of ~0.05 oC per decade in the lowest 8 km of atmosphere. In other words, the GCR-induced low and high cloud changes may explain why the Earth’s surface has warmed much more rapidly (~0.1 oC per decade) than the lowest 8 km of atmosphere during the last two decades. A piece of suggestive evidence to support such a claim is that the regions of large differences in surface and low atmospheric temperature trends remarkably coincide with the regions of high correlation between cosmic ray and low cloud top temperature as shown in Figures 6 and 7. The spatial correlation map in Figure 6 shows how ISCCP-D2 low cloud top temperatures co-vary with GCR flux. The correlation coefficients, r, are calculated from the 12-month running mean at each grid point (following Marsh and Svensmark 6,7). Figure 6 reveals a band of significantly high correlation centered around the tropics where stratocumulus and marine stratus clouds are dominant6,7. The fraction of Earth with r > 0.6 is ~ 30% which is significant at 99.9%. Figure 7 shows the global distribution of temperature trends in the lowest five miles of the atmosphere derived from MSU t2lt data for the period 1979-2000. If we compare MSU t2lt trends with similar temperature trends derived from surface observations (not shown), we see warming over the northern third of the globe both at the surface and in the five-mile-deep layer of air above. The largest differences in the trends of temperature between surface and low atmosphere measurements are over the tropical regions where the surface data show a significant warming while MSU t2lt data show a slightly cooling. The tropical radiosonde temperature data show the same patterns of surface warming and tropospheric cooling since 1979 as the independent surface and MSU observations16. As we have mentioned, the regions of large differences in surface and low atmospheric temperature trends remarkably coincide with the regions of high correlation between cosmic ray and low cloud top temperature (Figure 6). Such a nice coincidence may suggest that the differences in surface and low atmospheric temperature trends over tropics are associated with the solar indirect forcing via GCR-Cloud link. It may become necessary to include the solar indirect forcing via GCR-induced cloud change in the future climate model, as current models still cannot account fully for the apparent difference between observed surface and troposphere temperature trends since 1979 13-15. Unlike homogeneous greenhouse gases which warm both the surface and low troposphere, the potential influence of GCR variations on clouds are different at different regions of the atmosphere and the associated radiative forcing are spatial and temporal inhomogeneous. The observed rapid warming during the past two decades over the northern third of the globe both at the surface and in the air above is likely due to the greenhouse effect. Figure 6. Global correlation map of GCR with anomalies of IR low cloud top temperature. White pixels indicate regions with either no data or an incomplete monthly time series. Figure 7. Low-to-middle atmosphere decadal linear temperature trend ( oC/decade) for period 1979-2000 derived from MSU t2lt data. 6. SUMMARY AND DISCUSSION The dependence of ultrafine production rate on galactic cosmic ray ionization rate at different altitudes has been investigated. Our primary studies indicate that an increase in GCR ionization rate leads to an increase in CN production in the lower troposphere (>680 mb), but a decrease in CN production in the upper troposphere (<440 mb). In the lower troposphere, the ionization rate is low and the H2SO4 concentration is relatively high, the particle formation is limited by ionization and an increase in ionization rate leads to an increase in nucleation. In the upper troposphere, the ionization rate is very high and the H2SO4 concentration is relatively low, the particle formation is limited by H2SO4 concentration and an increase in ionization rate inhibit the nucleation by reducing the lifetime of ion clusters. The average change of CN production as the ionization rate increases is small in the middle troposphere (440-680 mb). Since an increase in ultrafine production rate is likely to increase the CCN abundance and cloudiness, we can expect that the correlation between GCR changes and global cloud cover (if any) should be positive for low cloud, negative for high cloud, and weak for the middle cloud. In addition to the reported positive correlation between GCR variations and low cloudiness, our analyses of ISCCP D2 IR cloud data further reveal that high cloudiness may be anti-correlated with GCR variations if volcano and El Niño impacts are excluded. During a solar cycle, the absolute change of high and low cloud amounts is opposite in sign but similar in magnitude (~1.5-2%). The fluctuations of middle cloud anomalies are small compared to that of low clouds, and no obvious correlation exists between middle cloudiness and GCR variations. Therefore, the observed different correlations between GCR variations and low, middle and high cloud anomalies seem to be consistent with the predicted dependence of CN production on GCR variations at different altitudes. Such a consistency suggests that solar activity might affect global cloudiness by modulating GCR fluxes. However, due to the limit of cloud cover data available and uncertainties in the volcano and El Niño impacts, our conclusions, especially with regard to the existence of anti-correlation between high cloudiness and cosmic ray variations, are not definitive. The climate implications associated with the possible GCR-induced cloud changes are discussed. Since cloud is critical to Earth radiation budget, opposite systematic variations of low and high clouds associated with solar activity, if confirmed, may represent an important mechanism to amplify the effect of solar variability on Earth’s climate. The decrease in GCR intensity during the last two decades might have led to a decrease in global mean low cloud amount and an increase in high cloud amount, which might have warmed the Earth’s surface and cooled the low troposphere. The potential GCR-induced change in cloud albedo and absorption may enhance such an impact (a decrease in cloud droplet concentration due to fewer CCN as a result of reduced GCR fluxes may also imply a low cloud albedo and absorption). We suggest that, the GCR-induced natural variability of global cloudiness, together with the greenhouse gases which warm both the surface and low troposphere, may reconcile the apparent differences in global mean temperature trends at Earth’s surface (rapidly warming, as recorded by thermometers) and in the lowest 8 km of atmosphere (little if any warming, as monitored by satellites and balloons). While this study provides additional evidence for the inferred correlation between variations in global cloud properties and the solar-modulated GCR fluxes, much more work is needed to understand how and how much the GCR variations will affect the global cloud properties. The first key process (i.e, influence of GCR variations on nucleation and CN abundance) in our proposed GCR-CN-CCN-Cloud hypothesis seems to be consistent the spatially dependent influence of GCR variations on cloud properties. However, we currently do not know how much the natural GCR variations will affect the CCN abundance and cloud properties. Laboratory and field measurements, as well as theoretical studies are needed to validate the predicted dependent- behaviors of nucleation on ionization rates at different altitudes, to investigate the effect of GCR variations on CCN abundance, and to clarify the complex microphysics of aerosol/cloud interactions. The current analyses of GCR-cloud correlations are limited by the uncertainties associated with the cloud data and short periods of cloud data available. Improved cloud cover data covering longer time periods will be very useful in studying GCR-cloud connections. ACKNOWLEDGEMENTS This work was supported by NSF grant and start-up fund from State University of New York, Albany. REFERENCES 1. IPCC, Climate Change 1995: The Science of Climate Change, Intergovernmental Panel on Climate Change. Ed. by J. T. Houghton et al., Cambridge Univ. Press, New York (1996). 2. Hartmann, D.L. Radiative effects of clouds on Earth’s climate, in Aerosol-Cloud-Climate Interactions. Ed. by P.V. Hobbs, Academic Press Inc., San Diego (1993). 3. Svensmark, H., and E. Friis-Christensen, Variation of cosmic ray flux and global cloud coverageÐA missing link in solar-climate relationships. J. Atmos. Sol. Terr. Phys., 59, 1225 (1997). 4. Rossow, W.B., and R.A. Schiffer, ISCCP Cloud Data Products. Bull. Am. Met. Soc., 72, 2- 20 (1991). 5. Rossow, W.B., and R.A. Schiffer, Advances in Understanding Clouds From ISCCP. Bull. Amer. Meteor. Soc., 80, 2261-2287 (1999). 6. Marsh, N., and H. Svensmark, Cosmic rays, clouds, and climate. Space Sci. Rev., 94, 215- 230 (2000a). 7. Marsh, N., and H. Svensmark, Low Cloud Properties influenced by Cosmic Rays, Phys. Rev. Lett., 85, 5004-5007 (2000b). 8. Kernthaler, S. C., R. Toumi, and J. D. Haigh, Some doubts concerning a link between cosmic ray fluxes and global cloudiness. Geophys. Res. Lett., 26, 863-865 (1999). 9. J¿rgensen, T. S., and A. W. Hansen, Comments on “Variation of cosmic ray flux and global cloud coverage Ð a missing link in solar-climate relationships” by Henrik Svensmark and Eigil Friis-Christensen. J. Atmos. Sol. Terr. Phys., 62, 73-77 (2000). 10. Kristjánsson, J. E., and J. Kristiansen, Is there a cosmic ray signal in recent variations in global cloudiness and cloud radiative forcing? J. Geophys. Res., 105, 11,851-11863 (2000). 11. IPCC, Climate Change 2001: The Scientific Basis, Working Group I contribution to the IPCC Third Assessment Report: Summary for policymakers. Intergovernmental Panel on Climate Change, Shanghai Draft 21-01-2001 (2001). 12. NRC, Reconciling Observations of Global Temperature Change. Panel on Reconciling Temperature Observations, National Research Council, National Academy Press, Washinton D. C. (2000). 13. Hansen, J. M., M. Sato, and R. Ruedy, Radiative forcing and climate response. J. Geophys. Res., 102, 6831-6864 (1997). 14. Bengtsson, L., Roeckner, E., and Stendel, M., Why is the global warming proceeding much slower than expected? J. Geophys. Res., 104, 3865-76 (1999). 15. Santer, B. D., et. al., Interpreting differential temperature trends at the surface and in the lower troposphere. Science, 287, 1227-1232 (2000). 16. Gaffen, D. J., et al., Multidecadal changes in the vertical temperature structure of the tropical troposphere. Science, 287, 1242-1245 (2000). 17. Millikan, R. A., H. V. Neher, and W. H. Pickering, Further studies on the origin of cosmic rays. Phys. Rev., 66, 295-302 (1944). 18. Neher, H.V., Cosmic-Rays at high latitudes and altitudes covering four solar maxima. J. Geophys. Res., 76, 1637-1651 (1971). 19. Yu, F., and R. P. Turco, Ultrafine aerosol formation via ion-mediated nucleation. Geophys. Res. Lett., 27, 883-886 (2000). 20. Yu, F. and R. P. Turco, From molecular clusters to nanoparticles: The role of ambient ionization in tropospheric aerosol formation. J. Geophys. Res., 106, 4797-4814 (2001). 21. Weber, R. J., et al., Measurement of new particle formation and ultrafine particle growth rates at a clean continental site. J. Geophys. Res., 102, 4375-4385 (1997). 22. Yu, F., and R. P. Turco, The role of ions in the formation and evolution of particles in aircraft plumes. Geophys. Res. Lett., 24, 1927-1930 (1997). 23. Yu, F., R. P. Turco and B. Kärcher, The possible role of organics in the formation and evolution of ultrafine aircraft particles. J. Geophys. Res., 104, 4079-4087 (1999). 24. Kärcher, B., F. Yu, F. P. Schroeder and R. P. Turco, Ultrafine aerosol particles in aircraft plumes: Analysis of growth mechanisms. Geophys. Res. Lett., 25, 2793-2796 (1998). 25. Weber, R. J., et al., New particle formation in the remote troposphere: A comparison of observations at various sites. Geophys. Res. Lett., 26, 307-310 (1999). 26. Clarke, A. D., Kapustin, V. N., Eisele, F. L., Weber, R. J., and P. H. McMurry, Particle production near marine clouds: Sulfuric acid and predictions from classical binary nucleation. Geophys. Res. Lett., 26, 2425-2428 (1999). 27. Twomey, S., The influence of pollution on the shortwave albedo of clouds. J. Atmos. Sci., 34, 1149-1152 (1977). 28. Albrecht, B.A., Aerosols, cloud microphysics, and fractional cloudiness. Science, 245, 1227-1230 (1989). 29. Hobbs, P.V., Aerosol-cloud interactions. In: Aerosol-cloud-climate interactions. Ed. by P.V. Hobbs, Academic Press Inc., San Diego (1993). 30. Stubenrauch, C. J., W. B. Rossow, F. Chéruy, A. Chédin, and N.A. Scott, Clouds as seen by satellite sounders (3I) and imagers (ISCCP), Part I, Evaluation of cloud parameters. J. Clim., 12, 2189-2213 (1999). 31. Haeffelin, M., R. Kandel, and C. Stubenrauch, Improved diurnal interpolation of reflected broadband shortwave observations using ISCCP data. J. Atmos. Oceanic Technol., 16, 38- 54 (1999). 32. Gates, W. L., et al., An overview of the results of the Atmospheric Model Intercomparison Project (AMIP I), Bull. Am. Meteorol. Soc., 80, 29-55 (1999). 33. Jensen, E.J., and O. B. Toon, The potential effects of volcanic aerosols on cirrus cloud microphysics. Geophys. Res. Lett., 19, 1759-1762 (1992). 34. Knollenberg, R.G., Kelly, K., and J. C. Wilson, Measurements of high number densities of ice crystals in the tops of tropical cumulonimbus. J. Geophys. Res., 98, 8639-8664 (1993). 35. Wylie, D. P., W. P. Menzel, H. M. Woolf, and K. I. Strabala, Four years of global cirrus cloud statistics using HIRS. J. Clim., 7, 1972-1986 (1994). 36. Trepte, C. R., and M. H. Hitchman, Tropical stratospheric circulation deduced from satellite aerosol data. Nature, 355, 626Ð628 (1992). 37. Trepte, C. R., R. E. Veiga, and M. P. McCormick, The pole-ward dispersal of Mount Pinatubo volcanic aerosol. J. Geo-phys. Res., 98, 18,563Ð18,573 (1993). 38. Holton, J. R., P. H. Haynes, M. E. McIntyre, A. R. Douglass, R. B. Rood, and L. Pfister, Stratosphere-troposphere ex-change. Rev. Geophys., 33, 403Ð439 (1995). 39. Kent, G. S., Wang, P. H., McCormick, M.P., and K. M. Skeens, Multiyear stratospheric aerosol and gas experiment II measurements of upper tropospheric aerosol characteristics. J. Geophys. Res., 100, 13,875-13,899 (1995). 40. Kent, G. S., Trepte, C. R., and P. L. Lucker, Long-term stratospheric aerosol and gas experiment I and II measurements of upper tropospheric aerosol extinction. J. Geophys. Res., 103, 28,863-28,874 (1998). 41. Wang et al., A 6-year climatology of cloud occurrence frequency from Stratospheric Aerosol and Gas Experiment II observations (1985-1990). J. Geo-phys. Res., 101, 29407- 29 (1996). 42. Liou, K.N., Radiation and Cloud Processes in the Atmosphere: Theory, Observation, and Modeling. p305, 393-395, Oxford University Press, New York, 487 pp (1992). 43. Lockwood, M., R. Stamper, and M.N. Wild, A doubling of the Sun’s coronal magnetic field during the past 100 years. Nature, 399, 437-439 (1999). 44. Ahluwalia, H.S., Galactic cosmic ray intensity variations at a high latitude sea level site 1937-1994. J. Geophys. Res., 102, 24229-24236 (1997). CLOUD DROPLET GROWTH R. Bingham Rutherford Appleton Laboratory, Chilton, Didcot,Oxon, U.K. Abstract We present a new mechanism for rapid cloud droplet growth. The proposed mechanism relies on an isotropic pressure introduced as a result of shadowing between two droplets and can occur in both a charged and neutral atmosphere. We consider the possibility of enhancing the growth of cloud droplets in the presence of charged particles in particular ions produced from cosmic rays. 1. INTRODUCTION The subject of raindrop formation in the presence of charged particles is not new, in C.T.R. Wilson 1899 reported on experiments promoting the formation of raindrops. Wilson 1899 pointed out that a “slight rain-like condensation takes place in which a supply of ions has been produced by the action of Ro¨ntgen rays or other ionising agent”. The model we propose for the rapid growth in rain droplet size was originally used to explain the formation of dust particulates in the plasma etching process (Bingham & Tsytovich, 2001). In these plasma experiments dust agglomeration has been shown to be important in laboratory etching experiments (Garscadden et al. 1994) where the growth of dust is extremely rapid and is due to dust-dust attraction by a plasma or neutral bombardment force known as the shadow force (Bingham & Tsytovich, 2001). In the atmosphere charged drops are the norm and exist in an atmosphere of charged ions and electrons whose densities are enhanced possibly by cosmic rays. The growth of the droplet is then influenced by the presence of the plasma particles. The shadow force is caused by both ion and neutral atom droplet collisions. The mechanism of the shadow attraction even for droplets of the same charge is described by a relatively simple expression. In a charged atmo- sphere Debye screening significantly reduces the repulsive Coulomb force while the attractive force due to the bombardment are not affected by screening and dominate. The nature of the attractive force is due to bombardment by charge plasma particles and neutral atoms. For a single droplet in a plasma atmosphere with an isotropic distribution of particles, direct bombardment and deflection transfers no net momentum to an isolated droplet and therefore no net force acts on a single drop. Another drop at distance r shadows the flux to the first with a solid angle a2/r2 where a is the drop radius. the net momentum transfer is proportional to the solid angle, to the surface 2 area of the drop a and to the neutral and plasma pressure nkBT where n is the density of neutrals or plasma, kB is Boltzmann’s constant and T is the temperature of neutral/plasma. The force imported by the anisotropic pressure is given by a2 F + nk T a2 (1) r2 B a and attractive potential Ua given by a2 U = nk T a2 (2) a a r B where the coefficient a consists of three parts = 11.1 ± 2.1 fSm-1, compared to day 170 when -1 = 14.8 ± 8.6 fSm . There was also a weak positive correlation (significant to > 90%) with Sg, suggesting that there was a different régime to the solar-induced turbulent mixing discussed above. Mechanical turbulence may have dominated, particularly in the morning when the wind was approximately constant at 3.5 ms-1, due to the remnants of a cold front. On day 170 the wind arises primarily from local convective mixing and shows the usual diurnal cycle with a maximum in the afternoon. During the afternoon, which was sunnier on day 165, both the s and wind traces for days 165 and 170 are more closely related, with the conductivity minimum occurring at similar times, associated with the maximum temperature. 2.4 Factors influencing conductivity for two cloud-free periods Here the predominantly clear period during day 169 is compared to the same time for day 170. Direct comparisons of two days are not normally possible because of the great variation in meteorological conditions. Days 169 and 170 were therefore rare because the synoptic conditions and incoming solar radiation were very similar, allowing detailed analysis of specific effects on s. Day 169 was almost cloud-free from 0530 until sunset. Some cirrus cloud was recorded manually at 1000, radar data† however showed no cloud over the whole day. There appears to have been – st. dev. (fSm-1) (ms ) RH(%) Dirn (deg)