Steady- and transient-state inversion in hydrogeology by successive flux estimation

A calibration method to solve the groundwater inverse problem under steady- and transient-state conditions is presented. The method compares kriged and numerical head field gradients to modify hydraulic conductivity without the use of non-linear optimization techniques. The process is repeated iteratively until a close match with piezometric data is reached. The approach includes a damping factor to avoid divergence and oscillation of the solution in areas of low hydraulic gradient and a weighting factor to account for temporal head variation in transient simulations. The efficiency of the method in terms of computing time and calibration results is demonstrated with a synthetic field. It is shown that the proposed method provides parameter fields that reproduce both hydraulic conductivity and piezometric data in few forward model solutions. Stochastic numerical experiments are conducted to evaluate the sensitivity of the method to the damping function and to the head field estimation errors.     

Bias in hydraulic head measurements from multilevel vibrating-wire piezometers with excessively permeable backfill

Hydrogeology Journal - An extensive network of multilevel vibrating-wire piezometers (VWP) was recently created to monitor the spatial and temporal variation of pore pressure (and hydraulic head)...     

Defining the Potential of Nanoscale Re‐Os Isotope Systematics Using Atom Probe Microscopy

Atom probe microscopy (APM) is a relatively new in situ tool for measuring isotope fractions from nanoscale volumes (< 0.01 μm³). We calculate the theoretical detectable difference of an isotope ratio measurement result from APM using counting statistics of a hypothetical data set to be ± 4δ or 0.4% (2s). However, challenges associated with APM measurements (e.g., peak ranging, hydride formation and isobaric interferences), result in larger uncertainties if not properly accounted for. We evaluate these factors for Re‐Os isotope ratio measurements by comparing APM and negative thermal ionisation mass spectrometry (N‐TIMS) measurement results of pure Os, pure Re, and two synthetic Re‐Os‐bearing alloys from Schwander et al. (2015, Meteoritics and Planetary Science, 50, 893) [the original metal alloy (HSE) and alloys produced by heating HSE within silicate liquid (SYN)]. From this, we propose a current best practice for APM Re‐Os isotope ratio measurements. Using this refined approach, mean APM and N‐TIMS ¹⁸⁷Os/¹⁸⁹Os measurement results agree within 0.05% and 2s (pure Os), 0.6–2% and 2s (SYN) and 5–10% (HSE). The good agreement of N‐TIMS and APM ¹⁸⁷Os/¹⁸⁹Os measurements confirms that APM can extract robust isotope ratios. Therefore, this approach permits nanoscale isotope measurements of Os‐bearing alloys using the Re‐Os geochronometer that could not be measured by conventional measurement principles.     

Methane-derived authigenic carbonates from active hydrocarbon seeps of the St. Lawrence Estuary,Canada

Nearly 2000 pockmarks with diameters ranging from a few tens of meters up to 700 m are present on the seafloor of the St. Lawrence Estuary in eastern Canada. Coring of some pockmarks resulted in the recovery of various-sized and shaped carbonate concretions in a predominantly silty mud matrix. Petrographic and geochemical data on four authigenic carbonate concretions are reported as well as data from shell material in the unconsolidated sediment. Video observations and echo-sounder images indicate that the sampled pockmarks are actively gas venting. The video images show significant look-alike microbial mats in areas where gas is venting. The carbonate concretions are primarily made up of carbonate cements with varying percentage of shell fragments, micrite particles and fine-grained clastics. Orthorhombic crystal morphology and diagenetic fabrics including isopachous layers and botryoids characterize the aragonite cement. Oxygen isotopes ratios for the cement crusts do not record any thermal anomaly at the site of precipitation with δ¹⁸O_VPDB ratios (+3‰) in equilibrium with cold (5 °C) deep marine waters, whereas significant negative δ¹³C_VPDB ratios (−9.9 to −33.5‰) for cement and shell material within concretions indicate that the carbonates largely derive from the microbial oxidation of methane. The δ¹³C_VPDB ratios of aragonite shells (−2.7 to −5.6‰) taken from unconsolidated sediments at some distance from the concretions/vents show variable dilution of HCO₃⁻ with negative δ¹³C_VPDB ratios derived from microbial oxidation of methane with isotopically normal (0‰) marine bicarbonate. These results are in agreement with other lines of evidence suggesting that pockmarks formed through the recent and still active release of gas from a reservoir within the Paleozoic sedimentary succession.     

Micro-fracture instabilities in granular solids

We study the aspect of unstable behavior (like strain localization bands) in elastic solids as a consequence of micro-fracturing. A two-scale approach of computational homogenization is considered. The macroscopic behavior is investigated by finite element computations on a unit cell. At the micro-level, we consider a granular structure with elastic grains. The inter-granular boundaries are modeled with cohesive laws, friction and unilateral contact. We show that decohesion between grains gives rise to macro-instabilities, indicated by the loss of ellipticity, typical for deformation localization bands. The relation between the microscopic softening on inter-granular boundaries and the onset of macro-instabilities is studied through numerical examples. The influence of the cohesive law and friction parameters is analyzed. For periodic distributions of granular structures, we prove the loss of periodicity by failure and the corresponding size dependence effect in the homogenized response. We present numerical examples of bifurcation of solutions for granular cell structures and of particular solutions specific to elementary volumes with periodic cell distribution. Size dependence appears in the unstable regime and is strongly influenced by cohesion and friction parameters.     

Chemical mass transfer in the Rio Tanama system,west-central Puerto Rico

Equilibrium relationships are defined between stream waters and weathering products, kaolinite and calcium montmorillonite, for the Rio Tanama system, west-central Puerto Rico. The major element composition of 46 water samples of springs and streams define a reaction path in the system CaO-Na₂O-MgO-Al₂O₃-SiO₂-H₂O between acid waters containing low concentrations of alkali cations and detrital reactant minerals. The principal reactant phases appear to be chlorite, plagioclase and orthoclase and occasionally anhydrite or calcite. Headward erosion by the Rio Tanama supplies the reactant phases to the stream silt load.The chemical denudation rate calculated for the Rio Tanama system is about 30 m/million yr. The chemical stream load appears to be buffered by the product phases in the main river over the 15–20 km river length sampled in this study.The silt and soil mineralogy and water compositions are used to define a log K for the hydrolysis of Ca-montmorillonite at 25°C of 35.0 ± 0.8. This value is in reasonable agreement with the value of 37.1 ± 1.0 defined by Garrels and Mackenzie (1967) in a similar manner for spring waters in the Sierra Nevada.     

Morphology,structure and kinematics of a rainfall controlled slow‐moving Andean landslide,Peru

The large slow‐moving landslide of Maca is located in the upper Colca valley (southern Peru), a region characterized by a well pronounced rainy period, and intense and recurrent sustained seismicity. The landslide, developed in deep lacustrine deposits, has recently accelerated, threatening the Maca village. This work aims at understanding the rupture mechanism and the causes of the recent landslide reactivation/acceleration. We present a multidisciplinary characterization of the Maca landslide that includes: (i) geological and morphological mapping in the field; (ii) remote sensing analysis using an historical aerial photograph of 1955 and the Pléiades satellite images (2013); (iii) global positioning system (GPS) including time‐series of surveys over 13 years, and continuous measurements over 14 months; (iv) a geophysical campaign with deep electrical resistivity tomography profiles acquired across the landslide mass. Our study shows that this 60 Mm³ landslide, which can be classified as a clay/silt compound landslide, moved by 15 m between 2001 and 2014 with a large inter‐annual velocity variation (up to a factor of 500) depending on the rainfall intensity. We suggest that these dramatic changes in velocity are the result of the combination of a threshold mechanism and the short intense rainy season in Peru. This study reveals three main driving factors acting at different timescales: (i) over several decades, the river course has significantly changed, causing the Maca landslide reactivation in the 1980s due to the erosion of its toe; (ii) at the year scale, a minimum amount of rainfall is required to trigger the motion and this amount controls the landslide velocity; (iii) transient changes in slide velocity may occur anytime due to earthquakes. This study particularly highlights the non‐linear behaviour of the motion with rainfall. Copyright © 2016 John Wiley & Sons, Ltd.