>',iA3A-CR-191955

Hans-F. Graf

Ingo Kirchner

Ingrid Schult

Max- PIe nck-In stitut

fur Meteorologie

/ _.j'

,J-;,,j /

Alan Robock Department of Meteorology

University of Maryland

College Park

MD 20742

USA

(NASA-CR-191955) PINATUBO ERUPTION

WINTER CLIMATE EFFECTS: MODEL

VERSUS OBSERVATIONS

(Max-Planck-Inst. fuer

Meteoro]ogie) 28 p

G3145

N93-21946

Unclas

0141487

_NNX_

sTi F_LI'I_

E] e._- _ i_A_l :_

MAX-PLANCK-INSTITUT

FOR METEOROLOGIE

BUNDESSTRASSE 55

D-2000 HAMBURG 13

F.R. GERMANY

Tel.:

Telex:

Telem(_il:

Telefox:

+49 (_o)4 I_ 73-0

211092 mpime d

MPI.METEOROLOGY

+49 (40) 4. 11 73-298 REPb 94
Pinatubo Eruption Winter Climate Effects:

Model Versus Observations

H.-F. Graf, I. Kirchner, A. Robock" and I. Schult

Max-Planck-lnstitut fur Meteorologie

Bundesstr. 55, 2000 Hamburg 13, Germany

Department of Meteorology, University of Maryland,

College Park, MD 20742, USA

October 1992

ISSN 0937-1060
Abstract:

Large volcanic eruptions, in addition to the weU-known effect of producing

global cooling for a year or two, havc been observed to produce shorterterm responses in the climate system involving non-linear dynamical processes. In this paper, we use the ECHAM2 general circulation model forced with

stratospheric aerosols to test some of these ideas. Run in a perpetualJanuary mode, with tropical stratospheric heating from the volcanic aerosols typical of the 1982 El Chich6n eruption or the 1991 Pinatubo eruption,

we find a dynamical response with an increased polar night jet in the Northem Hemisphere (NH) and stronger zonal winds which extend down into the

troposphere. The Azores High shifts northward with increased tropospheric

westerlies at 60°N and increased easterlies at 30"N. Surface temperatures

are higher both in northern Eurasia and North America, in agreement with

observations for the NH winters or 1982-83 and 1991-92 as well as the winters following the other 10 largest volcanic eruptions since 1883.
Introduction:

The eruptionof Philippine volcano Mr. Pinatubo on June 15, 199l, was possibly the most sulfur-rich of this century, putting about 20 Mt of S02 into

the Stratosphere (Bluth et al. 1992), which convened to 40 to 50 Mt of

sulfuric acid aerosols. Although conventional wisdom holds that volcanic

aerosols produce cooling at the surface due to reduction of incoming solar

radiation (and indeed that effect is to be expected on a two- to three-year

timescale (Robock 1991)), the winter of 1991-92 (DJF) was warmer than average over North America and the Eurasian middle latitudes both in satellite

derived temperatures of the lower half of the troposphere and in synoptic

surface observations at the same time. Negative temperature anomalies were

observed mainly over the eastern Mediterranean and over the northeastern

part of North America, the Davis strait and Greenland. The question arises

as to whether this anomaly pattern is produced by a deterministic process

other than the reduction of shortwave radiation or is simply a manifestation of the internal climate variability. The possibile effect of the coincidence of volcanic aerosol with an El Nifio leading to the higher continental

winter temperatures is discussed e.g. by Robock (1984) and Graf (1986).

Experiments to study this effect in more detail are currendy ongoing and

will be described in a forthcoming paper.

In the literature some contradictive results have been published about the

effect of volcanic aerosol on the climate of the cold season both from

model studies and observations. Recent numerical volcano experiments with

general circulation models (GCM) (Graf 1992 and Hansen et al. 1992), in

which the additional aerosol was introduced implicitly by reducing the

.solar radiation, shows the strongest climate signal during the warmer

seasons. This is not surprising, since the simple linear reduction of incoming solar radiation as used in these studies has only a minor absolute

effect in high latitudes during winter. Another experiment, mainly focussed

on the stratospheric processes (Rind et al. 1992) and including longwave

and shortwave heating, found decreased tropospheric midlatitude westerlies,

An energy-balance model (Robock, 1984) and an annual cycle GCM calculation

(Hansen et al. 1988, as analyzed by Robock and Liu, in preparation) actually find enhanced cooling in the winter polar regions due to the sea

ice / thermal inertia feedback. However there are observational studies

that show winter warming after large eruptions over western North America
(Lough and Fritts, 1987) and northern Europe (Groisman,I 1985, 1992). A

reanalysis of global temperature data (Robock and Mao 1992, submitted)

reveals winter warming over large parts of Eurasia after all major volcanic

eruptions of the last century since the eruption of Krak:ltau in 1883. The

question thus arises, whether there are processes other than the direct

effect of radiation which determine the climate response to volcanic aerosol during winter. One possibility is the dynamic coupling of stratosphere

and troposphere, as has been shown with simple linear models (GeUer and

Alpert 1980, Schrnitz and Grieger 1980). Numerical experiments (Boville

1986) showed for instance reduced meridional heat transport in the troposphere in the NCAR Climate Community Model for conditions typical for stratospheric westerly circulation (cold polar and warm tropical stratosphere).

We therefore conducted a GCM experiment in which, in addition to the reduction of solar radiation at the top of the model atmosphere, heating rate

anomalies were added in the stratosphere producing a rather realistic temperature anomaly field in the stratosphere. We want to study the atmospheric response to the volcanic disturbance keeping the ocean fixed. Therefore

any transitional effects from slowly varying pans of tlae climate system

can be neglected. The atmospheric memory is at best two weeks. Since our

main interest is in wintertime conditions, the model was integrated in a

"Perpetual January"-experiment mode. Cumulative hydrological effects like

the build-up of an enhanced snow cover and its effect on the Asian monsoon

system are discussed in Graf (1992).

Experimental Setup

The stratospheric aerosol distribution observed in January 1983 in midlatitudes of the northern hemisphere, ten months after the eruption of El Chich6n, is very similar to that seven months after Pinatubo in January 1992.

The transmission of direct solar radiation 10r wavelengths of L=0.3_m to

2.8_tm over Mauna Loa (see Climate Diagnostics Bull., June 1992) is somewhat

reduced in January 1992 ('_=0.86) compared with January 1983 (x=0.89), and

the column aerosol backscauer coefficient at ruby-wavelength (L=0.69 0.m)at

the same place is larger by a factor of two for January 1992. But, in middle latitudes of the northern hemisphere, lidar measurements of the aerosol

column backscatter are comparable for both Januaries (Jager, pets. comm.,

see Figure 1). Later in 1992 the Pinatubo aerosol became more effective in
0.005

"7.

7.

M. 0.004-

uu

I.-

0.003

,,t

(J

m 0.002

.m,

.s

_) 0.001

(J

O92HFGo

,,, i,,t,,ll,,i,,l'i'l''l''I''J''l''l''r[''l''

I I -

I ! I

I I I

69,.3 .- ,I ,, ', ..

o-e-e-e EL CHICHON 19821 I I |

PINATUBO 1901 It , I _

I I -

I I

I I

) I --

t !

0 I

,' ,' /1o / :

z t _ % /_i" 1,.' i _

,, TM ,i

t e t ! t t t I I1-1 I I ! f I. '" .r- , , _ , , , , , , , _ I _ I

i ! ;i i ! I ; i i i i

0 30 60 9_1 120 150 _180 210 240 270 300 330 36p 30 60

94 196 : 305

DAY OF THE YEAR

Figure !: Column back._:.ler IL)r _.=6¢.N.3 nm, Lid.'tr ob_rvations over

Garmisch-Panenkircllcn (47.5"N. II.(FE) after the eruption of El Chichon. April 4, 1982 (open circles) and Pi,mtubo. June 15, ICY)l (full

circles) (from J_lger, 1_)2. o;tendcd and rcdrawn). Day of Year 94

()flh) i._ ll_C d;w ¢_f Iht' rl'_lx't'tivt' L,n,_H(m

extinguishing solar radiation also in midlatitudes by a factor of 2 to 3

compared to the E1 Ghich6n aerosol For our mode] simulation we use.d a

prescribed aerosol distribution according to January 1992 observations,

i.e. similar aerosol like in January 1983 for midlatitudes and stronger

aerosol in low latitudes leading to a reduction in shortwave radiation as

shown in Figure 2 (dashed line). Lacis et al. (1992) found that the only

important parameter for the net radiative forcing due to volcanic stratospheric aerosol is the optical depth. Therefore, lacking further detailed

information about the optical properties of the January 1992 aerosol 0ike

SOLAR FORCING TOP O-" ATMOSPHERE

2{tO0

Figure 2: Reduction of .solar radia(ion IW/m21

due to volcanic aerosol lot Januz, ry 1983 (._)-

lid line) and !_.)2 (broken line). The Ja,_u:try

1992 values are used as solar forcing in |he

expenmems.

IC.00

E O.00

-_0._0

_ El Chichon

• _-- P;nOtuDO

-20.30 , , , . , , ' ,

-90 -70 -50 -30 -10 :O ' _O ' 50 70 90

MODEL _ATITUDE
e.g. effective radii of the panicles), we used model calculations based on

El Chichrn to prescribe these.

The climate model used in our computations is the ECHAM2 GCM, used at the

Max-Planck-lnstitute fdr Meteorologie for climate modeUing research (Roeckner et al. 1992). We used the T21 version (horizontal resolution approximately 5.6*) with 19 vertical levels, the full diurnal cycle and standard

physics, including cloud-radiation interactions. Sea surface temperature

(SST) and sea ice extent are prescribed from climatology. An investigation

of the relative effects of anomalous SST assiociated with the 1991-92 E1

Nifio and the volcanic aerosol is in preparation.

Since the radiation scheme of the ECHAM2 climate model is unable to deal

with variable aerosol concentrations, we decided to use the "anomaly forcing" technique which has previously been used to study the climate response

to the smoke from the burning Kuwait oil wells (Bakan et al. 1991). This

technique allows the inclusion of any radiational effects of trace gases or

aerosols in the model by the use of an additional radiation model, which

explicitely treats the species under investigation. In this case we used

the 8-Eddington model, Bakan (1982), with the extensions for aerosol and

ozone from Schult (1991). We ran this model parallel to the original ECHAM2

radiation code once with and once without the prescribed aerosol with the

atmospheric conditions from the climate model. The flux differences of the

two 8-Eddington computations arethen added to the results of the original

code. After six months of running this complete interactive radiation and

circulation code we obtained mean monthly anomalies of the radiative fluxes

and heating rates. The effcct of the aerosol longwave radiation on the

surface radiation balance was more than one order of magnitude less than

the total effect. It will therefore (also in accordance with Lacis et al.

1992) be neglected in the computations described below.

Since the fully interactive model increases the computation time for ECHAM2

by a factor of five, we avoided its explicit application in order to be

able to run a series of experiments long enough for statistical investigations. Thus, we used the externally .computed mean heating rate anomaly and

the mean reduction of the shortwave radiation (Figures 2a and 3) of the

fully interactive model for a series of perpetual January simulations. The
zo40-

_" 6o!

BO100 -

120

model was initiated with data from one arbitrarily chosen January l of the

control experiment and integrated for 30 days. Doing this 60 rimes (each

new 1st January being initialized with the fields of the 30th January of

the former inte_ation) we produced a series of 60 January simulations both

for the disturbed (Experiment) and the undisturbed (Control) case.

Heating Rates Kelvin/Day

-90 -60 -30 0 30 60 90

, , I , , t , . , . , )_. , . ' _/'

40

60

100

120

%40 . ,_.. _. . , . .._, , 14o

-gO -60 -30 0 30 60 90

Perpetual January Anomalies Temperature (31-80)

-90 -6o -30 o 3o 6o 9o

-go -oo -30 o 3o 60 IN)

Perpetual January Anomalies Zonal Flow (31-80)

-90 -60 -30 0 :10 60 9o

90t \ _f-'_ f-'_J-'_\\\\\\\\",_//_/,/, 90

ii i

I ' " I ' " I ' " I ' ' t

-IO -eO -30 0 30 80 90

lsUtude

Figure 3: Cro_ _ctions of the zonal mean vaJuc_ of:

a) Heating rate anomalies IK/dl, ax determined from the radiative

transfer model for January 1983 or 1992 (top).

b) Temperature ;momalics IKI h)r Perpetual Janu:,ry Experiment (monlh 31

IO 60) minus Control (middle).

c) Zonal wind ;mona;dies Im/sl for Pcrpctual January Experiment minus

control {month 31 to 6()) (tx)ttom).
Figure 4 gives an overview of the mean thermal and solar net radiation

anomalies at earth surface as the climate model produced for three different experiments for January conditions: Perm3 is our Pinatubo experiment,

SST-2 is an El Nifio experiment using the observed tropical (30°N to 30"S)

sea surface temperature anomalies for January 1983, and Perm6 is the combined El Nino / Pinatubo experiment which will be discussed in a succeeding

paper. Anomalies of surface iongwave as well as shortwave radiation are in

the same order of magnitude for all experiments. In all cases there is net

increase of longwave radiation in middle and high latitudes. The pure El

Nitro experiment only leads to a meridional shift of cloud areas and therefore to a global mean zero effect of the shortwave anomaly, while the volcanic aerosol reduces solar radiation in the order of 5 to 7 W/m".

NE" THERMAL RADIATION SFC. ANOMALY

20.00

C'q

E 0.00

RERM6 - CONTROL

• -- _ PERM,3 - CONTROL

..... SST2 - CONTROL

1o.oo

t

" -

- 'O.O0

- 20.30 , ,

-90 -70 -50 30 - 10 _0 30 50 70 90

MODEL _ATI,UDE

NET SOLAR RADIATION S;'C. ANOMALY

20.00

10.00

F: o.oo

- ',o.oo

J '1 ")"

PERM6 - CONTROL

"-I - CONTROL

..... CONTROL

-20.00 -9o'-_o'-_,o*-._o'-'_o ?o'3'o sb'Tb'

MODEL _ATI'UDE

Figure 4: Net anomaly IW/m"l at the surface for thermal (Fig. 4a) and solar

(Fig. 4b) radiation ?or difli:rent Permanent January Experiments. PER.M3 is

the actual volcano experiment. SST2 is forced by tropical sea surface

temperature anomalies .',s observed in January 1983, and PERM6 is a combination of volcanic and SST Ibrcing.s.

We inspected carefully the results for possible trends in the zonal mean of

temperatures and geopotential and found that there were no strong trends

detectable after few months. To check for the stability of any signal, we

perlbrmed analyses with different sub._ts consisting of 30 Januaries each.

So we chose blocks of 30 succeeding Januaries and sets of randomly extracted Januaries, always excluding the first two months which show initial

drift. All subsets gave the same principal patterns with only small varia.
tions in amplitudesand pattern. Here we will concentrateour discussionto

the last half of the ensemble.The permanent January simulations are not

independent from each other. This must be considered in the statistical

tests by reducing the number of degrees of freedom. The effective number of

degrees of freedom was computed for every gridpoint using the autocorrelation function. The resultant Bartels' number B (Barrels, 1935) has a global

mean value of B = t.8 (B = 1.4) lot the 2m temperature (500 hPa height).

This means that the effectivc degree of freedom is N_r = N / B. This is

consistent with the experience that the atmospheric memory is at most two

weeks, Trying to stay at the conservative side of estimates, we made the

assumption that every second experiment January is independent. Consequently for a set of 30 Januaries we use an effective degree of freedom of 15.

This leads to a slight underestimation of the significance of the anomaly

amplitudes over most parts of the gtobe. Only over Antarctica the Barrels'

number is larger, indicating some trend or the model in this region.

Model Results

1. Stratosphere

Although the model stratosphere has only a coarse resolution, the obtained

results are comparable with observations and other model studies.

The zonal mean cross section of the simulated temperature anomaly (,Figure

3b) is in good qualitative agreement with the corresponding observations

from December 1991 (,Figure 5).Positive temperature anomalies are simulated

and observed in low latitudes, and negative anomalies in the polar night

region of the stratosphere. Although the heating rate anomaly is positive

over the North pole (see Figure 3a), the pole cools due to a dynamic effect; because of the intensified polar stratospheric vortex the meridional

(poleward) heat transport decreases. The meridional temperature gradient

between the tropics and the winter pole is strengthened in the lower stratosphere, in January 1992 a minor stratospheric warming occured increasing

the temperature in high northern latitudes by 5 to 6 K. Although stratospheric warmings in principlc can bc simulated by ECHAM2, it is not a predictable feature and is not reproduced by the model in this case.
8

E  ̧

c" 2-

< .

tD 0

k0-2-

1- 4 _

J

,_RASO _RASO 12/82 12/91

uSU 12/82

MSU 12/91

-8 -9o'-_o'-_o'-_o'-_o ,b '3b 'sb 'Tb '90

La:itude

Figure 5: Observed slratosphcric December

lempcralurc anomaly from the 1982-91 mean

alter lhe cruptions of El C-"hich6n (1982)

and Pinatuho (1991), data from Labitzkc

(radiosonde, solid lilts) and Chrisly (MSU,

broken lines).

Our results are also consistent with the study of Rind et al. (1992), who

used the GISS model with an increased number of model layers (23 layers,

top at 0.001 hPa). They investigated the effect of a homogenous aerosol

layer of optical depth 1:=0.15 on the zonal mean conditions in the middle

atmosphere.

The dynamic effects in the stratosphere are primarily due to changes in the

meridional temperature distribution. The mean zonal wind anomaly (Figure

3c) shows an enhanced polamiglat jet in the stratosphere with anomalies up

to I0 m/s. These westerly anomalies penetrate also downward into the troposphere resulting in a i:x)leward shift of the upper tropospheric westerly jet

of the northern hemisphere.

2. Troposphere

The model simulates a significant reduction of the variance of the geopotential height of the 200. 5(X) and 850 hPa levels over the whole North

Atlantic region. This holds also for the variances of thc zonal wind field

and temperature (the latter mainly in the lower troposphere) for the complete set of the experiment (Figure 6) as well as for sub_ts. This is especiaUy clear for the geopotential over the Azores high, over Labrador and

Southern Greenland. No significant signal is detectable elsewhere applying

the F-test. or the variance is increased in the experiment (Figure ¢_).

I0
Geopotential Height 850 hPa gpm

,cow,5ow,_o,,o, cow _o, o _oE co,.OOE ,2o_.,50E,80E

, , I , , I , , I , , I,.I , . I. , I . . I , , I . , I . . I , ,

. . .....

30N . o._. +_, -_--,.,._u_ _ . - ,_ + "

........ T.... _---- "_,"_ '-" _ -_" .... t.44.,_"'-' .

o.as_ _ i/-_ - ,_.61 _ -

3os.... '. -----_,_'- .......... _4 • "

"1_-_ ....,-.... _: "-,.i_--c_._ _

_--_' ,.._.<::j,.+/ ) -- _,._

I'' [ ' ' [ ' ' r-_ ' I ' ' i ' ' I ' ' I ' ' I _ I _ I ' ' I ' '

!80WISOWI20W 901 601 30W

- 60N

- 30N

0

30S

I eos

0 30E 60E g0E I20E 150E 180E

Temperature 850 hPa Celsius

180W 15011' 120W 901 (101I 30W 0 30E 60E 90E I20E 150E 180E

-,-,, I,, I,, I,, I,, I,, I,, I i, I,, I,, I,, I , ,

U'.Igi__++_.._+._.- • _b,+,:. __ _-

, ;_..-'_" ___1_,_'-...._ ,-- _ .....

........ '_r'--+''_+" _':_"zz_ .,_ "J_"_L'['- +lZ_,.'_...,.."'_,_ _E" -.,_+". " "'" "_. ' _.

60N - :.: I':M _ " " ..,,';. ,t .._..__/ _.: o.es.,=,.. " " - 60N

.].64 ! _.em /,.:'_,;_- ., :,. :

3o_..... :.... _V:_)_-..;:_u+----."::- "_'_ "+ ; " . 3o.

_ _.a_ o. .,w4-._ll "t - ' " " " "*" _ " ''"_"

• _ _ u._ d-", k.W'_'" 0 .............. .,--_ ........... o.,_. ' o-.-.

o++ _'_ '-_ +,, ,,. ,4., ,....,4 : "" .'.-_ 0

' - o e- _ /,,_,.:_,.._, =o , .0.6_ <.7 "_'-" "' U

' .'.¢" __m.-.,,,--"J:/X_Jl- _ .'_. o:__ ;.. -

301_

, ,+ ,

301S

6o_-L . '--J_., ,_-__-_._L-L.... : ...... Z,/.. /:--,_ -eos

,__:_:: : ,.. _ _.. \ _.F_,

I ' ' I ' ' I ' ' I ' " I ' ' I ' ' I ' ' I ' ' I ' ' I ' ' ] ' ' I ' ' ,

180W!50WI20W O0W 60W 30W 0 30£ 80E 90E 120E ISOE 180E

U-Velocity 850 hPa m/s

180W 150W 120W 9011' 60W 30W 0 30E 60E OOE 120E 150E 1801

,, I , ,t, , I, , I,,I ,,l,, I , ,I _ _ I +, , I _ ,I, ,

:_._: o+4_ +...... :- .... o._-

• " '_;':'_'" ..... :_tN! ° _ " " ' ":" "

6oN _:_s_--_,"'-; --:__+: _ +:..-.' • - 6oN

_o. o._. --:.. :'-'__ ":::,,+-,+,%;.,,+ _ -._o.

-. I.._,+,:'__-_.0.1t4 _ _)" _.., _.,,. ,: ,_-

"_.---.'_'0_" ' . o . - _ :. ._" ._:s4_+ o.+s --

o.M, !,i4-- • . , ,._ _.s4. \;_% , _4

0 ;_ 0.1_ -:-" /_.'_-_ ; " ' _(;j_ .... ._AK-:K._.._ 0

o.oH'__ o.-'='+_.,- "+: , - '_.: '_---__ o

• '_ _ ," ; 'IIP_5.' ' ;

3oe- :........... _mm.. . ,, _ -

,.. _ .5, o.,, _. oos

• o._ 1"" _;_.+4. . ..._ .__x._.sY

6os .... ,_u_. •_2_. _ - eos

I ' ' ' ' I ' ' I ' ' ! ' ' ' I ' '

180W I50W 120W g0w 60W 30W 0 30E 60E g0E 120E I50E 180E

Figure 6: Changcs in local varinncc (F-tcst) a! Ihc 8..,¢0hPn Icvcl ['or:

Gcopolcnlia! hcigh[ (top), zonal vclocity (middlc), and [cmpcratum (bottom) for Ihc Pcrmancnl January ExpcrimcnL Shading indica[cs xigni_can!

diffcrcnccs bclwccn Exrx:rimcn! ;rod Control. Valucs I;]rgcr (smailcr) Lhan

I mcans [hat [hc v_r,;mcc in [hc cxl'x:rimcn! is lar_cr (smaJlcr) than in

thc Conlrol.
A local T-test was used to test for the simaificance of the amplitudes of

the zonal wind, geopotential height and 2m-temperature anomalies in different sets of 30 experiment months. Here we will focus on the results for

the Januaries 31 to 60, i.e. on the second half of the experiment. Principially, the anomaly patterns were similar in the other sets of each 30

Januaries.

In addition to the t-test, and in order to test for the difference of the

distribution functions of the Control and the Experiment, a recurrency

analysis (von Storch and Zwiers, 1988) was also applied to several sets of

Januaries. T-test and recurrency test give comparable results, especially

concerning the areas of statistically significant signals. In each set the

pattern of anomalies of the geopotential heights is similar throughout the

troposphere.

The largest amplitudes of the geopotential heights in all tropospheric

levels occur over the North Atlantic, where the pattern indicates a northward shift of the tropospheric circulation (Figure 7). The Azores high is

displaced northward by about 15° and the Icelandic low shifted towards

Greenland. This pattern is significant (15=0.05) following the T-test and

recurrent in all sets as well as are the negative tropical geopotential

anomalies, mainly in the middle troposphere. The most prominent signal is

the anticyclonic anomaly in the area of the northward shifted Azores high.

Thus, the area of strongest and recurrent circulation anomalies is the

Nonh Atlantic. This holds for the whole troposphere. The structure of the

anomaly patterns throughout the troposphere suggests a barotropic response.

It is suggested that the interaction of the strengthened stratospheric

polarnight jet with tropospheric lec effects of the Rocky Mountains is

responsible for the strongest response of the circulation to the aerosol

forcing over the North Atlantic. The primary effect may be a change in the

stationary planetary wave pattern.

The westerlies over the North Atlantic are shifted northwards resulting in

westerly wind anomalies around 60" N and easterly anomalies at 30°N (Figure 8). This pattern is again stable throughout the whole troposphere. It is

also combined with a significant reduction of the variance of the zonal

wind component as is obvious from the F-test. Accordingly the Nonh Atlan12
tic centers of action shifl northward.

Geopotential Height 200 hPa gpm

1B011'1601F 1JIOIF flOW 8011' _ O _ _ _ lJlOE 160E IBOE

Geopotential Height

1_ _ LI0W 901

,,|,, -- I

500 hPa gpm

BOW' 301F 0 _ _ IN_ J.ilOE150111BQE

,,I,,I,,I,, ,I,,I,,l,,I,,l, I

._. .',...-a_ .._.:'.---7 -.-.-.::.::..-.- ...... _ 3lira'

] I /:ii!

Geopotential Height 850 hPa gpm

LBOW'LBOIFLROTOOW O01f SOW' 0 301 _ flOE LSOELBOEIB01

I , , I , , I , , I , , I , , I , I , ! , , I , , I , i I., , I , ,

........i _ '- ! i')/'" ,

_.!:!':":.. - - ; ....._..._

I ' ' I '" i ' ' I" ' I ' ' I ' ' ' I ' ' J ' ' I ' ' I ' '

18011'1601F lll01r gOlF BOW ll011' 0 _ _ flOE UlOE 1501l IBOZ

Figure 7: Anom:dics of Ilyc gcorx)lcnlial hcighl at dif[crcnl Iroposphcric

Icvcls. Sh.'_din_ indicatcs statisLically significant (error chancc 5%)

differcnccs bclwecn Control :_nd Experiment for months 3J to 60 on thc

basis of a local t-lest, dcgrccs of freedom rcduccd Io 15.
U-Velocity 500 liPa m/s

I Jr_

• _ i g_--a.___=k '_ "-_-__1._ " --_

. " ._;i .......... I - ,4 _

_''1' '1''1''1''1'11 '1' '1' 'l''[''i''l

1801r J 1,110l Ilolr BoW _ 0 801 601 _ 111011lO01 tll_l

0

808-

8011-

U-Velocity 200 hPa m/s

1B011'150V 111011'gO11' _ _ 0 _ _ 9_ llOE UlGI 1BOI

18Olr1801ulov gol _ _ o _ 6ol _ Uloze 111o!_L801

Figvre 8: Anomalics of thc znn:ll wind componcni at diffcrcnt trt4_osphcnc

Icvcl._. Sh;=ding in(lic:uc._ ._tati._lic;dly significant (crror chancc 5%)

diffcrcnccs bClwccn Control :rod Exl_rinlcnl for months .31 Io 50 on thc

basis of a local t-tcst, dc_rccs of Irccdom rcduccd to 15.

14

OR_GIN_L P)=C_ IS

OF POOR QUALITY
Temperature 500 ]_a Celsius

IOGW__ _ _ 801F 0 _ _ _ lmll_lBOZ

I,,I,,I,,I,,I,,I, I ,I, ,I, ,I,,I,,I,,

.:..:..,-., o a@ ), ; '_ :_.,/T I''_ _ o

• .:._'._;.'._; ........ j..-_.;;;..;...::. ..::'_; o---,_ _... _...,.,._,_. - ._-_.-.-.,.,_-::.._:::::.: . -:..-.:-_ .t-_:..1 f ......":_",._-.-:_-_ _

Temperature 850 hPa Celsius

1BOW 1501F JJIOIF gOlF 6olr _ 0 3Oll _ _ 1sloe 15oE 1B_

, ,I, , l, , !,, I , I ,, I i ,

o g_l o o

I :'_. ..--i_:.--_, _-

'_-.-.-.._:--:i___2o_ I 7 • f ._

2m Temperature Celsius

1BOIF 1501F 1JlOIF BOW BOW 301W 0 _ _ 9011 l_O_ 150B 1BO_

_ ,I =, _ e , I .... I,, I,, I ,, I ,, I,,

I_, -, , ,

.

0 0

........ __.. __._-_--_. __.... __-_

_ | I i *

1BOW 15olr _ Bow _ _ 0 3011 IIOB 9011 180_ 150_ IBOle

Fig.re 9: Anonlalic._ of the air |eml-x:raturc in the lower tm[x).'_pherc. Shading indicatc,,; slatis|ic;lll), signii]¢ant (error chance 5%) di ITcmncc._

between Comrol and Ex[x_rimcn! [or n'mnth._ 31 to ()(} on thc ba_i,_ o[ a

local L-I.C._L.dc_rccs of Ircedom rcdu¢cd (l_ 1,5.
In the tropics, the anomalies of the zonal wind component change sign from

the lower to the upper troposphere. In both sets over the eastern equatorial Pacific in the lower (upper) troposphere easterly (westerly) wind anomalies are found. This indicates a slightly strengthened Walker circulation

over the eastern equatorial Pacific during winter.

The meteorological parameter which normally is of interest in connection

with violent volcanic eruptions is the surface air temperature. Figure 9

0ower panel) shows the mean anomaly of this parameter, again for the last

30 Januaries of our experiment. The main features, however, are similar

also for other subsets of the 60 January integration after few months of

initial drift. The T-test shows a major significant area of cooling over

North Africa and the Middle East, strong cooling over Greenland and parts

of China. In the .southern hemisphere the only significant signals are a

cold southern tip of South America and Antarctic cooling. Significant warming only occurs in the Barents Sea. Here one of the model's restrictions

is apparent, i.e. the fixed .sea ice coverage and tempentture. This in part

prevents the development of positive anomalies of the 2m temperature. In

the free troposphere (e.g. in 850 hPa, Figure 9, middle panel), where the

surface conditions are less important and the transport processes become

dominant, the warming is significant over large pans of the Arctic ocean

north of Siberia and the warm anomaly reaches further .south over Asia.

Tropospheric Observations

In order to compare our simulation of the impact of the Pinatubo (or an

other similar volcanic eruption) stratospheric aerosol on climate with

global observations data, we used the 500 hPa geopotential height anomalies

as analyzed by the US National Weather Service, the satellite observed

lower tropospheric tcmperature from the Microsounding Unit (MSU) (see Spencer and Christy, 1990), and (to a certain degree) synoptic observations

collected by the Seewettcramt of the Deutscher Wetterdicnst. The last data

set is based on a non-uniform distribution of stations with large gaps over

the oceans. This gives serious problems when the data are interpolated to a

coarse map such as the T21 grid. Thus these data are used only for qualitative comparisons.

The MSU temperature data use the tcmper:lture dependence of the microwave

]6
emission of oxygen molecules. For the tropospheric signal, calibrated MSU

channel 2 (53.74 GHz) brightness temperatures are combined with a weightrung function describing the vertical sensitivity to thermal emission by

molecular oxygen in the atmosphere. The main part of the signal we used is

from the layer between the earth surface and 400 hPa (with a peak at 700

hPa), i.e. from the lower pan of the troposphere rather than from the

earth surface. The spatial resolution of these observations is in the order

of 10°, even though the data are published on a 2.5 ° grid (Spencer et al.,

1990).

In order to compare with previous violent volcanic eruptions, the surface

temperature data set from Jones et at. (1988) was also used (see Robock and

Mao, 1992 submitted).

An inspection of the mean December 1991 through February 1992 500 hPa height anomaly (Fibre 10) shows some similarity between our simulation and

observed anomalies: The strongest anomalies occur in the Atlantic region.

9ON"

DJF 1991-92

Observed 500 hPa Oeopotential Anomaly (m)

-':---:"...... Z

30S

gos

180

Figure !11: Anomaly of the 500 hPa gcolx)tential height as analyzed by II_c US

NWS. mean of December 1991 Ihrough February 1992. Anomalies calculated

from mean 1982-1990.

17
The simulatedanomaliesare generally weaker than the observedones and a

tendency to a westward displacementis characteristic.Such zonal (westward) displacements(compared with observations)of the mean circulation

features are common in low resolution GCMs, resulting mainly from the

smooth orographic representation(Sausen 1992, pet's, comm.). The negative

anomalyover Greenlandexists in both fields (for the simulatedanomaly see

Figure 6), while the observedcenter of the positive geopotentialanomaly

is shifted eastwardsto England.Both, in the observationsand in the simulation, a trough develops over Eastern Europe/ Western Asia, which also

has a more westerly position than in the observations. The pattern anomaly

correlation (Table 1) is between 0.37 for the globe and 0.40 for the northern hemisphere if we use the mean pattern of the experiments 03-60. As

stated above, subsets have shown the same patterns. This means that only

about 15% of the total North hemispheric observed spatial variance of the

500 hPa field of the middle troposphere in winter 1991-92 can be explained

by the volcanic impact alone. Much better results we get for the western

hemisphere (180°W to 0°) where this value can be improved to 49% (r=0.70).

Table 1: Pattern correlation coefficients (area weighted, multiplied

by 100) between observed and simulated meteorological fields. A

1:1 correlation means that no spatial shift is allowed, "optimal" correlation allows for 2 gridpoints shift in all directions

to get a best fit. All fields were interpolated to the T21 grid

before being analyzed. Mean of the experiments 03-60 used.

Parameter global N-hcmisph. 90°W-90°E 0°-180"E 0°-180°W

M S U-Tern p.

1: 1 18 57 39 54 60

optimal 29 57 62 56 63

! 500 hPa

1:1 35 39 1 -29 68

i optimal 37 40 41 -46 70

Over Siberia the patterns of 5(X) hPa geopotential height anomalies do not

agree well because of a simulated high over North Siberia which is not

observed in winter 1991/92.

For thc lower tropospheric temperature (Figure 11), reasonable agreement

between simulations and observations in winter 1991-92 is found for the

18
northern hemisphereagain. The "optimal" pattern correlations are between

z=0.56 and x=0,63 for different subareas of the Northern hemisphere (Table

l), resulting in z--0.57 for the whole northern hemisphere. This means that

about 35 to 40% of the total spatial variance of the obscrved 2m temperature can be explained alone by the volcanic forcing. Thc model reproduced

especially well the cooling over Greenland and North Africa. Some of the

warming of the middle and high latitudes of the northern continents, which

is seen in both patterns, is also modelled at the fight place, although not

significant with respect to the amplitude because of the high variability

of the temperature in these latitudes in the control run.

Lower Troposphere Temperature Anomaly (K)

9_ / .... _._) ._ ._ ........................ ".., ---...._.. t ------- ,... " .,,.).- .......

Figure 11: Anomaly of Ihc air Icmpcralurc near surfacc I K I ,'Ls ob_rvcd

Dcccml:x:r 1991 through February 1992 from MSU _tcllilc data (Spencer and

Christy. 10()1). Anom;dics calculalcd from mean for 1982-1990.

The good agreement between the model simulation and the observations in

winter ]999]/92 still could be coincidental. This might be excluded if

similar anomaly patterns evolved after other volcanic eruptions as well.

An examination of the near surface temperature data set from Jones et al.

(1988 and 1991, updated, Jonc._, pcr_ona[ comm.) after the twelve most violent volcanic eruptions of the last century (Figure 12) also shows good

agreement with our simulations. For the six tropical volcanoes (Krakatau,

Sta Maria. Agung, Fuego. El Chich6n and Pinatubo) the anomalies are for the

]9
u_r 3urtGce tempero[ure _nornu_y k_)

After 12 Largest Volcanoes 1885-19£2

30N

60S,

Figure 12: Ob._rvcd December Ihrough February mean Icmpcralurc anomaly

aflcr the twelve largcsl volcanic eruptions of Ihc last century (Data from

Jones. It)88. UlXJ"alcd, Jones, i_rs. com,n.).

following boreal winter, for the six high latitude volcanoes (Novarupta,

Tarawera, Bandai, Ksudach, Cerro Azul and Bezymiarmy) anomalies are for the

second winter after the eruption to allow time for the aerosols to reach

the tropical stratosphere. The anomalies are calculated from the 1951 to

1980 meanas differences from low pass (cut-off 10 years) filtered data. The

pattern after every volcano is virtually the same as the mean effect shown

in Figure 12, except some diffcrences over North America for winters with

and without El Nifios. The best pattern correlation between the composite

signal of post eruption winter temperatures and our simulation is found for

the sector between 90°W and 90°E on the Northern hemisphere (x=0.35). There

is a zero pattern correlation for the western half of the North hemisphere

and only z=0.21 for the eastern half.

Conclusions

The perpetual January GCM simulation has shown that the stratospheric aerosol produced by a violent and sulfur-rich volcanic eruption can alter the

stratospheric and tropospheric climate conditions significantly. Using the

20
"anomaly forcing" technique the stratospheric forcing due to the aerosol

(which is not included in the original GCM) can be simulated in accordance

with other model simulations and observations. The stratospheric winter

signal can be described as a sigrtificant warming of the aerosol containing

layers in the tropics and midlatitudes due to absorption of radiative energy. This warming increases the poleward meridional temperature gradient,

leads to an increased polar night jet, reduced meridional heat transport

and therefore a cooling tendency in the stratosphere over the winter pole.

In this study the tropospheric signal was found to be meridionally nonuniform. It is strongest over the North Atlantic and consistsof enhanced

westwinds over high latitudes (around 60°N), strengthened zonal heat and

moisture transport in higher latitudes and a northward shift of the Azores

high. The areas of strongest anomalies also show significantly reduced

variance of the parameters.

A temperature anomaly field was found in the lower troposphere, which (in

the main features) is very similar to observations during the winter 199]-

92 Similar anomalies were observed also during other post volcano winters. There is a significant cooling over Greenland and North

Africa / Middle East. In agreement with observations some warming tendency

with a maximum amplitude of about 2 K is seen in the mean results over

North America and northern Eurasia. This positive temperature anomaly can

not be proved to be statistically significant because of the large natural

variability of surface temperature in these latitudes. This warm anomaly is

supported by observations following historic eruptions (Robock and Mao

1992) and it determines the high pattern correlation found between simulation and observations.

The simulated and observed temperature anomalies following a violent

volcanic eruption in part contradict the findings based on linear multiple

regression analysis (Cress and Sch0nwicsc 1991) and encrgy balance models

(Robock 1984). Thc results undcrlinc thc importance of nonlinear and

dynamic responses of the atmosphere to external forcing.

21
ACKNOWLEDGMENTS: We thank John Christy, Phil Jones and Karin Labitzke for

providing the temperature data. We are grateful to M. Grunert for preparing

Figure I, to Jiaoping Mao for producing the analyses in Figure 12 and to

Brian Doty for the GRADS software for Figures 10 to 12. We thank Huug van

den Dool and John Janowiak for the 500 hPa data. Data are from the Climate

Analysis Center of the National Meteorological Center of NOAA. Comments of

Lennart Bengtsson and Mojib Latif on a first version of the manuscript

helped clarify the presentation. This work was partially supported by the

Bundesministerium fi.ir Forschung und Technologie (No. 07-KFT-86_), NSF

grant ATM-8920590 and NASA grant NAG 5-1835.

22
Literature:

Bakan,S., A. Chlond, U. Cubasch, J. Feichter, H. Graf, H. Grassl,

K. Hasselmann, I. Kirchner, M. Latif, E. Roeckner, R. Sausen, U.

Schlese, D. Schriever, I. Schult, U. Schumann, F. Sielmann and W.

Welke, 1991: Climate response to smoke from the burning oil wells

in Kuwait. Nature, 351, 367-371

Bakan,S., 1982: Strahlungsgetriebene Zellularkonvektion

twolken. Dissertation, Univ. Hamburg, 99 S.

in SchichBartels,J., 1935: Random lluctuations, persistence, and

sistence in geophysical and cosmical periodicities.

Magnetism atmospher. Electr., 40, 1-60

quasiperTerrestr.

Boville,B.A., 1986: The influence of the polar night jet on the

tropospheric circulation in a GCM. Joum. Atmosph. Sci., 41, 1132-

1142

Bluth,G.J.S.; S.D.Doiron; A.J.Krueger; L.S. Walter and C.C.Schnetzler, 1991: Global Tracking of the SO2 clouds from the June, 1991

Mount Pinatubo Eruptions, Geophys.Res.Lett., 19, Jan. 1992

Cress, A. und Ch.-D. Sch0nwiese, 1990: Vulkanische Einflflsse auf

die bodenn_e und stratospN(rische Lufttemperatur der Erde. Ber.

Inst. f. Meteoroiogie und Geophysik Univ. Frankfurt/Main. Nr. 82,

148 S.

Geller,M.A. and J.C. Alpert, 1980: Planetary wave coupling between

the troposphere and the middle atmosphere as a possible sunweather mechanism. J. Atmos. Sci., 37, 1197-1215

Graf, H.-F., 1986: On El Nifio/Southcrn Oscillation and northern

hemispheric temperature. Gerl. Beitr. Geophys. 95, 63-75

Graf, H.-F., 1992: Arctic radiation

ty. Climate Dynamics, 7, 19-28

deficit and climate variabiliGroismann,P.Y., 1985: Regional climatic consequences ol volcanic

eruptions. Meteorology and Hydrology, No.4, 39-45. [in Russian]

Groisman, P.Y., 1992: Possible region.'ll consequences of the Pinatubo eruption: an empirical approach. Geophys. Res. Lett., 19,

1603- 1606

Hansen.J.: Lacis,A.; Ruedy,R. and M.Sato, 1992: Potential climate

impact of Mount Pinatubo eruption. Geophys.Res.Lett., 19, 215-218

Jager, H., 1992: The Pinatubo cloud observed by Lidar over Garmisch-Partenkirchen. Geophys. Rcs. Lett. 19, 191-194

Jones, P.D., 1988: Hemispheric surface air temperature variations:

recent trends and an updatc to 1987. J. Climate, I, 654-660

Jones, P.D., T.M.L. Wigley, C.K. Folland, D.E. Parker, J.K. An23
gell, S. Lebedeff and J.E. Hansen,1988: Evidencefor global warming in the pastdecade.Nature,332,790

Lacis, A., J. Hansenand M. Sato, 1992: Climate forcing by stratospheric aerosols. Geophys. Res. lett. ,19, 1607-1610

LoughJ.M, H.C. Fritts, 1987: An assessment of the possible effects of volcanic eruptions on North American climate using treering data, 1602-1900 A.D., Climatic Change, 10, 219-239.

Rind, D., N.K. Balachandran and R. Suozzo, 1992: Climate change

and the middle atmosphere. Part II: The impact of volcanic aerosols. J. Climate, 5, 189-208

Robock,A.D., 1984: Climate model simulations of the El Chich6n

eruption. Geofisica Intemacional, 23, 403-414

Robock, A.D., in Greenhouse-Gas-lnduced Climatic Change: A Critical Appraisal of Simulations and Observations (ed, M.E. Schlesinger), 429-444, (Elsevier, Amsterdam,1991).

Robock, A. and J. Mao, 1992: Winter warming from large volcanic

eruptions. Submitted to Geophys. Res. Lett.

Roeckner, E., K. Arpe, L. Bengtsson, S. Brinkop, L. D_imenil, M.

Esch, E. Kirk, F. Lunkeit, M. Ponater, B. Rockel, R. Sausen, U.

Schlese, S. Schubert and M. Windelband, 1992: Simulation of the

present-day climate with the ECHAM model: Impact of model physics

and resolution. MPI Report No. 95

Schmitz,G. and N. Grieger, 1980: Model calculations of the structure of planetary waves in the upper troposphere and lower stratosphere as a function of the wind field in the upper stratosphere.

Tellus, 32, 207-214

Schult,l., 1991: Bildung und Transport yon Aerosolteilchen in der

Stratosphere und ihre Bcdeutung f_ir den Stratdungshaushalt. Examensarbeit Nr. 11, Max-Planck-lnstitut f_ir Meteorologie Hamburg

Spencer.R.W.,J.R. Christy, and N.C. Grody 1990: Global amaospheric

temperature monitorin_ with satellite microwave measurements:

Method and results 1979-84, J. Climate, 3, 1111-1128.

Spencer,R.W. and J.R. Christy, 1991: General Temperature update to

31 October 1991. NASA Headquarters, Washington. (pers. comm.)

Storch, H. yon, F. W. Zwiers, 1988: Recurrenceanalysis of climate

sensitivity experiments. J. Climate, !, 157-171.

24
MPI Report-Reference: *Reprint,

Report No. 1-60

Report No. 61

March 1991

Report No. 62

April 1991

Report No. 63

May 1991

Report No. 64

June 1991

Report No. 65

June 1991

Report No. 66

June 1991

Report No. 67

July 1991

Report No. 68

July 1991

Report No. 69

August 1991

Report No. 70

September 1991

Report No. 71

September 1991

Report No. 72

September 1991

Report No. 73

September 1991

Please order the reference list from the MPI Hamburg

Three-Dimenslonal Simulation of Cloud Street Development du_ng

a Cold Air Outbreak

Anclreas Chlond

Modelling the low-frequency sea surface temperature varfabillty

In the Northern Pacific

Ute Luksch, Hans yon Storch

*Journal of Climate, 1992

How much predictive skill Is contained In the thermal structure of an OGCl

Mojib Laflf, Nicholas E. Graham

• Journal of Phys. Oceanography

Downscaling of Global Climate Change Estimates to Regional Scales:

An Application to Iberian RalnfaU in Wintertime

Hans von Storch, Eduardo Zorita, Ulrich Cubasch

*Journal of Climate, 1992

Southern Ocean Sea-lce Simulations forced with operationally derived

Atmospheric Analyses Data

Achim StOssel

The Phase of 30- to 60-Day Oscillation and the Genesis of Tropical Cyclone_

In the Western Pacific

Hans yon Storch, Ann Smallegange

Time-Dependent Greenhouse Warming Computations wlth a Coupled

Ocean-Atmosphere Model

Ulrich Cubasch, K]aus Hasselmann, Heinke H6ck, Ernst Maier-Reimer,

Uwe Mikolajewicz, Benjamin D. Santer, Robert Sausen

On the sensitivity of the Global Ocean Circulation to Changes in the Surface

Heat Flux Forcing

Ernst Maier-Reimer, Uwe Mikolajewicz, Klaus Hasse_mann

Acoustic Detection of Greenhouse-lnduced Climate Changes In the

Presence of Natural Variability

Uwe Mikolajewicz.ErnstMaier-Reimer.Tim P. Barnett

The Effect of a Regional Increase In Ocean Surface Roughness on the

Tropospheric Circulation: A GCM Experiment

Uwe Uibrich, Gerd BUrger. Dierk Schriever, Hans yon Storch,

Susanne L. Weber, Gerhard Schmitz

A differential absorption Ildar system for high resolution water vapor

measurements in the troposphere

Jens B6senberg

Coupling an Ocean Wave Model to an Atmospheric General

Circulation Model

Susanne L. Weber, Hans yon Storch, Petro V'n'erbo, Liana Zambresky

Climate Varlablllty In a Coupled GCM. Part I: The Tropical Paclflc

Mojib Laflf, Andreas Stert, Ernst Maier-Reimer, Martina M. Junge

-I-
MPI-Report-Reference: *Reprintedin

Report No. 74

October 1991

Report No. 75

October 1991

Report No. 76

January 1992

Report No. 77

January 1992

Report No. 78

January 1992

Report No. 79

January 1992

Report No. 80

March 1992

Report No. 81

March 1992

Report No. 82

May 1992

Report No. 83

May 1992

Report No. 84

May 1992

Report No. 85

June 1992

Report No. 86

June 1992

Report No. 87

August 1992

Report No. 88

August 1992

Report No. 89

September 1992

On the Space-Time .Structureon ENSO

Mojib Latlf,Andreas Sterl,ErnstMaier-Reimer

On the Navigation of Polar Orbiting Meteorological Scatellltes

JarusJav K]okocnik, Peter SchlL_ssel,Jan Kostelecky, Hartmut Grassl

A Hybrid Coupled Tropical Atmosphere Ocean Model:

Sensitivities and Hlndcast Skill

Richard K]eeman, Mojib Latlf, Moritz Fl0gel

A Direct Mathernatlcal Approach to Optimize the Age-Depth Relation

of Deep-Sea Sedlment Cores

Wolfgang BrOggemann

The Joint Normal Modes of the coupled Atmosphere-Ocean System

observed from 1967 to 1986

Jin-Song Xu

Orbital tuning of martne sedlmentary core= An automatic procedure

based on a general llnear model

Bj6m Grieger

Application of ocean models for the interpretation of AGCM experiments

on the climate of the last glacial maxlmum

Michael Lautenschlager, Uwe Mikolajewicz, Ernst Maier-Reimer,

Christoph Heinze

Normal Modes of the Atmosphere as estimated by Prlnclpal Oscillation Patterns

and derived from Quasi-Geostrophic Theory

Reiner Schnur, Gerhard Schmit-z. Norbert Grieger, Hans yon Storch

On regional surface fluxes over partly forested areas

Martin Claussen, Wim Klaassen

On the cold start problem In transient slmulatlons

with coupled atmosphere-ocean models

K]aus Hasselmann, Robert Sausen, Ernst Maier-Reimer, Reinhard Voss

Iron fertillzatlon of the Austral Ocean

- a model assessment

Katharina D. Kurz, Ernst Maier-Reimer

"Global Biogeochemical Cycles

A snow cover model for global climatic slmulatlons

Betlina Loth, Hans F. Graf, Josef M. Oberhuber

Complex principal oscillation pattern analysis

Gerd B0rger

Scale aggregatlon In sernl-smooth flow

Martin Claussen

Optimal fingerprints for the detection of time dependent climate change

K]aus Hasselmann

Blomes computed from simulated cllmatologles

Martin Claussen, Monika Esch

-2-
MPI-Report-Reference: °Reprinted in

Report No, 90

September 1992

Report No. 91

September 1992

Report No. 92

October 1992

Report No. 93

October 1992

Report No. 94

October 1992

New generation of radiation budget measurements from space

and their use in climate modelling and diagnostic studies

L.DQmenil,E,Raschke

Modal structure of vartatlom in the tropic climate system

Part I: Observations

Mojib Lallf, Timothy P. Barnett, Keisuke Mizuno

The K6ppen Cllmate Classification

as a Dlagnostlc Tool for General Circulation Models

Ulrike Lohmann, Robert Sausen, Lennart Bengtsson, Ulrich Cubasch

Jan PerlwJtz,ErichRoeckner

Simulation of the present-day climate wlth the ECHAM model:

Impact of model physics and resolution

E.Roeckner, K. Arpe, L.Bengtsson,S.B_inkop,L D_menll,

M. Esch,E.Kirk,F.Lunkeit,M. Ponater, B.Rockel, R.Sausen,

U.Schlese,S.Schubert, M. Windelband

Plnatubo Eruption Winter Cllmate Effects:

Model Versus Obr_vatlon$

Hans-F. Graf, Ingo Kirchner, Alan Robock, Ingrid Schult

-3-