Description: In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) and document the procedures that we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) Project mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo- and preindustrial (PI) time slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene, as simulated with our COSMOS setup and PRISM boundary conditions, is both warmer and wetter in the global mean than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.

Description: The stable water isotopes H218O and HDO are incorporated as passive tracers into the oceanic general circulation model MPI-OM, and a control simulation under present-day climate conditions is analyzed in detail. Both δ18O and δD distributions at the ocean surface and deep ocean are generally consistent with available observations on the large scale. The modelled δD-δ 18O relations in surface waters slightly deviates from the slope of the global meteoric water line in most basins, and a much steeper slope is detected in Arctic Oceans. The simulated deuterium excess of ocean surface waters shows small variations between 80° S and 55° N, and a strong decrease north of 55° N. The model is also able to capture the quasi-linear relationship between δ18O and salinity S, as well as δD and S, as seen in observational data. Both in the model results and observations, the surface δ−S relations show a steeper slope in extra-tropical regions than in tropical regions, which indicates relatively more addition of isotopically depleted water at high latitudes.

Description: Interpreting stable oxygen isotope (δ18O) records from stalagmites is still one of the complex tasks in speleothem research. Here, we present a novel model-based approach, where we force a model describing the processes and modifications of δ18O from rain water to speleothem calcite (Oxygen isotope Drip water and Stalagmite Model – ODSM) with the results of a state-of-the-art atmospheric general circulation model enhanced by explicit isotope diagnostics (ECHAM5-wiso). The approach is neither climate nor cave-specific and allows an integrated assessment of the influence of different varying climate variables, e.g. temperature and precipitation amount, on the isotopic composition of drip water and speleothem calcite. First, we apply and evaluate this new approach under present-day climate conditions using observational data from seven caves from different geographical regions in Europe. Each of these caves provides measured δ18O values of drip water and speleothem calcite to which we compare our simulated isotope values. For six of the seven caves modeled δ18O values of drip water and speleothem calcite are in good agreement with observed values. The mismatch of the remaining caves might be caused by the complexity of the cave system, beyond the parameterizations included in our cave model. We then examine the response of the cave system to mid-Holocene (6000yr before present, 6ka) climate conditions by forcing the ODSM with ECHAM5-wiso results from 6ka simulations. For a set of twelve European caves, we compare the modeled mid-Holocene-to-modern difference in speleothem calcite δ18O to available measurements. We show that the general European changes are simulated well. However, local discrepancies are found, and might be explained either by a too low model resolution, complex local soil-atmosphere interactions affecting evapotranspiration or by cave specific factors such as non-equilibrium fractionation processes. The mid-Holocene experiment pronounces the potential of the presented approach to analyse δ18O variations on a spatially large (regional to global) scale. Modelled as well as measured European δ18O values of stalagmite samples suggest the presence of a strong, positive mode of the North Atlantic Oscillation at 6 ka before present, which is supported by the respective modelled climate parameters.

Description: Modelling the evolution of the Earth system on long timescales requires the knowledge and understanding of driving mechanisms, such as the hydrological cycle. This is dominant in all components of the Earth’s system, such as atmosphere, ocean, land surfaces/vegetation and the cryosphere. Observations and measurements of stable water isotopes in climate archives can help to decipher and reconstruct climate change and its regional variations. For the cryosphere, the δ18O cycle in the current generation of Earth-System-Models is missing and an efficient and accurate tracer transport scheme is required. We describe ISOPOLIS 1.0 a modular semi-Lagrangian transport scheme of second-order accuracy, which is coupled to the polythermal and thermomechanical ice sheet model SICOPOLIS (version 2.9). Model skill is demonstrated by experiments with a simplified ice sheet geometry and by comparisons of simulated ice cores with data from Greenland (GRIP) and Antarctica (Vostok). The presented method is a valuable tool to investigate the transport of any kind of passive tracer inside a polythermal ice sheet and is an important step to model the whole δ18O cycle.

Description: The large-scale boreal winter climatic patterns associated with interannual variability in a coral oxygen isotope (δ18O) record from the southern Red Sea (Klein et al. [1997]) covering most of the last century are investigated. From the early 1930’s to the early 1960’s, the winter coral δ18O record, reflecting temperature and salinity variations in southern Red Sea surface waters, is associated with global (or large scale) sea surface temperature (SST) and 850mb geopotential height (Z850) anomalies which project on the corresponding patterns associated with the El Niño-Southern Oscillation (ENSO). In contrast, since the early 1960’s the winter coral δ18O record is related to a Z850 pattern that reflects the ENSO-independent part of the East Asian Winter Monsoon (EAWM), which includes the Siberian High, the East Asian through and the East Asian upper-tropospheric Jet. Our results indicate a weakening of the ENSO control on interannual temperature/salinity variations in southern Red Sea surface waters in the early 1960’s, due to the warming of the Indian Ocean, and suggest that information about the non-stationarity in the relationship between ENSO and two distinct modes of EAWM can be documented in southern Red Sea coral δ18O records.