Search USGSSearch

SEARCH · Search USGS

Results for “Surface Water Techniques”

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 397 records · Page 22Linked to original sources

Direct-current vertical electrical-resistivity soundings in the Lower Peninsula of Michigan

Ninety-three direct-current vertical electrical-resistivity soundings were conducted in the Lower Peninsula of Michigan from June through October 1987. These soundings were made to assist in mapping the depth to brine in areas where borehole resistivity logs and water-quality data are sparse or lacking. The Schlumberger array for placement of current and potential electrodes was used for each sounding. Vertical electrical-resistivity sounding field data, shifted and smoothed sounding data, and electric layers calculated using inverse modeling techniques are presented. Also included is a summary of the near-surface conditions and depths to conductors and resistors for each sounding location.

Michigan

Traveltime and reaeration characteristics for a reach of the Rio Grande, Albuquerque, New Mexico, October 1991

Traveltime characteristics were determined using stream-velocity data and tracer-dye data for a reach of the Rio Grande. Traveltimes determined by the stream-velocity method were virtually the same as those determined by the tracer-dye and tracer-gas technique. The mean velocity of the stream was 1.12 miles per hour at a flow of about 300 cubic feet per second. Reaeration characteristics were determined using a propane tracer gas and a tracer-dye (rhodamine WT). Reaeration coefficients were adjusted for water temperature and the effects of wind movement on the water surface. The peak method-adjusted reaeration-coefficient mean value for the reach was 7.0 per day and ranged from 4.6 to 8.3 per day. The area method-adjusted reaeration- coefficient mean value for the reach was 7.7 per day and ranged from 5.5 to 10.4 per day.

Water-Resources Investigations Report

The role of hydrates, competing chemical constituents, and surface composition on CLNO2 formation

Atomic chlorine (Cl • ) affects air quality and atmospheric oxidizing capacity. Nitryl chloride (ClNO 2 ) – a common Cl • source–forms when chloride-containing aerosols react with dinitrogen pentoxide (N 2 O 5 ). A recent study showed that saline lakebed (playa) dust is an inland source of particulate chloride (Cl – ) that generates high ClNO 2 . However, the underlying physiochemical factors responsible for observed yields are poorly understood. To elucidate these controlling factors, we utilized single particle and bulk techniques to determine the chemical composition and mineralogy of playa sediment and dust samples from the southwest United States. Single particle analysis shows trace highly hygroscopic magnesium and calcium Cl-containing minerals are present and likely facilitate ClNO 2 formation at low humidity. Single particle and mineralogical analysis detected playa sediment organic matter that hinders N 2 O 5 uptake as well as 10 Å-clay minerals (e.g., Illite) that compete with water and chloride for N 2 O 5 . Finally, we show that the composition of the aerosol surface, rather than the bulk, is critical in ClNO 2 formation. These findings underscore the importance of mixing state, competing reactions, and surface chemistry on N 2 O 5 uptake and ClNO 2 yield for playa dusts and, likely, other aerosol systems. Therefore, consideration of particle surface composition is necessary to improve ClNO 2 and air quality modeling.

Environmental Science Technology

Faulting and groundwater in a desert environment: constraining hydrogeology using time-domain electromagnetic data

Within the south-western Mojave Desert, the Joshua Basin Water District is considering applying imported water into infiltration ponds in the Joshua Tree groundwater sub-basin in an attempt to artificially recharge the underlying aquifer. Scarce subsurface hydrogeological data are available near the proposed recharge site; therefore, time-domain electromagnetic (TDEM) data were collected and analysed to characterize the subsurface. TDEM soundings were acquired to estimate the depth to water on either side of the Pinto Mountain Fault, a major east-west trending strike-slip fault that transects the proposed recharge site. While TDEM is a standard technique for groundwater investigations, special care must be taken when acquiring and interpreting TDEM data in a twodimensional (2D) faulted environment. A subset of the TDEM data consistent with a layered-earth interpretation was identified through a combination of three-dimensional (3D) forward modelling and diffusion time-distance estimates. Inverse modelling indicates an offset in water table elevation of nearly 40 m across the fault. These findings imply that the fault acts as a low-permeability barrier to groundwater flow in the vicinity of the proposed recharge site. Existing production wells on the south side of the fault, together with a thick unsaturated zone and permeable near-surface deposits, suggest the southern half of the study area is suitable for artificial recharge. These results illustrate the effectiveness of targeted TDEM in support of hydrological studies in a heavily faulted desert environment where data are scarce and the cost of obtaining these data by conventional drilling techniques is prohibitive.

Near Surface Geophysics

Analysis of water surface and flow distribution for the design flood at a proposed highway crossing of the Sabine River near Tatum, Texas

The hydraulic effects of the proposed Texas Highway 43 crossing of the Sabine River near Tatum, Texas, were determined on the basis of results from a two-dimensional finite-element surface-water-flow model. In planning the replacement crossing by the Texas State Department of Highways and Public Transportation, approximations of apportionment of flow among the openings and velocities within the openings were of concern. The model was used to simulate flow in the river floodplain system for the proposed design, an alternate design, and for the natural condition. The proposed bridge design by the State Department of Highways and Public Transportation consisted of a 320-foot main channel opening, four left overflow bridges with widths of 120, 320, 320, and 280 feet, and one right overflow opening with a width of 440 feet. The alternate design consisted of a 950-foot main channel opening, three left overflow bridges with widths of 320, 320, and 280 feet, and one right overflow opening with a width of 200 feet. Preliminary one-dimensional computations were used as an aid in establishing the boundary conditions for the two-dimensional analysis. The analysis evaluates the interaction of complex hydraulic characteristics of varied features found at many multiple-opening highway crossings. Valid applications of accepted standard one-dimensional techniques would not reveal some aspects of the overall hydraulics of the crossing. The simulation of the proposed design shows, among other things, which openings are the most productive per foot of bridge length. The results of the two-dimensional simulation of the proposed design indicate some differences in the apportionment of flow among the openings when compared to the one-dimensional proposed crossing computations. The alternate design results in computed water-surface altitudes which were slightly lower than those of the proposed design. The alternate design would require less modification to the existing embankment. Velocities computed within the openings and at the bridge abutments, using the two-dimensional model, were within the design specifications of the State Department of Highways and Public Transportation. The simulations of the proposed and alternate designs indicate a lateral component of the water-surface slope at the embankment. Redistribution of flow across the floodplain also is indicated in both simulations. Some of the differences in the response between the two designs are affected by geometric features of the floodplain other than the embankment-opening geometry.

Water-Resources Investigations Report

Analytical solutions for benchmarking cold regions subsurface water flow and energy transport models: one-dimensional soil thaw with conduction and advection

Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.

Advances in Water Resources

Flood-plain delineation for Bull Run, Little Rocky Run, Johnny Moore Creek, and Popes Head Creek basins, Fairfax County, Virginia

Water-surface profiles of the 25-year and 100-year floods and maps on which the 25-, 50-, and 100-year flood boundaries are delineated for Little Rocky Run Basin, Johnny Moore Creek Basin, Popes Head Creek Basin and that part of Bull Run basin within Fairfax County are presented in this report. The techniques used in the computation of the flood profiles and delineation of flood boundaries are presented, and specific hydraulic problems encountered within the study area are also included.

Virginia

Bathymetry and vegetation in isolated marsh and cypress wetlands in the northern Tampa Bay Area, 2000-2004

Wetland bathymetry and vegetation mapping are two commonly used lines of evidence for assessing the hydrologic and ecologic status of expansive coastal and riverine wetlands. For small isolated freshwater wetlands, however, bathymetric data coupled with vegetation assessments are generally scarce, despite the prevalence of isolated wetlands in many regions of the United States and the recognized importance of topography as a control on inundation patterns and vegetation distribution. In the northern Tampa Bay area of west-central Florida, bathymetry was mapped and vegetation was assessed in five marsh and five cypress wetlands. These 10 isolated wetlands were grouped into three categories based on the effects of ground-water withdrawals from regional municipal well fields: natural (no effect), impaired (drier than natural), and augmented (wetlands with artificially augmented water levels). Delineation of the wetland perimeter was a critical component for estimating wetland-surface area and stored water volume. The wetland perimeter was delineated by the presence of Serenoa repens (the 'palmetto fringe') at 9 of the 10 sites. At the 10th site, where the palmetto fringe was absent, hydric-soils indicators were used to delineate the perimeter. Bathymetric data were collected using one or more techniques, depending on the physical characteristics of each wetland. Wetland stage was measured hourly using continuous stage recorders. Wetland vegetation was assessed semiannually for 2 1/2 years in fixed plots located at three distinct elevations. Vegetation assessments were used to determine the community composition and the relative abundance of obligate, facultative wet, and facultative species at each elevation. Bathymetry maps were generated, and stage-area and stage-volume relations were developed for all 10 wetlands. Bathymetric data sets containing a high density of data points collected at frequent and regular spatial intervals provided the most useful stage-area and stage-volume relations. Bathymetric maps of several wetlands also were generated using a low density of data points collected along transect lines or contour lines. In a comparative analysis of the three mapping approaches, stage-area and stage-volume relations based on transect data alone underestimated (by 50-100 percent over certain ranges of stage) the wetland area and volume compared to results using a high density of data points. Adding data points collected along one elevation contour below the wetland perimeter to the transect data set greatly improved the agreement of the resulting stage-area and stage-volume relations to the high-density mapping approach. Stage-area relations and routinely monitored stage data were used to compare and contrast the average flooded area in a natural marsh and an impaired marsh over a 2-year period. Vegetation assessments used together with flooded-area information provided the potential for extrapolating vegetation results from points or transects to wetlands as a whole. A comparison of the frequency of flooding of different areas of the wetland and the species composition in vegetation plots at different elevations indicated the dependence of vegetation on inundation frequency. Because of the broad tolerances of many wetlands plants to a range of inundation conditions, however, vegetation assessments alone provided less definitive evidence of the hydrologic differences between the two sites, and hydrologic changes occurring during the 2 years, than the flooded-area frequencies. Combining flooded-area frequencies with vegetation assessments could provide a more versatile and insightful approach for determining the ecological status of wetlands than using vegetation and stage data alone. Flooded-area frequencies may further provide a useful approach for assessing the ecological status of wetlands where historical vegetation surveys and stage data are lacking. Comparing the contemporary flooded-area frequencies a

Scientific Investigations Report

Documentation of the dynamic parameter, water-use, stream and lake flow routing, and two summary output modules and updates to surface-depression storage simulation and initial conditions specification options with the Precipitation-Runoff Modeling System (PRMS)

This report documents seven enhancements to the U.S. Geological Survey (USGS) Precipitation-Runoff Modeling System (PRMS) hydrologic simulation code: two time-series input options, two new output options, and three updates of existing capabilities. The enhancements are (1) new dynamic parameter module, (2) new water-use module, (3) new Hydrologic Response Unit (HRU) summary output module, (4) new basin variables summary output module, (5) new stream and lake flow routing module, (6) update to surface-depression storage and flow simulation, and (7) update to the initial-conditions specification. This report relies heavily upon U.S. Geological Survey Techniques and Methods, book 6, chapter B7, which documents PRMS version 4 (PRMS-IV). A brief description of PRMS is included in this report.

Techniques and Methods

Applied Geochemistry Special Issue on Environmental geochemistry of modern mining

Environmental geochemistry is an integral part of the mine-life cycle, particularly for modern mining. The critical importance of environmental geochemistry begins with pre-mining baseline characterization and the assessment of environmental risks related to mining, continues through active mining especially in water and waste management practices, and culminates in mine closure. The enhanced significance of environmental geochemistry to modern mining has arisen from an increased knowledge of the impacts that historical and active mining can have on the environment, and from new regulations meant to guard against these impacts. New regulations are commonly motivated by advances in the scientific understanding of the environmental impacts of past mining. The impacts can be physical, chemical, and biological in nature. The physical challenges typically fall within the purview of engineers, whereas the chemical and biological challenges typically require a multidisciplinary array of expertise including geologists, geochemists, hydrologists, microbiologists, and biologists. The modern mine-permitting process throughout most of the world now requires that potential risks be assessed prior to the start of mining. The strategies for this risk assessment include a thorough characterization of pre-mining baseline conditions and the identification of risks specifically related to the manner in which the ore will be mined and processed, how water and waste products will be managed, and what the final configuration of the post-mining landscape will be. In the Fall 2010, the Society of Economic Geologists held a short course in conjunction with the annual meeting of the Geological Society of America in Denver, Colorado (USA) to examine the environmental geochemistry of modern mining. The intent was to focus on issues that are pertinent to current and future mines, as opposed to abandoned mines, which have been the focus of numerous previous short courses. The geochemical challenges of current and future mines share similarities with abandoned mines, but differences also exist. Mining and ore processing techniques have changed; the environmental footprint of waste materials has changed; environmental protection has become a more integral part of the mine planning process; and most historical mining was done with limited regard for the environment. The 17 papers in this special issue evolved from the Society of Economic Geologists’ short course. The relevant geochemical processes encompass the source, transport, and fate of contaminants related to the life cycle of a mine. Contaminants include metals and other inorganic species derived from geologic sources such as ore and solid mine waste, and substances brought to the site for ore processing, such as cyanide to leach gold. Factors, such as mine-waste mineralogy, hydrologic setting, mine-drainage chemistry, and microbial activity, that affect the hydrochemical risks from mining are reviewed by Nordstrom et al. In another paper, Nordstrom discusses baseline characterization at mine sites in a regulatory framework, and emphasizes the influence of mineral deposits in producing naturally elevated concentrations of many trace elements in surface water and groundwater. Surface water quality in mineralized watersheds is influenced by a number of processes that act on daily (diel) cycles and can produce dramatic variations in trace element concentrations as described by Gammons et al. Pre-mining baseline characterization studies should strive to capture the magnitude of these diel variations. Desbarats et al., using a case study of mine drainage from a gold mine, illustrate how elements that commonly occur as negatively charged species (anions) in solution, such as arsenic as arsenate, behave in an opposite fashion than most metals, which occur as positively charged species (cations). Significant improvement in the understanding of factors that influence the toxicity of metals to aquatic organisms in surface water has highlighted the importance of aqueous chemistry, particularly dissolved organic carbon, as described by Smith et al. Stream sediment contamination is another important pathway for affecting aquatic organisms, as reviewed by Besser et al. Understanding and predicting environmental consequences from mining begins with knowing the mineralogy and mineral reactivity of the ore, the wastes, and of secondary minerals formed later. Jamieson et al. review the importance of mineralogical studies in mine planning and remediation. A number of types of site-specific studies are needed to identify environmental risks related to individual mines. Lapakko reviews the general framework of mine waste characterization studies that are integral to the mine planning process. Hageman et al. present a comparative study of several static tests commonly used to characterize mine waste. The mining and ore processing practices employed at a specific mine site will vary on the basis of the commodities being targeted, the geology of the deposit, the geometry of the deposit, and the mining and ore processing methods used. Thus, these factors, in addition to the waste management practices used, can result in a variety of end-member mine waste features, each of which has its own set of challenges. Open pit mines and underground mines require waste rock to be removed to access ore. Waste rock presents unique problems because the rock is commonly mineralized at sub-economic grades and has not been processed to remove potentially problematic minerals, such as pyrite. Amos et al. examine the salient aspects of the geochemistry of waste rock. Mill tailings – the waste material after ore minerals have been removed – are a volumetrically important solid waste at many mine sites. Their fine grain size and the options for their management make their behavior in the environment distinct from that of waste rock. Lindsay et al. describe some of these differences through three case-study examples. Subaqueous disposal of tailings is another option described by Moncur et al. Cyanide leaching for gold extraction is a common method throughout the world. Johnson describes environmental aspects of cyanidation. Uranium mining presents unique environmental challenges, particularly since in-situ recovery has seen widespread use. Campbell et al. review the environmental geochemistry of uranium mining and current research on bioremediation. Ore concentrates from many types of metal mining undergo a pyrometallurgical technique known as smelting to extract the metal. Slag is the result of smelting, and it may be an environmental liability or a valuable byproduct, as described by Piatak et al. Finally, the open pits that result from surface mining commonly reach below the water table. At the end of mining, these pits may fill to form lakes that become part of the legacy of the mine. Castendyk et al., in two papers, review theoretical aspects of the environmental limnology of pit lakes. They also describe approaches that have been used to model pit lake water balance, wall-rock contributions to pit lake chemistry, pit lake water quality, and limnological processes, such as vertical mixing, through the use of three case studies.

Applied Geochemistry

Biological effects-based tools for monitoring impacted surface waters in the Great Lakes: a multiagency program in support of the Great Lakes Restoration Initiative

There is increasing demand for the implementation of effects-based monitoring and surveillance (EBMS) approaches in the Great Lakes Basin to complement traditional chemical monitoring. Herein, we describe an ongoing multiagency effort to develop and implement EBMS tools, particularly with regard to monitoring potentially toxic chemicals and assessing Areas of Concern (AOCs), as envisioned by the Great Lakes Restoration Initiative (GLRI). Our strategy includes use of both targeted and open-ended/discovery techniques, as appropriate to the amount of information available, to guide a priori end point and/or assay selection. Specifically, a combination of in vivo and in vitro tools is employed by using both wild and caged fish (in vivo), and a variety of receptor- and cell-based assays (in vitro). We employ a work flow that progressively emphasizes in vitro tools for long-term or high-intensity monitoring because of their greater practicality (e.g., lower cost, labor) and relying on in vivo assays for initial surveillance and verification. Our strategy takes advantage of the strengths of a diversity of tools, balancing the depth, breadth, and specificity of information they provide against their costs, transferability, and practicality. Finally, a series of illustrative scenarios is examined that align EBMS options with management goals to illustrate the adaptability and scaling of EBMS approaches and how they can be used in management decisions.

Great Lakes

Methods for quantifying interactions between groundwater and surface water

Driven by the need for integrated management of groundwater (GW) and surface water (SW), quantification of GW–SW interactions and associated contaminant transport has become increasingly important. This is due to their substantial impact on water quantity and quality. In this review, we provide an overview of the methods developed over the past several decades to investigate GW–SW interactions. These methods include geophysical, hydrometric, and tracer techniques, as well as various modeling approaches. Different methods reveal valuable information on GW–SW interactions at different scales with their respective advantages and limitations. Interpreting data from these techniques can be challenging due to factors like scale effects, heterogeneous hydrogeological conditions, sediment variability, and complex spatiotemporal connections between GW and SW. To facilitate the selection of appropriate methods for specific sites, we discuss the strengths, weaknesses, and challenges of each technique, and we offer perspectives on knowledge gaps in the current science.

Annual Review of Environment and Resources

Streamflow at selected gaging stations on the James River in North Dakota and South Dakota, 1953-82, with a section on climatology

Historic stream flow data were compiled and record extension techniques were used, when necessary, to develop a monthly stream flow record for 1953-82 at streamflow-gaging stations on the James River in North Dakota and South Dakota. The record extension techniques included both Maintenance of Variance Extension Type 1 and Ordinary Least Squares. In addition to the historic stream flow record, synthesized unregulated streamflow was computed for gaging stations on the James River for 1953-82 by eliminating the effects of Jamestown Reservoir, Pipestem Reservoir, Sand Lake National Wildlife Refuge, and the consumptive surface-water withdrawals. Maintenance of Variance Extension Type 1, Ordinary Least Squares regression, water-balance procedures, and drainage-area ratio methods were used to compute the unregulated streamflows. Mean annual historic stream flow of the James River at Jamestown is about 7,000 acre-feet less than the mean annual unregulated streamflow. Mean monthly unregulated streamflow is greater than mean monthly historic streamflow during March, April, and May; during all other months, unregulated stream flow is less than historic streamflow. The climate during 1953-82 was compared to the climate records, which range in length from 53 to 102 years, in the James River basin to determine if the climate during 1953-82 was similar to the climate that has occurred during the last 90 years. On the basis of the climate analysis, the data-development period (1953-82) does not represent an abnormally wet or dry period. Thus, the streamflow during the data-development period represents a range of streamflows that might be expected to occur during the last 90 years.

Water-Resources Investigations Report

Interactions between groundwater and surface water in a Virginia coastal plain watershed. 2. Acid-base chemistry

A field study of surface water and groundwater interactions during baseflow and stormflow conditions was performed at the Reedy Creek watershed in the Virginia Coastal Plain. Three estimates of the average saturated hydraulic conductivity ( K s ) of the unconfined aquifer were in reasonable agreement (ranging from 0.0033 to 0.010 cm/s), indicating that baseflow in the creek is entirely from the drainage of shallow groundwater from the relatively thin (1–6 m thick) unconfined aquifer. This relatively permeable surficial aquifer was found to be underlain by dark, olive grey, clay-silt and diatomaceous Miocene deposits of low permeability known as the Calvert Formation, which is believed to function as a confining bed in the area. A chemical hydrograph separation technique was used to resolve the contributions of [old] (pre-event) and [new] (event) water to stormflow. Results from a major rainstorm indicated that old water dominated the stormflow response of the watershed, although the new water contribution approached 40% at the hydrograph peak. Stormflow at Reedy Creek appears to result from saturation overland flow from variable source areas which include the stream channels and a significant part of the riparian wetland area. This response appears to be attributable to the transient dynamics of the shallow groundwater flow system and to the formation of localized groundwater mounds which raise the water-table to the wetland surface.

Hydrological Processes

Documentation of the Unsaturated-Zone Flow (UZF1) Package for modeling Unsaturated Flow Between the Land Surface and the Water Table with MODFLOW-2005

Percolation of precipitation through unsaturated zones is important for recharge of ground water. Rain and snowmelt at land surface are partitioned into different pathways including runoff, infiltration, evapotranspiration, unsaturated-zone storage, and recharge. A new package for MODFLOW-2005 called the Unsaturated-Zone Flow (UZF1) Package was developed to simulate water flow and storage in the unsaturated zone and to partition flow into evapotranspiration and recharge. The package also accounts for land surface runoff to streams and lakes. A kinematic wave approximation to Richards? equation is solved by the method of characteristics to simulate vertical unsaturated flow. The approach assumes that unsaturated flow occurs in response to gravity potential gradients only and ignores negative potential gradients; the approach further assumes uniform hydraulic properties in the unsaturated zone for each vertical column of model cells. The Brooks-Corey function is used to define the relation between unsaturated hydraulic conductivity and water content. Variables used by the UZF1 Package include initial and saturated water contents, saturated vertical hydraulic conductivity, and an exponent in the Brooks-Corey function. Residual water content is calculated internally by the UZF1 Package on the basis of the difference between saturated water content and specific yield. The UZF1 Package is a substitution for the Recharge and Evapotranspiration Packages of MODFLOW-2005. The UZF1 Package differs from the Recharge Package in that an infiltration rate is applied at land surface instead of a specified recharge rate directly to ground water. The applied infiltration rate is further limited by the saturated vertical hydraulic conductivity. The UZF1 Package differs from the Evapotranspiration Package in that evapotranspiration losses are first removed from the unsaturated zone above the evapotranspiration extinction depth, and if the demand is not met, water can be removed directly from ground water whenever the depth to ground water is less than the extinction depth. The UZF1 Package also differs from the Evapotranspiration Package in that water is discharged directly to land surface whenever the altitude of the water table exceeds land surface. Water that is discharged to land surface, as well as applied infiltration in excess of the saturated vertical hydraulic conductivity, may be routed directly as inflow to specified streams or lakes if these packages are active; otherwise, this water is removed from the model. The UZF1 Package was tested against the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model for a vertical unsaturated flow problem that includes evapotranspiration losses. This report also includes an example in which MODFLOW-2005 with the UZF1 Package was used to simulate a realistic surface-water/ground-water flow problem that includes time and space variable infiltration, evapotranspiration, runoff, and ground-water discharge to land surface and to streams. Another simpler problem is presented so that the user may use the input files as templates for new problems and to verify proper code installation.

Techniques and Methods

Extending the capture map concept to estimate discrete and risk-based streamflow depletion potential

A popular and contemporary use of numerical groundwater models is to estimate the discrete relation between groundwater extraction and surface-water/groundwater exchange. Previously, the concept of a “capture map” has been put forward as a means to effectively summarize this relation for decision-making consumption. While capture maps have enjoyed success in the environmental simulation industry, they are deterministic, ignoring uncertainty in the underlying model. Furthermore, capture maps are not typically calculated in a manner that facilitates analysis of varying combinations of extraction locations and/or reaches. That is, they are typically constructed with focus on a single reach or group of reaches. The former of these limitations is important for conveying risk to decision makers and stakeholders, while the latter is important for decision-making support related to surface-water management, where future foci may include reaches that were not the focus of the original capture analysis. Herein, we use the concept of a response matrix to generalize the theory of the capture-map approach to estimate spatially discrete streamflow depletion potential. We also use first-order, second-moment uncertainty estimation techniques with the concept of “risk shifting” to place capture maps and streamflow depletion potential in a stochastic, risk-based framework. Our approach is demonstrated for an integrated groundwater/surface-water model of the lower San Antonio River, Texas, USA.

Groundwater

Evaluation of land subsidence and ground failures at Bicycle Basin, Fort Irwin National Training Center, California, 1992–2017

Groundwater has been pumped in the Bicycle Basin at Fort Irwin National Training Center since the 1960s, and the amount pumped has generally increased since the 1990s. After a large crack (approximately 0.5-kilometer long) formed at the surface of Bicycle Lake playa during 2005–06 in the area used as an aircraft runway, a monitoring study was initiated by the U.S. Geological Survey, in cooperation with the U.S. Army Fort Irwin National Training Center, to help determine the cause. The extent and effect of groundwaterlevel declines and land-surface deformation in Bicycle Basin were evaluated using a number of approaches, including water-level measurements from December 2007 to June 2017, land surveys across the playa area, interferometric synthetic aperture radar (InSAR) analyses, geophysical surveys of the playa area, and numerical experiments to test hypotheses about soil mechanical processes. A specific objective of this study was to evaluate the recent development of ground failures in the form of large, surface cracks that pose a hazard to aircraft operations on the Bicycle Lake playa airstrip. Another large crack and set of interconnected cracks formed at the surface of Bicycle Lake playa in 2013 following a period of inundation similar to that of 2005–06; however, the 2013 cracks formed a network rather than a single feature. Groundwater pumping resulted in more than 27 meters (90 feet) of water-level declines in wells north of the Bicycle Lake playa from 1990 to 2017, and InSAR interferograms indicated more than 400 millimeters (16 inches) of subsidence during 1993–2015. Subsidence rates calculated from InSAR interferograms were variable, temporally and spatially. Results of leveling surveys indicated differential subsidence between 2009 and 2016 along a transect across the playa to the area of maximum subsidence; there was less subsidence south of the 2005–06 crack than north of it. The steepest subsidence gradient for this time was in the area of maximum subsidence. Repeat tape-extensometer measurements from April 2009 to November 2017 across the main-runway crack indicated slight opening along part of the crack. A baseline lidar (light detection and ranging) survey of the main-runway crack was done in January 2009 to track the development of this feature. During the 95-months from January 2009 to December 2016, 0.5 meters (19.7 inches) of subsidence was observed next to the subsidence pit in the crack along the western end of the scanned area, presumably due to erosion and backfilling when the lake was inundated. Standing water in the crack along the eastern end of the scanned area prevented determining change for that part of the crack. Time-series, shaded-relief images show the progression of change in the crack from January 2009 through December 2016, with the crack closing and opening, depending on the time of most recent inundation relative to when the scans were taken. Volumetric changes in the size of the crack, calculated between the sequential lidar surveys, show that from 2009 to 2016, the cumulative volumetric change was 1.5 cubic meters (54 cubic feet). This volume is a minimum estimate because the crack volume below the water surface in the pits was not included in the calculations. Subsurface imaging of the 2005–06 main-runway crack by a series of electrical resistivity tomography surveys in 2008 and 2017 indicated that the crack could extend 5 meters (16 feet) in the subsurface. Electromagnetic induction surveys in 2008, 2014, and 2015 evaluated the technique for ground-failure monitoring and showed that the technique was effective at revealing anomalies correlated with the features of concern in the playa area, such as the 2005–06 crack, numerous ”healed” macropolygon features, and scattered sink-like depressions. Results from numerical experiments simulating watertable decline at depth indicated that the material deposits in the desaturating capillary fringe zone might not be able to transmit large enough stresses up through the overlying soil to cause cracks at the land surface. Results from simulations of desiccation in the presence of a regional tectonically induced stress field, however, tended to support the hypothesis that the combined processes could control the formation of giant desiccation macropolygons in certain areas, such as Bicycle Lake playa.

California

Spatially consistent high-resolution land surface temperature mosaics for thermophysical mapping of the Mojave Desert

Daytime and nighttime thermal infrared observations acquired by the ASTER and MODIS instruments onboard the NASA Terra spacecraft have produced a dataset that can be used to map thermophysical properties across large regions, which have implications on surface processes, thermal environments and habitat suitability for desert species. ASTER scenes acquired between 2004 and 2012 are combined using new mosaicking and data-fusion techniques to produce a map of daytime and nighttime land surface temperature with coverage exclusive of the effects of clouds and weather. These data are combined with Landsat 7 visible imagery to generate a consistent map of apparent thermal inertia (ATI), which is related to the presence of exposed bedrock, rocks, fine-grained sediments and water on the surface. The resulting datasets are compared to known geomorphic units and surface types to generate an interpreted mechanical composition map of the entire Mojave Desert at 100 m per pixel that is most sensitive to large clast size distinctions in grain size distribution.

ARizone, California, Nevada, Utah