Depth‐dependent hydraulic conductivity distribution patterns of a streambed

Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (K_v) of the channel sediments and the aquifer materials. This method was applied to produce K_v profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the K_v profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of K_v in the channel sediments. This depth‐dependent pattern of K_v distribution for the channel sediments is considered to be produced by hyporheic processes. This K_v‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulations of non-ergodic solute transport in three-dimensional heterogeneous porous media

Numerical simulations of non-ergodic transport of a non-reactive solute plume by steady-state groundwater flow under a uniform mean velocity, , were conducted in a three-dimensional heterogeneous and statistically isotropic aquifer. The hydraulic conductivity, K(x), is modeled as a random field which is assumed to be log-normally distributed with an exponential covariance. Significant efforts are made to reduce the simulation uncertainties. Ensemble averages of the second spatial moments of the plume and the plume centroid variances were simulated with 1600 Monte Carlo (MC) runs for three variances of log K, _Y²=0.09, 0.23, and 0.46, and a square source normal to of three dimensionless lengths. It is showed that 1600 MC runs are needed to obtain stabilized results in mildly heterogeneous aquifers of _Y²0.5 and that large uncertainty may exist in the simulated results if less MC runs are used, especially for the transverse second spatial moments and the plume centroid variance in transverse directions. The simulated longitudinal second spatial moment and the plume centroid variance in longitudinal direction fit well to the first-order theoretical results while the simulated transverse moments are generally larger than the first-order values. The ergodic condition for the second spatial moments is far from reaching in all cases simulated and transport in transverse directions may reach ergodic condition much slower than that in longitudinal direction.     

The elastic properties of single-crystal fayalite as determined by dynamical measurement techniques

We present new elasticity measurements on single-crystal fayalite and combine our results with other data from resonance, pulse superposition interferometry, and Brillouin scattering to provide a set of recommended values for the adiabatic elastic moduliC _ij and their temperature variations. We use a resonance method (RPR) with specimens that were previously investigated by pulse superposition experiments. The nineC _ij of fayalite are determined from three new sets of measurements. One set of our newC _ij data is over the range 300–500 K. We believe that the relatively large discrep ncies found in someC _ij are due in large part to specimen inhomogeneities (chemical and microstructural) coupled with differences in the way various techniques sample, rather than only systematic errors associated with experimental procedures or in the preparations of the specimens.Our recommendeaC _ij's (GPa) and (C _ij/T) _p (GPa/K) are: The resulting values for the isotropic bulk and shear moduli,K _s and , and their temperature derivatives are:K _s=134(4) GPa; =50.7(0.3) GPa; (K _s/T) _p =–0.024(0.005) GPa/K; and (/T) _p =–0.013(0.001) GPa/K. An important conclusion is thatK _s increases as the Fe/(Fe+Mg) ratio in olivine is increased.     

Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. II: modelling the effects on groundwater flow

Anisotropy and heterogeneity of hydraulic conductivity (K) are seldom considered in models of mire hydrology. We investigated the effect of anisotropy and heterogeneity on groundwater flow in bog peat using a steady‐state groundwater model. In five model simulations, four sets of K data were used. The first set comprised measured K values from an anisotropic and heterogeneous bog peat. These data were aggregated to produce the following simplified data sets: an isotropic and heterogeneous distribution of K; an isotropic and homogeneous distribution; and an anisotropic and homogeneous distribution. We demonstrate that, where anisotropy and heterogeneity exist, groundwater flow in bog peat is complex. Fine‐scale variations in K have the potential to influence patterns and rates of groundwater flow. However, for our data at least, it is heterogeneity and not anisotropy that has the greater influence on producing complex patterns of groundwater flow. We also demonstrate that patterns and rates of groundwater flow are simplified and reduced when measured K values are aggregated to create a more uniform distribution of K. For example, when measured K values are aggregated to produce isotropy and homogeneity, the rate of modelled seepage is reduced by 28%. We also show that when measured K values are used, the presence of a drainage ditch can increase seepage through a modelled cross‐section. Our work has implications for the accurate interpretation of hydraulic head data obtained from peat soils, and also the understanding of the effect of drainage ditches on patterns and rates of groundwater flow. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of anisotropy on solute transport in degraded fen peat soils

Peat soils are heterogeneous, anisotropic porous media. Compared to mineral soils, there is still limited understanding of physical and solute transport properties of fen peat soils. In this study, we aimed to explore the effect of soil anisotropy on solute transport in degraded fen peat. Undisturbed soil cores, taken in vertical and horizontal direction, were collected from one drained and one restored fen peatland both in a comparable state of soil degradation. Saturated hydraulic conductivity (K _s) and chemical properties of peat were determined for all soil cores. Miscible displacement experiments were conducted under saturated steady state conditions using potassium bromide as a conservative tracer. The results showed that (1) the K _s in vertical direction (K _sv) was significantly higher than that in horizontal direction (K_sh), indicating that K _s of degraded fen peat behaves anisotropically; (2) pronounced preferential flow occurred in vertical direction with a higher immobile water fraction and a higher pore water velocity; (3) the 5% arrival time (a proxy for the strength of preferential flow) was affected by soil anisotropy as well as study site. A strong correlation was found between 5% arrival time and dispersivity, K _s and mobile water fraction; (4) phosphate release was observed from drained peat only. The impact of soil heterogeneity on phosphate leaching was more pronounced than soil anisotropy. The soil core with the strongest preferential flow released the highest amount of phosphate. We conclude that soil anisotropy is crucial in peatland hydrology but additional research is required to fully understand anisotropy effects on solute transport.     

Flow and transport through two-dimensional isotropic media of binary conductivity distribution. Part 1: NUMERICAL methodology and semi-analytical solutions

Flow and transport take place in a heterogeneous medium made up from inclusions of conductivity K submerged in a matrix of conductivity K ₀. We consider two-dimensional isotropic media, with circular inclusions of uniform radii, that are placed at random and without overlap in the matrix. The system is completely characterized by the conductivity contrast =K/K ₀ and by the volume fraction n. The flow is uniform in the mean, of velocity U=const. The derivation of the velocity field is achieved by a numerical method of high accuracy, based on analytical elements. Approximate analytical solutions are derived by a few methods: composite elements, effective medium, dilute systems and first-order approximation in logconductivity variance. The latter was employed by Rubin (1995), while the dilute system approximation was used by Eames and Bush (1999) and Dagan and Lessoff (2001). Transport is solved in a Lagrangean framework, with trajectories determined numerically from the velocity field, by particle tracking. Results for the velocity variance and for the longitudinal macrodispersivity, for a few values of and n, are presented in Part 2.     

Reservoir lifetime and heat recovery factor in geothermal aquifers used for urban heating

Simple models are discussed to evaluate reservoir lifetime and heat recovery factor in geothermal aquifers used for urban heating. By comparing various single well and doublet production schemes, it is shown that reinjection of heat depleted water greatly enhances heat recovery and reservoir lifetime, and can be optimized for maximum heat production. It is concluded that geothermal aquifer production should be unitized, as is already done in oil and gas reservoirs.Nomenclature a distance between doublets in multi-doublet patterns, meters - A area of aquifer at base temperature, m² drainage area of individual doublets in multidoublet patterns, m² - D distance between doublet wells, meters - h aquifer thickness, meters - H water head, meters - Q production rate, m³/sec. - r _e aquifer radius, meters - r _w well radius, meters - R _g heat recovery factor, fraction - S water level drawdown, meters - t producing time, sec. - T aquifer transmissivity, m²/sec. - v stream-channel water velocity, m/sec. - actual temperature change, °C - theoretical temperature change, °C - water temperature, °C - heat conductivity, W/m/°C - _r rock heat conductivity, W/m/°C - _aC_a aquifer heat capacity, J/m³/°C - _aC_r rock heat capacity, J/m³/°C - _WC_W water heat capacity, J/m³/°C - aquifer porosity, fraction     

Stochastic analysis of flux and head moments in a heterogeneous aquifer system

This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈q_x〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈q_y〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈q_x〉 field; finally, 〈q_z〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.     

Changes in Entrapped Gas Content and Hydraulic Conductivity with Pressure

Water table fluctuations continuously introduce entrapped air bubbles into the otherwise saturated capillary fringe and groundwater zone, which reduces the effective (quasi‐saturated) hydraulic conductivity, K_quasi, thus impacting groundwater flow, aquifer recharge and solute and contaminant transport. These entrapped gases will be susceptible to compression or expansion with changes in water pressure, as would be expected with water table (and barometric pressure) fluctuations. Here we undertake laboratory experiments using sand‐packed columns to quantify the effect of water table changes of up to 250 cm on the entrapped gas content and the quasi‐saturated hydraulic conductivity, and discuss our ability to account for these mechanisms in ground water models. Initial entrapped air contents ranged between 0.080 and 0.158, with a corresponding K_quasi ranging between 2 and 6 times lower compared to the K_s value. The application of 250 cm of water pressure caused an 18% to 26% reduction in the entrapped air content, resulting in an increase in K_quasi by 1.16 to 1.57 times compared to its initial (0 cm water pressure) value. The change in entrapped air content measured at pressure step intervals of 50 cm, was essentially linear, and could be modeled according to the ideal gas law. Meanwhile, the changes in K_quasi with compression–expansion of the bubbles because of pressure changes could be adequately captured with several current hydraulic conductivity models.     

Investigating peat hydrological properties using field and laboratory methods: application to the Lanoraie peatland complex (southern Quebec,Canada)

Flow dynamics within a peatland are governed by hydraulic parameters such as hydraulic conductivity, dispersivity and specific yield, as well as by anisotropy and heterogeneity. The aim of this study is to investigate hydraulic parameters variability in peat through the use of different field and laboratory methods. An experimental site located in the Lanoraie peatland complex (southern Quebec, Canada) was used to test the different approaches. Slug and bail tests were performed in piezometer standpipes to investigate catotelm hydraulic conductivity. Combined Darcy tests and tracer experiments were conducted on cubic samples using the modified cube method (MCM) to assess catotelm hydraulic conductivity, anisotropy and dispersivity. A new laboratory method is proposed for assessing acrotelm hydraulic conductivity and gravity drainage using a laboratory experimental tank. Most of slug tests' recovery curves were characteristic of compressible media, and important variability was observed depending on the initial head difference. The Darcy experiments on cubic samples provided reproducible results, and anisotropy (K_h > K_v) was observed for most of samples. All tracer experiments displayed asymmetrical breakthrough curves, suggesting the presence of retardation and/or dual porosity. Hydraulic conductivity estimates performed using the experimental tank showed K variations over a factor of 44 within the upper 40 cm of the acrotelm. The results demonstrate that the intrinsic variability associated with the different field and laboratory methods is small compared with the spatial variability of hydraulic parameters. It is suggested that a comprehensive assessment of peat hydrological properties can be obtained through the combined use of complementary field and laboratory investigations. Copyright © 2007 John Wiley & Sons, Ltd.     

Field Study of Subsurface Heterogeneity with Steady‐State Hydraulic Tomography

Remediation of subsurface contamination requires an understanding of the contaminant (history, source location, plume extent and concentration, etc.), and, knowledge of the spatial distribution of hydraulic conductivity (K) that governs groundwater flow and solute transport. Many methods exist for characterizing K heterogeneity, but most if not all methods require the collection of a large number of small‐scale data and its interpolation. In this study, we conduct a hydraulic tomography survey at a highly heterogeneous glaciofluvial deposit at the North Campus Research Site (NCRS) located at the University of Waterloo, Waterloo, Ontario, Canada to sequentially interpret four pumping tests using the steady‐state form of the Sequential Successive Linear Estimator (SSLE) ( Yeh and Liu 2000 ). The resulting three‐dimensional (3D) K distribution (or K‐tomogram) is compared against: ( 1 ) K distributions obtained through the inverse modeling of individual pumping tests using SSLE, and ( 2 ) effective hydraulic conductivity (K_eff) estimates obtained by automatically calibrating a groundwater flow model while treating the medium to be homogeneous. Such a K_eff is often used for designing remediation operations, and thus is used as the basis for comparison with the K‐tomogram. Our results clearly show that hydraulic tomography is superior to the inversions of single pumping tests or K_eff estimates. This is particularly significant for contaminated sites where an accurate representation of the flow field is critical for simulating contaminant transport and injection of chemical and biological agents used for active remediation of contaminant source zones and plumes.     

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Flow and transport through two-dimensional isotropic media of binary conductivity distribution. Part 2: NUMERICAL simulations and comparison with theoretical results

In the present part the results of numerical simulations of flow and transport in media made up from circular inclusions of conductivity K that are submerged in a matrix of conductivity K ₀, subjected to uniform mean velocity, are presented. This is achieved for a few values of =K/K ₀ (0.01, 0.1 and 10) and of the volume fraction n (0.05, 0.1 and 0.2). The numerical simulations (NS) are compared with the analytical approximate models presented in Part 1: the composite elements (CEA), the effective medium (EMA), the dilute system (DSA) and the first-order in the logconductivity variance (FOA). The comparison is made for the longitudinal velocity variance and for the longitudinal macrodispersivity. This is carried out for n<0.2, for which the theoretical and simulation models represent the same structure of random and independent inclusions distribution. The main result is that transport is quite accurately modeled by the EMA and CEA for low , for which _L is large, whereas in the case of =10, the EMA matches the NS for n<0.1. The first-order approximation is quite far apart from the NS for the values of examined . This material is based upon work supported by the National Science Foundation under Grant No. 0218914. Authors also wish to thank the Center of Computational Research, University at Buffalo for assistance in running numerical simulations.     

Stochastic analysis of groundwater flow in aquifers with leakage

 Stochastic analysis of one- and two-dimensional flow through a shallow semi-confined aquifer with spatially variable hydraulic conductivity K represented by a stationary (statistically homogeneous) random process is carried out by using the spectral technique. The hydraulic head covariance functions for flows in a semi-confined aquifer bounded by a leaky layer above and an impervious stratum below are derived by assuming that the randomness forcing the head variation to originate from the hydraulic conductivity field of the aquifer. The head covariance functions are studied using two convenient forms of the logarithmic hydraulic conductivity process. The results demonstrate the significant reduction in the head variances and covariances due to the presence of a leaky layer. The hydraulic head correlation distance is also reduced greatly due to the presence of the leaky layer.     

Coordinated Simultaneous Freja Satellite and Niemegk Observatory ULF Measurements

The standard FFT analysis was applied to thirteen pulsation events selected from the March–April 1993 simultaneous measurements in space (Freja satellite) and on the ground (Niemegk Observatory). The spectral processing of the six-minute sections of record in two coordinate components perpendicular to the magnetic field lines had two principal subjects in view: a) Frequency-amplitude satellite-ground relations. The average values of smoothed spectral amplitude transmission coefficients in the total Pc4-3 and, separately, in Pc4 and Pc3 bands were 0.60, 0.48 and 0.70, respectively. The maximum values were observed within the Pc3 band, near f 30 mHz. b) Estimation of frequency bands of field line resonances (FLRs) recorded on Freja. Freja's motion in a broad latitude range (±5° around Niemegk) during the measurements enabled two frequency bands of FLRs with a width of 10's of mHz to be identified on Freja. The weighted frequencies of the FLRs bands on Freja were f _wA 36 mHz and f _wB 18 mHz with fine structure separations of about 5 mHz and 4 mHz, respectively. The second band with f _wB 18 mHz was mostly not observed at Niemegk, and this could be the consequence of the satellite moving partly through regions of L greater than that of Niemegk.     

Incorporating spatially heterogeneous infiltration capacity into hydrologic models with applications for simulating post‐wildfire debris flow initiation

Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity (K_s) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of K_s, and the spatial correlation length scale associated with K_s. The proposed method for determining is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale K_s measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.     

Theory of the piezomagnetic effect in titanomagnetite-bearing rocks

Summary The stress sensitivitiesS _x andS _R of susceptibility and remanence for titanomagnetite-bearing rocks are calculated in terms of magnetostriction constants ₁₀₀ and ₁₁₁ and anisotropy constantsK ₁,K ₂ of the magnetic minerals.S _x andS _R are represented by quite different algebraic expressions but happen to have comparable numerical values over the whole range of titanomagnetite compositions. Both increase strongly with titanium content. This leads to more optimistic calculations of tectonomagnetic effects than with the previously assumed stress sensitivity for pure magnetite.     

Steady flows driven by sources of random strength in heterogeneous aquifers with application to partially penetrating wells

Average steady source flow in heterogeneous porous formations is modelled by regarding the hydraulic conductivity K(x) as a stationary random space function (RSF). As a consequence, the flow variables become RSFs as well, and we are interested into calculating their moments. This problem has been intensively studied in the case of a Neumann type boundary condition at the source. However, there are many applications (such as well-type flows) for which the required boundary condition is that of Dirichlet. In order to fulfill such a requirement the strength of the source must be proportional to K(x), and therefore the source itself results a RSF. To solve flows driven by sources whose strength is spatially variable, we have used a perturbation procedure similar to that developed by Indelman and Abramovich (Water Resour Res 30:3385–3393, 1994) to analyze flows generated by sources of deterministic strength. Due to the linearity of the mathematical problem, we have focused on the explicit derivation of the mean head distribution G _d (x) generated by a unit pulse. Such a distribution represents the fundamental solution to the average flow equations, and it is termed as mean Green function. The function G _d (x) is derived here at the second order of approximation in the variance σ² of the fluctuation (where K _A is the mean value of K(x)), for arbitrary correlation function ρ(x), and any dimensionality d of the flow domain. We represent G _d (x) as product between the homogeneous Green function G _d ⁽⁰⁾(x) valid in a domain with constant K _A , and a distortion term Ψ _d (x) = 1 + σ²ψ _d (x) which modifies G _d ⁽⁰⁾(x) to account for the medium heterogeneity. In the case of isotropic formations ψ _d (x) is expressed via one quadrature. This quadrature can be analytically calculated after adopting specific (e.g.. exponential and Gaussian) shape for ρ(x). These general results are subsequently used to investigate flow toward a partially-penetrating well in a semi-infinite domain. Indeed, we construct a σ²-order approximation to the mean as well as variance of the head by replacing the well with a singular segment. It is shown how the well-length combined with the medium heterogeneity affects the head distribution. We have introduced the concept of equivalent conductivity K ^eq(r,z). The main result is the relationship where the characteristic function ψ^(w)(r,z) adjusts the homogeneous conductivity K _A to account for the impact of the heterogeneity. In this way, a procedure can be developed to identify the aquifer hydraulic properties by means of field-scale head measurements. Finally, in the case of a fully penetrating well we have expressed the equivalent conductivity in analytical form, and we have shown that (being the effective conductivity for mean uniform flow), in agreement with the numerical simulations of Firmani et al. (Water Resour Res 42:W03422, 2006).     

Soil hydraulic conductivities of latosols under pasture,forest and teak in Rondonia,Brazil

We investigated the changes of saturated hydraulic conductivity, K_sat , with depth of latosols developed on Precambrian basement rocks under primary rainforest, pasture and teak. In all cases, K_sat decreased with depth, with most of the decrease occurring between the surface and a depth of 30 cm. In conjunction with prevailing rainfall intensities and frequencies, this anisotropy supports a pronounced lateral component of hillslope flow paths, and also of overland flow under pasture. Our results are at variance with data from other latosols where K_sat tends to increase with depth, and hence suggest that considerable restraint is needed in generalization and extrapolation until results from a co‐ordinated effort at hydrology‐oriented data collection become available. Copyright © 1999 John Wiley & Sons, Ltd.     

Sensitivity analysis of a no-crossflow model for the transient flowmeter test

 Logarithmic sensitivities and plausible relative errors are studied in a simple no-crossflow model of a transient flowmeter test (TFMT). This model is identical to the model of a constant-rate pumping test conducted on a fully penetrating well with wellbore storage, surrounded by a thick skin zone, and situated in a homogeneous confined aquifer. The sensitivities of wellbore drawdown and wellface flowrate to aquifer and skin parameters are independent of the pumping rate. However, the plausible relative errors in the aquifer and skin parameters estimated from drawdown and wellface flowrate data can be proportionally decreased by increasing the pumping rate. The plausible relative errors vary by many orders of magnitude from the beginning of the TFMT. The practically important flowrate and drawdown measurements in this test, for which the plausible relative errors vary by less than one order of magnitude from the minimum plausible relative errors, can begin approximately when the dimensionless wellface flowrate exceeds q _D =q/Q≈0.4. During most of this stage of the test, the plausible relative errors in aquifer hydraulic conductivity (K _a ) are generally an order of magnitude smaller than those in aquifer specific storativity. The plausible relative errors in the skin hydraulic conductivity (K _s ) are generally larger than the plausible relative errors in the aquifer specific storativity when the thick skin is normal (K _s >K _a ) and smaller when the thick skin is damaged (K _s <K _a ). The specific storativity of the skin zone would be so biased that one should not even attempt to estimate it from the TFMT. We acknowledge Wiebe H. van der Molen for recommending the De Hoog algorithm and sharing his code. This research was partially supported by the US Geological Survey, USGS Agreement #1434-HQ-96-GR-02689 and North Carolina Water Resources Research Institute, WRRI Project #70165.