Engineering Papers⌕ Search

Engineering topics

Russell, Gary L.

Publications and source records attributed to Russell, Gary L..

Symmetric Equations on the Surface of a Sphere as Used by Model GISS:IB

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.

Icosahedral grid↗

Configuration and Assessment of the GISS ModelE2 Contributions to the CMIP5 Archive

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.

ocean model↗

CMIP5 Historical Simulations (1850-2012) with GISS ModelE2

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.

ocean heat content↗

The Role of Long-Lived Greenhouse Gases as Principal LW Control Knob that Governs the Global Surface Temperature for Past and Future Climate Change

The climate system of the Earth is endowed with a moderately strong greenhouse effect that is characterized by non-condensing greenhouse gases (GHGs) that provide the core radiative forcing. Of these, the most important is atmospheric CO2. There is a strong feedback contribution to the greenhouse effect by water vapor and clouds that is unique in the solar system, exceeding the core radiative forcing due to the non-condensing GHGs by a factor of three. The significance of the non-condensing GHGs is that once they have been injected into the atmosphere, they remain there virtually indefinitely because they do not condense and precipitate from the atmosphere, their chemical removal time ranging from decades to millennia. Water vapor and clouds have only a short lifespan, with their distribution determined by the locally prevailing meteorological conditions, subject to Clausius-Clapeyron constraint. Although solar irradiance is the ultimate energy source that powers the terrestrial greenhouse effect, there has been no discernible long-term trend in solar irradiance since precise monitoring began in the late 1970s. This leaves atmospheric CO2 as the effective control knob driving the current global warming trend. Over geological time scales, volcanoes are the principal source of atmospheric CO2, and the weathering of rocks is the principal sink, with the biosphere participating as both a source and a sink. The problem at hand is that human industrial activity is causing atmospheric CO2, to increase by 2 ppm per year, whereas the interglacial rate has been 0.005 ppm per year. This is a geologically unprecedented rate to turn the CO2 climate control knob. This is causing the global warming that threatens the global environment.

Greenhouse effect↗

Fast Atmosphere-Ocean Model Runs with Large Changes in CO2

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.

ocean models↗

Comparison of Mean Climate Trends in the Northern Hemisphere Between N.C.E.P. and Two Atmosphere-Ocean Model Forced Runs

Results are presented for two greenhouse gas experiments of the Goddard Institute for Space Studies Atmosphere-Ocean Model (AOM). The computed trends of surface pressure, surface temperature, 850, 500 and 200 mb geopotential heights and related temperatures of the model for the time frame 1960-2000 are compared to those obtained from the National Centers for Environmental Prediction observations. A spatial correlation analysis and mean value comparison are performed, showing good agreement. A brief general discussion about the statistics of trend detection is presented. The domain of interest is the Northern Hemisphere (NH) because of the higher reliability of both the model results and the observations. The accuracy that this AOM has in describing the observed regional and NH climate trends makes it reliable in forecasting future climate changes.

Lucarini, Valerio↗

Modeling Feedbacks Between Water and Vegetation in the Climate System

Not only is water essential for life on earth, but life itself affects the global hydrologic cycle and consequently the climate of the planet. Whether the global feedbacks between life and the hydrologic cycle tend to stabilize the climate system about some equilibrium level is difficult to assess. We use a global climate model to examine how the presence of vegetation can affect the hydrologic cycle in a particular region. A control for the present climate is compared with a model experiment in which the Sahara Desert is replaced by vegetation in the form of trees and shrubs common to the Sahel region. A second model experiment is designed to identify the separate roles of two different effects of vegetation, namely the modified albedo and the presence of roots that can extract moisture from deeper soil layers. The results show that the presence of vegetation leads to increases in precipitation and soil moisture in western Sahara. In eastern Sahara, the changes are less clear. The increase in soil moisture is greater when the desert albedo is replaced by the vegetation albedo than when both the vegetation albedo and roots are added. The effect of roots is to withdraw water from deeper layers during the dry season. One implication of this study is that the insertion of vegetation into the Sahara modifies the hydrologic cycle so that the vegetation is more likely to persist than initially.

Miller, James R.↗

Projected Impact of Climate Change on the Energy Budget of the Arctic Ocean by a Global Climate Model

The annual energy budget of the Arctic Ocean is characterized by a net heat loss at the air-sea interface that is balanced by oceanic heat transport into the Arctic. The energy loss at the air-sea interface is due to the combined effects of radiative, sensible, and latent heat fluxes. The inflow of heat by the ocean can be divided into two components: the transport of water masses of different temperatures between the Arctic and the Atlantic and Pacific Oceans and the export of sea ice, primarily through Fram Strait. Two 150-year simulations (1950-2099) of a global climate model are used to examine how this balance might change if atmospheric greenhouse gases (GHGs) increase. One is a control simulation for the present climate with constant 1950 atmospheric composition, and the other is a transient experiment with observed GHGs from 1950 to 1990 and 0.5% annual compounded increases of CO2 after 1990. For the present climate the model agrees well with observations of radiative fluxes at the top of the atmosphere, atmospheric advective energy transport into the Arctic, and surface air temperature. It also simulates the seasonal cycle and summer increase of cloud cover and the seasonal cycle of sea-ice cover. In addition, the changes in high-latitude surface air temperature and sea-ice cover in the GHG experiment are consistent with observed changes during the last 40 and 20 years, respectively. Relative to the control, the last 50-year period of the GHG experiment indicates that even though the net annual incident solar radiation at the surface decreases by 4.6 W(per square meters) (because of greater cloud cover and increased cloud optical depth), the absorbed solar radiation increases by 2.8 W(per square meters) (because of less sea ice). Increased cloud cover and warmer air also cause increased downward thermal radiation at the surface so that the net radiation into the ocean increases by 5.0 Wm-2. The annual increase in radiation into the ocean, however, is compensated by larger increases in sensible and latent heat fluxes out of the ocean. Although the net energy loss from the ocean surface increases by 0.8 W (per square meters), this is less than the interannual variability, and the increase may not indicate a long-term trend. The seasonal cycle of heat fluxes is significantly enhanced. The downward surface heat flux increases in summer (maximum 2 of 19 W per square meters or 23% in June) while the upward heat flux increases in winter (maximum of 16 W per square meters or 28% in November). The increased downward flux in summer is due to a combination of increases in absorbed solar and thermal radiation and smaller losses of sensible and latent heat. The increased heat loss in winter is due to increased sensible and latent heat fluxes, which in turn are due to reduced sea-ice cover. On the other hand, the seasonal cycle of surface air temperature is damped, as there is a large increase in winter temperature but little change in summer.

Miller, James R.↗

Comparison of Model and Observed Regional Temperature Changes During the Past 40 Years

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.

Russell, Gary L.↗

Modeling the Effect of Wetlands, Flooding, and Irrigation on River Flow: Application to the Aral Sea

As the world's population continues to increase, additional stress is placed on water resources. This stress, coupled with future uncertainties regarding climate change, makes arid and semi-arid regions particularly vulnerable. One example is the Aral Sea where the freshwater inflow, which is dominated by snowmelt runoff, has decreased significantly since the expansion of intensive irrigation in the 1960s. The purpose of this paper is to use a river routing scheme from a global climate model to examine the flow of the Amu Dar'ya River into the Aral Sea. The river routing scheme is modified to include groundwater flow, flooding, and evaporative losses in the river's wetlands and floodplain, and anthropogenic withdrawals for irrigation. A set of scenarios is designed to test the sensitivity of river flow to the inclusion of these modifications into the river routing scheme. When riverine wetlands and floodplains are present, the river flow is reduced significantly and is similar to the observed flow. In addition the model results show that it is essential to incorporate human diversions to accurately represent the inflow to the Aral Sea, and they also indicate potential management strategies that might be appropriate to maintain a balance between inflow to the Sea and upstream diversions for irrigation.

Ferrari, Michael R.↗

Response to CO2 Transient Increase in the GISS Coupled Model: Regional Coolings in a Warming Climate

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.

Russell, Gary L.↗

Regional Sea Level Changes Projected by the NASA/GISS Atmosphere-Ocean Model

Sea level has been rising for the past century, and inhabitants of the Earth's coastal regions will want to understand and predict future sea level changes. In this study we present results from new simulations of the Goddard Institute for Space Studies (GISS) global atmosphere-ocean model from 1950 to 2099. Model results are compared with observed sea level changes during the past 40 years at 17 coastal stations around the world. Using observed levels of greenhouse gases between 1950 and 1990 and a compounded 0.5% annual increase in Co2 after 1990, model projections show that global sea level measured from 1950 will rise by 61 mm in the year 2000, by 212 mm in 2050, and by 408 mm in 2089. By 2089, two thirds of the global sea level rise will be due to thermal expansion and one third will be due to ocean mass changes. The spatial distribution of sea level rise is different than that projected by rigid lid ocean models.

Russell, Gary L.↗

Projected Impact of Climate Change on the Water and Salt Budgets of the Arctic Ocean by a Global Climate Model

The annual flux of freshwater into the Arctic Ocean by the atmosphere and rivers is balanced by the export of sea ice and oceanic freshwater. Two 150-year simulations of a global climate model are used to examine how this balance might change if atmospheric greenhouse gases (GHGs) increase. Relative to the control, the last 50-year period of the GHG experiment indicates that the total inflow of water from the atmosphere and rivers increases by 10% primarily due to an increase in river discharge, the annual sea-ice export decreases by about half, the oceanic liquid water export increases, salinity decreases, sea-ice cover decreases, and the total mass and sea-surface height of the Arctic Ocean increase. The closed, compact, and multi-phased nature of the hydrologic cycle in the Arctic Ocean makes it an ideal test of water budgets that could be included in model intercomparisons.

Miller, James R.↗

Stable water isotope behavior during the last glacial maximum: A general circulation model analysis

Global water isotope geochemisty during the last glacial maximum (LGM) is simulated with an 8 deg x 10 deg atmospheric general circulation model (GCM). The simulation results suggest that the spatial delta O-18/temperature relationships observed for the present day and LGM climates are very similar. Furthermore, the temporal delta O-18/temperature relationship is similar to the present-day spatial relationship in regions for which the LGM/present-day temperature change is significant. This helps justify the standard practice of applying the latter to the interpretation of paleodata, despite the possible influence of other factors, such as changes in the evaportive sources of precipitation or in the seasonality of precipitation. The model suggests, for example, that temperature shifts inferred from ice core data may differ from the true shifts by only about 30%.

Jouzel, Jean↗

Continental-scale river flow in climate models

The hydrologic cycle is a major part of the global climate system. There is an atmospheric flux of water from the ocean surface to the continents. The cycle is closed by return flow in rivers. In this paper a river routing model is developed to use with grid box climate models for the whole earth. The routing model needs an algorithm for the river mass flow and a river direction file, which has been compiled for 4 deg x 5 deg and 2 deg x 2.5 deg resolutions. River basins are defined by the direction files. The river flow leaving each grid box depends on river and lake mass, downstream distance, and an effective flow speed that depends on topography. As input the routing model uses monthly land source runoff from a 5-yr simulation of the NASA/GISS atmospheric climate model (Hansen et al.). The land source runoff from the 4 deg x 5 deg resolution model is quartered onto a 2 deg x 2.5 deg grid, and the effect of grid resolution is examined. Monthly flow at the mouth of the world's major rivers is compared with observations, and a global error function for river flow is used to evaluate the routing model and its sensitivity to physical parameters. Three basinwide parameters are introduced: the river length weighted by source runoff, the turnover rate, and the basinwide speed. Although the values of these parameters depend on the resolution at which the rivers are defined, the values should converge as the grid resolution becomes finer. When the routing scheme described here is coupled with a climate model's source runoff, it provides the basis for closing the hydrologic cycle in coupled atmosphere-ocean models by realistically allowing water to return to the ocean at the correct location and with the proper magnitude and timing.

Miller, James R.↗

The impact of global warming on river runoff

A global atmospheric model is used to calculate the annual river runoff for 33 of the world's major rivers for the present climate and for a doubled CO2 climate. The model has a horizontal resolution of 4 x 5 deg, but the runoff from each model grid box is quartered and added to the appropriate river drainage basin on a 2 x 2.5 deg resolution. The computed runoff depends on the model's precipitation, evapotranspiration, and soil moisture storage. For the doubled CO2 climate, the runoff increased for 25 of the 33 rivers, and in most cases the increases coincide with increased rainfall within the drainage basins. There were runoff increases in all rivers in high northern latitudes, with a maximum increase of 47 percent. At low latitudes there were both increases and decreases ranging from a 96 increase to a 43 percent decrease. The effect of the simplified model assumptions of land-atmosphere interactions on the results is discussed.

Miller, James R.↗

The global geochemistry of bomb-produced tritium - General circulation model compared to available observations and traditional interpretations

Observational evidence suggests that of the tritium produced during nuclear bomb tests that has already reached the ocean, more than twice as much arrived through vapor impact as through precipitation. In the present study, the Goddard Institute for Space Studies 8 x 10 deg atmospheric general circulation model is used to simulate tritium transport from the upper atmosphere to the ocean. The simulation indicates that tritium delivery to the ocean via vapor impact is about equal to that via precipitation. The model result is relatively insensitive to several imposed changes in tritium source location, in model parameterizations, and in model resolution. Possible reasons for the discrepancy are explored.

Koster, Randal D.↗

Comparison of horizontal difference schemes for the shallow water equations on a sphere

The accuracy of horizontal difference schemes used in the hydrodynamics parts of General Circulation Models are compared by means of numerical experiments for the shallow water equations on a sphere. As expected, the phase lag of moving waves decreases as the order of accuracy of a scheme increases or as the grid resolution increases. Overall, Takano and Wurtele's partial fourth order energy and potential enstrophy conserving scheme on the C grid is most accurate. It is clearly superior to the other schemes for the Rossby-Haurwitz wave number 6 initial conditions for coarse grid resolution.

Russell, Gary L.↗