Description: A significant influence of changes in the westerly winds over the Southern Ocean was proposed as a mechanism to explain a large portion of the glacial atmospheric pCO2 drawdown (Toggweiler et al., 2006). However, additional modelling studies with Earth System Models of Intermediate Complexity do not confirm the size and sometimes even the sign of the impact of southern hemispheric winds on the glacial pCO2 as suggested by Toggweiler (Men- viel et al., 2008; Tschumi et al., 2008, d’Orgeville et al., 2010). We here add to this discussion and explore the potential contribution of changes in the latitudinal position of the winds on Southern Ocean physics and the carbon cycle by using a state-of-the-art ocean general circulation model (MITgcm) in a spatial resolution increasing in the Southern Ocean (2◦ longitude; northern hemisphere: 2◦ latitude; southern hemisphere: 2◦cos(α)). We discuss how the change in carbon cycling is related to the upwelling strength and pattern in the Southern Ocean and how they depend on the changing wind fields and/or the sea ice coverage. While the previous studies explored the impact of the westlies starting from present day or pre-industrial back- ground conditions, we here perform simulations from LGM background climate. Ocean surface conditions are for reasons of consistency taken from output of the COSMOS Earth System model for a pre-industrial control and two LGM runs (Zhang et al., in preparation). Additionally, a northwards shift (by 10◦) of the westerly wind belt as proposed by Toggweiler is investigated.

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.