Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Water Resources Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17Linked to original sources
USGS2000

Timescales for migration of atmospherically derived sulphate through an alpine/subalpine watershed, Loch Vale, Colorado

Sulphur 35, a cosmogenically produced radioisotope with a short half‐life (87 days), was measured in snowpack during 1993–1997 and at four locations within the Loch Vale watershed during 1995–1997. The four sites include the two main drainages in the watershed, Andrews Creek and Icy Brook, a small south facing catchment flowing into Andrews Creek (Andrews Spring 1), and a similar north facing catchment flowing out of a scree field into Icy Brook (Spring 19). Concentrations ranged from a high of almost 50 mBq/L for a sample from Spring 19 in June 1996 to a concentration near the detection limit for a sample from Andrews Creek in April 1997. Sulphur 35 concentrations were normalized to sulphate (as mBq/mg SO 4 −2 ) and were decay‐corrected to a Julian day of 90 (April 1) for each year. Snowpack had the highest 35 S concentration with an average concentration of 53 mBq/mg SO 4 −2 . Concentrations in the streams were much lower, even when corrected for decay relative to JD 90. The large 35 S concentrations found in Spring 19 were the result of increases in concentration due to sublimation and/or evapotranspiration and were lower than snowpack when normalized to sulphate. Using 35 S concentrations found in snowpack as of JD 90 as a beginning concentration, the fraction of sulphate in streamflow that was derived from atmospheric deposition within the prior water year was estimated. For Icy Brook and Andrews Creek the fraction of the sulphate in streamflow derived from that year's snowpack and precipitation was low prior to the beginning of the main spring melt, reached a maximum during the period of maximum flow, and decreased as the summer progressed. A calculation of the seasonal flux indicated that about 40% of the sulphate that flowed out of the watershed was derived from atmospheric sulphate deposited during the previous year. This suggests that more than half of the sulphate deposited in the watershed by atmospheric processes during the previous year was removed during the following summer. Thus sulphate retention in alpine watersheds like Loch Vale is very limited, and changes in sulphate deposition should be quickly reflected in stream chemistry.

Colorado↗
USGS1999

Influence of various water quality sampling strategies on load estimates for small streams

Extensive streamflow and water quality data from eight small streams were systematically subsampled to represent various water‐quality sampling strategies. The subsampled data were then used to determine the accuracy and precision of annual load estimates generated by means of a regression approach (typically used for big rivers) and to determine the most effective sampling strategy for small streams. Estimation of annual loads by regression was imprecise regardless of the sampling strategy used; for the most effective strategy, median absolute errors were ∼30% based on the load estimated with an integration method and all available data, if a regression approach is used with daily average streamflow. The most effective sampling strategy depends on the length of the study. For 1‐year studies, fixed‐period monthly sampling supplemented by storm chasing was the most effective strategy. For studies of 2 or more years, fixed‐period semimonthly sampling resulted in not only the least biased but also the most precise loads. Additional high‐flow samples, typically collected to help define the relation between high streamflow and high loads, result in imprecise, overestimated annual loads if these samples are consistently collected early in high‐flow events.

Water Resources Research↗
USGS2000

Multispecies reactive tracer test in an aquifer with spatially variable chemical conditions

A field investigation of multispecies reactive transport was conducted in a well‐characterized, sand and gravel aquifer on Cape Cod, Massachusetts. The aquifer is characterized by regions of differing chemical conditions caused by the disposal of secondary sewage effluent. Ten thousand liters of groundwater with added tracers (Br, Cr(VI), and EDTA complexed with Pb, Zn, Cu, and Ni) were injected into the aquifer and distributions of the tracers were monitored for 15 months. Most of the tracers were transported more than 200 m; transport was quantified using spatial moments computed from the results of a series of synoptic samplings. Cr(VI) transport was retarded relative to Br; the retardation factor varied from 1.1 to 2.4 and was dependent on chemical conditions. At 314 days after the injection, dissolved Cr(VI) mass in the tracer cloud had decreased 85%, with the likely cause being reduction to Cr(III) in a suboxic region of the aquifer. Transport of the metal‐EDTA complexes was affected by aqueous complexation, adsorption, and dissolution‐precipitation reactions of Fe oxyhydroxide minerals in the aquifer sediments. Dissolved Pb‐EDTA complexes disappeared from the tracer cloud within 85 days, probably due to metal exchange reactions with Fe and adsorbed Zn (present prior to the injection from contamination by the sewage effluent). About 30% of the Cu‐EDTA complexes remained within the tracer cloud 314 days after injection, even though the thermodynamic stability of the Pb‐EDTA complex is greater than Cu‐EDTA. It is hypothesized that stronger adsorption of Pb 2+ to the aquifer sediments causes the Pb‐EDTA complex to disassociate to a greater degree than the Cu‐EDTA complex. The mass of dissolved Zn‐EDTA increased during the first 175 days of the tracer test to 140% of the mass injected, with the increase due to desorption of sewage‐derived Zn. Dissolved Ni‐EDTA mass remained nearly constant throughout the tracer test, apparently only participating in reversible adsorption reactions. The results of the field experiment provide a chemically complex data set that can be used in the testing of reactive transport models of flow coupled with chemical reactions.

Water Resources Research↗
USGS2000

Controls on nitrogen flux in alpine/subalpine watersheds of Colorado

High‐altitude watersheds in the Front Range of Colorado show symptoms of advanced stages of nitrogen excess, despite having less nitrogen in atmospheric deposition than other regions where watersheds retain nitrogen. In two alpine/subalpine subbasins of the Loch Vale watershed, atmospheric deposition of NO 3 − plus NH 4 + was 3.2–5.5 kg N ha −1 , and watershed export was 1.8–3.9 kg N ha −1 for water years 1992–1997. Annual N export increased in years with greater input of N, but most of the additional N was retained in the watershed, indicating that parts of the ecosystem are nitrogen‐limited. Dissolved inorganic nitrogen (DIN) concentrations were greatest in subsurface water of talus landscapes, where mineralization and nitrification augment high rates of atmospheric deposition of N. Tundra landscapes had moderately high DIN concentrations, whereas forest and wetland landscapes had low concentrations, indicating little export of nitrogen from these landscapes. Between the two subbasins the catchment of Icy Brook had greater retention of nitrogen than that of Andrews Creek because of landscape and hydrologic characteristics that favor greater N assimilation in both the terrestrial and aquatic ecosystems. These results suggest that export of N from alpine/subalpine watersheds is caused by a combination of direct flushing of N from atmospheric deposition and release of N from ecosystem biogeochemical processes (N cycling). Sensitivity of alpine ecosystems in the western United States to atmospheric deposition of N is a function of landscape heterogeneity, hydrologic flow paths, and climatic extremes that limit primary productivity and microbial activity, which, in turn, control retention and release of nitrogen. Conceptual and mechanistic models of N excess that have been developed for forested ecosystems need to be modified in order to predict the response of alpine ecosystems to future changes in climate and atmospheric deposition of N.

Water Resources Research↗
USGS2000

Colorado River sediment transport: 1. Natural sediment supply limitation and the influence of Glen Canyon Dam

Analyses of flow, sediment‐transport, bed‐topographic, and sedimentologic data suggest that before the closure of Glen Canyon Dam in 1963, the Colorado River in Marble and Grand Canyons was annually supply‐limited with respect to fine sediment (i.e., sand and finer material). Furthermore, these analyses suggest that the predam river in Glen Canyon was not supply‐limited to the same degree and that the degree of annual supply limitation increased near the head of Marble Canyon. The predam Colorado River in Grand Canyon displays evidence of four effects of supply limitation: (1) seasonal hysteresis in sediment concentration, (2) seasonal hysteresis in sediment grain size coupled to the seasonal hysteresis in sediment concentration, (3) production of inversely graded flood deposits, and (4∥ development or modification of a lag between the time of a flood peak and the time of either maximum or minimum (depending on reach geometry) bed elevation. Analyses of sediment budgets provide additional support for the interpretation that the predam river was annually supply‐limited with respect to fine sediment, but it was not supply‐limited with respect to fine sediment during all seasons. In the average predam year, sand would accumulate and be stored in Marble Canyon and upper Grand Canyon for 9 months of the year (from July through March) when flows were dominantly below 200–300 m 3 /s; this stored sand was then eroded during April through June when flows were typically higher. After closure of Glen Canyon Dam, because of the large magnitudes of the uncertainties in the sediment budget, no season of substantial sand accumulation is evident. Because most flows in the postdam river exceed 200–300 m 3 /s, substantial sand accumulation in the postdam river is unlikely.

Colorado River, Glen Canyon Dam↗
USGS2000

Colorado River sediment transport: 2. Systematic bed‐elevation and grain‐size effects of sand supply limitation

The Colorado River in Marble and Grand Canyons displays evidence of annual supply limitation with respect to sand both prior to [ Topping et al , this issue] and after the closure of Glen Canyon Dam in 1963. Systematic changes in bed elevation and systematic coupled changes in suspended‐sand concentration and grain size result from this supply limitation. During floods, sand supply limitation either causes or modifies a lag between the time of maximum discharge and the time of either maximum or minimum (depending on reach geometry) bed elevation. If, at a cross section where the bed aggrades with increasing flow, the maximum bed elevation is observed to lead the peak or the receding limb of a flood, then this observed response of the bed is due to sand supply limitation. Sand supply limitation also leads to the systematic evolution of sand grain size (both on the bed and in suspension) in the Colorado River. Sand input during a tributary flood travels down the Colorado River as an elongating sediment wave, with the finest sizes (because of their lower settling velocities) traveling the fastest. As the fine front of a sediment wave arrives at a given location, the bed fines and suspended‐sand concentrations increase in response to the enhanced upstream supply of finer sand. Then, as the front of the sediment wave passes that location, the bed is winnowed and suspended‐sand concentrations decrease in response to the depletion of the upstream supply of finer sand. The grain‐size effects of depletion of the upstream sand supply are most obvious during periods of higher dam releases (e.g., the 1996 flood experiment and the 1997 test flow). Because of substantial changes in the grain‐size distribution of the bed, stable relationships between the discharge of water and sand‐transport rates (i.e., stable sand rating curves) are precluded. Sand budgets in a supply‐limited river like the Colorado River can only be constructed through inclusion of the physical processes that couple changes in bed‐sediment grain size to changes in sand‐transport rates.

Colorado River, Grand Canyon↗
USGS2000

Relations between basin characteristics and stream water chemistry in alpine/subalpine basins in Rocky Mountain National Park, Colorado

Relations between stream water chemistry and topographic, vegetative, and geologic characteristics of basins were evaluated for nine alpine/subalpine basins in Rocky Mountain National Park, Colorado, to identify controlling parameters and to better understand processes governing patterns in stream water chemistry. Fractional amounts of steep slopes (≥30°), unvegetated terrain, and young surficial debris within each basin were positively correlated to each other. These terrain features, which commonly occur on steep valley side slopes underlain by talus, were negatively correlated with concentrations of base cations, silica, and alkalinity and were positively correlated with nitrate, acidity, and runoff. These relations might result from the short residence times of water and limited soil development in the talus environment, which limit chemical weathering and nitrogen uptake. Steep, unvegetated terrains also tend to promote high Ca/Na ratios in stream water, probably because physical weathering rates in those areas are high. Physical weathering exposes fresh bedrock that contains interstitial calcite, which weathers relatively quickly. The fractional amounts of subalpine meadow and, to a lesser extent, old surficial debris in the basins were positively correlated to concentrations of weathering products and were negatively correlated to nitrate and acidity. These relations may reflect more opportunities for silicate weathering and nitrogen uptake in the lower‐energy environments of the valley floor, where soils are finer‐grained, older, and better developed and slopes are relatively flat. These results indicate that in alpine/subalpine basins, slope, vegetation (or lack thereof), and distribution and age of surficial materials are interrelated and can have major effects on stream water chemistry.

Colorado↗
USGS2000

Development of a pore network simulation model to study nonaqueous phase liquid dissolution

A pore network simulation model was developed to investigate the fundamental physics of nonequilibrium nonaqueous phase liquid (NAPL) dissolution. The network model is a lattice of cubic chambers and rectangular tubes that represent pore bodies and pore throats, respectively. Experimental data obtained by Powers [1992] were used to develop and validate the model. To ensure the network model was representative of a real porous medium, the pore size distribution of the network was calibrated by matching simulated and experimental drainage and imbibition capillary pressure‐saturation curves. The predicted network residual styrene blob‐size distribution was nearly identical to the observed distribution. The network model reproduced the observed hydraulic conductivity and produced relative permeability curves that were representative of a poorly consolidated sand. Aqueous‐phase transport was represented by applying the equation for solute flux to the network tubes and solving for solute concentrations in the network chambers. Complete mixing was found to be an appropriate approximation for calculation of chamber concentrations. Mass transfer from NAPL blobs was represented using a corner diffusion model. Predicted results of solute concentration versus Peclet number and of modified Sherwood number versus Peclet number for the network model compare favorably with experimental data for the case in which NAPL blob dissolution was negligible. Predicted results of normalized effluent concentration versus pore volume for the network were similar to the experimental data for the case in which NAPL blob dissolution occurred with time.

Water Resources Research↗
USGS2000

Influence of net freshwater supply on salinity in Florida Bay

An annual water budget for Florida Bay, the large, seasonally hypersaline estuary in the Everglades National Park, was constructed using physically based models and long‐term (31 years) data on salinity, hydrology, and climate. Effects of seasonal and interannual variations of the net freshwater supply (runoff plus rainfall minus evaporation) on salinity variation within the bay were also examined. Particular attention was paid to the effects of runoff, which are the focus of ambitious plans to restore and conserve the Florida Bay ecosystem. From 1965 to 1995 the annual runoff from the Everglades into the bay was less than one tenth of the annual direct rainfall onto the bay, while estimated annual evaporation slightly exceeded annual rainfall. The average net freshwater supply to the bay over a year was thus approximately zero, and interannual variations in salinity appeared to be affected primarily by interannual fluctuations in rainfall. At the annual scale, runoff apparently had little effect on the bay as a whole during this period. On a seasonal basis, variations in rainfall, evaporation, and runoff were not in phase, and the net freshwater supply to the bay varied between positive and negative values, contributing to a strong seasonal pattern in salinity, especially in regions of the bay relatively isolated from exchanges with the Gulf of Mexico and Atlantic Ocean. Changes in runoff could have a greater effect on salinity in the bay if the seasonal patterns of rainfall and evaporation and the timing of the runoff are considered. One model was also used to simulate spatial and temporal patterns of salinity responses expected to result from changes in net freshwater supply. Simulations in which runoff was increased by a factor of 2 (but with no change in spatial pattern) indicated that increased runoff will lower salinity values in eastern Florida Bay, increase the variability of salinity in the South Region, but have little effect on salinity in the Central and West Regions.

Florida↗
USGS2001

Solute transport along preferential flow paths in unsaturated fractures

Laboratory experiments were conducted to study solute transport along preferential flow paths in unsaturated, inclined fractures. Qualitative aspects of solute transport were identified in a miscible dye tracer experiment conducted in a transparent replica of a natural granite fracture. Additional experiments were conducted to measure the breakthrough curves of a conservative tracer introduced into an established preferential flow path in two different fracture replicas and a rock‐replica combination. The influence of gravity was investigated by varying fracture inclination. The relationship between the travel times of the solute and the relative influence of gravity was substantially affected by two modes of intermittent flow that occurred: the snapping rivulet and the pulsating blob modes. The measured travel times of the solute were evaluated with three transfer function models: the axial dispersion, the reactors‐in‐series, and the lognormal models. The three models described the solute travel times nearly equally well. A mechanistic model was also formulated to describe transport when the pulsating blob mode occurred which assumed blobs of water containing solute mixed with residual pools of water along the flow path.

Water Resources Research↗
USGS2002

Using groundwater temperature data to constrain parameter estimation in a groundwater flow model of a wetland system

Parameter estimation is a powerful way to calibrate models. While head data alone are often insufficient to estimate unique parameters due to model nonuniqueness, flow‐and‐heat‐transport modeling can constrain estimation and allow simultaneous estimation of boundary fluxes and hydraulic conductivity. In this work, synthetic and field models that did not converge when head data were used did converge when head and temperature were used. Furthermore, frequency domain analyses of head and temperature data allowed selection of appropriate modeling timescales. Inflows in the Wilton, Wisconsin, wetlands could be estimated over periods such as a growing season and over periods of a few days when heads were nearly steady and groundwater temperature varied during the day. While this methodology is computationally more demanding than traditional head calibration, the results gained are unobtainable using the traditional approach. These results suggest that temperature can efficiently supplement head data in systems where accurate flux calibration targets are unavailable.

Water Resources Research↗
USGS2000

An analysis of the pilot point methodology for automated calibration of an ensemble of conditionally simulated transmissivity fields

An analysis of the pilot point method for automated calibration of an ensemble of conditionally simulated transmissivity fields was conducted on the basis of the simplifying assumption that the flow model is a linear function of log transmissivity. The analysis shows that the pilot point and conditional simulation method of model calibration and uncertainty analysis can produce accurate uncertainty measures if it can be assumed that errors of unknown origin in the differences between observed and model‐computed water pressures are small. When this assumption is not met, the method could yield significant errors from overparameterization and the neglect of potential sources of model inaccuracy. The conditional simulation part of the method is also shown to be a variant of the percentile bootstrap method, so that when applied to a nonlinear model, the method is subject to bootstrap errors. These sources of error must be considered when using the method.

Water Resources Research↗
USGS2000

Estimating the variance and integral scale of the transmissivity field using head residual increments

A modification of previously published solutions regarding the spatial variation of hydraulic heads is discussed whereby the semivariogram of increments of head residuals (termed head residual increments HRIs) are related to the variance and integral scale of the transmissivity field. A first‐order solution is developed for the case of a transmissivity field which is isotropic and whose second‐order behavior can be characterized by an exponential covariance structure. The estimates of the variance σ Y 2 and the integral scale λ of the log transmissivity field are then obtained via fitting a theoretical semivariogram for the HRI to its sample semivariogram. This approach is applied to head data sampled from a series of two‐dimensional, simulated aquifers with isotropic, exponential covariance structures and varying degrees of heterogeneity (σ Y 2 = 0.25, 0.5, 1.0, 2.0, and 5.0). The results show that this method provided reliable estimates for both λ and σ Y 2 in aquifers with the value of σ Y 2 up to 2.0, but the errors in those estimates were higher for σ Y 2 equal to 5.0. It is also demonstrated through numerical experiments and theoretical arguments that the head residual increments will provide a sample semivariogram with a lower variance than will the use of the head residuals without calculation of increments.

Water Resources Research↗
USGS2000

Mechanism and rate of denitrification in an agricultural watershed: Electron and mass balance along groundwater flow paths

The rate and mechanism of nitrate removal along and between groundwater flow paths were investigated using a series of well nests screened in an unconfined sand and gravel aquifer. Intensive agricultural activity in this area has resulted in nitrate concentrations in groundwater often exceeding drinking water standards. Both the extent and rate of denitrification varied depending on the groundwater flow path. While little or no denitrification occurred in much of the upland portions of the aquifer, a gradual redox gradient is observed as aerobic upland groundwater moves deeper in the aquifer. In contrast, a sharp shallow redox gradient is observed adjacent to a third‐order stream as aerobic groundwater enters reduced sediments. An essentially complete loss of nitrate concurrent with increases in excess N 2 provide evidence that denitrification occurs as groundwater enters this zone. Electron and mass balance calculations suggest that iron sulfide (e.g., pyrite) oxidation is the primary source of electrons for denitrification. Denitrification rate estimates were based on mass balance calculations using nitrate and excess N 2 coupled with groundwater travel times. Travel times were determined using a groundwater flow model and were constrained by chlorofluorocarbon‐based age dates. Denitrification rates were found to vary considerably between the two areas where denitrification occurs. Denitrification rates in the deep, upland portions of the aquifer were found to range from <0.01 to 0.14 m M of N per year; rates at the redoxcline along the shallow flow path range from 1.0 to 2.7 m M of N per year. Potential denitrification rates in groundwater adjacent to the stream may be much faster, with rates up to 140 m M per year based on an in situ experiment conducted in this zone.

Water Resources Research↗
USGS2000

Characterizing multiple timescales of stream and storage zone interaction that affect solute fate and transport in streams

The fate of contaminants in streams and rivers is affected by exchange and biogeochemical transformation in slowly moving or stagnant flow zones that interact with rapid flow in the main channel. In a typical stream, there are multiple types of slowly moving flow zones in which exchange and transformation occur, such as stagnant or recirculating surface water as well as subsurface hyporheic zones. However, most investigators use transport models with just a single storage zone in their modeling studies, which assumes that the effects of multiple storage zones can be lumped together. Our study addressed the following question: Can a single‐storage zone model reliably characterize the effects of physical retention and biogeochemical reactions in multiple storage zones? We extended an existing stream transport model with a single storage zone to include a second storage zone. With the extended model we generated 500 data sets representing transport of nonreactive and reactive solutes in stream systems that have two different types of storage zones with variable hydrologic conditions. The one storage zone model was tested by optimizing the lumped storage parameters to achieve a best fit for each of the generated data sets. Multiple storage processes were categorized as possessing I, additive; II, competitive; or III, dominant storage zone characteristics. The classification was based on the goodness of fit of generated data sets, the degree of similarity in mean retention time of the two storage zones, and the relative distributions of exchange flux and storage capacity between the two storage zones. For most cases (>90%) the one storage zone model described either the effect of the sum of multiple storage processes (category I) or the dominant storage process (category III). Failure of the one storage zone model occurred mainly for category II, that is, when one of the storage zones had a much longer mean retention time ( t s ratio > 5.0) and when the dominance of storage capacity and exchange flux occurred in different storage zones. We also used the one storage zone model to estimate a “single” lumped rate constant representing the net removal of a solute by biogeochemical reactions in multiple storage zones. For most cases the lumped rate constant that was optimized by one storage zone modeling estimated the flux‐weighted rate constant for multiple storage zones. Our results explain how the relative hydrologic properties of multiple storage zones (retention time, storage capacity, exchange flux, and biogeochemical reaction rate constant) affect the reliability of lumped parameters determined by a one storage zone transport model. We conclude that stream transport models with a single storage compartment will in most cases reliably characterize the dominant physical processes of solute retention and biogeochemical reactions in streams with multiple storage zones.

Water Resources Research↗
USGS2000

Identifying fracture‐zone geometry using simulated annealing and hydraulic‐connection data

A new approach is presented to condition geostatistical simulation of high‐permeability zones in fractured rock to hydraulic‐connection data. A simulated‐annealing algorithm generates three‐dimensional (3‐D) realizations conditioned to borehole data, inferred hydraulic connections between packer‐isolated borehole intervals, and an indicator (fracture zone or background‐ K bedrock) variogram model of spatial variability. We apply the method to data from the U.S. Geological Survey Mirror Lake Site in New Hampshire, where connected high‐permeability fracture zones exert a strong control on fluid flow at the hundred‐meter scale. Single‐well hydraulic‐packer tests indicate where permeable fracture zones intersect boreholes, and multiple‐well pumping tests indicate the degree of hydraulic connection between boreholes. Borehole intervals connected by a fracture zone exhibit similar hydraulic responses, whereas intervals not connected by a fracture zone exhibit different responses. Our approach yields valuable insights into the 3‐D geometry of fracture zones at Mirror Lake. Statistical analysis of the realizations yields maps of the probabilities of intersecting specific fracture zones with additional wells. Inverse flow modeling based on the assumption of equivalent porous media is used to estimate hydraulic conductivity and specific storage and to identify those fracture‐zone geometries that are consistent with hydraulic test data.

Water Resources Research↗
USGS2000

Tracer transport in fractured crystalline rock: Evidence of nondiffusive breakthrough tailing

Extended tailing of tracer breakthrough is often observed in pulse injection tracer tests conducted in fractured geologic media. This behavior has been attributed to diffusive exchange of tracer between mobile fluids traveling through channels in fractures and relatively stagnant fluid between fluid channels, along fracture walls, or within the bulk matrix. We present a field example where tracer breakthrough tailing apparently results from nondiffusive transport. Tracer tests were conducted in a fractured crystalline rock using both a convergent and weak dipole injection and pumping scheme. Deuterated water, bromide, and pentafluorobenzoic acid were selected as tracers for their wide range in molecular diffusivity. The late time behavior of the normalized breakthrough curves were consistent for all tracers, even when the pumping rate was changed. The lack of separation between tracers of varying diffusivity indicates that strong breakthrough tailing in fractured geologic media may be caused by advective transport processes. This finding has implications for the interpretation of tracer tests designed to measure matrix diffusion in situ and the prediction of contaminant transport in fractured rock.

Water Resources Research↗
USGS2000

Landslide triggering by rain infiltration

Landsliding in response to rainfall involves physical processes that operate on disparate timescales. Relationships between these timescales guide development of a mathematical model that uses reduced forms of Richards equation to evaluate effects of rainfall infiltration on landslide occurrence, timing, depth, and acceleration in diverse situations. The longest pertinent timescale is A / D 0 , where D 0 is the maximum hydraulic diffusivity of the soil and A is the catchment area that potentially affects groundwater pressures at a prospective landslide slip surface location with areal coordinates x , y and depth H . Times greater than A / D 0 are necessary for establishment of steady background water pressures that develop at ( x , y , H ) in response to rainfall averaged over periods that commonly range from days to many decades. These steady groundwater pressures influence the propensity for landsliding at ( x , y , H ), but they do not trigger slope failure. Failure results from rainfall over a typically shorter timescale H 2 / D 0 associated with transient pore pressure transmission during and following storms. Commonly, this timescale ranges from minutes to months. The shortest timescale affecting landslide responses to rainfall is √(H/g), where g is the magnitude of gravitational acceleration. Postfailure landslide motion occurs on this timescale, which indicates that the thinnest landslides accelerate most quickly if all other factors are constant. Effects of hydrologic processes on landslide processes across these diverse timescales are encapsulated by a response function, R(t*) = √(t*/π) exp (-1/t*) - erfc (1/√t*), which depends only on normalized time, t*. Use of R(t*) in conjunction with topographic data, rainfall intensity and duration information, an infinite-slope failure criterion, and Newton's second law predicts the timing, depth, and acceleration of rainfall-triggered landslides. Data from contrasting landslides that exhibit rapid, shallow motion and slow, deep-seated motion corroborate these predictions.

Water Resources Research↗