Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

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.

(Table S1) Tie-points used in age model construction for ODP Leg 172 sites

ix Ocean Drilling Program (ODP) sites, in the Northwest Atlantic have been used to investigate kinematic and chemical changes in the "Western Boundary Undercurrent" (WBUC) during the development of full glacial conditions across the Marine Isotope Stage 5a/4 boundary (~70,000 years ago). Sortable silt mean grain size(sort s) measurements are employed to examine changes in near bottom flow speeds, together with carbon isotopes measured in benthic foraminifera and % planktic foraminiferal fragmentation as proxies for changes in water-mass chemistry. A depth transect of cores, spanning 1.8-4.6 km depth, allows changes in both the strength and depth of the WBUC to be constrained across millennial scale events. Sort s measurements reveal that the flow speed structure of the WBUC during warm intervals ("interstadials") was comparable to modern (Holocene) conditions. However, significant differences are observed during cold intervals, with higher relative flow speeds inferred for the shallow component of the WBUC (~2 km depth) during all cold "stadial" intervals (including Heinrich Stadial 6), and a substantial weakening of the deep component (~3-4 km) during full glacial conditions. Our results therefore reveal that the onset of full glacial conditions was associated with a regime shift to a shallower mode of circulation (involving Glacial North Atlantic Intermediate Water) that was quantitatively distinct from preceding cold stadial events. Furthermore, our chemical proxy data show that the physical response of the WBUC during the last glacial inception was probably coupled to basin-wide changes in the water-mass composition of the deep Northwest Atlantic.