B and Li isotopes as intrinsic tracers for injection tests in aquifer storage and recovery systems

Boron and Li isotopes have been tested as environmental tracers of treated sewage injected into a sandy aquifer (Shafdan reclamation project, Israel). During a 38 days injection test in a newly dug injection well, a conservative artificial tracer (Br⁻) was monitored together with δ¹¹B and δ⁷Li in the injectate, in the unsaturated soil zone (porous cup) and an observation well in the aquifer. In spite of B and Li concentrations in the injectate close to background values, significant shifts of the isotope signatures could be observed over the duration of the injection test. Boron isotope ratios show a breakthrough curve delayed with respect to Br⁻ breakthrough due to some reversible sorption on the aquifer material. No isotope fractionation was observed in the unsaturated or the saturated zone so that B isotopes can be considered as conservative in the investigated part of the aquifer system. Lithium isotopes are strongly fractionated, probably due to sorption processes. Lithium concentrations point to a Li sink in the system, δ⁷Li values vary strongly with a tendency of ⁷Li depletion in the liquid phase over the duration of the experiment. This is opposite to the expected preferential sorption of ⁶Li onto clay minerals. Boron isotopes reveals a valuable tracer of artificial recharge of freshwaters derived from treated sewage, both for short term tracer tests and for long-term monitoring of artificial recharge, even if in aquifers with higher clay contents, sorption-linked isotope fractionation cannot be excluded. More data are needed on Li isotope fractionation in natural groundwater systems to assess the potential of this tracer as monitoring tool.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Deglaciation history at the Alpine-Mediterranean transition (Argentera-Mercantour,SW Alps) from 10Be dating of moraines and glacially polished bedrock

Estimating the extent and age of the last glacial maxima as well as the chronology of glacial recessions in various environmental contexts is key to source-to-sink studies and paleoclimate reconstructions. The Argentera-Mercantour massif is located at the transition between the Alps and the Mediterranean Sea, therefore, its deglaciation chronology can be compared to the sediment budget of the Var River basin. Based on 13 new cosmic-ray exposure (CRE) beryllium-10 (¹⁰Be) datings performed on moraines and polished crystalline bedrocks and 22 reassessed ¹⁰Be CRE ages from similar altitude nearby steep basement surfaces, and from a lake sediment core, we can constrain the deglaciation chronology of the Argentera-Mercantour massif. These data allow for the first time to fully reconstruct the deglaciation history at the scale of the entire massif in agreement with a major glacier recession at c. 15 ka, at the onset of Bølling transition between the Oldest and Older Dryas. Main deglaciation of the upper slopes [2700–2800 m above sea level (a.s.l.)] occurred after the Last Glacial Maximum (LGM) at 20.8–18.6 ka, followed by the main deglaciation of the lower slopes (2300 m a.s.l.) at 15.3–14.2 ka. Finally, the flat polished surfaces above 2600 m a.s.l. and the zones confined within narrow lateral valleys were likely affected by progressive ice melting of remaining debris covered glaciers and moraine erosion following the Younger Dryas re-advance stage between 12 and 8–9 ka. At lower elevations, the Vens Lake located at 2300 m a.s.l., allows evidence of the onset of lake sedimentation at c. 14 ka and a transition towards a vegetated environment that mainly occurred before 8 ka. Moraine final stabilization at 5 ka might reflect denudation acceleration during the Holocene humid phase. This contribution reveals a glacier–climate relationship more sensitive to warming phases in the southern Alps highlighted by a major decrease of glaciers after c. 15 ka. This major deglaciation is correlated with a 2.5-fold decrease of sediment discharge of rivers into the Mediterranean Sea. © 2019 John Wiley & Sons, Ltd.     

The April 2010 North African heatwave: when the water vapor greenhouse effect drives nighttime temperatures

Climate Dynamics - North Africa experienced a severe heatwave in April 2010 with daily maximum temperatures ( $$T_{max}$$ ) frequently exceeding $$40\,^{\circ }\mathrm{C}$$ and daily minimum...     

Geothermal waters from the Taupo Volcanic Zone,New Zealand: Li,B and Sr isotopes characterization

Chemical and isotopic data for 23 geothermal water samples collected in New Zealand within the Taupo Volcanic Zone (TVZ) are reported. Major and trace elements including Li, B and Sr and their isotopic compositions (δ⁷Li, δ¹¹B, ⁸⁷Sr/⁸⁶Sr) were determined in high temperature geothermal waters collected from deep boreholes in different geothermal fields (Ohaaki, Wairakei, Mokai, Kawerau and Rotokawa geothermal systems). Lithium concentrations are high (from 4.5 to 19.9 mg/L) and Li isotopic compositions (δ⁷Li) are homogeneous, ranging between −0.5‰ and +1.4‰. In particular, it is noteworthy that, except for the samples from the Kawerau geothermal field having slightly higher δ⁷Li values (+1.4%), the other geothermal waters have a near constant δ⁷Li signature around a mean value of 0‰ ± 0.6 (2σ, n = 21). Boron concentrations are also high and relatively homogeneous for the geothermal samples, falling between 17.5 and 82.1 mg/L. Boron isotopic compositions (δ¹¹B) are all negative, and display a range between −6.7‰ and −1.9‰. These B isotope compositions are in agreement with those of the Ngawha geothermal field in New Zealand. Lithium and B isotope signatures are in a good agreement with a fluid signature mainly derived from water/rock interaction involving magmatic rocks with no evidence of seawater input. On the other hand, Sr concentrations are lower and more heterogeneous and fall between 2 and 165 μg/L. The ⁸⁷Sr/⁸⁶Sr ratios range from 0.70549 to 0.70961. These Sr isotope compositions overlap those of the Rotorua geothermal field in New Zealand, confirming that some geothermal waters (with more radiogenic Sr) have interacted with bedrocks from the metasedimentary basement. Each of these isotope systems on their own reveals important information about particular aspects of either water source or water/rock interaction processes, but, considered together, provide a more integrated understanding of the geothermal systems from the TVZ in New Zealand.     

Relay zone evolution: a history of repeated fault propagation and linkage,central Corinth rift,Greece

Established models indicate that, before being breached, relay zones along rift borders can evolve either by lengthening and rotating during progressive overlap of growing fault segments (isolated fault model), or, by simply rotating without lengthening before breaching (coherent fault model). The spatio‐temporal distribution of vertical motions in a relay zone can thus be used to distinguish fault growth mechanisms. Depositional relay zones that develop at sea level and accommodate both deposition on the ramp itself as well as transfer of sediments from the uplifting footwall into the hangingwall depocentres and provide the most complete record of vertical motions. We examine the development of a depositional relay ramp on the border of the active Corinth rift, Greece to reconstruct fault interaction in time and space using both onshore and offshore (2D seismic lines) data. The Akrata relay zone developed over a period of ca. 0.5 Myr since the Middle Pleistocene between the newly forming East Helike Fault (EHF) that propagated towards the older, more established Derveni Fault (DF). The relay zone captured the Krathis River, which deposited prograding Gilbert‐type deltas on the sub‐horizontal ramp. Successive oblique faults record progressive linkage and basinward migration of accommodation along the ramp axis, whereas marine terraces record diachronous uplift in their footwalls. Although early linkage of the relay zone occurs, continuous propagation and linkage of the EHF onto the static DF is recorded before final beaching. Rotation on forced folds above the upward and laterally propagating normal faults at the borders of the relay zone represents the ramp hinges. The Akrata relay zone cannot be compared directly to a simple fault growth model because (1) the relay zone connects two fault segments of different generations; (2) multiple linkages during propagation was facilitated by the presence of pre‐existing crustal structures, inherited from the Hellenide fold and thrust belt. The linkage of the EHF to the DF contributed to the westward and northward propagation of the southern rift border.     

Nutrient source attribution: Quantitative typology distinction of active and legacy source contributions to waterborne loads

Distinction between active and legacy sources of nutrients is needed for effective reduction of waterborne nutrient loads and associated eutrophication. This study quantifies main typological differences in nutrient load behaviour versus water discharge for active and legacy sources. This quantitative typology is used for source attribution based on monitoring data for water discharge and concentrations of total nitrogen (TN) and total phosphorous (TP) from 37 catchments draining into the Baltic Sea along the coastline of Sweden over the period 2003–2013. Results indicate dominant legacy source contributions to the monitored loads of TN and TP in most (33 of the total 37) study catchments. Dominant active sources are indicated in 1 catchment for TN, and mixed sources are indicated in 3 catchments for TN, and 4 catchments for TP. The TN and TP concentration contributions are quantified to be overall higher from the legacy than the active sources. Legacy concentrations also correlate well with key indicators of human activity in the catchments, agricultural land share for TN (R² = 0.65) and population density for TP (R² = 0.56). Legacy-dominated nutrient concentrations also change more slowly than in catchments with dominant active or mixed sources. Various data-based results and indications converge in indicating legacy source contributions as largely dominant, mainly anthropogenic, and with near-zero average change trends in the present study of catchments draining into the Baltic Sea along the coastline of Sweden, as in other parts of the world. These convergent indications emphasize needs to identify and map the different types of sources in each catchment, and differentiate strategies and measures to target each source type for possible achievement of shorter- and longer-term goals of water quality improvement.