Particle camera and settling chamber measurements at site PC (GeoB11834-3) and SC (GeoB12914-4) off Cape Blanc/Mauritania

Particles sinking out of the euphotic zone are important vehicles of carbon export from the surface ocean. Most of the particles produce heavier aggregates by coagulating with each other before they sink. We implemented an aggregation model into the biogeochemical model of Regional Oceanic Modelling System (ROMS) to simulate the distribution of particles in the water column and their downward transport in the Northwest African upwelling region. Accompanying settling chamber, sediment trap and particle camera measurements provide data for model validation. In situ aggregate settling velocities measured by the settling chamber were around 55 m d**-1. Aggregate sizes recorded by the particle camera hardly exceeded 1 mm. The model is based on a continuous size spectrum of aggregates, characterised by the prognostic aggregate mass and aggregate number concentration. Phytoplankton and detritus make up the aggregation pool, which has an averaged, prognostic and size dependent sinking. Model experiments were performed with dense and porous approximations of aggregates with varying maximum aggregate size and stickiness as well as with the inclusion of a disaggregation term. Similar surface productivity in all experiments has been generated in order to find the best combination of parameters that produce measured deep water fluxes. Although the experiments failed to represent surface particle number spectra, in the deep water some of them gave very similar slope and spectrum range as the particle camera observations. Particle fluxes at the mesotrophic sediment trap site off Cape Blanc (CB) have been successfully reproduced by the porous experiment with disaggregation term when particle remineralisation rate was 0.2 d**-1. The aggregation-disaggregation model improves the prediction capability of the original biogeochemical model significantly by giving much better estimates of fluxes for both upper and lower trap. The results also point to the need for more studies to enhance our knowledge on particle decay and its variation and to the role that stickiness play in the distribution of vertical fluxes.

Magnesium isotope ratios of pore-fluids and dolomites

Magnesium concentrations in deep-sea sediment pore-fluids typically decrease down core due to net precipitation of dolomite or clay minerals in the sediments or underlying crust. To better characterize and differentiate these processes, we have measured magnesium isotopes in pore-fluids and sediment samples from Ocean Drilling Program sites (1082, 1086, 1012, 984, 1219, and 925) that span a range of oceanographic settings. At all sites, magnesium concentrations decrease with depth. At sites where diagenetic reactions are dominated by the respiration of organic carbon, pore-fluid d26Mg values increase with depth by as much as 2 per mil. Because carbonates preferentially incorporate 24Mg (low d26Mg), the increase in pore-fluid d26Mg values at these sites is consistent with the removal of magnesium in Mg-carbonate (dolomite). In contrast, at sites where the respiration of organic carbon is not important and/or weatherable minerals are abundant, pore-fluid d26Mg values decrease with depth by up to 2 per mil. The decline in pore-fluid d26Mg at these sites is consistent with a magnesium sink that is isotopically enriched relative to the pore-fluid. The identity of this enriched magnesium sink is likely clay minerals. Using a simple 1D diffusion-advection-reaction model of pore-fluid magnesium, we estimate rates of net magnesium uptake/removal and associated net magnesium isotope fractionation factors for sources and sinks at all sites. Independent estimates of magnesium isotope fractionation during dolomite precipitation from measured d26Mg values of dolomite samples from sites 1082 and 1012 are very similar to modeled net fractionation factors at these sites, suggesting that local exchange of magnesium between sediment and pore-fluid at these sites can be neglected. Our results indicate that the magnesium incorporated in dolomite is 2.0-2.7 per mil depleted in d26Mg relative to the precipitating fluid. Assuming local exchange of magnesium is minor at the rest of the studied sites, our results suggest that magnesium incorporated into clay minerals is enriched in d26Mg by 0 per mil to +1.25 per mil relative to the precipitating fluid. This work demonstrates the utility of magnesium isotopes as a tracer for magnesium sources/sinks in low-temperature aqueous systems.

Radium in Arctic surface water

The transpolar drift is strongly enriched in 228Ra accumulated on the wide Arctic shelves with subsequent rapid offshore transport. We present new data of Polarstern expeditions to the central Arctic and to the Kara and Laptev seas. Because 226Ra activities in Pacific waters are 30% higher than in Atlantic waters, we correct 226Ra for the Pacific admixture when normalizing 228Ra with 226Ra. The use of 228Ra decay as age marker critically depends on the constancy in space and time of the source activity, a condition that has not yet adequately been tested. While 228Ra decays during transit over the central basin, ingrowth of 228Th could provide an alternative age marker. The high 228Th/228Ra activity ratio (AR = 0.8-1.0) in the central basins is incompatible with a mixing model based on horizontal eddy diffusion. An advective model predicts that 228Th grows to an equilibrium AR, the value of which depends on the scavenging regime. The low AR over the Lomonosov Ridge (AR = 0.5) can be due to either rapid transport (minimum age without scavenging 1.1 year) or enhanced scavenging. Suspended particulate matter load (derived from beam transmission and particulate 234Th) and total 234Th depletion data show that scavenging, although extremely low in the central Arctic, is enhanced over the Lomonosov Ridge, making an age of 3 years more likely. The combined data of 228Ra decay and 228Th ingrowth confirm the existence of a recirculating gyre in the surface water of the eastern Eurasian Basin with a river water residence time of at least 3 years.

Model of interactive visualization of diffusion processes from different strains, and source routes and spatial diffusion of capuchin monkeys over the past 6 million years

Source routes and Spatial Diffusion of capuchin monkeys over the past 6 million years, rebuilt in the SPREAD 1.0.6 from the MCC tree. The map shows the 10 different regions to which distinctive samples were associated. The different transmission routes have been calculated from the average rate over time. Only rates with Bayes factor> 3 were considered as significantly different from zero. Significant diffusion pathways are highlighted with color varying from dark brown to red, being the dark brown less significant rates and deep red the most significant rates.

ODP Site 162-980 age model and color reflectance record of sediments from ODP Leg 162 sites

Sediment spectral reflectance measurements were generated aboard the JOIDES Resolution during Ocean Drilling Program Leg 162 shipboard operations. The large size of the raw data set (over 1.3 gigabytes) and limited computer hard disk storage space precluded detailed analysis of the data at sea, although broad band averages were used as aids in developing splices and determining lithologic boundaries. This data report describes the methods used to collect these data and their shipboard and postcruise processing. These initial results provide the basis for further postcruise research.

Salinization Model with links to source files

This computational model of irrigation agriculture is used to study the effects of salinization in Mesopotamia. Scholars have long suspected that central and southern Mesopotamia present environments which limited agricultural production over long-term periods. In regions such as central Mesopotamia, where salinization likely affected settlement in different periods but was more manageable than in more southern regions, fallowing regimes, natural and engineered leaching, and decisions made on when to crop were strategies applied in order to limit the effects of salinization. The model is used to assess the effectiveness of these coping strategies by incorporating projected climate, soil, and landscape conditions with agricultural practices. The simulation results not only demonstrate the effectiveness and limitations of techniques to inhibiting progressive salinization but can be compared with the archaeological record in order to determine if the results correspond to past events and help to interpret past settlement history.

Late glacial and deglacial age model and biomarker analyses of sediment core MSM05/5_712-2 from the western continental slope of Svalbard

The transition from last glacial to deglacial and subsequently to modern interglacial climate conditions was accompanied by abrupt shifts in the palaeoceanographic setting in the subpolar North Atlantic. Knowledge about the role that sea ice coverage played during these rapid climate reversals is limited since most marine sediment cores from the higher latitudes provide only a coarse temporal resolution and often poorly preserved microfossils. Here we present a highly resolved reconstruction of sea ice conditions that characterised the eastern Fram Strait - a key area for water mass exchange between the Arctic Ocean and the North Atlantic - for the past 30 ka BP. This reconstruction is based on the distribution of the sea ice biomarker IP25 and phytoplankton derived biomarkers in a sediment core from the continental slope of western Svalbard. During the late glacial (30 ka to 19 ka BP), recurrent advances and retreats of sea ice characterised the study area and point to a hitherto less considered oceanic (and/or atmospheric) variability. A long-lasting perennial sea ice coverage in eastern Fram Strait persisted only at the very end of the Last Glacial Maximum (i.e. from 19.2 to 17.6 ka BP) and was abruptly reduced at the onset of Heinrich Event 1 - coincident with or possibly even inducing the collapse of the Atlantic Meridional Overturning Circulation (AMOC). Further maximum sea ice conditions prevailed during the Younger Dryas cooling event and support the assumption of an AMOC reduction due to increased formation and export of Arctic sea ice through Fram Strait. A significant retreat of sea ice and sea surface warming are observed for the Early Holocene.

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.