Age models, accumulation rates, sand and carbonates at DSDP Leg 74 holes

Accumulation rates for the five sites drilled during Leg 74 of the Glomar Challenger are presented on a common timescale based on calibration of datum levels to paleomagnetic records in Leg 74 sediments for the Paleogene, and a new compilation by Berggren et al. (1985), for the Neogene, and using the seafloor-spreading magnetic anomaly timescale of Kent (1985). We present data on accumulation of total sediment, of foraminifers, of the noncarbonate portion, and of fish teeth that give a history of productivity, winnowing, carbonate dissolution, and nonbiogenic input to what was then a part of the South Atlantic at about 30 deg S.

Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.

Age models and foraminifer isotopes of four ODP Leg 202 sites

We present benthic isotope stratigraphies for Sites 1236, 1237, 1239, and 1241 that span the late Miocene-Pliocene time interval from 6 to 2.4 Ma. Orbitally tuned timescales were generated for Sites 1237 and 1241 by correlating the high-frequency variations in gamma ray attenuation density, percent sand of the carbonate fraction, and benthic d13C to variations in Earth's orbital parameters. The astronomical timescales for Sites 1237 and 1241 are in agreement with the one from Atlantic Site 925/926 (Ocean Drilling Program Leg 154). The comparison of benthic d18O and d13C records from the east Pacific sites and Atlantic Site 925/926 revealed a surprising clarity of the "41-k.y. signal" in d13C records and a remarkably good correlation between their d13C records. This suggests that the late Miocene-Pliocene amplitudes of obliquity-related d13C cycles reflect a magnitude of global response often larger than that provided by obliquity-related d18O cycles. At Site 1237, the orbitally derived ages of Pliocene magnetic reversal boundaries between the base of Réunion and the top of Thvera confirm astronomical datings of the generally accepted ATNTS2004 timescale, except for the top of Kaena and the base of Sidufjall. Our astronomical age for the top of Kaena is about one obliquity cycle older. The base of Sidufjall appears to be about one precession cycle younger. The age models of Sites 1236 and 1239 were established by correlating their benthic d18O and d13C records directly to the orbitally tuned isotope record of Site 1241.

Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.

A 5 million year reconstruction of sea level, temperature, and sea water d18O

Marine sediment records from the Oligocene and Miocene reveal clear 400,000-year (400-kyr) climate cycles related to variations in orbital eccentricity. These cycles are also observed in the Plio-Pleistocene records of the global carbon cycle. However they are absent in the Late Pleistocene ice-age record over the past 1.5 million years. Here, we present a simulation of global ice volume over the past 5 million years with a coupled system of four 3-D ice-sheet models. Our simulation shows that the 400-kyr long eccentricity cycles of Antarctica vary coherently with d13C records during the Pleistocene suggesting that they drive the long-term carbon cycle changes throughout the past 35 million years. The 400-kyr response of Antarctica is eventually suppressed by the dominant 100-kyr glacial cycles of the large ice sheets in the Northern Hemisphere (NH).