Estimating River Conductance from Prior Information to Improve Surface‐Subsurface Model Calibration

Most groundwater models simulate stream‐aquifer interactions with a head‐dependent flux boundary condition based on a river conductance (CRIV). CRIV is usually calibrated with other parameters by history matching. However, the inverse problem of groundwater models is often ill‐posed and individual model parameters are likely to be poorly constrained. Ill‐posedness can be addressed by Tikhonov regularization with prior knowledge on parameter values. The difficulty with a lumped parameter like CRIV, which cannot be measured in the field, is to find suitable initial and regularization values. Several formulations have been proposed for the estimation of CRIV from physical parameters. However, these methods are either too simple to provide a reliable estimate of CRIV, or too complex to be easily implemented by groundwater modelers. This paper addresses the issue with a flexible and operational tool based on a 2D numerical model in a local vertical cross section, where the river conductance is computed from selected geometric and hydrodynamic parameters. Contrary to other approaches, the grid size of the regional model and the anisotropy of the aquifer hydraulic conductivity are also taken into account. A global sensitivity analysis indicates the strong sensitivity of CRIV to these parameters. This enhancement for the prior estimation of CRIV is a step forward for the calibration and uncertainty analysis of surface‐subsurface models. It is especially useful for modeling objectives that require CRIV to be well known such as conjunctive surface water‐groundwater use.     

Simulating Stable Isotope Ratios in Plumes of Groundwater Pollutants with BIOSCREEN‐AT‐ISO

BIOSCREEN is a well‐known simple tool for evaluating the transport of dissolved contaminants in groundwater, ideal for rapid screening and teaching. This work extends the BIOSCREEN model for the calculation of stable isotope ratios in contaminants. A three‐dimensional exact solution of the reactive transport from a patch source, accounting for fractionation by first‐order decay and/or sorption, is used. The results match those from a previously published isotope model but are much simpler to obtain. Two different isotopes may be computed, and dual isotope plots can be viewed. The dual isotope assessment is a rapidly emerging new approach for identifying process mechanisms in aquifers. Furthermore, deviations of isotope ratios at specific reactive positions with respect to “bulk” ratios in the whole compound can be simulated. This model is named BIOSCREEN‐AT‐ISO and will be downloadable from the journal homepage.     

Simulating Head Recovery and Intermittent Discharge in Free Flowing Wells

Modeling of the head recovery (buildup) in and intermittent operation of flowing wells lacks rigorous analytical tools. Presently different methods are applied in discharging and recovery phases. The presented reliable approximate solution considers a continuous free flow process and time variant screen loss parameters. The latter are temporarily set to large value to prevent inflow into the well bore during the single recovery or in repeated shut down periods. This method is called as screen loss control or SLC technique and is verified by means of comparative analysis with the constant rate pumping‐recovery test simulated in leaky aquifer. The evaluation of the free flow and recovery test in the free flowing well 28 (Artesia Heights) concluded with calibrated parameters close to those obtained in the earlier analyses based on split well flow process. The simulation of intermittent flowing well operation also underwent successful validation resulting in reasonable values of the Q_w and s response functions.     

Occurrence and distribution of bacteria indicators,chemical tracers and pathogenic vibrios in Singapore coastal waters

Water quality in Singapore's coastal area was evaluated with microbial indicators, pathogenic vibrios, chemical tracers and physico-chemical parameters. Sampling sites were grouped into two clusters (coastal sites at (i) northern and (ii) southern part of Singapore). The coastal sites located at northern part of Singapore along the Johor Straits exhibited greater pollution. Principal component analysis revealed that sampling sites at Johor Straits have greater loading on carbamazepine, while turbidity poses greater influence on sampling sites at Singapore Straits. Detection of pathogenic vibrios was also more prominent at Johor Straits than the Singapore Straits. This study examined the spatial variations in Singapore's coastal water quality and provided the baseline information for health risk assessment and future pollution management.     

Paleomagnetism and magnetic fabrics of Mio-Pliocene hypabyssal rocks of the Combia event,Colombia: tectonic implications

Mio-Pliocene hypabyssal rocks of the Combia event in the Amagá basin (NW Andes-Colombia), contain a deformational record of the activity of the Cauca-Romeral fault system, and the interaction of terranes within the Choco and northern Andean blocks. Previous paleomagnetic studies interpreted coherent counterclockwise rotations and noncoherent modes of rotation about horizontal axes for the Combia intrusives. However, rotations were determined from in-situ paleomagnetic directions and the existing data set is small. In order to better understand the deformational features of these rocks, we collected new paleomagnetic, structural, petrographic and magnetic fabric data from well exposed hypabyssal rocks of the Combia event. The magnetizations of these rocks are controlled by a low-coercivity ferromagnetic phase. Samples respond well to alternatingfield demagnetization isolating a magnetization component of moderate coercivity. These rocks do not have ductile deformation features. Anisotropy of magnetic susceptibility and morphotectonic analysis indicate that rotation about horizontal axes is consistently to the south-east, suggesting the need to apply a structural correction to the paleomagnetic data. The relationships between magnetic foliations and host-rock bedding planes indicate tectonic activity initiated before ~10 Ma. We present a mean paleomagnetic direction (declination D = 342.8°, inclination I = 12.1°, 95% confidence interval α₉₅ = 12.5°, precision parameter k = 8.6, number of specimens n = 18) that incorporates structural corrections. The dispersion S = 27° of site means cannot be explained by secular variation alone, but it indicates a counterclockwise rotation of 14.8° ± 12.7° relative to stable South America. Paleomagnetic data within a block bounded by the Sabanalarga and Cascajosa faults forms a more coherent data set (D = 336.5°, I = 17.4°, α₉₅ = 11.7°, k = 12.5, n = 14), which differs from sites west of the Sabanalarga fault and shows a rotation about a vertical axis of 20.2° ± 10.7°. Deformation in the Amagá basin may be tentatively explained by the obduction of the Cañas Gordas terrane over the northwestern margin of the northern Andean block. However, it can also be related to the local effects of the Cauca-Romeral fault system.     

Geodynamic Evolution of Northeastern Tunisia During the Maastrichtian–Paleocene Time: Insights from Integrated Seismic Stratigraphic Analysis

The Maastrichtian–Paleocene El Haria formation was studied and defined in Tunisia on the basis of outcrops and borehole data; few studies were interested in its three-dimensional extent. In this paper, the El Haria formation is reviewed in the context of a tectono-stratigraphic interval using an integrated seismic stratigraphic analysis based on borehole lithology logs, electrical well logging, well shots, vertical seismic profiles and post-stack surface data. Seismic analysis benefits from appropriate calibration with borehole data, conventional interpretation, velocity mapping, seismic attributes and post-stack model-based inversion. The applied methodology proved to be powerful for charactering the marly Maastrichtian–Paleocene interval of the El Haria formation. Migrated seismic sections together with borehole measurements are used to detail the three-dimensional changes in thickness, facies and depositional environment in the Cap Bon and Gulf of Hammamet regions during the Maastrichtian–Paleocene time. Furthermore, dating based on their microfossil content divulges local and multiple internal hiatuses within the El Haria formation which are related to the geodynamic evolution of the depositional floor since the Campanian stage. Interpreted seismic sections display concordance, unconformities, pinchouts, sedimentary gaps, incised valleys and syn-sedimentary normal faulting. Based on the seismic reflection geometry and terminations, seven sequences are delineated. These sequences are related to base-level changes as the combination of depositional floor paleo-topography, tectonic forces, subsidence and the developed accommodation space. These factors controlled the occurrence of the various parts of the Maastrichtian–Paleocene interval. Detailed examinations of these deposits together with the analysis of the structural deformation at different time periods allowed us to obtain a better understanding of the sediment architecture in depth and the delineation of the geodynamic evolution of the region.     

A wind tunnel investigation of particle segregation,ripple formation and armouring within sand beds of systematically varied texture

This paper reports on a wind tunnel investigation of particle segregation, ripple formation and surface armouring within sand beds of systematically varied particle size distribution, from coarsely skewed to bimodal. By design, the system was closed with no external inputs of mass from an external particle feed. Particles too coarse to travel in saltation for the given range in wind speed were dyed red in order to distinguish them in optical images from finer sand particles, which could be entrained into the unidirectional airflow. A 3D laser scanner measured the changing bed topography at regular time intervals during 18 experiments involving varied combinations of wind speed and bed texture. Image classification techniques were used to investigate the coincident self‐organization of the two populations of particles, as distinguished by their colour. As soon as saltation commenced, some of the red particles segregated into thin discontinuous patches. Particle trapping and sheltering on these rough patches was strongly favoured, causing them to grow preferentially. During the earliest stages of formation, bedform growth coincided with: (i) rapid coarsening of the surface texture; and (ii) the merging of proto‐ripple ‘crests’ to generate larger rhythmic bedforms of lower frequency. Consistent with previous work, ripple size was observed to increase under stronger winds when not exceeding the threshold for entrainment of the coarse‐mode or red particles from the crest. With declining rates of mass transport and particle segregation as the bed surface armoured, and the consequent deceleration of ripple propagation through to the end of each experiment, all surfaces eventually attained a steady‐state morphometry. At saturation, the largest ripples developed on beds having the lowest initial concentration of red particles. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessment of soil factors controlling ephemeral gully erosion on agricultural fields

The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion. Copyright © 2018 John Wiley & Sons, Ltd.     

Transport of moderately sorted gravel at low bed shear stresses: The role of fine sediment infiltration

A reliable estimation of sediment transport in gravel‐bed streams is important for various practical engineering and biological studies (e.g., channel stability design, bed degradation/aggradation, restoration of spawning habitat). In the present work, we report original laboratory experiments investigating the transport of gravel particles at low bed shear stresses. The laboratory tests were conducted under unsteady flow conditions inducing low bed shear stresses, with detailed monitoring of the bed topography using a laser scanner. Effects of bed surface arrangements were documented by testing loose and packed bed configurations. Effects of fine sediments were examined by testing beds with sand, artificial fine sand or cohesive silt infiltrated in the gravel matrix. Analysis of the experimental data revealed that the transport of gravel particles depends upon the bed arrangement, the bed material properties (e.g., size and shape, consolidation index, permeability) and the concentration of fine sediments within the surface layer of moving grains. This concentration is directly related to the distribution of fine particles within the gravel matrix (i.e., bottom‐up infiltration or bridging) and their transport mode (i.e., bedload or suspended load). Compared to loose beds, the mobility of gravel is reduced for packed beds and for beds clogged from the bottom up with cohesive fine sediments; in both cases, the bed shear stress for gravel entrainment increases by about 12%. On the other hand, the mobility of gravel increases significantly (bed shear stress for particle motion decreasing up to 40%) for beds clogged at the surface by non‐cohesive sand particles. Copyright © 2017 John Wiley & Sons, Ltd.