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 (/glade/p/cesmdata/cseg/collections/cesm1_1_2_LENS). The Last Millennium Ensemble (LME) also has a special CESM tag (/glade/p/cesmdata/cseg/collections/cesm1_1_2_rel03.LME) which can be used to reproduce the Last Millennium 850AD control and 850AD-2005 control simulations.  The instructions for setting up the LME simulations are linked below. The input datasets are available on Yellowstone (/glade/p/cesm/cseg/inputdata/) and on NCAR's SVN input database (https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/).


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 preserved as a special CESM tag (/glade/p/cesmdata/cseg/collections/cesm1_1_2_rel03.LME). 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 Yellowstone (/glade/p/cesm/cseg/inputdata/) and on NCAR's SVN input database (https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/).

Instructions: 850AD Control
Instructions: 850AD-1849 Transient Last Millennium


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

The forcing files are available on yellowstone under (/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 yellowstone (/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 Yellowstone: /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.n