The origin of reaction textures in mantle peridotite xenoliths from Sal Island, Cape Verde: the case for “metasomatism” by the host lava

Reaction zones around minerals in mantle xenoliths have been reported from many localities worldwide. Interpretations of the origins of these textures fall into two groups: mantle metasomatic reaction or reaction during transport of the xenoliths to the surface. A suite of harzburgitic mantle xenoliths from Sal, Cape Verde show clear evidence of reaction during transport. The reactions resulted in the formation of olivine–clinopyroxene and Si- and alkali-rich glass reaction zones around orthopyroxene and sieve-textured clinopyroxene and sieve textured spinel, both of which are associated with a Si- and alkali-rich glass similar to that in the orthopyroxene reaction zones. Reaction occurred at pressures less than the mantle equilibration pressure and at temperatures close to the liquidus temperature of the host magma. In addition, there is a clear spatial relation of reaction with the host lava: reaction is most intense near the lava/xenolith contact. The residence time of the xenoliths in the host magma, determined from Fe–Mg interdiffusion profiles in olivine, was approximately 4 years. Our results cannot be reconciled with a recent model for the evolution of the mantle below the Cape Verde Archipelago involving mantle metasomatism by kimberlitic melt. We contend that alkali-rich glasses in the Sal xenoliths are not remnants of a kimberlitic melt, but rather they are the result of reaction between the host lava or a similar magma and xenolith minerals, in particular orthopyroxene. The formation of a Si- and alkali-rich glass by host magma–orthopyroxene reaction appears to be a necessary precursor to formation of sieve textured spinel and clinopyroxene.     

Rapid biotic molecular transformation of fulvic acids in a karst aquifer

The study of molecular transformation processes of dissolved organic carbon (DOC) in the environment significantly contributes to a better understanding of the global biogeochemical organic matter cycle. In an oxic karst groundwater system, in which the most powerful abiotic DOC degradative reactions, photodegradation and metal-mediated redox chemistry, are at best marginal contributors, a near complete turnover of fulvic acids (FAs) has been observed within decades (∼60 years). Depletion of oxygen for a very extensive range of aliphatic and aromatic carbon chemical environments has been confirmed as well as the formation of novel classes of compounds, suggesting a major contribution from biotic processes. From these results we infer that FAs must be perceived as a rather active participant in the global carbon cycle. Molecular-level alterations of such magnitude and rapidity on such short-time scales ought to be considered as widespread in the processing of “refractory” DOC in the environment.     

Estimating uncertainty associated with water stages from a single SAR image

It is the goal of remote sensing to infer information about objects or a natural process from a remote location. This invokes that uncertainty in measurement should be viewed as central to remote sensing. In this study, the uncertainty associated with water stages derived from a single SAR image for the Alzette (G.D. of Luxembourg) 2003 flood is assessed using a stepped GLUE procedure. Main uncertain input factors to the SAR processing chain for estimating water stages include geolocation accuracy, spatial filter window size, image thresholding value, DEM vertical precision and the number of river cross sections at which water stages are estimated. Initial results show that even with plausible parameter values uncertainty in water stages over the entire river reach is 2.8 m on average. Adding spatially distributed field water stages to the GLUE analysis following a one-at-a-time approach helps to considerably reduce SAR water stage uncertainty (0.6 m on average) thereby identifying appropriate value ranges for each uncertain SAR water stage processing factor. For the GLUE analysis a Nash-like efficiency criterion adapted to spatial data is proposed whereby acceptable SAR model simulations are required to outperform a simpler regression model based on the field-surveyed average river bed gradient. Weighted CDFs for all factors based on the proposed efficiency criterion allow the generation of reliable uncertainty quantile ranges and 2D maps that show the uncertainty associated with SAR-derived water stages. The stepped GLUE procedure demonstrated that not all field data collected are necessary to achieve maximum constraining. A possible efficient way to decide on relevant locations at which to sample in the field is proposed. It is also suggested that the resulting uncertainty ranges and flood extent or depth maps may be used to evaluate 1D or 2D flood inundation models in terms of water stages, depths or extents. For this, the extended GLUE approach, which copes with the presence of uncertainty in the observed data, may be adopted.     

Anisotropic dual-continuum representations for multiscale poroelastic materials: Development and numerical modelling

Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations.     

Causes and consequences of long-term groundwater overabstraction in Jordan

Hydrogeology Journal - In 2017, a comprehensive review of groundwater resources in Jordan was carried out for the first time since 1995. The change in groundwater levels between 1995 and 2017 was...     

A data assimilation approach to discharge estimation from space

River discharge is currently monitored by a diminishing network of gauges, which provide a spatially incomplete picture of global discharges. This study assimilated water level information derived from a fused satellite Synthetic Aperture Radar (SAR) image and digital terrain model (DTM) with simulations from a coupled hydrological and hydrodynamic model to estimate discharge in an un‐gauged basin scenario. Assimilating water level measurements led to a 79% reduction in ensemble discharge uncertainty over the coupled hydrological hydrodynamic model alone. Measurement bias was evident, but the method still provided a means of improving estimates of discharge for high flows. The study demonstrates the potential of currently available synthetic aperture radar imagery to reduce discharge uncertainty in un‐gauged basins when combined with model simulations in a data assimilation framework, where sufficient topographic data are available. The work is timely because in the near future the launch of satellite radar missions will lead to a significant increase in the volume of data available for space‐borne discharge estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

Thermal influence of urban groundwater recharge from stormwater infiltration basins

Groundwater warming below cities has become a major environmental issue; but the effect of distinct local anthropogenic sources of heat on urban groundwater temperature distributions is still poorly documented. Our study addressed the local effect of stormwater infiltration on the thermal regime of urban groundwater by examining differences in water temperature beneath stormwater infiltration basins (SIB) and reference sites fed exclusively by direct infiltration of rainwater at the land surface. Stormwater infiltration dramatically increased the thermal amplitude of groundwater at event and season scales. Temperature variation at the scale of rainfall events reached 3 °C and was controlled by the interaction between runoff amount and difference in temperature between stormwater and groundwater. The annual amplitude of groundwater temperature was on average nine times higher below SIB (range: 0·9–8·6 °C) than at reference sites (range: 0–1·2 °C) and increased with catchment area of SIB. Elevated summer temperature of infiltrating stormwater (up to 21 °C) decreased oxygen solubility and stimulated microbial respiration in the soil and vadose zone, thereby lowering dissolved oxygen (DO) concentration in groundwater. The net effect of infiltration on average groundwater temperature depended upon the seasonal distribution of rainfall: groundwater below large SIB warmed up (+0·4 °C) when rainfall occurred predominantly during warm seasons. The thermal effect of stormwater infiltration strongly attenuated with increasing depth below the groundwater table indicating advective heat transport was restricted to the uppermost layers of groundwater. Moreover, excessive groundwater temperature variation at event and season scales can be attenuated by reducing the size of catchment areas drained by SIB and by promoting source control drainage systems. Copyright © 2009 John Wiley & Sons, Ltd.     

Large-strain deformation and strain partitioning in polyphase rocks: Dislocation creep of olivine–magnesiowüstite aggregates

Aggregates composed of olivine and magnesiowüstite have been deformed to large strains at high pressure and temperature to investigate stress and strain partitioning, phase segregation and possible localization of deformation in a polyphase material. Samples with 20 vol.% of natural olivine and 80 vol.% of (Mg_0.7Fe_0.3)O were synthesized and deformed in a gas-medium torsion apparatus at temperatures of 1127 °C and 1250 °C, a confining pressure of 300 MPa and constant angular displacement rates equivalent to constant shear strain rates of 1–3.3 × 10^− 4 s^− 1. The samples deformed homogeneously to total shear strains of up to γ  15. During constant strain rate measurements the flow stress remained approximately stable at 1250 °C while it progressively decreased after the initial yield stress at the lower temperature. Mechanical data, microstructures and textures indicate that both phases were deforming in the dislocation creep regime. The weaker component, magnesiowüstite, controlled the rheological behavior of the bulk material and accommodated most of the strain. Deformation and dynamic recrystallization lead to grain refinement and to textures that were not previously observed in pure magnesiowüstite and may have developed due to the presence of the second phase. At 1127 °C, olivine grains behaved as semi-rigid inclusions rotating in a viscous matrix. At 1250 °C, some olivine grains remained largely undeformed while deformation and recrystallization of other grains oriented for a-slip on (010) resulted in a weak foliation and a texture typical for pure dry olivine aggregates. Both a-slip and c-slip on (010) were activated in olivine even though the nominal stresses were up to 2 orders of magnitude lower than those needed to activate these slip systems in pure olivine at the same conditions.