UNCERTAINTIES IN

GREENHOUSE GAS INDUCED

/4berta ENVIRONMENfT

UNCERTAINTIES IN

GREENHOUSE GAS INDUCED

CLIMATE CHANGE

by:

Madhav L. Khandekar

Consulting Meteorologist 52 Montrose Crescent Unionville, Ontario L3R 7Z5 E-mail: [email protected]

Prepared for

Science and Technology Branch Environmental Sciences Division Alberta Environment 9820-106 Street Edmonton, Alberta T5K 2J6

March 2000 Pub. No. T/522 ISBN: 0-7785-1051-4

Although prepared with funding from Alberta Environment (AENV), the contents of this report do not necessarily reflect the views or policies of AENV, nor does mention of trade names or commercial products constitute endorsement or recommendation for use.

For further infomriation regarding this report, contact:

Infonmation Centre Alberta Environment Main Floor, Great West Life Building 9920-108 Street Edmonton, Alberta T5K 2M4 Phone: (780)944-0313

This report may be cited as: Khandekar, M.L., 2000. Uncertainties in greenhouse gas induced climate change. Report prepared for Science and Technology Branch, Alberta Environment, ISBN 0-7785-1051-4, Edmonton, Alberta. FOREWORD

The increase in the atmospheric concentration of greenhouse gases is an important global issue because of the risk of climate change. Despite considerable uncertainties in climate change science, the nations of the world have decided on a global effort to reduce the emissions of greenhouse gases. However, for the assessment of future actions and the formulation of an adaptation strategy, an updated knowledge of the developments regarding the uncertain aspects of the science of climate change is essential.

The science of climate change has been advancing in a fast pace. The purpose of this study is to review the recent scientific literature with particular emphasis on the latest development in several key areas of climate science uncertainties. These key areas are noted in the report. Since climate change science has an extensive literature, this allows focus on the more important aspects of the issue.

The report was prepared by Dr. Madhav Khandekar, who is a former Research Scientist fi-om Environment , where he worked for over 25 years. He holds M.S. and Ph.D. degrees in Meteorology from the Florida State University, USA, and has been in the fields of cHmatology, meteorology, and oceanography for over 42 years. He has published over 75 papers, reports, and book reviews, and has written a book on Ocean Surface Wave Analysis and Modelling, published by Springer-Verlag in 1989. His current research interests are: El Nino/Southern Oscillation, Asian Monsoon, and global weather anomalies.

Raymond Wong, Ph.D Project Coordinator Science and Technology Branch Environmental Sciences Division Digitized by the Internet Archive

in 2015

https://archive.org/details/uncertaintiesingOOkhan ABSTRACT

The purpose of this study was to assess the uncertainty associated with the science of climate change and global warming and determine the present state of our knowledge. The study examines several issues (as defined in the terms of reference) in the light of a large number of publications in the recent literature and describes the present status and uncertainties associated with these issues. The findings of this study are summarized below.

Uncertainties Related to Observational Evidence and Data Analysis

It is now certain that mean global temperature has increased over the last 100 years; however, there is still uncertainty regarding the actual magnitude of this increase. For the recent 20-year period,

1978-1997, the mean surface temperature increase is estimated to be about 0.32°C, or about 0.1 6°C per decade; however, the impact of urbanization and land-use change, which has not been thoroughly established as yet, could contribute up to 0.1 °C (or more) towards the mean surface temperature increase. There is still some uncertainty about the temperature trend at the surface (as measured by surface network of stations) versus the temperature trend in the troposphere (as measured by earth-orbiting satellites and/or by radiosonde data). Based on the available satellite data (1979-1997), the mean temperature trend in the troposphere is now estimated to be about 0.10**C per decade; based on radiosonde data (1958-1998), the mean temperature trend for the troposphere (850-300 mb layer) is about 0.10**C per decade, the trend being stronger in the southern hemisphere than in the northern hemisphere. A causal and unequivocal link between the mean surface temperature increase and the anthropogenic greenhouse gas increase has not yet been established. The most probable cause of the mean surface temperature increase is considered to be a combination of intemally and externally forced natural variability and anthropogenic sources.

Significant uncertainty still exists relating the total (direct plus indirect) radiative forcing by anthropogenic aerosols (e.g., sulfate, black carbon, dust etc.). Recent studies suggest that the negative total radiative forcing by anthropogenic aerosols may offset the positive forcing by the greenhouse gases. Precipitation trends in different regions of the world do not present conclusive evidence about the intensification of the hydrologic cycle of the atmosphere-ocean system. There is still uncertainty relating trends in storm (tropical as well as extratropical) fi-equency in different parts of the world. Available climate data do not show any increasing trend in events (e.g., extreme precipitation, extreme drought, thunderstorms, winter bUzzards) in any part of the world.

Uncertainties Related to Climate Models and Other Issues

Global cloud climatology and its realistic simulation is perhaps the single largest uncertainty in climate models at present. An associated aspect is the modeling of the impact of anthropogenic aerosols on clouds and on the radiation budget. The presence of climate drift in simulations using coupled climate models and the associated flux adjustment problem represents a significant uncertainty regarding (model) climate sensitivity to natural and anthropogenic forcing of various origins. Sea ice cover is not well simulated in most climate models; fiulher, most climate models do not include adequate dynamics in simulating growth and decay of sea ice. This represents a significant weakness in climate models, many ofwhich project excessive depletion of ice cover fi*om the Arctic and Antarctic regions. There is considerable uncertainty in the simulation of mesoscale atmospheric and oceanic features and their impacts on fixture climate. Finally, the

Uncertainties in Greenhouse Gas Induced Climate Change iii strongest climatic signal, namely the El Nino southern oscillation (ENSO) phenomenon is still inadequately simulated by most climate models.

Possible Areas for Further Work

In addition to a number of broad areas of research, several research activities on a regional and

provincial scale are identified. In particular, it is suggested that a thorough analysis be made of precipitation and temperature trends over westem Canada to determine spatial and temporal

variability of these trends. It is further suggested that data on extreme events like blizzards, heat waves, intense precipitation, and floods, droughts, etc. be collected and carefully analyzed to establish trends (if any) in these extreme events

IV Uncertainties in Greenhouse Gas Induced Climate Change TABLE OF CONTENTS

Page

ABSTRACT iii

Uncertainties Related to Observational Evidence and Data Analysis iii

Uncertainties Related to Climate Models and Other Issues iii Recommendations for Further Work iv

LIST OF FIGURES vi

ACKNOWLEDGEMENTS vii

1.0 INTRODUCTION 1

1 . 1 Historical Notes and Overview 1

1 .2 Global Warming as Evidenced by Mean Temperature Increase 3

1.3 References , 5

2.0 ASSESSMENT OF UNCERTAINTIES IN OBSERVATIONS 7 2.1 Introduction 7 2.2 Cooling in Extratropical North Atlantic Ocean 7 2.3 Satellite and Land-Based Temperature Records 10

2.4 Solar Lradiance Change Affecting Last 1 00 Years 12 2.5 Change in Intensity or Frequency of Extreme Weather Events 13

2.6 Mean Temperature Increase and its Link to Greenhouse Gas Increase 19 2-7 References 21

3.0 ASSESSMENT OF UNCERTAINTIES RELATED TO CLIMATE MODELS AND OTHER ISSUES 28 3.1 Introduction 28 3.2 Specific Modeling Uncertainties 29 3.3 Other Issues 40 3.4 References 41

4.0 SUMMARY OF FINDINGS AND PRIORITY AREAS OF RESEARCH 46

4. 1 Summary of Important findings 46 4.2 Priority Areas of Research 48 , 4.3 References 49

Uncertainties in Greenhouse Gas Induced Climate Change V 1

LIST OF FIGURES

Page

Figure 1 . 1 Direct atmospheric measurements of CO2 concentrations at Mauna Loa and South Pole since about 1956 2

Figure 1 .2a Combined land and sea surface temperatures 1 869- 1 990 relative to the 1961-1990 mean 4 Figure 1.2b Hemispheric and global temperature averages 1856-1998 relative to 1961-1990 4 Figure 2.1a Changes of annual land-surface air temperature and sea surface temperature 8 Figure 2.1b Winter surface temperature trends 1961-1990 9 Figure 2.2a Interannual variation of Atlantic hxirricanes and hurricane days 14 Figure 2.2b Intense Atlantic hurricanes and strongest Atlantic hurricane wind speeds 15 Figure 2.3a Total number of tropical cyclones and number of intense systems in Australian region 1969-1995 16 Figure 2.3b Annual number of tropical storms and typhoons over western North Pacific.... 16 Figure 2.4 Changes in annual global-mean temperature and over the past 138 years relative to the 1961-1990 average 21 Figure 3.1a Global energy balance for annual mean conditions 3 Figure 3.1b Global average clear-sky radiation budget and average cloudy fluxes and cloud radiative forcing from ERBE data 32 Figure 3.2 Estimates of globally and annually averaged anthropogenic radiative forcing due to changes in concentrations of greenhouse gases and aerosols since pre-industrial times and due to natural changes in solar output since 1850 33 Figure 3.3a Simulated global annual mean warming from 1880 in two simulations with GHG forcing only and two simulations that include both GHG and direct sulfate aerosol forcing 35 Figure 3.3b Evolution of changes in annual global mean surface air temperature compared to the mean for observations 35 Figure 3.4 Evolution of decadal mean Arctic sea ice area for three experiments 38

VI Uncertainties in Greenhouse Gas Induced Climate Change ACKNOWLEDGEMENTS

I am grateful to Raymond Wong of the Environmental Sciences Division of Alberta Environment (Edmonton), for initiating this contract work and for making necessary arrangements to present the results of this report in a couple of seminars, one at Alberta Environment and one at the University of Alberta in Edmonton.

Several (e-mail) communications with a number of scientists have provided valuable input in the preparation of this report. I wish to gratefully acknowledge useful input and receipt of reprints and reports from the following scientists: Dr. John Christy, University of Alabama, Huntsville, USA; Dr. James Hansen, NASA/GISS, New York, USA; Dr. James Hurrell, NCAR (Boulder, USA) and Hadley Centre for Climate Prediction and Research, UK; Dr. David Karoly, Monash University, Clayton, Australia; Dr. Gerald Meehl, NCAR, Boulder, USA; Dr. J.F.B. Mitchell, Hadley Centre for Climate Prediction and Research, UK; Dr. C. Prabhakara, NASA/GSFC, Greenbelt, USA; Dr. Erich Roeckner, Max-Planck Institut fur Meteorologie, Hamburg, Germany; and Dr. Francis Zwiers, Canadian Climate Centre, Victoria, Canada.

During my attendance at the lUGG (Intemational Union of Geodesy and Geophysics) 1999 General

Assembly in Birmingham UK, I held extensive discussions with several scientists who made important presentations at the general assembly. I wish to acknowledge, in particular, fruitful discussions with Dr. Phil Jones (University of East Anglia, Norwich, UK), Dr. Joyce Penner (University of Michigan, Ann Arbor, USA), Dr. V. Ramanathan (University of CaHfomia, San Diego, USA), and Dr. Vin Saxena (North Carolina State University, Ralei^, USA).

The excellent facilities of Environment Canada's Library in Downsview, Ontario are gratefully acknowledged. Special thanks are due to Ms. Roberta McCarthy and Maria Latyszewskyj for providing valuable help in the literature search and for making special arrangements to provide access to the library facilities outside of normal working hours.

Finally, I wish to express my gratitude to my wife Shalan for her help on the home computer during the preparation of the report.

Uncertainties in Greenhouse Gas Induced Climate Change vii Uncertainties in Greenhouse Gas Induced Climate Change 1.0 INTRODUCTION

1 . 1 Historical Notes and Overview

The concept that the earth's atmosphere acts somewhat hke the glass of a greenhouse, letting through the sunlight (the shortwave light rays) while retaining a portion of the longwave radiation from the earth's surface was first introduced by French mathematician Joseph Fourier (1827). The idea was further elaborated by Pouillet (1838), Tyndall (1865), and Langley (1890). Tyndall was the first person to propose the "greenhouse" wanning hypothesis while Langley calculated the absorption coefficient for carbon dioxide and water vapour, while measuring the intensity of radiation from the moon.

The first estimates as to how changes in global concentration of "carbonic acid" (a primary greenhouse gas, now more commonly referred to as carbon dioxide) might affect mean global surface temperature were made by Swedish chemist Svante Arrhenius more than 100 years ago. Arrhenius' work was motivated by a desire to explain the temperature variation of the earth's surface during the quaternary glaciation cycles; accordingly he showed that an increase in the atmospheric concentration of CO2 by a factor of two would lead to a heating of the earth's temperature by 5 to 6°C (Arrhenius, 1896). He ftirther addressed the question whether such changes in the atmospheric CO2 concentration were plausible and made some calculations of the human impact on climate that were published in a couple of papers a few years after his landmark paper of

1 896. Arrhenius was influenced by the work of his geologist colleague Arvid Hogbom on the processes and possible changes in the global carbon cycle.

Arrhenius' work was followed by the work of the American geologist T. C. Chamberlin who supported the principal conclusion of Arrhenius and fiirther elaborated on how concentrations of atmospheric CO2 might change with time. Chamberlin's emphasis was on the processes of CO2 emissions into the atmosphere from volcanic eruptions, the absorption and outgassing of CO2 into and out of the world's oceans, and the role of rock formation and weathering in controlling terrestrial carbon reservoirs (Chamberlin, 1899). The studies of Arrhenius and Chamberlin did not receive much support among the atmospheric science community of that time, since the general consensus was that the absorption of longwave radiation (emanating from the earth) by water vapour was so strong that the absorption by carbon dioxide was negligible. This consensus was held in view of the fact that water vapour is the strongest absorber in the spectral region of 6 to 8 />tm (10*

^m) bands, while in the 1 5 to 20 fim bands, the absorption by water vapour and carbon dioxide overlapped so much that it was difficult to separate the two effects solely on the basis of surface observations (Ramanathan and Vogehnann, 1997).

In 1938, British engineer G.S. Callendar reintroduced Arrhenius' work and possible impacts of CO2 on the atmosphere and the earth's climate. On the basis of much improved analyses of absorption spectra for different gases, Callendar (1938, 1940) demonstrated that CO2 does indeed have important absorption bands outside of those dominated by water vapour and that increased CO2 concentration could have significant global effects on the surface temperature of the earth. Callendar suggested that the increase in concentrations of atmospheric CO2 may account for the observed slight rise of average temperature in northem latitudes during recent decades. He also speculated for the first time that humans could have significant influence on the atmospheric CO2

Uncertainties in Greenhouse Gas Induced Climate Change 1 — 1

concentrations, but estimated that it would take several centuries of continued industrial emissions to achieve a doubling of concentration. Callendar's work rekindled interest among atmospheric scientists about the role of CO2 on past as well as future climates, leading to the publication of several papers linking CO2 concentration with atmospheric temperature increase. A paper by Plass (1956) suggested that a 10% increase in atmospheric carbon dioxide would increase the average temperature by 0.36°C. In an important paper in 1957, the well-known American geophysicist Roger Revelle proposed that "humans are carrying out a large-scale geophysical experiment through world-wide industrial activity that could lead to a buildup of CO2 greater than the rate of CO2 production from volcanoes" (Revelle and Suess, 1957). Revelle was instrumental in establishing the first station for long-term monitoring of atmospheric CO2 at Mauna Loa (Hawaii) and in launching an accelerated intemational research program on the potential human influence on the climate system (Hengeveld, 1996). Another CO2 monitoring station was later established at the south pole and several global CO2 monitoring networks were established during the 1970s and 1980s. These networks have established the steady buildup of atmospheric CO2 as can be seen in Figure 1.1, which shows the well-publicized increase in concentration of CO2 at Mauna Loa and at South Pole during the last 40 years. The two curves (Mauna Loa and South Pole) clearly show seasonal cycles with opposite phases. The much larger amplitude in the Mauna Loa data is attributed to the larger land and vegetation areas of the northem hemisphere.

1 ; 1 1 1 1 1 . . 1 1 1 1 > 1 1 1 1 1 1 1 I 1 1 ? I I I 1 t 1 1 t t 1 1 1 1 t 1

Mill.

I -- tCOs South Pole - - —

ICO2 Mauna Loa, Hawaii 1

340 1 1 1 1 1 1 1 1 1 330 1 1 1 1 1 1 1 1 320 1 1 1 1 1 1 1 1 1 1)

1 [ 1 1 J I 1 r- II 1 I 1 I 1 1 1 1 1 1 1 1 1 I I ! 1 1 1 1 1 1 1 i — i— 1 1 i —— ———— \ 1960 1970 1980 1990

Figure 1.1 Direct atmospheric measurements of CO2 concentrations at Mauna Loa and South Pole since about 1956. Concentration is in units ofppmv (parts per million, volume basis). Note that seasonal oscillations in the interannual variation at South Pole are in opposite phase to those at Mauna Loa and considerably smaller.

Uncertainties in Greenhouse Gas Induced Climate Change The Intergovemmental Panel on Climate Change (IPCC) was jointly established by the World Meteorological Organization (WMO) and the United Nations Environment Program (UNEP) in 1988. The IPCC is now recognized as the prime source of scientific and technical information on climate

change and its environmental and socioeconomic impacts. The first assessment report on climate change was completed by the IPCC in 1990. In its second assessment report on climate change (IPCC, 1996), the IPCC identified the increasing concentration of CO2 and its possible link with the increasing mean surface temperature of the earth with the cautious statement that the balance of evidence suggests a discernible human influence on global climate. The IPCC has also recognized other greenhouse gases, namely, methane (CH4) and nitrous oxide (N2O). The atmospheric concentrations of these two gases have also grown significantly since pre-industrial times (i.e., since about 1750) by about 145% (methane) and by about 15% (nitrous oxide) based on 1992 values. However, the direct radiative forcing by these two gases is significantly smaller (0.47 W-m'"^ and 0.14

W-m'^) as compared to 1 .56 W-m""^ for carbon dioxide; consequently, the present emphasis is on the global carbon cycle and the fiiture growth of carbon dioxide in the earth-atmosphere system over the next 50 to 100 years. The sharply increasing concentration of methane and its consequent radiative impact on the fiiture climate has also been emphasized by the IPCC.

1.2 Global Warming as Evidenced by Mean Temperature Increase

The key question at present is whether the atmosphere-earth system has warmed in recent times and, if so, by how much? The answer to the first part of the question is definitely affirmative, while the answer to the second part of the question is still indeterminate because of a number of uncertainties that will be discussed later.

Determination of the mean temperature of the combined land/ocean surface of the earth for a given year is perhaps the most difficult and most elusive task. Several land and water surface processes (e.g., urbanization, land-use change, surface and subsurface water circulations in the oceans) can influence and distort mean surface temperature calculations. Nevertheless, painstaking efforts by Phil Jones of the Climate Research Unit, University of East Anglia (UK) and his coworkers, through carefiil analysis of a vast amount of land-ocean surface data, have helped elucidate mean temperature variations over the last 140 years. These mean temperature calculations published in a series of papers

(Jones et al.; 1986a, 1986b, 1999) have now become a benchmark and provide a firm basis for the scientific discourse on climate change. The mean surface temperature changes for the period 1861- 1995 as given in IPCC (1996) are shown in Figure 1.2a, while Figure 1.2b shows the mean temperature variations over the northern and southem hemispheres together with global variations based on a more recent paper by Jones et al. (1999). These temperature curves show mean temperature departures (**C) relative to 1961-1990. Some of the important variations in the mean temperature curves may be noted as follows:

1.2.1 The earth's mean temperature shows a steady increase of ahnost 0.5°C fi-om about 1 91 0 to about 1940. Mean temperature shows a decline of about 0.2*'C fi-om about 1950 to 1975. From 1980 through 1995, the mean temperature of the earth shows a relatively steep

increase of about 0.25°C; the recent study by Jones et al. indicates a warming of about 0.32°C over the period 1978-1997.

1.2.2 The temperature change over the southem hemisphere is sHghtly larger than over the northem

hemisphere, while the global temperature change over the entire period 1861-1997 is about 0.57°C.

Uncertainties in Greenhouse Gas Induced Climate Change 3 S -0-6

O -0.8 1860 1880 1900 1 920 1 940 I960 1 980 2000 Year

1860 1880 1900 1920 1940 1960 1980 2OO0

Figure 1.2a (top) Combined land/sea surface temperatures 1869-1990 relative to the 1961-1990 mean. The solid curve represents smoothing of the annual values, shown by the bars, to suppress sub-decadal time-scale variations. The dashed smooth curve is the corresponding result from the IPCC (1992) (from IPCC, 1996).

Figure 1.2b (bottom) Hemispheric and global temperature averages 1856-1998 relative to 1961-1990. The smooth line on the time series highlights variation on the decadal time scale (10-

year Gaussian filter) (from Jones et al., 1999).

4 Uncertainties in Greenhouse Gas Induced Climate Change . .

These temperature changes, although numerically small, are nevertheless significant from the standpoint of radiative equilibrium and warrant a close examination of their possible influence on the earth's climate. The questions facing the scientific community today are:

a. Are these temperature changes caused by anthropogenic greenhouse gases, natural variations in climate, or a combination of the two? b. Are these temperature changes inducing climate change in terms of increased variability or increased incidence of extreme weather events or both?

This report provides an assessment of uncertainties in greenhouse gas induced climate change based on recent studies and the present state of our knowledge. In the next two chapters, several key issues facing the science of climate change (as defined in the terms of reference of the project) are examined and a detailed analysis and assessment of these issues is presented. The final chapter presents a brief summary and identifies priority areas for fiirther research.

1.3 References

Arrhenius, S., 1 896: On the influence of carbonic acid in the air upon the temperature of the ground. Philosophical magazine and Journal ofScience, 41 : 237-276.

Callendar, G.S., 1938: The artificial production of carbon dioxide and its influence on temperature, g. J. Royal Met. Soc, 64, 223-

Callendar, G.S., 1940: Variation in the amount of carbon dioxide in different air currents. Q.J. Royal Met Soc, 66, 395-

Chamberlin, T.C., 1899: An attempt to frame a working hypothesis of the cause of glacial periods on an atmospheric basis. J. ofGeology, 7, pp.75 1-787.

Fourier, J., 1827: 569-604.moire sur les temperatures du globe terrestre et des espaces planetaires,

Mem. de I'Academie Royale des Sciences de Vlinstut de France, 1, 569-604.

Hengeveld, H., 1996: Arrehenius' hypothesis: One century later. CMOS (Can. Meteor, and Oceanogr. Soc.) Bulletin, 24, no.3, 63-65.

IPCC, 1996: Climate Change 1995: The Science of Climate Change. Houghton, J.T., Meira Filho, E.G., Callendar, B.A., Harris, N., Kattenberg, A. & Maskell, K. (eds.), Cambridge Univ. Press, 570p.

Jones, P. D., S.C.B. Raper, R.S.Bradley, H.F.Diaz, P.M. Kelly and T.M. L. Wigley, 1986a: Northem hemisphere surface air temperature variations: 1851-1984. /. Climate & Applied Meteor., 25, 161-179.

Jones, P. D., S.C.B. Raper and T.M. L. Wigley, 1986b: Southem Hemisphere surface air temperature variations: 1851-1984. J. Clim & AppHed Meteor. 25, 1213-1230.

Uncertainties in Greenhouse Gas Induced Climate Change 5 Jones, P. D., M. New, D. E. Parker, S. Martin and LG. Rigor, 1999: Surface air temperature and its change over the past 150 years. Reviews of Geophysics, 37, no.2, 173-199.

Langley, S.P., 1890: The temperature of the moon. Mem. of the Nat. Acad, of Sci., iv: 193.

Plass, G. N., 1956: The carbon dioxide theory of chmatic change. Tellus, 8, 140-154.

Poullet, C. 1837: Memoire sur la chaleur solaire, sur les pouvoirs rayonnonants et absorbants d'lair atmospherique, et sur la temperature de I'espace. Compte rendus hebdomadaires des

seances de 1, Academic des Sciences 7, 24-65.

Ramanathan, V. and A. M. Vogelmann, 1997: Greenhouse effect, atmospheric solar absorption and the Earth's radiation budget: from the Arrhenius-Langley era to the 1990s. AMBIO, Vol. xxxvi, no.l, 38-46.

Revelle, R. and H.E. Suess, 1957: Carbon dioxide exchange between atmosphere and ocean and the question of an increase of atmospheric C02 during the past decades. Tellus, 9, 18-27.

Tyndall, J., 1 861 : On the absorption and radiation of heat by gases and vapours and the physical connection of radiation, absorption and conduction. Phil. Magazine, 22, 169-194.

6 Uncertainties in Greenhouse Gas Induced Climate Change 2.0 ASSESSMENT OF UNCERTAINTIES IN OBSERVATIONS

2.1 Introduction

Since the publication of the Chmate Change 1995 report (IPCC, 1996), a large number of studies based on observed data have attempted to verify some of the conclusions regarding present and the future climate by analyzing the available observed data and by examining model simulations of future climate. The thrust of most of the observational studies is to assess if recent weather events and atmospheric data (e.g., precipitation, temperature, tropical and subtropical storms, etc.) indicate any variabiHty that can be associated with possible climate change that, in turn, can be linked to global warming. Many of these studies have identified areas where observations do not support projections of future climate made by computer models. There have also been studies reported in recent Hterature questioning the magnitude of the warming as well as the validity of the mean surface temperature calculations and associated impacts on future climate. (Lindzen, 1990; Michaels and Stookebury, 1992). The global warming debate began in the early 1980s, well before the publication of the IPCC Climate Change 1995 report. The IPCC report brought the debate into sharper focus and has led to publication of reports by a group of scientists (e.g., ESEF Report,

1996) or commentaries (e.g., Singer et al., 1997, Kondratyev, 1997) and books hke "Hot Talk, Cold Science" by Singer (1997). These publications have prompted many scientists and researchers to take a closer look at the available data to determine if the perceived climate variability should be attributed to warming due to greenhouse gas increases, natural climate variability, or a combination ofboth.

In this chapter, several issues related to observational uncertainties (as defined in the terms of

reference) are analyzed in detail. The various issues are first stated in italics and then analyzed in the light of recent studies and the present state of our knowledge:

2.2 Cooling in Extratropical North Atlantic Ocean

Cooling has been observed in the North Atlantic andparts ofeastern Canada for the last 25-30 years, while most coupled models do not simulate this cooling very well.

The spatial distribution of mean surface temperature relative to the 1961-1990 period shows considerable variation in the magnitude of warming as well as cooling in different geographical regions of the earth (see Figure 2.1a,b). Three significant areas of cooling are: north central Pacific, northwest Atlantic, eastern Canada, and parts of central and eastern Europe. The most notable among these regions is the northwest Atlantic together with eastern Canada, which show significant cooling in the winter months (Figure 2.1b), and another region extending fi-om southeastern USA to the Caribbean Sea and thence southward (Figure 2.1a). Many of the earlier

climate models (e.g., Washington and Meehl, 1989; Manabe et al., 1991,1992) do not reproduce this cooling in the northwest Atlantic very well (in fact some of the climate models show a warming of up to 2 to 4°C in this region). A recent model developed at the Max Planck Institute

for Meteorology, Hamburg (Germany) by Roeckner et al., (1995, 1998) shows a large area of cooling over the Canadian sector of the North Atlantic, in agreement with the observed temperature

Uncertainties in Greenhouse Gas Induced Climate Change 7 8 Uncertainties in Greenhouse Gas Induced Climate Change Uncertainties in Greenhouse Gas Induced Climate Change 9 distribution shown in Figure 2.1b. This improvement is attributed to improved simulation of the natural climatic variability by the model (see Bengtsson, 1997).

Further discussion on climate model performance and weaknesses will be deferred to the next

chapter; suffice it to say here that the cooling over the northwest Atlantic and over eastern

Canada is now attributed to the positive phase of the NAO (North Atlantic oscillation), according to recent studies by Hurrell (1995, 1998) and others. Another recent study by Shabbar et al., (1997) has developed an index called the BWA (Baffin Bay-west Atlantic) index to explain the cooling over eastern Canada; this BWA index is a reflection of the western structure of the NAO and is more closely related to the negative temperature anomalies over the northwest Atlantic. The NAO and other large-scale oscillations, for example, the recently identified Arctic oscillation (Thompson and Wallace, 1998; Kerr, 1999) are not yet well simulated by many climate models. The cooling over southeastern USA in the Atlantic seaboard states is now attributed to negative radiative forcing by sulfate aerosols (e.g., Saxena and Yu, 1998; Saxena and Menon, 1999). A large number of studies on the important topic of climate forcing by anthropogenic aerosols have been reported in the recent literature; a detailed discussion on this appears in the next chapter.

2.3 Satellite and Land-Based Temperature Records

Satellite records show a slight cooling instead ofwarming during their past 19 years ofexistence. Recent analysis shows that thefalling altitude ofthe satellites can be a reason for the apparent cooling trend. However, the results show great uncertainty about the impact ofthis effect on the satellite records. Land-based temperature records are also compared with satellite data. Uncertainties in land-based records include the effect ofurbanization and uneven distribution of stations. These have impacts on mean surface temperature calculations.

As mentioned in the first chapter, calculation of the mean surface temperature of the earth, free of sampling error associated with urbanization and other effects, is perhaps the most difficult task. The mean surface temperature calculation and its upward trend are being increasingly questioned because of strong disagreements with mean tropospheric temperature trends calculated using satellite-based radiometric data fi-om the past 19 years (1979-1997). Several important and painstaking studies, notably by Spencer and Christy (1990, 1992) and Christy et al. (1995, 1998) have merged MSU (microwave sounding unit) data fi^om several satellites to produce a mean daily temperature value for the troposphere and its trend over the years for which satellite data have been available. These studies show a small cooling trend for the mean tropospheric temperature, inconsistent with the warming trend for the mean surface temperature documented in several papers by Jones and his coworkers. This inconsistency in temperature trends has been analyzed by Hurrell and Trenberth (1996, 1998) and by Jones et al. (1997), among others. The difference in temperature trends has been attributed to a number of factors, such as merging of data fi-om different satellites, uncertainties regarding the level and character of natural differential variability between the surface and lower troposphere, uncertainties in differential response to anthropogenic and/or natural (e.g., volcanic) forcing, and errors due to satellite drift from original local-equator crossing time (LCT).

The recent paper by Christy et al. (1998) recalculates the mean temperature trend for the troposphere applying several corrections and using improved techniques for merging data from

10 Uncertainties in Greenhouse Gas Induced Climate Change different satellites. The revised calculations yield a smaller cooling trend (~ -0.046°C per decade) for the mean tropospheric temperature, still at odds with the surface temperature trend of +0.1 5°C per decade. Hansen et al. (1995) made a comprehensive analysis of the differences between surface and tropospheric temperature trends and have attributed part of the difference to stratospheric , which cools the troposphere more than the surface. Hansen et al. further suggest that urbanization would influence the mean surface temperature calculation by at least 0.1 °C over 100 years. The satellite-based MSU data have been reanalyzed by Prabhakara et al. (1998) and a mean temperature trend of about +0.1 TC per decade during the period 1980-1996 has been estimated for the troposphere; this trend is in better agreement with the surface temperature trend. Another recent study by Angell (1999) calculated the mean temperature trend for the troposphere (850-300 mb layer) using a 63-station radiosonde network and obtained a mean tropospheric temperature trend of +0.1 0°C per decade during the 41 -year period (1958-1998). This value appears to be in very good agreement with the temperature trend calculated by Prabhakara et al. for the period 1980-1996.

Assessing the impacts of urbanization and land-use change on the mean surface temperature calculation is a challenging task. The classical studies of Mitchell (1961) and Oke (1973) suggest that the effect could be significant even for towns with populations of few thousand people. For large cities with populations of two million or more, Oke's study estimates the urbanization impact could increase the surface air temperature at urban sites by 10°C or more.

According to the IPCC (1990), the urban heat island (and its potential impact on surface air temperature) is probably the most serious source of systematic error in land surface climatological measurements. Several studies published in the last 15 years or so have attempted to assess the effects of urbanization on local and regional climate (e.g., Cayan and Douglas, 1984; Balling and Idso, 1989; Karl et al.,1988; Goodrich, 1992; Portman,1993). In the context of global warming there have been several studies, notably by Karl and Jones (1989) and Jones et al. (1990) to assess the impact of urbanization on the calculation of the mean surface temperature of the earth as a whole and to estimate a numerical value for it. According to Jones et al. (1990), the impact of urbanization on the mean surface temperature would be no more than 0.05**C over 100 years. This value appears to be too small in view of similar calculations for individual large cities or local regions. For example, Fujibe (1995) obtains an average rising trend of 2-5°C per 100 years in the rtiinimum temperature at several large cities in Japan, while Hingane (1996) estimates rising trends of 0.84°C and L39°C per 100 years in the mean surface temperatures calculated for Bombay and

Calcutta, two of the largest cities in hidia. It is now recognized that urbanization and changing laud use influence minimum temperature, which, in the last 100 years has risen faster than maximum temperature at most locations. The more rapid increase in minimum temperature (than maximum temperature) has led to a decrease in the diurnal temperature range (DTR) at most locations on a worldwide basis (Karl, 1993) and an increase in the mean surface temperature, estimated to be up to

0.25°C or more (Easterling et al. 1997), thus accounting for up to half of the global warming estimated by Jones and coworkers. In a more recent study (Gallo et al., 1999), over 1200 weather observing stations in USA were designated as either urban, sub-urban, or rural using a night-light index from a satellite-based device, and temperature trends at rural and urban sites were determined for the period 1950-1996. The study found that the decreasing trend in DTR (primarily associated with increases in minimum temperature) was smaller at rural stations than at urban stations by about 0.45°C per 100 years. This value, although not statistically significant, suggests that urbanization (e.g., increased humidity and cloud cover ) tends to increase the minimum temperature more than the maximum temperature, thus reducing DTR and boosting mean

Uncertainties in Greenhouse Gas Induced Climate Change 11 temperature. Another recent study (Hansen et al., 1999) makes a comprehensive analysis of surface temperature change and estimates the anthropogenic urban contribution to the global mean temperature curve for the past century to be about 0. TC.

In summary, the satellite-based MSU data, when suitably corrected for decay of satellite altitude (Wentz and Schabel, 1998) and other corrections as mentioned above show a small but definite warming trend in mean temperature of the troposphere over the last 19 years (1979-1997); the

actual magnitude of the warming trend is still a subject of study (C. Prabhakara, personal communication). The MSU data and the mean surface temperature data appear to be reconciling towards a small but a definite increasing trend throughout the tropospheric depths of the

atmosphere (Hansen et al., 1998).

2.4 Solar Irradiance Change Affecting Last 1 00 Years

The influence ofsolar activities may contribute to the warming observed in the ground-level temperature record. This may be significant in the post-1900period. What are the recentfindings on the solar contribution to observed warming?

The IPCC 1996 report did consider solar irradiance change over the last 100 years and its possible

impact on the global warming. However, it was concluded that the solar irradiance variations of the past century are likely to have been considerably smaller than the anthropogenic radiative forcing,

and hence its impact on climate was considered to be insignificant. Several recent studies (e.g.,

Hoyt and Schatten,1993; Lean et al., 1995a, 1995b; Solanki and Fhgge, 1998) have attempted to reconstruct the solar irradiance variations over the last 300 years or more and these reconstructed solar irradiance values have been studied in conjunction with variations in the earth's mean surface

temperature by Lean et al. (1995) and Lean and Rind (1998), among others. These studies now suggest that solar forcing can be significant and may have contributed to some of the recent observed global warming. According to Lean et al.(1995), solar forcing may have contributed about half of the observed 0.55°C surface warming since 1860 and about one third of the warming since 1970. In a recent comprehensive study. Lean and Rind (1998) correlated the reconstructed solar irradiance values with the northern hemisphere surface temperature anomalies and found a correlation of 0.86 over the preindustrial period from 1610 to 1800, implying a predominant solar influence on the preindustrial climate. Extending this correlation analysis to the industrial era from

about 1850 to about 1960, during which the solar activity has generally increased, it is possible to attribute up to half of the observed global warming to the increased solar output. Since about 1970, greenhouse gas radiative forcing has been significantly larger than solar irradiance-change forcing. Consequently the solar forcing signal in the global warming data cannot be easily isolated; nevertheless. Lean and Rind attribute up to one third of the recent warming to increasing solar irradiance output.

In summary, recent studies on reconstruction of solar irradiance and its possible impact on the earth's climate has added a new dimension to the global warming debate. Further, these studies strongly suggest that at least part of the observed global wanning can be explained due to increased solar irradiance during the last 100 years.

12 Uncertainties in Greenhouse Gas Induced Climate Change 2.5 Change in Intensity or Frequency of Extreme Weather Events

Extreme weather is expected to increase in frequency and intensity with global warming. Recent analyses ofextreme weather intensities andfrequencies show conflicting results. What are the majorfindingsfor drought, storms, andfloods? What do GCMs say about extreme weather? A related aspect is the relationship between El Nino and warming. Are there new developments regarding the linkages between El Nino and warmer climate?

There are several aspects of this issue that require careful assessment. The projections of future climate according to the IPCC 1996 report suggest

a. increased frequency and intensity of storms (tropical as well as extratropical), b. increase in the area of the globe affected by droughts and floods, and

c. intensification of the global hydrologic cycle (e.g., vigorous Indian/Asian monsoon).

Let us consider these points separately:

Storms: An analysis of tropical storms (hurricanes, tropical cyclones, typhoons, etc.) and their interannual variability indicates that these storms do not show any increasing trend in terms of frequency or intensity. The North Atlantic tropical storms (or hurricanes as they are popularly known) have been studied extensively by Gray (Colorado State University, USA) and his coworkers, and the results of their studies have been reported in a number of papers (e.g.. Gray,

1984a, 1984b; Gray et al., 1994; Landsea, 1993; Landsea et al., 1999). Based on these and other studies it is now well established that the Atlantic hurricane frequency and intensity are governed by several large-scale features, including the phase of the ENSO, Caribbean basin sea level pressure and 200-mb zonal winds, west Sahel (Africa) rainfall, and Atlantic sea surface temperature and its variation. According to Landsea et al. (1999), the Atlantic hurricane variability can be characterized as lacking strong linear trends, but is comprised of robust multidecadal variations. Li particular, intense hurricanes—those reaching sustained winds of 50 m/s"' were very common in the 1940s through the 1960s, and were much reduced in occurrence from the 1970s through the early 1990s. The years 1995 and 1996 showed a dramatic increase in hurricane activity, most likely because of the La Nina (cold) phase of the ENSO. The 1997 hurricane season was remarkably quiet (Rappaport, 1999) with only six tropical storms formed in the Atlantic basin, compared to the long- term average of ten. For ready reference, interannual variation of hurricanes, hurricane days, strongest Atlantic hurricanes, and hurricane wind speed are shown in Figure 2.2a,b. These figures underscore the conclusion of Landsea et al., as stated earher. In the central Pacific, the total number of tropical cyclones and typhoons shows a decreasing trend from 1959 till the mid-1970s, followed by an increasing trend during recent years (Figure 2.3a). In the Australian region, there was a distinct reduction in the total number of tropical cyclones in recent years while a small increase is seen in the number of intense cyclones (Figure 2.3b). No study has been reported so far on the interannual variability of tropical cyclones in the Indian Ocean area, while the tropical cyclones

(and depressions) in the Bay of Bengal show no discernible trend. In simmiary, there is no increasing trend in the tropical cyclone frequency or intensity in any of the three ocean basins, and the interannual variability appears to be primarily governed by the ENSO phase and other large-scale features of the atmosphere-ocean system.

Uncertainties in Greenhouse Gas Induced Climate Change 13 Atlantic Hurricanes 1944-1996

12 1 -J

10 -

Atlantic Humcane Days 1944-1996 60 1— 1 —

Figure 2.2a Interaimual variation of Atlantic hurricanes and hmricane days (from Landsea et al, 1999)

14 Uncertainties in Greenhouse Gas Induced Climate Change " r

Atlantic Intense Humcanes 1944-1996

o Is ^4

J 22H

1 1 H

r r i 0 r r r r r r r r r r r r r r r r I rrrrri rrrruprrrrrr rrr rpr rrrrrrr> 1944 1949 19S4 1959 1984 1969 1974 1979 1984 1989 1994 Year

Strongest Atlantic Humcane WIndspeed 1944-1996 90

80

^ 70 - r S80 H r o 50 -1 340

|»30 -I ^20

10 -

r r P r y.£r^' r r r r r r r r rrrrrrrrrr rr pt r * " "r* tt rrrrrrrrrrprrr . 1944 1948 1954 1959 1984 1989 1974 1979 1984 1989 1994 Year

Figure 2.2b Intense Atlantic hurricanes and strongest Atlantic hurricane wind speeds (from LandseaetaL, 1999)

I

Uncertainties in Greenhouse Gas Induced Climate Change 15 — —

' 0 — — •—. . V 1969 1972 1975 1978 1981 1984 1987 1990 1993 Season

Figure 2.3a (top) Total number of tropical cyclones (solid lines) and number of intense systems (dotted lines) in Australian region (105-160°E) 1969-1995

Figure 2.3b (bottom) Annual number of tropical storms and typhoons over western North Pacific. Two second-order polynomial fits to the series are also shown (from Shan and Shi, 1996).

Uncertainties in Greenhouse Gas Induced Climate Change The interaimual variability of extratropical storms appears to show an increasing trend but the results are conflicting. An earlier study by von Storch (1993) concludes that there was no observational evidence of systematic change in either the frequency or the severity of extratropical storms in the North Atlantic over the past 100 years. Lambert (1996) analyzed intense cyclones in the North Pacific and North Atlantic and found evidence to support an increase in intensity during the past several decades. A more recent study by Bijl et al. (1999) analyzes long-term observational sea-level data sets in the context of a possible change in storminess over northwest Europe. The study concludes that, for sea-level stations in the southern North Sea, there has been a tendency towards a small weakening of the storm climate over the past 100 years. For the more northern stations, there was no indication of a weakening of the storm climate; a small but insignificant increasing trend was indicated, hi another recent study

(Hirsch et al., 1999), an analysis of winter storms climatology (based on 1951-1997 data) reveals a decreasing trend in the frequency of winter storms for the east coast of North America. Thus there is no increasing trend emerging out of the studies reported so far.

Floods and droughts: A study by Balling Jr. (1996) analyzes over 100 years of excellent drought records in the USA during a period of substantial buildup of greenhouse gases and concludes that there is no evidence of a drying trend in any climatological region of the USA. A similar study on drought in Canada (Maybank et al. 1995) finds that the Pahner Drought Severity Lidex (PDSI, a commonly used drought index) shows no significant trend during the last 80 years over the

Canadian prairies, the most drought-prone region in Canada. A more recent study by Dai et al. (1998) calculates the PDSI on a global scale using gridded monthly air temperature and precipitation data over the period 1900-1995. The study finds large multi-year to decadal variations in the percentage areas in severe drought (PDSI < -3.0) and severe moisture surplus (PDSI > +3.0) over many land areas, but no significant secular trends. Since the late 1970s, there have been some increases in the combined percentage areas in severe drought and severe moisture surpluses, however, these recent changes are attributed to the shift in the ENSO phase towards a warmer phase, which has resulted in more frequent El Nino events during the last 25 years or so. Another recent study (Plummer et al., 1999) on changes in climate extremes over the Australian region reveals that the percentage area of AustraUa experiencing extreme wet conditions has increased sHghtly while the area of extreme dryness has decreased slightly since 1910. A comprehensive study on interannual variability of seasonal rainfall over various subdivisions of India

(Parthasarathy et al., 1995) based on over 125 years of excellent data sets shows no increasing or decreasing trend in rainfall amounts over any of the subdivisions, including the desert region of Rajasthan in northwest India. The studies reported so far do not reveal any increase or decrease in the total areas of floods and droughts worldwide.

Intensification of global hydrologic cycle: This is perhaps the most contentious aspect of climate change as a result of global warming. Several recent studies have made detailed analyses of precipitation data on local, regional, and continental scales attempts to find an increasing trend in either the total precipitation or in the frequency of "heavy precipitation" events; the results appear to be conflicting at this time. Among the often-quoted studies is the painstaking analysis by Karl and Knight (1998) showing an increasing trend in the percentage contribution of the upper 10th percentile of daily precipitation events to the total annual precipitation, area-averaged across the conterminous USA. However, such an analysis for Canada does not show any definite increasing trend. According to Mekis and Hogg (1999), the fraction of annual precipitation falling in heavy

Uncertainties in Greenhouse Gas Induced Climate Change 17 —

(90th percentile) events has decreased by over 4% since 1910 in southern regions of Canada. At more northern stations (north of 55°N), this jfraction has increased by 5%, but this trend is based on a shorter data set (1940-1995). Another recent study on precipitation trends on the Canadian

prairies (Akinremi et al., 1999) finds an increase in the number of precipitation events on the Canadian prairies, but this increase is mainly due to increase in the number of low-intensity events (precipitation amounts in the 0.5 to 5.0 mm range); the study further concludes that "precipitation events are not getting more intense on the Canadian prairies". Outside of North America, precipitation trend studies show inconclusive results. Over the Australian region, heavy rainfall has

increased in some areas (based on data period 1910-1995), but this increase is not significant. The

study by Zhai et al. (1999) on climate extremes in China concludes that, for China as a whole, there were no significant trends in much greater than normal amounts of annual precipitation either in 1- day or 3-day maximum rainfall or in rainfalls derived firom daily precipitation amounts greater than 10, 50, or 100 mm. For the monsoonal climate of India, Rakhecha and Soman (1994) analyzed

annual extreme rainfall time series of 1 to 3 days duration at 316 rainfall stations, well-distributed

over India. Only 14 of the 316 stations showed significant trends in all the three durations stations on the west coast of peninsular India (north of 12°N) showing an increasing trend while stations in the southern peninsula showing a decreasing trend. It may be noted that several ofthe highest 3-day rainfall amounts were more than 500 mm, well over annualprecipitation amounts at

many Canadian stations! Another comprehensive study by Groisman et al., (1999) applies a statistical model of daily precipitation based on the gamma distribution and obtains a 20% increase in the probability of summer daily precipitation amounts over 25.4 mm in a few northern European countries (e.g., Russia, Norway, and Poland); however, this increasing trend has not been reported elsewhere in Europe. Further, the precipitation amounts for some of the northern European countries are much too small to be of any real significance; for example, the increased number of days with summer precipitation in Norway is based on threshold values of 0.2 mm and 1 .0 mm per day, much too small compared to daily values in the monsoonal climate of south Asia. No studies have been reported so far on increasing trends in "heavy" precipitation over southeast Asia. A few studies on interannual variability of rainfall over southeast Asia (Khandekar et al., 1999; Kriplani and Kulkami,1997) indicate that seasonal rainfall over Indonesia/Malaysia and vicinity is primarily govemed by the ENSO phase and does not reveal any increasing tendency in recent years. In summary, the hydrologic cycle, while showing variability over different regions, does not reveal any intensification on a global or continental scale. Over the tropical monsoon regions where the hydrologic cycle is the strongest, the annual and seasonal rainfalls appear to be primarily govemed by ENSO and other local and regional features.

A brief mention about the GCM (General Circulation Model) simulation of extreme climate in a warmer world: Typically GCM simulations obtain changes in mean temperature, daily precipitation, and wind speed in a doubled CO2 environment and examine these changes (against a "control" climate) in terms of increased/decreased probability of certain threshold values of these weather parameters being exceeded. Recent GCM simulation studies (e.g., Jones et al., 1997; Hennessy et al., 1997; Zwiers and ICharin, 1998) indicate an increase in global average precipitation by 10% or more and increased fi-equency of heavy precipitation, especially over the south Asian monsoon zone. According to Zwiers and Kharin (1998), strongest increase in precipitation probability (precipitation over 50 mm/day) is found over northwest and central India, due to intensification of the Asian summer monsoon under a doubled CO2 atmosphere. Recent observational studies of the Indian/Asian monsoon (reported earlier) do not support the findings of

18 Uncertainties in Greenhouse Gas Induced Climate Change Zwiers and Kharin. Further, the global wanning projected by Zwiers and Kharin is much larger than expected on the basis of recent and improved model projections (this will be discussed in the next chapter). Another study (Bhaskaran and Mitchell, 1998) simulates monsoon precipitation change over the Indian subcontinent in a doubled CO2 environment and obtains (with both greenhouse gas and aerosol forcing included) a 12% decrease in precipitation in the western region of India, and an increase of 20% in the eastern and northeastern regions of India. Once again, there is no observational support for these simulations (see Rupa Kumar et al., 1992).

Linkages between EI Nino and warmer climate: The tendency for more frequent El Nino events and fewer La Nina events since the late 1970s has prompted a question whether El Nino events are being influenced by the warmer climate due to increasing greenhouse gas concentration. The El Nino event of 1990-1995 is considered to be the longest on record and this has prompted Trenberth and Hoar (1996) to link El Nino events to decadal changes throughout the Pacific basin. Using a simple ARMA (autoregressive moving average) model to a measure of the southern oscillation index, Trenberth and Hoar obtain the probability of occurrence of a prolonged El Nifio (like the 1990-95 El Nino) to be extremely small; this has led to the suggestion that such a prolonged ENSO and the associated changes may have been partly caused by the observed increase in greenhouse gas concentrations. The suggestion of Trenberth and Hoar was closely examined by Rajgopalan et al, (1997) who used a nonparametric statistical technique to analyze 114 years of continuous data on Darwin sea-level pressure (a good measure of the southern oscillation index) and demonstrated that the probability of occurrence of an El Nifio event like that in 1990-95 is much higher than that suggested by Trenberth and Hoar. According to Rajgopal et al., such a prolonged (and unusual) El Nino event could occur through natural climatic variation. Nevertheless, the suggestion that El Nifio events could become more frequent in a warmer world is intriguing and has caught the attention of many atmospheric scientists. It is not clear at this time whether El Nifio events are becoming more frequent. Following the intense El Nifio of 1997-98 (from March 1997 to about May 1998), the SST (sea surface temperature) anomalies in the equatorial eastern and central Pacific have remained negative (La Nifia phase) for the past 16 months or so, and may remain negative until March 2000 according to recent projections (Climate Diagnostics Bulletin, 1999). More studies are needed to investigate interannual and decadal variability in the Pacific basin before we can link El Nifio events to warmer climate.

2.6 Mean Temperature Increase and its Link to Greenhouse Gas Increase

The relationship between mean global temperature increase and greenhouse gas increase has not been unequivocally established. What are the uncertainties in the relationship?

The increase in greenhouse gas concentration in the atmosphere over the last 100 years has been well established as discussed earlier. The increased concentration has led to an increase in the mean radiative forcing of 2.45 W-m'^ and this has also been well established by the IPCC (1996) and several other studies. What is not unequivocally established is the direct link between the mean surface temperature increase and the increase in the radiative forcing by the greenhouse gas increase. There are several other processes in the complex earth-atmosphere-ocean system that make this unequivocal link difficult to establish. Some of these processes have already been mentioned and discussed, for example, the role of large-scale atmospheric circulation features

Uncertainties in Greenhouse Gas Induced Climate Change 19 like ENSO and NAO in influencing the mean surface temperature patterns over different areas of landmasses and oceans. Considering the radiative equilibrium of the earth-atmosphere system, the uncertainty stems from several factors, including the role of increasing solar irradiance, the direct and indirect impact of aerosols and carbonaceous particles, and the role of cloud microphysics and fractional cloudiness on radiation balance and ultimately on mean surface temperature. A recent paper by Solanki and Fligge (1998) estimates an increase of 2-5 W-m'^ in solar irradiance over the past 100 years and suggests a possible link between solar irradiance increase and mean surface temperature. Further, direct and indirect negative radiative forcing by sulfate aerosols are now estimated to be comparable in magnitude to positive forcing by increased greenhouse gas concentrations. Consequently, a simple and direct link between the increasing concentration of CO2 (the most important of the greenhouse gases) and mean surface temperature increase is not easy to establish. This point is illustrated in Figure 2.4, which shows estimated changes in annual global-mean temperature together with changes in CO2 concentration over the past 138 years relative to thel 961 -1990 average. The variation in CO2 concentration is shown by a steadily increasing (smoothed) curve, whereas the mean temperature curve shows irregular variations with several time periods (each several years long) during which the mean global temperature has increased or decreased, either steadily or abruptly. Such increases and decreases of the mean temperature cannot be easily explained by a more-or-less linear increase in CO2 concentration. In particular, the steep increase in the mean temperature from 1910 to 1940 cannot be explained by anthropogenic forcing alone (see Wigley, 1999). Thus, an unequivocal and direct link between the mean temperature increase and the increase in greenhouse gas concentration remains elusive at this time. The most probable cause of the observed warming now appears to be a combination of internally and externally forced natural

variability and anthropogenic sources. Li a recent comprehensive study, Bamett et al. (1999) concludes that, due to uncertainties in climate models and limitations in observed data, "a reliable weighting ofthe differentfactors contributing to the observed climate (temperature) change cannot currently be given'\

20 Uncertainties in Greenhouse Gas Induced Climate Change Figure 2.4 Changes in annual global-mean temperature (thin line) and carbon dioxide (thick line) over the past 138 years relative to the 1961-1990 average. Carbon dioxide values for 1959- 1996 are from direct measurements at Mauna Loa, Hawaii; earlier values are estimated from ice core data (dashed line). The scale for carbon dioxide is in ppmv (parts per milhon, volume basis) relative to a mean value of 333.7 (from Hurrell, 1998).

2.7 References

Akinremi, O.O., S.M. McGinn and H.W. Cutforth, 1999: Precipitation trends on the Canadian prairies. J. of Climate, 12, 2996-3003.

Angell, J. A., 1999: Comparison of surface and tropospheric temperature trends estimated from a 63-station radiosonde network, 1958-1998. Geophy. Research Letters, 26, 2761-2764.

Balling Jr. R.C., 1996: Century-long variations in United States drought severity. J. ofAgri. and Forest Meteor., 82, 293-299.

" and S. Idso, 1989: Historical temperature trends in the United States and the effect of urban population growth. J. Geophysical Research, 94, D3, 3359-3363.

Bamett, T.P., K. Hasselmann, M. Chelliah, T. Delworth, G .Hegrel, P. Jones, E. Rasmusson, E. Roeckner, C. Ropelewski, B. Santer and S. Tett, 1999: Detection and attribution of recent climate change: A status report. Bull Amer. Met. Soc, 80, 12, 2631-2659.

Bengtsson, L. 1997: A numerical simulation of anthropogenic climate change. AMBIO, 26, 58-65.

Uncertainties in Greentiouse Gas Induced Climate Change 21 Bijl, W., R. Flather, J.G. de Ronde and T. Schmith, 1999: Changing storminess? An analysis of long-term sea level data sets. Climate Research, 11, 161-172.

Cayan, D.A. and A.V. Douglas, 1984: Urban influences on surface temperatures in southwestern

United States during recent decades. J. Climate and Applied Meteor., 23, 1520-1530.

Chan, J.C.L. and Jui-en Shi, 1996: Long-term trends and interannual variability in tropical cyclone activity over the western North Pacific. Geophy. Research Letters, 23, 2765-2767.

Christy, J. R., R.W. Spencer and R.T. McNider, 1995: Reducing noise in the MSU daily lower tropospheric global temperature data sets. J. of Climate, 8, 888-896.

" and E.S. Lobl, 1998: Analysis of the merging procedure for the MSU daily temperature time series. J. of Climate, 11, 2016-2041.

Climate Diagnostics Bulletin, 1999: National Oceanic and Atmospheric Administration, U.S. Department of Commerce, Camp Springs, USA, 75p.

Dai, A., K.E. Trenberth and T.R. Karl, 1998: Global variations in droughts and wet spells. Geophysical Res. Letters, 25, 3367-3370.

Easterling, D.R., Horton, Jones, Peterson, Karl, Parker, Salinger, Razuvayev, Plummer, Jamason and Folland, 1997: Maximum and minimum temperature trends for the globe. Science, 277, 364- 367.

ESEF, 1996: The global warming debate. Report ofthe European Science and environment Forum. John Emsley (ed.). Bourn Press Ltd. Great Britain, UK, 288p.

Fujibe, F. 1995: Temperature rising trends at Japanese cities during the last one hundred years and their relationships with population, increasing rates and daily temperature ranges. Papers in Meteor, and Geophysics, 46, 35-55.

Gallo, K.P., T.W. Owen, D.R. Easterling and P.M. Jamason, 1999: Temperature trends of the U.S. historical climatology network based on satellite-designated land use/land cover. J. ofClimate, 12, 1344-1348.

Goodrich, J.D. 1992: Urban bias influences on long-term California air temperature trends. Atmospheric Environment, 26B, 1-7.

Gray, W. 1984a: Atlantic seasonal hurricane frequency. Part I: El Nino and 30 mb quasi-biennial oscillation influences. Mow. Wea. Review, 112, 1649-1668.

Gray, W. 1984b: Atlantic seasonal hurricane frequency: Part 11; forecasting its variability. Mon. Wea. Review, 112, 1669-1683.

22 Uncertainties in Greenhouse Gas Induced Climate Change " C. Landsea, P. W. Mielke and KJ. Berry, 1994: Predicting Atlantic basin seasonal tropical cyclone activity by 1 June. Weather and Forecasting, 9, 103-1 15.

Groisman, P. Ya, Karl, Easterling, Knight, Jamason, Hennessy, Suppiah, Page, Wibig, Fortuniak, Razuvaev, Douglas, Forland and Zhai, 1999: Changes in the probabihty of heavy precipitation: important indicators of climate change. Climatic Change, 42, 243-283.

Hansen, J., H. Wilson, M. Sato, R. Ruedy, K. Shah and E. Hansen, 1995: Satellite and surface data at odds? reply to John Christy and Roy Spencer. Climatic Change, 30, 103-1 17.

" , M.Sato, R. Ruedy, A. Lacis and J. Glascoe, 1998: Global climate data and models: a reconcihation: Science, 281, 930-932.

" , R. Ruedy, J. Glascoe and M. Sato, 1999: GISS analysis of surface temperature change. J. Geophysical Research, 104, D24, 30997-31022.

Hennessy, K.J., J.M. Gregory and J.F.B. Mitchell, 1997: Changes in daily precipitation under enhanced greenhouse conditions. Climate Dynamics, 13,667-680.

Hingane, L.S. 1996: Is a signature of socio-economic impact written on the climate? Climatic Change, 32, 91-102.

Hirsch, M.E., A.T. DeGaetano and S.J. Colucci,1999: The development and analysis of an east coast winter storm climatology. Preprints, Eighth Conf. On Climate Variations, 13-17 Sept. 1999, Denver, USA, Amer. Meteor. Soc, p.10-15.

Hoyt, D.V. and K.H. Schatten, 1993: A discussion on plausible solar irradiance variations. J. Geophy. Research, 98, 18895-18906.

Hurrell, J.W., 1996: Influence of variations in extratropical wintertime teleconnections on Northem Hemisphere temperature. Geophy. Res. Letters, 23, 665-668.

" 1998: Relationships among recent atmospheric circulation changes, global warming, and satellite temperatures. Science Progress, 81, 205-224.

and K.E. Trenberth, 1996: Satelhte versus surface estimates of air temperature since 1979:

J. ofClimate, 9, 2222-2232.

" 1998: Difficulties in obtaining reliable temperature trends: Reconciling the surface and satellite microwave sounding unit records. J. ofClimate, 11, 945-967.

IPCC, 1990: Climate Change: The IPCC scientific assessment. Houghton, Jenkins, and Ephraums (eds.) Cambridge Univ. Press, UK, 365p.

Uncertainties in Greenhouse Gas Induced Climate Change 23 Jones, P.D., P.Ya. Groisman, M. Coughlan, N. Plummer, W.C. Wang and T.R. Karl, 1990: Assessment of urbanization effects in time series of surface air temperature over land. Nature, 2>A1, 169-172.

T.J. Osbom, T.M.L. Wigley, P.M. Kelly and B.D. Santer, 1997: Comparison between microwave sounding unit temperature record and the surface temperature record from 1979 to

1996: Real difference or potential discontinuities? J. Geophy. Research, 102, D25, 30135-30135.

Jones, R.G., J.M. Murphy, M. Noguer and A.B. Keen, 1997: Simulation of climate change over

Europe using a nested regional-climate model. EL. Comparison of driving and regional model

responses to a doubling of carbon dioxide. Quart. J. Roy. Meteor. Soc, 123, 265-292.

Karl, T.R., H.F. Diaz and G. Kukla, 1988: Urbanization: its detection and effect in the United States

climate record. /. of Climate, 1, 1099-1 123.

" 1993: A new perspective on global warming. EOS, Am. Geoph. Union, 74, no.2, p.25.

" and P.D. Jones, 1989: Urban bias in area-averaged surface air temperature trends. Bull. Am. Met. Soc, 70, 265-267.

" and R.W. Knight, 1998: Secular trends of precipitation amount, frequency, and intensity in the United States. Bull Am. Met. Soc, 79, 231-241.

Kerr, R. A., 1999: A new force in high-latitude climate. Science, 284, 241-242.

Khandekar, M.L., T.S. Murty, D. Scott and W. Baird, 2000: The 1997 El Nino, Indonesian forest fires and Malaysian smoke problem: a deadly combination of natural and man-made hazards. Natural Hazards, (in press).

Kripalani, R.H. and A. Kulkami, 1997: Rainfall variability over south-east Asia - Connections with

Indian Monsoon and ENSO extremes: New perspectives. Int. J. of Climatology, 17, 1 155-1 168. Lambert, S. 1996: Intense extratropical Northern Hemisphere winter cyclone events: 1899-1991. J.Geophy. Research, 101, 21319-21325.

Landsea, C. W., 1993: A climatology of intense (or major) Atlantic hurricanes. Mon. Wea. Review, 121, 1703-1713.

" R.A. Pielke, A.M. Mestas-Nunez and J.A. Knaff, 1999: Atlantic basin hurricanes: indices of climate change. Climatic Change, 42, 89-129.

Lean J., J. Beer and R. Bradley, 1995: Reconstruction of solar irradiance since 1610: Implications for climate change. Geophy. Res. letters, 22, 3195-3198.

" O.R. White and A. Skumanich, 1995: On the secular ultraviolet spectral irradiance during

the Maunder Minimum. Global Biogeochemical Cycles, 9, 171-182.

24 Uncertainties in Greenhouse Gas Induced Climate Change " and D. Rind, 1 998: Climate forcing by changing solar radiation. J. of Climate, 1 1 , 3069- 3094.

Lindzen, R. 1990: Some coolness concerning global warming. Bull. Am. Met. Soc, 71, 288-203.

Manabe, S., R.J. Stouffer, M.J. Spellman and K. Bryan, 1991: Transient response of a coupled ocean-atmosphere model to gradual changes to CO2, Part I: annual mean response. J. of Climate, 4, 785-818.

" M.J. Spellman and R.J. Stouffer, 1992: Transient response of a coupled ocean-atmosphere model to gradual changes to atmospheric CO2, Part 11: seasonal response. J. of Climate, 5, 105-126.

Maybank, J. B. Bonsai, K. Jones, R. Lawford, E.G. O'Brien, E.A. Ripley and E. Wheaton, 1995: Drought as a natural disaster, Atmosphere-Ocean, 33, 195-222.

Mekis, E. and W. Hogg, 1999: Rehabilitation and analysis of Canadian daily precipitation time series. Atmosphere-Ocean, 37, 53-85.

Michaels P. D. Stookesbury, 1992: Global warming: A reduced threat? Bull. Am. Met. Soc, 73, 1563-1577.

Mitchell, J.M., 1961: The temperature of cities. Weatherwise, 14, 224-229.

Oke, T. R., 1973: City size and the urban heat island. Atmospheric Environment, 7, 769-779.

Parthasarathy, B., A.A. Munot and D.R. Kothawale, 1995: Monthly and seasonal rainfall series for all-India homogeneous regions and meteorological subdivisions: 1871-1994. Contr. Indian Inst, of Trop. Meteor., August 1995, Pune, India, 1 13p.

Plummer, N., Salinger, Nicholls, Suppiah, Hennessy, Leighton, Trewin, Page and Lough, 1999: Changes in climate extremes over the Australian region and New Zealand during the twentieth century. Climatic Change, 42, 183-202.

Prabhakara, C, R. lacovazzi, J.M. Yoo and G. Dalu, 1998: Global warming deduced from MSU. Geophy. Res. Letters, 25, 1927-1930.

Rajgopalan, B., U. Lall and M.A. Cane, 1997: Anomalous ENSO occurrences: An alternate view. J. ofClimate, 10, 2351-2357.

Rakhecha, P.R. and M.K. Soman, 1994: Trends in the annual extreme rainfall events of 1 to 3 days duration over India. Theor. and Appl Climatology, 48, 227-237.

Rappaport, E.N., 1999: Atlantic hurricane season of 1997. Mon. Wea. Review, 111, 2012-2026.

Uncertainties in Greenhouse Gas Induced Climate Change 25 Roeckner, E., T. Siebert and J. Feichter, 1995: Climatic response to anthropogenic sulfate forcing

simulated with a general circulation model. Aerosolforcing of Climate, RJ. Charlson and J. Heintzenberg (eds). John Wiley, New York, pp.349-362.

" , L, Bengston, J. Feichter, J. Leliveld and H. Rodhe, 1998: Transient climate change simulations with a coupled atmosphere-ocean GCM including the tropospheric sulphur cycle. Rep.266, Max Planck Inst, fur Meteor. Hamburg, Germany, 48p.

Saxena, V. and S. Yu, 1998: Searching for a regional fingerprint of aerosol radiative forcing in the southeastem USA. Geophy. Res. Letters, 25, 2833-2836.

" and S. Menon, 1999: Sulfate-induced cooling in the southeastem USA: An observational assessment. Geophy. Res. Letters, 26, 2489-2492.

Shabbar, A., K. Higuchi, W. Skinner and J. Knox, 1997: The association between the BWA index

and winter surface temperature variability over eastem Canada and west Greenland. Int. J. Climatology, 17, 1195-1210.

Singer, F., B.A. Boe, F.W. Decker, N. Frank, T. Gold, W. Gray, H.R. Linden, R. Lindzen, P.J. Michaels, W.A. Nirenberg, W. Porch and R. Stevenson, 1997: Comments on "Open Letter to Ben Santer" Bull. Am. Met. Soc, 78, 81-82.

" 1997: Hot Talk, Cold Science: Global warming's unfinished business. The Independent Inst, USA, lOOp.

Solanki, S.K. and M. FHgge, 1998: Solar irradiance since 1874 revisited. Geophy. Res. Letters, 25, 341-344.

Spencer, R.W. and J.R. Christy, 1990: Precise monitoring of global temperature trends fi-om satellites. Science, 247, 1558-1562.

" 1 992: Precision and radiosonde vahdation of satellite grid point temperature anomalies, n, tropospheric retrieval and trends during 1979-1990. J. ofClimate, 5, 858-866.

Thompson, D.W. and J.M. Wallace, 1998: The Arctic oscillation signature in the wintertime geopotential height and temperature fields. Geophy. Res. Letters, 25, 1297-1300.

Trenberth, K. E. and T.J. Hoar, 1996: The 1990-1995 El Nino-Southem Oscillation event: longest on record. Geophy. Res. Letters, 23, 57-60.

Von Storch, H. 1993: Changing statistics of storms in the North Atlantic. Climate trends and fixture

offshore design and operation criteria. Workshop no. 2, Reykjavik, Iceland, March 29-30, 1993.

Wentz, F. and M. Schabel, 1998: Effects of orbital decay on satellite-derived lower-tropospheric temperature trends. Nature, 394, 661-664.

26 Uncertainties in Greenhouse Gas Induced Climate Change Wigley, T.M.L., 1999: The science of climate change: Global and U.S. perspective. Report preparedfor the Pew Center on Global Climate Change, Arlington, VA, USA, June 1999, 48 pp.

Zhai, P., Sun, Ren, Liu, Gao, and Zhang, 1999: Changes of climate extremes in China. Climatic Change, 42, 203-218.

Zwiers, F. and V.V. Kharin, 1998: Changes in extremes of the climate simulated by CCC GCM2 under CO2 doubling. / of Climate, 11, 2200-2222.

Uncertainties in Greenhouse Gas Induced Climate Change 27 3.0 ASSESSMENT OF UNCERTAINTIES RELATED TO CLIMATE MODELS AND OTHER ISSUES

3.1 Introduction

Attempts at estimating the impact of increasing greenhouse gas concentration on climate using numerical general circulation models of the atmosphere have been going on for ahnost 30 years.

One of the first such studies was made almost a quarter of a century ago by Manabe and Weatherald (1975). The model they used was very simple with a Hmited computational domain, nine vertical levels in the atmosphere, a model ocean represented by an area of wet land (an area possessing an infinite source of soil moisture for evaporation), no heat transport by ocean currents, an idealized land topography, and fixed cloudiness. Despite these simplifications, the model produced results that are remarkably similar, in some ways, to some of the recent simulation studies. Specifically, the model obtained mean temperature changes (due to doubling of CO2 concentration) ofbetween 2 and 3°C in the troposphere fi-om equator to about 60**N, a pronounced warming of several degrees at high latitudes in the lower troposphere, cooling in the stratosphere from about 18 km to 30 km, and a more active hydrologic cycle as indicated by the greater rates of total precipitation and evaporation computed by the model. The pronounced warming at high latitudes was attributed to decrease in the area of snow/ice cover due to increased downward radiation from increased CO2 and consequently a pronounced warming due to a net increase in the amount of solar radiation absorbed by the earth. Similar high-latitude warming in the lower troposphere has been predicted by other climate models and reported in recent literature (e.g., Washington and Meehl, 1989; Boer et

al., 1992). Since the 1975 study, Manabe and coworkers have produced a series of papers in the recent literature detailing a variety of simulations of climate change using coupled atmosphere- ocean models of increasing complexity; some of these papers will be discussed later in this chapter.

Coupled atmosphere-ocean models are a powerfiil tool for assessing the impact of increasing greenhouse gas concentrations on present and ftiture climate. Typically, these models, which consist of several vertical levels in the atmosphere as well as in the depths of the ocean, generate changes in mean values of weather elements like pressure, temperature, and precipitation over the next 50 to 100 years by allowing the concentration of greenhouse gases to increase to two or more times their present value. With the availability of better and faster computer technology, GCMs have become more comprehensive and are able to include in their computations a variety of atmosphere-ocean processes such as heat transport by ocean currents, a number land/ocean surface processes, clouds and their spatial and temporal variation, and a detailed atmospheric chemistry involving changes in

concentration of carbon dioxide, methane, nitrous oxide, ozone, etc. Despite these advances, the

present GCMs cannot realistically include or simulate all the atmosphere-ocean processes and their ftiture changes. There are several weaknesses in the GCMs, some of which are examined below.

28 Uncertainties in Greenhouse Gas Induced Climate Change 3.2 Specific Modeling Uncertainties

3.2.1 Transient modeling: The gradual inclusion ofCO2 as the model integrates is the latest approach to climate modeling. Wl^at are the advantages and disadvantages?

The model by Manabe and Weatherald (referred to earlier) first simulated the present climate by integrating the model equations for several hundred days using the present concentration of CO2. A doubled amount of CO2 was then introduced in the model, which was integrated again for several hundred days until an equilibrium stage was reached. The differences between the two climate states were analyzed to assess the climate change. This instantaneous doubling of CO2 concentration in the model is, of course, unrealistic and has, in general, produced relatively large increases in predicted surface air temperature. Since about 1990 most climate models have been simulating climate change by gradually increasing the greenhouse gas concentrations (e.g., Manabe et al.,1990, 1991, 1992; Manabe and Stouffer, 1996; Washington and Meehl, 1989; Murphy, 1995, Murphy and Mitchell, 1995). Some of these studies (e.g., Washington and Meehl, 1989) have reported climate change results using two parallel runs, one in which the CO2 concentration was doubled instantaneously and the other in which CO2 concentration was increased gradually at the rate of 1% per year. The rate of 1% per year has been chosen for most climate models investigating the transient response because the total C02-equivalent radiative forcing of various greenhouse gases other than water vapour is currently increasing at the rate of approximately 1% per year (Hansen et al, 1988). According to Washington and Meehl (1989) the transient response (gradual increase of CO2) produces an increase in mean surface temperature of 0.7°C, while the equilibrium response (instantaneous doubling of CO2) produces an increase of 1 .6°C in the mean temperature. The difference in transient and equilibrium response has been analyzed by Manabe and Weatherald who find that the equilibrium response produces much larger increases in the mean surface temperature over higher-latitude land areas and over circumpolar oceans of the southern hemisphere. Over the northern hemispheric land areas at higher latitudes, the reduction in sea ice and snow in response to increased greenhouse gas concentrations leads to increased exchange of heat between oceanic water surfaces and overlying air and this results in larger increases in mean temperature over land areas. Over the northern regions of the North Atlantic and the circumpolar southern oceans, the transient response produces very little warming because of the location of the sinking branches of the thermohaline circulation, which carries large amount of heat (trapped by increasing atmospheric CO2) to greater depths in the oceans. The equilibrium response to doubling of CO2 is particularly large (as much as +9°C) along the coast of Antarctica due to poleward retreat of sea ice, which has high surface albedo. The changes in precipitation rates due to enhanced CO2 concentration are generally larger at high latitudes than at low latitudes because of increased poleward transport of water vapour, the equilibrium response simulating a larger poleward transport than the transient response. Consequently, the precipitation rate change at higher latitudes is larger for the equilibrium response than for the transient response.

In summary, the equilibrium response (with instantaneous doubling of CO2) produces greater temperature increases than the transient response, as well as greater precipitation rates in high- latitude land areas of the northern hemisphere. These changes are not realistic when compared with observed temperature and precipitation changes of recent years. On the plus side, the equilibrium response experiment is simpler to design and is less computer-intensive than the transient response

Uncertainties in Greenhouse Gas Induced Climate Change 29 experiment. Most climate models at present use the procedure of gradual increase in CO2 concentration to predict the climate change response.

3.2.2 Cloud feedback and aerosols: The modeling ofclouds has been a major uncertainty in climate modeling. What are the latest results in this respect? What kinds ofaerosol modeling have been included in the current GCMs? What are the strengths and weaknesses?

These two aspects will be analyzed and discussed in tandem since the two are intimately connected in the context of climate modeling. The appropriate representation of clouds and their spatial and

temporal distribution in climate models is perhaps the single most challenging task for the climate modeling community at present. The presence of clouds, even at a single level, complicates the radiation budget of the earth-atmosphere system. With the presence of multilevel clouds, absorption and transmission of solar radiation at various levels can complicate the radiation budget

considerably. Further, the presence of anthropogenic aerosols is now recognized as causing a significant perturbation on the atmospheric radiation field through direct and indirect interaction

with solar radiation. There is considerable uncertainty associated with each of these effects, making the radiation budget calculation extremely complex and often indeterminate.

hi climate modeling, the simulated cloudiness field, because of its simple dependence on relative humidity, reflects the broad features of the simulated moisture field. However, accurate simulation of the moisture field depends on the model's ability to simulate general circulation patterns in the atmosphere. Most climate models are able to simulate general circulation patterns in a broad sense only. Accordingly, they simulate large-scale cloud features like maximum cloudiness at the ITCZ (inter-tropical convergence zone at around 10-12°N) and maximum mid-latitude cloudiness at around 60°N during the boreal summer. Many regional and local features of cloud climatology are not well estimated by climate models. For example, the Canadian Climate Centre second-

generation GCM (McFarlane et al., 1992) underpredicts cloudiness in the polar regions of both hemispheres and does not produce relative maxima associated with coastal stratus clouds near the western coasts of South America and Afiica in the southern hemisphere and near the California coast in the northem hemisphere. These differences in the simulation of cloud climatology can affect regional and global radiation balances. In an important short communication, Albrecht (1989) discusses the effects of increasing aerosol concentrations (natural or man-made) on low-level cloudiness and on global albedo. It is estimated that marine stratocumulus clouds contribute about a third of the earth's global albedo (Charlson et al, 1987). Inaccurate simulation of marine stratocumulus clouds can significantly alter fiiture climate projections by GCMs; for example, Randall et al. (1984) argue that a 4% change in the amount ofmarine stratocumulus could offset the warming resulting fi-om a doubling of CO2 concentration. Further, the impact of Arctic cloud cover on the radiation budget and consequently on regional and global climate is poorly understood because of the sparsity of observations. In a comprehensive study. Curry et al. (1996) concluded that basic annual cycles of clouds and radiation characteristics in the Arctic remain uncertain enough that the precise role of cloud feedbacks cannot be determined. An analysis of several GCM products over the Arctic (Tao et al, 1995; Chen et al, 1995) reveals several discrepancies in modeled cloud cover and surface temperatures. Incomplete knowledge of cloud climatology in the

Arctic (as well as in the Antarctic) is a source of major uncertainty in climate modeling at present.

30 Uncertainties in Greenhouse Gas induced Climate Change Do clouds heat or cool the planet earth? This seemingly simple question has been a challenge to climate modelers and was settled only recently through a comprehensive analysis of satellite data jfrom ERBE (Earth Radiation Budget Experiment). The ERBE data separated the clear-sky radiation budget from the radiation budget of average cloudy skies and the difference in the two budgets is attributed to cloud radiative forcing which, according to Ramanathan (1998) is about -18 W-m"^, much larger than the 2.45 W-m""^ forcing from greenhouse gases. For ready reference, the global radiation balance for the annual mean condition together with average clear-sky and cloudy-sky radiation balances are shown in Figure 3.1a,b. It is not known at this time how this net forcing of - 18 W-m"^ (shown in Fig.S.lb) would change in response to global warming. Further, the cloud radiative forcing effects are concentrated regionally along storm-track cloud systems, over mid- latitude oceans, and in deep convective cloud systems associated with the ITCZ in the tropics and its seasonal migration (Ramanathan et al., 1998). Multilevel cloudiness and changing cloud types can reduce solar radiation significantly on a regional scale, as observed by Abakumova et al.. (1996) and Leipert and Kukla (1997).

GLOBAL ENERGY BALANCE (Unit=Wrn^)

237

Figure 3.1a Global energy balance for annual mean conditions. The top-of-the-atmosphere estimates of solar insolation (343±2 W-m"^} reflected solar radiation (106±3 W-m'^), and outgoing longwave radiation (237±3 W-m" ) are obtained from satellite data. The other quantities include atmospheric absorption of solar radiation (70-95 W-m"^), surface absorption of solar radiation (142-167 W-m"^), downward longwave emission by the atmosphere (327±15 W-m"^), upward longwave emission by the surface (390±15 W-m'^), and latent (H) and sensible (S) fluxes from the surface (from Ramanathan, 1998).

Uncertainties in Greeniiouse Gas Induced Climate Change 31 Figure 3.1b Global average clear-sky radiation budget (left panel) and average cloudy fluxes and cloud radiative forcing (right panel) from ERBE data. Outgoing arrows denote outgoing longwave radiation and incoming arrows denote net incoming solar radiation. Uncertainties in these fluxes are

±5 W-m"^ and the values are for 5-year averages between 1985 and 1989 (from Ramanathan et al., 1989).

The anthropogenic aerosols (sulfate, black carbon, dust, etc.) are now recognized as providing a significant and yet uncertain perturbation on the global radiation balance in terms of an overall

cooling effect (Charlson et al., 1990, 1992; Boucher and Lohmann, 1995). Through the direct

effect, aerosols can scatter and absorb solar radiation in cloud-free air. In terms of indirect effect, an increase in concentrations of aerosols composed of soluble substances that act as cloud condensation nuclei (CCN) would increase the cloud droplet number density, thereby reducing average cloud droplet size. This, in turn, would enhance cloud albedo (Twomey, 1977) and could suppress drizzle production and increase cloud cover and persistence (Albrecht, 1989; Leaitch et

al.,1992; Gultepe et al., 1996). There is considerable uncertainty in quantifying the impact of indirect forcing of anthropogenic aerosols because of the uncertainty in the relationship between

32 Uncertainties in Greenhouse Gas Induced Climate Change aerosol number distribution and anthropogenic pollution (Charlson et al, 1992; Penner et al., 1994). The IPCC 1996 report assigns "low" to "very low" confidence levels in the estimates of radiative forcing (direct as well as indirect) by aerosols as shown in the schematic reproduced in Figure 3.2. Several recent studies have attempted to reduce the uncertainty associated with the estimates of these forcings (Penner et al., 1998, 1999). These and other studies (Lohman and Feichter, 1997) have now provided estimates that are significantly higher than the IPCC estimates. These new values have significant implications in the global warming debate with regard to radiation balance. A few observational studies reported recently (Saxena and Yu, 1998; Saxena and Menon, 1999) have assessed these forcings to be as high as -4.8 and -4.0 W-m"^ for the southeastem USA, a region where the effect of anthropogenic aerosols could supersede that of anthropogenic greenhouse gases (IPCC, 1996) and could even produce a (small) cooling effect on a regional scale. Such a cooling effect has been documented in recent studies by Saxena and his coworkers

3 -

Halocarbons - N2O k 2 CH4 Tropospheric aerosols — direct effect CO2 Solar soot Blomass Sulphate burning lii -J- Tropospheric -p Stratospheric ozone CO ozone T

Tropospheric aerosols -2 — indirect effect

Confidence level High Low Low Low Very Very Very Very low low low low

Figure 3.2 Estimates of globally and annually averaged anthropogenic radiative forcing (W-m"^) due to changes in concentrations of greenhouse gases and aerosols since pre-industrial times and to natural changes in solar output since 1850. The height of the rectangular bar indicates a mid-range estimate of the forcing while the error bars show an estimate of the uncertainty range. Contributions of individual gases to direct greenhouse forcing are indicated in the first bar. Indirect forcing associated with depletion of stratospheric ozone and increased concentration of tropospheric ozone are shown in the second and third bars. The next three bars indicate the direct contributions of individual tropospheric aerosol components. Indirect effects are shown in the second last bar and the last bar shows an estimate of changes in radiative forcing due to variation in solar output. Confidence levels are indicated below each bar. The value and confidence level of indirect forcing has increased significantly according to recent studies (from IPCC, 1996).

Uncertainties in Greenliouse Gas Induced Climate Change 33 (referenced above). At this time, the best estimate for aerosol direct forcing on a global scale

appears to be somewhere between -0.5 and -0.9 W-m""^ (Penner et al., 1998). The estimate for

indirect aerosol forcing still appears to be associated with considerable imcertainty because of cloud dynamics and the impact of aerosols on cloud microphysics. The most recent estimate for indirect

aerosol forcing on a global scale is about -1.6 W m" (Penner et al., 1999).

Several modelers have attempted to include aerosols in climate modeling in the last few years. One

of the first attempts appears to have been by Mitchell et al. (1995a, 1995b) who included the radiative effects of both greenhouse gases and sulfate aerosols in their model. The concentrations of

greenhouse gases and sulfate aerosols used by Mitchell et al. in their model was based on the IPCC scenario IS92a (see IPCC 1996), which assumes a slow reduction in the rate ofworld economic

growth and a gradual increase in conservation measures. Mitchell et al. included only the direct radiative forcing of sulfate aerosol by increasing surface albedo. The resulting simulation produces significantly reduced global warming, about 0.2°C per decade compared to about 03°C per decade with only greenhouse gas (CO2) inclusion. Simulated mean surface temperature increase as modeled by Mitchell et al. (1995a) using the UKMO (United Kingdom Meteorological Office) model is shown in Figure 3.3a, together with similar projections using the MPI (Max Planck Institute) model. With the inclusion of direct forcing due to sulfate aerosol, the simulated global warming for the period 1900-1990 appears to be closer to observed mean warming as per the data

of Jones et al. The projected future warming for the year 2050 is about 1 .9°C, which is smaller by ahnost 1°C than the projected global warming in the presence of CO2 only. In a more recent study,

Roeckner et al., (1999) used the ECHAM4 model, which is an upgraded version of the MPI model referred to in Figure 3.3a, to simulate future global warming in an experiment that included direct as well as indirect radiative forcing of sulfate aerosol. The ECHAM4 model experiment includes

the direct radiative forcing (by sulfate aerosol) as was done by Mitchell et al., while, for the indirect effect, the cloud albedo is recalculated using empirical relationships between the cloud droplet number concentration and the sulfate aerosol mass mixing ratio. The model obtains future warming under three scenarios: 1. inclusion of GHG (CO2, CH4, N2O plus other industrial gases) only, 2. Inclusion of GHG plus sulfur cycle with direct effects only, and 3. Inclusion of GHG plus sulfur cycle with direct as well as indirect effects. Projected global warming for the year 2050 and beyond

under the three scenarios (GHG, GSD, and GSDIO) is shown in Figure 3.3b. It can be seen that the projected warming under scenario GSDIO is in closest agreement with observed warming from 1900 to the present time. This scenario projects additional warming over the next 50 years of about

1**C, a value presently considered to be probable. It may be noted that the GSDIO scenario does not include the impact of sulfate aerosol on cloud lifetime, which is likely to increase the cloud albedo further and may lead to further reduction in the projected global warming.

In summary, the impact of anthropogenic aerosols on present and future climate is still fraught with considerable uncertainty and this uncertainty has a direct impact on cloud representation in climate models. Present observational and theoretical studies may help reduce some of the uncertainties as discussed above and this may help produce more realistic simulations of future global warming and

associated climate change. Inclusion of additional aerosol and cloud processes in climate models is

likely to reduce the projected global warming. Further, if the hydrologic cycle is expected to

intensify in a warmer world, the global cloud cover is expected to increase producing additional (negative) forcing in general.

34 Uncertainties in Greenhouse Gas Induced Climate Change 1860 1900 1950 2000 2050 2100

Figure 3.3a Simulated global annual mean warming from 1 880 in two simulations with GHG forcing only (MPI(x) and UKMO(w)) and two simulations that include both GHG and direct sulfate aerosol forcing (MPI(y) and UKMO{z)) (from IPCC, 1996)

Figure 3.3b Evolution of changes in annual global mean surface air temperature compared to the mean for observations (1860-1990). Simulations are shown for three experiments: GHG,

GSD, and GSDIO. See text for details (from Roeckner et al., 1999).

Uncertainties in Greenhouse Gas Induced Climate Change 35 3.2.3 Oceans and sea ice: What are the major difficulties in incorporating oceans and sea ice in GCM's? What are the latestfindings?

A variety of oceanic processes ranging jfrom turbulent transfer of heat, momentum, and water vapour at the air-sea interface to the thermohahne circulation through the depths of the oceans are now identified as important controls of present and future climate of the earth. Also, interannual variations of sea ice over the Arctic and Antarctic regions are identified as having a significant influence on present and future climate. A comprehensive discussion on each of these processes and their inclusion in GCMs is beyond the scope of the present report. A brief discussion on some of the important processes as they relate to the issue of uncertainty is presented in the following:

One of the earliest GCMs (Manabe and Weatherald, 1975) represented the ocean by aii area of wet land or an area possessing an infinite source of soil moisture for evaporation. According to Manabe

and Weatherald, their model ocean resembles an actual ocean in the sense that it is wet but it lacks the effects of heat transport by ocean currents. This highly simphfied model ocean (often known as "swamp ocean model") was later replaced by a mixed-layer ocean, which was represented by a slab of sea water 50 m thick over which various fluxes based on climatological (or seasonal) values were prescribed. Some of the first-generation climate models developed in the mid-1980s and early

1990s (e.g., Washington and Meehl, 1984; McFarlane et al., 1992) used a mixed-layer model ocean over which various fluxes were prescribed. These fluxes were further adjusted so as to yield a realistic simulation of observed climatological SST (sea surface temperature) field over the oceanic regions of the climate model. Such a mixed-layer model ocean could not generate deep water circulation in the northern oceans of the subarctic region nor in the circumpolar oceans of the

southem hemisphere. The models of Manabe et al, (1991, 1992) were among the first to use several levels in the vertical for their model ocean which was coupled with an atmospheric model consisting of several vertical levels through the troposphere and lower stratosphere. Such coupled models were able to generate deep-water circulation, which helped extract a large amount of heat (trapped by increased CO2) from the Arctic as well as Antarctic regions and thus were able to produce a much smaller temperature increase in the higher latitudes in a doubled CO2 scenario. Most climate models today include several vertical levels for their model ocean and also consider the impact of ocean bottom topography. A climate model coupled with an ocean consisting of several vertical levels presents a formidable challenge from a computational point of view. Due to different time scales for the atmosphere and ocean (few weeks for the atmosphere versus hundreds of years for the ocean, because of its thermal inertia), the atmosphere and the ocean models are separately "spun-up" for initialization before they can be coupled using prescribed fluxes and a flux adjustment procedure. The "spin-up" of ocean models for initialization presents special computational problems requiring special "acceleration" techniques as described by Manabe et al., (1991). As the two components (atmosphere and ocean) of the climate models are coupled, the coupled system is found to drift into a new climate equilibrium state that is far removed from the observed climate. The model climate drift is a sensitive indicator of model imperfections (Sausen et al., 1988) and has to be fixed by a flux adjustment process which is briefly discussed below:

The flux adjustment problem (called a "modeling dilemma" by the IPCC 1996) arises because the values of various exchange coefficients (momentum, heat, and moisture) and their temporal and spatial variations are not completely known. The problem is fiirther compounded because of cloud parameterization in the mid-latitudes where shortwave planetary albedo is too high in the storm

36 Uncertainties in Greenhouse Gas Induced Climate Change track region (Murphy, 1995). The flux adjustment terms are calculated from the difference between the modeled surface fluxes and those required to keep the model close to current climate. After running the model for a period suitable for the calculation of average flux adjustments, these terms are applied throughout the control and the anomaly (increased GHG concentration) experiment. The main purpose of the flux adjustment is to ensure that any perturbation (due to increased CO2 for example) is applied to a reahstic reference climate so that distortion of the major climate feedback process is minimized. Strictly speaking, the flux adjustment should be relatively small in comparison to the magnitude of the flux terms; however, the flux adjustment terms in some models are almost as large as the flux terms. In the tropical and subtropical latitudes, the flux terms are of

the order of about 50 W-m'^ ; however, values in excess of 100 W-m"'^ occur in many extratropical regions, especially over the North Atlantic where many climate models reveal deficiencies in terms of insufficient heat transport by the Gulf Stream. The need for flux adjustment is expected to decrease in the future as the climate models are further improved. In a recent paper, Janssen and Viterbo (1996) investigate the role of ocean surface waves and associated momentum transfer on the atmospheric circulation. The paper further documents the impact of enhanced sea surface roughness caused by young wind seas (growing wind waves) on the mid-tropospheric atmosphere and on the storm track in the southern hemisphere in particular. The authors (Janssen and Viterbo) suggest the use of a coupled wind-wave model to ameliorate the flux adjustment problem in climate models. Another technique of applying small zonally averaged flux terms to eliminate drift of the thermohaline circulation has been suggested by Weaver and Hughes (1996).

The presence of sea ice in the Arctic and Antarctic regions and its evolution and interannual variation presents another challenging problem for climate modelers. The initial formation of sea ice depends upon the ocean surface temperature and most climate models employ the thermodynamic sea ice growth algorithm of Semtner (1976). The evolution of sea ice (growth and decay of ice mass) in most climate models is govemed primarily by a thermodynamic process and is based on a siuface energy budget and heat conduction through the ice mass. The simulation of sea ice cover in the initialization and control experiment ofmany climate models is found to be significantly less than observed sea ice cover and distribution (IPCC, 1996). Also the sea ice thickness is significantly underestimated by many climate models. For example, the CCC

(Canadian Climate Centre) model described by McFarlane et al., (1992) obtains maximum sea ice thickness of only about 2.5 m as compared to observed thicknesses of 5 m, and greater than 5 m in the Canadian archipelago (Flato, 1995). This underestimation of sea ice thickness in climate models is attributed to the exclusion of sea ice dynamics (e.g., advection and deformation) and intemal mechanics, which have a significant impacts on the sea ice distribution, ice thickness, and formation of ridges and leads. Several studies on sea ice dynamics have been reported in literature in the last 20 years (e.g., Hibler, 1979, 1980; Lemke et al., 1990; Pollard and Thompson, 1994). The study by Pollard and Thompson (1994) showed that inclusion of sea ice dynamics produced a more realistic simulation of ice extent and that the ice cover was less sensitive to C02-induced warming in the higher latitudes. The recent study by Roeckner et al., (1999) employs a sea ice model that includes ice dynamics and a viscous-plastic rheology to parameterize the stress tensor.

Their simulation of the evolution of sea ice area is much more realistic than other models (see Figure 3.4).

Uncertainties in Greenhouse Gas induced Climate Cliange 37 ^

r-T I V

till -

C 10 ^\ -

- 8 i ^. GSD

V B 6 1 — 1 I'i GHG:

' 1 i I 1 1 1 t 1 t 1 i i M « 1 1 1 t.r.. 1860 1900 1950 2000 2050 2100

Figure 3.4 Evolution of decadal mean Arctic sea ice area for three experiments: GHG, GSD, and GSDIO. The evolution of sea ice area from the control experiment is denoted by CTL (from Roeckner et al, 1999).

In summary, the atmospheric and oceanic components of climate models have been coupled with considerable success in recent years in terms of more realistic simulations of present climate;

however, there are still a number of areas of uncertainty that remain unresolved at this time. The most important issue appears to be the large flux adjustment in the North Atlantic and insufficient transport of heat from southem latitudes. These two aspects are intimately linked to the

thermohaline circulation in the North Atlantic and its possible breakdown (Rahmstorf, 1997; Weaver, 1993) as a result of increased GHG concentrations and associated wanning. Another area of major uncertainty appears to be in the inadequate simulation of sea ice cover and, in particular,

its thickness, which is significantly undersimulated. Developing appropriate numerical techniques

for coupling the atmospheric and oceanic components of the climate models is also an area that needs fiirther attention.

3.2.4 Interactive biology and chemistry: These have not been incorporated in the GCMs but are consideredpotential areas ofrefinementforfuture climate models. What do we know about the relationship ofthese to climate processes?

The impact of mcreasing concentrations of CO2 on vegetation growth has been discussed in a number of studies by Idso and others (Idso and Kimball, 1993; ESEF, 1996). A recent study (Tian

38 Uncertainties in Greenhouse Gas Induced Climate Change et al., 1998) analyzes the effect of interannual climate variability on carbon storage in Amazonian ecosystems. The study finds that during El Nino years, the Amazon ecosystems act as a source of carbon to the atmosphere, while in other years these ecosystems act as a carbon sink. These carbon fluxes are comparable in magnitude to the fluxes associated with deforestation in the Amazon basin in the early 1990s. The increasing growth of vegetation and its interannual variability could very well influence the carbon cycle and this, in turn, could influence global warming potential and projections of future climate change. The future role of the terrestrial biosphere in controlling atmospheric CO2 is difficult to predict at this time. No modeling studies that include interactive biology appear to have been reported so far.

In a recent study, Roeckner et al. (1999) include in their atmospheric chemistry model the three main GHGs, namely, CO2, CH4, and N2O plus several industrial gases such as (CFC-11, 12, 113, 114, and 115) and hydrochlorofluorocarbons (HCFC-22, 123, and 141b). The evolution of these gases is based on simple chemistry and the IS92a scenario as defined by IPCC 1996. None of the models reported so far incorporates interactive chemistry.

3.2.5 Regional climate: There have been various attempts to do climatic downscaling to estimate climate change scenarios on a regional scale. Some are statistical downscaling techniques and some are regional climate models. How much confidence can we place in the results? What are the limitations in doing this?

In recent years a number of studies have been reported using regional climate models developed on a relatively smaller computational domain covering a regional area over which climate change studies may have special particular socio-economic relevance or benefit. A review paper by Giorgi and Meams (1991) discusses various approaches to the simulation of regional climate change while a recent paper by Hewitson and Crane (1996) discusses techniques and applications of downscaling. These regional models, often identified as RCM (regional climate models) are typically developed as nested models so that the lateral boundary conditions for these models are obtained fi-om large-scale GCMs. Such a procedure can be computationally very expensive, hence a statistical downscaling approach is often used in which the local values of weather elements like temperature and precipitation are inferred fi'om either observed atmospheric predictor variables or by dynamical downscaling techniques in which the temperature and precipitation simulated at the nearest grid point of a GCM are used. There are several issues regarding the downscaling and the regional climate that are beyond the scope of this report. A very brief assessment on the utility and the limitations of downscaling will be presented here:

A recent comprehensive paper (Murphy, 1999) evaluates statistical and dynamical techniques for downscaling of local climate(temperature and precipitation) at over 900 European stations. The dynamical and statistical methods are compared in terms of correlation between the estimated and observed time series of monthly anomalies. According to Murphy, a high degree of skill is found for estimates of temperature, especially over westem, central, and northem Europe; for precipitation the skill is lower with average correlation ranging from 0.4 in summer to 0.7 in winter. Overall, the dynamical and statistical methods show similar levels of skill, although the statistical method is better for summertime estimates of temperature while the dynamical method gives slightly better estimates of wintertime precipitation. In another recent paper (Charles et al., 1999), a non- homogeneous hidden Markov model was used for downscaling of precipitation (from a hmited-area

Uncertainties in Greenhouse Gas Induced Climate Change 39 model) for southwestern Australia. The statistical technique was found to provide a credible downscaling procedure for assessing climate change in that region (southwest Australia) provided a

variable characterizing the absolute moisture content is included in the predictor set. The results further highlight that the validation of a statistical downscaling technique for present-day conditions does not necessarily imply validity for changed climate conditions. A statistical downscaling approach for estimating local changes in surface climate in the central Argentina region (South America) has been reported in another recent paper by Sohnan and Nunez (1999). The statistical approach appears to provide reasonable results, however, the prediction skill varied between the

station and the season. A paper by Laprise et al. (1998) reports on the development of a regional climate model (with grid spacing of 45 km) that is nested with a second-generation Canadian Global Climate Model (GCMII). The current climate simulation by this regional climate model appears to show some improvement over the large-scale GCMII; however, validating the current climate simulation presents difficulty due to lack of available climatology on the comparable regional scale. Further, the model simulation cannot be given too much credence because of several problems with the physical parameterization at the mesoscale.

In summary, the statistical downscaling approach appears to provide reasonable skill for assessing

climate change on a regional or a local scale. The statistical approach is certainly cost-effective as it does not involve costly development of a regional model. The utility of a regional climate model

appears to be limited at this time due to the problems noted by Laprise et al. (1998).

3.2.6 Additional modeling issues: There are a couple of additional issues regarding climate models that are not defined in the terms of reference of this project but have relevance in the present discussion. Most climate models at present are not able to adequately simulate the ENSO phenomenon and associated global weather anomalies. A recent simulation study (Timmermann et al., 1999) using the ECHAM3/LSG model obtains an ENSO-like oscillation with a characteristics period of 5 to 8 years. However, the ENSO amplitude in the model simulation is considerably underestimated. Simulating the ENSO phenomenon in all its detail is still a major challenge for the

climate modeling community. Another aspect of climate modeling that has not been discussed is the adequate representation of mesoscale oceanic and atmospheric features. Most climate models use a grid spacing that is much too coarse for simulating oceanic features like mesoscale eddies or atmospheric features like squall lines and associated convective clouds. These mesoscale features can have a significant impact on the energy and radiation budget of the atmosphere and inadequate representation of these features in a model atmosphere could lead to inaccurate projections of future climate.

3.3 Other Issues

3.3.1 Heat from fuel combustion: Burning offossilfuelproduces heat. David Taylor, an emeritus professor ofchemistry in England calculated that the heat producedfrom the burning offossilfuel may heat the atmosphere significantly. What are the latest developments on this issue?

This issue does not seem to be of serious concem to the climate change science community at this time. Only a few studies appear to have been reported in recent literature. A study by Fischer (1990)

40 Uncertainties in Greenhouse Gas Induced Climate Change makes an interesting calculation of the amount of heat generated by worldwide human activity and its relative importance with respect to global warming due to CO2. Fischer obtains a value of about 10^^ 13.8 TW (1 TW= W), which is equivalent to about 0.027 W-m"^ over the earth's surface, hi a doubled CO2 scenario, Fischer obtains a reduction of about 0.57 W-m'^ in net outgoing radiation due to absorption of longwave radiation by increased concentration of CO2 (an assumption of half- cloudy sky is used here to allow for absorption by water vapour). The heat generated directly by human activity at the time of doubled CO2 is estimated to be between 0.034 to 0.053 W-m'"^ over the earth's surface. Fischer concludes that the heat pollution by combustion of fossil fuel, biomass burning etc. can be significant compared to the CO2 absorption and should be taken into account in the global warming debate. Li a similar study by Kaufinan et al. (1991), the heating from buming of fossil fiiel and biomass (associated with deforestation) is analyzed relative to the cooling due to anthropogenic aerosols. The study concluded that the cooling effect from coal and oil buming may presently range from 0.4 to 8 times the heating effect. Within this large uncertainty, Fraser and

Mahoney hypothesize that it is presently more likely that the fossil fiiel buming causes cooling of the atmosphere rather than heating. A recent paper by Stanhill and Kalma (1995) analyzes the solar radiation depletion at Hong Kong over the last 35 years in relation to anthropogenic heating (due to explosive population growth and associated human activity) and compare the impact of this anthropogenic or urban heating on the mean annual minimum temperature increase at Hong Kong, which is about 0.6°C. Stanhill and Kalma ftirther analyze the urbanization impact on other large cities in North America and conclude that the increase in minimum temperature (of 0.6°C) was about half the urban-rural temperature difference at other similar large cities and that this was due to solar depletion at Hong Kong because of increased cloudiness and interaction of anthropogenic aerosols with clouds.

In summary, the impact of anthropogenic heating due to fossil friel buming appears to be more than offset by anthropogenic aerosol cooling at present time. According to Fraser and Mahoney (1991), future increase in coal and oil burning and the resulting increase in concentration of cloud condensation nuclei, which may saturate the cooling effect, may allow the (fossil fiiel) heating effect to dominate. The evolution of clouds and associated radiation budget may be critical in determining the impact of fossil fiiel buming in a fiiture climate.

3.4 References

Abakumova, G.M., E.M. Fiegelson, V. Russak and V.V. Stadnik, 1996: Evaluation of long-term changes in radiation, cloudiness and surface temperature on the territory of the former Soviet Union. J. of Climate, 9, 1319-1327.

Albrecht, B., 1989: Aerosols, cloud microphysics and fractional cloudiness. Science, 245, 1227- 1230.

Boer, GJ., N.A. McFarlane and M. Lozare, 1992: Greenhouse gas-induced climate change simulated with the CCC second-generation general circulation model. J. of Climate, 5, 1045-1077.

Bucher, O. and U. Lohmann, 1995: The sulfate-CCN-cloud albedo effect: A sensitivity study with two general circulation models. Tellus, 47B, 281-300.

Uncertainties in Greenhouse Gas Induced Climate Change 41 .

Charles, S.P., B.C. Bates, P.H. Whetton and J.P. Hughes, 1999: Validation of downscaling models for changed climate conditions: case study of southwestern Australia. Climate Research, 12, 1-14.

Charlson, R.J., J.E. Lovelock, M.O. Andreae and S.G. Warren, 1989: Oceanic phytoplankton, atmospheric sulphur, cloud albedo and climate. Nature, 326, 655-661.

S.E. Schwartz, J.M. Hales, R.D. Cess, J.A. Coakley, J.E. Hansen and D.J. Hoffmann, 1992: Climate forcing by anthropogenic aerosols. Science, 255, 423-430.

Chen, B., D.H. Bromwich, K.M. Hines and X. Pan, 1995: Simulation of the 1979-1988 polar climates by global climate models. Annals of Glaciology, 21, 83-

Curry, J.A., W.B. Rossow, D.A. Randall J.L. Schramm, 1996: Overview of Arctic cloud and

radiation characteristics. J. of Climate, 9, 1731-1764.

Flato, G., 1995: Spatial and temporal variability of Arctic ice thickness. Annals of Glaciology, 21, 323-32.

Fischer, G., 1990: Heat pollution and global warming. Environ. Conservation, 17, 1 17-122.

Giorgi, F. and L.O. Meams, 1991 : Approaches to the simulation of regional climate change: A review. Rev. Geophysics, 29, 191-216.

Hansen, J., Fung, Lacis, Rind, Lebedeff, Ruedy, Russell and P. Stone, 1988: Global climate change

as forecast by the Goddard Institute for Space Studies three-dimensional model. J. of Geophysical Research, 93, 9341-9364.

Hasselmann, K., L. Bengtsson, U. Cubasch, G.S. Hegrel, H. Rodhe, E. Roeckner, H.V. Storch, R.

Voss and J. Waszkewitz, 1995: Detection of anthropogenic climate change using a fingerprint method. Max Planck Institutfur Meteorologie, Report no. 168, Hamburg, Germany, 20 pp.

Hibler, W.D. HI, 1979: A dynamic thermodynamic sea ice model. J. Phy. Oceanography, 9, 817- 846.

" 1980: Modelling a variable thickness sea ice cover. Mon. Wea. Rev., 108, 1943-1973.

Hewitson, B.C. and R.G. Crane, 1996: Climate downscaling: techniques and appHcations. Climate

Research, 1, 85-95.

Idso, S.B. and B.A. Kimball, 1993: Tree growth in carbon dioxide enriched air and its implications for global carbon cycling and maximum levels for atmospheric CO2. Global Biogeochem. Cycles, 7, 537-555.

Janssen, P.A.E.M. and P. Viterbo, 1996: Ocean waves and the atmospheric climate. J. of Climate, 9, 1269-1287.

42 Uncertainties in Greenhouse Gas Induced Climate Change Kaufinan, YJ., R.S. Fraser and R.L. Mahoney, 1991 : Fossil fuel and biomass burning effect on climate- Heating or cooling? J. of Climate, 4, 578-588.

Leaitch, W.R., G.A. Issac, J.W. Strapp, C.W. Banic and H.A. Wiebe, 1992: The relationship between cloud droplet number concentration and anthropogenic pollution: Observations and climatic implications. J. of Geophysical Research, 97, 2463-2474.

Lemke, P., W.B. Owens and W.D. Hibler in, 1990: A coupled sea-ice mixed-layer pycnocline model for the Weddell Sea. J. Geophy. Res., 95, 9513-9525.

Leprise, R., D. Caya, M. Giguere, G. Bereron, H. Cote, J-P. Blanchet, G.J. Boer and N.A. McFarlane, 1998: Climate and climate change in western Canada as simulated by the Canadian Regional Climate Model. Atmosphere-Ocean, 36, 119-167.

Liepert, B.G. and G.J. Kukla, 1997: Decline in global solar radiation with increased horizontal visibility in Gemiany between 1964 and 1990. J. of Climate, 10, 2391-2401.

McFarlane, N.A., G.J. Boer, J.-P. Blanchet and M. Lazare, 1992: The Canadian Climate Centre second-generation general circulation model and its equilibrium climate. J. of Climate, 5, 1013- 1044.

Manabe, S., K. Bryan and M.J. Spelhnan, 1990: Transient response of a global ocean-atmosphere model to a doubling of atmospheric carbon dioxide. J. Physical Oceanography, 20, 722-749.

" R.J. Stouffer, M.J. Spebnan and K. Bryan, 1991 : Transient response of a coupled ocean- atmosphere to gradual changes of atmospheric CO2. Part I: Annual mean response. J. ofClimate, 4, 785-818.

" M.J. Spelman and R.J. Stouffer, 1992: Transient response of a coupled ocean-atmosphere model to gradual changes of atmospheric C02.Part 11: Seasonal response. J. of Climate, 5, 105-126.

" and R.J.Stouffer, 1 996: Low-frequency variability of surface air temperature in a 1 000 year integration of a coupled atmosphere-ocean-land surface model. J. ofClimate, 9, 376-393.

" and R.T. Weatherald, 1975: The effect of doubling CO2 concentration on the climate of GCM. J. Atmos. Sciences, 32, 3-15.

Mitchell, J.F.B., T.J. Johns, J.M. Gregory and S.B.F. Tett, 1995: Climate response to increasing levels of greenhouse gases and sulfate aerosols. Nature, 376, 501-504.

" , R.A. Davis, W.J. Ingram and C.A. Senior, 1995: On surface temperature, greenhouse gases and aerosols: models and observations. J. of Climate, 10, 2364-2386.

Murphy, J.M.,1995: Transient response of the Hadley Centre coupled ocean-atmosphere model to increasing carbon dioxide. Part I: Control climate and flux adjustment. J. ofClimate, 8, 35-56.

Uncertainties in Greenhouse Gas Induced Climate Change 43 " 1999: An evaluation of statistical and dynamical techniques for downscaling local climate. J. ofClimate, 12, 2256-2284.

" and J.F.B. Mitchell, 1 995 : Transient response of the Hadley Centre coupled ocean- atmosphere model to increasing carbon dioxide. Part 11: Spatial and temporal structure of response. J. of Climate, 8, 57-80.

Penner, J.E., R.J. Charlson, J.M. Hales, L. Laulainen, R. Leifer, T. Novakov, J. Ogren, L.F. Radke, S.E. Schwartz and L. Travis, 1994: Quantifying and minimizing uncertainty of climate forcing by anthropogenic aerosols. Bull. Am. Meteor. Soc, 75, 375-400.

" C.C. Chung and K. Grant, 1998: Climate forcing by carbonaceous and sulfate aerosols. Climate Dynamics, 14, 839-851.

" 1999: Climate change and radiative forcing by anthropogenic aerosols: a review of research during the last five years. Paper presented at the La Jolla International School of Science, The Institute for advanced physics studies. La Jolla CA, March 29-30, 1999.

Rahmstorf, S. 1997: Risk of sea-change in the Arctic. Nature, 388, 825-826.

Ramanathan, V., 1998: Trace-gas greenhouse effect and global warming. AMBIO, 27, 187-197.

R.D. Cess, E.F. Harrison, P. Minnis, B.R. Barkstrom, E. Ahmed and D. Hartmann, 1989: Cloud-radiative forcing and climate: Results fi'om the earth radiation budget experiment. Science, 243, 57-63.

Randall, D.A., J. A. Coakley, C.W. Fairall, R.A. Kropfli and D.H. Lenschow, 1984: Outlook for research on subtropical marine stratiform clouds. Bull. Amer. Meteor. Soc, 65, 1290-1301

Roeckner, E., L. Bengtsson, J. Feichter, J. LeHveld and H. Rodhe, 1999: Transient climate change simulation with a coupled atmosphere-ocean GCM including the tropospheric sulfur cycle. J. of Climate, 12, 3004-3032.

Sausen, R., K. Barthel and K. Hassehnann, 1988: Coupled ocean-atmosphere models with flux correction. Climate Dynamics, 2, 145-163.

Saxena, V. K. and S.-C. Yu, 1998: Searching for a regional fingerprint of aerosol radiative forcing in the southeastem USA. Geophy. Res. Letters, 25, 2833-2836.

" and S. Menon, 1999: Sulfate-induced cooling in the southeastem USA: An observational assessment. Geophy. Res. Letters, (in press).

Soloman, S.A. and M.N.Nunez, 1999: Local estimates of global climate change: a statistical downscaling approach. Int. J. of Climatology, 19, 835-861. Tao, X., W.L. Chapman and J.E. Walsh, 1995: Intercomparison of global climate simulations of Arctic temperature. Preprints, Fourth Conf. on Polar Meteorology and Oceanography, Dallas TX, Amer. Met. Soc. 138-143.

Tian, H., J.M. Melillo, D.W. Kicklighter, A.D. McGuire, J.V.K. Helfrich, B. Moore and C.J. Vorosmarty, 1998: Effect of interannual climate variability on carbon storage in Amazonian ecosystems. Nature, 396, 664-667.

Timmermann, A., M. Latif, A. Grotzner and R. Voss, 1999: Modes of climate variability as simulated by a coupled general circulation model. Part I: ENSO-like climate variability and its low- frequency modulation. Climate Dynamics, 15, 605-618.

Twomey, S. 1977: The influence of pollution on the shortwave albedo of clouds. J. ofAtmospheric Sciences, 34, 1149-1152.

Washington, W.M. and G.A. Meehl, 1989: Climate sensitivity due to increased CO2: experiments with a coupled atmosphere and ocean general circulation model. Climate Dynamics, 4, 1-38.

Weaver A. J., 1993: The oceans and global warming. Nature, 364, 192-193.

" and T.M.C. Hughes, 1996: On the incompatibility of ocean and atmosphere models and the need for flux adjustments. Climate Dynamics, 12, 141-170.

Uncertainties in Greenhouse Gas Induced Climate Change 45 4.0 SUMMARY OF FINDINGS AND PRIORITY AREAS OF RESEARCH

The present report examines a large number of publications reported in literature since the IPCC 1996 study, and attempts to assess the present status of several key areas of uncertainty in the science of global warming and climate change.

4.1 Summary of Important Findings

4.1.1 The annual global mean surface temperature warmed by 0.57X over the period 1861-1997 and by 0.62*'C over 1901-1997. For the recent 20-year period (1978-1997), the mean surface

temperature increase is estimated to be about 0.32°C or about 0.1 6^C per decade. The spatial distribution of the temperature change also reveals a cooling ofbetween 0.5 to 1.0°C during the last 30 years over the northwest Atlantic, parts of eastern Canada and southeastem United States, and a sizable area over central/eastern Europe and westem Russia.

4.1.2 The mean temperature trend of the troposphere as determined by the satellite-based

radiometric data over the last 19 years (1979-1997) is now estimated to be about O.TC per decade. Recent studies attempt to explain the difference in mean temperature trend between surface and troposphere in terms of cooling of the troposphere due to stratospheric ozone depletion and warming of the surface due to urbanization and land-use change.

4.1.3 The increase in mean surface temperature of the earth is attributable to a number of factors according to several recent studies. Differential changes in daily maximum and minimum temperatures and associated narrowing of the diumal temperature range (DTR) appear to provide at least a partial explanation for the mean temperature increase. Mean maximum temperature has increased over most areas of the earth with the notable exception of eastem Canada, the southem United States, portions of eastem Europe, southem China, and parts of southem South America. Mean minimum temperature has increased ahnost everywhere, except over eastem Canada and small areas of eastem Europe, the Middle East, and at most locations over Lidia. Since daily

minimimi temperature has been increasing at a higher rate (or decreasing at a lower rate) than daily maximum temperature, DTR has decreased in most areas and mean temperature has increased, by 0.3°C or more.

4.1.4 Slowly varying large-scale atmospheric oscillations like the SO (southem oscillation) and the NAO (North Atlantic oscillation), which were first identified by Sir Gilbert Walker in the 1920s, and the more recently identified Arctic oscillation (AO) can explain part of the observed increase in the mean surface temperature of the earth. The SO and the associated El Nino phenomenon have contributed to positive surface temperature anomalies over westem Canada and the northwestem United States, while the strong positive phase of the NAO has, since 1980, contributed to the positive temperature anomalies over westem Europe and parts of Eurasia.

4.1.5 Changing solar radiation since 1900 and its forcing on the earth's climate may have contributed to about half of the observed warming at the earth's surface during the last 100 years.

4.1.6 No studies reported so far have directly and unequivocally linked increased concentration of greenhouse gases to the recent increase of the mean surface temperature of the earth. The most

46 Uncertainties in Greenhouse Gas Induced Climate Change probable cause of the mean surface temperature increase is now considered to be a combination of internally and externally forced natural variability and anthropogenic sources.

4.1.7 Radiative forcing by anthropogenic aerosols (sulfate in particular) is now identified as an important perturbation that can significantly influence present and fiiture climate. The direct and indirect (cloud-mediated) radiative forcing effects due to sulfate aerosols are presently assessed to be between about -4.8 and -4.0 W-m"'^ over the southeastern USA. A few recent studies have identified a small cooling trend in the surface temperature records of the southeastern USA.

4.1.8 An analysis of precipitation data (1910-1995) over the conterminous United States shows an increasing trend in percent contribution of the upper ten percentile of daily precipitation events to the total annual precipitation, thus suggesting an intensification of the hydrologic cycle as per the IPCC 1996 conclusion. However, a similar analysis for Canada shows a decreasing trend from 1910 through 1995 for stations located in southem Canada; for stations located north of 55°N in

Canada, there is a 10% increase in precipitation based on a shorter (1940-1995) data set. Over the

Canadian prairies, there is no increase in heavy precipitation events. Elsewhere, an analysis of summer (June-September) monsoon rainfall over China and India does not show any increasing trend in extreme precipitation events. Over Indonesia/Malaysia, available data do not reveal any increasing trend in precipitation. For Australia, the percentage area of the land mass experiencing extreme wet conditions appears to have increased slightly while the area of extreme dryness has decreased slightly since 1910; however, this variability in extreme climate appears to be governed primarily by the frequency of El Nino and La Nina events. On a global or continental scale, there is no evidence of intensification of the hydrologic cycle in recent years.

4.1.9 There is no increasing trend in the total number of hurricanes occurring in the North

Atlantic, nor is there any increasing trend in the intensity of the Atlantic hurricanes as measured by the wind speed. The central Pacific region appears to show an increase in the total number of hurricanes while the Australian region shows a distinct reduction in the total number of tropical cyclones in recent years. The interannual variation of hurricanes in the Atlantic as well as in the Pacific is found to be primarily govemed by the phase of the ENSO (El Niiio southem oscillation) phenomenon.

4.1.10 Recent studies using sea-level data sets over northwest Europe indicate a small weakening of the storm climate in the southem North Sea, while, for the northem North Sea region, a small

(but insignificant) increasing trend in the storm climate is indicated.

4.1.11 There is no evidence of an increased fi-equency of El Nino (warm phase of SO) events in the equatorial central and eastem Pacific. The present La Nina (or the cold phase) of SO, which began in June 1998, is expected to continue until March 2000 according to recent analysis.

4.1.12 The predicted impact of global warming in terms of increased frequency of extreme weather events (e.g., extreme precipitation, extreme drought, increased number of tropical and extratropical storms) is not home out by recent observational studies at this time.

4.1.13 There are some issues related to fiiture scenarios that have not been discussed so far. The IPCC 1996 report developed a number of scenarios based on a wide range of assumptions

Uncertainties in Greenhouse Gas Induced Climate Change 47 regarding future economic development, population growth, and energy usage. From these scenarios, projected future atmospheric concentrations of GHG are calculated and are then used in climate models to develop projections of future climate. A recent paper by Gray (1998) makes a critical analysis of these scenarios and concludes that the IPCC scenarios exaggerate the extent of one or more factors determining future emissions. The paper further shows that the inaccuracies involved in calculating future atmospheric concentrations of GHG are so great as to render the calculations highly unreliable. This study, once again, points towards the uncertainty

involved in projecting the future evolution of atmospheric chemistry and its possible impact on future climate.

4.2 Priority Areas of Research

Based on the discussion and analysis in Chapters 2 and 3, several areas of research relating modeling and observational aspects of the science of global warming and climate change can be

identified. From an observational perspective, there are still large gaps in our observations of the atmosphere-ocean system. The Arctic and the Antarctic are probably the two best locations to exemplify the observational gaps and the uncertainty in appropriately defining the atmosphere- ocean system in these regions (see Curry et al, 1996). Other regions with observational gaps are land areas of tropical AJQica and South America, tropical southeast Asia with particular reference to the maritime continent of Indonesia, and vast areas of open ocean in the southern hemisphere. From a modeling perspective, several uncertainties exist in adequate simulation of the ENSO phenomenon and in suitable coupling of atmosphere and ocean components of the climate models.

Taking into account all these uncertainties, the following priority areas of research are identified:

4.2.1 A cooperative effort to build a data base and to construct temperature and precipitation trends on a regional, national, continental, and global basis. A detailed analysis of these trends to identify and remove any local and/or regional bias. An excellent example is a recent study (Skinner and Majorowicz,1999) on regional climatic warming and associated land-cover changes in northwest North America. The study finds an interesting correlation between the surface air temperature increase and ground surface temperature (as determined fi-om well data) and further

assesses the impact of land-cover changes on these temperatures. It would be useful to make such detailed analyses of temperature and precipitation trends in regions where local and regional influences are known to exist.

4.2.2 Improving the observational network, especially in the high-latitude regions (Arctic and Antarctic) and over tropical oceans and land areas. The recent GCOS (Global Climate Observing System) program, launched through the WMO, should help define climate and its variability in a more precise manner. This program should enable more precise detection of future climate change and more accurate attribution to natural and/or anthropogenic forcing. In a comprehensive paper Barry (1995) emphasizes the need for an observing system and data sets related to the

cryosphere in Canada. Barry further describes the cryospheric data and its geographical variety and extent in Canada as providing a valuable window on the stability of climate in the Arctic and in middle latitudes. In a recent study on climate extremes, Karl and Easterling (1999) emphasized the need for an improved data set to establish a link between climate extremes and GHG- induced climate change. As Karl and Easterling have succinctly summarized. There are a

48 Uncertainties in Greenhouse Gas Induced Climate Change number ofimpediments preventing usfrom more effectively understanding the linkages between changes in climate extremes and natural hazards to anthropogenically induced climate change.

Certainly model deficiencies are high on the list, butjust as important is our lack ofreliable long-term climatic data. Time and time again, wefind that our observing systems and data sets often have large systematic biases ofuncertain magnitude casting doubt on our ability to detect multi-decadal changes. This is why efforts like the GCOS are so critical

4.2.3 While acknowledging significant progress in climate modeling, it must be realized that most climate models have grid spacing that is too coarse to adequately represent mesoscale oceanic and atmospheric features. There is a definite need to develop finer-resolution models and develop suitable parameterization techniques to represent mesoscale features in the models. An associated problem will be to develop revised calculations of radiative forcing firom natural and anthropogenic aerosols and their interaction with mesoscale features.

4.2.4 The ENSO phenomenon is now identified as the strongest climatic signal in the coupled atmosphere-ocean system outside of the diurnal cycle. It is imperative that climate models must be improved to realistically simulate the ENSO cycle and the associated global weather anomalies. Climate models with inadequate simulation of the ENSO cycle are unlikely to produce a realistic response to perturbations fi-om increasing GHG concentration. An associated aspect is the inadequate simulation of summer monsoons fi-om Afiica to Indonesia (Gadgil and Sajnani, 1998). The ENSO-monsoon connection has been amply demonstrated in several recent studies (e.g., Webster and Yang, 1992; Khandekar, 1996). The summer monsoon and ENSO are two important features of the atmosphere-ocean system that must be adequately simulated in fixture climate models.

4.2.5 Finally, a few research areas that can be identified on a regional (provincial) scale as described in the terms of reference of this project:

a. Thoroughly analyze precipitation and temperature trends over western Canada with particular reference to Alberta and determine spatial and temporal variability of these trends.

b. Collect and carefiilly analyze data on extreme events (e.g., blizzards, heat waves, intense precipitation and associated floods, thunderstorms, and tornadoes) over western Canada and determine if any trends exist.

c. Collect data on local and regional land-use and land-cover changes and assess their influence on precipitation and temperature in particular.

4.3 References

Barry, R.G., 1995: Observing systems and data sets related to the cryosphere in Canada: A contribution to planning for the Global Climate Observing System. Atmosphere-Ocean, 33, 771- 807.

Curry, J.A., W.B. Rossow, D.A. Randall and X. Pan, 1996: Overview of Arctic cloud and radiation characteristics. J. of Climate, 9, 1731-1764.

Uncertainties in Greenhouse Gas Induced Climate Change 49 Gadgil, S. and S. Sajani, 1998: Monsoon precipitation in the AMIP runs. Climate Dynamics, 14, 659-689.

Gray, V. 1998: ThelPCC future projections: are they plausible? Climate Research, 10, 155-162.

Karl, T.R. and D.R. Easterling, 1999: Climate extremes: selected review and future research directions. Climatic Change, 42, 309-325.

Khandekar, M.L., 1996: El Nino/Southern Oscillation, Indian Monsoon and world grain yields-A synthesis. Land-based and Marine Hazards, M.I.El-Sabh et al.(eds.)., Kluwer Academic Pub. pp.79-95.

Skinner, W. R. and J.A. Majorowicz, 1999: Regional climatic warming and associated twentieth century land-cover changes in north-westem North America. Climate Research, 12, 39-52.

Webster, P. and S.Yang, 1992: Monsoon and ENSO: selectively interactive systems. Q.J.R.

Meteorol. Soc, 1 1 8, 877-926.

50 Uncertainties in Greenhouse Gas Induced Climate Change