Transient simulations of Heinrich stadial 1 with a carbon-isotope enabled model and impact of poleward intensified southern hemisphere westerlies in a global eddy permitting ocean, sea-ice, carbon cycle model v1-0

This dataset is a collection of outputs of numerical simulations performed with LOVECLIM, an Earth System Model and the eddy permitting global MOM5 model coupled to sea-ice model SYS and carbon cycle model WOMBAT. 
Transient experiments performed with LOVECLIM started from a LGM state described in Menviel et al., 2016. The model was then forced with evolving boundary conditions (orbital parameters, ice-sheet topography and albedo) but with prognostic CO2, d13CO2 and D14C.
During the experiment meltwater was added into the North Atlantic to simulate a weakening of NADW during H1.
Four experiments were performed:
LH1 - freshwater addition in the Southern Ocean and weak southern hemispheric westerlies
LH1SO - salt addition in the Southern Ocean, weak southern hemispheric westerlies
LH1SHW - freshwater addition in the Southern Ocean, stronger southern hemispheric westerlies
LH1SOSHW - salt addition in the Southern Ocean, stronger southern hemispheric westerlies

The dataset include AABW and NPIW strength, atmospheric CO2, d13C and D14C as well as oceanic DIC and d13C values for a subset of experiments. d13C and D14C values are the isotopic signatures of oceanic DIC and respectively represent a measure of the ratio of stable isotopes  13C  and radioactive isotope 14C with respect to the standard . Both are reported in parts per thousand. 

LOVECLIM is an acronym made from the names of the five different models that have been coupled to built the Earth system model: LOch-Vecode-Ecbilt-CLio-agIsm Model (LOVECLIM). LOVECLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. 
The atmospheric component is ECBilt2, a T21, 3-level quasi-geostrophic model. The oceanic component is CLIO3, which is made up of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is 3° y 3°, and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetationmodel that simulates 
the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the oceanic ca
rbon cycle is represented in LOCH, a comprehensive model that takes into account both the solubility and biological pumps.
LOVECLIM description is an extract from http://www.academia.edu/12279222/Description_of_the_Earth_system_model_of_intermediate_complexity_LOVECLIM_version_1.2

The experiments performed with MOM5-SYS-WOMBAT (http://cosima.org.au/index.php/models/mom-sis-025-wombat/) include a control run and a simulation with poleward intensified southern hemisphere westerlies.
The outputs from these simulations include dissolved inorganic carbon and CO2 flux out of the ocean. 

Reference: Menviel, L., P. Spence, J. Yu, M. Chamberlain, R. Matear, K. Meissner, M. England, "Abrupt early deglacial atmospheric CO2 increase driven by Southern Hemisphere westerlies", 2017, in review.

This dataset was produced by Laurie Menviel of the ARC Centre of Excellence for Climate Systems Sciences, as part of the center "Mechanisms and attribution of past and future ocean circulation change" research program.
More information: https://www.climatescience.org.au/content/398-research-program-mechanisms-and-attribution-past-and-future-ocean-circulation-change

Files included are representing:
Figure 1:
 LH1SOSHW-Fig1.dat
 LH1SOSHW-Fig2.txt

Figure 2:
Timeseries:
 LH1-Fig2.txt
 LH1SHW-Fig2.txt
 LH1SO-Fig2.txt
 LH1SOSHW-Fig2.txt
 LH1-Fig2_DICALK.dat
 LH1SHW-Fig2_DICALK.dat
 LH1SO-Fig2_DICALK.dat
 LH1SOSHW-Fig2_DICALK.dat

 3D DIC content for 19k:LGM_DIC_19k.nc
 3D DIC content for 15k in LH1:LH1_DIC_15k.nc
 3D DIC content for 15k in LH1-SO-SHW:LH1SOSHW_DIC_15k.nc

Figure 3:
 3D DIC Control run yr 50:  CTR_DIC_2019.nc
 3D DIC poleward intensified winds yr 50:  W_DIC_2019.nc
 CO2 flux Control run yr 50: CTR_FDIC_2019.nc
 CO2 flux poleward intensified winds yr 50:W_FDIC_2019.nc

Figure 4
 3D d13C content for 19k:LGM_d13C_19k.nc
 3D d13C content for 16k in LH1-SO-SHW:LH1SOSHW_d13C_16k.nc


Other relevant information
ARCCSS CMS wiki: http://climate-cms.unsw.wikispaces.net/ARCCSS+published+datasets

File organisation:
/g/data1/ua8/ARCCSS_Data/T1c13/v1-0/netcdf/....   contains the CF - ACDD compliant netcdf output
/g/data1/ua8/ARCCSS_Data/T1c13/v1-0/ascii/....   contains the ascii output

Contact: laurie.menviel@mq.edu.au for any question on the dataset content and provenance
         paola.petrelli@utas.edu.au for questions or issues with file accessibility

Citation: