Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.izmiran.rssi.ru/~obridko/papers/dorman.pdf Дата изменения: Thu Feb 10 16:13:56 2005 Дата индексирования: Sun Apr 10 01:05:01 2016 Кодировка: Поисковые слова: quasar

PREDICTION OF EXPECTED GLOBAL CLIMATE CHANGE BY FORECASTING OF GALACTIC COSMIC RAY INTENSITY TIME VARIATION IN NEAR FUTURE BASED ON SOLAR MAGNETIC FIELD DATA
A.V. Belov 1, L.I. Dorman1,2, R.T. Gushchina1, V.N. Obridko1, B.D. Shelting1, and V.G. Yanke
1

(1) IZMIRAN of Russian Academy of Science, RUSSIA; (2) Israel Cosmic Ray & Space Weather Center and Emilio Segre' Observatory affiliated to Tel Aviv University, Technion and Israel Space Agency, ISRAEL;E-mail: lid@physics.technion.ac.il, Fax: 972-4-6964952 ABSTRACT A method of prediction of expected part of global climate change caused by cosmic ray (CR) by forecasting of galactic cosmic ray intensity time variation in near future based on solar activity data prediction and determined parameters of convection-diffusion and drift mechanisms is presented. This gave possibility to made prediction of expected part of global climate change, caused by long-term cosmic ray intensity variation. In this paper we use the model of cosmic ray modulation in the heliosphere, which considers a relation between long-term cosmic ray variations with parameters of the solar magnetic field. The later now can be predicted with good accuracy. By using this prediction, the expected cosmic ray variations in the near Earth space also can be estimated with a good accuracy. It is shown that there are two possibilities: 1) to predict cosmic ray intensity for 1-6 months by using a delay of long-term cosmic ray variations relatively to effects of the solar activity and 2) to predict cosmic ray intensity for the next solar cycle. For the second case the prediction of the global solar magnetic field characteristics is crucial. For both cases reliable long-term cosmic ray and solar activity data as well as solar magnetic field are necessary. For solar magnetic field we used results of two magnetographs (from Stanford and Kitt Peak Observatories). The obtained by described method prediction on long-term cosmic ray intensity variation we used for estimation of expected part of global climate change caused by cosmic rays. SOLAR-HELIOSPHERIC PARAMETERS FOR MULTI-PARAMETRIC MODEL OF COSMIC RAY MODULATION

Long-term variations of galactic cosmic rays were compared with the behavior of various solar-activity indices and heliospheric parameters many times. A close relation of this parameter to the behavior of cosmic rays in the past decades was substantiated theoretically by (see e.g. Jokipii and Thomas, 1981; Kota and Jokipii, 1983) and proved experimentally by (see e.g. Simpson, 1963, Dorman and Dorman, 1967, Smith and Thomas, 1986; Webber and Lockwood, 1988; Smith 1990; Belov et al., 1995, 1999a,b, 2003; Mavromichalaki et al., 1995, 1998, Bazilevskaya and Svirzhevskaya, 1998, McCracken and McDonald, 2001, Usoskin et al., 2001). In this paper a model of long-term modulation of cosmic rays has been constructed using a set of parameters of the global solar magnetic field. A reliable model of CR modulation should incorporate at least several solar-heliospheric parameters. The choice of the parameters is explained in detail by Belov et al. (2002). We obtain the best results by using the heliospheric current sheet tilt , the intensity of the magnetic field radial component Br averaged over the entire source surface - BSS of the source surface magnetic field (SSMF), the polarity p (qA>0, qA<0) of the general solar magnetic field.


Expected CR variations were calculated basing on characteristics of the solar wind source surface field determined from observations on the Stanford and Kitt Peak Observatories. An existing difference between characteristics and observational results of the two magnetographs didn't allow choosing the best one among them for prediction purpose. These magnetographs register fields of different spatial scale. Mean intensity Bss of the source surface magnetic field and HCS tilt are successfully complement to each other. By the polarity of the large-scale solar magnetic field p, we mean the polarity of its dipole component at the north pole of the Sun.
We suppose that CR intensity long-term variation may be determined as following (taking into account direct action of polarity effect on CR as well as actions with different times-lag for other parameters): bp up b u bB uB =a+ (1) Е (1 +b p p f (t - )) (t - ) + + 1 Е Bss (t - ) + + 1 Е ps (t - ) , u + 1 =0 =0 =0 uB up where ps (t - ) and p f (t - ) are the periods of inversion of the solar magnetic field were determined on the source surface from optical observations and magnetic observations. We have applied this model to describe CR variations (10 GV rigidity) during a rather long period, when reliable information on solar magnetic field is available. The calculations were performed for parameters: a = 10.9±0.6, b= - 0.205±0.001 %/є, bp = - 0.65±0.08, bB = - 1.42±0.09 %/nT, bp = - 3.2±0.2 %. Were determined also u= 7 months, uB = 10 months, up= 3 months. The obtained correlation coefficient is 0.94. Results are shown in Fig. 1.

Fig. 1. Observed and simulated CR variations for 1976-2003 years (lower curves). Contributions to simulated CR variation from of mean source surface magnetic field intensity BSS, the heliospheric current sheet tilt , heliomagnetic polarity changes p are shown as upper curves.


The analysis of results shown in Fig. 2, leads to following conclusions. Two characteristics of the source surface magnetic field ­ the structural (HCS tilt angle) and quantitative (mean field BSS) ones ­ well supplement each other in describing CR variations. The changing of HCS tilt angle controls long-term variations (11-year cycles and their basic features), while BSS is responsible mainly for shorter period variations. Correspondingly, the HCS tilt angle plays the leading part in the periods of low and moderate solar activity, yielding to BSS in the vicinity of the cycle maxima. Results of modeling have revealed differences of CR modulation during periods with different direction of the solar global magnetic field that is described by influence of cyclic changes of local and large scale fields of the Sun registered by the observatories mentioned above on long-term CR variations. Both the Stanford and Kitt Peak observations could be used for CR forecasting Calculations performed by Obridko and Shelting (2001) show that the magnetic field at the source surface changes its sign much earlier than in the photosphere. The time boundaries in the photosphere and at the source surface were determined using both the line-of-sight observations of the polar field and the field calculated in radial direction. The behavior of CR correlated with the polarity changes at the source surface much better than in the photosphere. The CR variation is most closely related to the sign of magnetic fields obtained from H observations for the solar wind source surface- pHBss. The obtained boundaries are as follows: 09.1979­03.1981, 10.1989­ 03.1991, and 04.1999-06.2001. THE COSMIC RAY FORECASTING MODEL As the source of the CR variation data, we are using the Moscow neutron monitor. Since 1997, these data have been available through Internet in real time. The delay of CR variations reflects complicated individual features of the activity cycles and is related to the solar magnetic cycle. We suggest using this particularity of CR variations to predict the CR flux in the near-Earth space. In Fig. 2 are shown results of CR intensity forecasting bases on the method described in the previous Section for several months ahead (up to 12 months). The results (the were obtained for three periods: 1977-2003 ­ all time with 3 reversals of the global solar magnetic field; 1981­1989 ­ negative polarity of the global solar magnetic field (qA<0); 1991­2000 ­ positive polarity of the global solar magnetic field period (qA>0).

correlation coefficient and standard r.m.s. deviation as a function of time )

Fig. 2. The correlation coefficient ( three low curves, the left ordinate) and the standard r.m.s. deviation (upper three curves, the right ordinate) between the observed CR variations (Moscow neutron monitor) and the predicted CR variations with , Bss and p as modulating parameters for different time of the forecast (in months).


From Fig. 3 can be seen that the analysis of data for 1977­2003 has revealed a good correlation ( = 0.82 В 0.94) between the simulated CR variations with , Bss and p as modulating parameters of the multi-parameter model and galactic CR intensity behavior during long period, spanning several cycles of solar activity. The model gives the good prediction for 1-5 months in advance. Beginning from ~ 6 months forecasting time the quality of forecasting is not so good. Fig. 2 shows also that the using data for the period with one polarity (qA < 0 or qA > 0) gives much better results (higher correlation coefficient), and for longer number of months ahead (up to about 10 months). FORECAST OF COSMIC RAY INTENSITY MODULATION FOR 1 AND 12 MONTHS In Fig. 3 are shown results of CR intensity forecasting for 1 month ahead (the upper panel) in comparison with observed CR intensity variation. In the bottom panel are shown the same but for CR intensity forecasting for 12 months ahead.

Fig. 3. Observed on Moscow neutron monitor and predicted for 1 and 12 months ahead - the long-term CR variations.

We must admit that the model applied herein meets difficulties when forecasting for half a year or more. The longer the term of the forecast, the larger the fraction of the inner heliosphere excluded from consideration. More promising seems the approach that involves prediction of the main parameters of the solar magnetic field. Then, the forecast of CR variations would be based both on the measured solar indices and on their expected future behavior.
THE PREDICTION OF COSMIC RAY INTENSITY VARIATION FOR THE NEXT SOLAR CYCLE Predictions of global solar magnetic field characteristics are crucial for predictions long-term CR variations. CR variations for the next solar cycle are predicted by solar magnetic field data of two magnetographs (Stanford and Kitt Peak) using the statistical method. According to our predictions of CR behavior till 2013 the next CR maximum should be between the end of 2006 and the beginning of 2007 by observations in the Kitt Peak observatory, but it should be in the middle of 2006 by data of the Stanford observatory (see Fig. 4).


Fig. 4. The forecast of the CR behavior based on the predicted values of the global characteristics of the solar magnetic field, thick line ­ data of CR intensity observations (Moscow NM), thin line ­ the predicted CR variation up to 2013 based on data of Kitt Peak Observatory (upper panel) and based on data of Stanford Observatory (bottom panel). ON THE CONNECTION OF CR SOLAR CYCLE VARIATION WITH VARIATION OF PLANETARY CLOUD COVERAGE A very important result for an understanding of the mechanism of the influence of solar activity cycle on the Earth's climate has recently been obtained: it was found that the Earth's cloud coverage (observed by satellites) is strongly correlated with CR intensity (Svensmark and Friis-Christensen, 1997; Svensmark, 1998, 2000; Marsh and Svensmark, 2000a,b). Clouds influence irradiative properties of the atmosphere by both cooling through reflection of incoming short wave solar radiation, and heating through trapping of outgoing long wave radiation. The total result depends mostly on the height of the clouds. According to Hartmann (1993), high optically thin clouds tend to heat while low optically thick clouds tend to cool. It was shown that low clouds give a cooling of about 17 W.m-2 , so they play an important role in the Earth's radiation budget (Ohring and Clapp, 1980; Ramanathan et al., 1989; Ardanuy et al., 1991). So even small changes in the lower cloud coverage can give important changes in the radiation budget and considerably influence the Earth's climate (let us remember that the solar irradiance changes during solar cycle by only about 0.3 W.m-2 ). Fig. 5 shows the composite of satellite observations of the Earth's total cloud coverage in comparison with CR intensity (according to Climax NM) and solar activity data (intensity of 10.7 cm solar radio flux). From Fig.5 it can be seen that the correlation of global cloud coverage with CR intensity is much better than with solar activity. Marsh and Svensmark (2000a) came to conclusion that CR intensity connects very well with low global cloud coverage, but not with high and middle clouds. It is important to note that low clouds lead, as rule, to the cooling of the atmosphere. It means that with increasing CR intensity and cloud coverage (see Fig. 5), the surface temperature is expected to decrease. It is in good agreement with the situation for the last 1000 years, and with direct measurements of the surface temperature for the last four solar cycles.


Fig. 5. Changes in the Earth's cloud coverage: triangles - from satellite Nimbus 7, CMATRIX project, (Stowe et al., 1988); squares - from the International Satellite Cloud Climatology Project, ISCCP, (Rossow and Shiffer, 1991); diamonds ­ from the Defense Meteorological Satellite Program, DMSP (Weng and Grody, 1994, Ferraro et al., 1996). Solid curve ­ CR intensity variation according to Climax NM, normalized to May 1965. Broken curve ­ solar radio flux at 10.7 cm (in units 10 -22 W.m -2 .Hz -1 ). All data are smoothed using 12 months running mean. According to Svensmark (2000). EXPECTED PART OF CLIMATE CHANGE CAUSED BY CR INTENSITY VARIATION From Fig. 5 can be seen that about 20% of CR intensity increase in Climax NM for solar cycle corresponds to about 4% increase of global cloud covering, what can give sufficient change in radiation balance influenced on climate change. From other side, from Fig.4 can be seen that the accounting of data on general solar magnetic field give a possibility to predict with a good accuracy expected CR intensity variation. For example, for the period of few years (up to about 2006) the CR intensity expected to be increase on about 10%, so it is expected some global climate cooling and increasing of precipitation corresponded to increase of the global cloud covering on about 22.5%. Of course, this cooling can be compensated with the process of global warming caused by increasing of green gases, but in any case it is necessary to take into account all processes influenced on global climate change. CONCLUSIONS Our conclusions are as follows: 1. On the basis of solar magnetic field data by taking into account tilt angle data and data on reversals periods, can be made very good determination of CR intensity change and made prediction for the future with correlation coefficient between observed and predicted intensities about 0.97 for 1 month prediction, 0.91 for 6 months prediction, and 0.80 for 12 months prediction (see Fig. 2 and 3). 2. For the period of about 2 years ahead (up to 2006) the CR intensity expected to be increase on about 10%, so it is expected some global climate cooling and increasing of precipitation corresponded to increase of the global cloud covering on about 2.0-2.5% during these two years.


ACKNOWLEDGEMENTS This research is partly supported by EU Project COST 724, RFBR grant 02-02-16992 and 02-02-16199. We thank all people responded for the internet accessible data on the solar data: http://sun.stanford.edu/~wso/wso.html, http://www.nso.noao.edu/nsokp/dataarch.html.

REFERENCES Ardanuy P.E., L.L. Stowe, A. Gruber, and M. Weiss, Shortwave, longwave and net cloud-radiative forcing as determined from NIMBUS-7 observations, J. Geophys. Res., 96, 18537-18549, 1991. Bazilevskaya, G.A. and Svirzhevskaya, A.., On the Stratospheric Measurements of the Cosmic Rays, Space Sci. Rev. 85, 431-521, 1998. Belov A.V., Gushchina R.T., and Sirotina I.V., Long term cosmic ray variations and thier relation with solar activity parameters, Proc. 24-th ICRC, 4, 542-546. 1995. Belov, A.V., Gushchina, R.T., Yanke, V.G. et. al., Relation of long-term variations of cosmic rays to the magnetic field in the Sun and solar wind, Izv. RAN, ser.phys., 63, N8, 1606-1610, 1999a. Belov, A.V., Gushchina, R.T., & Yanke, V.G., On connection of cosmic ray long term variations with solar ­ heliospheric parameters, Proc. 26-th ICRC, 7, 175-179, 1999b. Belov, A.V., Shelting, B.D., Gushchina, R.T., Obridko, V.N., Yanke, V.G., The relation of the long-term cosmic rays modulation with the solar magnetic field characteristics, Geomag. and aeronomy, 42 (6), 727-735, 2002. Belov A.V., R.T. Gushchina, V.N.,Obridko, B.D.Shelting, V.G.Yanke, On the possibility to forecast longterm variations of galactic cosmic rays, Izv. RAN, ser.phys., 67 (4), .508-510, 2003.

Dorman L.I., Estimation of long-term cosmic ray intensity variation in near future and prediction of expected global climate change, Paper D2.1/C2.2/E3.1-0098-04, 2004.
Dorman, I.V., Dorman, L.I., On the nature of the lag of the cosmic ray intensity changes with respect to the solar activity changes, Cosmic rays. Articles, N8, 100-109, 1967. Ferraro R.R., F. Weng, N.C. Grody, and A. Basist, An eight-year (1987-1994) time series of rainfall, clouds, water vapor, snow cover, and sea ice derived from SSM/I measurements, Bull. Am. Met. Soc., 77, 891-905, 1996. Hartmann D.L., Radiative effects of clouds on the Earth's climate in aerosol-cloud-climate interactions, in P.V. Hobbs (ed) Aerosol-Cloud-Climate Interactions, Academic Press, 151, 1993. Jokipii, J.R., Thomas, B.T., Effects of drift on the transport of cosmic rays. Modulation by a wavy interplanetary current sheet, Astrophys. J. 243, 1115, 1981. Kota, J. and Jokipii, J.R., Effects of Drift on the Transport of Cosmic Rays. A Three-Dimensional Model Including Diffusion, Astrophys. J. 265, 573-581, 1983. Marsh N.D. and H. Svensmark , Low cloud properties influenced by cosmic rays, Phys. Rev. Lett., 85, 5004-5007 2000a. Marsh N. and H. Svensmark, Cosmic rays, clouds, and climate, Space Sci. Rev., 94, No. 1-2, 215-230, 2000b. Mavromichalaki H. Belehaki A. Rafios X. Tsagouri I.,Hale-Cycle Effects in Cosmic-Ray Intensity During the Last Four Cycles, Astrophysics and Space Science, 246, N1, 7-14, 1996. Mavromichalaki H., A. Belehaki and X. Rafios, Simulated effects at neutron monitor energies: evidence for a 22-year cosmic-ray variation, Astronomy and Astrophysics, 330, 764, 1998. McCracken K. G., McDonald F.B., The long term modulation of the galactic cosmic radiation, 1500-2000, Proc. 27-th ICRC, 9, 3753-3756, 2001. Obridko, V.N. and Shelting, B.D., Polar Global Magnetic Field of the Sun Reversals, Proc. of the Symposium "The Sun at the era of the magnetic field reversal", May 28 - June 2001, St ­ Petersburg, 2001. Ohring G. and P.F. Clapp, The effect of changes in cloud amount on the net radiation on the top of the atmosphere, J. Atmospheric Sci., 37, 447-454, 1980.


Ramanathan V., R.D. Cess, E.F. Harrison, P. Minnis, B.R. Barkstrom, E. Ahmad, and D. Hartmann, Cloudradiative forcing and climate: results from the Earth Radiation Budget Experiment, Science, 243, No. 4887, 5763, 1989. Rossow W.B. and R. Shiffer, ISCCP cloud data products, Bull. Am. Met. Soc., 72, 2-20. 1991. Simpson, J.A., Recent investigation of the low energy cosmic and solar particle radiations. Preprint N25, Vatican, 323-352, 1963. Smith, E.J., Thomas, B.T., Latitudinal Extent of the Heliospheric Current Sheet and Modulation of Galactic Cosmic Rays, J. Geophys. Res., 91, A3, 2933-2942, 1986. Stowe L.L., C.G. Wellemayer, T.F. Eck, H.Y.M. Yeh, and the Nimbus-7 Team, Nimbus-7 global cloud climatology, Part 1: Algorithms and validation, J. Climate, 1, 445-470. 1988. Svensmark H.,Influence of cosmic rays on the Earth's climate, Phys. Rev. Lett., 81, 5027-5030. 1998. Svensmark H., Cosmic rays and Earth's climate, Space Sci. Rev., 93, 175-185, 2000. Svensmark H. and E. Friis-Christiansen, Variation of cosmic ray flux and global cloud coverage ­ a missing link in solar-climate relationships, J. Atmos. Solar-Terr. Phys., 59, 1225-1232, 1997. Usoskin I.G., Mursula K., and Kovaltsov G.A., Odd and even cycles in cosmic rays and solar activity. Proc. 27-th ICRC, 9, 3791-3794, 2001. Webber, W.R., Lockwood, J.A, Charachteristics of the 22-year Modulation of Cosmic Rays as Seen by Neutron Monitors, J. Geophys. Res. 93, 8735-8740, 1988.

Weng F. and N.C. Grody, Retrieval of cloud liquid water using the special sensor microwave imager (SSM/I), J. Geophys. Res., 99, 25535-25551, 1994.
FIGURE CAPTIONS Fig. 1. Observed and simulated CR variations for 1976-2003 years (lower curves). Contributions in simulated CR variation from SSMF mean intensity Bss, HCS tilt angle, p helio-magnetic polarity changes are shown as upper curves. Fig. 2. The correlation coefficient ( three low curves, the left ordinate) and the standard r.m.s. deviation (upper three curves, the right ordinate) between the observed CR variations (Moscow neutron monitor) and the predicted CR variations with , Bss and p as modulating parameters for different time of the forecast (in months). Fig. 3. Observed (red) on Moscow neutron monitor and predicted (blue) for 1 and 12 months ahead - the long-term CR variations. Fig. 4. The forecast of the CR behavior based on the predicted values of the global characteristics of the solar magnetic field, red line ­ data of CR intensity observations (Moscow NM), black line ­ the predicted CR variation up to 2013 based on data of Kitt Peak Observatory (upper panel) and based on data of Stanford Observatory (bottom panel). Fig. 5. Changes in the Earth's cloud coverage: triangles - from satellite Nimbus 7, CMATRIX project, (Stowe et al., 1988); squares - from the International Satellite Cloud Climatology Project, ISCCP, (Rossow and Shiffer, 1991); diamonds ­ from the Defense Meteorological Satellite Program, DMSP (Weng and Grody, 1994, Ferraro et al., 1996). Solid curve ­ CR intensity variation according to Climax NM, normalized to May 1965. Broken curve ­ solar radio flux at 10.7 cm (in units 10 -22 W.m -2 .Hz -1 ). All data are smoothed using 12 months running mean. According to Svensmark (2000).