Environmental tracer transport (H and SF6) in the saturated and unsaturated zones and its use in nitrate pollution management

An important quantity in groundwater protection is the residence time of water in an aquifer. It relates to both the travel time of a pollutant to arrive at a well and the time span required for self-purification of a polluted aquifer after removal of pollutant inputs. Time scales for aquifers can be gained from artificial tracer experiments or from environmental tracer data, the latter offering the only realistic alternative if time scales of years or decades have to be taken into account.

Different tracers show different time scales due to their different transport mechanisms especially in the unsaturated zone. While solute tracers are moved advectively with the seepage water, gas tracers pass the unsaturated zone diffusively through the air phase. Depending on the properties of the unsaturated zone (hydraulic properties, thickness) this difference in behavior can be used to separate the subsurface transport process into the unsaturated and the saturated parts.

In a field study in Germany, SF₆ and ³H were used as environmental tracers. Both have a relatively well-known input function. Interpretation of data from observation wells by a box model approach led to spatially and temporally varying residence times. This was an indication that the influence of the unsaturated zone could not be neglected. While the gas tracer SF₆ shows only residence times in the saturated zone, the tracer ³H reflects the whole travel time of water including both the unsaturated and saturated zones. Using a one-dimensional plug-flow model for the unsaturated zone combined with a detailed two-dimensional flow and transport model for the saturated zone leads to a holistic and consistent interpretation of the measured tracer concentrations. The observed pattern of old water under thick loess cover and younger water under areas where the fractured basalt aquifer crops out is reproduced after adjusting only two parameters: the effective porosity of the saturated aquifer and the product of field capacity and thickness of the unsaturated zone. While the effective porosity of the saturated zone is adjusted by means of the SF₆ data, the field capacity of the loess layer is adjusted by means of the ³H observations. The thickness of the unsaturated zone is deduced from geological and pedological maps. All flow data are obtained from a calibrated flow model, which is based on geological data, observed heads and pumping tests only.

The transport model for the saturated zone was calibrated by fitting the porosity by means of gaseous tracer concentrations (SF₆). The combined saturated–unsaturated zone model was then calibrated by fitting the field capacity of the unsaturated zone by means of ³H concentrations. With this model it was possible to verify the observed NO₃ concentrations at the drinking water wells and to develop predictions for their future development under various scenarios of fertilizer input reduction in specific areas.     

Comparison of bryozoan assemblages from two contrasting Arctic shelf regions

The structure of bryozoan assemblages from two Arctic regions (East Greenland and West Spitsbergen shelf) was compared. Both areas are located at the same latitude and the samples were taken from similar depths; however, the regions differed in water temperature and oceanography. East Greenland, which has lower mean annual water temperatures, was found to be one third richer in taxa (86 species) than West Spitsbergen (59 species). Diversity (Shannon–Wiener index – H′) and abundance were also higher on average in East Greenland (e.g. H′ = 2.49) than West Spitsbergen (e.g. H′ = 2.12). However for species richness, diversity and abundance there were no significant statistical differences between means (ANOVA) from the two regions. In spite of these similarities the investigated assemblages differed to a large extent in both species composition and dominance structure. There were 55 species that occurred only in East Greenland and 28 species that were present only in West Spitsbergen. There was higher proportion of species with an Arctic distribution in East Greenland (46%) than in West Spitsbergen (24%). Observed dissimilarities were concluded to be due to different hydrological conditions between the two regions.     

Instrument intercomparison study on cloud droplet size distribution measurements: Holography vs. laser optical particle counter

During the EUROTRAC Ground Based Cloud Experiment (GCE) 1990, a newly developed HODAR (Holographic Droplet and Acrosol Recording) was operated for the first time to measure cloud droplet size distributions by recording Fraunhofer in-line holograms of small cloud sample volumes in situ and analyzing the holographic images in the laboratory.This technical note compares the resulting size distributions with those obtained from two FSSP-100 laser optical particle counters. For all holograms analyzed during the GCE90 field experiment, the size distributions obtained from the two different methods agree well. Additionally, the liquid water contents (LWC) were measured directly by a Gerber particulate volume monitor PVM-100. The LWC calculated from the measured droplet size distributions deviate from the PVM-100 data.     

Geochemical and petrological constraints on rear-arc magma genesis processes in Ecuador: The Puyo cones and Mera lavas volcanic formations

The Puyo scoria cones and the Mera lava flows, two newly recognized volcanic formations dated between Late Pliocene to Middle Pleistocene, extend the limits of the Ecuadorian rear-arc volcanic province some 100 km to the south. The Puyo scoria cones have erupted K-rich absarokites containing olivine, diopside and phlogopite, whereas the Mera lava flows display a basic andesite composition, with olivine and minor augite phenocrysts. In addition to high contents in LILE, LREE and HFSE, the Puyo absarokites exhibit many characteristics of primitive melts, namely high Cr (590–310 ppm) and Ni (330–154 ppm) contents, high Mg# (64–70) and they contain forsteritic olivine (Fo_82–89). The composition of the most primary Puyo absarokite was used in petrogenetic models, in order to constrain the genesis of these high-K magmas. Major and trace elements models, as well as isotopic data, indicate that the source of Puyo magmas is a hydrated phlogopite- and garnet-bearing lherzolite. Phlogopite crystallization in the mantle wedge is triggered by the metasomatism by 3–5% of a SiO₂-, H₂O-rich liquid generated by slab melting. Partial melting of the subducted oceanic crust beneath Ecuador is allowed by the subduction of the young and warm Carnegie Ridge, which modifies the thermal regime of the Benioff zone. A low degree (1–4%) of partial melting of the metasomatized mantle wedge, leaving a variable garnet (4–7%) ± phlogopite (0–4%) lherzolitic residual assemblage, leads to the compositions of the entire Puyo absarokite series and is consistent with previous petrogenetic models developed for the Ecuadorian volcanic arc. Indeed, the homogeneity of isotopic data across the arc suggests a similar source for the whole Ecuadorian magmas.     

Hydroxyl lines towards Puppis A

The 1667-MHz ground-state hydroxyl (OH) line has been observed towards positions on Puppis A as well as in the immediate vicinity of this supernova remnant (SNR). The spectra at two positions on Puppis A show absorption lines at velocities less than +7.6 km s¹, and emission lines above this velocity. Analysis of all four 18 cm OH lines indicates that Puppis A must lie between the OH clouds seen in absorption and emission, giving it a v _lsr of +7.6 km s¹. This implies that the kinematic distance to Puppis A is The OH clouds in front of Puppis A are slightly anomalously excited, as shown by analysis of the four 18 cm OH lines at the two observed positions on Puppis A. There is no evidence of any interaction of the SNR with the clouds in front of and behind it.     

Constraining the prograde and retrograde P-T paths of granulites using decomposition of initially zoned garnets: an example from the Central Indian Tectonic Zone

The present study from the Sausar Mobile Belt (SMB) in the southern part of the Central Indian Tectonic zone (CITZ) demonstrates how microdomainal compositional variation of a single garnet porphyroblast in a metapelite granulite sample records the different segments of a near complete P-T path of metamorphic evolution. The microdomainal variation is ascribed to the preservation of growth zoning and heterogeneous distribution of diverse inclusion mineral assemblages. Subsequent mineral reactions under changing P/T conditions were controlled by this compositional heterogeneity. Four stages of metamorphic evolution have been deciphered. An early prograde stage (M_o) is implied by the rare presence of staurolite-biotite-quartz and in places of kyanite inclusion assemblages in other metapelite samples, together with the growth zoning preserved in garnet. The peak metamorphism (M₁) at ~9.5 kbar, ~850 °C is consistent with the biotite dehydration melting that produced garnet-K-feldspar and granitic leucosomes. This was followed by near isothermal decompression (M₂) at ~6 kbar, ~825 °C, during which different garnet segments behaved as separate microscale bulk compositions and decomposed both internally and externally to produce different retrograde mineral assemblages. In the quartz-bearing domain of almandine-rich and grossular-rich garnet core, grossular components in garnet reacted with included sillimanite and quartz to produce coronal plagioclase (X_An=0.90). By contrast, grossular-rich garnet in quartz-absent domain reacted with included sillimanite to produce layered spinel_ss {X_Mg (Mg/Mg+Fe²⁺) = 0.23–0.26}, X_Al (Al/Al+Fe³⁺)=0.71–0.81}-plagioclase (X_An=0.91)-cordierite {X_Mg (Mg/Mg+Fe²⁺) = 0.80–0.83} coronas both in the core and inner rim region of garnet. During post-decompression cooling, reactions occurred at about 600 °C (M₃), whereby quartz-bearing, sillimanite-absent microdomains of pyrope-rich, grossular-poor garnet outer rim decomposed to form relatively magnesian assemblages of cordierite-anthophyllite and cordierite-biotite-quartz. M₂ spinel_ss decomposed to polyphase domains of spinel-magnetite±högbomite at this stage. Collating the textural and geothermobarometric results, a clockwise P-T path has been deduced. The deduced P-T loop is consistent with a model of crustal thickening due to continental collision, followed by rapid vertical thinning, which appears to be the general feature of the Sausar Mobile Belt. Using model calculations of the preserved growth and diffusion zoning in garnet, we demonstrate rather short-lived nature of this collision orogeny (in the order of 40–60 Ma).Editorial responsibility: W. Schreyer     

The vertical distribution of climate forcings and feedbacks from the surface to top of atmosphere

The radiative forcings and feedbacks that determine Earth’s climate sensitivity are typically defined at the top-of-atmosphere (TOA) or tropopause, yet climate sensitivity itself refers to a change in temperature at the surface. In this paper, we describe how TOA radiative perturbations translate into surface temperature changes. It is shown using first principles that radiation changes at the TOA can be equated with the change in energy stored by the oceans and land surface. This ocean and land heat uptake in turn involves an adjustment of the surface radiative and non-radiative energy fluxes, with the latter being comprised of the turbulent exchange of latent and sensible heat between the surface and atmosphere. We employ the radiative kernel technique to decompose TOA radiative feedbacks in the IPCC Fourth Assessment Report climate models into components associated with changes in radiative heating of the atmosphere and of the surface. (We consider the equilibrium response of atmosphere-mixed layer ocean models subjected to an instantaneous doubling of atmospheric CO₂). It is shown that most feedbacks, i.e., the temperature, water vapor and cloud feedbacks, (as well as CO₂ forcing) affect primarily the turbulent energy exchange at the surface rather than the radiative energy exchange. Specifically, the temperature feedback increases the surface turbulent (radiative) energy loss by 2.87 W m^?2 K^?1 (0.60 W m^?2 K^?1) in the multimodel mean; the water vapor feedback decreases the surface turbulent energy loss by 1.07 W m^?2 K^?1 and increases the surface radiative heating by 0.89 W m^?2 K^?1; and the cloud feedback decreases both the turbulent energy loss and the radiative heating at the surface by 0.43 and 0.24 W m^?2 K^?1, respectively. Since changes to the surface turbulent energy exchange are dominated in the global mean sense by changes in surface evaporation, these results serve to highlight the fundamental importance of the global water cycle to Earth’s climate sensitivity.     

Adakite-like volcanism of Ecuador: lower crust magmatic evolution and recycling

In the Northern Andes of Ecuador, a broad Quaternary volcanic arc with significant across-arc geochemical changes sits upon continental crust consisting of accreted oceanic and continental terranes. Quaternary volcanic centers occur, from west to east, along the Western Cordillera (frontal arc), in the Inter-Andean Depression and along the Eastern Cordillera (main arc), and in the Sub-Andean Zone (back-arc). The adakite-like signatures of the frontal and main arc volcanoes have been interpreted either as the result of slab melting plus subsequent slab melt–mantle interactions or of lower crustal melting, fractional crystallization, and assimilation processes. In this paper, we present petrographic, geochemical, and isotopic (Sr, Nd, Pb) data on dominantly andesitic to dacitic volcanic rocks as well as crustal xenolith and cumulate samples from five volcanic centers (Pululagua, Pichincha, Ilalo, Chacana, Sumaco) forming a NW–SE transect at about 0° latitude and encompassing the frontal (Pululagua, Pichincha), main (Ilalo, Chacana), and back-arc (Sumaco) chains. All rocks display typical subduction-related geochemical signatures, such as Nb and Ta negative anomalies and LILE enrichment. They show a relative depletion of fluid-mobile elements and a general increase in incompatible elements from the front to the back-arc suggesting derivation from progressively lower degrees of partial melting of the mantle wedge induced by decreasing amounts of fluids released from the slab. We observe widespread petrographic evidence of interaction of primary melts with mafic xenoliths as well as with clinopyroxene- and/or amphibole-bearing cumulates and of magma mixing at all frontal and main arc volcanic centers. Within each volcanic center, rocks display correlations between evolution indices and radiogenic isotopes, although absolute variations of radiogenic isotopes are small and their values are overall rather primitive (e.g., ε_Nd = +1.5 to +6, ⁸⁷Sr/⁸⁶Sr = 0.7040–0.70435). Rare earth element patterns are characterized by variably fractionated light to heavy REE (La/Yb_N = 5.7–34) and by the absence of Eu negative anomalies suggesting evolution of these rocks with limited plagioclase fractionation. We interpret the petrographic, geochemical, and isotopic data as indicating open-system evolution at all volcanic centers characterized by fractional crystallization and magma mixing processes at different lower- to mid-crustal levels as well as by assimilation of mafic lower crust and/or its partial melts. Thus, we propose that the adakite-like signatures of Ecuadorian rocks (e.g., high Sr/Y and La/Yb values) are primarily the result of lower- to mid-crustal processing of mantle-derived melts, rather than of slab melts and slab melt–mantle interactions. The isotopic signatures of the least evolved adakite-like rocks of the active and recent volcanoes are the same as those of Tertiary ”normal” calc-alkaline magmatic rocks of Ecuador suggesting that the source of the magma did not change through time. What changed was the depth of magmatic evolution, probably as a consequence of increased compression induced by the stronger coupling between the subducting and overriding plates associated with subduction of the aseismic Carnegie Ridge.     

Influence of atmospheric circulation patterns on local cloud and solar variability in Bergen,Norway

In a previous paper, we have shown that long-term cloud and solar observations (1965–2013) in Bergen, Norway (60.39°N, 5.33°E) are compatible with a largely cloud dominated radiative climate. Here, we explicitly address the relationship between the large scale circulation over Europe and local conditions in Bergen, identifying specific circulation shifts that have contributed to the observed cloud and solar variations. As a measure of synoptic weather patterns, we use the Grosswetterlagen (GWL), a daily classification of European weather for 1881–2013. Empirical models of cloud cover, cloud base, relative sunshine duration, and normalised global irradiance are constructed based on the GWL frequencies, extending the observational time series by more than 70 years. The GWL models successfully reproduce the observed increase in cloud cover and decrease in solar irradiance during the 1970s and 1980s. This cloud-induced dimming is traced to an increasing frequency of cyclonic and decreasing frequency of anticyclonic weather patterns over northern Europe. The changing circulation patterns in winter can be understood as a shift from the negative to the positive phase of the North Atlantic and Arctic Oscillation. A recent period of increasing solar irradiance is observed but not reproduce by the GWL models, suggesting this brightening is associated with factors other than large scale atmospheric circulation, possibly decreasing aerosol loads and local cloud shifts.