An efficient climate model with a 3D ocean and statistical–dynamical atmosphere*

We describe a coupled climate model of intermediate complexity designed for use in global warming experiments. The atmospheric component is a two-dimensional (zonally averaged) statistical-dynamical model based on the Goddard Institute for Space Study's atmospheric general circulation model (GCM). In contrast with energy-balance models used in some climate models of intermediate complexity, this model includes full representation of the hydrological and momentum cycles. It also has parameterizations of the main physical processes, including a sophisticated radiation code. The ocean component is a coarse resolution ocean GCM with simplified global geometry based on the Geophysical Fluid Dynamics Laboratory modular ocean model. Because of the simplified geometry the resolution in the western boundary layers can be readily increased compared to conventional coarse resolution models, without increasing the model's computational requirements in a significant way. The ocean model's efficiency is also greatly increased by using a moderate degree of asynchronous coupling between the oceanic momentum and tracer fields. We demonstrate that this still allows an accurate simulation of transient behavior, including the seasonal cycle. A 100 years simulation with the model requires less than 8 hours on a state-of the art workstation. The main novelty of the model is therefore a combination of computational efficiency, statistical-dynamical atmosphere and 3D ocean. Long-term present-day climate simulations are carried out using the coupled model with and without flux adjustments, and with either the Gent-McWilliams (GM) parametrization scheme or horizontal diffusion (HD) in the ocean. Deep ocean temperatures systematically decrease in the runs without flux adjustment. We demonstrate that the mismatch between heat transports in the uncoupled states of two models is the main cause for the systematic drift. In addition, changes in the circulation and sea-ice formation also contribute to the drift. Flux adjustments in the freshwater fluxes are shown to have a stabilizing effect on the thermohaline circulation in the model, whereas the adjustments in the heat fluxes tend to weaken the global "conveyor". To evaluate the model's response to transient external forcing global warming simulations are also carried out with the flux-adjusted version of the coupled model. The coupled model reproduces reasonably well the behavior of more sophisticated coupled GCMs for both current climate and for the global warming scenarios.     

Evaluating the performance of CMIP3 and CMIP5 global climate models over the north-east Atlantic region

One of the main sources of uncertainty in estimating climate projections affected by global warming is the choice of the global climate model (GCM). The aim of this study is to evaluate the skill of GCMs from CMIP3 and CMIP5 databases in the north-east Atlantic Ocean region. It is well known that the seasonal and interannual variability of surface inland variables (e.g. precipitation and snow) and ocean variables (e.g. wave height and storm surge) are linked to the atmospheric circulation patterns. Thus, an automatic synoptic classification, based on weather types, has been used to assess whether GCMs are able to reproduce spatial patterns and climate variability. Three important factors have been analyzed: the skill of GCMs to reproduce the synoptic situations, the skill of GCMs to reproduce the historical inter-annual variability and the consistency of GCMs experiments during twenty-first century projections. The results of this analysis indicate that the most skilled GCMs in the study region are UKMO-HadGEM2, ECHAM5/MPI-OM and MIROC3.2(hires) for CMIP3 scenarios and ACCESS1.0, EC-EARTH, HadGEM2-CC, HadGEM2-ES and CMCC-CM for CMIP5 scenarios. These models are therefore recommended for the estimation of future regional multi-model projections of surface variables driven by the atmospheric circulation in the north-east Atlantic Ocean region.     

A regional response in mean westerly circulation and rainfall to projected climate warming over Tasmania, Australia

Coupled ocean–atmosphere general circulation models (GCMs) lack sufficient resolution to model the regional detail of changes to mean circulation and rainfall with projected climate warming. In this paper, changes in mean circulation and rainfall in GCMs are compared to those in a variable resolution regional climate model, the Conformal Cubic Atmospheric Model (CCAM), under a high greenhouse gas emissions scenario. The study site is Tasmania, Australia, which is positioned within the mid-latitude westerlies of the southern hemisphere. CCAM projects a different response in mean sea level pressure and mid-latitude westerly circulation to climate warming to the GCMs used as input, and shows greater regional detail of the boundaries between regions of increasing and decreasing rainfall. Changes in mean circulation dominate the mean rainfall response in western Tasmania, whereas changes to rainfall in the East Coast are less related to mean circulation changes. CCAM projects an amplification of the dominant westerly circulation over Tasmania and this amplifies the seasonal cycle of wet winters and dry summers in the west. There is a larger change in the strength than in the incidence of westerly circulation and rainfall events. We propose the regional climate model displays a more sensitive atmospheric response to the different rates of warming of land and sea than the GCMs as input. The regional variation in these results highlight the need for dynamical downscaling of coupled general circulation models to finely resolve the influence of mean circulation and boundaries between regions of projected increases and decreases in rainfall.     

Evaluation of water and energy budgets in regional climate models applied over Europe

This study presents a model intercomparison of four regional climate models (RCMs) and one variable resolution atmospheric general circulation model (AGCM) applied over Europe with special focus on the hydrological cycle and the surface energy budget. The models simulated the 15 years from 1979 to 1993 by using quasi-observed boundary conditions derived from ECMWF re-analyses (ERA). The model intercomparison focuses on two large atchments representing two different climate conditions covering two areas of major research interest within Europe. The first is the Danube catchment which represents a continental climate dominated by advection from the surrounding land areas. It is used to analyse the common model error of a too dry and too warm simulation of the summertime climate of southeastern Europe. This summer warming and drying problem is seen in many RCMs, and to a less extent in GCMs. The second area is the Baltic Sea catchment which represents maritime climate dominated by advection from the ocean and from the Baltic Sea. This catchment is a research area of many studies within Europe and also covered by the BALTEX program. The observed data used are monthly mean surface air temperature, precipitation and river discharge. For all models, these are used to estimate mean monthly biases of all components of the hydrological cycle over land. In addition, the mean monthly deviations of the surface energy fluxes from ERA data are computed. Atmospheric moisture fluxes from ERA are compared with those of one model to provide an independent estimate of the convergence bias derived from the observed data. These help to add weight to some of the inferred estimates and explain some of the discrepancies between them. An evaluation of these biases and deviations suggests possible sources of error in each of the models. For the Danube catchment, systematic errors in the dynamics cause the prominent summer drying problem for three of the RCMs, while for the fourth RCM this is related to deficiencies in the land surface parametrization. The AGCM does not show this drying problem. For the Baltic Sea catchment, all models similarily overestimate the precipitation throughout the year except during the summer. This model deficit is probably caused by the internal model parametrizations, such as the large-scale condensation and the convection schemes.     

Soil moisture prediction over the Australian continent

Summary This paper describes an attempt to model soil moisture over the Australian continent with an integrated system of dynamic models and a Geographic Information System (GIS) data base. A land surface scheme with improved treatment of soil hydrological processes is described. The non-linear relationships between soil hydraulic conductivity, matric potential and soil moisture are derived from the Broadbridge and White soil model. For a single location, the prediction of the scheme is in good agreement with the measurements of the Hydrological and Atmospheric Pilot Experiment (HAPEX). High resolution atmospheric and geographic data are used in soil moisture prediction over the Australian continent. The importance of reliable land surface parameters is emphasized and details are given for deriving the parameters from a GIS. Predicted soil moisture patterns over the Australian continent in summer, with a 50 km spatial resolution, are found to be closely related to the distribution of soil types, apart from isolated areas and times under the influence of precipitation. This is consistent with the notion that the Australian continent in summer is generally under water stress. In contrast, predicted soil temperatures are more closely related to radiation patterns and changes in atmospheric circulation. The simulation can provide details of soil moisture evolution both in space and time, that are very useful for studies of land use sustainability, such as plant growth modelling and soil erosion prediction.With 12 Figures     

The sensitivity of the tropical circulation and Maritime Continent precipitation to climate model resolution

The dependence of the annual mean tropical precipitation on horizontal resolution is investigated in the atmospheric version of the Hadley Centre General Environment Model. Reducing the grid spacing from about 350 km to about 110 km improves the precipitation distribution in most of the tropics. In particular, characteristic dry biases over South and Southeast Asia including the Maritime Continent as well as wet biases over the western tropical oceans are reduced. The annual-mean precipitation bias is reduced by about one third over the Maritime Continent and the neighbouring ocean basins associated with it via the Walker circulation. Sensitivity experiments show that much of the improvement with resolution in the Maritime Continent region is due to the specification of better resolved surface boundary conditions (land fraction, soil and vegetation parameters) at the higher resolution. It is shown that in particular the formulation of the coastal tiling scheme may cause resolution sensitivity of the mean simulated climate. The improvement in the tropical mean precipitation in this region is not primarily associated with the better representation of orography at the higher resolution, nor with changes in the eddy transport of moisture. Sizeable sensitivity to changes in the surface fields may be one of the reasons for the large variation of the mean tropical precipitation distribution seen across climate models.     

MRI模式对华南春雨气候态及年际变率的模拟：不同模式分辨率的比较

本文利用日本气象研究所（MRI）参加第五次国际耦合模式比较计划（CMIP5）的大气环流模式在高、中、低三种分辨率下的AMIP试验结果,评估了其对华南春雨气候态和年际变率的模拟能力,比较了不同分辨率的模拟结果。结果表明,三种不同水平分辨率（120 km、60 km和20 km）的模式均能再现北半球春季位于中国东南部的降水中心。相较于120 km模式,20 km模式能够更为合理地模拟出华南春雨位于南岭—武夷山脉的降水中心。水汽收支分析表明,60 km、20 km模式高估了水汽辐合,使得华南春雨的降水强度被高估。在年际变率方面,在三种分辨率下,模式均能较好地再现观测中El Ni?o衰减年春季的西北太平洋反气旋以及华南春雨降水正异常。较之120 km模式,60 km、20 km模式模拟的降水正异常的空间分布和强度更接近观测,原因是后者模拟的El Ni?o衰减年春季华南地区的水平水汽平流异常更接近观测。本研究表明,发展高分辨率气候模式是提高华南春雨的气候态和年际变率模拟水平的有效途径之一。     

Diagnostic metrics for evaluation of annual and diurnal cycles

Two sets of diagnostic metrics are proposed for evaluation of global models?? simulation of annual and diurnal cycles of precipitation. The metrics for the annual variation include the annual mean, the solstice and equinoctial asymmetric modes of the annual cycle (AC), and the global monsoon precipitation domain and intensity. The metrics for the diurnal variation include the diurnal range, the land?Csea contrast and transition modes of the diurnal cycle (DC), and the diurnal peak propagation in coastal regions. The proposed modes for the AC and DC represent faithfully the first two leading empirical orthogonal functions and explain, respectively, 82% of the total annual variance and 87% of the total diurnal variance over the globe between 45°S and 45°N. The simulated AC and DC by the 20-km-mesh MRI/JMA atmospheric general circulation model (AGCM) are in a wide-ranging agreement with observations; the model considerably outperforms any individual AMIP II GCMs and has comparable performance to 12-AMIP II model ensemble simulation measured by Pearson??s pattern correlation coefficient. Comparison of four versions of the MRI/JMA AGCM with increasing resolution (180, 120, 60, and 20?km) reveals that the 20-km version reproduces the most realistic annual and diurnal cycles. However, the improved performance is not a linear function of the resolution. Marked improvement of the simulated DC (AC) occurs at the resolution of 60?km (20?km). The results suggest that better represented parameterizations that are adequately tuned to increased resolutions may improve models?? simulation on the forced responses. The common deficiency in representing the monsoon domains suggests the models having difficulty in replicating annual march of the Subtropical Highs that is largely driven by prominent east-west land?Cocean thermal contrast. Note that the 20-km model reproduces realistic diurnal cycle, but fails to capture realistic Madden-Julian Oscillation.     

Albedo enhancement of marine clouds to counteract global warming: impacts on the hydrological cycle

Recent studies have shown that changes in solar radiation affect the hydrological cycle more strongly than equivalent CO₂ changes for the same change in global mean surface temperature. Thus, solar radiation management ??geoengineering?? proposals to completely offset global mean temperature increases by reducing the amount of absorbed sunlight might be expected to slow the global water cycle and reduce runoff over land. However, proposed countering of global warming by increasing the albedo of marine clouds would reduce surface solar radiation only over the oceans. Here, for an idealized scenario, we analyze the response of temperature and the hydrological cycle to increased reflection by clouds over the ocean using an atmospheric general circulation model coupled to a mixed layer ocean model. When cloud droplets are reduced in size over all oceans uniformly to offset the temperature increase from a doubling of atmospheric CO₂, the global-mean precipitation and evaporation decreases by about 1.3% but runoff over land increases by 7.5% primarily due to increases over tropical land. In the model, more reflective marine clouds cool the atmospheric column over ocean. The result is a sinking motion over oceans and upward motion over land. We attribute the increased runoff over land to this increased upward motion over land when marine clouds are made more reflective. Our results suggest that, in contrast to other proposals to increase planetary albedo, offsetting mean global warming by reducing marine cloud droplet size does not necessarily lead to a drying, on average, of the continents. However, we note that the changes in precipitation, evaporation and P-E are dominated by small but significant areas, and given the highly idealized nature of this study, a more thorough and broader assessment would be required for proposals of altering marine cloud properties on a large scale.     

The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection

The LMDZ4 general circulation model is the atmospheric component of the IPSL–CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley–Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke’s convection scheme, used in previous versions, the Emanuel’s scheme improves the representation of the Hadley–Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke’s parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.     

Regional‐scale surface hydrologic simulations from global climate models: A case study

Abstract

Current understanding of the regional nature of global changes in the climato‐logical regime of the earth is limited. General circulation climate models (GCMs) cannot provide consistent and detailed information on the regional patterns of precipitation, soil moisture and runoff that are required by water resource planners. A case study is presented that couples a limited area model, with high spatial resolution and realistic land‐surface parametrization, to a global climate model. Results for July are presented for the continent of Australia.     

High-resolution simulations of global climate, part 2: effects of increased greenhouse cases

We report results from the highest-resolution simulations of global warming yet performed with an atmospheric general circulation model. We compare the climatic response to increased greenhouse gases of the National Center for Atmospheric Research (NCAR) climate model, CCM3, at T42 and T170 resolutions (horizontal grid spacing of 300 and 75 km respectively). All simulations use prescribed sea surface temperatures (SST). Simulations of the climate of 2100 ad use SSTs based on those from NCAR coupled model, Climate System Model (CSM). We find that the global climate sensitivity and large-scale patterns of climate change are similar at T42 and T170. However, there are important regional scale differences that arise due to better representation of topography and other factors at high resolution. Caution should be exercised in interpreting specific features in our results both because we have performed climate simulations using a single atmospheric general circulation model and because we used with prescribed sea surface temperatures rather than interactive ocean and sea-ice models.     

不同分辨率对淮河流域连续暴雨过程影响的中尺度模拟试验

利用NCAR 和NOAA 发展的新一代中尺度模式WRF(Weather Research and Forecast),对2003年7 月上旬发生在我国淮河流域三个连续暴雨过程(7 月1—11 日)进行了数值模拟试验,研究的重点是了解不同水平分辨率(45、30、20 和10 km)对WRF 模拟结果的影响。模拟结果与观测的比较表明,WRF 模式能够合理地模拟不同时段的降水带以及平均环流形势的分布特征,对于区域平均等压面上的物理量也有较好的模拟性能。不同分辨率的模拟结果比较表明:不同分辨率对降水的模拟效果影响较大,提高模式水平分辨率有助于预报效果的改善,但高分辨率模拟的降水强度偏强,空报偏多;不同分辨率对环流形势的模拟效果影响不大,各个分辨率的低层风场误差都存在一个5～6 天的波动,并向模式的中高层传播,传播速度约为3 天。      

GRAPES气象-水文模式在一次洪水预报中的应用

尝试将GRAPES (Global-Regional Assimilation and PrEdiction System) 模式与水文模型结合,构建GRAPES气象-水文单向耦合模式,进行洪水预报。气象模式选取GRAPES_Meso模式,分别采用15 km×15 km和5 km×5 km水平分辨率,15 km×15 km的GRAPES模式由NCEP全球预报场提供初始场和侧边界条件;5 km×5 km的GRAPES模式由15 km×15 km GRAPES模式提供初始场和侧边界条件,将GRAPES_Meso模式的定量降水预报分辨率统一降尺度到5 km×5 km分辨率,用于驱动水文模式。水文模型选取新安江模型与分布式新安江模型。以淮河王家坝站以上流域和息县流域为试验流域,将GRAPES降水预报场驱动水文模型进行单向耦合,构建GRAPES气象-水文单向耦合模式,选择2009年8月28日08:00(北京时,下同)—9月9日14:00汛期一次洪水过程,进行实际预报试验。结果表明:15 km×15 km和5 km×5 km的GRAPES模式预报降水与实况降水分布相一致;与水文站观测降水驱动水文模型洪水模拟结果相比,GRAPES气象-水文模式对洪水预报的预见期延长效果明显,对洪水模拟精度也较高,与水文模型输入场分辨率要求相匹配的降水产品对洪水模拟的精度更高。     

Tropical Pacific internal atmospheric dynamics and resolution in a coupled GCM

A noise reduction technique, namely the interactive ensemble (IE) approach is adopted to reduce noise at the air–sea interface due to internal atmospheric dynamics in a state-of-the-art coupled general circulation model (CGCM). The IE technique uses multiple realization of atmospheric general circulation models coupled to a single ocean general circulation model. The ensembles mean fluxes from the atmospheric simulations are communicated to the ocean component. Each atmospheric simulation receives the same SST coming from the ocean component. The only difference among the atmospheric simulations comes from perturbed initial conditions, thus the atmospheric states are, in principle synoptically independent. The IE technique can be used to better understand the importance of weather noise forcing of natural variability such as El Niño Southern Oscillation (ENSO). To study the impact of weather noise and resolution in the context of a CGCM, two IE experiments are performed at different resolutions. Atmospheric resolution is an important issue since the noise statistics will depend on the spatial scales resolved. A simple formulation to extract atmospheric internal variability is presented. The results are compared to their respective control cases where internal atmospheric variability is left unchanged. The noise reduction has a major impact on the coupled simulation and the magnitude of this effect strongly depends on the horizontal resolution of the atmospheric component model. Specifically, applying the noise reduction technique reduces the overall climate variability more effectively at higher resolution. This suggests that “weather noise” is more important in sustaining climate variability as resolution increases. ENSO statistics, dynamics, and phase asymmetry are all modified by the noise reduction, in particular ENSO becomes more regular with less phase asymmetry when noise is reduced. All these effects are more marked for the higher resolution case. In contrast, ENSO frequency is unchanged by the reduction in the weather noise, but its phase-locking to the annual cycle is strongly dependent on noise and resolution. At low resolution the noise structure is similar to the signal, whereas the spatial structure of the noise deviates from the spatial structure of the signal as resolution increases. It is also suggested that event-to-event differences are largely driven by atmospheric noise as opposed to chaotic dynamics within the context of the large-scale coupled system, suggesting that there is a well-defined “canonical” event.     

Importance of instantaneous radiative forcing for rapid tropospheric adjustment

To better understand CFMIP/CMIP inter-model differences in rapid low cloud responses to CO₂ increases and their associated effective radiative forcings, we examined the tropospheric adjustment of the lower tropospheric stability (LTS) in three general circulation models (GCMs): HadGEM2-A, MIROC3.2 medres, and MIROC5. MIROC3.2 medres showed a reduction in LTS over the sub-tropical ocean, in contrast to the other two models. This reduction was consistent with a temperature decrease in the mid-troposphere. The temperature decrease was mainly driven by instantaneous radiative forcing (RF) caused by an increase in CO₂. Reductions in radiative and latent heating, due to clouds, and in adiabatic and advective heating, also contribute to the temperature decrease. The instantaneous RF in the mid-troposphere in MIROC3.2 medres is inconsistent with the results of line-by-line (LBL) calculations, and thus it is considered questionable. These results illustrate the importance of evaluating the vertical profile of instantaneous RF with LBL calculations; improved future model performance in this regard should help to increase our confidence in the tropospheric adjustment in GCMs.     

Dynamical downscaling: Fundamental issues from an NWP point of view and recommendations

Dynamical downscaling has been recognized as a useful tool not only for the climate community, but also for associated application communities such as the environmental and hydrological societies. Although climate projection data are available in lower-resolution general circulation models (GCMs), higher-resolution climate projections using regional climate models (RCMs) have been obtained over various regions of the globe. Various model outputs from RCMs with a high resolution of even as high as a few km have become available with heavy weight on applications. However, from a scientific point of view in numerical atmospheric modeling, it is not clear how to objectively judge the degree of added value in the RCM output against the corresponding GCM results. A key factor responsible for skepticism is based on the fundamental limitations in the nesting approach between GCMs and RCMs. In this article, we review the current status of the dynamical downscaling for climate prediction, focusing on basic assumptions that are scrutinized from a numerical weather prediction (NWP) point of view. Uncertainties in downscaling due to the inconsistencies in the physics packages between GCMs and RCMs were revealed. Recommendations on how to tackle the ultimate goal of dynamical downscaling were also described.     

Using moisture conservation to evaluate oceanic surface freshwater fluxes in climate models

We apply a diagnostic based on moisture conservation in the atmosphere, integrated over planetary-scale ocean basins and drainage areas to evaluate freshwater fluxes over the ocean surface to three generations of the Hadley Centre climate model (HadCM3, HadGEM1 and HadGEM2-AO). The coherent inclusion of runoff by the diagnostic enables model surface freshwater fluxes to be compared directly with observational estimates of precipitation, evaporation and river discharge. We also introduce a normalised metric, based on model-observation RMS differences, to assess the representation of the fluxes by the model. This methodology could be a powerful tool for evaluating model performance during future model development and model intercomparison exercises. Using this diagnostic, and defining the drainage areas from the global river routing model TRIP, we obtain large-scale surface oceanic fluxes from ERA40 and NCAR-NCEP reanalysis data, which we compare with analogous budgets computed from a set of individual observational estimates of evaporation, precipitation and river discharge. The sum of errors in the Hadley Centre climate model in all ocean basins suggests a steady improvement over the three generations of the model. However, an analysis of sources and sinks of water vapour shows common errors in the models, like an excess of evaporation in the tropical-subtropical Atlantic, and a surplus of water vapour export from tropical-subtropical areas to the mid-latitude regions, making the oceanic surface fluxes too fresh at mid latitudes. Errors in the models are consistent with an excessively strong hydrological cycle.     

Global 60?km simulations with CCAM: evaluation over the tropics

A six-member ensemble of 60?km resolution global atmospheric simulations has been performed for studying future climate scenarios of Pacific island nations. The simulations were performed using the CSIRO Conformal Cubic Atmospheric Model (CCAM), driven by bias-corrected sea surface temperatures (SSTs) provided by six Coupled Model Intercomparison Project phase 3 global climate models (GCMs) from the Intergovernmental Panel on Climate Change Fourth Assessment Report for the period 1971–2100. This paper focuses on results for the representation of the current climate in the tropical region, a region where the “cold tongue” problem is apparent in all host GCMs. The SST bias-correction and the fine horizontal resolution employed in the CCAM simulations produce a significant improvement over the host GCMs in the rainfall patterns for the transient seasons March–April–May and September–October–November, and a moderate improvement for December–January–February and June–July–August. CCAM also simulates improved rainfall patterns over the South Pacific Convergence Zone. The performance of other tropical features, such as El Ni?o Southern Oscillation and the Walker circulation, is also evaluated.     

Laurentian Great Lakes double-CO2 climate change hydrological impacts

The Great Lakes Environmental Research Laboratory has developed conceptual daily models for simulating moisture storages in and runoff from the 121 watersheds draining into the Laurentian Great Lakes, over-lake precipitation into each lake, and the heat storages in and evaporation from each lake. We combine these components as net basin supplies for each lake to consider climate change scenarios developed from atmospheric general circulation models (GCMs). Recent scenarios of a doubling of atmospheric Co₂, available from the Goddard Institute for Space Studies, the Geophysical Fluid Dynamics Laboratory, and Oregon State University are considered by making changes in historical meteorological data similar to the changes observed in the GCMs, observing the impact of the changed data in the model outputs, and comparing outputs to model results using unchanged data, representing comparison to an unchanged atmosphere. This study indicates a 23 to 51% reduction in net basin supplies to all the Great Lakes; there is significant variation in the components of these supplies among the three GCMs. The basins various moisture storages become dryer and the lakes are warmer with associated hydrological impacts.GLERL Contribution NO. 646.