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.

Isotopic signature of specific biomarkers derived from anaerobic methanotrophic archaea (ANME groups) and sulphate-reducing bacteria (SRB) at different cold seep provinces

Anaerobic methane-oxidizing microbial communities in sediments at cold methane seeps are important factors in controlling methane emission to the ocean and atmosphere. Here, we investigated the distribution and carbon isotopic signature of specific biomarkers derived from anaerobic methanotrophic archaea (ANME groups) and sulphate-reducing bacteria (SRB) responsible for the anaerobic oxidation of methane (AOM) at different cold seep provinces of Hydrate Ridge, Cascadia margin. The special focus was on their relation to in situ cell abundances and methane turnover. In general, maxima in biomarker abundances and minima in carbon isotope signatures correlated with maxima in AOM and sulphate reduction as well as with consortium biomass. We found ANME-2a/DSS aggregates associated with high abundances of sn-2,3-di-O-isoprenoidal glycerol ethers (archaeol, sn-2-hydroxyarchaeol) and specific bacterial fatty acids (C16:1omega5c, cyC17:0omega5,6) as well as with high methane fluxes (Beggiatoa site). The low to medium flux site (Calyptogena field) was dominated by ANME-2c/DSS aggregates and contained less of both compound classes but more of AOM-related glycerol dialkyl glycerol tetraethers (GDGTs). ANME-1 archaea dominated deeper sediment horizons at the Calyptogena field where sn-1,2-di-O-alkyl glycerol ethers (DAGEs), archaeol, methyl-branched fatty acids (ai-C15:0, i-C16:0, ai-C17:0), and diagnostic GDGTs were prevailing. AOM-specific bacterial and archaeal biomarkers in these sediment strata generally revealed very similar d13C-values of around -100 per mill. In ANME-2-dominated sediment sections, archaeal biomarkers were even more 13C-depleted (down to -120 per mill), whereas bacterial biomarkers were found to be likewise 13C-depleted as in ANME-1-dominated sediment layers (d13C: -100 per mill). The zero flux site (Acharax field), containing only a few numbers of ANME-2/DSS aggregates, however, provided no specific biomarker pattern. Deeper sediment sections (below 20 cm sediment depth) from Beggiatoa covered areas which included solid layers of methane gas hydrates contained ANME-2/DSS typical biomarkers showing subsurface peaks combined with negative shifts in carbon isotopic compositions. The maxima were detected just above the hydrate layers, indicating that methane stored in the hydrates may be available for the microbial community. The observed variations in biomarker abundances and 13C-depletions are indicative of multiple environmental and physiological factors selecting for different AOM consortia (ANME-2a/DSS, ANME-2c/DSS, ANME-1) along horizontal and vertical gradients of cold seep settings.

The impact of ocean acidification and warming on the skeletal mechanical properties of the sea urchin Paracentrotus lividus from laboratory and field observations

Increased atmospheric CO2 concentration is leading to changes in the carbonate chemistry and the temperature of the ocean. The impact of these processes on marine organisms will depend on their ability to cope with those changes, particularly the maintenance of calcium carbonate structures. Both a laboratory experiment (long-term exposure to decreased pH and increased temperature) and collections of individuals from natural environments characterized by low pH levels (individuals from intertidal pools and around a CO2 seep) were here coupled to comprehensively study the impact of near-future conditions of pH and temperature on the mechanical properties of the skeleton of the euechinoid sea urchin Paracentrotus lividus. To assess skeletal mechanical properties, we characterized the fracture force, Young's modulus, second moment of area, material nanohardness, and specific Young's modulus of sea urchin test plates. None of these parameters were significantly affected by low pH and/or increased temperature in the laboratory experiment and by low pH only in the individuals chronically exposed to lowered pH from the CO2 seeps. In tidal pools, the fracture force was higher and the Young's modulus lower in ambital plates of individuals from the rock pool characterized by the largest pH variations but also a dominance of calcifying algae, which might explain some of the variation. Thus, decreases of pH to levels expected for 2100 did not directly alter the mechanical properties of the test of P. lividus. Since the maintenance of test integrity is a question of survival for sea urchins and since weakened tests would increase the sea urchins' risk of predation, our findings indicate that the decreasing seawater pH and increasing seawater temperature expected for the end of the century should not represent an immediate threat to sea urchins vulnerability