The Climate Change Research Section (CCR) is part of the Climate and
Global Dynamics (CGD) Division at the National Center for Atmospheric Research (NCAR) in
Boulder, Colorado.
The Climate Change Research (CCR) Section uses appropriate climate system models to study
the sensitivity and stability of the Earth system to a variety of forcings, including
changes of greenhouse gases, aerosols, solar irradiance, volcanic forcing, land
characteristics, and land use change.
CCR provides a focal point for NCAR and university paleoclimate research and serves as a
resource to the paleoclimatic and climate change research community in the use of the
Community Climate System Model (CCSM).
CCR is involved with DOE Laboratories in developing and using high-performance coupled
climate models.
Paleoclimate
Bette Otto-Bliesner (CCR), Esther Brady (CCR),
and Christine Shields (CCR) and Sang-Ik Shin and Zhengyu Liu (visitors, University of
Wisconsin) used PaleoCSM1 to investigate the coupled climate system in the tropical
Pacific region over the last glacial-interglacial cycle. CSM reproduces recent
estimates, based on alkenones and Mg/Ca ratios, of sea surface temperature (SST) changes
and gradients in the tropical Pacific and predicts weaker El
Carrie Morrill (ASP) and
Otto-Bliesner used this series of PaleoCSM1 simulations to also study the influence of El
Otto-Bliesner and Amy Clement (University of Miami)
compared the
Liu, Brady, and J. Lynch-Stieglitz (LDEO) studied the global upper
ocean evolution in the Holocene due to changed orbital forcing. They found that the
annual mean SST changes in the early to mid-Holocene are forced primarily by the annual
mean insolation forcing with an overall symmetric response of colder equator (<0.5
Otto-Bliesner collaborating with Jonathan Overpeck (University of
Arizona) and Gifford Miller (CU) are using simulations with CCSM2 to study Arctic climate
change during the Last Interglacial (130 ky BP). Proxy data for the Last
Interglacial (LIG) indicates that terrestrial sites warmed rapidly and early (most by 4-6
This figure shows June-July-August surface temperature change, 130 ky BP (pccsm.08b) versus present (pccsm.08), simulated by CCSM2.
Otto-Bliesner and Brady show that the significant improvements in the component models of CCSM2 now result in improved simulation of the "Greening of the Sahara" during the early to mid-Holocene. Mid-Holocene proxies for northern Africa suggest that lake levels were significantly higher and steppe and xerophytic vegetation grew in present-day desert regions. PMIP atmosphere-only GCMs and previous atmosphere-ocean GCMs (including CSM1) underestimated both the northward shift and the magnitude of the precipitation increase required to maintain steppe vegetation. A CCSM2 simulation for 8.5 ky BP gives a more realistic northward expansion of precipitation into the Saharan-Sahel region with both a longer African summer monsoon season and greater rainfall during the monsoon season.
This figure show the annual precipitation change (cm) over Africa for 8.5 ky BP compared to present simulated by CCSM2. The insets give details of the monthly precipitation changes simulated for 8.5 ky BP.
Brady and
Otto-Bliesner investigated the response of circulation in the Mediterranean Sea and the
exchange with the Atlantic Ocean at Gibraltar to orbital and greenhouse gas forcing in
CCSM2 simulations for the mid-Holocene (8.5 ky BP) and Last Interglacial (130 ky
BP). Sediment cores taken from the Eastern Mediterranean show periodic depositions of dark
meter-thick layers of rich organic sediments called sapropels. Compared to the
present-day simulation, results from these paleo-integrations, which have perihelion in
boreal summer, are consistent with the hypothesis that sapropels may have been deposited
because an enhanced African monsoon increased Nile river discharge into the Eastern
Mediterranean weakening deep water formation. Enhanced Nile river discharge may have
increased the nutrient supply enhancing productivity in the Eastern basin and the weaker
deep water formation may have allowed for anoxic benthic conditions that would have led to
better preservation of the organic matter. Although this model lacks the biogeochemistry
component that would be able to model these conditions definitively, the results of Brady
and Otto-Bliesner are consistent with the physical changes hypothesized. A transient
increasing CO2(by 1%/yr) simulation shows a weak increase in the Mediterranean
Overflow Water (MOW) salinity
and transport at the Strait of Gibraltar due to an increase in evaporation and decrease in
precipitation over the basin. This is consistent with the observed trend over the last 50
years.
This figure shows the zonal overturning streamfunction calculated for the Mediterranean Sea. Units are 106m3/s and positive streamlines exhibit clockwise flow. The streamlines show surface Atlantic water flowing eastward, sinking in the eastern basin and returning as Mediterranean overflow water. In the Paleo-simulations, a reversed circulation is found in the Eastern basin.
Climate of the Past 1000 Years
Caspar Ammann (ASP) and Fortunat Joos (University of Bern, Switzerland) completed transient simulations with the coupled PaleoCSM1 covering roughly one millennium of natural climate variability prior to significant anthropogenic influence. In collaboration with David Schimel (TSS), Otto-Bliesner and Robert Tomas (CCR), experiments were designed to test the fidelity of the previously chosen approach of implementing natural external forcing factors in the coupled simulations of the 20th century. A new new multi-ice core based volcanic forcing series for the period 850-2000 AD and the latest atmospheric composition data from Antarctica were used as forcings. Emphasis of the simulations was focused on the magnitude of potential solar irradiance changes on climate. A best guess scenario was run that included a background trend of solar irradiance back to the Maunder Minimum (when sunspots almost entirely vanished for decades) of roughly double the magnitude of the instrumentally confirmed 11-year cycle range. This was compared with a simulation forced by a solar series scaled by a significantly larger trend (~6x). The base-history for generating past irradiance was based on 10-Beryllium isotopes recovered from polar ice cores. The temporal statistics of this series is consistent with commonly employed sunspot based forcing series. The best guess scenario, for reference, is also consistent with the forcing used in previously reported 20th century simulations. Therefore, the long pre-anthropogenic millennium simulations test the coupled models ability to reproduce naturally forced climate variations.
The two simulations (small solar and large solar forcing) are compared to different proxy-based climate reconstructions. Remarkable agreement of the best-guess run (low solar) is found with most climate reconstructions. Although the experiment remains on the cooler side, it stays almost entirely within the uncertainty bounds. Only very large explosive eruptions cause the hemispheric and global surface temperature to drop out of the expected range, a problem commonly encountered in climate simulations. The large solar experiment, on the other hand, is significantly colder and is thus less consistent with the climate data. In the 20th century, the simulation that was quite successful in generating pre-anthropogenic natural variations remains 0.5 to 0.7 degrees Celsius too cool if anthropogenic forcings are held constant at 1870 AD conditions. Yet if the greenhouse gas changes and tropospheric sulfate aerosol (only direct effect) are included, the simulations reproduce observations very closely.
This figure shows global annual average surface temperature anomalies from 1850-2000 AD simulated by PaleoCSM compared instrumental record (black) and Jones et al. reconstruction (yellow). PaleoCSM1 simulations are small solar (red), large solar (blue), and large solar without anthropogenic forcing after 1870 (green).
Given the good match over the last 1000 years including the instrumental period, Ammann then extended the 'best-guess' simulations in collaboration with Warren Washington (CCR)and Gerald Meehl (CCR) to the year 2100 using the IPCC scenario A2 for future anthropogenic forcing. Natural forcings were held constant at 2000 AD conditions.
This figure shows the comparison of this run with a multi-proxy reconstruction and the instrumental record. Based on the good agreement of the low-frequency magnitude in the past confirming, at least to a first order, the climate sensitivity of the model, the projected future warming of +2.5 degrees Celsius to the end of the 21st century can now be seen in a different, more striking perspective.
Under the umbrella of the Weather and Climate Change Assessment Science Initiative, several projects investigated aspects of the Millennium simulations from the standpoint of uncertainty. Eugene Wahl (ESIG) and Ammann tested stationarity of teleconnections to ENSO throughout the millennium. Many paleoclimate reconstruction techniques rely on the assumptions that stationarity is preserved. Results based on the model output largely support this notion, but some exceptions carry potentially important implications. Additionally, Ammann, Otto-Bliesner and Linda Mearns (ESIG) have worked with the UCAR Education and Outreach Office (Roberta Johnson, Linda Carbone, Susan Foster and Doug Haller) to generate the framework for a teachers guide to the new NCAR Climate Discovery exhibit focusing on the Little Ice Age.
Philippe Naveau (University of Colorado), Hee-Seok Oh (University of Alberta, Canada) and Ammann have generated and applied new statistical tools to extract external forcing fingerprints from climate time series. Using discrete wavelet decomposition on decadal to century time scales, the detection was successful in isolating potential solar influence in different multi-proxy reconstructions and simplify the detection in the coupled model simulations. For the very short-lived volcanic effects, a state-space model approach was chosen that allows quantification of the volcanic cooling with associated posterior probability.
Alan Robock (Rutgers University) and Ammann, with support from the global ice-coring community, started to collect the latest high-resolution ice core sulfate data from around the world. A meeting at the Niels Bohr Institute in Denmark launched this effort to generate a new and improved volcanic forcing dataset specifically designed for climate modeling. This record will be made available to the community through the NOAA Paleoclimate website.
Brad Adams and Michael Mann (University of
Virginia) and Ammann have analyzed the effects of past explosive volcanism on the El
Future Climate Change
As part of the DOE/UCAR Cooperative Agreement, Washington and Meehl have carried out the largest set of publicly available climate change simulations with the Parallel Climate Model (PCM) that we are aware of anywhere in the world. The data from these simulations are available from the Program for Climate Model Diagnosis and Intercomparison (PCMDI) at Lawrence Livermore National Laboratory (LLNL), National Energy Research Scientific Computing Center (NERSC) located at Lawrence Berkeley Laboratory, and at NCAR. This archive of simulations includes a 1000-year control run, SRES scenario runs for the 21st century, business as usual, and stabilization of greenhouse gas runs to the end of the 22nd century, special simulations for the demonstration project of the Advanced Climate Prediction Initiative (ACPI), a large number of single and combined forcing runs for 20th and other simulations. Information on the simulations can be found at www.cgd.ucar.edu/pcm/. Gary Strand has been the coordinator of the large data archive. and he has been the lead CGD person involved in the DOE Earth System Grid (ESG) project that will provide web access to PCM and CCSM that are distributed at several supercomputing centers. This archive complements the smaller number of CSM simulations with some of the same external climate forcings. The computer time for these simulations was obtained from the Scientific Computing Division of NCAR, NERSC, Oak Ridge National Laboratory (ORNL), and Los Alamos National Laboratory (LANL).
Higher Resolution Versions of CCSM
Both the CCAWG and the DOE Climate Change Prediction Program (CCPP) are particularly interested in regional aspects of climate change. We have seen that by increasing horizontal resolution with the uncoupled version of CCSM Atmosphere Model (CAM) at T85 (developed by the Climate Modeling Section (CMS) and the CCSM Atmospheric Working Group) that there is substantial improvement in the regional aspects of precipitation especially in mountainous areas. This is particularly important for the western North American regions. Also, there appears to be improvement in the surface winds in the Arctic Ocean region, which in turn will affect the flow of sea ice. The present T42 version does not drive the sea ice correctly in several areas of the Arctic Ocean. During the last year we have performed some limited tests using T85 coupled to the CCSM ocean, sea ice, and land components. The tests indicate several improvements. This is very encouraging since these have been persistent systematic errors in previous lower resolution versions.
PCM Climate Change Simulations of the 20th and 21st Centuries: Sensitivity and Decadal Variability
Gerald Meehl (joint with CCR and CAS) collaborated with Aiguo Dai,
(CCR) Warren Washington, (CCR) Tom Wigley (CAS) and Julie Arblaster (CCR) to run the
model and analyze the influences of various anthropogenic (GHGs, ozone, sulfate aerosols)
and natural (volcanic and solar) forcings over the 20th century in the PCM. Four member
ensembles of the single forcings and an additional eight experiments with combinations of
forcings in four member ensembles produced the largest number of forced 20th century
climate experiments ever run with a global coupled climate model. Early century warming
was shown to be mainly due to solar forcing, with enhanced tropical convection
contributing to amplifying the solar forcing. The global temperature signal from the
forcings was shown to be roughly additive.
Meehl, Washington, Thomas Bettge (SCD), Arblaster and Wigley collaborated with Ben Santer
(LLNL) and a team of scientists to analyze tropopause height as an indicator of climate
change in models and observations, with a warming troposphere associated with increased
tropopause height in the model and observations. Through this work they were able to
quantify the contributions of the combination of natural and anthropogenic forcings on
changes in tropopause height, and confirm that the troposphere actually warmed during the
second half of the 20th century due to anthropogenic forcing factors, mainly GHGs and
ozone.
In collaborative work with Arblaster, Meehl documented the mechanisms involved with the
model result of increased mean south Asian monsoon precipitation and increased interannual
variability in a future warmer climate. Sensitivity
experiments with CCM3 showed that increases in mean monsoon precipitation were caused
mainly by warmer Indian Ocean SSTs, but increased variability was associated with forcing
from the tropical Pacific via the Walker Circulation.
Meehl collaborated with Ping Liu (University of Hawaii) and analyzed results from the
global coupled models in the Coupled Model Intercomparison Project (CMIP) to study
possible future changes of the Sahara Desert. Five of the eighteen CMIP models were chosen
for analysis based on their ability to simlulate a reasonable present-day climatology of
the Sahara Desert with similar rainfall distributions and meridional boundaries as in the
observational data. When CO2 concentration was increased at one percent per
year for 80 years in these models, the Sahara moved north, became hotter, and dried. The
local enhanced greenhouse effect from increased CO2 increased the net surface
sensible heat flux, which in turn contributed to the warming trend.
Meehl also collaborated with Washington and Arblaster to analyze four recent NCAR global
coupled models (PCM, CSM, PCTM and CCSM) to show that the relevant feedbacks for climate
system response to increasing CO2, (ice/albedo. water vapor, and clouds) are
managed by the atmospheric model. The ocean, sea ice and land surface play secondary roles
for the globally averaged response. Two models with identical atmospheres but different
ocean and sea ice components (PCM and PCTM) have the most similar resonse to increasing CO2,
followed closely by the CSM with comparable atmosphere and different ocean and sea ice.
The model that is most dissimilar, the CCSM, has a different response from either of the
other three, and in particular is different from PCTM in spite of very similar ocean and
sea ice but different atmospheric model components. Analysis of ocean heat uptake in these
models showed that processes and
response were model dependent, and must be studied more thoroughly in each model to
understand the global climate sensitivity of the coupled system.
Meehl and Arblaster collaborated with David Karoly (University of Oklahoma) and a team of
scientists to document a set of simple indices for climate change. These indices capture
aspects of regional climate such as land-ocean temperature contrast for North American
climate. The PCM agrees with observations for 20th century climate change for these
indices. This suggests that such indices can be used for climate change/detection to
provide more detail than globally averaged temperature, while still capturing the
essential aspects of regional scale climate change. This can be done without going to
pattern-based techniques that may have problems with data coverage in observations, or
processes at unresolved spatial scales in the current medium resolution climate models.
Aiguo Dai (CCR), with others in CCR, analyzed trends and variability
in the surface air temperature, precipitation, ENSO and Arctic Oscilliation variability,
and Atlantic Ocean circulations in the 1200-year control run of the PCM. The PCM control
run is among the first few millennial simulations by a coupled GCM without flux
adjustments. The PCM shows small drifts in surface temperature and other fields that
are comparable to those of flux-corrected GCMs. The Atlantic thermohaline
circulation (THC) in the PCM control run shows many interesting features, including
a multidecadal (~23 yr) oscillation, a reversal of the THC trends, the influence of the
Antarctic Bottom Water on the THC, and a horizontal current changes in the North Atlantic
associated with the THC changes.
Dai also analyzed a PCM climate change simulation in which CO2 and other
greenhouse gases and sulfate aerosols are projected to increase in the 21st century and
stablize in the 22nd century (CO2 at 550ppmv), and thereafter the greenhouse
gases decrease so that their concentrations are symmetric around year 2200 (i.e., CO2
concentration in 2450 equals that in 1950). The results show that the response of the
Atlantic thermohaline circulation (THC) is fairly linear, i.e., the THC strength decreases
with the CO2 increases and increases with the CO2 decline with very small phase
lag.
Dai (with others in CCR) analyzed the PCM simulations for DOE-sponsored Accelerated
Climate Prediction Initiative (ACPI) Program In these runs, the oceans were initialized to
1995 conditions by a group from the Scripps Institution of Oceanography and other
institutions. An ensemble of model runs was then carried out to the year 2099 using the
projected forcing. It is shown that the initialization to1995 conditions removes a large
part of the un-forced oceanic temperature and salinity drifts that occurred in the
standard 20th century integration.
Using
a coupled climate system model, the Parallel Climate Model (PCM), regional responses of
the THC in the North Atlantic to increased CO2 and the underlying physical
processes were studied by Aixue Hu (CCR). The Atlantic THC shows a 20-year cycle in the
control run, qualitatively consistent with observations. Compared with the
Hu
also performed CMIP coordinated xperiments Using CCSM2.0.
Land Cover Change
Feddema (U. Kansas), J. Dobson
(U. Kansas), and S. Freire (ORNL), have made a great deal of progress on developing data
sets and changes to the models specifying anthropogenic land surface changes from 1750 to
2050 that can be used in PCM and CCSM. This will provide a new forcing for anthropogenic
changes in the climate system.
More information may be found on the Parallel Climate Model and CCSM Climate Change and Assessment Working Group (CCAWG) web pages.
The following are scientists and programming staff involved in the PCM effort or its components in alphabetical order: J. Arblaster (NCAR), T. Bettge (NCAR), L. Buja (NCAR), A. Craig (NCAR), A. Dai (NCAR), J. Dennis (NCAR), J. Dukowicz (LANL), J. Hack (NCAR), A. Hu, (NCAR), E. Hunke (LANL), R. James (NCAR), P. Jones (LANL), R. Loft (NCAR), R. Malone (LANL), M. Maltrud (LANL), W. Maslowski (NPS), G. Meehl (NCAR), A. Middleton (NCAR), A. Semtner (NPS), R. Smith (LANL), Ilana Stern (NCAR), G. Strand (NCAR), Robert Tomas (NCAR), W. Washington (NCAR), V. Wayland (NCAR), D. Williamson (NCAR), and Y. Zhang (NPS).
Climate Change Sensitivity Studies
Jeffrey Kiehl (CCR) in collaboration with Christine Shields (CCR) have carried out an analysis of the equilibrium slab ocean simulations of the CAM2 and CAMX atmospheric versions of the CCSM. These simulations show the climate sensitivity of the new model is 2.6 °K for a doubling of CO2, while the sensitivity of the CAM2 model is 2.2 °K. Analysis of these two equilibrium simulations and simulations of the fully coupled models forced with a 1% increase per year in CO2, indicate that changes in the low cloud response to warming is the major reason for the increase in sensitivity. In CAM2 low clouds were diagnosed as a function of relative humidity for all cloud types including convective clouds. In the new model convective cloud amount is diagnosed as a function of cloud mass flux. Studies from high resolution cloud resolving models have shown that cloud mass flux is a better predictor of convective cloud amount than relative humidity. This change in the formulation of low cloud amount leads to a smaller increase in low cloud amount for a warmer climate state.
This figure shows the time series of global annual change in shortwave cloud forcing (Wm-2) due to a 1% per increase in atmospheric carbon dioxide. The blue curve is for the CCSM2 and the red curve is for the CCSM3.
The more negative the cloud forcing the greater the negative feedback on the climate system. The CCSM2 has a much stronger negative cloud feedback than the CCSM3 model.
This figure shows the correlation between longwave cloud forcing and 500 hPa vertical velocity from the new version of the CAM (red) and observations (black).
These analyses have been carried out for tropical, subtropical and extratropical synoptic regimes. The largest differences between model and observations are in the simulation of high cloud amount and high latitude predicted cloud water.
A comparison of the CAMX with the latest version of the NOAA Geophysical Fluid Dynamics Laboratory (GFDL) atmospheric model has also been carried out. This project falls under the new memorandum of understanding between NCAR and NOAA GFDL. The goal of this collaboration is to understand model sensitivities and the sources of agreement or difference in climate sensitivity. It is hoped that the regional sensitivity analysis will provide a useful means to understand inter-model differences.
Another method of analysis is to force the atmosphere model with sea surface temperature perturbations from coupled CO2 simulations. Two sea surface temperature perturbation patterns have been used to force the atmospheric model. The first pattern comes from the slab ocean 2XCO2 of the CAMX, the second pattern was created by B. Soden (NOAA GFDL) from a composite of CMIP2 1% simulations. Both the CAMX and the GFDL AM2 model have been forced with the CMIP sea surface temperature perturbation data. These simulations are being used to study the response of regional climates to a prescribed identical forcing.