A. Lacis

Efficient Three-Dimensional Global Models for Climate Studies: Models I and II

Abstract A global atmospheric model is developed with a computational efficiency which allows long-range climate experiments. The model solves the simultaneous equations for conservation of mass, energy and momentum, and the equation of state on a grid. Differencing schemes for the dynamics are based on work of Arakawa; the schemes do not need any viscosity for numerical stability, and can thus yield good results with coarse resolution. Radiation is computed with a semi-implicit spectral integration, including all significant atmospheric gases, aerosols and cloud particles. Cloud cover and vertical distribution are computed. Convection mixes moisture, heat and momentum, with buoyant air allowed to penetrate to a height determined by its buoyancy. Ground temperature calculations include diurnal variation and seasonal heat storage. Ground hydrology incorporates a water-holding capacity appropriate for the root zone of local vegetation. Snow depth is computed. Snow albedo includes effects of snow age and mas...

Present-Day Atmospheric Simulations Using GISS ModelE: Comparison to In Situ, Satellite, and Reanalysis Data

Abstract A full description of the ModelE version of the Goddard Institute for Space Studies (GISS) atmospheric general circulation model (GCM) and results are presented for present-day climate simulations (ca. 1979). This version is a complete rewrite of previous models incorporating numerous improvements in basic physics, the stratospheric circulation, and forcing fields. Notable changes include the following: the model top is now above the stratopause, the number of vertical layers has increased, a new cloud microphysical scheme is used, vegetation biophysics now incorporates a sensitivity to humidity, atmospheric turbulence is calculated over the whole column, and new land snow and lake schemes are introduced. The performance of the model using three configurations with different horizontal and vertical resolutions is compared to quality-controlled in situ data, remotely sensed and reanalysis products. Overall, significant improvements over previous models are seen, particularly in upper-atmosphere te...

Climate forcings in Goddard Institute for Space Studies SI2000 simulations

[1] We define the radiative forcings used in climate simulations with the SI2000 version of the Goddard Institute for Space Studies (GISS) global climate model. These include temporal variations of well-mixed greenhouse gases, stratospheric aerosols, solar irradiance, ozone, stratospheric water vapor, and tropospheric aerosols. Our illustrations focus on the period 1951–2050, but we make the full data sets available for those forcings for which we have earlier data. We illustrate the global response to these forcings for the SI2000 model with specified sea surface temperature and with a simple Q-flux ocean, thus helping to characterize the efficacy of each forcing. The model yields good agreement with observed global temperature change and heat storage in the ocean. This agreement does not yield an improved assessment of climate sensitivity or a confirmation of the net climate forcing because of possible compensations with opposite changes of these quantities. Nevertheless, the results imply that observed global temperature change during the past 50 years is primarily a response to radiative forcings. It is also inferred that the planet is now out of radiation balance by 0.5 to 1 W/m 2 and that additional global warming of about 0.5� C is already ‘‘in the pipeline.’’ INDEX TERMS: 1620 Global Change: Climate dynamics (3309); 1635 Global Change: Oceans (4203); 1650 Global Change: Solar variability;

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.

Nonsphericity of dust‐like tropospheric aerosols: Implications for aerosol remote sensing and climate modeling

T-matrix computations of light scattering by polydispersions of randomly oriented nonspherical aerosols and Mie computations for equivalent spheres are compared. Findings show that even moderate nonsphericity results in suubstantial errors in the retrieved aerosol optical thickness if satellite reflectance measurements are analyzed using Mie theory. On the other hand, the use of Mie theory for nonspherical aerosols produces negligible errors in the computation of albedo and flux related quantities, provided that the aerosol size distribution and optical thickness are known beforehand. The first result can be explained by large nonspherical-spherical differences in scattering phase function, while the second result follows from small nonspherical-spherical differences in single-scattering albedo and asymmetry parameter. No cancellation of errors occurs if one consistently uses Mie theory in the retrieval algorithm and then in computing the albedo for the retrieved aerosol optical thickness.

Forcings and chaos in interannual to decadal climate change

We investigate the roles of climate forcings and chaos (unforced variability) in climate change via ensembles of climate simulations in which we add forcings one by one. The experiments suggest that most interannual climate variability in the period 1979–1996 at middle and high latitudes is chaotic. But observed SST anomalies, which themselves are partly forced and partly chaotic, account for much of the climate variability at low latitudes and a small portion of the variability at high latitudes. Both a natural radiative forcing (volcanic aerosols) and an anthropogenic forcing (ozone depletion) leave clear signatures in the simulated climate change that are identified in observations. Pinatubo aerosols warm the stratosphere and cool the surface globally, causing a tendency for regional surface cooling. Ozone depletion cools the lower stratosphere, troposphere and surface, steepening the temperature lapse rate in the troposphere. Solar irradiance effects are small, but our model is inadequate to fully explore this forcing. Well-mixed anthropogenic greenhouse gases cause a large surface wanning that, over the 17 years, approximately offsets cooling by the other three mechanisms. Thus the net calculated effect of all measured radiative forcings is approximately zero surface temperature trend and zero heat storage in the ocean for the period 1979–1996. Finally, in addition to the four measured radiative forcings, we add an initial (1979) disequilibrium forcing of +0.65 W/m2. This forcing yields a global surface warming of about 0.2°C over 1979–1996, close to observations, and measurable heat storage in the ocean. We argue that the results represent evidence of a planetary radiative imbalance of at least 0.5° W/m2; this disequilibrium presumably represents unrealized wanning due to changes of atmospheric composition prior to 1979. One implication of the disequilibrium forcing is an expectation of new record global temperatures in the next few years. The best opportunity for observational confirmation of the disequilibrium is measurement of ocean temperatures adequate to define heat storage.

The GISS Global Climate-Middle Atmosphere Model. Part I: Model Structure and Climatology

Abstract The GISS global climate model (Hansen et al.) has been extended to include the middle atmosphere up to an altitude of approximately 85 km. The model has the full array of processes used for climate research, i.e., numerical solutions of the primitive equations, calculation of radiative and surface fluxes, a complete hydrologic cycle with convective and cloud cover parameterizations, etc. In addition, a parameterized gravity wave drag formulation has been incorporated, in which gravity-wave momentum fluxes due to flow over topography, wind shear and convection are calculated at each grid box, using theoretical relationships between the grid-scale variables and expected source strengths. The parameterized waves then propagate vertically upward depending on the instantaneous wind and temperature profiles, with waves breaking at levels in which their momentum flux exceed the background saturation value. Radiative damping is also calculated, and the total momentum convergence in each layer is used to ...

Interpretation of Snow-Climate Feedback as Produced by 17 General Circulation Models

Snow feedback is expected to amplify global warming caused by increasing concentrations of atmospheric greenhouse gases. The conventional explanation is that a warmer Earth will have less snow cover, resulting in a darker planet that absorbs more solar radiation. An intercomparison of 17 general circulation models, for which perturbations of sea surface temperature were used as a surrogate climate change, suggests that this explanation is overly simplistic. The results instead indicate that additional amplification or moderation may be caused both by cloud interactions and longwave radiation. One measure of this net effect of snow feedback was found to differ markedly among the 17 climate models, ranging from weak negative feedback in some models to strong positive feedback in others.

Intercomparison of models representing direct shortwave radiative forcing by sulfate aerosols

The importance of aerosols as agents of climate change has recently been highlighted. However, the magnitude of aerosol forcing by scattering of shortwave radiation (direct forcing) is still very uncertain even for the relatively well characterized sulfate aerosol. A potential source of uncertainty is in the model representation of aerosol optical properties and aerosol influences on radiative transfer in the atmosphere. Although radiative transfer methods and codes have been compared in the past, these comparisons have not focused on aerosol forcing (change in net radiative flux at the top of the atmosphere). Here we report results of a project involving 12 groups using 15 models to examine radiative forcing by sulfate aerosol for a wide range of values of particle radius, aerosol optical depth, surface albedo, and solar zenith angle. Among the models that were employed were high and low spectral resolution models incorporating a variety of radiative transfer approximations as well as a line-by-line model. The normalized forcings (forcing per sulfate column burden) obtained with the several radiative transfer models were examined, and the discrepancies were characterized. All models simulate forcings of comparable amplitude and exhibit a similar dependence on input parameters. As expected for a non-light-absorbing aerosol, forcings were negative (cooling influence) except at high surface albedo combined with small solar zenith angle. The relative standard deviation of the zenith-angle-averaged normalized broadband forcing for 15 models was 8% for particle radius near the maximum in this forcing (∼0.2 μm) and at low surface albedo. Somewhat greater model-to-model discrepancies were exhibited at specific solar zenith angles. Still greater discrepancies were exhibited at small particle radii, and much greater discrepancies were exhibited at high surface albedos, at which the forcing changes sign; in these situations, however, the normalized forcing is quite small. Discrepancies among the models arise from inaccuracies in Mie calculations, differing treatment of the angular scattering phase function, differing wavelength and angular resolution, and differing treatment of multiple scattering. These results imply the need for standardized radiative transfer methods tailored to the direct aerosol forcing problem. However, the relatively small spread in these results suggests that the uncertainty in forcing arising from the treatment of radiative forcing of a well-characterized aerosol at well-specified surface albedo is smaller than some of the other sources of uncertainty in estimates of direct forcing by anthropogenic sulfate aerosols and anthropogenic aerosols generally.

Errors induced by the neglect of polarization in radiance calculations for rayleigh-scattering atmospheres

Abstract Although neglecting polarization and replacing the rigorous vector radiative transfer equation by its approximate scalar counterpart has no physical background, it is a widely used simplification when the incident light is unpolarized and only the intensity of the reflected light is to be computed. In this paper we employ accurate vector and scalar multiple-scattering calculations to perform a systematic study of the errors induced by the neglect of polarization in radiance calculations for a homogeneous, plane-parallel Rayleigh-scattering atmosphere (with and without depolarization) above a Lambertian surface. Specifically, we calculate percent errors in the reflected intensity for various directions of light incidence and reflection, optical thickness of the atmosphere, single-scattering albedos, depolarization factors, and surface albedos. The numerical data displayed can be used to decide whether or not the scalar approximation may be employed depending on the parameters of the problem. We show that the errors decrease with increasing depolarization factor and/or increasing surface albedo. For conservative or nearly conservative scattering and small surface albedos, the errors are maximum at optical thicknesses of about 1. The calculated errors may be too large for some practical applications, and, therefore, rigorous vector calculations should be employed whenever possible. However, if approximate scalar calculations are used, we recommend to avoid geometries involving phase angles equal or close to 0° and 90°, where the errors are especially significant. We propose a theoretical explanation of the large vector/scalar differences in the case of Rayleigh scattering. According to this explanation, the differences are caused by the particular structure of the Rayleigh scattering matrix and come from lower-order (except first-order) light scattering paths involving right scattering angles and right-angle rotations of the scattering plane.