Glacier mass balance reconstruction by sublimation induced enrichment of chemical species on Cerro Tapado (Chilean Andes)

A 36 m long ice core down to bedrock from the Cerro Tapado glacier (5536 m a.s.l, 30°08' S, 69°55' W) was analyzed to reconstruct past climatic conditions for Northern Chile. Because of the marked seasonality in the precipitation (short wet winter and extended dry summer periods) in this region, major snow ablation and related post-depositional processes occur on the glacier surface during summer periods. They include predominantly sublimation and dry deposition. Assuming that, like measured during the field campaign, the enrichment of chloride was always related to sublimation, the chemical record along the ice core may be applied to reconstruct the history of such secondary processes linked to the past climatic conditions over northern Chile. For the time period 1962–1999, a mean annual net accumulation of 316 mm water equivalent (weq) and 327 mm weq loss by sublimation was deduced by this method. This corresponds to an initial total annual accumulation of 539 mm weq. The annual variability of the accumulation and sublimation is related with the Southern Oscillation Index (SOI): higher net-accumulation during El-Niño years and more sublimation during La Niña years. The deepest part of the ice record shows a time discontinuity; with an ice body deposited under different climatic conditions: 290 mm higher precipitation but with reduced seasonal distribution (+470 mm in winter and –180 mm in summer) and –3°C lower mean annual temperature. Unfortunately, its age is unknown. The comparison with regional proxy data however let us conclude that the glacier buildup did most likely occur after the dry mid-Holocene.

The role of climate-driven chemical weathering on soil production

Climate plays an important role in controlling rates of weathering and weathered regolith production. Regolith production functions, however, seldom take climate parameters into account. Based on a climate-dependent weathered regolith production model, at low denudation rates, relative regolith thicknesses are less sensitive to changes in precipitation rates, while at high denudation rates, small changes in climatic parameters can result in complete stripping of hillslopes. This pattern is compounded by the long residence times and system response times associated with low denudation rates, and vice versa. As others have shown, the transition between regolith-mantled and bedrock slopes is dependent on the ratio of denudation to production. Here, we further suggest that this is itself a function of precipitation rate and temperature. We suggest that climatic parameters can be easily incorporated into existing soil production models and that such additions improve the predictive power of soil production models. (C) 2013 Elsevier B.V. All rights reserved.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

Temperature and precipitation effects on δ13C depth profiles in SOM under temperate beech forests

Enrichment of 13C in SOM with soil depth is related to interacting processes influenced by temperature and precipitation. Our objectives were to derive climate effects on patterns of vertical δ13C values of soil organic matter (SOM) while minimizing the effect of confounding variables. We investigated vertical changes in δ13C values of SOM in 1-cm depth intervals in silvicultural mature beech (Fagus sylvatica L.) forest ecosystems in northern Rhineland-Palatinate across gradients of MAT (7.9 to 9.7 °C mean annual temperature) and MAP (607 to 1085 mm mean annual precipitation) in winter 2011. Forest stands (n = 10) were chosen based on data sets provided by the Rhineland-Palatinate Forest Administration so that variations in these gradients occurred while other environmental factors like physico-chemical soil properties, tree species, stand age, exposition and precipitation (for the temperature gradient) or temperature (for the precipitation gradient) did not differ among study sites. From litter down to the mineral soil at 10 cm depth, soil organic carbon (SOC) content decreased (47.5 ± SE 0.1% to 2.5 ± 0.1%) while the δ13C values increased (− 29.4 ± 0.1‰ to − 26.1 ± 0.1‰). Litter of sites under higher MAP/lower MAT had lower δ13C values which was in line with literature data on climate driven plant physiological process. To compare the dimension of the vertical 13C enrichment, δ13C values were regressed linearly against log-transformed carbon contents yielding absolute values of these slopes (beta). Beta values ranged between 0.6 and 4.5 (range of r from − 0.7 to − 1.0; p < 0.01). Due to an assumed decay continuum and similar variations of δ13C values in litter and in 10 cm depth, we conclude that effects on isotope composition in the Oi layer continue vertically and therefore, δ13C values in litter do not solely control beta values. Beta values decreased with increasing MAT (r = − 0.83; p < 0.05). Reduced soil moisture and therefore both, reduced microbial activity and reduced downward transport of microbial cycled DOM (=13C enriched) might be responsible for less pronounced δ13C depth profiles in case of high temperatures. Greater C:N ratios (lower degradability) of the litter under higher temperatures likely contributed to these depth trends. Beta values increased with increasing MAP (r = 0.73; p < 0.05). We found decreasing C:N ratios in the mineral soil that possibly indicates higher decomposition under higher precipitation. Exclusion of the organic layers from linear regressions indicated a stronger impact of MAP on the development of δ13C depth profiles. Our results confirm temperature and precipitation effects on δ13C depth profiles and indicate stronger 13C enrichment under lower MAT/higher MAP. Therefore, time series of vertical δ13C depth profiles might provide insights into climate change effects.

Dissolution-Precipitation in Porous Media: Experiments and Modelling

The evolution of porosity due to dissolution/precipitation processes of minerals and the associated change of transport parameters are of major interest for natural geological environments and engineered underground structures. We designed a reproducible and fast to conduct 2D experiment, which is flexible enough to investigate several process couplings implemented in the numerical code OpenGeosys-GEM (OGS-GEM). We investigated advective-diffusive transport of solutes, effect of liquid phase density on advective transport, and kinetically controlled dissolution/precipitation reactions causing porosity changes. In addition, the system allowed to investigate the influence of microscopic (pore scale) processes on macroscopic (continuum scale) transport. A Plexiglas tank of dimension 10 × 10 cm was filled with a 1 cm thick reactive layer consisting of a bimodal grain size distribution of celestite (SrSO4) crystals, sandwiched between two layers of sand. A barium chloride solution was injected into the tank causing an asymmetric flow field to develop. As the barium chloride reached the celestite region, dissolution of celestite was initiated and barite precipitated. Due to the higher molar volume of barite, its precipitation caused a porosity decrease and thus also a decrease in the permeability of the porous medium. The change of flow in space and time was observed via injection of conservative tracers and analysis of effluents. In addition, an extensive post-mortem analysis of the reacted medium was conducted. We could successfully model the flow (with and without fluid density effects) and the transport of conservative tracers with a (continuum scale) reactive transport model. The prediction of the reactive experiments initially failed. Only the inclusion of information from post-mortem analysis gave a satisfactory match for the case where the flow field changed due to dissolution/precipitation reactions. We concentrated on the refinement of post-mortem analysis and the investigation of the dissolution/precipitation mechanisms at the pore scale. Our analytical techniques combined scanning electron microscopy (SEM) and synchrotron X-ray micro-diffraction/micro-fluorescence performed at the XAS beamline (Swiss Light Source). The newly formed phases include an epitaxial growth of barite micro-crystals on large celestite crystals (epitaxial growth) and a nano-crystalline barite phase (resulting from the dissolution of small celestite crystals) with residues of celestite crystals in the pore interstices. Classical nucleation theory, using well-established and estimated parameters describing barite precipitation, was applied to explain the mineralogical changes occurring in our system. Our pore scale investigation showed limits of the continuum scale reactive transport model. Although kinetic effects were implemented by fixing two distinct rates for the dissolution of large and small celestite crystals, instantaneous precipitation of barite was assumed as soon as oversaturation occurred. Precipitation kinetics, passivation of large celestite crystals and metastability of supersaturated solutions, i.e. the conditions under which nucleation cannot occur despite high supersaturation, were neglected. These results will be used to develop a pore scale model that describes precipitation and dissolution of crystals at the pore scale for various transport and chemical conditions. Pore scale modelling can be used to parameterize constitutive equations to introduce pore-scale corrections into macroscopic (continuum) reactive transport models. Microscopic understanding of the system is fundamental for modelling from the pore to the continuum scale.

Barium isotope fractionation during witherite (BaCO3) dissolution, precipitation and at equilibrium

This study examines the behavior of Ba isotope fractionation between witherite and fluid during mineral dissolution, precipitation and at chemical equilibrium. Experiments were performed in batch reactors at 25 oC in 10-2 M NaCl solution where the pH was adjusted by continuous bubbling of a water saturated gas phase of CO2 or atmospheric air. During witherite dissolution no Ba isotope fractionation was observed between solid and fluid. In contrast, during witherite precipitation, caused by a pH increase, a preferential uptake of the lighter 134Ba isotopomer in the solid phase was observed. In this case, the isotope fractionation factor αwitherite-fluid is calculated to be 0.99993 ± 0.00004 (or Δ137/134Bawitherite-fluid ≈ -0.07 ± 0.04 ‰, 2sd). The most interesting feature of this study, however, is that after the attainment of chemical equilibrium, the Ba isotope composition of the aqueous phase is progressively becoming lighter, indicating a continuous exchange of Ba2+ ions between witherite and fluid. Mass balance calculations indicate that the detachment of Ba from the solid is not only restricted to the outer surface layer of the solid, but affects several (~7 unit cells) subsurface layers of the crystal. This observation comes in excellent agreement with the concept of a dynamic system at chemical equilibrium in a mineral-fluid system, denoting that the time required for the achievement of isotopic equilibrium in the witherite-fluid system is longer compared to that observed for chemical equilibrium. Overall, these results indicate that the isotopic composition of Ba bearing carbonates in natural environments may be altered due to changes in fluid composition without a net dissolution/precipitation to be observed.