Inversion of interference hydraulic pumping tests in both homogeneous and fractal dual media

This work proposes a complete method for automatic inversion of data from hydraulic interference pumping tests based on both homogeneous and fractal dual-medium approaches. The aim is to seek a new alternative concept able to interpret field data, identify macroscopic hydraulic parameters and therefore enhance the understanding of flow in porous fractured reservoirs. Because of its much contrasted sensitivities to parameters, the dual-medium approach yields an ill-posed inverse problem that requires a specific optimization procedure including the calculation of analytical sensitivities and their possible re-scaling. Once these constraints are fulfilled, the inversion proves accurate, provides unambiguous and reliable results. In the fractal context inverting several drawdown curves from different locations at the same time reveals more accurate. Finally, hydraulic parameters drawn from inversion should be taken into account to improve in various situations the conditioning of up-scaled flow in fractured rocks.     

Combining pump-and-treat and physical barriers for contaminant plume control

A detailed analysis is presented of the hydraulic efficiency of plume management alternatives that combine a conventional pump-and-treat system with vertical, physical hydraulic barriers such as slurry walls or sheet piles. Various design settings are examined for their potential to reduce the pumping rate needed to obtain a complete capture of a given contaminated area. Using established modeling techniques for flow and transport, those barrier configurations (specified by location, shape, and length) that yield a maximum reduction of the pumping rate are identified assuming homogeneous aquifer conditions. Selected configurations are further analyzed concerning their hydraulic performance under heterogeneous aquifer conditions by means of a stochastic approach (Monte Carlo simulations) with aquifer transmissivity as a random space function. The results show that physical barriers are an appropriate means to decrease expected (mean) pumping rates, as well as the variance of the corresponding pumping rate distribution at any given degree of heterogeneity. The methodology presented can be transferred easily to other aquifer scenarios, provided some basic premises are fulfilled, and may serve as a basis for reducing the pumping rate in existing pump-and-treat systems.     

Measuring the flow resistance of submerged grass

Through laboratory experiments conducted in a grass‐lined flume, the hydraulic resistance of grass is measured and quantified. For the grass examined, it is found that Manning's n value is greater than those recommended in well‐established texts such as Chow (1959. Open Channel Flow. McGraw‐Hill: Singapore), relatively lower than those predicted by n‐UR methods, but corresponds well with the value found from calibration studies of two‐ and three‐dimensional numerical models. The assumption of a uniform Manning's n value with flow depth, which is often made in numerical modelling, may be invalid depending on the relative submergence of the vegetation. Drag coefficients are evaluated for a method applicable to three‐dimensional numerical models. Further detailed experimental investigation and application of these approaches within a numerical modelling framework is now recommended. Copyright © 2002 John Wiley & Sons, Ltd.     

Experimental and Numerical Investigation of the Flow Hydraulic in Gradual Transition Open Channels

Channel expansions are common in both natural and artificial open channels. With increasing cross-sectional dimensions in an expansion, the flow decelerates. Due to separation of flow and subsequent eddy formation, a significant head loss is occurred along the transition. This study presents the results of experimental investigations on subcritical flow along the expansive transition of rectangular to trapezoidal channels. Also, a numerical simulation was developed using the finite volume method with Reynolds Stress turbulent model. Water surface profiles and velocity distributions of flow through the transition were measured experimentally and compared with the numerical results. Also, hydraulic efficiency of the transition and coefficient of energy head loss were calculated. The results show that with increasing the upstream Froude number, hydraulic efficiency of the transition and coefficient of energy head loss are decreased and increased, respectively. The results also showed the ability of numerical simulation for simulating the flow separation zones and bed shear stress along the transition for different inlet discharges and inflow Froude numbers.     

On estimating the hydraulic properties of soil,Part 3. Parameters of the Philip infiltration equation

Using a large number of experimental data, the paper estimates the parameters of the Philip infiltration equation. The values of these parameters are obtained using the Philip approach, which employs the hydraulic properties of soil. In the absence of extensive measurements, the tabulated results can be used in modelling infiltration. The paper introduces approximate sorptivity expressions that use a limited number of parameters. Sorptivity may be estimated reasonably well from a knowledge of moisture deficit, a representative capillary pressure, and saturated hydraulic conductivity. The accuracy of prediction is highest for coarser materials. For coarse materials (sands), second-power sorptivity can be approximated as a linear function of initial effective saturation. The proposed expressions for sorptivity can be used as a substitute for the lengthy Philip procedure.     

Variability of bedload transport and channel morphology in a braided river hydraulic model

This paper investigates variability in bedload transport and channel morphology for 11 replicate experimental runs in an approximately 1:50 braided river model. The experiments, each of 90 h duration, were carried out in a 20 × 3m tilting flume. All the experiments started with the same initial conditions. Bedload transport was measured at 5 min intervals in a collection drum at the exit from the flume. The model showed reasonable hydraulic similarity when compared to prototype rivers. Results show that mean bedload transport rates for the 11 runs vary in the range 0·98 to 1·49gs^?1 (mean + 1·21, coefficient of variation 11 per cent). Within-run transport rates commonly vary from close to zero, to two and occasionally three or four times the mean rate. Within the bedload series, several irregular phases of transport intensity can be observed, but time series analysis of the data show little underlying serial structure (an AR(2) autoregressive model is appropriate). Channel patterns are narrow/braided, are established quickly and remain relatively stable throughout the runs, although channel widths increase between 20 and 103 per cent over the 11 runs. Channel behaviour varies from aggradational to transitional between aggradation and degradation. Time-averaged bedload transport rate is weakly correlated with braiding intensity. In general, these results demonstrate that for a given set of controlling variables, bedload transport and channel morphology can be approximately replicated.     

An Analytical Model for Hydraulic Fracturing in Shallow Bedrock Formations

A theoretical method is proposed to estimate post‐fracturing fracture size and transmissivity, and as a test of the methodology, data collected from two wells were used for verification. This method can be employed before hydrofracturing in order to obtain estimates of the potential hydraulic benefits of hydraulic fracturing. Five different pumping test analysis methods were used to evaluate the well hydraulic data. The most effective methods were the Papadopulos‐Cooper model (1967), which includes wellbore storage effects, and the Gringarten‐Ramey model (1974), known as the single horizontal fracture model. The hydraulic parameters resulting from fitting these models to the field data revealed that as a result of hydraulic fracturing, the transmissivity increased more than 46 times in one well and increased 285 times in the other well. The model developed by dos Santos (2008) , which considers horizontal radial fracture propagation from the hydraulically fractured well, was used to estimate potential fracture geometry after hydrofracturing. For the two studied wells, their fractures could have propagated to distances of almost 175 m or more and developed maximum apertures of about 2.20 mm and hydraulic apertures close to 0.30 mm. Fracturing at this site appears to have expanded and propagated existing fractures and not created new fractures. Hydraulic apertures calculated from pumping test analyses closely matched the results obtained from the hydraulic fracturing model. As a result of this model, post‐fracturing geometry and resulting post‐fracturing well yield can be estimated before the actual hydrofracturing.     

Model studies of dense water overflows in the Faroese Channels

The overflow of dense water from the Nordic Seas through the Faroese Channel system was investigated through combined laboratory experiments and numerical simulations using the Massachusetts Institute of Technology General Circulation Model. In the experimental study, a scaled, topographic representation of the Faroe-Shetland Channel, Wyville-Thomson Basin and Ridge and Faroe Bank Channel seabed bathymetry was constructed and mounted in a rotating tank. A series of parametric experiments was conducted using dye-tracing and drogue-tracking techniques to investigate deep-water overflow pathways and circulation patterns within the modelled region. In addition, the structure of the outflowing dense bottom water was investigated through density profiling along three cross-channel transects located in the Wyville-Thomson Basin and the converging, up-sloping approach to the Faroe Bank Channel. Results from the dye-tracing studies demonstrate a range of parametric conditions under which dense water overflow across the Wyville-Thomson Ridge is shown to occur, as defined by the Burger number, a non-dimensional length ratio and a dimensionless dense water volume flux parameter specified at the Faroe-Shetland Channel inlet boundary. Drogue-tracking measurements reveal the complex nature of flow paths and circulations generated in the modelled topography, particularly the development of a large anti-cyclonic gyre in the Wyville-Thompson Basin and up-sloping approach to the Faroe Bank Channel, which diverts the dense water outflow from the Faroese shelf towards the Wyville-Thomson Ridge, potentially promoting dense water spillage across the ridge itself. The presence of this circulation is also indicated by associated undulations in density isopycnals across the Wyville-Thomson Basin. Numerical simulations of parametric test cases for the main outflow pathways and density structure in a similarly-scaled Faroese Channels model domain indicate excellent qualitative agreement with the experimental observations and measurements. In addition, the comparisons show that strong temporal variability in the predicted outflow pathways and circulations have a strong influence in regulating the Faroe Bank Channel and Wyville-Thomson Ridge overflows, as well as in determining the overall response in the Faroese Channels to changes in the Faroe-Shetland Channel inlet boundary conditions.     

Spatial variability of the poroelastic response of the heterogeneous medium

In this work, a stochastic methodology is applied to analyze the variability of the poroelastic response of the heterogeneous medium at the field scale. To solve the problem analytically, we restrict our attention to the one-dimensional models, where fluid flow as well as deformation occurs in one direction only under a constant applied stress. Assuming statistic homogeneity, the closed-form solutions that describe the variability of fluid pressure head, and a solid's strain and displacement are developed using a spectral approach based on Fourier–Stieltjes representations for the perturbed quantities. The influence of the correlation length of the log hydraulic conductivity on these results is investigated. It is found that the variances of the solid's strain and displacement increase with the correlation length of the log hydraulic conductivity, while the correlation length of the log hydraulic conductivity plays the role in reducing the variability of the specific discharge.     

LiDAR-based fluvial remote sensing to assess 50–100-year human-driven channel changes at a regional level: The case of the Piedmont Region,Italy

Remote Sensing (RS) technology has recently offered new and promising opportunities to analyze river systems. In this paper, we present a calibration of characteristic Hydraulic Scaling Law (HSL) using a regional database of river geomorphic features. We consistently linked discharge with channel geometry features for estimated Bankfull Channel Depth (eBCD), Active Channel Width (ACW), and Low Flow water Channel Width (LFCW), which are continuously available from RS data along the river course. We then used historical information and external sources of information on channel reaches that were relatively unaffected by human pressure over periods ranging from a few decades to a century (measured in comparable geographical areas) to infer relatively Unaltered HSLs (rUHSLs). Adopting rUHSL validated with available local historical evidence on channel geometry, we were able to assess historical changes in channel geometry consistently over the entire region and within the studied temporal window. The case study was conducted for the Po basin in the Piedmont Region, north-west Italy. From our analysis, it emerges that regionally 74% of the river network has riverbed incisions exceeding 1 m, while 66% of channels have halved their historical widths with a total of 617 ha of land subtracted from the active channel. LFCW is, on average, wider in Alpine rivers compared with those located in the North Apennines. Although it is currently not possible to measure the accuracy of these estimates, the evidence generated is coherent with available historical information, characteristic hydraulic scaling laws, evidence from relatively unaltered reaches and the available literature on local fluvial systems. This methodology provides robust, novel and quantitative information regarding decadal to secular channel changes that have occurred on a regional scale. This new layer of information enriches our ability to rationally address assessments of large-scale past and future channel trajectories. © 2018 John Wiley & Sons, Ltd.     

Dissolved and particulate trace metal fluxes through the central English Channel, and the influence of coastal gyres

Measurements of dissolved Cd, Co, Cu, Mn, Ni, Pb, and Zn have been made on a seasonal basis at five stations on a north–south transect across the central English Channel between Cherbourg and the Isle of Wight. Vertical and horizontal distributions of dissolved Cd, Pb, Cu and Zn are relatively uniform except for sampling sites near the English coast. Dissolved Mn and Co show increased concentrations in the English coastal waters, and for Mn the seasonal trend in concentration follows the pattern seen in the Strait of Dover with higher values in the late summer. Ni and Cu are higher in concentration on the English side, which reflects mainly riverine sources. Measurements were also made of particulate forms of the metals above plus particulate Al, Ca, Fe, Mg, Sr and Ti. Water column concentrations of particulate metals broadly follow the distribution of suspended particulate matter, with highest concentrations near the UK coast. Trace metal concentrations have been integrated with modelled data on fluxes of water to provide estimates of fluxes for these elements into the eastern Channel, and an initial comparison is made with data for fluxes of metals through the Strait of Dover obtained during an earlier study. A major influence on the fluxes of particulate metals through the Isle of Wight-Cherbourg transect is the gyre system to the South east to the Isle of Wight, which has important east to west as well as west to east transport components. For those elements where the dissolved form of the metal dominates, the large flow of water in the central Channel waters leads to major fluxes of the metals towards the east and the Strait of Dover. However, the high suspended particulate matter loadings in the coastal waters and impact of the gyre system lead to net east to west fluxes of particulate Al, Fe, Mn and Ti. Comparison of these fluxes with data on the net west to east transport of these materials through the Strait of Dover infers that there must be a significant supply of these particulate metals to the eastern Channel.     

Quantitative assessment of dissolved radiotracers in the English Channel: Sources, average impact of la Hague reprocessing plant and conservative behaviour (1983, 1986, 1988, 1994)

Knowledge of long-term movements of water-masses in the English Channel has been substantially improved using hydrodynamic modelling coupled with radio-tracers studies; nevertheless, the precision of results so obtained is still largely dependent on measurement precision. New tools are now available to make more accurate determinations of radio-tracer distribution: (1) Repositioning of station locations at the same tide reference-time, giving a homogeneous spatial data set, coupled with the possibility of interpolating and quantifying the amounts of dissolved radioactivity flowing through the English Channel; (2) the first measurements of tritium (³H) in seawater on a large scale in the English Channel demonstrate that this fully conservative radionuclide is a clearly identifiable marker of industrial releases; (3) recent campaigns carried out during the FLUXMANCHE II CCE (1994) programme show the general distribution of dissolved radionuclides ¹³⁷Cs, ¹³⁴Cs, ⁶⁰Co, ¹²⁵Sb, ¹⁰⁶Ru and ³H in the English Channel and the Irish Sea; and (4) the re-utilisation of data from previous campaigns (1983, 1986, 1988) provides indications, at any given location in the English Channel, about the average dilution and distribution of releases derived from the La Hague reprocessing plant. Excesses and losses of radionuclides are now quantified with respect to known source terms; estimates of losses are provided for non-conservative radionuclides, while an excess of ¹³⁷Cs was observed in the English Channel during the period 1983–1994. This excess, which has the same order of magnitude as the quantities released from La Hague plant in the English Channel, could be explained by about 1% of the Sellafield reprocessing plant releases entering the Channel. These results confirm and give a more detailed picture of the previously known distribution of water masses in the English Channel. They lead to clear information about transit times and dilution at this scale, and provide directly comparable data for the validation of hydrodynamic models.     

Projected entrainment of fish resulting from aggregate dredging

Previous research to assess impacts from aggregate dredging has focussed on infaunal species with few studies made of fish entrainment. Entrainment evidence from hydraulic dredging studies is reviewed to develop a sensitivity index for benthic fish. Environmental monitoring attendant with the granting of new licences in the Eastern Channel Region (ECR) in 2006 offers a unique opportunity to assess the effects of dredging upon fish. Projected theoretical fish entrainment rates are calculated based upon: abundance data from 4m beam trawl sampling of fish species over the period 2005-2008; sensitivity data; and dredging activity and footprint derived from Electronic monitoring System (EMS) data. Results have been compared with actual entrainment rates and also against summary results from independent analysis of the changes in fish population over the period 2005-2008 (Drabble, 2012). The case is made for entrainment surveys to form part of impact monitoring for marine aggregate dredging.     

Conceptualization and Calibration of Anisotropic Alluvial Systems: Pitfalls and Biases

Physical properties of alluvial environments typically feature a high degree of anisotropy and are characterized by dynamic interactions between the surface and the subsurface. Hydrogeological models are often calibrated under the assumptions of isotropic hydraulic conductivity fields and steady-state conditions. We aim at understanding how these simplifications affect predictions of the water table using physically based models and advanced calibration and uncertainty analysis approaches based on singular value decomposition and Bayesian analysis. Specifically, we present an analysis of the information content provided by steady-state hydraulic data compared to transient data with respect to the estimation of aquifer and riverbed hydraulic properties. We show that assuming isotropy or fixed anisotropy may generate biases both in the estimation of aquifer and riverbed parameters as well as in the predictive uncertainty of the water table. We further demonstrate that the information content provided by steady-state hydraulic heads is insufficient to jointly estimate the aquifer anisotropy together with the aquifer and riverbed hydraulic conductivities and that transient data can help to reduce the predictive uncertainty to a greater extent. The outcomes of the synthetic analysis are applied to the calibration of a dynamic and anisotropic alluvial aquifer in Switzerland (The Rhône River). The results of the synthetic and real world modeling and calibration exercises documented herein provide insight on future data acquisition as well as modeling and calibration strategies for these environments. They also provide an incentive for evaluation and estimation of commonly made simplifying assumptions in order to prevent underestimation of the predictive uncertainty.     

Hydraulic Performance of the Horizontal Reactive Media Treatment Well: Pilot and Numerical Study

This study is focused on a passive treatment system known as the horizontal reactive treatment well (HRX Well®) installed parallel to groundwater flow, which operates on the principle of flow focusing that results from the hydraulic conductivity (K) ratio of the well and aquifer media. Passive flow and capture in the HRX Well are described by simplified equations adapted from Darcy's Law. A field pilot-scale study (PSS) and numerical simulations using a finite element method (FEM) were conducted to verify the HRX Well concept and test the validity of the HRX Well-simplified equations. The hydraulic performance results from both studies were observed to be within a close agreement to the simplified equations and their hydraulic capture width approximately five times greater than the well diameter (0.20 m). Key parameters affecting capture included the aquifer thickness, well diameter, and permeability ratio of the HRX Well treatment media and aquifer material. During pilot testing, the HRX Well captured 39% of flow while representing 0.5% of the test pit cross-sectional volume, indicating that the well captures a substantial amount of surrounding groundwater. While uncertainty in the aquifer and well properties (porosity, K, well losses), including the effects of boundary conditions, may have caused minor differences in the results, data from this study indicate that the simplified equations are valid for the conceptual design of a field study. A full-scale HRX Well was installed at Site SS003 at Vanderberg Air Force Base, California, in July/August 2018 based on outcomes from this study.     

Modelling larval dispersal of the king scallop (Pecten maximus) in the English Channel: examples from the bay of Saint-Brieuc and the bay of Seine

The king scallop (Pecten maximus) is one of the most important benthic species of the English Channel as it constitutes the first fishery in terms of landings in this area. To support strategies of spatial fishery management, we develop a high-resolution biophysical model to study scallop dispersal in two bays along the French coasts of the English Channel (i.e. the bay of Saint-Brieuc and the bay of Seine) and to quantify the relative roles of local hydrodynamic processes, temperature-dependent planktonic larval duration (PLD) and active swimming behaviour (SB). The two bays are chosen for three reasons: (1) the distribution of the scallop stocks in these areas is well known from annual scallop stock surveys, (2) these two bays harbour important fisheries and (3) scallops in these two areas present some differences in terms of reproductive cycle and spawning duration. The English Channel currents and temperature are simulated for 10 years (2000–2010) with the MARS-3D code and then used by the Lagrangian module of MARS-3D to model the transport. Results were analysed in terms of larval distribution at settlement and connectivity rates. While larval transport in the two bays depended both on the tidal residual circulation and the wind-induced currents, the relative role of these two hydrodynamic processes varied among bays. In the bay of Saint-Brieuc, the main patterns of larval dispersal were due to tides, the wind being only a source of variability in the extent of larval patch and the local retention rate. Conversely, in the bay of Seine, wind-induced currents altered both the direction and the extent of larval transport. The main effect of a variable PLD in relation to the thermal history of each larva was to reduce the spread of dispersal and consequently increase the local retention by about 10 % on average. Although swimming behaviour could influence larval dispersal during the first days of the PLD when larvae are mainly located in surface waters, it has a minor role on larval distribution at settlement and retention rates. The analysis of the connectivity between subpopulations within each bay allows identifying the main sources of larvae which depend on both the characteristics of local hydrodynamics and the spatial heterogeneity in the reproductive outputs.     

Hydraulic modeling of flood pulses in the Middle Bussento Karst System (MBSKS), UNESCO Cilento Global Geopark,southern Italy

Karst systems provide water for domestic and industrial uses and for generating hydropower, but they can also create fluvial hazards, such as upstream back‐flooding and downstream karst flash‐flood events. However, these hazards are difficult to foresee due to the complex recharge‐discharge processes as well as the lack of information on the inside of the system, which has often not been completely surveyed by speleologists or explored by boreholes. To overcome these difficulties, hydro‐chemical data from the monitoring system in the Middle Bussento Karst System (MBSKS), one of the first Experimental Karst Systems in southern Italy, were recorded and previously discussed. Based on shared background in flood karst hydraulic modeling, this paper describes the conceptual premises and rationale of a general‐purpose hydraulic model that is suitable both for the MBSKS and for other Mediterranean, multi‐recharge, mature, conduit‐dominated karst systems. To test the reliability of the model, simulations of time–space behavior and response are performed using natural and artificial flood pulses “as tracers”, considering a “pulse” as a significant variation in water quantity and/or quality. The results of the model explain the interactions between allogenic, autogenic, and anthropogenic recharges from differentiated sources and phreatic conduit systems. These results also clarify the overall response of karst springs at typical time scales of flood pulses. Table acronym name     

Reducing Groundwater Pumping Cost Using Alternate Pulsed Pumping

This paper deals with water pumping cost minimization, in a confined infinite aquifer, proposing an alternate pulsed pumping schedule. The transient flow analysis is conducted for two wells with equal pumping rates. Specifically, two pumping schedules are analytically compared. In the first case, well users pump simultaneously, and in the second one they cooperate so that they pump alternately. This paper proves that the proposed alternate pumping schedule works as a stabilizer, reducing the high hydraulic drawdowns values, regardless of the aquifer characteristics. Moreover, pumping alternately is better in terms of pumping cost, compared to simultaneous pumping, though benefit become negligible as distance between wells becomes large. Two simplified equations are proposed, one to find the difference of the hydraulic drawdowns between the two pumping schedules and the other one to find the economic benefit of each well from cooperation. The equations are finally applied to a number of cases and their results are compared to the results obtained from an algorithm created to calculate the hydraulic drawdowns and the pumping cost, using the Theis equation. The results can be very useful in irrigation scheduling, as they can be applied to systems of well users/farmers, to reduce pumping cost.     

Distribution of trace metals and Pb isotopes in bottom sediments of the Murucupi River,North Brazil

The Murucupi River belongs to the hydrographic network of the Pará River estuary, at the southern portion of the Amazon River mouth, which consists of a fluvial-marine transitional zone under strong impact of both tidal and fluvial currents. The geochemical results obtained for bottom sediments from the Murucupi River, the Arrozal Channel, and the Pará River indicate a natural variation of Pb, Cr, Cu, Zn, and Ni content among these water ways with no significant anthropogenic influence. According to the threshold effects level (TEL), the contents of trace metals do not offer risk to the local biota. By contrast, the differences in the Pb isotopic composition of sediments in the Murucupi River, the Arrozal Channel, and the Pará River are significant. These isotopic signatures indicate an anthropogenic contribution of Pb in the Murucupi River originating from the domestic effluents of urban centers; industrial waste represented by red mud is not included. These results demonstrate that the Pb isotopic signature is a prospective indicator for future contamination of bottom sediments by trace metals and is useful for identifying contaminants among the possible anthropogenic sources.     

Dissolved nutrient distributions in the Central English Channel

Samples for analysis of nitrate plus nitrite, dissolved reactive phosphate and silicate were collected in the English Channel along a transect between Cherbourg and the Isle of Wight, within the framework of the FLUXMANCHE II programme. A total of 130 samples were obtained during 15 cruises between December 1994 and December 1995, thus giving a full seasonal coverage. Nutrient concentrations were generally highest on the UK side of the Channel, where riverine inputs were greatest. Nutrient fluxes across this transect from west to east for the winter period, when the nutrients exhibited quasi-conservative behaviour, were calculated using the N, P and Si data, and estimates of water fluxes derived from a two-dimensional hydrodynamic numerical model. A gyre system off the Isle of Wight led to localised east to west fluxes of nutrients. Although nutrient concentrations are lower in the central Channel than in the coastal zones, the nutrient fluxes are greatest through this central zone because of the higher fluxes of water through this part of the transect. The fluxes obtained through this mid-Channel transect are greater than fluxes estimated through the Straits of Dover for a similar winter period as part of the earlier FLUXMANCHE I (December 1990– March 1991) programme, suggesting significant variability from year to year, and/or removal of nutrients in the eastern Channel.