Polychlorinated Biphenyls in Glaciers. 2. Model Results of Deposition and Incorporation Processes

In previous work, Alpine glaciers have been identified as a secondary source of persistent organic pollutants (POPs). However, detailed understanding of the processes organic chemicals undergo in a glacial system was missing. Here, we present results from a chemical fate model describing deposition and incorporation of polychlorinated biphenyls (PCBs) into an Alpine glacier (Fiescherhorn, Switzerland) and an Arctic glacier (Lomonosovfonna, Norway). To understand PCB fate and dynamics, we investigate the interaction of deposition, sorption to ice and particles in the atmosphere and within the glacier, revolatilization, diffusion and degradation, and discuss the effects of these processes on the fate of individual PCB congeners. The model is able to reproduce measured absolute concentrations in the two glaciers for most PCB congeners. While the model generally predicts concentration profiles peaking in the 1970s, in the measurements, this behavior can only be seen for higher-chlorinated PCB congeners on Fiescherhorn glacier. We suspect seasonal melt processes are disturbing the concentration profiles of the lower-chlorinated PCB congeners. While a lower-chlorinated PCB congener is mainly deposited by dry deposition and almost completely revolatilized after deposition, a higher-chlorinated PCB congener is predominantly transferred to the glacier surface by wet deposition and then is incorporated into the glacier ice. The incorporated amounts of PCBs are higher on the Alpine glacier than on the Arctic glacier due to the higher precipitation rate and aerosol particle concentration on the former. Future studies should include the effects of seasonal melt processes, calculate the quantities of PCBs incorporated into the entire glacier surface, and estimate the quantity of chemicals released from glaciers to determine the importance of glaciers as a secondary source of organic chemicals to remote aquatic ecosystems.

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.

Wet deposition of fossil and non-fossil derived particulate carbon: insights from radiocarbon measurement

Radiocarbon (14C) measurements of both organic carbon (OC) and elemental carbon (EC) allow a more detailed source apportionment, leading to a full and unambiguous distinction and quantification of the contributions from non-fossil and fossil sources. A thermal-optical method with a commercial OC/EC analyzer to isolate water-insoluble OC (WIOC) and EC for their subsequent 14C measurement was applied for the first time to filtered precipitation samples collected at a costal site in Portugal and at a continental site in Switzerland. Our results show that WIOC in precipitation is dominated by non-fossil sources such as biogenic and biomass-burning emissions regardless of rain origins and seasons, whereas EC sources are shared by fossil-fuel combustion and biomass burning. In addition, monthly variation of WIOC in Switzerland was characterized by higher abundance in warm than in cold seasons, highlighting the importance of biogenic emissions to particulate carbon in rainwater. Samples with high particulate carbon concentrations in Portugal were found to be associated with increased biogenic input. Despite the importance of non-fossil sources, fossil emissions account for approximately 20% of particulate carbon in wet deposition for our study, which is in line with fossil contribution in bulk rainwater dissolved organic carbon as well as aerosol WIOC and EC estimated by the 14C approach from other studies.

Flood triggering in Switzerland: the role of daily to monthly preceding precipitation

Determining the role of different precipitation periods for peak discharge generation is crucial for both projecting future changes in flood probability and for short- and medium-range flood forecasting. In this study, catchment-averaged daily precipitation time series are analyzed prior to annual peak discharge events (floods) in Switzerland. The high number of floods considered – more than 4000 events from 101 catchments have been analyzed – allows to derive significant information about the role of antecedent precipitation for peak discharge generation. Based on the analysis of precipitation times series, a new separation of flood-related precipitation periods is proposed: (i) the period 0 to 1 day before flood days, when the maximum flood-triggering precipitation rates are generally observed, (ii) the period 2 to 3 days before flood days, when longer-lasting synoptic situations generate "significantly higher than normal" precipitation amounts, and (iii) the period from 4 days to 1 month before flood days when previous wet episodes may have already preconditioned the catchment. The novelty of this study lies in the separation of antecedent precipitation into the precursor antecedent precipitation (4 days before floods or earlier, called PRE-AP) and the short range precipitation (0 to 3 days before floods, a period when precipitation is often driven by one persistent weather situation like e.g., a stationary low-pressure system). A precise separation of "antecedent" and "peak-triggering" precipitation is not attempted. Instead, the strict definition of antecedent precipitation periods permits a direct comparison of all catchments. The precipitation accumulating 0 to 3 days before an event is the most relevant for floods in Switzerland. PRE-AP precipitation has only a weak and region-specific influence on flood probability. Floods were significantly more frequent after wet PRE-AP periods only in the Jura Mountains, in the western and eastern Swiss plateau, and at the outlet of large lakes. As a general rule, wet PRE-AP periods enhance the flood probability in catchments with gentle topography, high infiltration rates, and large storage capacity (karstic cavities, deep soils, large reservoirs). In contrast, floods were significantly less frequent after wet PRE-AP periods in glacial catchments because of reduced melt. For the majority of catchments however, no significant correlation between precipitation amounts and flood occurrences is found when the last 3 days before floods are omitted in the precipitation amounts. Moreover, the PRE-AP was not higher for extreme floods than for annual floods with a high frequency and was very close to climatology for all floods. The fact that floods are not significantly more frequent nor more intense after wet PRE-AP is a clear indicator of a short discharge memory of Pre-Alpine, Alpine and South Alpine Swiss catchments. Our study poses the question whether the impact of long-term precursory precipitation for floods in such catchments is not overestimated in the general perception. The results suggest that the consideration of a 3–4 days precipitation period should be sufficient to represent (understand, reconstruct, model, project) Swiss Alpine floods.

Wet and dry deposition of 129I in Seville (Spain) measured by accelerator mass spectrometry

Iodine-129 (Full-size image (<1 K)) concentrations have been determined by accelerator mass spectrometry in rainwater samples taken at Seville (southwestern Spain) in 1996 and 1997. This technique allows a reduction in the detection limits for this radionuclide in comparison to radiometric counting and other mass spectrometric methods such as ICP-MS. Typical 129I concentrations range from 4.7×107129I atoms/l (19.2%) to 4.97×109129I atoms/l (5.9%), while 129I depositions are normally in the order of 108–1010 atoms/m2 d. These values agree well with other results obtained for recent rainwater samples collected in Europe. Apart from these, the relationship between 129I deposition and some atmospheric factors has been analyzed, showing the importance of the precipitation rate and the concentration of suspended matter in it.

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.