Experiment and Simulation of Non-Reactive Solute Transport in Porous Media

To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater.     

Closed Analytic Elements with Flexible Geometry

A simple and fast treatment of hydrogeologic features with irregularly shaped boundaries in two‐dimensional analytic element groundwater flow models is presented. The star domain shapes of the features are restricted to closed shapes represented as smooth and continuous single‐valued functions of distance from a focus point, . The element can be used to treat a variety of boundary and continuity conditions, including those of irregularly shaped lakes or heterogeneities in hydraulic conductivity. The new element is demonstrated via some simple illustrative test cases and shown to be efficient, accurate, and much simpler to implement than existing solutions for irregular shapes.     

A Physically Based Approach for Estimating Hydraulic Conductivity from HPT Pressure and Flowrate

The hydraulic profiling tool (HPT) is widely used to generate profiles of relative permeability vs. depth. In this work, prior numerical modeling results are used to develop a relationship between probe advance rate V (cm/s), probe diameter D (cm), water injection rate Q (mL/min), corrected pressure P_c (psi), and hydraulic conductivity K (feet/d) ((1)) where E is an empirically derived hydraulic efficiency factor. The relationship is validated by 23 HPT profiles that, after averaging K vertically, were similar to slug test results in adjoining monitoring wells. The best fit value of E for these profiles was 2.02. This equation provides a physically based approach for generating hydraulic conductivity profiles with HPT tooling.     

A Simple Analytical Formula for the Leakage Flux Through a Perforated Aquitard

In this methods note, we present a simple analytical formula to quantify the steady‐state leakage flux (Q) over a perforated aquitard. The flux depends on the aquitard thickness (D), the radius of the perforation (R), the hydraulic conductivity of the material inside the perforation (k_fill), the conductivities of the overlying and underlying aquifers (k₁ and k₂, respectively), and the head difference between the two aquifers (ΔH): This equation assumes an aquitard separating two homogeneous and infinite aquifers (R ? aquifer thickness) in which radial flow to and from the perforation occurs, with no other recharge or discharge boundaries near the perforation. The flux through a perforation in a hypothetical case study with D = 10 m, k₁ = 10 m/d, k₂ = 20 m/d, R = 0.072 m, and ΔH = 1 m ranges between less than 1 mL/d if the hole is backfilled with bentonite (k_fill = 10^?4 m/d), to several liters per day if the perforation is backfilled with sand from the overlying aquifer (k_fill = 10 m/d), to several m³/d if the perforation forms an open conduit (k_fill = 10⁵ m/d). The leakage fluxes calculated with this model agree well with those calculated using a numerical model (MODFLOW).     

Comparison of machine learning models for predicting fluoride contamination in groundwater

Groundwater is an especially important freshwater source for water supplies in the Maku area of northwest Iran. The groundwater of the area contains high concentrations of fluoride and is, therefore, important in predicting the fluoride contamination of the groundwater for the purpose of planning and management. The present study aims to evaluate the ability of the extreme learning machine (ELM) model to predict the level of fluoride contamination in the groundwater in comparison to multilayer perceptron (MLP) and support vector machine (SVM) models. For this purpose, 143 water samples were collected in a five-year period, 2004–2008. The samples were measured and analyzed for electrical conductivity, pH, major chemical ions and fluoride. To develop the models, the data set—including Na⁺, K⁺, Ca²⁺ and HCO₃ ^? concentrations as the inputs and fluoride concentration as the output—was divided into two subsets; training/validation (80% of data) and testing (20% of data), based on a cross-validation technique. The radial basis-based ELM model resulted in an R ² of 0.921, an NSC of 0.9071, an RMSE of 0.5638 (mg/L) and an MABE of 0.4635 (mg/L) for the testing data. The results showed that the ELM models performed better than MLP and SVM models for prediction of fluoride contamination. It was observed that ELM models learned faster than the other models during model development trials and the SVM models had the highest computation time.     

A Correction on Coastal Heads for Groundwater Flow Models

We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)h_s ? B/2α, in which h_s is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models.     

PIG: a numerical index for dissemination of groundwater contamination zones

A pollution index of groundwater (PIG) is proposed for quantification of water contamination. PIG quantifies the status of concentrations of water quality measures with respect to their drinking water quality standards. The validity of the proposed index is verified by choosing the data of groundwater quality of the Varaha River Basin (Visakhapatnam District, Andhra Pradesh, India) as a case study. The computed index from the study area varies from 0.83 to 2.55. The index disseminates the area into zones of insignificant (PIG <1.0), low (PIG: 1.0 to 1.5), moderate (PIG: 1.5 to 2.0), high (PIG 2.0 to 2.5) and very high (PIG >2.5) pollution. Insignificant pollution zone is observed from the upstream area, where the groundwater is dominated by , and very high pollution zone from the downstream area, where the groundwater is associated with Cl^?. This indicates that the quality of groundwater in the study area is mainly influenced by the source of geogenic origin, but it is subsequently modified by the effects of anthropogenic and marine sources. Geochemical ratios (Na⁺ : Cl^?, : Cl^?, Na⁺ : Ca²⁺ and Mg²⁺ : Ca²⁺) also form the quantitative basis of the index. The present study paves the way to implement appropriate management strategies at a specific site to circumvent the pollution. As the classification of the pollution zones with PIG depends upon the drinking water quality standards, it becomes a universal assessment tool for groundwater contamination at any test area. Copyright © 2011 John Wiley & Sons, Ltd.     

Approximation of van Genuchten Parameter Ranges from Hydraulic Conductivity Data

The use of retention function and relative conductivity function is essential for the calculation of flow in a variably saturated media using the Richards equation. A widely used mathematical model for this is the Mualem-van Genuchten model which requires the shape parameters and . These, however, are difficult to obtain. When data is scarce, and are often taken from literature and may deviate largely from actual values. The current study presents a novel mathematical model for the approximation of $$ and $$ and for the further estimation of realistic value ranges, which may be used as parameter space, for example, for the calibration of a numerical model. The model was developed for cases where data is scarce and only values of saturated hydraulic conductivity are available. It is based on a large data set from literature and it is demonstrated that the model estimates mean values from that data set with a good accuracy. In order to show the applicability of the model, a second data set has been compiled anew (provided as Supporting Information). The model is incorporated into the current version of the freeware computer program HYPAGS, which enables easy usage.     

Spatio-temporal characteristics of small-scale wave–current ripples on the Ameland ebb-tidal delta

Ebb-tidal deltas are highly dynamic environments affected by both waves and currents that approach the coast under various angles. Among other bedforms of various scales, these hydrodynamics create small-scale bedforms (ripples), which increase the bed roughness and will therefore affect hydrodynamics and sediment transport. In morphodynamic models, sediment transport predictions depend on the roughness height, but the accuracy of these predictors has not been tested for field conditions with strongly mixed (wave–current dominated) forcing. In this study, small-scale bedforms were observed in the field with a 3D Profiling Sonar at five locations on the Ameland ebb-tidal delta, the Netherlands. Hydrodynamic conditions ranged from wave dominated to current dominated, but were mixed most of the time. Small-scale ripples were found on all studied parts of the delta, superimposed on megaripples. Even though a large range of hydrodynamic conditions was encountered, the spatio-temporal variations in small-scale ripple dimensions were relatively small (height 0.015 m, length 0.11 m). Also, the ripples were always highly three-dimensional. These small dimensions are probably caused by the fact that the bed consists of relatively fine sediment. Five bedform height predictors were tested, but they all overestimated the ripple heights, partly because they were not created for small grain sizes. Furthermore, the predictors all have a strong dependence on wave- and current-related velocities, whereas the ripple heights measured here were only related to the near-bed orbital velocity. Therefore, ripple heights and lengths in wave–current-dominated, fine-grained coastal areas ( mm) may be best estimated by constant values rather than values dependent on the hydrodynamics. In the case of the Ameland ebb-tidal delta, these values were found to be m and m. ©2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Investigating young water fractions in a small Mediterranean mountain catchment: Both precipitation forcing and sampling frequency matter

The proportion of water younger than 2–3 months (young water fraction, F_yw) has become increasingly investigated in catchment hydrology. F_yw is typically estimated by comparing seasonal tracer cycles in precipitation and streamflow, through water sampling. However, some open research questions remain, such as: (i) whether part of the summer precipitation should be discarded because the high evapotranspiration demand, (ii) how well F_yw serves as a metric to compare catchments, and (iii) how sampling frequency affects F_yw estimates. To address these questions, we investigated F_yw in soil-, ground- and stream waters for the small Mediterranean Can Vila catchment. Rainfall was sampled at 5-mm intervals. Mobile soil water and groundwater were sampled fortnightly. Stream water was sampled depending on flow at variable time intervals (30 min to 1 week). Over 58 months, this sampling provided 1,529 δ¹⁸O determinations. Isotopic analyses results led us to include summer precipitation in the input signal. We found the highest F_yw in mobile soil waters (34%), while this was almost zero for groundwater except during wet periods. For stream waters, F_yw depended on the discharge variations, so that the flow-weighted young water fraction () was 22.6%, whereas the time-weighted F_yw was just 6.2%. Both and its discharge sensitivity (S_d) varied when different 12-month sampling periods were investigated. The young water fraction that would be obtained from a virtual thorough sampling () was estimated from the S_d and the observed stream flow. This showed an underestimation of by 25% for the frequent dynamic sampling and 66% for weekly sampling, due to missing high flows. Our results confirm that F_yw and its discharge sensitivity are metrics very sensitive to meteorological forcing during the analysed period. Thus, comparisons between catchments need long-term mean annual values and their variability. Our findings also support the dependence of F_yw estimates on the sampling rate and show the advantages of flow-weighted analysis. Finally, catchment water turnover investigations should be accompanied by the analysis of flow duration curves.     

Geomorphology and climate interact to control organic carbon stock and age in mountain river valley bottoms

Organic carbon (OC) in valley bottom downed wood and soil that cycles over short to moderate timescales (10¹ to 10⁵ years) represents a large, dynamic, and poorly quantified pool of carbon whose distribution and residence time affects global climate. We sought to quantify this potentially important OC pool at the watershed scale to estimate its magnitude and age, as well as determine the controls on its variability within watersheds. To do this, we compared four disparate mountain river basins to show that mountain river valley bottoms store substantial estimated OC stocks in floodplain soil and downed wood (median OC of MgC/ha, n = 178). Although soil OC is generally young (exhibiting a median radiocarbon fraction modern value of , n = 121), geomorphic processes regulate soil burial and processes that limit microbial respiration, preserving aged OC in especially deep, unconfined, wet, and/or high-elevation floodplain soils. We statistically modeled OC stocks to show that valley bottom morphology and hydrology regulate variability in floodplain soil retention and resulting variability in OC stock and age in floodplain soil throughout river networks. Comparing the distribution of OC stocks between wood and soil, we find that where floodplain soils are present, their OC stocks are generally greater than OC stocks stored in wood. Our results suggest that although mountain rivers may accumulate large OC stocks relatively rapidly, those stocks are highly sensitive to alterations in soil and wood retention, implying that human alterations to either disturb or restore floodplain wood and soil storage may have substantial impacts on OC storage in river corridors. © 2020 John Wiley & Sons, Ltd.     

A computational workflow for rupture‐to‐structural‐response simulation and its application to Istanbul

Scenario‐based earthquake simulations at regional scales hold the promise in advancing the state‐of‐the‐art in seismic risk assessment studies. In this study, a computational workflow is presented that combines (i) a broadband Green's function‐based fault‐rupture and ground motion simulation—herein carried out using the “UCSB (University of California at Santa Barbara) method”, (ii) a three‐dimensional physics‐based regional‐scale wave propagation simulation that is resolved at  Hz, and (iii) a local soil‐foundation‐structure finite element analysis model. These models are interfaced with each other using the domain reduction method. The innermost local model—implemented in ABAQUS—is additionally enveloped with perfectly matched layer boundaries that absorb outbound waves scattered by the structures contained within it. The intermediate wave propagation simulation is carried out using Hercules , which is an explicit time‐stepping finite element code that is developed and licensed by the CMU‐QUAKE group. The devised workflow is applied to a  km region on the European side of Istanbul, which was modeled using detailed soil stratigraphy data and realistic fault rupture properties, which are available from prior microzonation surveys and earthquake scenario studies. The innermost local model comprises a chevron‐braced steel frame building supported by a shallow foundation slab, which, in turn, rests atop a three‐dimensional soil domain. To demonstrate the utility of the workflow, results obtained using various simplified soil‐structure interaction analysis techniques are compared with those from the detailed direct model. While the aforementioned demonstration has a limited scope, the devised workflow can be used in a multitude of ways, for example, to examine the effects of shallow‐layer soil nonlinearities and surface topography, to devise site‐ and structure‐specific seismic fragilities, and for calibrating regional loss models, to name a few.     

Separation of multi-mode surface waves by supervised machine learning methods

Logistic regression, neural networks and support vector machines are tested for their effectiveness in isolating surface waves in seismic shot records. To distinguish surface waves from other arrivals, we train the algorithms on three distinguishing features of surface-wave dispersion curves in the domain: spectrum coherency of the trace's magnitude spectrum, local dip and the frequency range for a fixed wavenumber k in the spectrum. Numerical tests on synthetic data show that the kernel-based support vector machines algorithm gives the highest accuracy in predicting the surface-wave window in the domain compared to neural networks and logistic regression. This window is also used to automatically pick the fundamental dispersion curve. The other two methods correctly pick the low-frequency part of the dispersion curve but fail at higher frequencies where there is interference with higher-order modes.     

Isotopic and hydrochemical insights into the groundwater characteristics along an arid to semi-humid climate gradient in China

In arid to semi-arid regions, groundwater is a critical water resource heavily relied upon, with the recharge sources and patterns being predominantly shaped by climate change and regional disparities. To compare the characteristics of groundwater in the endorheic and exorheic river basins with the climate transition zone of Gansu Province, this study uses isotopic hydrochemical analyses. This study summarizes the differences in regional groundwater recharge and evolutionary patterns. The results shows that the distribution patterns of precipitation isotopes in endorheic and exorheic river basins are opposite to those of groundwater isotopes. Specifically, the precipitation in the endorheic areas is more depleted in heavy isotopes, whereas the groundwater is more enriched. Both endorheic areas and exorheic areas exhibit similar characteristics of groundwater hydrochemical evolution, evolving from low-mineralization Mg²⁺ HC ${\mathrm{O}}_3^{-}$ recharge water to Na⁺ Cl⁻ type water with saline characteristics. The former is primarily replenished by surface water, whereas the latter is primarily replenished by precipitation. Variations in recharge patterns along with the differences in climatic conditions lead to distinct groundwater conditions in the two regions.     

Starved versus alluvial river bedforms: an experimental investigation

Laboratory experiments were conducted to investigate the formation of river bedforms under sediment supply-limited conditions, i.e. when a motionless substratum is bared by the dynamics of the mobile sediments. Three series of experiments were organized in a laboratory flume by fixing all the hydrodynamic and morphodynamic parameters but varying the thickness of the initial layer of mobile sediments which covers the rigid bottom of the flume. At the end of all the experiments, which lasted for the same amount of time, the formation of transverse sand dunes was observed. For decreasing , the rigid bottom of the flume was bared progressively earlier during the experiment and the measurements showed a clear tendency of the bedforms to lengthen, i.e. to increase their crest-to-crest distance. Moreover, under strong supply limitation, the two-dimensional transverse dunes turned into three-dimensional barchanoid forms and into isolated barchan dunes characterized by an abrupt reduction in bedform heights. A two-dimensional Fourier analysis of the bottom profile was performed, providing the amplitude of the main streamwise and spanwise harmonic components of the bottom morphology as a function of . © 2019 John Wiley & Sons, Ltd.     

Orthotropic anisotropy: on contributions of elasticity parameters to a difference in quasi-P-wave-squared velocities resulted from propagation in two orthogonal symmetry planes

We investigate the dependence of quasi P-wave phase velocity propagating in orthotropic media on particular elasticity parameters. Specifically, due to mathematical facilitation, we consider the squared-velocity difference, , resulted from propagation in two mutually perpendicular symmetry planes. In the context of the effective medium theory, may be viewed as a parameter evaluating the influence of cracks – embedded in the background medium – parallel to one or both aforementioned planes. Our investigation is both theoretical and numerical. Based on Christoffel's equations, we propose two accurate approximations of . Due to them, we interpret the aforementioned squared-velocity difference as being twice more dependent on , than on . To describe the magnitude of the dependence, we consider the proportions between the partial derivatives of . Further, it occurs that is influenced by the ratio of vertically propagating quasi P-wave to vertically propagating quasi S-wave. Anomalously high might be caused by the low P/S ratio, which in turn can be an indicator of the presence of gas in natural fractures or aligned porosity. Also, we carry out numerical sensitivity study, according to which is approximately twice more dependent on than on , twice more sensitive to than to , and equally dependent on as on . The dependence on and can be neglected, especially for small phase angles. We verify the approximations and perform the sensitivity study, using eight examples of the elasticity tensors.     

Post-fire management effects on sediment (dis)connectivity in Mediterranean forest ecosystems: Channel and catchment response

Forest fires and post-fire practices influence sediment connectivity (SC). In this study, we use the ‘aggregated index of connectivity’ (AIC) to assess SC in five Mediterranean catchments (198–1090 ha) affected by a wildfire in 2012 in south-eastern Spain. Two temporal scenarios were considered, immediately after the fire and before post-fire management, and 2 years after the fire including all practices (hillslope barriers, check-dams, afforestation, salvage logging and skid trails). One LiDAR (light detection and ranging)-derived digital elevation model (DEM, 2 m × 2 m resolution) was generated, per scenario. The five catchment outlets were established as the computation target (AIC_OUT), and structural and functional SC were calculated. Index outputs were normalized to make the results of the non-nested catchments comparable (AIC_N-OUT). The output analysis includes the SC distribution along the catchments and at local scale (929 sub-catchments, 677 in the burned area), the hillslope and channel measures' effect on SC, and a sedimentological analysis using observed area-specific sediment yield (SSY) at 10 new (built after post-fire practices) concrete check-dams located in the catchments (SSY = 1.94 Mg ha⁻¹ yr⁻¹; σ = 1.22). The catchments with more circular shapes and steeper slopes were those with higher AIC_N-OUT. The structural SC maps – removing the rainfall erosivity influence – allowed evaluating the actual role played by the post-fire practices that reduced SC ( $\overline{x}$= − 1.19%; σ = 0.41); while functional SC was linked to the actual change of SC ( $\overline{x}$= + 5.32%; σ = 0.62). Hillslope treatments resulted in significant changes on AIC_N-OUT at sub-catchment scale with certain disconnectivity. A good and positive correlation was found between the SSY and the changes of AIC_N-OUT. However, the coarse DEM resolution explained the lack of effect of the rock check-dams – located on the secondary channels – on AIC_N-OUT. AIC_N-OUT proved to be a useful tool for decision making in post-fire restoration, but an optimal input data is still necessary to refine calculations.     

Multi-parameter reflection waveform inversion for acoustic transversely isotropic media with a vertical symmetry axis

Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures.