Search

Standard vector calculus formulas of Cartesian three space are projected onto the surface of a sphere. This produces symmetric equations with three nonindependent horizontal velocity components. Each orthogonal axis has a velocity component that rotates around its axis (eastward velocity rotates around the north–south axis) and a specific angular momentum component that is the product of the velocity component multiplied by the cosine of axis’ latitude. Angular momentum components align with the fixed axes and simplify several formulas, whereas the rotating velocity components are not orthogonal and vary with location. Three symmetric coordinates allow vector resolution and calculus operations continuously over the whole spherical surface, which is not possible with only two coordinates. The symmetric equations are applied to one-layer shallow water models on cubed-sphere and icosahedral grids, the latter being computationally simple and applicable to an ocean domain. Model results are presented for three different initial conditions and five different resolutions.

The key uncertainties in the climate sensitivity to the increasing greenhouse gases lie in the behavior and impact of short-lived species, such as tropospheric aerosols and ozone, and secondly, in the response and impact of the ocean circulation.

We use the new version of NASA Goddard Institute for Space Studies (GISS) climate model, modelE2 with 2º by 2.5º horizontal resolution and 40 vertical layers, with the model top at 0.1 hPa [Schmidt et al., 2014]. We use two different treatments of the atmospheric composition and aerosol indirect effect: (1) TCAD(I) version has fully interactive Tracers of Aerosols and Chemistry in both the troposphere and stratosphere. This model predicts total aerosol number and mass concentrations [Shindell et al., 2013]; (2) TCAM is the aerosol microphysics and chemistry model based on the quadrature methods of moments [Bauer et al., 2008]. Both TCADI and TCAM models include the first indirect effect of aerosols on clouds [Menon et al., 2010]; the TCAD model includes only the direct aerosol effect. We consider the results of the TCAD, TCADI and TCAM models coupled to "Russell ocean model" [Russell et al., 1995], E2-R. We examine the climate response for the "historical period" that include the natural and anthropogenic forcings for 1850 to 2012. The effect of clouds, their feedbacks, as well as the aerosol-cloud interactions are assessed for the transient climate change.

Extreme monsoon rainfall in India has disastrous consequences, including significant socio- economic impacts. However, little is known about the overall trends and climate factors associated with extreme rainfall because rainfall greatly varies across India and because few appropriate methods are available to measure extreme rainfall in the context of such heterogeneity. To provide a comprehensive assessment of extreme monsoon rainfall, we developed a metric using record rainfall data to measure the changes in the likelihood of extreme high and extreme low rainfall over time; this metric is independent of the characteristics of the underlying rainfall distributions. Hence, the metric is ideally suited to aggregate extreme rainfall information across heterogeneous regions covering India. We found that from 1930 to 2013, the likelihood of extreme high and extreme low rainfall increases 2-fold and 4-fold, respectively. These overall trend increases are driven by anomalous increases, particularly in the early 2000s; the likelihood of extreme high and extreme low rainfall increases 5-fold and 18-fold in 2005 and 2002, respectively. These findings imply a broadening of the underlying monsoon rainfall distribution over the past century. We also show that the time patterns of the likelihood of extreme rainfall in recent decades are correlated with the El Nino Southern Oscillation, Indian Ocean Dipole, and surface air temperature in the Northern Hemisphere.

The adjustment of parameterized gravity waves associated with model convection and finer vertical resolution has made possible the generation of the quasi-biennial oscillation (QBO) in two Goddard Institute for Space Studies (GISS) models, GISS Middle Atmosphere Global Climate Model III and a climate/middle atmosphere version of Model E2. Both extend from the surface to 0.002 hPa, with 2deg × 2.5deg resolution and 102 layers. Many realistic features of the QBO are simulated, including magnitude and variability of its period and amplitude. The period itself is affected by the magnitude of parameterized convective gravity wave momentum fluxes and interactive ozone (which also affects the QBO amplitude and variability), among other forcings. Although varying sea surface temperatures affect the parameterized momentum fluxes, neither aspect is responsible for the modeled variation in QBO period. Both the parameterized and resolved waves act to produce the respective easterly and westerly wind descent, although their effect is offset in altitude at each level. The modeled and observed QBO influences on tracers in the stratosphere, such as ozone, methane, and water vapor are also discussed. Due to the link between the gravity wave parameterization and the models' convection, and the dependence on the ozone field, the models may also be used to investigate how the QBO may vary with climate change.

Simulations of the preindustrial and doubled CO2 climates are made with the GISS Global Climate Middle Atmosphere Model 3 using two different estimates of the absolute solar irradiance value: a higher value measured by solar radiometers in the 1990s and a lower value measured recently by the Solar Radiation and Climate Experiment. Each of the model simulations is adjusted to achieve global energy balance; without this adjustment the difference in irradiance produces a global temperature change of 0.48C, comparable to the cooling estimated for the Maunder Minimum. The results indicate that by altering cloud cover the model properly compensates for the different absolute solar irradiance values on a global level when simulating both preindustrial and doubled CO2 climates. On a regional level, the preindustrial climate simulations and the patterns of change with doubled CO2 concentrations are again remarkably similar, but there are some differences. Using a higher absolute solar irradiance value and the requisite cloud cover affects the model's depictions of high-latitude surface air temperature, sea level pressure, and stratospheric ozone, as well as tropical precipitation. In the climate change experiments it leads to an underestimation of North Atlantic warming, reduced precipitation in the tropical western Pacific, and smaller total ozone growth at high northern latitudes. Although significant, these differences are typically modest compared with the magnitude of the regional changes expected for doubled greenhouse gas concentrations. Nevertheless, the model simulations demonstrate that achieving the highest possible fidelity when simulating regional climate change requires that climate models use as input the most accurate (lower) solar irradiance value.

We present a description of the ModelE2 version of the Goddard Institute for Space Studies (GISS) General Circulation Model (GCM) and the configurations used in the simulations performed for the Coupled Model Intercomparison Project Phase 5 (CMIP5). We use six variations related to the treatment of the atmospheric composition, the calculation of aerosol indirect effects, and ocean model component. Specifically, we test the difference between atmospheric models that have noninteractive composition, where radiatively important aerosols and ozone are prescribed from precomputed decadal averages, and interactive versions where atmospheric chemistry and aerosols are calculated given decadally varying emissions. The impact of the first aerosol indirect effect on clouds is either specified using a simple tuning, or parameterized using a cloud microphysics scheme. We also use two dynamic ocean components: the Russell and HYbrid Coordinate Ocean Model (HYCOM) which differ significantly in their basic formulations and grid. Results are presented for the climatological means over the satellite era (1980-2004) taken from transient simulations starting from the preindustrial (1850) driven by estimates of appropriate forcings over the 20th Century. Differences in base climate and variability related to the choice of ocean model are large, indicating an important structural uncertainty. The impact of interactive atmospheric composition on the climatology is relatively small except in regions such as the lower stratosphere, where ozone plays an important role, and the tropics, where aerosol changes affect the hydrological cycle and cloud cover. While key improvements over previous versions of the model are evident, these are not uniform across all metrics.

Observations of climate change during the CMIP5 extended historical period (1850-2012) are compared to trends simulated by six versions of the NASA Goddard Institute for Space Studies ModelE2 Earth System Model. The six models are constructed from three versions of the ModelE2 atmospheric general circulation model, distinguished by their treatment of atmospheric composition and the aerosol indirect effect, combined with two ocean general circulation models, HYCOM and Russell. Forcings that perturb the model climate during the historical period are described. Five-member ensemble averages from each of the six versions of ModelE2 simulate trends of surface air temperature, atmospheric temperature, sea ice and ocean heat content that are in general agreement with observed trends, although simulated warming is slightly excessive within the past decade. Only simulations that include increasing concentrations of long-lived greenhouse gases match the warming observed during the twentieth century. Differences in twentieth-century warming among the six model versions can be attributed to differences in climate sensitivity, aerosol and ozone forcing, and heat uptake by the deep ocean. Coupled models with HYCOM export less heat to the deep ocean, associated with reduced surface warming in regions of deepwater formation, but greater warming elsewhere at high latitudes along with reduced sea ice. All ensembles show twentieth-century annular trends toward reduced surface pressure at southern high latitudes and a poleward shift of the midlatitude westerlies, consistent with observations.

How does climate sensitivity vary with the magnitude of climate forcing? This question was investigated with the use of a modified coupled atmosphere-ocean model, whose stability was improved so that the model would accommodate large radiative forcings yet be fast enough to reach rapid equilibrium. Experiments were performed in which atmospheric CO2 was multiplied by powers of 2, from 1/64 to 256 times the 1950 value. From 8 to 32 times, the 1950 CO2, climate sensitivity for doubling CO2 reaches 8 C due to increases in water vapor absorption and cloud top height and to reductions in low level cloud cover. As CO2 amount increases further, sensitivity drops as cloud cover and planetary albedo stabilize. No water vapor-induced runaway greenhouse caused by increased CO2 was found for the range of CO2 examined. With CO2 at or below 1/8 of the 1950 value, runaway sea ice does occur as the planet cascades to a snowball Earth climate with fully ice covered oceans and global mean surface temperatures near 30 C.

The authors simulate transient twentieth-century climate in the Goddard Institute for Space Studies (GISS) GCM, with aerosol and ozone chemistry fully coupled to one another and to climate including a full dynamic ocean. Aerosols include sulfate, black carbon (BC), organic carbon, nitrate, sea salt, and dust. Direct and BC snow-albedo radiative effects are included. Model BC and sulfur trends agree fairly well with records from Greenland and European ice cores and with sulfur deposition in North America; however, the model underestimates the sulfur decline at the end of the century in Greenland. Global BC effects peak early in the century (1940s); afterward the BC effects decrease at high latitudes of the Northern Hemisphere but continue to increase at lower latitudes. The largest increase in aerosol optical depth occurs in the middle of the century (1940s-80s) when sulfate forcing peaks and causes global dimming. After this, aerosols decrease in eastern North America and northern Eurasia leading to regional positive forcing changes and brightening. These surface forcing changes have the correct trend but are too weak. Over the century, the net aerosol direct effect is -0.41 Watts per square meter, the BC-albedo effect is -0.02 Watts per square meter, and the net ozone forcing is +0.24 Watts per square meter. The model polar stratospheric ozone depletion develops, beginning in the 1970s. Concurrently, the sea salt load and negative radiative flux increase over the oceans around Antarctica. Net warming over the century is modeled fairly well; however, the model fails to capture the dynamics of the observedmidcentury cooling followed by the late century warming.Over the century, 20% of Arctic warming and snow ice cover loss is attributed to the BC albedo effect. However, the decrease in this effect at the end of the century contributes to Arctic cooling. To test the climate responses to sulfate and BC pollution, two experiments were branched from 1970 that removed all pollution sulfate or BC. Averaged over 1970-2000, the respective radiative forcings relative to the full experiment were +0.3 and -0.3 Watts per square meter; the average surface air temperature changes were +0.2 degrees and -0.03 C. The small impact of BC reduction on surface temperature resulted from reduced stability and loss of low-level clouds.

Ample physical evidence shows that carbon dioxide (CO2) is the single most important climate-relevant greenhouse gas in Earth s atmosphere. This is because CO2, like ozone, N2O, CH4, and chlorofluorocarbons, does not condense and precipitate from the atmosphere at current climate temperatures, whereas water vapor can and does. Noncondensing greenhouse gases, which account for 25% of the total terrestrial greenhouse effect, thus serve to provide the stable temperature structure that sustains the current levels of atmospheric water vapor and clouds via feedback processes that account for the remaining 75% of the greenhouse effect. Without the radiative forcing supplied by CO2 and the other noncondensing greenhouse gases, the terrestrial greenhouse would collapse, plunging the global climate into an icebound Earth state.

The USCLI VAR working group on drought recently initiated a series of global climate model simulations forced with idealized SST anomaly patterns, designed to address a number of uncertainties regarding the impact of SST forcing and the role of land-atmosphere feedbacks on regional drought. Specific questions that the runs are designed to address include: What are the mechanisms that maintain drought across the seasonal cycle and from one year to the next? What is the role of the leading patterns of SST variability, and what are the physical mechanisms linking the remote SST forcing to regional drought, including the role of land-atmosphere coupling? The runs were carried out with five different atmospheric general circulation models (AGCM5), and one coupled atmosphere-ocean model in which the model was continuously nudged to the imposed SST forcing. This paper provides an overview of the experiments and some initial results focusing on the responses to the leading patterns of annual mean SST variability consisting of a Pacific El Nino/Southern Oscillation (ENSO)-like pattern, a pattern that resembles the Atlantic Multi-decadal Oscillation (AMO), and a global trend pattern. One of the key findings is that all the AGCMs produce broadly similar (though different in detail) precipitation responses to the Pacific forcing pattern, with a cold Pacific leading to reduced precipitation and a warm Pacific leading to enhanced precipitation over most of the United States. While the response to the Atlantic pattern is less robust, there is general agreement among the models that the largest precipitation response over the U.S. tends to occur when the two oceans have anomalies of opposite sign. That is, a cold Pacific and warm Atlantic tend to produce the largest precipitation reductions, whereas a warm Pacific and cold Atlantic tend to produce the greatest precipitation enhancements. Further analysis of the response over the U.S. to the Pacific forcing highlights a number of noteworthy and to some extent unexpected results. These include a seasonal dependence of the precipitation response that is characterized by signal-to-noise ratios that peak in spring, and surface temperature signal-to-noise ratios that are both lower and show less agreement among the models than those found for the precipitation response. Another interesting result concerns what appears to be a substantially different character in the surface temperature response over the U.S. to the Pacific forcing by the only model examined here that was developed for use in numerical weather prediction. The response to the positive SST trend forcing pattern is an overall surface warming over the world's land areas with substantial regional variations that are in part reproduced in runs forced with a globally uniform SST trend forcing. The precipitation response to the trend forcing is weak in all the models.

This report critically reviews current knowledge about global distributions and properties of atmospheric aerosols, as they relate to aerosol impacts on climate. It assesses possible next steps aimed at substantially reducing uncertainties in aerosol radiative forcing estimates. Current measurement techniques and modeling approaches are summarized, providing context. As a part of the Synthesis and Assessment Product in the Climate Change Science Program, this assessment builds upon recent related assessments, including the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR4, 2007) and other Climate Change Science Program reports. The objectives of this report are (1) to promote a consensus about the knowledge base for climate change decision support, and (2) to provide a synthesis and integration of the current knowledge of the climate-relevant impacts of anthropogenic aerosols for policy makers, policy analysts, and general public, both within and outside the U.S government and worldwide.

The simulation of major midwinter stratospheric sudden warmings (SSWs) in six stratosphere-resolving general circulation models (GCMs) is examined. The GCMs are compared to a new climatology of SSWs, based on the dynamical characteristics of the events. First, the number, type, and temporal distribution of SSW events are evaluated. Most of the models show a lower frequency of SSW events than the climatology, which has a mean frequency of 6.0 SSWs per decade. Statistical tests show that three of the six models produce significantly fewer SSWs than the climatology, between 1.0 and 2.6 SSWs per decade. Second, four process-based diagnostics are calculated for all of the SSW events in each model. It is found that SSWs in the GCMs compare favorably with dynamical benchmarks for SSW established in the first part of the study. These results indicate that GCMs are capable of quite accurately simulating the dynamics required to produce SSWs, but with lower frequency than the climatology. Further dynamical diagnostics hint that, in at least one case, this is due to a lack of meridional heat flux in the lower stratosphere. Even though the SSWs simulated by most GCMs are dynamically realistic when compared to the NCEP-NCAR reanalysis, the reasons for the relative paucity of SSWs in GCMs remains an important and open question.

The current state of our understanding of solar effects on climate is reviewed. As an example of the relevant issues, the climate during the Maunder Minimum is compared with current conditions in GCM simulations that include a full stratosphere and parameterized ozone response to solar spectral irradiance variability and trace gas changes. The GISS Global Climate/Middle Atmosphere Model coupled to a q-flux/mixed layer model is used for the simulations, which begin in 1500 and extend to the present. Experiments were made to investigate the effect of total versus spectrally-varying solar irradiance changes; spectrally-varying solar irradiance changes on the stratospheric ozone/climate response with both pre-industrial and present trace gases; and the impact on climate and stratospheric ozone of the preindustrial trace gases and aerosols by themselves. The results showed that: (1) the Maunder Minimum cooling relative to today was primarily associated with reduced anthropogenic radiative forcing, although the solar reduction added 40% to the overall cooling. There is no obvious distinguishing surface climate pattern between the two forcings. (2)The global and tropical response was greater than 1 C, in a model with a sensitivity of 1.2 C per W m-2. To reproduce recent low-end estimates would require a sensitivity 1/4 as large. (3) The global surface temperature change was similar when using the total and spectral irradiance prescriptions, although the tropical response was somewhat greater with the former, and the stratospheric response greater with the latter. (4) Most experiments produce a relative negative phase of the NAO/AO during the Maunder Minimum, with both solar and anthropogenic forcing equally capable, associated with the tropical cooling and relative poleward EP flux refraction. (5) A full stratosphere appeared to be necessary for the negative AO/NAO phase, as was the case with this model for global warming experiments, unless the cooling was very large, while the ozone response played a minor role and did not influence surface temperature significantly. (6) Stratospheric ozone was most affected by the difference between present day and preindustrial atmospheric composition and chemistry, with increases in the upper and lower stratosphere during the Maunder Minimum. While the estimated UV reduction led to ozone decreases, this was generally less important than the anthropogenic effect except in the upper middle stratosphere, as judged by two different ozone photochemistry schemes. (7) The effect of the reduced solar irradiance on stratospheric ozone and on climate was similar in Maunder Minimum and current atmospheric conditions.

How the elements of the global hydrologic cycle may respond to climate change is reviewed, first from a discussion of the physical sensitivity of these elements to changes in temperature, and then from a comparison of observations of hydrologic changes over the past 100 million years. Observations of current changes in the hydrologic cycle are then compared with projected future changes given the prospect of global warming. It is shown that some of the projections come close to matching the estimated hydrologic changes that occurred long ago when the earth was very warm.

We examine the climate response to solar irradiance changes between the late 17th century Maunder Minimum and the late 18th century. Global average temperature changes are small (about 0.3 to 0.4 C) in both a climate model and empirical reconstructions. However, regional temperature changes are quite large. In the model, these occur primarily through a forced shift toward the low index state of the Arctic Oscillation/North Atlantic Oscillation. This leads to colder temperatures over the Northern Hemisphere continents, especially in winter (1 to 2 C), in agreement with historical records and proxy data for surface temperatures.

The Goddard Institute for Space Studies (GISS) climate/middle atmosphere model has been used to study the impacts of increasing greenhouse gases, polar ozone depletion, volcanic eruptions, and solar cycle variability. We focus on the projection of the induced responses onto Northern Hemisphere winter surface climate. Changes in the model's surface climate take place largely through enhancement of existing variability patterns, with greenhouse gases, polar ozone depletion and volcanic eruptions primarily affecting the Arctic Oscillation (AO) pattern. Perturbations descend from the stratosphere to the surface in the model by altering the propagation of planetary waves coming up from the surface, in accord with observational evidence. Models lacking realistic stratospheric dynamics fail to capture these wave flux changes. The results support the conclusion that the stratosphere plays a crucial role in recent AO trends. We show that in our climate model, while ozone depletion has a significant effect, greenhouse gas forcing is the only one capable of causing the large, sustained increase in the AO observed over recent decades. This suggests that the AO trend, and a concurrent strengthening of the stratospheric vortex over the Arctic, are very likely anthropogenic in origin.

A tropospheric chemistry module has been developed for use within the Goddard Institute for Space Studies (GISS) general circulation model (GCM) to study interactions between chemistry and climate change. The model uses a simplified chemistry scheme based on CO-NOx-CH4 chemistry, and also includes a parameterization for emissions of isoprene, the most important non-methane hydrocarbon. The model reproduces present day annual cycles and mean distributions of key trace gases fairly well, based on extensive comparisons with available observations. Examining the simulated change between present day and pre-industrial conditions, we find that the model has a similar response to that seen in other simulations. It shows a 45% increase in the global tropospheric ozone burden, within the 25% - 57% range seen in other studies. Annual average zonal mean ozone increases by more than 125% at Northern Hemisphere middle latitudes near the surface. Comparison of model runs that allow the calculated ozone to interact with the GCM's radiation and meteorology with those that do not shows only minor differences for ozone. The common usage of ozone fields that are not calculated interactively seems to be adequate to simulate both the present day and the pre-industrial ozone distributions. However, use of coupled chemistry does alter the change in tropospheric oxidation capacity, enlarging the overall decrease in OH concentrations from the pre-industrial to the present by about 10% (-5.3% global annual average in uncoupled mode, -5.9% in coupled mode). This indicates that there may be systematic biases in the simulation of the pre-industrial to present day decrease in the oxidation capacity of the troposphere (though a 10% difference is well within the total uncertainty). Global annual average radiative forcing from pre-industrial to present day ozone change is 0.32 W/sq m. The forcing seems to be increased by about 10% when the chemistry is coupled to the GCM. Forcing values greater than 0.8 W/sq m are seen over large areas of the United States, Southern Europe, North Africa, the Middle East, Central Asia, and the Arctic. Radiative forcing is greater than 1.5 W/sq m over parts of these areas during Northern summer Though there are local differences, the radiative forcing is overall in good agreement with the results of other modeling studies in both its magnitude and spatial distribution, demonstrating that the simplified chemistry is adequate for climate studies.

The GISS Global Climate Middle Atmosphere Model (GCMAM) is used to investigate the effect of estimated solar irradiance changes on climate for the past 500 years. This model is employed to allow the impact of UV variations on the stratosphere to affect the troposphere via wave-mean flow interactions. Multiple experiments are done with only a total solar irradiance change (peaking at 0.2 percent from the Maunder Minimum to today); with estimated spectrally-varying irradiance changes (i.e., peak changes of 0.7 percent in the UV, 0.2 percent in the visible and near IR; and 0.07 percent in the IR greater than 1 micron); and the spectrally-varying changes in conjunction with model calculated ozone responses in the stratosphere. Results of the varying temperature patterns and radiation response will be discussed. Of interest is whether the different methods of forcing the solar-induced climate change produce different spatial surface temperature signatures, particularly ones that can be differentiated from greenhouse gas warming. In preliminary tests, spectrally-varying solar forcing with induced ozone changes for solar maximum minus solar minimum conditions results in a temperature signal that is primarily at high latitudes.The high latitude response arises due to solar/ozone-induced alterations in the stratospheric wind field that affect planetary wave propagation from the troposphere, and alter tropospheric advection patterns. In contrast, forcing by total solar irradiance changes produces significant response at low and subtropical latitudes as well, driven by water vapor and cloud feedbacks to the radiative perturbation.

To investigate the effects of changes in the latitudinal temperature gradient and the global mean temperature on dust concentration in the Northern Hemisphere, experiments with the GISS GCM (Goddard Institute for Space Studies General Circulation Model) are performed. The dust concentration over Greenland is calculated from sources in central and eastern Asia, which are integrated on-line in the model. The results show that an increase in the latitudinal temperature gradient increases both the Asian dust source strength and the concentration over Greenland. The source increase is the result of increased surface winds, and to a minor extent, the increase in Greenland dust is also associated with increased northward transport. Cooling the climate in addition to this increased gradient leads to a decrease in precipitation scavenging, which helps produce a further (slight) increase in Greenland dust in this experiment. Reducing the latitudinal gradient reduces the surface wind and hence the dust source, with a subsequent reduction in Greenland dust concentrations. Warming the climate in addition to this reduced gradient leads to a further reduction in Greenland dust due to enhanced precipitation scavenging. These results can be used to evaluate the relationship of Greenland ice core temperature changes to changes in the latitudinal and global temperatures.

To investigate the effects of changes in the latitudinal temperature gradient and the global mean temperature on dust concentration in the Northern Hemisphere, experiments with the Goddard Institute for Space Studies General Circulation Model (GISS GCM) are performed. The dust concentration over Greenland is calculated from sources in central and eastern Asia, which are integrated on-line in the model. The results show that an increase in the latitudinal temperature gradient increases both the Asian dust source strength and the concentration over Greenland. The source increase is the result of increased surface winds, and to a minor extent, the increase in Greenland dust is also associated with increased northward transport. Cooling the climate in addition to this increased gradient leads to a decrease in precipitation scavenging, which helps produce a further (slight) increase in Greenland dust in this experiment. Reducing the latitudinal gradient reduces the surface wind and hence the dust source, with a subsequent reduction in Greenland dust concentrations. Warming the climate in addition to this reduced gradient leads to a further reduction in Greenland dust due to enhanced precipitation scavenging. These results can be used to evaluate the relationship of Greenland ice core temperature changes to changes in the latitudinal and global temperatures.

The possibility of North Atlantic Deep Water (NADW) changes in both past and future climates has raised the issue of how the Southern Ocean would respond. Recent experiments with the GISS coupled atmosphere-ocean model have shown that a "bipolar see-saw" between NADW production and Antarctic Bottom Water (AABW) production in the Weddell Sea can occur in conjunction with freshening of the North Atlantic. However, this effect operates not through a slow ocean response but via a rapid atmospheric mechanism. As NADW reduces, colder temperatures in the North Atlantic, and Northern Hemisphere in general, are associated with higher surface pressure (increased atmospheric mass). Reduced mass in the Southern Hemisphere occurs in response, with lower pressure over the South Pole (an EOF #1 effect, the "high phase" of the Antarctic Oscillation).The lower pressure is associated with stronger west winds that generate an intensified Antarctic Circumpolar Current (ACC), which leads to longitudinal heat divergence in the South Atlantic (and heat convergence in the Southern Indian Ocean). Colder temperatures in the Weddell Sea region lead to sea ice growth, increased salinity and surface water density, and greater Weddell Sea Bottom Water production. Increased poleward transport of heat occurs in the South Atlantic in conjunction with increased bottom water production, but its convergence at high latitudes is not sufficient to offset the longitudinal heat divergence due to the intensified ACC. The colder temperatures at high latitudes in the South Atlantic increase the latitudinal temperature gradient, baroclinic instability, eddy energy and eddy poleward transport of momentum, helping to maintain the lower pressure over the pole in an interactive manner. The heat flux convergence in the Indian Ocean provides a warming tendency in that region, and overall global production of AABW remains unchanged. These results have implications for the interpretation of the ice core records of the last deglaciation, but may also be relevant for changes during the Holocene and perhaps even in response to increased CO2 forcing

The Goddard Institute for Space Studies global climate model (GISS GCM) is used to examine the sensitivity of the simulated climate to sea ice concentration specifications in the type of simulation done in the Atmospheric Modeling Intercomparison Project (AMIP), with specified oceanic boundary conditions. Results show that sea ice concentration uncertainties of +/- 7% can affect simulated regional temperatures by more than 6 C, and biases in sea ice concentrations of +7% and -7% alter simulated annually averaged global surface air temperatures by -0.10 C and +0.17 C, respectively, over those in the control simulation. The resulting 0.27 C difference in simulated annual global surface air temperatures is reduced by a third, to 0.18 C, when considering instead biases of +4% and -4%. More broadly, least-squares fits through the temperature results of 17 simulations with ice concentration input changes ranging from increases of 50% versus the control simulation to decreases of 50% yield a yearly average global impact of 0.0107 C warming for every 1% ice concentration decrease, i.e., 1.07 C warming for the full +50% to -50% range. Regionally and on a monthly average basis, the differences can be far greater, especially in the polar regions, where wintertime contrasts between the +50% and -50% cases can exceed 30 C. However, few statistically significant effects are found outside the polar latitudes, and temperature effects over the non-polar oceans tend to be under 1 C, due in part to the specification of an unvarying annual cycle of sea surface temperatures. The +/- 7% and 14% results provide bounds on the impact (on GISS GCM simulations making use of satellite data) of satellite-derived ice concentration inaccuracies, +/- 7% being the current estimated average accuracy of satellite retrievals and +/- 4% being the anticipated improved average accuracy for upcoming satellite instruments. Results show that the impact on simulated temperatures of imposed ice concentration changes is least in summer, encouragingly the same season in which the satellite accuracies are thought to be worst. Hence the impact of satellite inaccuracies is probably less than the use of an annually averaged satellite inaccuracy would suggest.

MM II defined a series of experiments to better understand and characterize model transport and to assess the realism of this transport by comparison to observations. Measurements from aircraft, balloon, and satellite, not yet available at the time of MM I [Prather and Remsberg, 1993], provide new and stringent constraints on model transport, and address the limits of our transport modeling abilities. Simulations of the idealized tracers the age spectrum, and propagating boundary conditions, and conserved HSCT-like emissions probe the relative roles of different model transport mechanisms, while simulations of SF6 and C02 make the connection to observations. Some of the tracers are related, and transport diagnostics such as the mean age can be derived from more than one of the experiments for comparison to observations. The goals of the transport experiments are: (1) To isolate the effects of transport in models from other processes; (2) To assess model transport for realistic tracers (such as SF6 and C02) for comparison to observations; (3) To use certain idealized tracers to isolate model mechanisms and relationships to atmospheric chemical perturbations; (4) To identify strengths and weaknesses of the treatment of transport processes in the models; (5) To relate evaluated shortcomings to aspects of model formulation. The following section are included:Executive Summary, Introduction, Age Spectrum, Observation, Tropical Transport in Models, Global Mean Age in Models, Source-Transport Covariance, HSCT "ANOY" Tracer Distributions, and Summary and Conclusions.

Results are presented for six simulations of the Goddard Institute for Space Studies (GISS) global atmosphere-ocean model for the years 1950 to 2099. There are two control simulations with constant 1950 atmospheric composition from different initial states, two GHG experiments with observed greenhouse gases up to 1990 and compounded .5% CO2 annual increases thereafter, and two GHG+SO4 experiments with the same varying greenhouse gases plus varying tropospheric sulfate aerosols. Surface air temperature trends in the two GHG experiments are compared between themselves and with the observed temperature record from 1960 and 1998. All comparisons show high positive spatial correlation in the northern hemisphere except in summer when the greenhouse signal is weakest. The GHG+SO4 experiments show weaker correlations. In the southern hemisphere, correlations are either weak or negative which in part are due to the model's unrealistic interannual variability of southern sea ice cover. The model results imply that temperature changes due to forcing by increased greenhouse gases have risen above the level of regional interannual temperature variability in the northern hemisphere over the past 40 years. This period is thus an important test of reliability of coupled climate models.

The 6-year wind archives from the Goddard Institute for Space Studies/Global Climate-Middle Atmosphere Model (GISS/GCMAM) were in- put to the GISS/Harvard/Irvine Chemical Transport Model (G/H/I CTM) to study the seasonal and interannual variability of the budgets and distributions of nitrous oxide (N2O) and trichlorofluoromethane (CCl3F), with the corresponding chemical loss frequencies recycled and boundary conditions kept unchanged from year to year. The effects of ozone feedback and quasi-biennial oscillation (QBO) were not included. However, the role of circulation variation in driving the lifetime variability is investigated. It was found that the global loss rates of these tracers are related to the extratropical planetary wave activity, which drives the tropical upward mass flux. For N2O, a semiannual signal in the loss rate variation is associated with the interhemispheric asymmetry in the upper stratospheric wave activity. For CCl3F, the semiannual signal is weaker, associated with the comparatively uniform wave episodes in the lower stratosphere. The loss rates lag behind the wave activity by about 1-2 months. The interannual variation of the GCM generated winds drives the interannual variation of the annually averaged lifetime. The year-to-year variations of the annually averaged lifetimes can be about 3% for N2O and 4% for CCl3F.

Using GCM (General Circulation Model) simulations of the troposphere-stratosphere circulation in different climate regimes, the sensitivity of the Brewer-Dobson Circulation to climate change will be discussed. Included will be simulations of the stratosphere with increasing/doubled CO2 as done with different versions of the GISS (Goddard Institute for Space Studies) Global Climate Middle Atmosphere Model, with and without ozone response. To put the results in a broader perspective, stratospheric circulation changes from several paleoclimate simulations will also be presented, including the Last Glacial Maximum and the Paleocene, each with several different (plausible) latitudinal sea surface temperature gradients. Hence the climate changes investigated will include strong variations in baroclinic and topographic tropospheric planetary wave and gravity wave forcing, as well as global mean temperature and atmospheric CO2 levels. Results will be analyzed in terms of wave driving via E-P flux convergences and gravity wave effects, and will be shown to be strongly dependent on wave propagation characteristics.

Increased droughts are to be expected in a warmer world, and so are increased floods. A warmer atmosphere can hold more moisture, and evaporate more water from the surface. Thus, when it is not raining, available soil water should be reduced. When it is raining, it could very well rain harder. Most researchers agree then that a warmer world will have greater hydrologic extremes. In addition, there is a basic imbalance that develops as climate warms, between the loss of moisture from the soil by evaporation and replenishment via precipitation. The land has a smaller heat capacity than the ocean, so it should warm faster. Evaporation from the land proceeds at the rate of its warming, while precipitation derives primarily from evaporation at the ocean surface. As the latter is increasing more slowly, in a warmer world, precipitation will not increase as rapidly as evaporation due to the fact that the oceans warm more slowly than the land surface (evaporation over the ocean is slower than over the land). Hence, more droughts are anticipated in a warmer world, but the specific location of such droughts is somewhat uncertain. To address the question of where droughts are likely to occur, one first needs to have a reasonable sense of what the future magnitude of warming will be, and what the latitudinal distribution of warming will be. For example, the greater the warming at high latitudes relative to low latitudes, the more likely there will be increased drought over the U.S. in summer. In contrast, substantial tropical warming could give us El Nino-like precipitation, with intensified flooding along the southern tier of the U.S. All of these conditions are likely to intensify as the global temperature rises.

Using GCM simulations of the troposphere-stratosphere circulation in different climate regimes, the sensitivity of the Brewer-Dobson Circulation to climate change will be discussed. Included will be simulations of the stratosphere with increasing/doubled C02 as done with different versions of the GISS Global Climate Middle Atmosphere Model, with and without ozone response. To put the results in a broader perspective, stratospheric circulation changes from several paleoclimate simulations will also be presented, including the Last Glacial Maximum and the Paleocene, each with several different (plausible) latitudinal sea surface temperature gradients. Hence the climate changes investigated will include strong variations in baroclinic and topographic tropospheric planetary wave and gravity wave forcing, as well as global mean temperature and atmospheric C02 levels. Results will be analyzed in terms of wave driving via E-P flux convergences and gravity wave effects, and will be shown to be strongly dependent on wave propagation characteristics.

The response of water vapor to climate change is investigated through a series of model studies with varying latitudinal temperature gradients, mean temperatures, and ultimately, actual climate change configurations. Questions to be addressed include: what role does varying convection have in water vapor feedback; do Hadley Circulation differences result in differences in water vapor in the upper troposphere; and, does increased eddy energy result in greater eddy vertical transport of water vapor in varying climate regimes?

A key component of climate/chemistry modeling is how to handle the influx into (and egress from) the troposphere. This is especially important when considering tropospheric ozone, and its precursors (e.g., NO(x) from aircraft). A study has been conducted with various GISS models to determine the minimum requirements necessary for producing realistic troposphere-stratosphere exchange. Four on-line tracers are employed: CFC-11 and SF6 for mixing from the troposphere into the stratosphere, Rn222 for vertical mixing within the troposphere, and 14C for mixing from the stratosphere into the troposphere. Four standard models are tested, with varying vertical resolution, gravity wave drag and location of the model top, and additional subsidiary models are employed to examine specific features. The results show that proper vertical transport between the troposphere and stratosphere in the GISS models requires lifting the top of the model considerably out of the stratosphere, and including gravity wave drag in the lower stratosphere. Increased vertical resolution without these aspects does not improve troposphere-stratosphere exchange. The transport appears to be driven largely by the residual circulation within the stratosphere; associated E-P flux convergences require both realistic upward propagating energy from the troposphere, and realistic pass-through possibilities. A 23 layer version with a top at the mesopause and incorporating gravity wave drag appears to have reasonable stratospheric-tropospheric exchange, in terms of both the resulting tracer distributions and atmospheric mass fluxes.

The potential impact of contrails and alterations in the lifetime of background cirrus due to subsonic airplane water and aerosol emissions has been investigated in a set of experiments using the GISS GCM connected to a q-flux ocean. Cirrus clouds at a height of 12-15km, with an optical thickness of 0.33, were input to the model "x" percentage of clear-sky occasions along subsonic aircraft flight paths, where x is varied from .05% to 6%. Two types of experiments were performed: one with the percentage cirrus cloud increase independent of flight density, as long as a certain minimum density was exceeded; the other with the percentage related to the density of fuel expenditure. The overall climate impact was similar with the two approaches, due to the feedbacks of the climate system. Fifty years were run for eight such experiments, with the following conclusions based on the stable results from years 30-50 for each. The experiments show that adding cirrus to the upper troposphere results in a stabilization of the atmosphere, which leads to some decrease in cloud cover at levels below the insertion altitude. Considering then the total effect on upper level cloud cover (above 5 km altitude), the equilibrium global mean temperature response shows that altering high level clouds by 1% changes the global mean temperature by 0.43C. The response is highly linear (linear correlation coefficient of 0.996) for high cloud cover changes between 0. 1% and 5%. The effect is amplified in the Northern Hemisphere, more so with greater cloud cover change. The temperature effect maximizes around 10 km (at greater than 40C warming with a 4.8% increase in upper level clouds), again more so with greater warming. The high cloud cover change shows the flight path influence most clearly with the smallest warming magnitudes; with greater warming, the model feedbacks introduce a strong tropical response. Similarly, the surface temperature response is dominated by the feedbacks, and shows little geographical relationship to the high cloud input. Considering whether these effects would be observable, changing upper level cloud cover by as little as 0.4% produces warming greater than 2 standard deviations in the Microwave Sounding Unit (MSU) channels 4, 2 and 2r, in flight path regions and in the subtropics. Despite the simplified nature of these experiments, the results emphasize the sensitivity of the modeled climate to high level cloud cover changes, and thus the potential ability of aircraft to influence climate by altering clouds in the upper troposphere.

The (GISS) Goddard Institute for Space Studies coupled atmosphere-ocean model is used to investigate the effect of increased atmospheric CO2 by comparing a compounded 1 percent CO2 increase experiment with a control simulation. After 70 years of integration, the global surface air temperature in the 1 percent CO2 experiment is 1.43 C warmer. In spite of this global warming, there are two distinct regions, the northern Atlantic Ocean and the southern Pacific Ocean, where the surface air temperature is up to 4 C cooler. This situation is maintained by two positive feedbacks: a local effect on convection in the South Pacific and a non-local impact on the meridional circulation in the North Atlantic. The poleward transport of latent energy and dry static energy by the atmosphere is greater in the 1 percent CO2 experiment, caused by warming and therefore increased water vapor and greater greenhouse capacity at lower latitudes. The larger atmospheric transports tend to reduce upward vertical fluxes of heat and moisture from the ocean surface at high latitudes, which has the effect of stabilizing the ocean, reducing both convection and the thermohaline circulation. With less convection, less warm water is brought up from below, and with a reduced North Atlantic thermohaline circulation (by 30 percent at time of CO2 doubling), the poleward energy transport by the oceans decreases. The colder water then leads to further reductions in evaporation, decreases of salinity at high latitudes, continued stabilization of the ocean, and maintenance of reduced convection and meridional overturning. Although sea ice decreases globally, it increases in the cooling regions which reduces the overall climate sensitivity; its effect is most pronounced in the Southern Hemisphere. Tropical warming has been observed over the past several decades; if modeling studies such as this and others which have produced similar effects are valid, these processes may already be beginning.

Understanding tropical sensitivity is perhaps the major concern confronting researchers, for both past and future climate change issues. Tropical data has been beset by contradictions, and many techniques applicable to the extratropics are either unavailable or fraught with uncertainty when applied at low latitudes. Paleoclimate data, if interpreted within the context of the latitudinal temperature gradient data they imply, can be used to estimate what happened to tropical temperatures in the past, and provide a first guess for what might happen in the future. The approach is made possible by the modeling result that atmospheric dynamical changes, and the climate impacts they produce, respond primarily to temperature gradient changes. Here we review some "rules" obtained from GCM (General Circulation Model) experiments with different sea surface temperature gradients and different forcing, that can be used to relate paleoclimate reconstructions to the likely temperature gradient changes they suggest.

The ability to forecast seasonal climate is of great practical interest. One of the most obvious benefits would be agriculture, for which various preparations (planting, machinery, irrigation, manpower) would be enabled. The expectation of being able to make such forecasts far enough in advance (on the order of 9 months) hinges on components of the system with the longest persistence or predictability. The mixed results of El Nino forecasts has raised the hope that tropical Pacific sea surface temperatures (SST) fall into this category. For agriculturally-relevant forecasts to be made, and utilized, requires several conditions. The SST in the regions that affect agricultural areas must be forecast successfully, many months in advance. The climate response to such sea surface temperatures must then be ascertained, either through the use of historical empirical studies or models (e.g., GCMS). For practical applications, the agricultural production must be strongly influenced by climate, and farmers on either the local level or through commercial concerns must be able to adjust to using such forecasts. In a continuing series of papers, we will explore each of these components. This article concerns the question of utilizing SST to forecast the climate in several regions of agricultural production. We optimize the possibility of doing so successfully by using observed SST in a hindcast mode (i.e., a perfect forecast), and we also use the globally observed values (rather than just those from the tropical Pacific, for which predictability has been shown). This then is the ideal situation; in subsequent papers we will explore degrading the results by using only tropical Pacific SSTs, and then using only

Multiple realizations of the 1979-1998 time period have been simulated by the Goddard Institute for Space Studies Atmospheric General Circulation Model (GISS AGCM) to explore its responsiveness to accumulated forcings, particularly over sensitive agricultural regions. A microwave radiative transfer postprocessor has produced the AGCM's lower tropospheric, tropospheric and lower stratospheric brightness temperature (Tb) time series for correlations with the various Microwave Sounding Unit (MSU) time series available. MSU maps of monthly means and anomalies were also used to assess the AGCM's mean annual cycle and regional variability. Seven realizations by the AGCM were forced by observed sea surface temperatures (sst) through 1992 to gather rough standard deviations associated with internal model variability. Subsequent runs hindcast January 1979 through April 1998 with an accumulation of forcings: observed ssts, greenhouse gases, stratospheric volcanic aerosols. stratospheric and tropospheric ozone and tropospheric sulfate and black carbon aerosols. The goal of narrowing gaps between AGCM and MSU time series was complicated by MSU time series, by Tb simulation concerns and by unforced climatic variability in the AGCM and in the real world. Lower stratospheric Tb correlations between the AGCM and MSU for 1979-1998 reached as high as 0.91 +/-0.16 globally with sst, greenhouse gases, volcanic aerosol, stratospheric ozone forcings and tropospheric aerosols. Mid-tropospheric Tb correlations reached as high as 0.66 +/-.04 globally and 0.84 +/-.02 in the tropics. Oceanic lower tropospheric Tb correlations similarly reached 0.61 +/-.06 globally and 0.79 +/-.02 in the tropics. Of the sensitive agricultural areas considered, Nordeste in northeastern Brazil was simulated best with mid-tropospheric Tb correlations up to 0.75 +/- .03. The two other agricultural regions, in Africa and in the northern mid-latitudes, suffered from higher levels of non-sst variability. Zimbabwe had a maximum mid-tropospheric correlation of 0.54 +/- 0.11 while the U.S. Cornbelt had only 0.25 +/- .10. Precipitation and surface temperature performance are also examined over these regions. Correlations of MSU and AGCM time series mostly improved with addition of explicit atmospheric forcings in zonal bands but not in agricultural regional bins each encompassing only six AGCM gridcells.

A mesoscale atmospheric model was used to evaluate the impact of subgrid-scale landscape discontinuities on the vertical profiles of resolved temperature, moisture, and moist static energy in the planetary boundary layer (PBL) of general circulation models (GCMs). These profiles were produced with a 3D version of the model (using a horizontal grid resolution of 7.5 km and 13 vertical layers in the PBL) by averaging horizontally the various atmospheric variables over a 180 x 180 sq km domain-about the size of the horizontal domain represented by a single grid element in a GCM. They were compared to corresponding vertical profiles produced with a 1D version of the model, which simulates the PBL, as in a GCM, over a single horizontal grid element. Differences obtained between the horizontally averaged atmospheric variables produced with the 3D simulations and the 1D simulations emphasize the impact of subgrid-scale landscape discontinuities on GCM-resolved variables. Various types of landscape discontinuities, characterized by horizontal contrasts of surface wetness and size of land patches, were simulated under various background-wind conditions. Differences of temperature, specific humidity, and moist static energy as large as 4 K, 6 g/kg, and 10 kJ/kg were obtained in some cases. These differences were not affected significantly by moderate winds but were sensitive to the spatial distribution of surface wetness. These results emphasize the need to parameterize mesoscale processes induced by landscape discontinuities in GCMss.

A comparison is made of the vertical distribution of high-level cloud tops derived from the Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements and from the International Satellite Cloud Climatology Project (ISCCP) for all Julys and Januarys in 1985 to 1990. The results suggest that ISCCP overestimates the pressure of high-level clouds by up to 50-150 mbar, particularly at low latitudes. This is caused by the frequent presence of clouds with diffuse tops (greater than 50% time when cloudy events are observed). The averaged vertical extent of the diffuse top is about 1.5 km. At midlatitudes where the SAGE II and ISCCP cloud top pressure agree best, clouds with distinct tops reach a maximum relative proportion of the total level cloud amount (about 30-40%), and diffuse-topped clouds are reduced to their minimum (30-40%). The ISCCP-defined cloud top pressure should be regarded not as the material physical height of the clouds but as the level which emits the same infrared radiance as observed. SAGE II and ISCCP cloud top pressures agree for clouds with distinct tops. There is also an indication that the cloud top pressures of optically thin clouds not overlying thicker clouds are poorly estimated by ISCCP at middle latitudes. The average vertical extent of these thin clouds is about 2.5 km.

Global high-level clouds identified in Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements for January and July in the period 1985 to 1990 are compared with near-nadir-looking observations from the International Satellite Cloud Climatology Project (ISCCP). Global and zonal mean high-level cloud amounts from the two data sets agree very well, if clouds with layer extinction coefficients of less than 0.008/km at 1.02 micrometers wavelength are removed from the SAGE II results and all detected clouds are interpreted to have an average horizontal size of about 75 km along the 200 km transimission path length of the SAGE II observations. The SAGE II results are much more sensitive to variations of assumed cloud size than to variations of detection threshold. The geographical distribution of cloud fractions shows good agreement, but systematic regional differences also indicate that the average cloud size varies somewhat among different climate regimes. The more sensitive SAGE II results show that about one third of all high-level clouds are missed by ISCCP but that these clouds have very low optical thicknesses (less than 0.1 at 0.6 micrometers wavelength). SAGE II sampling error in monthly zonal cloud fraction is shown to produce no bias, to be less than the intraseasonal natural variability, but to be comparable with the natural variability at longer time scales.

Future climate change could have significant repercussions for lightning-caused wildfires. Two empirical fire models are presented relating the frequency of lightning fires and the area burned by these fires to the effective precipitation and the frequency of thunderstorm activity. One model deals with the seasonal variations in lightning fires, while the second model deals with the interannual variations of lightning fires. These fire models are then used with the Goddard Institute for Space Studies General Circulation Model to investigate possible changes in fire frequency and area burned in a 2 X CO2 climate. In the United States, the annual mean number of lightning fires increases by 44%, while the area burned increases by 78%. On a global scale, the largest increase in lightning fires can be expected in untouched tropical ecosystems where few natural fires occur today.

A general circulation model (GCM) is used to model global lightning distributions and frequencies. Both total and cloud-to-ground lightning frequencies are modeled using parameterizations that relate the depth of convective clouds to lightning frequencies. The model's simulations of lightning distributions in time and space show good agreement with available observations. The model's annual mean climatology shows a global lightning frequency of 77 flashes per second, with cloud-to-ground lightning making up 25% of the total. The maximum lightning activity in the GCM occurs during the Northern Hemisphere summer, with approximately 91% of all lightning occurring over continental and coastal regions.

The Goddard Institute for Space Studies (GISS) general circulation model (GCM) is used to study the possible implications of past and future climate change on global lightning frequencies. Two climate change experiments were conducted: one for a 2 x CO2 climate (representing a 4.2 degs C global warming) and one for a 2% decrease in the solar constant (representing a 5.9 degs C global cooling). The results suggest at 30% increase in global lightning activity for the warmer climate and a 24% decrease in global lightning activity for the colder climate. This implies an approximate 5-6% change in global lightning frequencies for every 1 degs C global warming/cooling. Both intracloud and cloud-to-ground frequencies are modeled, with cloud-to-ground lightning frequencies showing larger sensitivity to climate change than intracloud frequencies. The magnitude of the modeled lightning changes depends on season, location, and even time of day.

Solar variability is examined in search of implications for global change. The topics covered include the following: solar variation modification of global surface temperature; the significance of solar variability with respect to future climate change; and methods of reducing the uncertainty of the potential amplitude of solar variability on longer time scales.

The response of the climate system to a temporally and spatially constant amount of volcanic particles is simulated using a general circulation model (GCM). The optical depth of the aerosols is chosen so as to produce approximately the same amount of forcing as results from doubling the present CO2 content of the atmosphere and from the boundary conditions associated with the peak of the last ice age. The climate changes produced by long-term volcanic aerosol forcing are obtained by differencing this simulation and one made for the present climate with no volcanic aerosol forcing. The simulations indicate that a significant cooling of the troposphere and surface can occur at times of closely spaced multiple sulfur-rich volcanic explosions that span time scales of decades to centuries. The steady-state climate response to volcanic forcing includes a large expansion of sea ice, especially in the Southern Hemisphere; a resultant large increase in surface and planetary albedo at high latitudes; and sizable changes in the annually and zonally averaged air temperature.

A general circulation model was used to investigate the sensitivity of Asian and African climate to prescribed changes in boundary conditions with the objective of identifying the relative importance of individual high-latitude glacial boundary conditions on seasonal climate and providing a physical basis for interpreting the paleoclimate record. The circulation model is described and results are presented. Insolation forcing increased summer Asian monsoon winds, while increased high-latitude ice cover strengthened winter Asian trade winds causing decreased precipitation. These factors had little effect on African climate. Cooler North Atlantic sea surface temperatures enhanced winter trade winds over North Africa, southern Asian climate was relatively unaffected. Reducing Asian orography enhanced Asian winter circulation while decreasing the summer monsoon. These model results suggest that African and southern Asian climate respond differently to separate elements of high-latitude climate variability.

The ratio of intracloud lightning (IC) to cloud-to-ground lightning (CG) in thunderstorms is observed to vary with latitude. The reason frequently given for this behavior is that the height of the freezing level, and hence the negative charge center in thunderstorms, varies with latitude, resulting in higher ratios in the tropics compared with those in midlatitudes. This study shows that this hypothesis could well be incorrect. Analysis of cloud and lightning data indicates that the IC/CG ratio is linked to the thickness of the cold cloud region in thunderstorms (0 C to cloud top), rather than to the freezing level height.

The purpose of this project was to couple a regional (Southern Ocean) ocean/sea ice model to the existing Goddard Institute for Space Science (GISS) atmospheric general circulation model (GCM). This modification recognizes: the relative isolation of the Southern Ocean; the need to account, prognostically, for the significant air/sea/ice interaction through all involved components; and the advantage of translating the atmospheric lower boundary (typically the rapidly changing ocean surface) to a level that is consistent with the physical response times governing the system evolution (that is, to the base of the fast responding ocean surface layer). The deeper ocean beneath this layer varies on time scales several orders of magnitude slower than the atmosphere and surface ocean, and therefore the boundary between the upper and deep ocean represents a more reasonable fixed boundary condition.

The International Satellite Cloud Climatology Project dataset is used to correlate variations of cloud optical thickness and cloud temperature in today's atmosphere. The analysis focuses on low clouds in order to limit the importance of changes in cloud vertical extent, particle size, and water phase. Coherent patterns of change are observed on several time and space scales. On the planetary scale, clouds in colder, higher latitudes are found to be optically thicker than clouds in warmer, lower latitudes. On the seasonal scale, winter clouds are, for the most part, optically thicker than summer clouds. The logarithmic derivative of cloud optical thickness with temperature is used to describe the sign and magnitude of the optical thickness-temperature correlation. The seasonal, latitudinal, and day-to-day variations of this relation are examined for Northern Hemisphere clouds in 1984. In cold continental clouds, optical thickness increases with temperature, consistent with the temperature variation of the adiabatic cloud water content. In warm continental and in almost all maritime clouds, however, optical thickness decreases with temperature.

A simple parameterization has been developed to simulate global lightning distributions. Convective cloud top height is used as the variable in the parameterization, with different formulations for continental and marine thunderstorms. The parameterization has been validated using two lightning data sets: one global and one regional. In both cases the simulated lightning distributions and frequencies are in very good agreement with the observed lightning data. This parameterization could be used for global studies of lightning climatology; the earth's electric circuit; in general circulation models for modeling global lightning activity, atmospheric NO(x) concentrations, and perhaps forest fire distributions for both the present and future climate; and, possibly, even as a short-term forecasting aid.