Search USGS⌕ Search

SEARCH · Search USGS

Results for “Techniques and Methods”

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 1,261 records · Page 70Linked to original sources
USGS2022

Evaluation of electrical and electromagnetic geophysical techniques to inspect earthen dam and levee structures in Arkansas

Within the state of Arkansas there is an increasing number of aging dams and levees that have little to no documentation concerning their construction or composition. Surface geophysical surveys offer a non-intrusive method for investigating these structures: To describe their lithologic makeup, to evaluate the materials that they were constructed upon, and to identify potential flow paths through them. Techniques such as electrical resistivity tomography, seismic refraction, and electromagnetic induction have all been used to image dams and levees and require additional information from geologic outcrops, geotechnical borings, or drill cores in order to make informed geologic interpretations of the geophysical models. These geologic models then allow the owners of these structures to make more informed decisions about their operation and maintenance. Between 2011 and 2018, the U.S. Geological Survey conducted geophysical and geotechnical investigations of three earthen structures within the state of Arkansas. Electrical and electromagnetic geophysical data were used to develop lithologic models of these structures and the underlying geology. Self-potential surveys were utilized to detect the movement of water through these structures indicating possible seepage pathways. Geotechnical methods such as electric and hydraulic direct-push well logs and cores acted as both a control on the geophysical interpretations and a confirmation of anomalies. This integrated approach detected the lack of an impermeable core within a levee, imaged a change in lithology of the bedrock forming the seal beneath a gravity dam, and identified a potential seepage feature within the core of an earthen dam. These results further support that this method of extending known lithologic features via surface and borehole geophysics is a useful approach for characterizing earthen water control structures.

Arkansas↗
USGS2011

A preliminary assessment of the spatial sources of contemporary suspended sediment in the Ohio River basin, United States, using water quality data from the NASQAN programme in a source tracing procedure

Reliable information on catchment scale suspended sediment sources is required to inform the design of management strategies for helping abate the numerous environmental issues associated with enhanced sediment mobilization and off‐site loadings. Since sediment fingerprinting techniques avoid many of the logistical constraints associated with using more traditional indirect measurement methods at catchment scale, such approaches have been increasingly reported in the international literature and typically use data sets collected specifically for sediment source apportionment purposes. There remains scope for investigating the potential for using geochemical data sets assembled by routine monitoring programmes to fingerprint sediment provenance. In the United States, routine water quality samples are collected as part of the US Geological Survey's revised National Stream Quality Accounting Network programme. Accordingly, the geochemistry data generated from these samples over a 10‐year period (1996–2006) were used as the basis for a fingerprinting exercise to assess the key tributary sub‐catchment spatial sources of contemporary suspended sediment transported by the Ohio River. Uncertainty associated with the spatial source estimates was quantified using a Monte Carlo approach in conjunction with mass balance modelling. Relative frequency weighted means were used as an alternative way of summarizing the spatial source contributions, thereby avoiding the need to use confidence limits. The results should be interpreted in the context of the routine, but infrequent nature, of the suspended sediment samples used to assemble geochemistry as a basis for the sourcing exercise. Nonetheless, the study demonstrates how routine monitoring samples can be used to provide some preliminary information on sediment provenance in large drainage basins.

Tennessee, Ohio, West Virginia, Pennsylvania, Indi↗
USGS1989

Artificial recharge of groundwater and its role in water management

This paper summarizes and discusses the various aspects and methods of artificial recharge with particular emphasis on its uses and potential role in water management in the Arabian Gulf region. Artificial recharge occurs when man's activities cause more water to enter an aquifer, either under pumping or non-pumping conditions, than otherwise would enter the aquifer. Use of artificial recharge can be a practical means of dealing with problems of overdraft of groundwater. Methods of artificial recharge may be grouped under two broad types: (a) water spreading techniques, and (b) well-injection techniques. Successful use of artificial recharge requires a thorough knowledge of the physical and chemical characteristics of the aquifier system, and extensive onsite experimentation and tailoring of the artificial-recharge technique to fit the local or areal conditions. In general, water spreading techniques are less expensive than well injection and large quantities of water can be handled. Water spreading can also result in significant improvement in quality of recharge waters during infiltration and movement through the unsaturated zone and the receiving aquifer. In comparison, well-injection techniques are often used for emplacement of fresh recharge water into saline aquifer zones to form a manageable lens of fresher water, which may later be partially withdrawn for use or continue to be maintained as a barrier against salt-water encroachment. A major advantage in use of groundwater is its availability, on demand to wells, from a natural storage reservoir that is relatively safe from pollution and from damage by sabotage or other hostile action. However, fresh groundwater occurs only in limited quantities in most of the Arabian Gulf region; also, it is heavily overdrafted in many areas, and receives very little natural recharge. Good use could be made of artificial recharge by well injection in replenishing and managing aquifers in strategic locations if sources of freshwater could be made available for the artificial-recharge operations.

Desalination↗
USGS1979

An algol program for dissimilarity analysis: a divisive-omnithetic clustering technique

Clustering techniques are used properly to generate hypotheses about patterns in data. Of the hierarchical techniques, those which are divisive and omnithetic possess many theoretically optimal properties. One such method, dissimilarity analysis, is implemented here in ALGOL 60, and determined to be competitive computationally with most other methods. ?? 1979.

Computers & Geosciences↗
USGS2007

Determining metal assimilation efficiency in aquatic invertebrates using enriched stable metal isotope tracers

We employ a novel approach that combines pulse-chase feeding and multi-labelled stable isotopes to determine gut passage time (GPT), gut retention time (GRT), food ingestion rate (IR) and assimilation efficiency (AE) of three trace elements for a freshwater gastropod. Lettuce isotopically enriched in 53 Cr, 65 Cu and 106 Cd was fed for 2 h to Lymnaea stagnalis . The release of tracers in feces and water was monitored for 48 h, during which unlabelled lettuce was provided ad libidum . The first defecation of 53 Cr occurred after 5 h of depuration (GPT), whereas 90% of the ingested 53 Cr was recovered in the feces after 22.5 h of depuration (GRT). 53 Chromium was not significantly accumulated in the soft tissues upon exposure. In contrast, 65 Cu and 106 Cd assimilation was detectable for most experimental snails, i.e., 65/63 Cu and 106/114 Cd ratios in exposed snails were higher than those for controls. Food IR during the labelled feeding phase was 0.16 ± 0.07 g g −1 d −1 . IR was inferred from the amount of 53 Cr egested in the feces during depuration and the concentration of 53 Cr in the labelled lettuce. Assimilation efficiencies (±95% CI) determined using mass balance calculations were 84 ± 4% for Cu and 85 ± 3% for Cd. The ratio method yields similar AE estimates. Expanding the application of this novel stable isotope tracer technique to other metals in a wide variety of species will provide unique opportunities to evaluate the interplay between digestive processes and dietary influx of metals. Understanding the biological processes that modulate dietborne metal uptake is crucial to assess the toxicity of dietborne metals.

Aquatic Toxicology↗
USGS2022

Evaluating the effect of expert elicitation techniques on population status assessment in the face of large uncertainty

Population projection models are important tools for conservation and management. They are often used for population status assessments, for threat analyses, and to predict the consequences of conservation actions. Although conservation decisions should be informed by science, critical decisions are often made with very little information to support decision-making. Conversely, postponing decisions until better information is available may reduce the benefit of a conservation decision. When empirical data are limited or lacking, expert elicitation can be used to supplement existing data and inform model parameter estimates. The use of rigorous techniques for expert elicitation that account for uncertainty can improve the quality of the expert elicited values and therefore the accuracy of the projection models. One recurring challenge for summarizing expert elicited values is how to aggregate them. Here, we illustrate a process for population status assessment using a combination of expert elicitation and data from the ecological literature. We discuss the importance of considering various aggregation techniques, and illustrate this process using matrix population models for the wood turtle ( Glyptemys insculpta ) to assist U.S. Fish and Wildlife Service decision-makers with their Species Status Assessment. We compare estimates of population growth using data from the ecological literature and four alternative aggregation techniques for the expert-elicited values. The estimate of population growth rate based on estimates from the literature (λ mean = 0.952, 95% CI: 0.87–1.01) could not be used to unequivocally reject the hypotheses of a rapidly declining population nor the hypothesis of a stable, or even slightly growing population, whereas our results for the expert-elicited estimates supported the hypothesis that the wood turtle population will decline over time. Our results showed that the aggregation techniques used had an impact on model estimates, suggesting that the choice of techniques should be carefully considered. We discuss the benefits and limitations associated with each method and their relevance to the population status assessment. We note a difference in the temporal scope or inference between the literature-based estimates that provided insights about historical changes, whereas the expert-based estimates were forward looking. Therefore, conducting an expert-elicitation in addition to using parameter estimates from the literature improved our understanding of our species of interest.

Journal of Environmental Management↗
USGS2013

Radiocarbon dating late Quaternary loess deposits using small terrestrial gastropod shells

Constraining the ages and mass accumulation rates of late Quaternary loess deposits is often difficult because of the paucity of organic material typically available for 14 C dating and the inherent limitations of luminescence techniques. Radiocarbon dating of small terrestrial gastropod shells may provide an alternative to these methods as fossil shells are common in loess and contain ∼12% carbon by weight. Terrestrial gastropod assemblages in loess have been used extensively to reconstruct past environmental conditions but have been largely ignored for dating purposes. Here, we present the results of a multi-faceted approach to understanding the potential for using small terrestrial gastropod shells to date loess deposits in North America. First, we compare highly resolved 14 C ages of well-preserved wood and gastropod shells (Succineidae) recovered from a Holocene loess section in Alaska. Radiocarbon ages derived from the shells are nearly identical to wood and plant macrofossil ages throughout the section, which suggests that the shells behaved as closed systems with respect to carbon for at least the last 10 ka (thousands of calibrated 14 C years before present). Second, we apply 14 C dating of gastropod shells to late Pleistocene loess deposits in the Great Plains using stratigraphy and independent chronologies for comparison. The new shell ages require less interpretation than humic acid radiocarbon ages that are commonly used in loess studies, provide additional stratigraphic coverage to previous dating efforts, and are in correct stratigraphic order more often than their luminescence counterparts. Third, we show that Succineidae shells recovered from historic loess in the Matanuska River Valley, Alaska captured the 20th century 14 C bomb spike, which suggests that the shells can be used to date late Holocene and historic-aged loess. Finally, results from Nebraska and western Iowa suggest that, similar to other materials, shell ages approaching ∼40 ka should be viewed with caution as they may reflect trace amounts of contamination. In sum, our results show that small terrestrial gastropod shells, especially from the Succineidae family, provide reliable ages for late Quaternary loess deposits in North America.

Quaternary Science Reviews↗
USGS2001

Partitioning sources of variation in vertebrate species richness

Aim: To explore biogeographic patterns of terrestrial vertebrates in Maine, USA using techniques that would describe local and spatial correlations with the environment. Location: Maine, USA. Methods: We delineated the ranges within Maine (86,156 km2) of 275 species using literature and expert review. Ranges were combined into species richness maps, and compared to geomorphology, climate, and woody plant distributions. Methods were adapted that compared richness of all vertebrate classes to each environmental correlate, rather than assessing a single explanatory theory. We partitioned variation in species richness into components using tree and multiple linear regression. Methods were used that allowed for useful comparisons between tree and linear regression results. For both methods we partitioned variation into broad-scale (spatially autocorrelated) and fine-scale (spatially uncorrelated) explained and unexplained components. By partitioning variance, and using both tree and linear regression in analyses, we explored the degree of variation in species richness for each vertebrate group that Could be explained by the relative contribution of each environmental variable. Results: In tree regression, climate variation explained richness better (92% of mean deviance explained for all species) than woody plant variation (87%) and geomorphology (86%). Reptiles were highly correlated with environmental variation (93%), followed by mammals, amphibians, and birds (each with 84-82% deviance explained). In multiple linear regression, climate was most closely associated with total vertebrate richness (78%), followed by woody plants (67%) and geomorphology (56%). Again, reptiles were closely correlated with the environment (95%), followed by mammals (73%), amphibians (63%) and birds (57%). Main conclusions: Comparing variation explained using tree and multiple linear regression quantified the importance of nonlinear relationships and local interactions between species richness and environmental variation, identifying the importance of linear relationships between reptiles and the environment, and nonlinear relationships between birds and woody plants, for example. Conservation planners should capture climatic variation in broad-scale designs; temperatures may shift during climate change, but the underlying correlations between the environment and species richness will presumably remain.

Journal of Biogeography↗
USGS1958

Mass control of insects: The effects on fish and wildlife

The mass control of insects carried on during the past ten years has made possible the economical suppression and, in a few instances, the near eradication of pest insect populations over widespread areas. These large operations, usually featuring the use of the airplane for applying insecticides quickly and cheaply, have proven to be the means of effectuating control programs which never could have been carried on through use of the techniques and equipment of early-day economic entomology. With the development of the methods employed in modern mass control came the introduction of new toxicants, well-adapted for use against insects because of their potency, low cost and wide availability. The present-day entomologist is thus armed with the means to carry on a control task of almost any magnitude.

Bulletin of the Entomological Society of America↗
USGS2014

Forest Ecosystem respiration estimated from eddy covariance and chamber measurements under high turbulence and substantial tree mortality from bark beetles

Eddy covariance nighttime fluxes are uncertain due to potential measurement biases. Many studies report eddy covariance nighttime flux lower than flux from extrapolated chamber measurements, despite corrections for low turbulence. We compared eddy covariance and chamber estimates of ecosystem respiration at the GLEES Ameriflux site over seven growing seasons under high turbulence (summer night mean friction velocity (u*) = 0.7 m s −1 ), during which bark beetles killed or infested 85% of the aboveground respiring biomass. Chamber-based estimates of ecosystem respiration during the growth season, developed from foliage, wood and soil CO 2 efflux measurements, declined 35% after 85% of the forest basal area had been killed or impaired by bark beetles (from 7.1 ±0.22 μmol m −2 s −1 in 2005 to 4.6 ±0.16 μmol m −2 s −1 in 2011). Soil efflux remained at ~3.3 μmol m −2 s −1 throughout the mortality, while the loss of live wood and foliage and their respiration drove the decline of the chamber estimate. Eddy covariance estimates of fluxes at night remained constant over the same period, ~3.0 μmol m −2 s −1 for both 2005 (intact forest) and 2011 (85% basal area killed or impaired). Eddy covariance fluxes were lower than chamber estimates of ecosystem respiration (60% lower in 2005, and 32% in 2011), but the mean night estimates from the two techniques were correlated within a year (r 2 from 0.18-0.60). The difference between the two techniques was not the result of inadequate turbulence, because the results were robust to a u* filter of > 0.7 m s −1 . The decline in the average seasonal difference between the two techniques was strongly correlated with overstory leaf area (r 2 =0.92). The discrepancy between methods of respiration estimation should be resolved to have confidence in ecosystem carbon flux estimates.

Wyoming↗
USGS2006

Application of continuous seismic reflection methods to hydrologic studies

Oil and gas exploration and engineering studies in water‐covered areas routinely use continuous seismic profiling techniques to obtain subsurface geologic information. Such profiling also can be used effectively in hydrologic studies to define the geologic framework of aquifer systems, to locate hydrologic boundaries, and in some places, to interpret the lithologic character of aquifers and confining beds. High‐resolution continuous seismic profiling, through the use of nonexplosive sound sources, can be used to produce continuous records that require little data processing before hydrogeologic interpretation. High‐resolution tuned transducer, minisparker, Uniboom, and small airgun systems operating from small boats in shallow water are capable of transmitting energy that can penetrate up to a hundred meters of earth materials. The resulting analog records of the reflected seismic signal closely resemble geologic sections. Surveys on the Housatonic River in Connecticut and Massachusetts, and on the Connecticut River, in water from 1 to 10 meters deep, have defined the bedrock surface beneath 60 meters of stratified drift. Seismic‐reflection profiling also was used to determine the extent and thickness of recent lake‐bottom deposits in two Connecticut lakes. Surveys along 90 kilometers of river channels in the Sarasota‐port Charlotte, Florida, area defined in detail the stratigrapahy and continuity of the shallow aquifers.

Connecticut, Massachusetts↗
USGS1999

Estimation of methanogen biomass via quantitation of coenzyme M

Determination of the role of methanogenic bacteria in an anaerobic ecosystem often requires quantitation of the organisms. Because of the extreme oxygen sensitivity of these organisms and the inherent limitations of cultural techniques, an accurate biomass value is very difficult to obtain. We standardized a simple method for estimating methanogen biomass in a variety of environmental matrices. In this procedure we used the thiol biomarker coenzyme M (CoM) (2-mercaptoethanesulfonic acid), which is known to be present in all methanogenic bacteria. A high-performance liquid chromatography-based method for detecting thiols in pore water (A. Vairavamurthy and M. Mopper, Anal. Chim. Acta 78:363–370, 1990) was modified in order to quantify CoM in pure cultures, sediments, and sewage water samples. The identity of the CoM derivative was verified by using liquid chromatography-mass spectroscopy. The assay was linear for CoM amounts ranging from 2 to 2,000 pmol, and the detection limit was 2 pmol of CoM/ml of sample. CoM was not adsorbed to sediments. The methanogens tested contained an average of 19.5 nmol of CoM/mg of protein and 0.39 ± 0.07 fmol of CoM/cell. Environmental samples contained an average of 0.41 ± 0.17 fmol/cell based on most-probable-number estimates. CoM was extracted by using 1% tri-( N )-butylphosphine in isopropanol. More than 90% of the CoM was recovered from pure cultures and environmental samples. We observed no interference from sediments in the CoM recovery process, and the method could be completed aerobically within 3 h. Freezing sediment samples resulted in 46 to 83% decreases in the amounts of detectable CoM, whereas freezing had no effect on the amounts of CoM determined in pure cultures. The method described here provides a quick and relatively simple way to estimate methanogenic biomass.

Applied and Environmental Microbiology↗
USGS1990

Seismic site effects and the spatial interpolation of earthquake seismograms: Results using aftershocks of the 1986 North Palm Springs, California, earthquake

We address the following two questions. Can a microearthquake's ground motions be modeled by incident P and S waves that excite a site transfer-function that is a smooth function of incidence angle? Given recorded ground motions from a set of earthquakes having known locations and mechanisms, can we derive such a site transfer-function and use it to obtain the ground motions that would result from an earthquake source occurring somewhere in the same volume but having a location and mechanism that are different from the recorded events? Although many factors will cause two distinct microearthquake sources to have different seismograms at a common station, in this paper we concentrate only upon the differences caused by source mechanisms, P- and S-wave travel-time variations and by variations in the site transfer-function. We specifically exclude the effects of waves scattered from heterogeneities in the geologic structure away from the seismic site. We express the site transfer-function as a sum of several terms having simple dependences upon incidence angle and azimuth. Each term is an independent function of time. Given a set of seismograms observed at the site, we solve a linear system of equations for the time dependences of each term. These time series may be used to calculate the seismograms that would have resulted from an earthquake having arbitrary mechanism and location. This step is an interpolation. We have applied this technique to seismograms after aftershocks of the 1986 North Palm Springs earthquake. Our interpolation technique works fairly well within the volume occupied by the recorded events, but the method is not very successful at providing accurate seismograms for sources located outside the aftershock volume. The primary causes of the inaccuracy are the inadequacy of our chosen angular functions to model the site response fully and the likely scattering of seismic waves by geological heterogeneities (in this case, the Banning and Mission Creek faults) near the seismic stations. Our methods could be used to determine the effects of single scattering from lateral heterogeneities in geologic structure.

California↗
USGS2013

Reconnaissance investigation of the placer gold deposits in the Zarkashan Area of Interest, Ghazni Province, Afghanistan

This study is a reconnaissance investigation of the placer gold deposits in the Zarkashan Area of Interest (AOI) in Ghazni Province, Afghanistan. Detailed investigations of the Zarkashan gold deposits were conducted by Soviet and Afghan geologists in the 1960s and 1970s, prior to the development of satellite-based remote-sensing platforms and new methods of geomorphic mapping. The purpose of this study was to integrate new mapping techniques with previously collected concentration and borehole sampling data and geomorphologic interpretations to reassess the placer gold deposits in the Zarkashan AOI. A methodology combining the collection and analysis of historical sampling data, digital database development, hydrologic analysis, and geomorphic modeling was used. The analysis led to the reinterpretation of four gold-bearing seams along the Zarkashan River, and the calculation of an estimated gold reserve of approximately 3,000 kilograms (kg). This estimate is approximately 1,500 kg greater than the Soviet estimate. The result differs in large part due to the reinterpretation of the seams based on a much lower cutoff grade of 100 mg/m 3 . Because cutoff grade is dependent in part on the price of gold, the sevenfold increase in the price of gold since the undertaking of the Soviet investigation warranted our re-evaluation of their 500 mg/m 3 cutoff grade.

Ghazni Province↗
USGS2013

Reconnaissance investigation of the alluvial gold deposits in the North Takhar Area of Interest, Takhar Province, Afghanistan

This study is a reconnaissance assessment of the alluvial gold deposits of the North Takhar Area of Interest (AOI) in Takhar Province, Afghanistan. Soviet and Afghan geologists collected data and calculated the gold deposit reserves in Takhar Province in the 1970s, prior to the development of satellite-based remote-sensing platforms and new methods of geomorphic mapping. The purpose of this study was to integrate new mapping techniques with previously collected borehole sampling and concentration sampling data and geomorphologic interpretations to reassess the alluvial gold placer deposits in the North Takhar AOI. Through a combination of historical borehole and cross-section data and digital terrain modeling, the Samti, Nooraba-Khasar-Anjir, and Kocha River placer deposits were reassessed. Resource estimates were calculated to be 20,927 kilograms (kg) for Samti, 7,626 kg for Nooraba-Khasar-Anjir, 160 kg for the mouth of the Kocha, 1,047 kg for the lower Kocha, 113 kg for the middle Kocha, and 168 kg for the upper Kocha. Previous resource estimates conducted by the Soviets for the Samti and Nooraba-Khasar-Anjir deposits estimated 30,062 kg and 802 kg of gold, respectively. This difference between the new estimates and previous estimates results from the higher resolution geomorphic model and the interpretation of areas outside of the initial work zone studied by Soviet and Afghan geologists.

Takhar Province↗
USGS2017

An evaluation of the efficacy of using environmental DNA (eDNA) to detect giant gartersnakes (Thamnophis gigas)

Detecting populations of rare or cryptic species is essential for their conservation. For species like giant gartersnakes ( Thamnophis gigas ), conventional survey methods can be expensive and inefficient. These sampling difficulties might be overcome by modern techniques that detect deoxyribonucleic acid (DNA) shed by organisms into the environment (eDNA). We evaluated the efficacy of detecting giant gartersnake eDNA in water samples from the laboratory and at locations with known giant gartersnake populations in the Sacramento Valley of California, and failed to detect giant gartersnake DNA in most laboratory and all field samples. Aspects of giant gartersnake biology—such as highly keratinized skin and spending extensive time in the terrestrial environment, as well as hot, sunny, and turbid conditions in wetlands and canals of the Sacramento Valley—likely contributed to low detection probabilities. Although detection of eDNA shows promise under many conditions, further development is needed before sampling for eDNA is a viable option for detecting giant gartersnake populations.

California↗
USGS2009

Topographic Change Detection at Select Archeological Sites in Grand Canyon National Park, Arizona, 2006-2007

Topographic change of archeological sites within the Colorado River corridor of Grand Canyon National Park (GCNP) is a subject of interest to National Park Service managers and other stakeholders in the Glen Canyon Dam Adaptive Management Program. Although long-term topographic change resulting from a variety of natural processes is typical in the Grand Canyon region, a continuing debate exists on whether and how controlled releases from Glen Canyon Dam, located immediately upstream of GCNP, are impacting rates of site erosion, artifact transport, and the preservation of archeological resources. Continued erosion of archeological sites threatens both the archeological resources and our future ability to study evidence of past cultural habitation. Understanding the causes and effects of archaeological site erosion requires a knowledge of several factors including the location and magnitude of the changes occurring in relation to archeological resources, the rate of the changes, and the relative contribution of several potential causes, including sediment depletion associated with managed flows from Glen Canyon Dam, site-specific weather patterns, visitor impacts, and long-term climate change. To obtain this information, highly accurate, spatially specific data are needed from sites undergoing change. Using terrestrial lidar data collection techniques and novel TIN- and GRID-based change-detection post-processing methods, we analyzed topographic data for nine archeological sites. The data were collected using three separate data collection efforts spanning 16 months (May 2006 to September 2007). Our results documented positive evidence of erosion, deposition, or both at six of the nine sites investigated during this time interval. In addition, we observed possible signs of change at two of the other sites. Erosion was concentrated in established gully drainages and averaged 12 cm to 17 cm in depth with maximum depths of 50 cm. Deposition was concentrated at specific locations outside of drainages but generally was spread over larger areas (tens to hundreds of square meters). Maximum depths of deposition averaged 12 cm to 15 cm and reached as much as 35 cm. Overall, we found that the spatial distribution and magnitudes of surface change are specific to each site and that a thorough understanding of the geomorphology, weather, and sand supply is requisite for a complete understanding of the data. Additional work in combining these results with site-specific weather, hydrology, and geomorphology data will assist in the development of working models for determining the causes of the documented topographic changes.

Scientific Investigations Report↗
USGS2014

Numerical simulation of groundwater flow in the Columbia Plateau Regional Aquifer System, Idaho, Oregon, and Washington

A three-dimensional numerical model of groundwater flow was constructed for the Columbia Plateau Regional Aquifer System (CPRAS), Idaho, Oregon, and Washington, to evaluate and test the conceptual model of the system and to evaluate groundwater availability. The model described in this report can be used as a tool by water-resource managers and other stakeholders to quantitatively evaluate proposed alternative management strategies and assess the long‑term availability of groundwater. The numerical simulation of groundwater flow in the CPRAS was completed with support from the Groundwater Resources Program of the U.S. Geological Survey Office of Groundwater. The model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model, MODFLOW-NWT. The model uses 3-kilometer (9,842.5 feet) grid cells that subdivide the model domain by 126 rows and 131 columns. Vertically, the model domain was subdivided into six geologic model units. From youngest to oldest, the units are the Overburden, the Saddle Mountains Basalt, the Mabton Interbed, the Wanapum Basalt, the Vantage Interbed, and the Grande Ronde Basalt. Natural recharge was estimated using gridded historical estimates of annual precipitation for the period 1895–2007. Pre-development recharge was estimated to be the average natural recharge for this period. Irrigation recharge and irrigation pumping were estimated using a remote-sensing based soil-water balance model for the period 1985–2007. Pre-1985 irrigation recharge and pumping were estimated using previously published compilation maps and the history of large-scale irrigation projects. Pumping estimates for municipal, industrial, rural, residential, and all other uses were estimated using reported values and census data. Pumping was assumed to be negligible prior to 1920. Two models were constructed to simulate groundwater flow in the CPRAS: a steady-state predevelopment model representing conditions before large-scale pumping and irrigation altered the system, and a transient model representing the period 1900–2007. Automated parameter-estimation techniques (steady-state predevelopment model) and traditional trial-and-error (transient model) methods were used for calibration. To calibrate the steady-state and transient models, 10,525 and 46,460 water level measurements, respectively, and 50 base-flow estimates were used. The steady-state model simulated the shape, slope, and trends of a potentiometric surface that was generally consistent with mapped water levels. For the transient model, the mean and median difference between simulated and measured hydraulic heads is -10 and 4 ft, respectively, with a standard deviation of 164 ft over a 5,648 ft range of measured heads. The residuals for the simulation period show that 52 percent of the simulated heads exceeded measured heads with a median residual value of 43 ft, and 48 percent were less than measured heads with a median residual value of -76 ft. The CPRAS model was constructed to derive components of the groundwater budget and help understand the interactions of stresses, such as recharge, groundwater pumping, and commingling wells on the groundwater and surface-water system. Through these applications, the model can be used to identify trends in groundwater storage and use, and quantify groundwater availability. The annual groundwater budgets showed several patterns of change over the simulation period. Groundwater pumping was negligible until the 1950s and began to increase significantly during the 1970s and 1980s. Recharge was highly variable due to the interannual variability of precipitation, but began to increase in the late 1940s due to the increase in surface-water irrigation projects. Groundwater contributions to streamflow (base flow) followed recharge closely. However, in areas of significant groundwater-level decline, base flow is reduced. Groundwater pumping had the greatest effect on water levels, followed by irrigation enhanced recharge. Commingling was a larger factor in structurally complex upland areas where hydraulic-head gradients are naturally high. Groundwater pumping has increased substantially over the past 40–50 years; this increase resulted in declining water levels at depth and decreased base flows over much of the study area. The effects of pumping are mitigated somewhat by the increase of surface-water irrigation, especially in the shallow Overburden unit, and commingling wells in some areas. During dry to average years, groundwater pumping causes a net loss of groundwater in storage and current condition (2000–2007) groundwater pumping exceeds recharge in all but the wettest of years.

Idaho, Oregon, Washington↗