Instructions for Reproducing - Protocol and Forcing Information

CESM Code Tag

The Large Ensemble Community Project (LENS) and the Last Millennium Ensemble (LME) simulations used a common CESM1.1 code base.

The LENS tag was preserved as a special CESM tag

cesm1_1_2_LENS

The Last Millennium Ensemble (LME) was not maintained as a dedicated CESM tag. The instructions linked below outline how to reproduce the Last Millennium 850AD control and 850AD-2005 control simulations from a CESM1.1 or CESM1.2 release tag.

Instructions

850AD Control
850AD-1849 Transient Last Millennium

Initial Conditions

Initial conditions among ensemble members were varied using round off error introduced through the CAM namelist. OCN/ICE/CLM/RTM initial conditions were held constant.

Setup instructions

The Last Millennium code base was not preserved as a special CESM tag. However, the instructions listed below can be used to reproduce the Last Millennium 850AD control and 850AD-1849 control simulations. We also include instructions for extending the Last Millennium through the 20th Century.


The input datasets needed for these simulations are available on the NCAR glade filesystem:

/glade/p/cesm/cseg/inputdata/

And on NCAR's SVN input database:

https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/

Forcing files 850AD; 850-1849; 1850-2005

The forcing files are available on NCAR's glade

/glade/p/cesm/cseg/inputdata

Many of the filenames listed in the tables below are also linked to the SVN input data repository [https://svn-ccsm-inputdata.cgd.ucar.edu/] or to other data centers. To gain access the NCAR files you may register for a guest SVN account, or use an existing SVN account and password to download these files directly to your machine. Glade directory paths are relative to

/glade/p/cesm/cseg/inputdata/

SOLAR FORCING: The solar forcing uses the reconstruction of Vieira, Solanki, Krivova, and Usoskin (Vieira et al., Astron. Astrophys.,2011) for total solar irradiance with the spectral variations and '11-year' solar cycle from Schmidt et al., (Geosci. Model Dev., 2012). The input datasets are available on NCAR glade filesystem

/glade/p/cesmdata/cseg/inputdata/atm/cam/solar/SOLAR_SPECTRAL_VK_Lean_849-2008_annual_c130909.nc

VOLCANIC FORCING: For the LM volcanic forcing, we adopt the ice core–based index of Gao et al. (J. Geophys. Res, 2008), which is based on volcanic deposition signals from 54 ice core records from both the Arctic and Antarctica. This reconstruction provides estimates of the stratospheric sulfate aerosol loadings from volcanic eruptions, modeled as a function of latitude, altitude, and month. The twentieth-century simulation uses an eruption dataset of similar temporal and spatial resolution as described by Ammann et al.(Geophys. Res. Lett.,2003) and used in CMIP5 twentieth-century simulations

GREENHOUSE GAS FORCING: The LME evolution of the concentrations of the principal long-lived greenhouse gases (GHGs; CO2, CH4, and N2O) are imposed based on estimates from highresolution ice cores in Antarctica [see Schmidt et al. (Geosci Model Dev., 2012) for discussion and references]. Small variations from 850 to about 1800 are related to natural feedbacks in the carbon and nitrogen cycles.

LAND USE/LAND COVER BOUNDARY CONDITIONS: In the LM simulation, we merged the Pongratz et al. (2008) reconstruction of land use, proposed by PMIP3, with that of Hurtt et al. (2009) used in the CCSM4 twentieth-century simulations to give a consistent and smoothly evolving land use change. The datasets are spliced together by scaling the Pongratz et al. (Global Biogeochem. Cycles. 2008) dataset to match the Hurtt et al. (iLeaps Newsletter, 2009) dataset at 1500 for every land model grid point. The only plant functional types (PFTs) that are changed are those for crops and pasture; all other PFTs remain at their 1850 control prescription.

Note: The surface dataset used for the 20th Century extensions (1850-2005) of the LM simulations was modified for members 002-008

  • Ensemble member 001: fsurdat = surfdata_1.9x2.5_simyr1850_c110114.nc
  • Ensemble members 002-008: fsurdat = surfdata_1.9x2.5_simyr1850_c140303.nc

ORBITAL FORCING: The CESM model adjusts orbital position (eccentricity, obliquity and precession) automatically every model year through the coupler namelist variable (orb_iyear) based on Berger et al., (J. Geophys. Res.,1993).

AEROSOL FORCING: Stratospheric aerosols are prescribed in the model as a fixed single-size distribution in the three layers in the lower stratosphere above the tropopause

OZONE FORCING: The ozone forcing is from WACCM, not SPARC. For the 1850 control run ozone was specified as a repeating annual cycle based on an existing WACCM 300-yr control run. For the transient historical period (1850-2005) we use the average of 2 existing WACCM ensemble members and apply a 10-year running mean to each month separately. The input datasets are available on NCAR glade filesystem

/glade/p/cesm/cseg/inputdata/atm/cam/ozone

Table 1: Reference Case

RUN_REFCASE
b.e11.B1850C5CN.f19_g16.0850cntl.001
RUN_REFDATE
0850-01-01
Restart Files
SVN Download


Table 2: 850 Control Forcing Files

Atmosphere

Filename

Setting

/atm/cam/topo/ consistent-topo-fv1.9x2.5_c130424.nc N/A
/atm/cam/volc/ N/A N/A
/atm/cam/solar/ SOLAR_SPECTRAL_VK_Lean_849-2008_annual_c130909.nc FIXED
/atm/cam/ozone/ ozone_1.9x2.5_L66_1850clim_c130408.nc CYCLICAL
/atm/cam/ggas/ ghg_pmip3_850-2007_annual_c100517.v2.nc FIXED

srf_emis_specifier

/atm/cam/chem/trop_mozart_aero/emis/

DMS: aerocom_mam3_dms_surf_2000_c090129.nc
SO2: ar5_mam3_so2_surf_1850_c090726.nc
SOAG: ar5_mam3_soag_1.5_surf_1850_c100217.nc
bc_a1: ar5_mam3_bc_surf_1850_c090726.nc
num_a1: ar5_mam3_num_a1_surf_1850_c090726.nc
num_a2: ar5_mam3_num_a2_surf_1850_c090726.nc
pom_a1: ar5_mam3_oc_surf_1850_c090726.nc
so4_a1: ar5_mam3_so4_a1_surf_1850_c090726.nc
so4_a2: ar5_mam3_so4_a2_surf_1850_c090726.nc
srf_emis_type = 'CYCLICAL'
CYCLICAL

ext_frc_specifier

(/atm/cam/chem/trop_mozart_aero/emis/)

SO2: ar5_mam3_so2_elev_1850_c090726.nc
bc_a1: ar5_mam3_bc_elev_1850_c090726.nc
num_a1: ar5_mam3_num_a1_elev_1850_c090726.nc
num_a2: ar5_mam3_num_a2_elev_1850_c090726.nc
pom_a1: ar5_mam3_oc_elev_1850_c090726.nc
so4_a1: ar5_mam3_so4_a1_elev_1850_c090726.nc
so4_a2: ar5_mam3_so4_a2_elev_1850_c090726.nc
ext_frc_type = 'CYCLICAL'
CYCLICAL


Table 3: 850AD-1849 Transient Forcing Files

Atmosphere

Filename

Setting

/atm/cam/topo/ consistent-topo-fv1.9x2.5_c130424.nc N/A
/atm/cam/volc/ IVI2LoadingLatHeight501-2000_L18_c20100518.nc N/A
/atm/cam/solar/ SOLAR_SPECTRAL_VK_Lean_849-2008_annual_c130909.nc SERIAL
/atm/cam/ozone/ ozone_1.9x2.5_L66_1850clim_c130408.nc CYCLICAL
/atm/cam/ggas/ ghg_pmip3_850-2007_annual_c100517.v2.nc RAMPED

srf_emis_specifier

/atm/cam/chem/trop_mozart_aero/emis/

DMS: aerocom_mam3_dms_surf_2000_c090129.nc
SO2: ar5_mam3_so2_surf_1850_c090726.nc
SOAG: ar5_mam3_soag_1.5_surf_1850_c100217.nc
bc_a1: ar5_mam3_bc_surf_1850_c090726.nc
num_a1: ar5_mam3_num_a1_surf_1850_c090726.nc
num_a2: ar5_mam3_num_a2_surf_1850_c090726.nc
pom_a1: ar5_mam3_oc_surf_1850_c090726.nc
so4_a1: ar5_mam3_so4_a1_surf_1850_c090726.nc
so4_a2: ar5_mam3_so4_a2_surf_1850_c090726.nc
srf_emis_type = 'CYCLICAL'
CYCLICAL

ext_frc_specifier

(/atm/cam/chem/trop_mozart_aero/emis/)

SO2: ar5_mam3_so2_elev_1850_c090726.nc
bc_a1: ar5_mam3_bc_elev_1850_c090726.nc
num_a1: ar5_mam3_num_a1_elev_1850_c090726.nc
num_a2: ar5_mam3_num_a2_elev_1850_c090726.nc
pom_a1: ar5_mam3_oc_elev_1850_c090726.nc
so4_a1: ar5_mam3_so4_a1_elev_1850_c090726.nc
so4_a2: ar5_mam3_so4_a2_elev_1850_c090726.nc
ext_frc_type = 'CYCLICAL'
CYCLICAL


Table 4: 1850-2005 Transient Forcing Files

Atmosphere

Filename

Setting

/atm/cam/topo/ consistent-topo-fv1.9x2.5_c130424.nc N/A
/atm/cam/volc/ CCSM4_volcanic_1850-2008_prototype1.nc N/A
/atm/cam/solar/ SOLAR_SPECTRAL_VK_Lean_849-2008_annual_c130909.nc SERIAL
/atm/cam/ozone/ ozone_1.9x2.5_L66_1849-2006_c130613.nc INTERP
/atm/cam/ggas/ ghg_pmip3_850-2007_annual_c100517.v2.nc RAMPED

srf_emis_specifier

/atm/cam/chem/trop_mozart_aero/emis/

DMS: aerocom_mam3_dms_surf_1849-2006_c090804.nc
SO2: ar5_mam3_so2_surf_1850-2005_c090804.nc.nc
SOAG: ar5_mam3_soag_1.5_surf_1850-2005_c100429.nc
bc_a1: ar5_mam3_bc_surf_1850-2005_c090804.nc
num_a1: ar5_mam3_num_a1_surf_1850-2005_c090804.nc
num_a2: ar5_mam3_num_a2_surf_1850-2005_c090804.nc
pom_a1: ar5_mam3_oc_surf_1850-2005_c090804.nc
so4_a1: ar5_mam3_so4_a1_surf_1850-2005_c090804.nc
so4_a2: ar5_mam3_so4_a2_surf_1850-2005_c090804.nc
srf_emis_type = 'INTERP_MISSING_MONTHS'
INTERP

ext_frc_specifier

(/atm/cam/chem/trop_mozart_aero/emis/)

SO2: ar5_mam3_so2_elev_1850-2005_c090804.nc
bc_a1: ar5_mam3_bc_elev_1850-2005_c090804.nc
num_a1: ar5_mam3_num_a1_elev_1850-2005_c090804.nc
num_a2: ar5_mam3_num_a2_elev_1850-2005_c090804.nc
pom_a1: ar5_mam3_oc_elev_1850-2005_c090804.nc
so4_a1: ar5_mam3_so4_a1_elev_1850-2005_c090804.nc
so4_a2: ar5_mam3_so4_a2_elev_1850-2005_c090804.nc
ext_frc_type = 'INTERP_MISSING_MONTHS'
INTERP


Table 5: Transient Land Forcing Files 850AD-2005 :: fsurdat and fpftdyn

<

Path

Filename (fsurdat)

Time Period

/lnd/clm2/surfdata/ surfdata_1.9x2.5_simyr0850_c130710.nc Transient: 850-1099
/lnd/clm2/surfdata/ surfdata_1.9x2.5_simyr1100_c131015.nc Transient: 1100-1349
/lnd/clm2/surfdata/ surfdata_1.9x2.5_simyr1350_c131018.nc Transient: 1350-1599
/lnd/clm2/surfdata/ surfdata_1.9x2.5_simyr1600_c131018.nc Transient: 1600-1849
/lnd/clm2/surfdata/ surfdata_1.9x2.5_simyr1850_c140303.nc Transient: 1850-2005

Path

Filename (fpftdyn)

Time Period

/lnd/clm2/surfdata/ surfdata.pftdyn_1.9x2.5_hist_simyr0850-1100_c130710.nc Transient: 850-1099
/lnd/clm2/surfdata/ surfdata.pftdyn_1.9x2.5_hist_simyr1100-1350_c131015.nc Transient: 1100-1349
/lnd/clm2/surfdata/ surfdata.pftdyn_1.9x2.5_hist_simyr1350-1600_c131018.nc Transient: 1350-1599
/lnd/clm2/surfdata/ surfdata.pftdyn_1.9x2.5_hist_simyr1600-1850_c131018.nc Transient: 1600-1849
/lnd/clm2/surfdata/ surfdata.pftdyn_1.9x2.5_hist_simyr1850-2005_c140303.nc Transient: 1850-2005


Table 6: Land forcing files: 850 Control :: fsurdat and fpftdyn

Namelist

Path

Filename

fsurdat /lnd/clm2/surfdata/ surfdata_1.9x2.5_simyr0850_c130710.nc
fpftdyn /lnd/clm2/surfdata/ surfdata.pftdyn_1.9x2.5_hist_simyr0850-1100_c130710.nc