Effects of variations of river stage and hydraulic conductivity on temporal scaling of groundwater levels: numerical simulations

It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).     

Radially convergent groundwater flow in sloping terrain

Abstract

The groundwater flow equation governing the elevation (h) of the steady-state phreatic surface in a sloping aquifer fed by constant recharge over a bi-circular sector is rhh′ ? r ² Bh′ + Pr ² ? PR ² = 0, where r is the radial coordinate, P is a constant involving recharge and aquifer properties, and B is the slope of the aquifer—bedrock boundary. The derived flow equation describes radially convergent flow through a sloping aquifer that discharges to a water body of fixed head. One important simplification is that in which the width of the bi-circular sector is constant, and the draining land becomes a rectangular aquifer. The bi-circular sector and rectangular-strip groundwater flow problems are solved in terms of implicit equations. The solutions for the steady-state phreatic surfaces depend on the ratio of recharge to hydraulic conductivity, the slope of the aquifer-bedrock, and the downstream constant-head boundary. Computational examples illustrate the application of the solutions.     

Stochastically conditioned flow paths and travel time distribution in highly heterogeneous synthesized aquifers / Chemins d’écoulement et distribution des temps de parcours au sein d’aquifères synthétiques hautement hétérogènes sous conditions stochastiques

Abstract

Abstract Characterization of heterogeneity at the field scale generally requires detailed aquifer properties such as transmissivity and hydraulic head. An accurate delineation of these properties is expensive and time consuming, and for many if not most groundwater systems, is not practical. As an alternative approach, stochastic representation of random fields is used and presented in this paper. Specifically, an iterative stochastic conditional simulation approach was applied to a hypothetical and highly heterogeneous pre-designed aquifer system. The approach is similar to the classical co-kriging technique; it uses a linear estimator that depends on the covariance functions of transmissivity (T), and hydraulic head (h), as well as their cross-covariances. A linearized flow equation along with a conditional random field generator constitutes the iterative process of the conditional simulation. One hundred equally likely realizations of transmissivity fields with pre-specified geostatistical parameters were generated, and conditioned to both limited transmissivity and head data. The successful implementation of the approach resulted in conditioned flow paths and travel-time distribution under different degrees of aquifer heterogeneity. This approach worked well for fields exhibiting small variances. However, for random fields exhibiting large variances (greater than 1.0), an iterative procedure was used. The results show that, as the variance of the ln[T] increases, the flow paths tend to diverge, resulting in a wide spectrum of flow conditions, with no direct discernable relationship between the degree of heterogeneity and travel time. The applied approach indicates that high errors may result when estimation of particle travel times in a heterogeneous medium is approximated by an equivalent homogeneous medium.     

A numerical comparison between two upscaling techniques: non-local inverse based scaling and simplified renormalization

In this paper, we face the problem of upscaling transmissivity from the macroscopic to the megascopic scale; here the macroscopic scale is that of the continuous flow equations, whereas the megascopic scale is that of the flow models on a coarse grid. In this paper, we introduce the non-local inverse based scaling (NIBS) and compare it with the simplified renormalization (SR). The latter is a classical technique that we adapt to compute internode transmissivities for a finite differences flow model in a direct way. NIBS is implemented in three steps: in the first step, the macroscopic transmissivity, together with arbitrarily chosen auxiliary boundary conditions and sources, is used to solve forward problems (FPs) at the macroscopic scale; in the second step, the resulting heads are sampled at the megascopic scale; in the third step, the upscaled internode transmissivities are obtained by solving an inverse problem with the differential system method (DS) for which the heads resulting from the second step are used. NIBS is a non-local technique, because the computation of the internode transmissivities relies upon the whole transmissivity field at the macroscopic scale. We test NIBS against SR in the case of synthetic, isotropic, confined aquifers under the assumptions of two-dimensional (2D) and steady-state flow; the aquifers differ for the degree of heterogeneity, which is represented by a normally distributed uncorrelated component of lnT. For the comparison, the reference heads and fluxes at the megascopic scale are computed from the solution of FPs at the macroscopic scale. These reference values are compared with the heads and the fluxes predicted from models at the megascopic scale using the upscaled parameters of SR and NIBS. For the class of aquifers considered in this paper, the results of SR are better than those of NIBS, which hints that non-local effects can be disregarded at the megascopic scale. The two techniques provide comparable results when the heterogeneity increases, when the megascopic scale is large with respect to the heterogeneity length scale, or when the source terms are relevant.     

Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. I: laboratory measurements

Anisotropy and heterogeneity of hydraulic conductivity (K) are suspected of greatly affecting rates and patterns of ground‐water seepage in peats. A new laboratory method, termed here the modified cube method, was used to measure horizontal and vertical hydraulic conductivity (K_h and K_v) of 400 samples of bog peat. The new method avoids many of the problems associated with existing field and laboratory methods, and is shown to give relatively precise measurements of K. In the majority of samples tested, K_h was much greater than K_v, indicating that the bog peat was strongly anisotropic. Log₁₀K_h, log₁₀K_v, and log₁₀ (K_h/K_v) were found to vary significantly with depth, although none of the relationships was simple. We comment on the scale dependency of our measurements. Copyright © 2002 John Wiley & Sons, Ltd.     

A vortex-tube model of eddies in the inertial range

Abstract

A vortex-tube geometry of the cascade of energy to small-scale eddies, in the inertial range of fully-developed turbulence, is proposed. The model is a special case of the beta model of Frisch, Sulem and Nelkin (1978). We require that the cascade conserve the principal invariants of inviscid, incompressible flow, namely volume, topological knottedness, circulation, and, at discrete times marking the termination of steps in the cascade, energy. The process terminates in a finite time, as in any beta model, leaving behind a self-similar network of “inactive” tubes. We associate a self-similar scaling dimension D with the structure, equal to the Hausdorff dimension of the set of “active” tubes at the termination of the cascade. Because circulation Λ plays a key role in the analysis of the cascade, we refer to these vortex-tube geometries as “gamma models”. The viewpoint throughout is entirely deterministic.

We describe two examples of gamma models. In the ring geometry, an eddy is a vortex ring, and the cascade produces “rings upon rings”, so we allow cutting and fusing of tubes while conserving total helicity. In the preferred helical model, no cutting is needed, and the cascade produces an infinite progression of braided “coils upon coils”. We suggest that latter geometry as a candidate for the topology of a singularity of the inviscid limit of a Navier-Stokes flow, when modeled by discrete vortex tubes.

A crucial ingredient of a gamma model, not explicitly present in a beta model, is the possibility of “splitting” a vortex tube into sub-tubes carrying smaller circulation. We suggest a dynamical basis for this process, as an instability of tubes whose cores violate the Rayleigh criterion.

The parameters describing a gamma model are not uniquely determined by our study, but there is a “simplest” helical gamma model, involving minimal splitting and distortion of tubes. The dimension D of the structure is 13/5, with a scale factor Λ = 2^?5/4. This value of D agrees with that suggested by Hentschel and Procaccia (1982), by analogy with established results for certain branched polymers.     

Aquifer parameter estimation from surface resistivity data

This paper is devoted to the additional use, other than ground water exploration, of surface geoelectrical sounding data for aquifer hydraulic parameter estimation. In a mesoscopic framework, approximated analytical equations are developed separately for saline and for fresh water saturations. A few existing useful aquifer models, both for clean and shaley sandstones, are discussed in terms of their electrical and hydraulic effects, along with the linkage between the two. These equations are derived for insight and physical understanding of the phenomenon. In a macroscopic scale, a general aquifer model is proposed and analytical relations are derived for meaningful estimation, with a higher level of confidence, of hydraulic parameter from electrical parameters. The physical reasons for two different equations at the macroscopic level are explicitly explained to avoid confusion. Numerical examples from existing literature are reproduced to buttress our viewpoint.     

Comparison of approaches for predicting solute transport: sandbox experiments

The main purpose of this paper was to compare three approaches for predicting solute transport. The approaches include: (1) an effective parameter/macrodispersion approach (Gelhar and Axness 1983); (2) a heterogeneous approach using ordinary kriging based on core samples; and (3) a heterogeneous approach based on hydraulic tomography. We conducted our comparison in a heterogeneous sandbox aquifer. The aquifer was first characterized by taking 48 core samples to obtain local-scale hydraulic conductivity (K). The spatial statistics of these K values were then used to calculate the effective parameters. These K values and their statistics were also used for kriging to obtain a heterogeneous K field. In parallel, we performed a hydraulic tomography survey using hydraulic tests conducted in a dipole fashion with the drawdown data analyzed using the sequential successive linear estimator code (Yeh and Liu 2000) to obtain a K distribution (or K tomogram). The effective parameters and the heterogeneous K fields from kriging and hydraulic tomography were used in forward simulations of a dipole conservative tracer test. The simulated and observed breakthrough curves and their temporal moments were compared. Results show an improvement in predictions of drawdown behavior and tracer transport when the K tomogram from hydraulic tomography was used. This suggests that the high-resolution prediction of solute transport is possible without collecting a large number of small-scale samples to estimate flow and transport properties that are costly to obtain at the field scale.     

Density-dependent dispersion in heterogeneous porous media Part II: Comparison with nonlinear models

The results of a series of high-resolution numerical experiments are used to test and compare three nonlinear models for high-concentration-gradient dispersion. Gravity stable miscible displacement is considered. The first model, introduced by Hassanizadeh, is a modification of Fick’s law which involves a second-order term in the dispersive flux equation and an additional dispersion parameter β. The numerical experiments confirm the dependency of β on the flow rate. In addition, a dependency on travelled distance is observed. The model can successfully be applied to nearly homogeneous media (σ² = 0.1), but additional fitting is required for more heterogeneous media.The second and third models are based on homogenization of the local scale equations describing density-dependent transport. Egorov considers media that are heterogeneous on the Darcy scale, whereas Demidov starts at the pore-scale level. Both approaches result in a macroscopic balance equation in which the dispersion coefficient is a function of the dimensionless density gradient. In addition, an expression for the concentration variance is derived. For small σ², Egorov’s model predictions are in satisfactory agreement with the numerical experiments without the introduction of any new parameters. Demidov’s model involves an additional fitting parameter, but can be applied to more heterogeneous media as well.     

Spatial variability in streambed hydraulic conductivity of contrasting stream morphologies: channel bend and straight channel

Streambed hydraulic conductivity is one of the main factors controlling variability in surface water‐groundwater interactions, but only few studies aim at quantifying its spatial and temporal variability in different stream morphologies. Streambed horizontal hydraulic conductivities (K_h) were therefore determined from in‐stream slug tests, vertical hydraulic conductivities (K_v) were calculated with in‐stream permeameter tests and hydraulic heads were measured to obtain vertical head gradients at eight transects, each comprising five test locations, in a groundwater‐dominated stream. Seasonal small‐scale measurements were taken in December 2011 and August 2012, both in a straight stream channel with homogeneous elevation and downstream of a channel meander with heterogeneous elevation. All streambed attributes showed large spatial variability. K_h values were the highest at the depositional inner bend of the stream, whereas high K_v values were observed at the erosional outer bend and near the middle of the channel. Calculated K_v values were related to the thickness of the organic streambed sediment layer and also showed higher temporal variability than K_h because of sedimentation and scouring processes affecting the upper layers of the streambed. Test locations at the channel bend showed a more heterogeneous distribution of streambed properties than test locations in the straight channel, whereas within the channel bend, higher spatial variability in streambed attributes was observed across the stream than along the stream channel. Copyright © 2014 John Wiley & Sons, Ltd.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydraulic Tomography: Continuity and Discontinuity of High‐K and Low‐K Zones

Hydraulic tomography is an emerging field and modeling method that provides a continuous hydraulic conductivity (K) distribution for an investigated region. Characterization approaches that rely on interpolation between one‐dimensional (1D) profiles have limited ability to accurately identify high‐K channels, juxtapositions of lenses with high K contrast, and breaches in layers or channels between such profiles. However, locating these features is especially important for groundwater flow and transport modeling, and for design and operation of in situ remediation in complex hydrogeologic environments. We use transient hydraulic tomography to estimate 3D K in a volume of 15‐m diameter by 20‐m saturated thickness in a highly heterogeneous unconfined alluvial (clay to sand‐and‐gravel) aquifer with a K range of approximately seven orders of magnitude at an active industrial site in Assemini, Sardinia, Italy. A modified Levenberg‐Marquardt algorithm was used for geostatistical inversion to deal with the nonlinear nature of the highly heterogeneous system. The imaging results are validated with pumping tests not used in the tomographic inversion. These tests were conducted from three of five clusters of continuous multichannel tubing (CMTs) installed for observation in the tomographic testing. Locations of high‐K continuity and discontinuity, juxtaposition of very high‐K and very low‐K lenses, and low‐K “plugs” are evident in regions of the investigated volume where they likely would not have been identified with interpolation from 1D profiles at the positions of the pumping well and five CMT clusters. Quality assessment methods identified a suspect high‐K feature between the tested volume and a lateral boundary of the model.     

NEW WAYS FOR THE INTERPRETATION OF GEOELECTRICAL RESISTIVITY MEASUREMENTS IN THE SEARCH FOR AND DELIMITATION OF AQUIFERS

Abstract

Geoelectric resistivity measurements by means of direct current for solving hydrogeological problems have become increasingly significant in recent years. Measurements on the surface according to the four-point-method (SCHLUMBERGER or WENNER arrangement) result in so-called “apparent” resistivities ? α as a function of the electrode distance L. The evaluation of these measuring data, here in form of sounding graphs ? α(L/2), consists of the determination of true resistivities as a function of the depth z. Since a direct computation of ? (z) from α (L/2) is not possible in practice, theoretically computed master curves constitute the essential auxiliary means for the evaluation.

New simplified calculation techniques allow to establish accurately computed master curves for an underground consisting of more than three layers. By means of such standard graphs special problems of hydrogeology can quantitatively be solved by applying geoelectrical methods. The procedure is demonstrated on hand of complicated cases of aquifers devided into several storeys.     

Sediment and solute transport on soil slope under simultaneous influence of rainfall impact and scouring flow

Soil erosion and nutrient losses with surface runoff in the loess plateau in China cause severe soil quality degradation and water pollution. It is driven by both rainfall impact and runoff flow that usually take place simultaneously during a rainfall event. However, the interactive effect of these two processes on soil erosion has received limited attention. The objectives of this study were to better understand the mechanism of soil erosion, solute transport in runoff, and hydraulic characteristics of flow under the simultaneous influence of rainfall and shallow clear‐water flow scouring. Laboratory flume experiments with three rainfall intensities (0, 60, and 120 mm h⁻¹) and four scouring inflow rates (10, 20, 30, and 40 l min⁻¹) were conducted to evaluate their interactive effect on runoff. Results indicate that both rainfall intensity and scouring inflow rate play important roles on runoff formation, soil erosion, and solute transport in the surface runoff. A rainfall splash and water scouring interactive effect on the transport of sediment and solute in runoff were observed at the rainfall intensity of 60 mm h⁻¹ and scouring inflow rates of 20 l min⁻¹. Cumulative sediment mass loss (Ms) was found to be a linear function of cumulative runoff volume (Wr) for each treatment. Solute transport was also affected by both rainfall intensity and scouring inflow rate, and the decrease in bromide concentration in the runoff with time fitted to a power function well. Reynolds number (Re) was a key hydraulic parameter to determine erodability on loess slopes. The Darcy–Weisbach friction coefficients (f) decreased with the Reynolds numbers (Re), and the average soil and water loss rate (M_l) increased with the Reynolds numbers (Re) on loess slope for both scenarios with or without rainfall impact. Copyright © 2010 John Wiley & Sons, Ltd.     

The Navier–Stokes-α equations revisited

The Navier–Stokes-α equation is a regularised form of the Euler equation that has been employed in representing the sub-grid scales in large-eddy simulations. Determined efforts have been made to place it on a secure deductive foundation. This requires two steps to be completed. The first is fundamental and consists of establishing from the equations governing the fluid flow, a relationship between two velocities called by Holm (Chaos, 2002a, 12, 518) the “filtered” and “unfiltered” velocities. The second consists of the relation between these two velocities. Until now, the preferred route to the first objective has been variational, by varying the action using Hamilton's principle. Soward and Roberts (J. Fluid Mech., 2008, 604, 297) followed that variational route and established the existence of an important but unwelcome term omitted by Holm in his derivation. It is shown here that the Soward and Roberts result may be derived from Euler's equation by a direct approach with considerably greater efficiency. Holm achieved the second objective by making a “Taylor hypothesis”, which we use here to evaluate the unwelcome term missing from his analysis of the first step. The resulting model equations differ from those of Holm's α model, and the attractive mean Kelvin's circulation theorem that follows from his α equations is no longer valid. For that reason, we call the term omitted by Holm unwelcome.     

Rotating flow over long shallow ridges

Abstract

The flow of a rotating homogeneous, incompressible fluid past a long ridge is investigated. An analysis is presented for flows in which E ? 1, Ro ～ E^½, H/D ～ E⁰, h/D ～ E^½ and cosα ～ E⁰ where E is the Ekman number, Ro the Rossby number, H/D the fluid depth to ridge width ratio, h/D the ridge height to ridge width ratio and α the angle between the free stream flow and a line perpendicular to the ridge axis. The analysis includes effects of the nonlinear inertial terms. Particular examples of a ridge of triangular cross section and a sinusoidal topography are investigated in some detail. Experiments are presented for a triangular ridge which are in good agreement with the theory.     

Applicability of spectral analysis to determine hydraulic diffusivity

This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ^′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S _y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S _y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ^′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ^′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.     

A Correction on Coastal Heads for Groundwater Flow Models

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

Accounting for Aquifer Heterogeneity from Geological Data to Management Tools

A nested workflow of multiple‐point geostatistics (MPG) and sequential Gaussian simulation (SGS) was tested on a study area of 6 km² located about 20 km northwest of Quebec City, Canada. In order to assess its geological and hydrogeological parameter heterogeneity and to provide tools to evaluate uncertainties in aquifer management, direct and indirect field measurements are used as inputs in the geostatistical simulations to reproduce large and small‐scale heterogeneities. To do so, the lithological information is first associated to equivalent hydrogeological facies (hydrofacies) according to hydraulic properties measured at several wells. Then, heterogeneous hydrofacies (HF) realizations are generated using a prior geological model as training image (TI) with the MPG algorithm. The hydraulic conductivity (K) heterogeneity modeling within each HF is finally computed using SGS algorithm. Different K models are integrated in a finite‐element hydrogeological model to calculate multiple transport simulations. Different scenarios exhibit variations in mass transport path and dispersion associated with the large‐ and small‐scale heterogeneity respectively. Three‐dimensional maps showing the probability of overpassing different thresholds are presented as examples of management tools.     

Rain properties controlling soil splash detachment

Laboratory experiments have been conducted to study the effects of various rain properties on sand detachment resulting from raindrop impact. Splash cups were exposed to simulated rainfall intensities ranging between 10 and 140 mm h⁻¹. The detached sand was collected and weighed whereas rain intensity, equivalent drop diameter and fall velocity of raindrops were measured with an optical spectro‐pluviometer (OSP). The properties of the simulated rain (i.e. median volume diameter and kinetic energy) were compared with those observed in natural conditions. Statistical analyses have been undertaken in order to evaluate which rain variable best predicts the mass of sand detached. Linear and non‐linear correlations between the mass of detached sediment and the product of drop size (d) by drop velocity (v), i.e. D^αV^β, with values of α varying between 1 and 6 and β between 0 and 3, have been computed. The results indicate that the coefficient of determination (R²) for α ranging between 3 and 5 and β lower or equal to 2 are satisfying. Although kinetic energy (D³V²) described splash detachment relatively well, the product of momentum by drop diameter (D⁴V) was slightly superior in describing splash detachment. Therefore, the momentum multiplied by the drop diameter is recommended as the best rain variable to describe splash detachment. Copyright © 2000 John Wiley & Sons, Ltd.