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.

Dinoflagellate cyst counts and age model from sediment samples of ODP Hole 175-1085A

Sediment samples from ODP Site 1085 were investigated in order to obtain more information on the initiation and development of the Benguela upwelling system during the middle and upper Miocene. In particular, our intent was to establish the causes of the upwelling as well as the response of the upwelling regime to the development of the Antarctic Circumpolar Current. Based on changes in the calcareous dinoflagellate cyst association, we found an initial increase of the dinoflagellate cyst productivity, probably related to the initiation of upwelling about 11.8 Ma ago. Two distinct increases in cyst productivity in conjunction with temperature decreases of the upper water masses reflect upwelling pulses off Namibia and occur at the end of the Miocene cooling events Mi5 (about 11.5 Ma) and Mi6 (about 10.5 Ma). Both cooling events are associated with an ice volume increase in Antarctica and are thought to have led to an increase in southeasterly winds, possibly causing these two upwelling pulses. We demonstrate a decrease in dinoflagellate cyst productivity and enhanced terrigenous input via the Orange River after the Mi5 event. At about 11.1 Ma, the dinoflagellate cyst productivity increases again. The polar cyst species Caracomia arctica occurs here for the first time. This implies an influence of subantarctic mode water and therefore a change in the quality of the upwelling water which allowed the Benguela upwelling to develop into modern conditions. From about 10.4 Ma, C. arctica forms a permanent part of the association, pointing to an establishment of the upwelling regime.

Stable isotope record, Mg/Ca ratios and age model of St. Stephens Quarry, Alabama

In the largest global cooling event of the Cenozoic Era, between 33.8 and 33.5 Myr ago, warm, high-CO2 conditions gave way to the variable 'icehouse' climates that prevail today. Despite intense study, the history of cooling versus ice-sheet growth and sea-level fall reconstructed from oxygen isotope values in marine sediments at the transition has not been resolved. Here, we analyse oxygen isotopes and Mg/Ca ratios of benthic foraminifera, and integrate the results with the stratigraphic record of sea-level change across the Eocene-Oligocene transition from a continental-shelf site at Saint Stephens Quarry, Alabama. Comparisons with deep-sea (Sites 522 (South Atlantic) and 1218 (Pacific)) d18O and Mg/Ca records enable us to reconstruct temperature, ice-volume and sea-level changes across the climate transition. Our records show that the transition occurred in at least three distinct steps, with an increasing influence of ice volume on the oxygen isotope record as the transition progressed. By the early Oligocene, ice sheets were ~25% larger than present. This growth was associated with a relative sea-level decrease of approximately 105 m, which equates to a 67 m eustatic fall.

Stable carbon and oxygen isotope ratios and age model of ODP Hole 182-1127B

The aim of this study was to evaluate the potential of constructing an oxygen and carbon isotope stratigraphy for the late Pleistocene succession from Hole 1127B drilled on the Great Australian Bight. Stable isotope analyses were performed on bulk- and fine-fraction (<38 µm) sediment samples. The oxygen isotope variations are generally smaller in magnitude than expected from global pelagic records. This is most likely due to the neriticly dominated sediment composition. Correlation of the oxygen isotope data with carbonate mineralogy and downhole logging data shows simultaneous variations and trends, which are particularly evident in the mid-Pleistocene sediments. Correlation of the oxygen isotope data with the classic SPECMAP curve is used to evaluate the stratigraphic potential of the Site 1127 sediments. This study indicates that an isotope stratigraphy based on planktonic and benthic foraminifers is needed to fully evaluate the response of cool-water carbonates deposited in a margin setting to global ice-volume fluctuations and, hence, the associated sea level variations.

Static GPS measurements and elevation model in the hinterland of Neumayer

As part of the CryoSat Cal/Val activities and the pre-site survey for an ice core drilling contributing to the International Partnerships in Ice Core Sciences (IPICS), ground based kinematic GPS measurements were conducted in early 2007 in the vicinity of the German overwintering station Neumayer (8.25° W and 70.65° S). The investigated area comprises the regions of the ice tongues Halvfarryggen and Søråsen, which rise from the Ekströmisen to a maximum of about 760 m surface elevation, and have an areal extent of about 100 km x 50 km each. Available digital elevation models (DEMs) from radar altimetry and the Antarctic Digital Database show elevation differences of up to hundreds of meters in this region, which necessitated an accurate survey of the conditions on-site. An improved DEM of the Ekströmisen surroundings is derived by a combination of highly accurate ground based GPS measurements, satellite derived laser altimetry data (ICESat), airborne radar altimetry (ARA), and radio echo sounding (RES). The DEM presented here achieves a vertical accuracy of about 1.3 m and can be used for improved ice dynamic modeling and mass balance studies.

Static GPS measurements and elevation model in the vicinity of the EDML deep-drilling site

Interpretation of ice-core records requires accurate knowledge of the past and present surface topography and stress-strain fields. The European Project for Ice Coring in Antarctica (EPICA) drilling site (0.0684° E and 75.0025° S, 2891.7 m) in Dronning Maud Land, Antarctica, is located in the immediate vicinity of a transient and splitting ice divide. A digital elevation model is determined from the combination of kinematic GPS measurements with the GLAS12 data sets from the ICESat satellite. Based on a network of stakes, surveyed with static GPS, the velocity field around the EDML drilling site is calculated. The annual mean velocity magnitude of 12 survey points amounts to 0.74 m/a. Flow directions mainly vary according to their distance from the ice divide. Surface strain rates are determined from a pentagon-shaped stake network with one center point, close to the drilling site. The strain field is characterised by along flow compression, lateral dilatation, and vertical layer thinning.

Microcharcoal concentration and elongation, and age model of sediment core MD96-2098

Although grassland and savanna occupy only a quarter of the world's vegetation, burning in these ecosystems accounts for roughly half the global carbon emissions from fire. However, the processes that govern changes in grassland burning are poorly understood, particularly on time scales beyond satellite records. We analyzed microcharcoal, sediments, and geochemistry in a high-resolution marine sediment core off Namibia to identify the processes that have controlled biomass burning in southern African grassland ecosystems under large, multimillennial-scale climate changes. Six fire cycles occurred during the past 170,000 y in southern Africa that correspond both in timing and magnitude to the precessional forcing of north-south shifts in the Intertropical Convergence Zone. Contrary to the conventional expectation that fire increases with higher temperatures and increased drought, we found that wetter and cooler climates cause increased burning in the study region, owing to a shift in rainfall amount and seasonality (and thus vegetation flammability). We also show that charcoal morphology (i.e., the particle's length-to-width ratio) can be used to reconstruct changes in fire activity as well as biome shifts over time. Our results provide essential context for understanding current and future grassland-fire dynamics and their associated carbon emissions.