Geology Reports⌕ Search

Geology topics

C. Amanda Garcia

Publications and source records attributed to C. Amanda Garcia.

At least 19 recordsLinked to original sources

Water-level change from a multiple-well aquifer test in volcanic rocks, Umatilla Indian Reservation near Mission, northeastern Oregon, 2016

The U.S. Geological Survey, in cooperation with the Confederated Tribes of the Umatilla Indian Reservation (CTUIR), (1) estimated water-level change from a multiple-well aquifer test centered on CTUIR well number 422 and (2) evaluated hydraulic connections between the pumping and observation wells on the Umatilla Indian Reservation near Mission, northeastern Oregon to improve the understanding of aquifer characteristics and hydrologic flow boundaries. Water-level changes, or pumping responses, were determined by distinguishing the pumping signal from environmental fluctuations in groundwater levels using analytical water-level models. The pumping well produces water from basalt units from a depth of 450 to 1,057 feet below land surface and was intermittently pumped during February 1–April 18, 2016. Water-level responses to pumping were estimated in the pumping well and in seven observation wells within 4 miles (mi) of the pumping well. The observation wells are open to basalt and some observation wells are either separated from the pumping well by faults and other structural features, within structural zones, or adjacent to structural features. Pumping responses at the observation wells were classified as detected in two wells, ambiguous in one well, and not detected in four wells. Observation-well open-interval elevations overlapped with the pumping-well open interval in both wells with detected pumping responses. Observation wells with detections are 1.8 mi east of the pumping well and across a fault, and 1.4 mi south of the pumping well. The pumping response was classified as ambiguous in an observation well located 1.4 mi west of the pumping well, where the dip of the basalt unit steepens, and adjacent to the Agency syncline. Pumping responses were not detected in observation wells within 0.3 mi of the pumping well where observation-well open-interval elevations are above the top of the pumping well open interval. Analysis of pumping responses indicates (1) a more permeable zone of basalt is adjacent to the lower portion of the pumping-well open interval and extends eastward, (2) basalt adjacent to the upper portion of the pumping-well open-interval is less permeable than the lower portion or separated from the lower portion by a less permeable zone, and (or) (3) a less permeable zone limits vertical hydraulic connectivity between the pumping well and the overlying basalt.

Oregon↗

Groundwater resources of the Harney Basin, southeastern Oregon

In response to increasing groundwater demand and declining groundwater levels in the Harney Basin of southeastern Oregon, the U.S. Geological Survey and the Oregon Water Resources Department conducted a cooperative groundwater-availability study during 2016–22. This Fact Sheet summarizes the results of this study. Full details of the study are provided in Gingerich and others (2022a, 2022b), Garcia and others (2022), and the other supporting documents listed on the last page of this Fact Sheet.

Oregon↗

Groundwater resources of the Harney Basin, southeastern Oregon

Groundwater development has increased substantially in southeastern Oregon’s Harney Basin since 2010, mainly for the purpose of large-scale irrigation. Concurrently, some areas of the basin experienced groundwater-level declines of more than 100 feet, and some shallow wells have gone dry. The Oregon Water Resources Department has limited new groundwater development in the basin until an improved understanding of the groundwater-flow system is available. This report describes the results of a hydrologic investigation undertaken to provide that understanding. The investigation encompasses the groundwater hydrology of the entire 5,240-square-mile Harney Basin. Most of the precipitation in the Harney Basin falls in the higher-elevation areas of the Blue Mountains and Steens Mountain. Although considerable groundwater recharge occurs in these upland areas, most (83 percent) re-emerges as streams and springs in the uplands. Groundwater recharge in the lowlands is provided through infiltration of surface water flowing onto the lowlands from rivers and streams leaving the uplands and as groundwater flow from the surrounding upland rocks. Water-balance calculations indicate that the rate of groundwater recharge to the Harney Basin lowlands (where most groundwater is withdrawn) averages 173,000 acre-feet per year (acre-ft/yr). Groundwater in the Harney Basin lowlands mainly discharges through evapotranspiration from groundwater-irrigated (supplied from wells) crops or from natural vegetation drawing groundwater from the shallow water table and capillary fringe. Groundwater discharge in the lowlands is estimated to be about 283,000 acre-ft/yr, which exceeds the estimated groundwater recharge to the lowlands by about 110,000 acre-ft/yr. This imbalance results in removal of groundwater from storage in the aquifer system and is evidenced by the large declines observed in groundwater levels in the areas of greatest groundwater pumpage. To a large degree, the location and depth of pumpage dictate the timing and distribution of the effects of groundwater use in the Harney Basin. Pumpage is commonly greatest in the areas where higher-permeability geologic units allow for higher well yields. However, many of these higher-permeability units are bounded by lower-permeability units that cannot supply groundwater at a sufficient rate to replenish the areas of greatest pumpage, resulting in groundwater-level declines. Three Harney Basin areas with a combined area exceeding 140 square miles have experienced groundwater-level declines exceeding 40 feet compared to pre-development conditions: near the Weaver Spring/Dog Mountain area, in the northeastern floodplains along Highway 20, and near Crane. Areas of more modest groundwater-level decline (about 10 feet) were identified in the Virginia Valley area and the Silver Creek floodplain north of Riley. Smaller localized areas of groundwater-level depression have also formed around individual wells or groups of wells throughout the Harney Basin lowlands. Most groundwater being pumped from the Harney Basin lowlands, including all three areas experiencing large groundwater-level declines, was recharged more than 12,000 years ago, near the end of the last glacial period when the climate in the basin was cooler and wetter than today. Geochemical evidence indicates that modern recharge generally circulates to a depth no greater than 100 feet below the floodplains of major rivers and streams in the lowlands. Away from the major river and stream corridors, pre-modern water commonly is found at the water table. Recharge to groundwater and recovery of groundwater levels in the most heavily pumped areas in the Harney Basin lowlands are restricted by the limited spatial extent and depth of modern recharge in the Harney Basin lowlands and the relatively fine-grained deposits underlying most of the lowland areas.

Oregon↗

Hydrologic budget of the Harney Basin groundwater system, southeastern Oregon

Groundwater-level declines and limited quantitative knowledge of the groundwater-flow system in the Harney Basin prompted a cooperative study between the U.S. Geological Survey and the Oregon Water Resources Department to evaluate the groundwater-flow system and budget. This report provides a hydrologic budget of the Harney Basin groundwater system that includes separate groundwater budgets for upland and lowland areas to avoid double counting water that recharges in the uplands, discharges to streams and springs in the uplands, flows downstream to the lowlands, and recharges the lowland groundwater system. Lowlands generally represent the conterminous valleys within the center of the basin, including floodplains of the major streams and uplands represent all other areas in the basin. The upland groundwater budget is minimally affected by groundwater development and generally represents the budget of the natural system. In upland areas during 1982–2016, mean-annual recharge totaled 288,000 acre-feet (acre-ft) and mean-annual discharge totaled 239,000 acre-ft, resulting in a net recharge of 49,000 acre-ft. Upland groundwater recharge occurs as infiltration of precipitation and snowmelt and was estimated using the USGS Soil-Water-Balance model calibrated to estimates of runoff, evapotranspiration (ET), base flow, and snow-water equivalent. Groundwater discharge to streams is the predominant discharge mechanism in upland areas and was estimated as 225,000 acre-feet per year (acre-ft/yr) during 1982–2016 using hydrograph separation and summer low-flow estimates in streamgaged watersheds and a linear relation between estimated streamflow and base flow in ungaged watersheds. The remaining upland discharge occurs through springs (14,000 acre-ft/yr) that either emerge downgradient of locations where groundwater discharge to streams was estimated or are routed to irrigated areas. Spring discharge was estimated as a compilation of current and historical measurements. The net upland recharge, which is 17 percent of total upland recharge, ultimately recharges lowland areas as groundwater flow from uplands to lowlands. The lowland groundwater budget for the Harney Basin represents a combination of natural conditions and human activity as more than 99 percent of groundwater development has occurred either inside or within 2 miles of the lowland boundary. In lowland areas during 1982–2016, mean annual groundwater recharge totaled 173,000 acre-ft and groundwater discharge totaled 283,000 acre-ft, indicating discharge exceeded recharge by more than 60 percent. Excluding groundwater pumping, the lowland groundwater budget is more in balance with a mean annual recharge of 165,000 acre-ft and a mean annual discharge of 131,000 acre-ft during 1982–2016. The 23-percent difference between non-pumping recharge and discharge mostly represents the cumulative uncertainty in the estimates of the various groundwater budget components but also likely includes a small reduction in natural groundwater discharge captured by pumping. Lowland groundwater is predominantly recharged by infiltration of surface water (116,000 acre-ft/yr) through streams, floodwater, and irrigation, with a lesser amount as groundwater inflow from uplands and minimal recharge beneath Malheur and Harney Lakes. Recharge from streams and floodwater (natural and irrigation) was estimated using a balance of measured and estimated surface-water inflow to and outflow from lowland areas including streamflow, springflow, and ET where a portion of surface-water inflow to lowland areas is comprised of upland discharge to streams and springs. Groundwater ET (119,000 acre-ft/yr) is the predominant natural discharge mechanism in lowland areas and was estimated as the mean from two remote-sensing based approaches incorporating groundwater ET measurements from other similar basins and 23 years (1987–2015) of Landsat imagery. Discharge of lowland groundwater into Malheur and Harney Lakes is about 700 acre-ft/yr and is represented in groundwater ET estimates. The remaining natural groundwater discharge from lowland areas issues from Sodhouse Spring (8,900 acre-ft/yr) and as groundwater flow to the Malheur River Basin through Virginia Valley (3,100 acre-ft/yr). The relatively large amount of groundwater discharged to springs in Warm Springs Valley (25,000 acre-ft/yr) is accounted for in groundwater ET estimates. Natural groundwater discharge in lowland areas of the Harney Basin has remained relatively constant during the last 80 years based on comparisons with estimates north of Malheur Lake and west of Harney Lake published in the 1930s. Annual net amount of groundwater pumped (pumpage) from the Harney Basin during 2017–18 averaged 144,000 acre-ft. The net value is the difference between pumpage (about 152,000 acre-ft/yr) and reinfiltration of groundwater pumped for irrigation and non-irrigation purposes (about 8,000 acre-ft/yr). Net pumpage was estimated in concurrent studies that compiled groundwater-use data and coupled reported groundwater pumpage data from wells with remote-sensing-based ET estimates from groundwater-irrigated fields. Total pumpage for irrigation has increased from about 54,000 acre-ft/yr during 1991–92 to 145,000 acre-ft/yr during 2017–18. Presently, pumpage is greatest in the lowland region north of Malheur Lake (81,000 acre-ft/yr), with lesser amounts to the north and northwest of Harney Lake (41,000 acre-ft/yr) and to the south and east of Malheur Lake (22,000 acre-ft/yr). During this study, mean annual lowland groundwater discharge (including pumpage) exceeded mean annual recharge, indicating that the lowland hydrologic budget is out of balance. Net groundwater pumpage during 2017–18 is similar to groundwater discharge from all other sources in the lowlands and is four times the imbalance between non-pumping lowland recharge and discharge (34,000 acre-ft/yr). Declining groundwater levels at depth across many parts of the Harney Basin lowlands indicate that pumpage is depleting aquifer storage and is likely capturing a small amount of natural groundwater discharge to springs and ET in some lowland areas. If pumping continues, aquifer storage depletion will continue until the capture rate of natural discharge to springs and ET is equal to the pumping rate. If groundwater development occurs in upland areas and reduces either the streamflow or groundwater inflow to lowland areas, the deficit in the lowland water budget will increase.

Oregon↗

Hydraulic characterization of carbonate-rock and basin-fill aquifers near Long Canyon, Goshute Valley, northeastern Nevada

Understanding groundwater flow and pumping effects near pending mining operations requires accurate subsurface hydraulic characterization. To improve conceptual models of groundwater flow and development in the complex hydrogeologic system near Long Canyon Mine, in northwestern Goshute Valley, northeastern Nevada, the U.S. Geological Survey characterized the hydraulic properties of carbonate rocks and basin-fill aquifers using an integrated analysis of steady-state and stressed aquifer conditions informed by water chemistry and aquifer-test data. Hydraulic gradients and groundwater-age data in northern Goshute Valley indicate carbonate rocks in the Pequop Mountains just west and south of the Long Canyon Mine project area constitute a more permeable and active flow system than saturated rocks in the northern Pequop Mountains, western Toano Range, and basin fill. Permeable carbonate rocks in the northern Pequop Mountains, in part, discharge to the Johnson Springs wetland complex (JSWC), where mean groundwater ages range from 500 to 2,400 years and samples all contain a small fraction of modern waters, relative to mean ages of 8,600 to more than 22,000 years for most groundwater sampled to the north and east. Recharge to the JSWC occurs from a roughly 27-square-mile area in the upgradient Pequop Mountains to the west, composed mostly of permeable carbonate rock and fractured quartzite, and bounded by low-permeability shales and marbleized and siliclastic rocks. Single-well aquifer-test analyses provided transmissivity estimates at pumped wells. Transmissivity estimates ranged from 7,000 to 400,000 feet squared per day (ft 2 /d) in carbonate rocks and from 2,000 to 80,000 ft 2 /d in basin fill near the Long Canyon Mine. Water-level drawdown from multiple-well aquifer testing and rise from unintentional leakage into the overlying basin-fill aquifer were estimated and distinguished from natural fluctuations in 93 pumping and monitoring sites using analytical water-level models. Leakage of disposed aquifer-test pumpage occurred south of the aquifer test area through an unlined irrigation ditch. Drawdown was detected at distances of as much as 3 miles (mi) from pumping wells at all but one carbonate-rock site, at basin-fill sites on the alluvial fan immediately downgradient from pumping wells, and in Big Spring and spring NS-05. Similar drawdowns in carbonate rocks within the drawdown detection area suggest all wells penetrate a highly transmissive zone (HTZ) that is bounded by low-permeability rocks. Drawdown was not detected in carbonate rocks to the west of Canyon fault, in any basin-fill sites on the valley floor east of the Hardy fault, or at volcanic sites to the north, indicating that these major fault structures and (or) permeability contrasts between hydrogeologic units impeded groundwater flow or obscured pumping signals. Alternatively, unintentional leakage might have obscured drawdown at basin-fill sites on the valley floor, where water-level rise was detected at nine sites over 3 mi. Consistent hydraulic properties were estimated by simultaneously interpreting steady-state flow during predevelopment conditions and changes in groundwater levels and springflows from the 2016 carbonate-rock aquifer test with an integrated groundwater-flow model. Hydraulic properties were distributed across carbonate rocks, basin fill, volcanic rocks, and siliciclastic rocks with a hydrogeologic framework developed from geologic mapping and hydraulic testing. Estimated transmissivity distributions spanned at least three orders of magnitude in each rock unit. In the HTZ, simulated transmissivities ranged from 10,000 to 23,000,000 ft 2 /d, with the most transmissive areas occurring around Big Spring. Comparatively low carbonate-rock transmissivities of less than 10,000 ft 2 /d were estimated in the northern Pequop Mountains and poorly defined values of less than 1,000 ft 2 /d were estimated in the western Toano Range. Transmissivities in basin fill ranged from less than 10 to 80,000 ft 2 /d and were minimally constrained by the 2016 carbonate-rock aquifer test because poorly quantified leakage affected water levels more so than pumping. The most transmissive areas were informed by single-well aquifer tests along the eastern edge of the Pequop Mountains near Long Canyon Mine and could be indicative of a hydraulic connection between basin fill and more transmissive underlying carbonate rocks. Simulated transmissivities of volcanic and low-permeability rocks mostly are less than 1,000 ft 2 /d. The estimated hydraulic-property distributions and informed interpretation of hydraulic connections among hydrogeologic units improved the characterization and representation of groundwater flow near the Long Canyon Mine.

Nevada↗

Evaluation of hydrologic impact of an irrigation curtailment program in the Upper Klamath Lake Basin using Landsat satellite data

Upper Klamath Lake (UKL) is the source of the Klamath river that flows through southern Oregon and northern California. The UKL basin is home to two endangered species and provides water for 81,000+ ha (200,000+ acres) of irrigation on the United States Bureau of Reclamation (USBR) Klamath Project located downstream of the UKL basin. Irrigated agriculture also occurs along the tributaries to UKL. During 2013–2016, water right calls resulted in various levels of curtailment of irrigation diversions from the tributaries to UKL. However, information on the extent of curtailment, how much irrigation water was saved, and its impact on the UKL is unknown. In this study, we combined Landsat-based actual evapotranspiration (ETa) data obtained from the Operational Simplified Surface Energy Balance (SSEBop) model with gridded precipitation and USGS station discharge data to evaluate the hydrologic impact of the curtailment program. Analysis was performed for five base years (2004, 2006, 2008-2010) and four target years (2013-2016) over irrigated areas above UKL. Our results indicated that the impact of the curtailment program over the June to September time-period was highest during 2013 and declined in each of the following years. The total on-field water savings were approximately 60 hm3 in 2013 and 2014, 44 hm3 in 2015, and 32 hm3 in 2016. The instream water flow change or extra water available (EWA) were found at 92, 68, 45, and 26 hm3 respectively for 2013, 2014, 2015 and 2016. Most water savings came from pasture and wetlands. Alfalfa showed the most decline in water use among grain crops. The resulting EWA from the curtailment contributed to a maximum of 19% of the lake inflows and 50% of the lake volume. This study presents the use of Landsat-based ETa and other remote sensing datasets for evaluating water-related impacts of the irrigation curtailment program.

Upper Klamath Lake Basin↗

A process to estimate net infiltration using a site-scale water-budget approach, Rainier Mesa, Nevada National Security Site, Nevada, 2002–05

This report documents a process used to estimate net infiltration from precipitation, evapotranspiration (ET), and soil data acquired at two sites on Rainier Mesa. Rainier Mesa is a groundwater recharge area within the Nevada National Security Site where recharged water flows through bedrock fractures to a deep (450 meters) water table. The U.S. Geological Survey operated two ET stations on Rainier Mesa from 2002 to 2005 at sites characterized by pinyon-juniper and scrub-brush vegetative cover. Precipitation and ET data were corrected to remove measurement biases and gap-filled to develop continuous datasets. Net infiltration (percolation below the root zone) and changes in root-zone water storage were estimated using a monthly water-balance model. Site-scale water-budget results indicate that the heavily-fractured welded-tuff bedrock underlying thin (<40 centimeters) topsoil is a critical water source for vegetation during dry periods. Annual precipitation during the study period ranged from fourth lowest (182 millimeters [mm]) to second highest (708 mm) on record (record = 55 years). Annual ET exceeded precipitation during dry years, indicating that the fractured-bedrock reservoir capacity is sufficient to meet atmospheric-evaporative demands and to sustain vegetation through extended dry periods. Net infiltration (82 mm) was simulated during the wet year after the reservoir was rapidly filled to capacity. These results support previous conclusions that preferential fracture flow was induced, resulting in an episodic recharge pulse that was detected in nearby monitoring wells. The occurrence of net infiltration only during the wet year is consistent with detections of water-level rises in nearby monitoring wells that occur only following wet years.

Nevada↗

Groundwater discharge by evapotranspiration, flow of water in unsaturated soil, and stable isotope water sourcing in areas of sparse vegetation, Amargosa Desert, Nye County, Nevada

This report documents methodology and results of a study to evaluate groundwater discharge by evapotranspiration (GWET) in sparsely vegetated areas of Amargosa Desert and improve understanding of hydrologic-continuum processes controlling groundwater discharge. Evapotranspiration and GWET rates were computed and characterized at three sites over 2 years using a combination of micrometeorological, unsaturated zone, and stable-isotope measurements. One site (Amargosa Flat Shallow [AFS]) was in a sparse and isolated area of saltgrass ( Distichlis spicata ) where the depth to groundwater was 3.8 meters (m). The second site (Amargosa Flat Deep [AFD]) was in a sparse cover of predominantly shadscale ( Atriplex confertifolia ) where the depth to groundwater was 5.3 m. The third site (Amargosa Desert Research Site [ADRS]), selected as a control site where GWET is assumed to be zero, was located in sparse vegetation dominated by creosote bush ( Larrea tridentata ) where the depth to groundwater was 110 m. Results indicated that capillary rise brought groundwater to within 0.9 m (at AFS) and 3 m (at AFD) of land surface, and that GWET rates were largely controlled by the slow but relatively persistent upward flow of water through the unsaturated zone in response to atmospheric-evaporative demands. Greater GWET at AFS (50 ± 20 millimeters per year [mm/yr]) than at AFD (16 ± 15 mm/yr) corresponded with its shallower depth to the capillary fringe and constantly higher soil-water content. The stable-isotope dataset for hydrogen (δ 2 H) and oxygen (δ 18 O) illustrated a broad range of plant-water-uptake scenarios. The AFS saltgrass and AFD shadscale responded to changing environmental conditions and their opportunistic water use included the time- and depth-variable uptake of unsaturated-zone water derived from a combination of groundwater and precipitation. These results can be used to estimate GWET in other areas of Amargosa Desert where hydrologic conditions are similar.

Nevada↗

Hydraulic characterization of volcanic rocks in Pahute Mesa using an integrated analysis of 16 multiple-well aquifer tests, Nevada National Security Site, 2009–14

An improved understanding of groundwater flow and radionuclide migration downgradient from underground nuclear-testing areas at Pahute Mesa, Nevada National Security Site, requires accurate subsurface hydraulic characterization. To improve conceptual models of flow and transport in the complex hydrogeologic system beneath Pahute Mesa, the U.S. Geological Survey characterized bulk hydraulic properties of volcanic rocks using an integrated analysis of 16 multiple-well aquifer tests. Single-well aquifer-test analyses provided transmissivity estimates at pumped wells. Transmissivity estimates ranged from less than 1 to about 100,000 square feet per day in Pahute Mesa and the vicinity. Drawdown from multiple-well aquifer testing was estimated and distinguished from natural fluctuations in more than 200 pumping and observation wells using analytical water-level models. Drawdown was detected at distances greater than 3 miles from pumping wells and propagated across hydrostratigraphic units and major structures, indicating that neither faults nor structural blocks noticeably impede or divert groundwater flow in the study area. Consistent hydraulic properties were estimated by simultaneously interpreting drawdown from the 16 multiple-well aquifer tests with an integrated groundwater-flow model composed of 11 well-site models—1 for each aquifer test site. Hydraulic properties were distributed across volcanic rocks with the Phase II Pahute Mesa-Oasis Valley Hydrostratigraphic Framework Model. Estimated hydraulic-conductivity distributions spanned more than two orders of magnitude in hydrostratigraphic units. Overlapping hydraulic conductivity ranges among units indicated that most Phase II Hydrostratigraphic Framework Model units were not hydraulically distinct. Simulated total transmissivity ranged from 1,600 to 68,000 square feet per day for all pumping wells analyzed. High-transmissivity zones exceeding 10,000 square feet per day exist near caldera margins and extend along the northern and eastern Pahute Mesa study area and near the southwestern edge of the study area. The estimated hydraulic-property distributions and observed hydraulic connections among geologic structures improved the characterization and representation of groundwater flow at Pahute Mesa.

Nevada↗

Budgets and chemical characterization of groundwater for the Diamond Valley flow system, central Nevada, 2011–12

The Diamond Valley flow system consists of six hydraulically connected hydrographic areas in central Nevada. The general down-gradient order of the areas are southern and northern Monitor Valleys, Antelope Valley, Kobeh Valley, Stevens Basin, and Diamond Valley. Groundwater flow in the Diamond Valley flow system terminates at a large playa in the northern part of Diamond Valley. Concerns relating to continued water-resources development of the flow system resulted in a phased hydrologic investigation that began in 2005 by the U.S. Geological Survey in cooperation with Eureka County. This report presents the culmination of the phased investigation to increase understanding of the groundwater resources of the basin-fill aquifers in the Diamond Valley flow system through evaluations of groundwater chemistry and budgets. Groundwater chemistry was characterized using major ions and stable isotopes from groundwater and precipitation samples. Groundwater budgets accounted for all inflows, outflows, and changes in storage, and were developed for pre-development (pre-1950) and recent (average annual 2011&ndash;12) conditions. Major budget components include groundwater discharge by evapotranspiration and groundwater withdrawals; groundwater recharge by precipitation, and interbasin flow; and storage change. Groundwater in the basin-fill aquifer of the Diamond Valley flow system was mostly a calcium or sodium bicarbonate water type and generally within acceptable drinking-water standards. The general water type was similar among the individual hydrographic areas. Stable isotopes of oxygen-18 and deuterium from precipitation varied seasonally, such that enrichment from evaporation was greater during warmer months than cooler months. The isotopic signature of shallow groundwater was similar to cool season precipitation, indicating recharge was relatively recent (similar to recent climatic conditions) and was derived from cool season precipitation. Site-scale groundwater evapotranspiration was estimated from eddy-covariance and micrometeorological measurements collected at four sites and ranged from 0.15 feet per year in sparse, undisturbed shrubland to 1.13 feet per year in a grassland meadow. Vegetation indices calculated from satellite imagery and field mapping were used to define three evapotranspiration units (shrubland, grassland, and playa) and to extrapolate site-scale groundwater evapotranspiration rates to basin-scale estimates. Annual pre-development groundwater evapotranspiration for individual hydrographic areas ranged from 2,900 acre-feet per year (acre-ft/yr) in northern Monitor Valley to 35,000 acre-ft/yr in Diamond Valley. Total groundwater evapotranspiration from the Diamond Valley flow system under pre-development conditions was about 70,000 acre-ft/yr. Areas of irrigated land in the Diamond Valley flow system increased from less than 5,000 acres in the early 1960s to more than 25,000 acres in 2012 and are mostly for growing alfalfa in southern Diamond Valley. Annual (2011&ndash;12) net groundwater withdrawals for irrigation, assumed to be the volume of groundwater consumed by crops and pastureland, ranged from about 420 acre-ft/yr in Antelope Valley to 67,000 acre-ft/yr in Diamond Valley. Total net groundwater withdrawals for irrigation in the Diamond Valley flow system were about 69,000 acre-ft/yr (2011&ndash;12). Groundwater recharge, the largest inflow component to the Diamond Valley flow system, was determined as the sum of groundwater evapotranspiration and net subsurface outflow (subsurface outflow minus subsurface inflow). Annual groundwater recharge estimates ranged from 200 acre-ft/yr in Stevens Basin to 35,000 acre-ft/yr in Diamond Valley. Subsurface flow between hydrographic basins was evaluated using estimated transmissivity, groundwater-flow sections derived from remotely sensed imagery, and hydraulic gradients determined from 2012 water-level data. Subsurface outflow ranged from 0 acre-ft/yr for Diamond Valley to 3,400 acre-ft/yr for northern Monitor Valley into western Kobeh Valley. Subsurface inflow ranged from 0 acre-ft/yr for southern Monitor Valley to 4,200 acre-ft/yr for Kobeh Valley from northern Monitor and Antelope Valleys. The pre-development, steady state, groundwater budget for the Diamond Valley flow system was estimated at about 70,000 acre-ft/yr of inflow and outflow. During years 2011&ndash;12, inflow components of groundwater recharge from precipitation and subsurface inflow from adjacent basins totaled 70,000 acre-ft/yr for the DVFS, whereas outflow components included 64,000 acre-ft/yr of groundwater evapotranspiration and 69,000 acre-ft/yr of net groundwater withdrawals, or net pumpage. Spring discharge in northern Diamond Valley declined about 6,000 acre-ft/yr between pre-development time and years 2011&ndash;12. Assuming net groundwater withdrawals minus spring flow decline is equivalent to the storage change, the 2011&ndash;12 summation of inflow and storage change was balanced with outflow at about 133,000 acre-ft/yr.

Nevada↗

Groundwater discharge by evapotranspiration, Dixie Valley, west-central Nevada, March 2009-September 2011

With increasing population growth and land-use change, urban communities in the desert Southwest are progressively looking toward remote basins to supplement existing water supplies. Pending applications by Churchill County for groundwater appropriations from Dixie Valley, Nevada, a primarily undeveloped basin east of the Carson Desert, have prompted a reevaluation of the quantity of naturally discharging groundwater. The objective of this study was to develop a revised, independent estimate of groundwater discharge by evapotranspiration (ET g ) from Dixie Valley using a combination of eddy-covariance evapotranspiration (ET) measurements and multispectral satellite imagery. Mean annual ET g was estimated during water years 2010 and 2011 at four eddy-covariance sites. Two sites were in phreatophytic shrubland dominated by greasewood, and two sites were on a playa. Estimates of total ET and ET g were supported with vegetation cover mapping, soil physics considerations, water‑level measurements from wells, and isotopic water sourcing analyses to allow partitioning of ET g into evaporation and transpiration components. Site-based ET g estimates were scaled to the basin level by combining remotely sensed imagery with field reconnaissance. Enhanced vegetation index and brightness temperature data were compared with mapped vegetation cover to partition Dixie Valley into five discharging ET units and compute basin-scale ET g . Evapotranspiration units were defined within a delineated groundwater discharge area and were partitioned as (1) playa lake, (2) playa, (3) sparse shrubland, (4) moderate-to-dense shrubland, and (5) grassland. Groundwater ET is influenced primarily by phreatophytic vegetative cover, salinity of soil and groundwater within the playa, depth to groundwater, solar radiation, and air temperature. The annual groundwater contribution to site‑scale ET ranged from 24 to 61 percent of total ET at vegetated sites and 4 to 15 percent of total ET at playa sites. Mean annual ET g from vegetated sites ranged from 53 millimeters (mm) (0.17 foot [ft], 7.3 percent vegetative cover) to 225 mm (0.74 ft, 24.8 percent vegetative cover). Cumulative liquid‑water fluxes in the unsaturated zone indicate that ET g at vegetated sites was influenced primarily by plant transpiration. Binary mixing analyses of oxygen-18 isotopes in groundwater and shallow soil water indicate that plants predominantly use groundwater throughout the year. Groundwater fractions in greasewood stem water varied seasonally and ranged from 0.63 to 1.0. Mean annual playa ET g ranged from about 11 mm (0.04 ft) at the inner playa site (near-surface volumetric water content of 37–53 percent) to about 20 mm (0.07 ft) at the outer playa site located within 2 kilometers of the playa edge (near-surface volumetric water content of 25–38 percent), but playa ET g estimates were within the probable error (plus or minus [±] 20–23 mm; 0.06–0.08 ft). Varying playa ET g was influenced predominantly by salinity rather than depth to groundwater. Osmotic resistance and physical impediments to ET (such as surface salt crusts and salt precipitate in the soil pore space) increased with increasing salinity toward the playa center, whereas vapor pressure decreased. Mean annual basin-scale ET g totaled about 28 million cubic meters (Mm 3 ) (23,000 acre-feet [acre-ft]), and represents the sum of ET g from all ET units. Annual groundwater ET from vegetated areas totaled about 26 Mm 3 (21,000 acre-ft), and was dominated by the moderate-to-dense shrubland ET unit (54 percent), followed by sparse shrubland (37 percent) and grassland (9 percent) ET units. Senesced grasses observed in the northern most areas of the moderate-to-dense ET unit likely confounded the vegetation index and led to an overestimate of ET g for this ET unit. Therefore, mean annual ET g for moderate-to-dense shrubland presented here is likely an upper bound. Annual groundwater ET from the playa ET unit was 2.2 Mm 3 (1,800 acre-ft), whereas groundwater ET from the playa lake ET unit was 0–0.1 Mm 3 (0–100 acre-ft). Oxygen-18 and deuterium data indicate discharge from the playa center predominantly represents removal of local precipitation-derived recharge. The playa lake estimate, therefore, is considered an upper bound. Mean annual ET g estimates for Dixie Valley are assumed to represent the pre‑development, long-term ET g rates within the study area.

Nevada↗

Hydrogeologic framework and occurrence, movement, and chemical characterization of groundwater in Dixie Valley, west-central Nevada

Dixie Valley, a primarily undeveloped basin in west-central Nevada, is being considered for groundwater exportation. Proposed pumping would occur from the basin-fill aquifer. In response to proposed exportation, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation and Churchill County, conducted a study to improve the understanding of groundwater resources in Dixie Valley. The objective of this report is to characterize the hydrogeologic framework, the occurrence and movement of groundwater, the general water quality of the basin-fill aquifer, and the potential mixing between basin-fill and geothermal aquifers in Dixie Valley. Various types of geologic, hydrologic, and geochemical data were compiled from previous studies and collected in support of this study. Hydrogeologic units in Dixie Valley were defined to characterize rocks and sediments with similar lithologies and hydraulic properties influencing groundwater flow. Hydraulic properties of the basin-fill deposits were characterized by transmissivity estimated from aquifer tests and specific-capacity tests. Groundwater-level measurements and hydrogeologic-unit data were combined to create a potentiometric surface map and to characterize groundwater occurrence and movement. Subsurface inflow from adjacent valleys into Dixie Valley through the basin-fill aquifer was evaluated using hydraulic gradients and Darcy flux computations. The chemical signature and groundwater quality of the Dixie Valley basin-fill aquifer, and potential mixing between basin-fill and geothermal aquifers, were evaluated using chemical data collected from wells and springs during the current study and from previous investigations. Dixie Valley is the terminus of the Dixie Valley flow system, which includes Pleasant, Jersey, Fairview, Stingaree, Cowkick, and Eastgate Valleys. The freshwater aquifer in the study area is composed of unconsolidated basin-fill deposits of Quaternary age. The basin-fill hydrogeologic unit can be several orders of magnitude more transmissive than surrounding and underlying consolidated rocks and Dixie Valley playa deposits. Transmissivity estimates in the basin fill throughout Dixie Valley ranged from 30 to 45,500 feet squared per day; however, a single transmissivity value of 0.1 foot squared per day was estimated for playa deposits. Groundwater generally flows from the mountain range uplands toward the central valley lowlands and eventually discharges near the playa edge. Potentiometric contours east and west of the playa indicate that groundwater is moving eastward from the Stillwater Range and westward from the Clan Alpine Mountains toward the playa. Similarly, groundwater flows from the southern and northern basin boundaries toward the basin center. Subsurface groundwater flow likely enters Dixie Valley from Fairview and Stingaree Valleys in the south and from Jersey and Pleasant Valleys in the north, but groundwater connections through basin-fill deposits were present only across the Fairview and Jersey Valley divides. Annual subsurface inflow from Fairview and Jersey Valleys ranges from 700 to 1,300 acre-feet per year and from 1,800 to 2,300 acre-feet per year, respectively. Groundwater flow between Dixie, Stingaree, and Pleasant Valleys could occur through less transmissive consolidated rocks, but only flow through basin fill was estimated in this study. Groundwater in the playa is distinct from the freshwater, basin-fill aquifer. Groundwater mixing between basin-fill and playa groundwater systems is physically limited by transmissivity contrasts of about four orders of magnitude. Total dissolved solids in playa deposit groundwater are nearly 440 times greater than total dissolved solids in the basin-fill groundwater. These distinctive physical and chemical flow restrictions indicate that groundwater interaction between the basin fill and playa sediments was minimal during this study period (water years 2009–11). Groundwater in Dixie Valley generally can be characterized as a sodium bicarbonate type, with greater proportions of chloride north of the Dixie Valley playa, and greater proportions of sulfate south of the playa. Analysis of major ion water chemistry data sampled during the study period indicates that groundwater north and south of Township 22N differ chemically. Dixie Valley groundwater quality is marginal when compared with national primary and secondary drinking-water standards. Arsenic and fluoride concentrations exceed primary drinking water standards, and total dissolved solids and manganese concentrations exceed secondary drinking water standards in samples collected during this study. High concentrations of boron and tungsten also were observed. Chemical comparisons between basin-fill and geothermal aquifer water indicate that most basin-fill groundwater sampled could contain 10–20 percent geothermal water. Geothermal indicators such as high temperature, lithium, boron, chloride, and silica suggest that mixing occurs in many wells that tap the basin-fill aquifer, particularly on the north, south, and west sides of the basin. Magnesium-lithium geothermometers indicate that some basin-fill aquifer water sampled for the current study likely originates from water that was heated above background mountain-block recharge temperatures (between 3 and 15 degrees Celsius), highlighting the influence of mixing with warm water that was possibly derived from geothermal sources.

Nevada↗

Advanced methods for modeling water-levels and estimating drawdowns with SeriesSEE, an Excel add-in

Water-level modeling is used for multiple-well aquifer tests to reliably differentiate pumping responses from natural water-level changes in wells, or &ldquo;environmental fluctuations.&rdquo; Synthetic water levels are created during water-level modeling and represent the summation of multiple component fluctuations, including those caused by environmental forcing and pumping. Pumping signals are modeled by transforming step-wise pumping records into water-level changes by using superimposed Theis functions. Water-levels can be modeled robustly with this Theis-transform approach because environmental fluctuations and pumping signals are simulated simultaneously. Water-level modeling with Theis transforms has been implemented in the program SeriesSEE, which is a Microsoft&reg; Excel add-in. Moving average, Theis, pneumatic-lag, and gamma functions transform time series of measured values into water-level model components in SeriesSEE. Earth tides and step transforms are additional computed water-level model components. Water-level models are calibrated by minimizing a sum-of-squares objective function where singular value decomposition and Tikhonov regularization stabilize results. Drawdown estimates from a water-level model are the summation of all Theis transforms minus residual differences between synthetic and measured water levels. The accuracy of drawdown estimates is limited primarily by noise in the data sets, not the Theis-transform approach. Drawdowns much smaller than environmental fluctuations have been detected across major fault structures, at distances of more than 1 mile from the pumping well, and with limited pre-pumping and recovery data at sites across the United States. In addition to water-level modeling, utilities exist in SeriesSEE for viewing, cleaning, manipulating, and analyzing time-series data.

Techniques and Methods↗

Assessing hydraulic connections across a complex sequence of volcanic rocks - Analysis of U-20 WW multiple-well aquifer test, Pahute Mesa, Nevada National Security Site, Nevada

Groundwater beneath Pahute Mesa flows through a complexly layered sequence of volcanic rock aquifers and confining units that have been faulted into distinct structural blocks. Hydraulic property estimates of rocks and structures in this flow system are necessary to assess radionuclide migration near underground nuclear testing areas. The U.S. Geological Survey (USGS) used a 12 month (October 1, 2008&mdash; October 1, 2009) intermittent pumping schedule of well U-20 WW and continuously monitored water levels in observation wells ER-20-6 #3, UE-20bh 1, and U-20bg as a multi-well aquifer test to evaluate hydraulic connections across structural blocks, bulk hydraulic properties of volcanic rocks, and the hydraulic significance of a major fault. Measured water levels were approximated using synthetic water levels generated from an analytical model. Synthetic water levels are a summation of environmental water-level fluctuations and a Theis (1935) transform of the pumping signal from flow rate to water-level change. Drawdown was estimated by summing residual differences between measured and synthetic water levels and the Theis-transformed pumping signal from April to September 2009. Drawdown estimates were used in a three‑dimensional numerical model to estimate hydraulic properties of distinct aquifers, confining units, and a major fault. A maximum water-level drawdown of nearly 0.4 foot in well UE-20bh 1, which is more than 1 mile from the pumping well, was detected across a major fault. Drawdown estimates in the observation well nearest to (ER-20-6 #3, less than 1 mile) and within the same structural block as the pumping well were less than detection (<0.1 foot). Minimal drawdown within the same structural block indicates that lava units are likely separated by bedded tuff confining units. Hydraulic property estimates indicate that wells U-20 WW, UE-20bh 1, and ER-20-6 #3 produce water from moderately permeable fractured lava, as hydraulic conductivity and specific storage estimates average 4.8 feet per day and 2.1&times;10 &ndash;6 per foot, respectively, and transmissivity estimates range from 1,200 to 3,600 feet squared per day. Sensitivity analyses indicate that the major fault is hydraulically similar to the permeable host rock and connects flow between structural blocks.

Nevada↗

Interpretation of Flow Logs from Nevada Test Site Boreholes to Estimate Hydraulic Conductivity Using Numerical Simulations Constrained by Single-Well Aquifer Tests

Hydraulic conductivities of volcanic and carbonate lithologic units at the Nevada Test Site were estimated from flow logs and aquifer-test data. Borehole flow and drawdown were integrated and interpreted using a radial, axisymmetric flow model, AnalyzeHOLE. This integrated approach is used because complex well completions and heterogeneous aquifers and confining units produce vertical flow in the annular space and aquifers adjacent to the wellbore. AnalyzeHOLE simulates vertical flow, in addition to horizontal flow, which accounts for converging flow toward screen ends and diverging flow toward transmissive intervals. Simulated aquifers and confining units uniformly are subdivided by depth into intervals in which the hydraulic conductivity is estimated with the Parameter ESTimation (PEST) software. Between 50 and 150 hydraulic-conductivity parameters were estimated by minimizing weighted differences between simulated and measured flow and drawdown. Transmissivity estimates from single-well or multiple-well aquifer tests were used to constrain estimates of hydraulic conductivity. The distribution of hydraulic conductivity within each lithology had a minimum variance because estimates were constrained with Tikhonov regularization. AnalyzeHOLE simulated hydraulic-conductivity estimates for lithologic units across screened and cased intervals are as much as 100 times less than those estimated using proportional flow-log analyses applied across screened intervals only. Smaller estimates of hydraulic conductivity for individual lithologic units are simulated because sections of the unit behind cased intervals of the wellbore are not assumed to be impermeable, and therefore, can contribute flow to the wellbore. Simulated hydraulic-conductivity estimates vary by more than three orders of magnitude across a lithologic unit, indicating a high degree of heterogeneity in volcanic and carbonate-rock units. The higher water transmitting potential of carbonate-rock units relative to volcanic-rock units is exemplified by the large difference in their estimated maximum hydraulic conductivity; 4,000 and 400 feet per day, respectively. Simulated minimum estimates of hydraulic conductivity are inexact and represent the lower detection limit of the method. Minimum thicknesses of lithologic intervals also were defined for comparing AnalyzeHOLE results to hydraulic properties in regional ground-water flow models.

Scientific Investigations Report↗

Transport of tritium contamination to the atmosphere in an arid environment

Soil–plant–atmosphere interactions strongly influence water movement in desert unsaturated zones, but little is known about how such interactions affect atmospheric release of subsurface water-borne contaminants. This 2-yr study, performed at the U.S. Geological Survey's Amargosa Desert Research Site in southern Nevada, quantified the magnitude and spatiotemporal variability of tritium (3H) transport from the shallow unsaturated zone to the atmosphere adjacent to a low-level radioactive waste (LLRW) facility. Tritium fluxes were calculated as the product of 3H concentrations in water vapor and respective evaporation and transpiration water-vapor fluxes. Quarterly measured 3H concentrations in soil water vapor and in leaf water of the dominant creosote-bush [ Larrea tridentat a (DC.) Coville] were spatially extrapolated and temporally interpolated to develop daily maps of contamination across the 0.76-km2 study area. Maximum plant and root-zone soil concentrations (4200 and 8700 Bq L−1, respectively) were measured 25 m from the LLRW facility boundary. Continuous evaporation was estimated using a Priestley–Taylor model and transpiration was computed as the difference between measured eddy-covariance evapotranspiration and estimated evaporation. The mean evaporation/transpiration ratio was 3:1. Tritium released from the study area ranged from 0.12 to 12 μg d−1 and totaled 1.5 mg (8.2 × 1010 Bq) over 2 yr. Tritium flux variability was driven spatially by proximity to 3H source areas and temporally by changes in 3H concentrations and in the partitioning between evaporation and transpiration. Evapotranspiration removed and limited penetration of precipitation beneath native vegetation and fostered upward movement and release of 3H from below the root zone.

Vadose Zone Journal↗

Portable chamber measurements of evapotranspiration at the Amargosa Desert Research Site near Beatty, Nye County, Nevada, 2003-06

Portable chamber measurements of evapotranspiration (ET) were made at the U.S. Geological Survey's Amargosa Desert Research Site in southern Nevada to help quantify component- and landscape-scale contributions to ET in an arid environment. Evapotranspiration data were collected approximately every 3 months from 2003 to 2006. Chamber measurements of ET were partitioned into bare-soil evaporation and mixed-species transpiration components. The component-scale ET fluxes from native shrubs typically surpassed those from bare soil by as much as a factor of four. Component-scale ET fluxes were extrapolated to landscape-scale ET using a one-layer, multi-component canopy model. Landscape-scale ET fluxes predominantly were controlled by bare-soil evaporation. Bare soil covered 94 percent of the landscape on average and contributed about 70 percent of the landscape-scale vapor flux. Creosote bush, an evergreen shrub, accounted for about 90 percent of transpiration on average due to its dominance across the landscape (80 percent of the 6 percent shrub cover) and evergreen character.

Nevada↗

Selected micrometeorological, soil-moisture, and evapotranspiration data at Amargosa Desert Research Site in Nye County near Beatty, Nevada, 2001-05

Selected micrometeorological and soil-moisture data were collected at the Amargosa Desert Research Site adjacent to a low-level radio-active waste and hazardous chemical waste facility near Beatty, Nevada, 2001-05. Evapotranspiration data were collected from February 2002 through the end of December 2005. Data were col-lected in support of ongoing research to improve the understanding of hydrologic and contaminant-transport processes in arid environments. Micrometeorological data include solar radiation, net radiation, air temperature, relative humidity, saturated and ambient vapor pressure, wind speed and direction, barometric pressure, precipitation, near-surface soil temperature, soil-heat flux and soil-water content. All micrometeorological data were collected using a 10-second sampling interval by data loggers that output daily and hourly mean values. Daily maximum and minimum values are based on hourly mean values. Precipitation data output includes daily and hourly totals. Selected soil-moisture profiles at depth include periodic measurements of soil volumetric water-content measurements at nine neutron-probe access tubes to depths ranging from 5.25 to 29.25 meters. Evapotranspiration data include measurement of daily evapotranspiration and 15-minute fluxes of the four principal energy budget components of latent-heat flux, sensible-heat flux, soil-heat flux, and net radiation. Other data collected and used in equations to determine evapotranspiration include temperature and water content of soil, temperature and vapor pressure of air, and covariance values. Evapotranspiration and flux estimates during 15-minute intervals were calculated at a 0.1-second execution interval using the eddy covariance method. Data files included in this report contain the complete micrometeorological, soil-moisture, and evapotranspiration field data sets. These data files are presented in tabular Excel spreadsheet format. This report highlights selected data contained in the computer generated data files using figures, tables, and brief discussions. Instrumentation used for data collection also is described. Water-content profiles are shown to demonstrate variability of water content with depth. Time-series data are plotted to illustrate temporal variations in micrometeorological, soil-water content, and evapotranspiration data.

Nevada↗