Modeling the marine aragonite cycle: changes under rising carbon dioxide and its role in shallow water CaCO3 dissolution

The marine aragonite cycle has been included in the global biogeochemical model PISCES to study the role of aragonite in shallow water CaCO3 dissolution. Aragonite production is parameterized as a function of mesozooplankton biomass and aragonite saturation state of ambient waters. Observation-based estimates of marine carbonate production and dissolution are well reproduced by the model and about 60% of the combined CaCO3 water column dissolution from aragonite and calcite is simulated above 2000 m. In contrast, a calcite-only version yields a much smaller fraction. This suggests that the aragonite cycle should be included in models for a realistic representation of CaCO3 dissolution and alkalinity. For the SRES A2 CO2 scenario, production rates of aragonite are projected to notably decrease after 2050. By the end of this century, global aragonite production is reduced by 29% and total CaCO3 production by 19% relative to pre-industrial. Geographically, the effect from increasing atmospheric CO2, and the subsequent reduction in saturation state, is largest in the subpolar and polar areas where the modeled aragonite production is projected to decrease by 65% until 2100.

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

Mycobacterium ulcerans persistence at a village water source of Buruli ulcer patients

Buruli ulcer (BU), a neglected tropical disease of the skin and subcutaneous tissue, is caused by Mycobacterium ulcerans and is the third most common mycobacterial disease after tuberculosis and leprosy. While there is a strong association of the occurrence of the disease with stagnant or slow flowing water bodies, the exact mode of transmission of BU is not clear. M. ulcerans has emerged from the environmental fish pathogen M. marinum by acquisition of a virulence plasmid encoding the enzymes required for the production of the cytotoxic macrolide toxin mycolactone, which is a key factor in the pathogenesis of BU. Comparative genomic studies have further shown extensive pseudogene formation and downsizing of the M. ulcerans genome, indicative for an adaptation to a more stable ecological niche. This has raised the question whether this pathogen is still present in water-associated environmental reservoirs. Here we show persistence of M. ulcerans specific DNA sequences over a period of more than two years at a water contact location of BU patients in an endemic village of Cameroon. At defined positions in a shallow water hole used by the villagers for washing and bathing, detritus remained consistently positive for M. ulcerans DNA. The observed mean real-time PCR Ct difference of 1.45 between the insertion sequences IS2606 and IS2404 indicated that lineage 3 M. ulcerans, which cause human disease, persisted in this environment after successful treatment of all local patients. Underwater decaying organic matter may therefore represent a reservoir of M. ulcerans for direct infection of skin lesions or vector-associated transmission.

Toward explaining the Holocene carbon dioxide and carbon isotope records: Results from transient ocean carbon cycle-climate simulations

[1] The Bern3D model was applied to quantify the mechanisms of carbon cycle changes during the Holocene (last 11,000 years). We rely on scenarios from the literature to prescribe the evolution of shallow water carbonate deposition and of land carbon inventory changes over the glacial termination (18,000 to 11,000 years ago) and the Holocene and modify these scenarios within uncertainties. Model results are consistent with Holocene records of atmospheric CO2 and δ13C as well as the spatiotemporal evolution of δ13C and carbonate ion concentration in the deep sea. Deposition of shallow water carbonate, carbonate compensation of land uptake during the glacial termination, land carbon uptake and release during the Holocene, and the response of the ocean-sediment system to marine changes during the termination contribute roughly equally to the reconstructed late Holocene pCO2 rise of 20 ppmv. The 5 ppmv early Holocene pCO2 decrease reflects terrestrial uptake largely compensated by carbonate deposition and ocean sediment responses. Additional small contributions arise from Holocene changes in sea surface temperature, ocean circulation, and export productivity. The Holocene pCO2 variations result from the subtle balance of forcings and processes acting on different timescales and partly in opposite direction as well as from memory effects associated with changes occurring during the termination. Different interglacial periods with different forcing histories are thus expected to yield different pCO2 evolutions as documented by ice cores.

Sediment core fossils in ancient Lake Ohrid: testing for faunal change since the Last Interglacial

Ancient Lake Ohrid is probably of early Pleistocene or Pliocene origin and amongst the few lakes in the world harbouring an outstanding degree of endemic biodiversity. Although there is a long history of evolutionary research in Lake Ohrid, particularly on molluscs, a mollusc fossil record has been missing up to date. For the first time, gastropod and bivalve fossils are reported from the basal, calcareous part of a 2.6 m long sediment succession (core Co1200) from the north-eastern part of Lake Ohrid. Electron spin resonance (ESR) dating of mollusc shells from the same stratigraphic level yielded an age of 130 ± 28 ka. Lithofacies III sediments, i.e. a stratigraphic subdivision comprising the basal succession of core Co1200 between 181.5–263 cm, appeared solid, greyish-white, and consisted almost entirely of silt-sized endogenic calcite (CaCO3>70%) and intact and broken mollusc shells. Here we compare the faunal composition of the thanatocoenosis with recent mollusc associations in Lake Ohrid. A total of 13 mollusc species (9 gastropod and 4 bivalve species) could be identified within Lithofacies III sediments. The value of sediment core fossils for reconstructing palaeoenvironmental settings was evaluated and the agreement between sediment and palaeontological proxies was tested. The study also aims at investigating major faunal changes since the Last Interglacial and searching for signs of extinction events. The combined findings of the ecological study and the sediment characteristics suggest deposition in a shallow water environment during the Last Interglacial. The fossil fauna exclusively included species also found in the present fauna, i.e. no extinction events are evident for this site since the Last Interglacial. The thanatocoenosis showed the highest similarity with recent Intermediate Layer (5–25 m water depth) mollusc assemblages. The demonstrated existence of a mollusc fossil record in Lake Ohrid sediment cores also has great significance for future deep drilling projects. It can be hoped that a more far reaching mollusc fossil record will then be obtained, enabling insight into the early evolutionary history of Lake Ohrid.

High-resolution late-glacial chronology for the Gerzensee lake record (Switzerland): δ18O correlation between a Gerzensee-stack and NGRIP

Oxygen-isotope variations were analyzed on bulk samples of shallow-water lake marl from Gerzensee, Switzerland, in order to evaluate major and minor climatic oscillations during the late-glacial. To highlight the overall signature of the Gerzensee δ18O record, δ18O records of four parallel sediment cores were first correlated by synchronizing major isotope shifts and pollen abundances. Then the records were stacked with a weighting depending on the differing sampling resolution. To develop a precise chronology, the δ18O-stack was then correlated with the NGRIP δ18O record applying a Monte Carlo simulation, relying on the assumption that the shifts in δ18O were climate-driven and synchronous in both archives. The established chronology on the GICC05 time scale is the basis for (1) comparing the δ18O changes recorded in Gerzensee with observed climatic and environmental fluctuations over the whole North Atlantic region, and (2) comparing sedimentological and biological changes during the rapid warming with smaller climatic variations during the Bølling/Allerød period. The δ18O record of Gerzensee is characterized by two major isotope shifts at the onset and at the termination of the Bølling/Allerød warm period, as well as four intervening negative shifts labeled GI-1e2, d, c2, and b, which show a shift of one third to one fourth of the major δ18O shifts at the beginning and end of the Bølling/Allerød. Despite some inconsistency in terminology, these oscillations can be observed in various climatic proxies over wide regions in the North Atlantic region, especially in reconstructed colder temperatures, and they seem to be caused by hemispheric climatic variations.

Spectroscopy and detectability of liquid brines on mars

Recent geomorphological observations as well as chemical and thermodynamic studies demonstrate that liquid water should be stable today on the Martian surface at some times of the day. In Martian conditions, brines would be particularly more stable than pure water because salts can depress the freezing point and lower the evaporation rate of water. Despite this evidence, no clear spectral signature of liquid has been observed so far by the hyperspectral imaging spectrometers OMEGA and CRISM. However, past spectral analysis lacks a good characterization of brines׳ spectral signatures. This study thus aims to determine how liquid brines can be detected on Mars by spectroscopy. In this way, laboratory experiments were performed for reproducing hydration and dehydration cycles of various brines while measuring their spectral signatures. The resulting spectra first reveal a very similar spectral evolution for the various brine types and pure water, with the main difference observed at the end of the dehydration with the crystallization of various hydrated minerals from brines. The main characteristic of this spectral behavior is an important decoupling between the evolution of albedo and hydration bands depths. During most of the wetting/drying processes, spectra usually display a low albedo associated with shallow water absorption band depths. Strong water absorption band depth and high albedo are respectively only observed when the surface is very wet and when the surface is very dry. These experiments can thus explain why the currently active Martian features attributed to the action of a liquid are only associated with low albedo and very weak spectral signatures. Hydration experiments also reveal that deliquescence occurs easily even at low temperature and moderate soil water vapor pressure and could thus cause seasonal darkening on Mars. These experiments demonstrate that the absence of water absorptions in CRISM in the middle afternoon does not rule out water activity and suggest future spectral investigations to identify water on the Martian surface.

Constraints on barium isotope fractionation during aragonite precipitation by corals

We present a barium (Ba) isotope fractionation study of marine biogenic carbonates (aragonitic corals). The major aim is to provide first constraints on the Ba isotope fractionation between modern surface sea water and coral skele- ton. Mediterranean surface sea water was found to be enriched in the heavy Ba isotopes compared to previously reported values for marine open ocean authi- genic and terrestrial minerals. In aquarium experiments with a continuous sup- ply of Mediterranean surface water, the Ba isotopic composition of the bulk sample originating from cultured, aragonitic scleractinian corals (d137/134Ba between +0.16 +/- 0.12permil and +0.41 +/-0.12permil) were isotopically identical or lighter than that of the ambient Mediterranean surface sea water (d137/134Ba = +0.42 +/- 0.07permil, 2SD), which corresponds to an empirical maximum value of Ba isotope fractionation of D137/134Bacoral-seawater = -0.26 +/- 0.14permil at 25°C. This maximum Ba isotope fractionation is close and identical in direction to previous results from inorganic precipitation experiments with aragonite- structured pure BaCO3 (witherite). The variability in measured Ba concentrations of the cultured corals is at odds with a uniform distribution coefficient, DBa/Ca, thus indicating stronger vital effects on isotope than element discrimination. This observation supports the hypothesis that the Ba isotopic compositions of these corals do not result from simple equilibrium between the skeleton and the bulk sea water. Complementary coral samples from natural settings (tropical shallow-water corals from the Bahamas and Florida and cold- water corals from the Norwegian continental shelf) show an even wider range in d137/134Ba values (+0.14 +/- 0.08permil and +0.77 +/- 0.11permil), most probably due to additional spatial and/or temporal sea water heterogeneity, as indicated by recent publications.

Comparative carbon cycle dynamics of the present and last interglacial

Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO₂ remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO₂ fertilization, land use, wild fire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO₂ dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO₂ and δ¹³ CO₂ changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO₂ dynamics from 8ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO₂ changes after 122 ka BP. This failure to simulate late-Eemian CO₂ dynamics could be a result of the imposed forcings such as prescribed CaCO₃ accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO₂ dynamics eshallow water CaCO₃ accumulation, peat and permafrost carbon dynamics are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO₂ dynamics.

Importance of vegetation, topography and flow paths for water transit times of base flow in alpine headwater catchments

The mean transit time (MTT) of water in a catchment gives information about storage, flow paths, sources of water and thus also about retention and release of solutes in a catchment. To our knowledge there are only a few catchment studies on the influence of vegetation cover changes on base flow MTTs. The main changes in vegetation cover in the Swiss Alps are massive shrub encroachment and forest expansion into formerly open habitats. Four small and relatively steep headwater catchments in the Swiss Alps (Ursern Valley) were investigated to relate different vegetation cover to water transit times. Time series of water stable isotopes were used to calculate MTTs. The high temporal variation of the stable isotope signals in precipitation was strongly dampened in stream base flow samples. MTTs of the four catchments were 70 to 102 weeks. The strong dampening of the stable isotope input signal as well as stream water geochemistry points to deeper flow paths and mixing of waters of different ages at the catchments' outlets. MTTs were neither related to topographic indices nor vegetation cover. The major part of the quickly infiltrating precipitation likely percolates through fractured and partially karstified deeper rock zones, which increases the control of bedrock flow paths on MTT. Snow accumulation and the timing of its melt play an important role for stable isotope dynamics during spring and early summer. We conclude that, in mountainous headwater catchments with relatively shallow soil layers, the hydrogeological and geochemical patterns (i.e. geochemistry, porosity and hydraulic conductivity of rocks) and snow dynamics influence storage, mixing and release of water in a stronger way than vegetation cover or topography do.