Numerical inverse interpretation of single-hole pneumatic tests in unsaturated fractured tuff

A numerical inverse method was used to interpret simultaneously multirate injection and recovery data from single-hole pneumatic tests in unsaturated fractured tuff at the Apache Leap Research Site near Superior, Arizona. Our model represents faithfully the three-dimensional geometry of boreholes at the site, and accounts directly for their storage and conductance properties by treating them as high-permeability and high-porosity cylinders of finite length and radius. It solves the airflow equations in their original nonlinear form and yields information about air permeability, air-filled porosity and dimensionless borehole storage coefficient. Some of this is difficult to accomplish with analytical type-curves. Air permeability values obtained by our inverse method agree well with those obtained by steady-state and type-curve analyses.     

Slug test analysis to evaluate permeability of compressible materials

The line-fitting methods such as the Hvorslev method and the Bouwer and Rice method provide a rapid and simple means to analyze slug test data for estimating in situ hydraulic conductivity ( k ) of geologic materials. However, when analyzing a slug test in a relatively compressible geologic formation, these conventional methods may have difficulties fitting a straight line to the semilogarithmic plot of the test data. Data from relatively compressible geologic formations frequently show a concave-upward curvature because of the effect of the compressibility or specific storage ( S _s). To take into account the compressibility of geologic formations, a modified line-fitting method is introduced, which expands on Chirlin's (1989) approach to the case of a partially penetrating well with the basic-time-lag fitting method. A case study for a compressible till is made to verify the proposed method by comparing the results from the proposed methods with those obtained using a type-curve method (Kansas Geological Survey method [ Hyder et al. 1994 ]).     

Applying Parameter-Estimation Methods to Recovery-Test and Slug-Test Analyses

Parameter-estimation methods, including an exhaustive-search method and PEST (Parameter ESTimation) software, were applied to recovery-test data and slug-test data to obtain best estimates of transmissivity (T) by minimizing the sums of residuals. Each residual represents the difference between the field-measured water-level value and the value calculated by the appropriate non-linear equation. The exhaustive-search method in both cases involves computing the sums of residuals for an array of transmissivity and storativity values selected by the user for testing. Two new Fortran programs are presented that employ the exhaustive-search method. They utilize Picking's method for analyzing recovery-test data and the analytical equation for analyzing slug-test data derived by Cooper, Bredehoeft, and Papadopulos. Picking's method involves application of the Papadopulos and Cooper's equation for drawdown in finite-diameter wells. Utilizing field data reported in the literature, the estimated transmissivity values from the exhaustive-search methods were compared to the literature values obtained by type-curve matching techniques. The exhaustive-search values corresponded closely to the curve-matching values. Estimates for T were also obtained from recovery-test and slug-test data from two sites in southeastern Pennsylvania. For these sites, the PEST program was also applied to the data to evaluate the accuracy of the exhaustive-search methods. The results from the two methods were generally in good agreement. The two new Fortran programs are practical tools for the hydrogeologist, as they require less time compared to type-curve matching and the PEST method, and they yield accurate estimates of transmissivity.     

Field study of hydrogeologic characterization methods in a heterogeneous aquifer

Hydraulic conductivity (K) and specific storage (S(s)) are required parameters when designing transient groundwater flow models. The purpose of this study was to evaluate the ability of commonly used hydrogeologic characterization approaches to accurately delineate the distribution of hydraulic properties in a highly heterogeneous glaciofluvial deposit. The metric used to compare the various approaches was the prediction of drawdown responses from three separate pumping tests. The study was conducted at a field site, where a 15 m × 15 m area was instrumented with four 18-m deep Continuous Multichannel Tubing (CMT) wells. Each CMT well contained seven 17 cm × 1.9 cm monitoring ports equally spaced every 2 m down each CMT system. An 18-m deep pumping well with eight separate 1-m long screens spaced every 2 m was also placed in the center of the square pattern. In each of these boreholes, cores were collected and characterized using the Unified Soil Classification System, grain size analysis, and permeameter tests. To date, 471 K estimates have been obtained through permeameter analyses and 270 K estimates from empirical relationships. Geostatistical analysis of the small-scale K data yielded strongly heterogeneous K fields in three-dimensions. Additional K estimates were obtained through slug tests in 28 ports of the four CMT wells. Several pumping tests were conducted using the multiscreen and CMT wells to obtain larger scale estimates of both K and S(s). The various K and S(s) estimates were then quantitatively evaluated by simulating transient drawdown data from three pumping tests using a 3D forward numerical model constructed using HydroGeoSphere (Therrien et al. 2005). Results showed that, while drawdown predictions generally improved as more complexity was introduced into the model, the ability to make accurate drawdown predictions at all CMT ports was inconsistent.     

Modeling time dependence in a tidal boundary layer

Although the rate of acceleration and deceleration in a depth-limited tidal flow can be quite large, it has been found that quasi-steady logarithmic models or steady eddy diffusivitics can explain the turbulent field quite well (Gross andNowell, 1983;Long, 1981). A simple, time-dependent model was devised based on kinetic energy closure (Bradshaw et al., 1967) to simulate a depth-limited tidal boundary layer without the assumption of a time-independent eddy diffusivity. Computer solutions revealed that the parameterS = u_s/δΩ (u_s is a turbulent velocity scale, δ is depth, Ω is tidal frequency) governs the phasing of the velocity to the forcing pressure gradient. However, for values ofS reasonable for depth-limited estuarine tidal flow the ratio of the turbulent energy production-dissipation balance to rate of change of turbulent energy is always quite large, indicating that the turbulent adjustment rate is fast enough for the flow to be modeled with quasi-steady techniques.     

In-Well Hydraulics of the Electromagnetic Borehole Flowmeter: Further Studies

The electromagnetic borehole flowmeter (EBF) is finding increasing application as a method for measuring hydraulic conductivity (K) distributions. A recent paper details an experimental/theoretical study of the effect of in-well hydraulics on calculated K distributions based on EBF measurements (Dinwiddie et al. 1999). Results showed that minimizing head loss associated with flow through the meter was the key to producing accurate K values. Using the same experimental procedures, the previous study has been extended to develop data from a larger diameter (1 inch) EBF, and to determine if an EBF can be calibrated effectively without using an inflatable packer to force all flow through the meter annulus. Both experiments were aimed at producing low head loss conditions. Results show that overall calibration can be accomplished in the absence of a packer, which reduces head losses to nonmeasurable levels, and use of the 1-inch EBF with a packer reduces head losses by a factor of 16 when compared with the 0.5-inch EBF studied by Dinwiddie et al. (1999).     

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.     

Selection and ranking of ground motion models for seismic hazard analysis in the Pyrenees

The issue addressed in this paper is the objective selection of appropriate ground motion models for seismic hazard assessment in the Pyrenees. The method of Scherbaum et al. (2004a) is applied in order to rank eight published ground motion models relevant to intraplate or to low deformation rate contexts. This method is based on a transparent and data-driven process which quantifies the model fit and also measures how well the underlying model assumptions are met. The method is applied to 15 accelerometric records obtained in the Pyrenees for events of local magnitude between 4.8 and 5.1, corresponding to moment magnitudes ranging from 3.7 to 3.9. Only stations at rock sites are considered. A total of 720 spectral amplitudes are used to rank the selected ground motion models. Some control parameters of these models, such as magnitude and distance definitions, may vary from one model to the other. It is thus important to correct the selected models for their difference with respect to the magnitude and distance definitions used for the Pyrenean data. Our analysis shows that, with these corrections, some of the ground motion models successfully fit the data. These are the Lussou et al. (2001) and the Berge-Thierry et al. (2003) models. According to the selected ground motion models, a possible scenario of a magnitude 6 event is proposed; it predicts response spectra accelerations of 0.08–0.1 g at 1 Hz at a hypocentral distance of 10 km.     

Fully Automated Objective-Based Method for Master Recession Curve Separation

The fully automated objective-based method for master recession curve (MRC) separation was developed by using Microsoft Excel spreadsheet and Visual Basic for Applications (VBA) code. The core of the program code is used to construct an MRC by using the adapted matching strip method (Posavec et al. 2006). Criteria for separating the MRC into two or three segments are determined from the flow-duration curve and are represented as the probable range of percent of flow rate duration. Successive separations are performed automatically on two and three MRCs using sets of percent of flow rate duration from selected ranges and an optimal separation model scenario, having the highest average coefficient of determination R², is selected as the most appropriate one. The resulting separated master recession curves are presented graphically, whereas the statistics are presented numerically, all in separate sheets. Examples of field data obtained from two springs in Istria, Croatia, are used to illustrate its application. The freely available Excel spreadsheet and VBA program ensures the ease of use and applicability for larger data sets.     

Source mechanism of the Alaskan earthquake of 1964 from amplitudes of free oscillations and surface waves — comments

Recently Ben-Menahem et al. (1972) (hereafter abbreviated to B.R.I.) proposed, on the basis of long-period surface-wave and free-oscillation data, a steeply dipping fault model for the Alaskan earthquake of 1964. This result differs from that obtained by Kanamori (1970) (hereafter abbreviated to K.70) who also used long-period Love and Rayleigh waves, but preferred a low-angle fault model. I would like to make several comments on this difference.     

Multimodel Validity Assessment of Groundwater Flow Simulation Models Using Area Metric Approach

We demonstrate the application of the Area Metric developed by Ferson et al. (2008) for multimodel validity assessment. The Area Metric quantified the degree of models' replicative validity: the degree of agreement between the observed data and the corresponding simulated outputs represented as their empirical cumulative distribution functions (ECDFs). This approach was used to rank multiple representations of a case study groundwater flow model of a landfill by their Area Metric scores. A multimodel approach allows to account for uncertainties that may either be epistemic (from lack of knowledge) or aleatory (from variability inherent in the system). The Area Metric approach enables explicit incorporation of model uncertainties, epistemic as well as aleatory, into validation assessment. The proposed approach informs understanding of the collected data and that of the model domain. It avoids model overfitting to a particular system state, and in fact is a blind assessment of the models' validity: models are not adjusted, or updated, to achieve a better numerical fit. This approach assesses the degree of models' validity, in place of the typical binary model validation/invalidation process. Collectively, this increases confidence in the model's representativeness that, in turn, reduces risk to model users.     

Evaluation of Pt and Au pressure scales based on MgO absolute pressure scale

The equations of state (EOSs) of MgO produced by two independent scale-free methods, (1) the simultaneous elastic wave velocity and in situ synchrotron X-ray measurements (Kono et al., 2010; Li et al., 2006) and (2) the first-principles calculations (Wu et al., 2008), agree well with each other to at least 150 GPa and 2000 K. Furthermore, the EOS from first-principles calculations also agrees well with shock wave data, another pressure-scale-free data. These agreements strongly support that these EOSs provide reliable absolute pressure scales. Here we evaluate Au and Pt EOSs based on the EOS of Wu et al. (2008) using the simultaneously measured volume data of MgO, Au, and Pt from the literature. The primary pressure scales developed by Tange et al. (2009) and Yokoo et al. (2009) using only pressure-scale-free experimental data of MgO, Au, and Pt produce internal consistent pressure and agree with EOS of Wu et al. (2008). The Au EOS by Tsuchiya (2003) works well at room temperature but underestimates pressure at high temperature. The Au EOS by Fei et al. (2007) can well describe thermal pressure. The EOSs of Pt by Holmes et al. (1989) and Ono et al. (2011) work well at both room temperature and high temperature. The results also suggest that the discrepancy between bulk modulus of iron from experiments (Mao et al., 1990) and those from Earth’s core (Dziewonski and Anderson, 1981) is not originated from the overestimation of pressure by the EOS of Holmes et al. (1989). At high pressure and temperature, pressure uncertainty resulted from volume error becomes similarly important as the accuracy of the pressure scale.     

Regional Flow and Deformation Analysis of Basin‐Fill Aquifer Systems Using Stress‐Dependent Parameters

Changes in effective stress due to water pressure variations modify the intrinsic hydrodynamic properties of aquifers and aquitards. Overexploited groundwater systems, such as basins with heavy pumping, are subject to nonrecoverable modifications. This results in loss of permeability, porosity, and specific storage due to system consolidation. This paper presents (1) the analytical development of model functions relating effective stress to hydrodynamic parameters for aquifers and aquitards constituted of unconsolidated granular sediments, and (2) a modeling approach for the analysis of aquifer systems affected by effective stress variations, taking into account the aforementioned dependency. The stress‐dependent functions were fit to laboratory data, and used in the suggested modeling approach. Based on only few unknowns, this approach is computationally simple, efficiently captures the hydromechanical processes that are active in regional aquifer systems under stress, and readily provides an estimate of their consolidation.     

Simulation of breaking waves using the high-order spectral method with laboratory experiments: wave-breaking energy dissipation

We examine the implementation of a wave-breaking mechanism into a nonlinear potential flow solver. The success of the mechanism will be studied by implementing it into the numerical model HOS-NWT, which is a computationally efficient, open source code that solves for the free surface in a numerical wave tank using the high-order spectral (HOS) method. Once the breaking mechanism is validated, it can be implemented into other nonlinear potential flow models. To solve for wave-breaking, first a wave-breaking onset parameter is identified, and then a method for computing wave-breaking associated energy loss is determined. Wave-breaking onset is calculated using a breaking criteria introduced by Barthelemy et al. (J Fluid Mech https://arxiv.org/pdf/1508.06002.pdf, submitted) and validated with the experiments of Saket et al. (J Fluid Mech 811:642–658, 2017). Wave-breaking energy dissipation is calculated by adding a viscous diffusion term computed using an eddy viscosity parameter introduced by Tian et al. (Phys Fluids 20(6): 066,604, 2008, Phys Fluids 24(3), 2012), which is estimated based on the pre-breaking wave geometry. A set of two-dimensional experiments is conducted to validate the implemented wave breaking mechanism at a large scale. Breaking waves are generated by using traditional methods of evolution of focused waves and modulational instability, as well as irregular breaking waves with a range of primary frequencies, providing a wide range of breaking conditions to validate the solver. Furthermore, adjustments are made to the method of application and coefficient of the viscous diffusion term with negligible difference, supporting the robustness of the eddy viscosity parameter. The model is able to accurately predict surface elevation and corresponding frequency/amplitude spectrum, as well as energy dissipation when compared with the experimental measurements. This suggests the model is capable of calculating wave-breaking onset and energy dissipation successfully for a wide range of breaking conditions. The model is also able to successfully calculate the transfer of energy between frequencies due to wave focusing and wave breaking. This study is limited to unidirectional waves but provides a valuable basis for future application of the wave-breaking model to a multidirectional wave field. By including parameters for removing energy due to wave-breaking into a nonlinear potential flow solver, the risk of developing numerical instabilities due to an overturning wave is decreased, thereby increasing the application range of the model, including calculating more extreme sea states. A computationally efficient and accurate model for the generation of a nonlinear random wave field is useful for predicting the dynamic response of offshore vessels and marine renewable energy devices, predicting loads on marine structures, and in the study of open ocean wave generation and propagation in a realistic environment.     

Modelling stream flow for use in ecological studies in a large,arid zone river,central Australia

Australian arid zone ephemeral rivers are typically unregulated and maintain a high level of biodiversity and ecological health. Understanding the ecosystem functions of these rivers requires an understanding of their hydrology. These rivers are typified by highly variable hydrological regimes and a paucity, often a complete absence, of hydrological data to describe these flow regimes. A daily time‐step, grid‐based, conceptual rainfall–runoff model was developed for the previously uninstrumented Neales River in the arid zone of northern South Australia. Hourly, logged stage data provided a record of stream‐flow events in the river system. In conjunction with opportunistic gaugings of stream‐flow events, these data were used in the calibration of the model. The poorly constrained spatial variability of rainfall distribution and catchment characteristics (e.g. storage depths) limited the accuracy of the model in replicating the absolute magnitudes and volumes of stream‐flow events. In particular, small but ecologically important flow events were poorly modelled. Model performance was improved by the application of catchment‐wide processes replicating quick runoff from high intensity rainfall and improving the area inundated versus discharge relationship in the channel sections of the model. Representing areas of high and low soil moisture storage depths in the hillslope areas of the catchment also improved the model performance. The need for some explicit representation of the spatial variability of catchment characteristics (e.g. channel/floodplain, low storage hillslope and high storage hillslope) to effectively model the range of stream‐flow events makes the development of relatively complex rainfall–runoff models necessary for multisite ecological studies in large, ungauged arid zone catchments. Grid‐based conceptual models provide a good balance between providing the capacity to easily define land types with differing rainfall–runoff responses, flexibility in defining data output points and a parsimonious water‐balance–routing model. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydraulic conductivity prediction for sandy soils

The potential for petroleum-contaminated soils to impact ground water is evaluated using the soil leachability model in South Carolina and Georgia. In this model, the Green and Ampt (1911) equation is used to estimate unsaturated flow with saturated hydraulic conductivity (K(S)) values obtained using the Rawls and Brakensiek (1989) equation from inputs of percent sand- and clay-sized particles. However, many soils have > 70% sand-sized particles, which is the maximum amount for which the Rawls and Brakensiek equation is valid. Therefore, 70 sets of K(S) and particle-size data from the literature for southeastern United States sandy soils were analyzed to develop a new equation for estimating K(S). A multiple linear regression model with an adjusted R2 = 0.65 (p < 0.0001) was developed from percent clay- and sand-sized particle data. Eight additional sets of data were used to validate the model. The root mean square deviation and maximum squared difference, (yi - yi(2)), values for the new model are smaller than those obtained using the Rawls and Brakensiek equation, or Rosetta, a neural network-based model (Schaap 1999). The new model provides estimates that are within an order of magnitude for all but one of the test data sets. Rosetta predicts K(S) within one order of magnitude of measured values for six test data, and predicts K(S) within two orders of magnitude for the other two. The Rawls and Brakensiek predictions are within one order of magnitude for four test data, but are two or more orders of magnitude too high for the remaining four points. The new equation is recommended for estimating K(S) for sandy southeastern United States soils for which inputs are limited to percent sand and clay.     

EXPERIMENTAL STUDY ON AGGRADATION

１ＩＮＴＲＯＤＵＣＴＩＯＮＨｕｍａｎａｃｔｉｖｉｔｉｅｓｍａｙｉｎｄｕｃｅｓｉｇｎｉｆｉｃａｎｔｃｈａｎｇｅｓｉｎｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔｃａｐａｃｉｔｙｏｆａｌｕｖｉａｌｒｉｖｅｒｃｈａｎｎｅｌｓａｎｄｃａｕｓｅｃｈａｎｇｅｓｉｎｔｈｅ...     

PUMPING TESTS ON OPEN WELLS IN PALAR ALLUVIUM,NEAR MADRAS CITY,INDIA AN APPLICATION OF THE PAPADOPULOS-COOPER METHOD

Abstract

The non-equilibrium formula of Theis has hitherto been found to be inadequate to analyse drawdown data from open wells, due to the influence of effects of well-storage on early drawdown data. Recently, a reliable type-curve method was proposed by Papadopulos and Cooper to analyse drawdowns in wells of large diameter. In order to study its practical validity, the Papadopulos-Cooper method was applied to drawdown data from six open wells, tapping the alluvial sands of the river Palar, near Madras city, India. The study has indicated that the method has considerable practical utility in analysing drawdown data from large diameter open wells. Besides, the method also offers an easy way of computing the duration over which the effects of well-storage are likely to influence drawdown data, in terms of aquifer transmissibility and well-radius. If a rough idea of aquifer trans-missibility is available from existing data, then, such a computation, to determine the duration over which wellstorage affects drawdown, could be initially made to advantage, so that proper judgement could be exercised in using other conventional methods of analysis, such as that of Jacob or that of Theis. (Key Words: Pumping tests; Open wells; Papadopulos-Cooper method.)     

A Type-Curve Method for the Analysis of Pumping Tests with Piecewise-Linear Pumping Rates

Aquifer hydraulic parameters are commonly inferred from constant-rate pumping tests, while variable pumping rates are frequently encountered in actual field conditions. In this study, we propose a generally applicable dimensionless form of the analytical solution for variable-rate pumping tests in confined aquifers. In particular, we adopt a piecewise-linear fitting of variable pumping rates and propose a new type-curve method for estimating the hydraulic conductivity (K ) and specific storage (S_s ) of the investigated confined aquifer. For each test, a series of type curves, which depend on the variable pumping rates, the location of observation wells and the introduced first dimensionless inflection time, need to be provided for matching the observed drawdown data on a log-log graph. We first demonstrate the applicability and robustness of this method through a synthetic pumping test. Subsequently, we apply this method to analyze drawdown data from four pumping tests conducted within a multilayered aquifer/aquitard system in Wuxi city, Jiangsu Province, China. The parameter estimates are then compared with those reported by PEST. The K and S_s values estimated by the new type-curve method are found to be quite close to PEST-based estimates. Parameter estimation results demonstrate the difference in K and S_s values between observation wells. The difference could be attributed to the spatial heterogeneity in K and S_s . A future research topic may focus on the characterization of K and S_s heterogeneity with the currently available drawdown data from variable-rate pumping tests.