Search NASA⌕ Search

SEARCH · Search NASA

Results for “climate model testing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Model Performance in Simulating the Mid-Holocene Green Sahara

The Green Sahara Period, spanning about 11,500 to 5,000 years ago, offers an opportunity to test the ability of climate models to simulate large-scale changes in northern African climate through the strengthening of the West African Monsoon. In this study, we evaluate the performance of four models in simulating the mid-Holocene (6,000 BP), namely – EC-Earth, iCESM, CCSM4-Toronto, and the GISS ModelE2.1-G. Two scenarios are considered for each model – a standard PMIP scenario simulated with the mid-Holocene orbital parameters and greenhouse gas concentrations with vegetation prescribed to pre-industrial conditions, as well as a Green-Sahara scenario which additionally considers factors such as enhanced vegetation, reduced dust, presence of lakes, and land and soil feedbacks. All mid-Holocene scenarios capture an increase in monsoonal precipitation in northern Africa. However, a comparison of the two mid-Holocene scenarios reveals significantly higher precipitation in northern Africa for all the Green-Sahara scenarios relative to the PMIP scenarios – an observation consistent across all models. Accompanied by a strengthened Saharan Heat Low, these changes in the West African Monsoon are also linked to polar amplification, a stronger Indian Summer Monsoon and alterations to the Walker circulation. Model results are in agreement with pollen-based SAT records, multi-proxy SST records and African lake level records. This comparison indicates that a realistic simulation of the mid-Holocene Green Sahara requires consideration of multiple factors in addition to orbital and greenhouse gas forcings.

Green Sahara Period↗
NASA NTRS2009

The GLACE-2 Experiment

A major motivation for the study of the coupled land-atmosphere system is the idea that soil moisture anomalies may affect future meteorological variables through their effects on future surface energy and water budgets. If true, the accurate initialization of soil moisture in a subseasonal or seasonal forecast system may improve forecast skill, making the forecast products more valuable to society. The Global Land-Atmosphere Coupling Experiment (GLACE-2) project is examining, with a wide variety of models, the degree to which subseasonal (out to two months) precipitation and air temperature forecasts improve through the realistic initialization of soil moisture. For the first time ever, a global consensus should emerge regarding the value of land initialization for forecasts, perhaps motivating national forecast centers to make full use of land moisture initialization in their operations. Participants in GLACE-2 perform two series of forecasts, each consisting of 100 2-month forecast ensembles (10 members per ensemble) covering ten boreal spring and summer start-dates in each of the years 1986-1995. Series 1 utilizes realistic land surface state initialization, provided through a decadal offline simulation using realistic meteorological forcing, as provided by the Global Soil Wetness Project - Phase 2 (GSWP-2), a research activity of the Global Land-Atmosphere System Study (GLASS) of GEWEX. Series 2 is identical to Series 1 in every way except for the fact that it does not benefit from realistic land state initialization. Through the comparison of Series 1 and 2, we isolate the impact of land initialization on the forecasts. Optional extensions to these base runs include forecasts covering additional years, using alternative meteorological forcing for the land initialization. To date, GLACE-2 has garnered participation from eleven modeling groups, covering 13 atmospheric models. Analysis of available results is already well underway, both at the individual modeling institutions and at the NASA Goddard Space Flight Center, which is coordinating the project. Analysis focuses on two elements of the forecast problem: (i) the quantification of model-specific "predictability" (i.e., the degree to which simulated atmospheric chaos will foil a forecast, even under the assumption of"perfect" model physics, initialization data, and validation data) and its decay with time; and (ii) the quantification of forecast skill, determined through a comparison of predicted precipitation and air temperature against observations. Indeed, the specific contribution of land initialization to both these elements is isolated through a comparison of the Series 1 and 2 forecasts. We examine the two elements at four different forecast leads: 1-15, 16-30, 31-45, and 46-60 days. Statistics-based approaches for enhancing skill (essentially using observational statistics to reduce the impact of the models' climatic biases) are also tested. In the present talk, we provide an update of progress in GLACE-2, featuring quantifications of predictability and forecast skill for a number of the participating models and providing a "first look" at the desired consensus view of land impacts on subseasonal forecasts.

Koster, Randal↗
NASA NTRS1997

An Ecosystem Simulation Model for Methane Production and Emission from Wetlands

Previous experimental studies suggest that methane emission from wetland is influenced by multiple interactive pathways of gas production and transport through soil and sediment layers to the atmosphere. The objective of this study is to evaluate a new simulation model of methane production and emission in wetland soils that was developed initially to help identify key processes that regulate methanogenesis and net flux of CH4 to the air, but which is designed ultimately for regional simulation using remotely sensed inputs for land cover characteristics. The foundation for these computer simulations is based on a well-documented model (CASA) of ecosystem production and carbon cycling in the terrestrial blaspheme. Modifications to represent flooded wetland soils and anaerobic decomposition include three new sub-models for: (1) layered soil temperature and water table depth (WTD) as a function of daily climate drivers, (2) CH4 production within the anoxic soil layer as a function of WTD and CO2 production under poorly drained conditions, and (3) CH4 gaseous transport pathways (molecular diffusion, ebullition, and plant vascular transport) as a function of WTD and ecosystem type. The model was applied and tested using climate and ecological data to characterize tundra wetland sites near Fairbanks, Alaska studied previously by Whalen and Reeburgh. Comparison of model predictions to measurements of soil temperature and thaw depth, water-table depth, and CH4 emissions over a two year period suggest that inter-site differences in soil physical conditions and methane fluxes could be reproduced accurately for selected periods. Day-to-day comparison of predicted emissions to measured CH4 flux rates reveals good agreement during the early part of the thaw season, but the model tends to underestimate production of CH4 during the months of July and August in both test years. Important seasonal effects, including that of falling WTD during these periods, are apparently overlooked in the model formulation. Nevertheless, reasonably close agreement was achieved between the model's mean daily and seasonal estimates of CH4 flux and observed emission rates for northern wetland ecosystems. Several features of the model are identified as crucial to more accurate prediction of wetland methane emission, including the capacity to incorporate influences of localized topographic and hydrologic features on site-specific soil temperature and WTD dynamics, and mechanistic simulation of methane emission transport pathways from within the soil profile.

Potter, C. S.↗
NASA NTRSDate not supplied

Introduction to our COSPAR-Sponsored Session: Space-Observation Contributions Supporting Climate Action

• The challenge for us at COSPAR in putting this program together was that we are primarily a scientifically oriented organization. “Climate Action,” as it is often expressed publicly, tends to be the domain of “activists.” As scientists, on the other hand, to retain our credibility – and for the contributions we make to continue carrying the weight of scientific truth – it is critical that we remain even-handed in our analysis and our reporting of results. • COSPAR is the Committee on Space Research, and fortunately, much of what we observe from space bears objectively on the attributes of the Earth System and how it is evolving. We can measure from space changes in the cryosphere, in the oceans, the land surface and the atmosphere, frequently, on a global scale, to a degree that is generally not possible any other way. And we have been making such observations, in increasing detail, scope, and precision, for about 60 years. • These data help us characterize the current state and identify patterns of change in the Earth System that allow us to infer processes and to improve climate models. With modeling we can test and refine our inferences, and make some predictions. The results can then inform policy and decision-makers, the public in general, as well as those taking direct action. • So, for today’s session we have invited 7 15-minute presentations spanning a broad range of climate-related topics, covering observations from space that can inform climate-related decision-making – from sea-ice extent and sea-level rise, to land surface properties and land-based water resources, to interactions between land or water surfaces and the atmosphere, to severe events in the atmosphere itself. • We begin with an overview of the observations contributed by one of the premier space agencies, given by NASA Chief Scientist and Senior Climate Advisor Dr. Kate Calvin, which will be followed by talks on what we are learning about specific attributes of the Earth System from space-based measurements that can inform climate-related policy and action.

Ralph Kahn↗
NASA NTRS1997

Direct Aerosol Radiative Forcing: Calculations and Measurements from the Tropospheric

Radiative forcing is defined as the change in the net (downwelling minus upwelling) radiative flux at a given level in the atmosphere. This net flux is the radiative power density available to drive climatic processes in the earth-atmosphere system below that level. Recent research shows that radiative forcing by aerosol particles is a major source of uncertainty in climate predictions. To reduce those uncertainties, TARFOX was designed to determine direct (cloud-free) radiative forcing by the aerosols in one of the world's major industrial pollution plumes--that flowing from the east coast of the US over the Atlantic Ocean. TARFOX measured a variety of aerosol radiative effects (including direct forcing) while simultaneously measuring the chemical, physical, and optical properties of the aerosol particles causing those effects. The resulting data sets permit a wide variety of tests of the consistency, or closure, among the measurements and the models that link them. Because climate predictions use the same or similar model components, closure tests help to assess and reduce prediction uncertainties. In this work we use the TARFOX-determined aerosol, gas, and surface properties to compute radiative forcing for a variety of aerosol episodes, with inadvisable optical depths ranging from 0.07 to 0.6. We calculate forcing by several techniques with varying degrees of sophistication, in part to test the range of applicability of simplified techniques--which are often the only ones feasible in climate predictions by general circulation models (GCMs). We then compare computed forcing to that determined from: (1) Upwelling and downwelling fluxes (0.3-0.7 mm and 0.7-3.0 mm) measured by radiometers on the UK MRF C-130. and (2) Daily average cloud-free absorbed solar and emitted thermal radiative flux at the top of the atmosphere derived from the AVHRR radiometer on the NOAA- 14 satellite. The calculations and measurements all yield aerosol direct radiative forcing in the range -50 to -190 W sq m per unit inadvisable optical depth. The magnitudes are about 15 to 100 times larger than the global-average direct forcing expected for the global-average sulfate aerosol optical depth of 0.04. The reasons for the larger forcing in TARFOX include the relatively large optical depths and the focus on cloud-free, daytime conditions over the dark ocean surface. These are the conditions that produce the actual major radiative forcing events that contribute to any global-average climate effect. Detailed comparisons of calculated and measured forcings for specific events are used for more refined tests of closure.

Russell, P. B.↗
DOE OSTI2021

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

Mt. Agung eruption as a confirmation of the effect of atmospheric radiative perturbations on climate

The effects of the 1963 Mt. Agung (Bali) volcanic eruption are studied with respect to the response of the climate system to global radiative perturbations. Volcanic aerosols, spread globally by stratospheric winds, cause a perturbation characterized by a reasonably well known forcing function. Since the climatic response to a large eruption should be pronounced, these phenomena might provide a means to test the validity of climate models. The Mt. Agung eruption is described, and the atmospheric thermal response is calculated. The magnitude, sign, and time delay of temperature changes computed with a simple one-dimensional climate model for both the stratosphere and troposphere are in agreement with the observed temperature changes.

Hansen, J. E.↗
NASA NTRS1987

On requirements for a satellite mission to measure tropical rainfall

Tropical rainfall data are crucial in determining the role of tropical latent heating in driving the circulation of the global atmosphere. Also, the data are particularly important for testing the realism of climate models, and their ability to simulate and predict climate accurately on the seasonal time scale. Other scientific issues such as the effects of El Nino on climate could be addressed with a reliable, extended time series of tropical rainfall observations. A passive microwave sensor is planned to provide information on the integrated column precipitation content, its areal distribution, and its intensity. An active microwave sensor (radar) will define the layer depth of the precipitation and provide information about the intensity of rain reaching the surface, the key to determining the latent heat input to the atmosphere. A visible/infrared sensor will provide very high resolution information on cloud coverage, type, and top temperatures and also serve as the link between these data and the long and virtually continuous coverage by the geosynchronous meteorological satellites. The unique combination of sensor wavelengths, coverages, and resolving capabilities together with the low-altitude, non-Sun synchronous orbit provide a sampling capability that should yield monthly precipitation amounts to a reasonable accuracy over a 500- by 500-km grid.

Thiele, Otto W.↗
NASA NTRS2003

Earth System Dynamics: The Determination and Interpretation of the Global Angular Momentum Budget using the Earth Observing System

The objective of this investigation has been to examine the mass and momentum exchange between the atmosphere, oceans, solid Earth, hydrosphere, and cryosphere. The investigation has focused on changes in the Earth's gravity field, its rotation rate, atmospheric and oceanic circulation, global sea level change, ice sheet change, and global ground water circulation observed by contemporary sensors and models. The primary component of the mass exchange is water. The geodetic observables provided by these satellite sensors are used to study the transport of water mass in the hydrological cycle from one component of the Earth to another, and they are also used to evaluate the accuracy of models. As such, the investigation is concerned with the overall global water cycle. This report provides a description of scientific, educational and programmatic activities conducted during the period July 1, 1999 through June 30,2000. Research has continued into measurements of time-varying gravity and its relationship to Earth rotation. Variability of angular momentum and the related excitation of polar motion and Earth rotation have been examined for the atmosphere and oceans at time-scales of weeks to several years. To assess the performance of hydrologic models, we have compared geodetic signals derived from them with those observed by satellites. One key component is the interannual mass variability of the oceans obtained by direct observations from altimetry after removing steric signals. Further studies have been conducted on the steric model to quantify its accuracy at global and basin-scales. The results suggest a significant loss of water mass from the Oceans to the land on time-scales longer than 1-year. These signals are not reproduced in any of the models, which have poorly determined interannual fresh water fluxes. Output from a coupled atmosphere-ocean model testing long-term climate change hypotheses has been compared to simulated errors from the Gravity Recovery and Climate Experiment (GRACE) mission. Results indicate that GRACE will be able to observe seasonal signals at half-wavelengths ranging from 1000 to 10000 km, and may be able to observe secular trends at half- wavelengths of greater than 2000-3000 km for soil moisture and snow depth if they are as large as some of the climate experiments predict.

Source record↗
NASA NTRS2009

Spatial and Temporal Inter-Relationship between Anomalies and Trends of Temperature, Moisture, Cloud Cover and OLR as Observed by AIRS/AMSU on Aqua

AIRS/AMSU is the advanced IR/MW atmospheric sounding system launched on EOS Aqua in May 2002. Products derived from AIRS/AMSU by the AIRS Science Team include surface skin temperature and atmospheric temperature profiled; atmospheric humidity profiles, fractional cloud clover and cloud top pressure, and OLR. Products covering the period September 2002 through the present have been derived from AIRS/AMSU using the AIRS Science Team Version 5 retrieval algorithm. In this paper, we will show results covering the time period September 2006 - November 2008. This time period is marked by a substantial warming trend of Northern Hemisphere Extra-tropical land surface skin temperatures, as well as pronounced El Nino - La Nina episodes. These both influence the spatial and temporal anomaly patterns of atmospheric temperature and moisture profiles, as well as of cloud cover and Clear Sky and All Sky OLR. The relationships between temporal and spatial anomalies of these parameters over this time period, as determined from AIRS/AMSU observations, will be shown with particular emphasis on which contribute significantly to OLR anomalies in each of the tropics and extra-tropics. Results will also be shown to evaluate the anomalies and trends of temperature profiles and OLR as determined from analysis of AIRS/AMSU data. Global and regional trends during the 6 1/3 year time period are not necessarily indicative of what has happened in the past, or what may happen in the future. Nevertheless, the inter-relationships of spatial and temporal anomalies of atmospheric geophysical parameters with those of surface skin temperature are indicative of climate processes, and can be used to test the performance of climate models when driven by changes in surface temperatures.

Susskind, Joel↗
NASA NTRS2008

Spatial and Temporal Inter-Relationships between Anomalies and Trends of Temperature, Moisture, Cloud Cover, and OLR as Observed by AIRS/AMSU on Aqua

AIRS/AMSU is the advanced IR/MW atmospheric sounding system launched on EOS Aqua in May 2002. Products derived from AIRS/AMSU by the AIRS Science Team include surface skin temperature and atmospheric temperature profiles; atmospheric humidity profiles, fractional cloud cover and cloud top pressure, and OLR. Products covering the period September 2002 through the present have been derived from AIRS/AMSU using the AIRS Science Team Version 5 retrieval algorithm. In this paper, we will show results covering the time period September 2006 - November 2008. This time period is marked by a substantial warming trend of Northern Hemisphere Extratropical land surface skin temperatures, as well as pronounced El Nino - La Nina episodes. These both influence the spatial and temporal anomaly patterns of atmospheric temperature and moisture profiles, as well as of cloud cover and Clear sky and All Sky OLR. The relationships between temporal and spatial anomalies of these parameters over this time period, as determined from AIRS/AMSU observations, will be shown, with particular emphasis on which contribute significantly to OLR anomalies in each of the tropics and extra-tropics. Results will also be shown to validate the anomalies and trends of temperature profiles and OLR as determined from analysis of AIRS/AMSU data. Global and regional trends during the 6 1/3 year period are not necessarily indicative of what has happened in the past, or what may happen in the future. Nevertheless, the inter-relationships of spatial and temporal anomalies of atmospheric geophysical parameters with those of surface skin temperature are indicative of climate processes, and can be used to test the performance of climate models when driven by changes in surface temperatures.

Susskind, Joel↗
NASA NTRS1981

Climate model postprocessing

The development of new postprocessing software of the climate modeling group is summarized. Code, test, and perform simulations with global general circulation models are described. The models improve understanding and ability to predict the vagaries of weather and climate. To learn from and utilize the model results, it is necessary to create elaborate postprocessing software to allow analysis of the large volume of data produced. The models produce sigma history tapes. The sigma history records are interpolated to pressure history records, which are written on a pressure history tape. The model results are analyzed on pressure surfaces, with snap shots or time averages.

Abeles, J.↗
NASA NTRS1983

Radiation budget measurement/model interface

This final report includes research results from the period February, 1981 through November, 1982. Two new results combine to form the final portion of this work. They are the work by Hanna (1982) and Stevens to successfully test and demonstrate a low-order spectral climate model and the work by Ciesielski et al. (1983) to combine and test the new radiation budget results from NIMBUS-7 with earlier satellite measurements. Together, the two related activities set the stage for future research on radiation budget measurement/model interfacing. Such combination of results will lead to new applications of satellite data to climate problems. The objectives of this research under the present contract are therefore satisfied. Additional research reported herein includes the compilation and documentation of the radiation budget data set a Colorado State University and the definition of climate-related experiments suggested after lengthy analysis of the satellite radiation budget experiments.

Vonderhaar, T. H.↗
NASA NTRS2014

The Impact of the Atlantic Cold Tongue on West African Monsoon Onset in Regional Model Simulations for 1998-2002

The Atlantic cold tongue (ACT) develops during spring and early summer near the Equator in the Eastern Atlantic Ocean and Gulf of Guinea. The hypothesis that the ACT accelerates the timing of West African monsoon (WAM) onset is tested by comparing two regional climate model (RM3) simulation ensembles. Observed sea surface temperatures (SST) that include the ACT are used to force a control ensemble. An idealized, warm SST perturbation is designed to represent lower boundary forcing without the ACT for the experiment ensemble. Summer simulations forced by observed SST and reanalysis boundary conditions for each of five consecutive years are compared to five parallel runs forced by SST with the warm perturbation. The article summarizes the sequence of events leading to the onset of the WAM in the Sahel region. The representation of WAM onset in RM3 simulations is examined and compared to Tropical Rainfall Measuring Mission (TRMM), Global Precipitation Climatology Project (GPCP) and reanalysis data. The study evaluates the sensitivity of WAM onset indicators to the presence of the ACT by analysing the differences between the two simulation ensembles. Results show that the timing of major rainfall events and therefore theWAM onset in the Sahel are not sensitive to the presence of the ACT. However, the warm SST perturbation does increase downstream rainfall rates over West Africa as a consequence of enhanced specific humidity and enhanced northward moisture flux in the lower troposphere.

Regional model simulation↗
NASA NTRS2014

Reply to Comment by Laprise on 'the Added Value to Global Model Projections of Climate Change by Dynamical Downscaling: a Case Study over the Continental U.S. Using the GISS-ModelE2 and WRF Models'

In his comment, Laprise raises several points that we agree merit consideration. His primary critique is that our study [Racherla et al., 2012] tested the ability of the WRF regional climate model to reproduce historical temperature and precipitation change relative to the driving global climate model (GCM) using only a single simulation rather than an ensemble. He asserts that the observed changes are smaller than the internal variability in the climate system (i.e., not statistically significant) and that thus a single simulation should not necessarily be able to capture the observations. Laprise points out that the statistical signal is reduced for a multi-decadal trend such as the one we analyzed in comparison with mean climatology and cites two studies showing that for particular climate parameters it can take any years for a signal to be discerned over internal variability. He states that The results of theexperiment as designed were strongly influenced by the presence of internal variability and sampling errors,which masked the rather small climate changes that may have occurred as a consequence of changes inforcing during the period considered. While Laprise discusses statistics in general terms at some length, for the actual climate trends examined in our study, he offers no evidence that the forced signal was smallcompared with internal variability. The two studies he cites [de Ela et al., 2013; Maraun, 2013] do not provide convincing evidence as they concern climate variables averaged over different times and areas. One in fact examines extreme precipitation events, which by definition are rare and thus have a lower significance level. We accept the general point that it is important to consider internal variability, and as noted in our paper we agree that an ensemble of simulations is in principle an optimal, though computationally expensive, approach. While we did not present the statistical significance of the observations in our original paper, we have now evaluated those for the regional temperature trends used in our study to evaluate the added value of WRF and thus can analyze data as to the magnitude of the trends with respect to internal variability.

climate↗
NASA NTRS1999

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.↗
NASA NTRS1978

Internally and externally caused climate change

A numerical climate model is used to simulate climate change forced only by random fluctuations of the atmospheric heat transport. This short-term natural variability of the atmosphere is shown to be a possible 'cause' not only of the variability of the annual world average temperature about its mean, but also long-term excursions from the mean. Various external causes of climate change are also tested with the model and the results compared with observations for the past 100 years. Volcanic dust is shown to have been an important cause of climate change, while the effects of sunspot-related solar constant variation and anthropogenic forcing are not evident.

Robock, A.↗