Search USGS⌕ Search

SEARCH · Search USGS

Results for “Algorithms”

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 847 records · Page 47Linked to original sources
USGS2002

A statistical evaluation of non-ergodic variogram estimators

Geostatistics is a set of statistical techniques that is increasingly used to characterize spatial dependence in spatially referenced ecological data. A common feature of geostatistics is predicting values at unsampled locations from nearby samples using the kriging algorithm. Modeling spatial dependence in sampled data is necessary before kriging and is usually accomplished with the variogram and its traditional estimator. Other types of estimators, known as non-ergodic estimators, have been used in ecological applications. Non-ergodic estimators were originally suggested as a method of choice when sampled data are preferentially located and exhibit a skewed frequency distribution. Preferentially located samples can occur, for example, when areas with high values are sampled more intensely than other areas. In earlier studies the visual appearance of variograms from traditional and non-ergodic estimators were compared. Here we evaluate the estimators' relative performance in prediction. We also show algebraically that a non-ergodic version of the variogram is equivalent to the traditional variogram estimator. Simulations, designed to investigate the effects of data skewness and preferential sampling on variogram estimation and kriging, showed the traditional variogram estimator outperforms the non-ergodic estimators under these conditions. We also analyzed data on carabid beetle abundance, which exhibited large-scale spatial variability (trend) and a skewed frequency distribution. Detrending data followed by robust estimation of the residual variogram is demonstrated to be a successful alternative to the non-ergodic approach.

Environmental and Ecological Statistics↗
USGS2003

Preliminary investigation of submerged aquatic vegetation mapping using hyperspectral remote sensing

The use of airborne hyperspectral remote sensing imagery for automated mapping of submerged aquatic vegetation (SAV) in the tidal Potomac River was investigated for near to realtime resource assessment and monitoring. Airborne hyperspectral imagery and field spectrometer measurements were obtained in October of 2000. A spectral library database containing selected ground-based and airborne sensor spectra was developed for use in image processing. The spectral library is used to automate the processing of hyperspectral imagery for potential real-time material identification and mapping. Field based spectra were compared to the airborne imagery using the database to identify and map two species of SAV (Myriophyllum spicatum and Vallisneria americana). Overall accuracy of the vegetation maps derived from hyperspectral imagery was determined by comparison to a product that combined aerial photography and field based sampling at the end of the SAV growing season. The algorithms and databases developed in this study will be useful with the current and forthcoming space-based hyperspectral remote sensing systems.

Environmental Monitoring and Assessment↗
USGS2003

A hydrologic network supporting spatially referenced regression modeling in the Chesapeake Bay watershed

The U.S. Geological Survey has developed a methodology for statistically relating nutrient sources and land-surface characteristics to nutrient loads of streams. The methodology is referred to as SPAtially Referenced Regressions On Watershed attributes (SPARROW), and relates measured stream nutrient loads to nutrient sources using nonlinear statistical regression models. A spatially detailed digital hydrologic network of stream reaches, stream-reach characteristics such as mean streamflow, water velocity, reach length, and travel time, and their associated watersheds supports the regression models. This network serves as the primary framework for spatially referencing potential nutrient source information such as atmospheric deposition, septic systems, point-sources, land use, land cover, and agricultural sources and land-surface characteristics such as land use, land cover, average-annual precipitation and temperature, slope, and soil permeability. In the Chesapeake Bay watershed that covers parts of Delaware, Maryland, Pennsylvania, New York, Virginia, West Virginia, and Washington D.C., SPARROW was used to generate models estimating loads of total nitrogen and total phosphorus representing 1987 and 1992 land-surface conditions. The 1987 models used a hydrologic network derived from an enhanced version of the U.S. Environmental Protection Agency's digital River Reach File, and course resolution Digital Elevation Models (DEMs). A new hydrologic network was created to support the 1992 models by generating stream reaches representing surface-water pathways defined by flow direction and flow accumulation algorithms from higher resolution DEMs. On a reach-by-reach basis, stream reach characteristics essential to the modeling were transferred to the newly generated pathways or reaches from the enhanced River Reach File used to support the 1987 models. To complete the new network, watersheds for each reach were generated using the direction of surface-water flow derived from the DEMs. This network improves upon existing digital stream data by increasing the level of spatial detail and providing consistency between the reach locations and topography. The hydrologic network also aids in illustrating the spatial patterns of predicted nutrient loads and sources contributed locally to each stream, and the percentages of nutrient load that reach Chesapeake Bay.

Environmental Monitoring and Assessment↗
USGS2004

Sampling issues affecting accuracy of likelihood-based classification using genetical data

We demonstrate the effectiveness of a genetic algorithm for discovering multi-locus combinations that provide accurate individual assignment decisions and estimates of mixture composition based on likelihood classification. Using simulated data representing different levels of inter-population differentiation (F st ~ 0.01 and 0.10), genetic diversities (four or eight alleles per locus), and population sizes (20, 40, 100 individuals in baseline populations), we show that subsets of loci can be identified that provide comparable levels of accuracy in classification decisions relative to entire multi-locus data sets, where 5, 10, or 20 loci were considered. Microsatellite data sets from hatchery strains of lake trout, Salvelinus namaycush, representing a comparable range of inter-population levels of differentiation in allele frequencies confirmed simulation results. For both simulated and empirical data sets, assignment accuracy was achieved using fewer loci (e.g., three or four loci out of eight for empirical lake trout studies). Simulation results were used to investigate properties of the 'leave-one-out' (L1O) method for estimating assignment error rates. Accuracy of population assignments based on L1O methods should be viewed with caution under certain conditions, particularly when baseline population sample sizes are low (<50).

Environmental Biology of Fishes↗
USGS2011

Workshop discusses community models for coastal sediment transport

Numerical models of coastal sediment transport are increasingly used to address problems ranging from remediation of contaminated sediments, to siting of sewage outfalls and disposal sites, to evaluating impacts of coastal development. They are also used as a test bed for sediment-transport algorithms, to provide realistic settings for biological and geochemical models, and for a variety of other research, both fundamental and applied. However, there are few full-featured, publicly available coastal sediment-transport models, and fewer still that are well tested and have been widely applied. This was the motivation for a workshop in Woods Hole, Massachusetts, on June 22–23, 2000, that explored the establishment of community models for coastal sediment-transport processes.

Eos Science News↗
USGS1999

Dynamic 3D simulations of earthquakes on en echelon faults

One of the mysteries of earthquake mechanics is why earthquakes stop. This process determines the difference between small and devastating ruptures. One possibility is that fault geometry controls earthquake size. We test this hypothesis using a numerical algorithm that simulates spontaneous rupture propagation in a three-dimensional medium and apply our knowledge to two California fault zones. We find that the size difference between the 1934 and 1966 Parkfield, California, earthquakes may be the product of a stepover at the southern end of the 1934 earthquake and show how the 1992 Landers, California, earthquake followed physically reasonable expectations when it jumped across en echelon faults to become a large event. If there are no linking structures, such as transfer faults, then strike-slip earthquakes are unlikely to propagate through stepover s >5 km wide.

Geophysical Research Letters↗
USGS2000

The crustal thickness of Australia

We investigate the crustal structure of the Australian continent using the temporary broadband stations of the Skippy and Kimba projects and permanent broadband stations. We isolate near-receiver information, in the form of crustal P-to-S conversions, using the receiver function technique. Stacked receiver functions are inverted for S velocity structure using a Genetic Algorithm approach to Receiver Function Inversion (GARFI). From the resulting velocity models we are able to determine the Moho depth and to classify the width of the crust-mantle transition for 65 broadband stations. Using these results and 51 independent estimates of crustal thickness from refraction and reflection profiles, we present a new, improved, map of Moho depth for the Australian continent. The thinnest crust (25 km) occurs in the Archean Yilgarn Craton in Western Australia; the thickest crust (61 km) occurs in Proterozoic central Australia. The average crustal thickness is 38.8 km (standard deviation 6.2 km). Interpolation error estimates are made using kriging and fall into the range 2.5–7.0 km. We find generally good agreement between the depth to the seismologically defined Moho and xenolith-derived estimates of crustal thickness beneath northeastern Australia. However, beneath the Lachlan Fold Belt the estimates are not in agreement, and it is possible that the two techniques are mapping differing parts of a broad Moho transition zone. The Archean cratons of Western Australia appear to have remained largely stable since cratonization, reflected in only slight variation of Moho depth. The largely Proterozoic center of Australia shows relatively thicker crust overall as well as major Moho offsets. We see evidence of the margin of the contact between the Precambrian craton and the Tasman Orogen, referred to as the Tasman Line.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Fault structure and mechanics of the Hayward Fault, California from double-difference earthquake locations

The relationship between small-magnitude seismicity and large-scale crustal faulting along the Hayward Fault, California, is investigated using a double-difference (DD) earthquake location algorithm. We used the DD method to determine high-resolution hypocenter locations of the seismicity that occurred between 1967 and 1998. The DD technique incorporates catalog travel time data and relative P and S wave arrival time measurements from waveform cross correlation to solve for the hypocentral separation between events. The relocated seismicity reveals a narrow, near-vertical fault zone at most locations. This zone follows the Hayward Fault along its northern half and then diverges from it to the east near San Leandro, forming the Mission trend. The relocated seismicity is consistent with the idea that slip from the Calaveras Fault is transferred over the Mission trend onto the northern Hayward Fault. The Mission trend is not clearly associated with any mapped active fault as it continues to the south and joins the Calaveras Fault at Calaveras Reservoir. In some locations, discrete structures adjacent to the main trace are seen, features that were previously hidden in the uncertainty of the network locations. The fine structure of the seismicity suggests that the fault surface on the northern Hayward Fault is curved or that the events occur on several substructures. Near San Leandro, where the more westerly striking trend of the Mission seismicity intersects with the surface trace of the (aseismic) southern Hayward Fault, the seismicity remains diffuse after relocation, with strong variation in focal mechanisms between adjacent events indicating a highly fractured zone of deformation. The seismicity is highly organized in space, especially on the northern Hayward Fault, where it forms horizontal, slip-parallel streaks of hypocenters of only a few tens of meters width, bounded by areas almost absent of seismic activity. During the interval from 1984 to 1998, when digital waveforms are available, we find that fewer than 6.5% of the earthquakes can be classified as repeating earthquakes, events that rupture the same fault patch more than one time. These most commonly are located in the shallow creeping part of the fault, or within the streaks at greater depth. The slow repeat rate of 2–3 times within the 15-year observation period for events with magnitudes around M = 1.5 is indicative of a low slip rate or a high stress drop. The absence of microearthquakes over large, contiguous areas of the northern Hayward Fault plane in the depth interval from ∼5 to 10 km and the concentrations of seismicity at these depths suggest that the aseismic regions are either locked or retarded and are storing strain energy for release in future large-magnitude earthquakes.

California↗
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

Regional interdisciplinary paleoflood approach to assess extreme flood potential

In the past decade, there has been a growing interest of dam safety officials to incorporate a risk‐based analysis for design‐flood hydrology. Extreme or rare floods, with probabilities in the range of about 10 −3 to 10 −7 chance of occurrence per year, are of continuing interest to the hydrologic and engineering communities for purposes of planning and design of structures such as dams [ National Research Council , 1988]. The National Research Council stresses that as much information as possible about floods needs to be used for evaluation of the risk and consequences of any decision. A regional interdisciplinary paleoflood approach was developed to assist dam safety officials and floodplain managers in their assessments of the risk of large floods. The interdisciplinary components included documenting maximum paleofloods and a regional analyses of contemporary extreme rainfall and flood data to complement a site‐specific probable maximum precipitation study [ Tomlinson and Solak , 1997]. The cost‐effective approach, which can be used in many other hydrometeorologic settings, was applied to Elkhead Reservoir in Elkhead Creek (531 km 2 ) in northwestern Colorado; the regional study area was 10,900 km 2 . Paleoflood data using bouldery flood deposits and noninundation surfaces for 88 streams were used to document maximum flood discharges that have occurred during the Holocene. Several relative dating methods were used to determine the age of paleoflood deposits and noninundation surfaces. No evidence of substantial flooding was found in the study area. The maximum paleoflood of 135 m 3 s −1 for Elkhead Creek is about 13% of the site‐specific probable maximum flood of 1020 m 3 s −1 . Flood‐frequency relations using the expected moments algorithm, which better incorporates paleoflood data, were developed to assess the risk of extreme floods. Envelope curves encompassing maximum rainfall (181 sites) and floods (218 sites) were developed for northwestern Colorado to help define maximum contemporary and Holocene flooding in Elkhead Creek and in a regional frequency context. Study results for Elkhead Reservoir were accepted by the Colorado State Engineer for dam safety certification.

Water Resources Research↗
USGS2002

Fault structure and kinematics of the Long Valley Caldera region, California, revealed by high-accuracy earthquake hypocenters and focal mechanism stress inversions

We have determined high-resolution hypocenters for 45,000+ earthquakes that occurred between 1980 and 2000 in the Long Valley caldera area using a double-difference earthquake location algorithm and routinely determined arrival times. The locations reveal numerous discrete fault planes in the southern caldera and adjacent Sierra Nevada block (SNB). Intracaldera faults include a series of east/west-striking right-lateral strike-slip faults beneath the caldera's south moat and a series of more northerly striking strike-slip/normal faults beneath the caldera's resurgent dome. Seismicity in the SNB south of the caldera is confined to a crustal block bounded on the west by an east-dipping oblique normal fault and on the east by the Hilton Creek fault. Two NE-striking left-lateral strike-slip faults are responsible for most seismicity within this block. To understand better the stresses driving seismicity, we performed stress inversions using focal mechanisms with 50 or more first motions. This analysis reveals that the least principal stress direction systematically rotates across the studied region, from NE to SW in the caldera's south moat to WNW-ESE in Round Valley, 25 km to the SE. Because WNW-ESE extension is characteristic of the western boundary of the Basin and Range province, caldera area stresses appear to be locally perturbed. This stress perturbation does not seem to result from magma chamber inflation but may be related to the significant (???20 km) left step in the locus of extension along the Sierra Nevada/Basin and Range province boundary. This implies that regional-scale tectonic processes are driving seismic deformation in the Long Valley caldera.

California↗
USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Time‐lapse imaging of saline‐tracer transport in fractured rock using difference‐attenuation radar tomography

Accurate characterization of fractured‐rock aquifer heterogeneity remains one of the most challenging and important problems in groundwater hydrology. We demonstrate a promising strategy to identify preferential flow paths in fractured rock using a combination of geophysical monitoring and conventional hydrogeologic tests. Cross‐well difference‐attenuation ground‐penetrating radar was used to monitor saline‐tracer migration in an experiment at the U.S. Geological Survey Fractured Rock Hydrology Research Site in Grafton County, New Hampshire. Radar data sets were collected every 10 min in three adjoining planes for 5 hours during each of 12 tracer tests. An innovative inversion method accounts for data acquisition times and temporal changes in attenuation during data collection. The inverse algorithm minimizes a combination of two functions. The first is the sum of weighted squared data residuals. Second is a measure of solution complexity based on an a priori space‐time covariance function, subject to constraints that limit radar‐attenuation changes to regions of the tomograms traversed by high difference‐attenuation ray paths. The time series of tomograms indicate relative tracer concentrations and tracer arrival times in the image planes; from these we infer the presence and location of a preferential flow path within a previously identified zone of transmissive fractures. These results provide new insights into solute channeling and the nature of aquifer heterogeneity at the site.

New Hampshire↗
USGS2003

Comparisons of two moments‐based estimators that utilize historical and paleoflood data for the log Pearson type III distribution

The expected moments algorithm (EMA) [ Cohn et al. , 1997 ] and the Bulletin 17B [ Interagency Committee on Water Data , 1982 ] historical weighting procedure (B17H) for the log Pearson type III distribution are compared by Monte Carlo computer simulation for cases in which historical and/or paleoflood data are available. The relative performance of the estimators was explored for three cases: fixed‐threshold exceedances, a fixed number of large floods, and floods generated from a different parent distribution. EMA can effectively incorporate four types of historical and paleoflood data: floods where the discharge is explicitly known, unknown discharges below a single threshold, floods with unknown discharge that exceed some level, and floods with discharges described in a range. The B17H estimator can utilize only the first two types of historical information. Including historical/paleoflood data in the simulation experiments significantly improved the quantile estimates in terms of mean square error and bias relative to using gage data alone. EMA performed significantly better than B17H in nearly all cases considered. B17H performed as well as EMA for estimating X 100 in some limited fixed‐threshold exceedance cases. EMA performed comparatively much better in other fixed‐threshold situations, for the single large flood case, and in cases when estimating extreme floods equal to or greater than X 500 . B17H did not fully utilize historical information when the historical period exceeded 200 years. Robustness studies using GEV‐simulated data confirmed that EMA performed better than B17H. Overall, EMA is preferred to B17H when historical and paleoflood data are available for flood frequency analysis.

Water Resources Research↗
USGS2003

An analytical formulation of two‐dimensional groundwater dispersion induced by surficial recharge variability

A predominant cause of dispersion in groundwater is advective mixing due to variability in seepage rates. Hydraulic conductivity variations have been extensively researched as a cause of this seepage variability. In this paper the effect of variations in surface recharge to a shallow surficial aquifer is investigated as an important additional effect. An analytical formulation has been developed that relates aquifer parameters and the statistics of recharge variability to increases in the dispersivity. This is accomplished by solving Fourier transforms of the small perturbation forms of the groundwater flow equations. Two field studies are presented in this paper to determine the statistics of recharge variability for input to the analytical formulation. A time series of water levels at a continuous groundwater recorder is used to investigate the temporal statistics of hydraulic head caused by recharge, and a series of infiltrometer measurements are used to define the spatial variability in the recharge parameters. With these field statistics representing head fluctuations due to recharge, the analytical formulation can be used to compute the dispersivity without an explicit representation of the recharge boundary. Results from a series of numerical experiments are used to define the limits of this analytical formulation and to provide some comparison. A sophisticated model has been developed using a particle‐tracking algorithm (modified to account for temporal variations) to estimate groundwater dispersion. Dispersivity increases of 9 percent are indicated by the analytical formulation for the aquifer at the field site. A comparison with numerical model results indicates that the analytical results are reasonable for shallow surficial aquifers in which two‐dimensional flow can be assumed.

Water Resources Research↗
USGS2003

High precision earthquake locations reveal seismogenic structure beneath Mammoth Mountain, California

In 1989, an unusual earthquake swarm occurred beneath Mammoth Mountain that was probably associated with magmatic intrusion. To improve our understanding of this swarm, we relocated Mammoth Mountain earthquakes using a double difference algorithm. Relocated hypocenters reveal that most earthquakes occurred on two structures, a near-vertical plane at 7–9 km depth that has been interpreted as an intruding dike, and a circular ring-like structure at ∼5.5 km depth, above the northern end of the inferred dike. Earthquakes on this newly discovered ring structure form a conical section that dips outward away from the aseismic interior. Fault-plane solutions indicate that in 1989 the seismicity ring was slipping as a ring-normal fault as the center of the mountain rose with respect to the surrounding crust. Seismicity migrated around the ring, away from the underlying dike at a rate of ∼0.4 km/month, suggesting that fluid movement triggered seismicity on the ring fault.

California↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 1. Theory and computation

To establish a theoretical basis for predicting and interpreting the behavior of rapid mass movements on Earth's surface, we develop and test a new computational model for gravity-driven motion of granular avalanches across irregular, three-dimensional (3-D) terrain. The principles embodied in the model are simple and few: continuum mass and momentum conservation and intergranular stress generation governed by Coulomb friction. However, significant challenges result from the necessity of satisfying these principles when deforming avalanches interact with steep and highly variable 3-D terrain. We address these challenges in four ways. (1) We formulate depth-averaged governing equations that are referenced to a rectangular Cartesian coordinate system (with z vertical) and that account explicitly for the effect of nonzero vertical accelerations on depth-averaged mass and momentum fluxes and stress states. (2) We compute fluxes of mass and momentum across vertical cell boundaries using a high-resolution finite volume method and Roe-type Riemann solver. Our algorithm incorporates flux difference splitting, an entropy correction for the flux, and eigenvector decomposition to embed the effects of driving and resisting forces in Riemann solutions. (3) We use a finite element method and avalanche displacements predicted by Riemann solutions to compute Coulomb stresses conjugate to the displacements in 3-D stress space. (4) We test the model output against analytical solutions, a sand cone conceptual experiment, and (in a companion paper) data from detailed laboratory experiments. Model results illustrate a complex interplay of basal traction and internal stress, and they successfully predict not only the gross behavior but also many details of avalanche motion from initiation to deposition.

Journal of Geophysical Research Earth Surface↗
USGS2004

Log Pearson type 3 quantile estimators with regional skew information and low outlier adjustments

The recently developed expected moments algorithm (EMA) [ Cohn et al. , 1997 ] does as well as maximum likelihood estimations at estimating log‐Pearson type 3 (LP3) flood quantiles using systematic and historical flood information. Needed extensions include use of a regional skewness estimator and its precision to be consistent with Bulletin 17B . Another issue addressed by Bulletin 17B is the treatment of low outliers. A Monte Carlo study compares the performance of Bulletin 17B using the entire sample with and without regional skew with estimators that use regional skew and censor low outliers, including an extended EMA estimator, the conditional probability adjustment (CPA) from Bulletin 17B , and an estimator that uses probability plot regression (PPR) to compute substitute values for low outliers. Estimators that neglect regional skew information do much worse than estimators that use an informative regional skewness estimator. For LP3 data the low outlier rejection procedure generally results in no loss of overall accuracy, and the differences between the MSEs of the estimators that used an informative regional skew are generally modest in the skewness range of real interest. Samples contaminated to model actual flood data demonstrate that estimators which give special treatment to low outliers significantly outperform estimators that make no such adjustment.

Water Resources Research↗