Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi‐objective genetic algorithm

Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.     

Sensitivity of interfacial hydraulics to the microtopographic roughness of water‐lain gravels

Flow within the interfacial layer of gravel‐bed rivers is poorly understood, but this zone is important because the hydraulics here transport sediment, generate flow structures and interact with benthic organisms. We hypothesized that different gravel‐bed microtopographies generate measurable differences in hydraulic characteristics within the interfacial layer. This was tested using a high density of spatially and vertically distributed, velocity time series measured in the interfacial layers above three surfaces of contrasting microtopography. These surfaces had natural water‐worked textures, captured in the field using a casting procedure. Analysis was repeated for three discharges, with Reynolds numbers between 165000 and 287000, to evaluate whether discharge affected the impact of microtopography on interfacial flows. Relative submergence varied over a small range (3.5 to 8.1) characteristic of upland gravel‐bed rivers. Between‐surface differences in the median and variance of several time‐averaged and turbulent flow parameters were tested using non‐parametric statistics. Across all discharges, microtopographic differences did not affect spatially averaged (median) values of streamwise velocity, but were associated with significant differences in its spatial variance, and did affect spatially averaged (median) turbulent kinetic energy. Sweep and ejection events dominated the interfacial region above all surfaces at all flows, but there was a microtopographic effect, with Q2 and Q4 events less dominant and structures less persistent above the surface with the widest relief distribution, especially at the highest Reynolds number flow. Results are broadly consistent with earlier work, although this analysis is unique because of the focus on interfacial hydraulics, spatially averaged ‘patch scale’ metrics and a statistical approach to data analysis. An important implication is that observable differences in microtopography do not necessarily produce differences in interfacial hydraulics. An important observation is that appropriate roughness parameterizations for gravel‐bed rivers remain elusive, partly because the relative contributions to flow resistance of different aspects of bed microtopography are poorly constrained. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Organohalogen contaminants and trace metals in North-East Atlantic porbeagle shark (Lamna nasus)

The North-East Atlantic porbeagle (Lamna nasus) population has declined dramatically over the last few decades and is currently classified as ‘Critically Endangered’. As long-lived, apex predators, they may be vulnerable to bioaccumulation of contaminants. In this study organohalogen compounds and trace elements were analysed in 12 specimens caught as incidental bycatch in commercial gillnet fisheries in the Celtic Sea in 2011. Levels of organohalogen contaminants were low or undetectable (summed CB and BDE concentrations 0.04–0.85 mg kg⁻¹ wet weight). A notably high Cd concentration (7.2 mg kg⁻¹ wet weight) was observed in one mature male, whereas the range observed in the other samples was much lower (0.04–0.26 mg kg⁻¹ wet weight). Hg and Pb concentrations were detected only in single animals, at 0.34 and 0.08 mg kg⁻¹ wet weight, respectively. These contaminant levels were low in comparison to other published studies for shark species.     

A Statistical Framework for Calculating and Assessing Compositional Linear Trends Within Fault Zones: A Case Study of the NE Block of the Clark Segment,San Jacinto Fault,California, USA

Utilizing chemical data derived from the various fault zone architectural components of the Clark strand of the San Jacinto fault, southern California, USA, we apply for the first time non-central principal component analysis to calculate a compositional linear trend within molar A–CN–K space. In this procedure A–CN–K are calculated as the molar proportions of Al₂O₃ (A), CaO* + Na₂O (CN), and K₂O (K) in the sum of molar Al₂O₃, Na₂O, CaO*, and K₂O. CaO* is the molar CaO after correction for apatite. We then derive translational invariant chemical alteration intensity factors, t, for each architectural component through orthogonal projection of analyzed samples onto the compositional linear trend. The chemical alteration intensity factor t determines the relative change in composition compared to the original state (i.e., the composition of the altered wall rocks). It is dependent on the degree of intensity to which the process or processes responsible for the change in composition of each architectural component has been active. These processes include shearing, fragmentation, fluid flow, and generation of frictional heat. Non-central principal component analysis indicates that principal component 1 explains 99.7 % of the spread of A–CN–K data about the calculated compositional linear trend (i.e., the variance). The significance level for the overall one-way analysis of variance (ANOVA) is 0.0001. Such a result indicates that at least one significant difference across the group of means of t values is different at the 95 % confidence level. Following completion of the overall one-way ANOVA, the difference in means t test indicated that the mean of the t values for the fault core are different than the means obtained from the transition and damage zones. In contrast, at the 95 % confidence level, the means of the t values for the transition and damage zones are not statistically distinguishable. The results of XRD work completed during this study revealed that the <2 µm fraction is composed primarily of illite/smectite with ~15 % illite in the damage zone, of illite/smectite with ~30 % illite in the transition zone, and of discreet illite with very minor smectite in the fault core. These changes parallel the increasing values of the chemical alteration intensity factors (i.e., t). Based on the above results, it is speculated that when fault zones are derived from tonalitic wall rocks at depths of ~400 ± 100 m, the onset of the illite/smectite to illite conversion will occur when t values exceed 0.20 ± 0.12, the average chemical alteration intensity factor calculated for the transition zone. Under such conditions during repeated rupturing events, frictional heat is produced and acidic fluids with elevated temperatures (≥ ~125 °C) are flushed through the fault core. Over time, the combination of shearing, fragmentation, and frictionally elevated temperatures eventually overcomes the kinetic barrier for the illite/smectite to illite transition. Such settings and processes are unique to fault zones, and as a result, they represent an underappreciated setting for the development of illite from illite/smectite. The success of non-central principal component analysis in this environment offers the first statistically rigorous methodology for establishing the existence of compositional linear trends in fault zones. This method also derives quantifiable alteration intensity factors that could potentially be used to compare the intensity of alteration at different segments of a fault, as well as offer a foundation to interpret the potential driving forces for said alteration and differences therein.     

Riparian responses to reduced flood flows: comparing and contrasting narrowleaf and broadleaf cottonwoods

Abstract

To enable assessment of risks of water management to riparian ecosystems at a regional scale, we developed a quantile-regression model of abundance of broadleaf cottonwoods (Populus deltoides and P. fremontii) as a function of flood flow attenuation. To test whether this model was transferrable to narrowleaf cottonwood (Populus angustifolia), we measured narrowleaf abundance along 39 river reaches in northwestern Colorado, USA. The model performed well for narrowleaf in all 32 reaches where reservoir storage was <75% of mean annual flow. Field data did not fit the model at four of seven reaches where reservoir storage was >90% of mean annual flow. In these four reaches, narrowleaf was abundant despite peak flow attenuation of 45–61%. Poor model performance in these four reaches may be explained in part by a pulse of narrowleaf cottonwood expansion as a response to channel narrowing and in part by differences between narrowleaf and broadleaf cottonwood response to floods and drought.

Editor Z.W. Kundzewicz; Guest editor M. Acreman

Citation Wilding, T.K., Sanderson, J.S., Merritt, D.M., Rood, S.B., and Poff, N.L., 2014. Riparian responses to reduced flood flows: comparing and contrasting narrowleaf and broadleaf cottonwoods. Hydrological Sciences Journal, 59 (3–4), 605–617.     

On‐Site Vapor‐Phase Analysis as a Novel Approach for Monitoring Groundwater Wells

The results of comprehensive field testing of on‐site vapor‐phase‐based groundwater monitoring methods are presented to demonstrate their utility as a robust and cost‐effective approach for rapidly obtaining volatile organic compounds (VOCs) concentration data from a monitoring well. These methods—which rely on sensitive, commercially available field equipment to analyze vapor in equilibrium with groundwater—proved easy to implement and can be tailored to site‐specific needs, including multilevel sampling. During field testing, low‐flow groundwater concentrations could be reasonably estimated using submerged passive vapor diffusion samplers or field equilibration of collected groundwater (R² = 0.85 to 0.96). These two methods are not as reliant on in‐well mixing to overcome vertical stratification within wells as simpler headspace methods. The importance of well and aquifer‐specific factors on concentration data (and therefore method selection) is highlighted, including the effect of changing in‐well patterns due to seasonal temperature gradients. Results indicated that vertical stratification was relatively limited within the set of wells included in these studies, resulting in similar performance for short depth‐discrete passive vapor diffusion samplers (constructed from 40‐mL vials) and longer samplers (2.5 to 5 feet in length) designed to cover a larger portion of the screened interval. A year‐long, multi‐event evaluation demonstrated that vapor‐phase‐based monitoring methods are no more variable than conventional groundwater monitoring methods, with both types subject to similar spatial and temporal variability that can be difficult to reduce. Vapor sampling methods represent a promising approach for estimation of groundwater concentrations by reducing the cost liabilities associated with monitoring while providing a more sustainable approach.     

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.     

基于分布式大流域径流模型的中国西北黑河流域水文模拟

水资源短缺是中国西北干旱地区长期的问题,区域人口增加、城市化扩张,加之气候变化的影响进一步加剧了西北地区水资源短缺,也使生活用水、灌溉用水、工业用水和维持生态系统稳定的用水危险加剧.采用分布式大流域径流模型(DLBRM)模拟黑河流域水文(中国第二大内陆河,流域面积128 000 km2)来理解区域的冰川和积雪融化水、地下水、地表水、蒸散发等方面的分布,评估气候变化对水文的影响和冰川退缩对中游和下游来水量的影响.模拟结果表明,黑河流域的大部分产流那源于黑河上游地区的祁连山.模拟1990-2000年黑河河流日流量变化结果认为,黑河中游正义峡给下游的供水为10×108m3,其中地表径流占51%,层间流占49%.中游地区沙土具有较高的蒸腾发能力,近一半的地表水被蒸发掉.模拟实践证明,分布式大流域径流模型可以结合气候变化、水资源管理方面的成果,改进流域水文模拟的精度.