Search USGS⌕ Search

SEARCH · Search USGS

Results for “Scientific Investigations Map”

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

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

At least 1,531 records · Page 85Linked to original sources
USGS2006

Geohydrological characterization, water-chemistry, and ground-water flow simulation model of the Sonoma Valley area, Sonoma County, California

The Sonoma Valley, located about 30 miles north of San Francisco, is one of several basins in Sonoma County that use a combination of ground water and water delivered from the Russian River for supply. Over the past 30 years, Sonoma Valley has experienced rapid population growth and land-use changes. In particular, there has been a significant increase in irrigated agriculture, predominantly vineyards. To provide a better understanding of the ground-water/surface-water system in Sonoma Valley, the U.S. Geological Survey compiled and evaluated existing data, collected and analyzed new data, and developed a ground-water flow model to better understand and manage the ground-water system. The new data collected include subsurface lithology, gravity measurements, groundwater levels, streamflow gains and losses, temperature, water chemistry, and stable isotopes. Sonoma Valley is drained by Sonoma Creek, which discharges into San Pablo Bay. The long-term average annual volume of precipitation in the watershed is estimated to be 269,000 acre-feet. Recharge to the ground-water system is primarily from direct precipitation and Sonoma Creek. Discharge from the ground-water system is predominantly outflow to Sonoma Creek, pumpage, and outflow to marshlands and to San Pablo Bay. Geologic units of most importance for groundwater supply are the Quaternary alluvial deposits, the Glen Ellen Formation, the Huichica Formation, and the Sonoma Volcanics. In this report, the ground-water system is divided into three depth-based geohydrologic units: upper (less than 200 feet below land surface), middle (between 200 and 500 feet), and lower (greater than 500 feet). Synoptic streamflow measurements were made along Sonoma Creek and indicate those reaches with statistically significant gains or losses. Changes in ground-water levels in wells were analyzed by comparing historical contour maps with the contour map for 2003. In addition, individual hydrographs were evaluated to assess temporal changes by region. In recent years, pumping depressions have developed southeast of Sonoma and southwest of El Verano. Water-chemistry data for samples collected from 75 wells during 2002-04 indicate that the ground-water quality in the study area generally is acceptable for potable use. The water from some wells, however, contains one or more constituents in excess of the recommended standards for drinking water. The chemical composition of water from creeks, springs, and wells sampled for major ions plot within three groups on a trilinear diagram: mixed-bicarbonate, sodium-mixed anion, and sodium-bicarbonate. An area of saline ground water in the southern part of the Sonoma Valley appears to have shifted since the late 1940s and early 1950s, expanding in one area, but receding in another. Sparse temperature data from wells southwest of the known occurrence of thermal water suggest that thermal water may be present beneath a larger part of the valley than previously thought. Thermal water contains higher concentrations of dissolved minerals than nonthermal waters because mineral solubilities generally increase with temperature. Geohydrologic Characterization, Water-Chemistry, and Ground-Water Flow Simulation Model of the Sonoma Valley Area, Sonoma County, California Oxygen-18 (d18 O) and deuterium (dD) values for water from most wells plot along the global meteoric water line, indicating that recharge primarily is derived from the direct infiltration of precipitation or the infiltration of seepage from creeks. Samples from shallow- and intermediate-depth wells located near Sonoma Creek and (or) in the vicinity of Shellville plot to the right of the global meteoric water line, indicating that these waters are partly evaporated. The d18 O and dD composition of water from sampled wells indicates that water from wells deeper than 200 feet is isotopically lighter (more negative) than water from wells less than 200 feet deep, possibly indicating that older ground wate

California↗
USGS2013

Simulation of groundwater flow, effects of artificial recharge, and storage volume changes in the Equus Beds aquifer near the city of Wichita, Kansas well field, 1935–2008

The Equus Beds aquifer is a primary water-supply source for Wichita, Kansas and the surrounding area because of shallow depth to water, large saturated thickness, and generally good water quality. Substantial water-level declines in the Equus Beds aquifer have resulted from pumping groundwater for agricultural and municipal needs, as well as periodic drought conditions. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project to store and later recover groundwater, and to form a hydraulic barrier to the known chloride-brine plume near Burrton, Kansas. In October 2009, the U.S. Geological Survey, in cooperation with the city of Wichita, began a study to determine groundwater flow in the area of the Wichita well field, and chloride transport from the Arkansas River and Burrton oilfield to the Wichita well field. Groundwater flow was simulated for the Equus Beds aquifer using the three-dimensional finite-difference groundwater-flow model MODFLOW-2000. The model simulates steady-state and transient conditions. The groundwater-flow model was calibrated by adjusting model input data and model geometry until model results matched field observations within an acceptable level of accuracy. The root mean square (RMS) error for water-level observations for the steady-state calibration simulation is 9.82 feet. The ratio of the RMS error to the total head loss in the model area is 0.049 and the mean error for water-level observations is 3.86 feet. The difference between flow into the model and flow out of the model across all model boundaries is -0.08 percent of total flow for the steady-state calibration. The RMS error for water-level observations for the transient calibration simulation is 2.48 feet, the ratio of the RMS error to the total head loss in the model area is 0.0124, and the mean error for water-level observations is 0.03 feet. The RMS error calculated for observed and simulated base flow gains or losses for the Arkansas River for the transient simulation is 7,916,564 cubic feet per day (91.6 cubic feet per second) and the RMS error divided by (/) the total range in streamflow (7,916,564/37,461,669 cubic feet per day) is 22 percent. The RMS error calculated for observed and simulated streamflow gains or losses for the Little Arkansas River for the transient simulation is 5,610,089 cubic feet per day(64.9 cubic feet per second) and the RMS error divided by the total range in streamflow (5,612,918/41,791,091 cubic feet per day) is 13 percent. The mean error between observed and simulated base flow gains or losses was 29,999 cubic feet per day (0.34 cubic feet per second) for the Arkansas River and -1,369,250 cubic feet per day (-15.8 cubic feet per second) for the Little Arkansas River. Cumulative streamflow gain and loss observations are similar to the cumulative simulated equivalents. Average percent mass balance difference for individual stress periods ranged from -0.46 to 0.51 percent. The cumulative mass balance for the transient calibration was 0.01 percent. Composite scaled sensitivities indicate the simulations are most sensitive to parameters with a large areal distribution. For the steady-state calibration, these parameters include recharge, hydraulic conductivity, and vertical conductance. For the transient simulation, these parameters include evapotranspiration, recharge, and hydraulic conductivity. The ability of the calibrated model to account for the additional groundwater recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project was assessed by using the U.S. Geological Survey subregional water budget program ZONEBUDGET and comparing those results to metered recharge for 2007 and 2008 and previous estimates of artificial recharge. The change in storage between simulations is the volume of water that estimates the recharge credit for the aquifer storage and recovery system. The estimated increase in storage of 1,607 acre-ft in the basin storage area compared to metered recharge of 1,796 acre-ft indicates some loss of metered recharge. Increased storage outside of the basin storage area of 183 acre-ft accounts for all but 6 acre-ft or 0.33 percent of the total. Previously estimated recharge credits for 2007 and 2008 are 1,018 and 600 acre-ft, respectively, and a total estimated recharge credit of 1,618 acre-ft. Storage changes calculated for this study are 4.42 percent less for 2007 and 5.67 percent more for 2008 than previous estimates. Total storage change for 2007 and 2008 is 0.68 percent less than previous estimates. The small difference between the increase in storage from artificial recharge estimated with the groundwater-flow model and metered recharge indicates the groundwater model correctly accounts for the additional water recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project. Small percent differences between inflows and outflows for all stress periods and all index cells in the basin storage area, improved calibration compared to the previous model, and a reasonable match between simulated and measured long-term base flow indicates the groundwater model accurately simulates groundwater flow in the study area. The change in groundwater level through recent years compared to the August 1940 groundwater level map has been documented and used to assess the change of storage volume of the Equus Beds aquifer in and near the Wichita well field for three different areas. Two methods were used to estimate changes in storage from simulation results using simulated change in groundwater levels in layer 1 between stress periods, and using ZONEBUDGET to calculate the change in storage in the same way the effects of artificial recharge were estimated within the basin storage area. The three methods indicate similar trends although the magnitude of storage changes differ. Information about the change in storage in response to hydrologic stresses is important for managing groundwater resources in the study area. The comparison between the three methods indicates similar storage change trends are estimated and each could be used to determine relative increases or decreases in storage. Use of groundwater level changes that do not include storage changes that occur in confined or semi-confined parts of the aquifer will slightly underestimate storage changes; however, use of specific yield and groundwater level changes to estimate storage change in confined or semi-confined parts of the aquifer will overestimate storage changes. Using only changes in shallow groundwater levels would provide more accurate storage change estimates for the measured groundwater levels method. The value used for specific yield is also an important consideration when estimating storage. For the Equus Beds aquifer the reported specific yield ranges between 0.08 and 0.35 and the storage coefficient (for confined conditions) ranges between 0.0004 and 0.16. Considering the importance of the value of specific yield and storage coefficient to estimates of storage change over time, and the wide range and substantial overlap for the reported values for specific yield and storage coefficient in the study area, further information on the distribution of specific yield and storage coefficient within the Equus Beds aquifer in the study area would greatly enhance the accuracy of estimated storage changes using both simulated groundwater level, simulated groundwater budget, or measured groundwater level methods.

Kansas↗
USGS2011

Water budgets and groundwater volumes for abandoned underground mines in the Western Middle Anthracite Coalfield, Schuylkill, Columbia, and Northumberland Counties, Pennsylvania: Preliminary estimates with identification of data needs

This report, prepared in cooperation with the Pennsylvania Department of Environmental Protection (PaDEP), the Eastern Pennsylvania Coalition for Abandoned Mine Reclamation, and the Dauphin County Conservation District, provides estimates of water budgets and groundwater volumes stored in abandoned underground mines in the Western Middle Anthracite Coalfield, which encompasses an area of 120 square miles in eastern Pennsylvania. The estimates are based on preliminary simulations using a groundwater-flow model and an associated geographic information system that integrates data on the mining features, hydrogeology, and streamflow in the study area. The Mahanoy and Shamokin Creek Basins were the focus of the study because these basins exhibit extensive hydrologic effects and water-quality degradation from the abandoned mines in their headwaters in the Western Middle Anthracite Coalfield. Proposed groundwater withdrawals from the flooded parts of the mines and stream-channel modifications in selected areas have the potential for altering the distribution of groundwater and the interaction between the groundwater and streams in the area. Preliminary three-dimensional, steady-state simulations of groundwater flow by the use of MODFLOW are presented to summarize information on the exchange of groundwater among adjacent mines and to help guide the management of ongoing data collection, reclamation activities, and water-use planning. The conceptual model includes high-permeability mine voids that are connected vertically and horizontally within multicolliery units (MCUs). MCUs were identified on the basis of mine maps, locations of mine discharges, and groundwater levels in the mines measured by PaDEP. The locations and integrity of mine barriers were determined from mine maps and groundwater levels. The permeability of intact barriers is low, reflecting the hydraulic characteristics of unmined host rock and coal. A steady-state model was calibrated to measured groundwater levels and stream base flow, the latter at many locations composed primarily of discharge from mines. Automatic parameter estimation used MODFLOW-2000 with manual adjustments to constrain parameter values to realistic ranges. The calibrated model supports the conceptual model of high-permeability MCUs separated by low-permeability barriers and streamflow losses and gains associated with mine infiltration and discharge. The simulated groundwater levels illustrate low groundwater gradients within an MCU and abrupt changes in water levels between MCUs. The preliminary model results indicate that the primary result of increased pumping from the mine would be reduced discharge from the mine to streams near the pumping wells. The intact barriers limit the spatial extent of mine dewatering. Considering the simulated groundwater levels, depth of mining, and assumed bulk porosity of 11 or 40 percent for the mined seams, the water volume in storage in the mines of the Western Middle Anthracite Coalfield was estimated to range from 60 to 220 billion gallons, respectively. Details of the groundwater-level distribution and the rates of some mine discharges are not simulated well using the preliminary model. Use of the model results should be limited to evaluation of the conceptual model and its simulation using porous-media flow methods, overall water budgets for the Western Middle Anthracite Coalfield, and approximate storage volumes. Model results should not be considered accurate for detailed simulation of flow within a single MCU or individual flooded mine. Although improvements in the model calibration were possible by introducing spatial variability in permeability parameters and adjusting barrier properties, more detailed parameterizations have increased uncertainty because of the limited data set. The preliminary identification of data needs includes continuous streamflow, mine discharge rate, and groundwater levels in the mines and adjacent areas. Data collected whe

Pennsylvania↗
USGS2011

Groundwater conditions in the Brunswick-Glynn County area, Georgia, 2009

The Upper Floridan aquifer is contaminated with saltwater in a 2-square-mile area of downtown Brunswick, Georgia. The presence of this saltwater has limited the development of the groundwater supply in the Glynn County area. Hydrologic, geologic, and water-quality data are needed to effectively manage water resources. Since 1959, the U.S. Geological Survey (USGS) has conducted a cooperative water program with the City of Brunswick and Glynn County to monitor and assess the effect of groundwater development on saltwater intrusion within the Floridan aquifer system. The potential development of alternative sources of water in the Brunswick and surficial aquifer systems also is an important consideration in coastal areas. During calendar year 2009, the cooperative water program included continuous water-level recording of 13 wells completed in the Floridan, Brunswick, and surficial aquifer systems; collecting water levels from 46 wells to map the potentiometric surface of the Upper Floridan aquifer in Glynn County during August 2009; and collecting and analyzing water samples from 55 wells completed in the Floridan aquifer system, of which 27 wells were used to map chloride concentrations in the upper water-bearing zone of the Upper Floridan aquifer in the Brunswick area during August 2009. Periodic water-level measurements also were collected from two wells completed in the Upper Floridan aquifer and four wells completed in the Brunswick aquifer system on Jekyll Island. Equipment was installed on one well to enable real-time specific conductance monitoring in the area surrounding the chloride plume. During 2008-2009, water levels in 30 of the 32 wells monitored in the Brunswick-Glynn County area rose at a rate of 0.24 to 7.58 feet per year (ft/yr). The largest rise of 7.58 ft/yr was in the Upper Floridan aquifer. These rises corresponded to a period of above normal precipitation and decreased pumping. Declines during 2008-2009 were recorded in wells completed in the Brunswick aquifer system (0.37 ft/yr) and Lower Floridan aquifer (0.83 ft/yr). Chloride data collected by two local industrial groundwater users at their well fields since 1958 were compiled and compared with data collected by the USGS during the same period. The results indicate that chloride concentrations at the two well fields have continued to rise despite modification of production wells to eliminate deep saline zones and decreases in pumpage at both facilities. One of the industrial users, Pinova Inc., plugged the lower portions of nine production wells in the mid to late 1960s, which generally decreased chloride concentrations to less than 100 milligrams per liter (mg/L) for a period of 10 to 20 years. However, chloride concentrations eventually returned to previous levels despite decreases in pumpage. During 1990-2009, chloride concentrations at the other industrial user's well field (Georgia-Pacific Cellulose LLC) generally increased despite a 16 million gallon per day decrease in pumpage during this period. Data from the Georgia-Pacific Cellulose well field and additional chloride data from USGS observation wells located to the east indicate continued movement of chloride from the source area located southeast of the site toward the well field.

Georgia↗
USGS2014

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

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

Idaho, Oregon, Washington↗
USGS2018

Hydrogeologic framework and groundwater conditions of the Ararat Basin in Armenia

Armenia is a landlocked country located in the mountainous Caucasus region between Asia and Europe. It shares borders with the countries of Georgia on the north, Azerbaijan on the east, Iran on the south, and Turkey and Azerbaijan on the west. The Ararat Basin is a transboundary basin in Armenia and Turkey. The Ararat Basin (or Ararat Valley) is an intermountain depression that contains the Aras River and its tributaries, which also form the border between Armenia and Turkey and divide the basin into northern and southern regions. The Ararat Basin also contains Armenia’s largest agricultural and fish farming zone that is supplied by high-quality water from wells completed in the artesian aquifers that underlie the basin. Groundwater constitutes about 40 percent of all water use, and groundwater provides 96 percent of the water used for drinking purposes in Armenia. Since 2000, groundwater withdrawals and consumption in the Ararat Basin of Armenia have increased because of the growth of aquaculture and other uses. Increased groundwater withdrawals caused decreased springflow, reduced well discharges, falling water levels, and a reduction of the number of flowing artesian wells in the southern part of Ararat Basin in Armenia. In 2016, the U.S. Geological Survey and the U.S. Agency for International Development (USAID) began a cooperative study in Armenia to share science and field techniques to increase the country’s capabilities for groundwater study and modeling. The purpose of this report is to describe the hydrogeologic framework and groundwater conditions of the Ararat Basin in Armenia based on data collected in 2016 and previous hydrogeologic studies. The study area includes the Ararat Basin in Armenia. This report was completed through a partnership with USAID/Armenia in the implementation of its Science, Technology, Innovation, and Partnerships effort through the Advanced Science and Partnerships for Integrated Resource Development program and associated partners, including the Government of Armenia, Armenia’s Hydrogeological Monitoring Center, and the USAID Global Development Lab and its GeoCenter. The hydrogeologic framework of the Ararat Basin includes several basin-fill stratigraphic units consisting of interbedded dense clays, gravels, sands, volcanic basalts, and andesite deposits. Previously published cross sections and well lithologic logs were used to map nine general hydrogeologic units. Hydrogeologic units were mapped based on lithology and water-bearing potential. Water-level data measured in the water-bearing hydrogeologic units 2, 4, 6, and 8 in 2016 were used to create potentiometric surface maps. In hydrogeologic unit 2, the estimated direction of groundwater flow is from the west to north in the western part of the basin (away from the Aras River) and from north to south (toward the Aras River) in the eastern part of the basin. In hydrogeologic unit 4, the direction of groundwater flow is generally from west to east and north to south (toward the Aras River) except in the western part of the basin where groundwater flow is toward the north or northwest. Hydrogeologic unit 6 has the same general pattern of groundwater flow as unit 4. Hydrogeologic unit 8 is the deepest of the water-bearing units and is confined in the basin. Groundwater flow generally is from the south to north (away from the Aras River) in the western part of the basin and from west to east and north to south (toward the Aras River) elsewhere in the basin. In addition to water levels, personnel from Armenia’s Hydrogeological Monitoring Center also measured specific conductance at 540 wells and temperature at 2,470 wells in the Ararat Basin using U.S. Geological Survey protocols in 2016. The minimum specific conductance was 377 microsiemens per centimeter (μS/cm), the maximum value was 4,000 μS/cm, and the mean was 998 μS/cm. The maximum water temperature was 24.2 degrees Celsius. An analysis between water temperature and well depth indicated no relation; however, spatially, most wells with cooler water temperatures were within the 2016 pressure boundary or in the western part of the basin. Wells with generally warmer water temperatures were in the eastern part of the basin. Samples were collected from four groundwater sites and one surface-water site by the U.S. Geological Survey in 2016. The stable-isotope values were similar for all five sites, indicating similar recharge sources for the sampled wells. The Hrazdan River sample was consistent with the groundwater samples, indicating the river could serve as a source of recharge to the Ararat artesian aquifer.

Ararat Basin↗
USGS2006

Analysis of the Magnitude and Frequency of Peak Discharges for the Navajo Nation in Arizona, Utah, Colorado, and New Mexico

Estimates of the magnitude and frequency of peak discharges are necessary for the reliable flood-hazard mapping in the Navajo Nation in Arizona, Utah, Colorado, and New Mexico. The Bureau of Indian Affairs, U.S. Army Corps of Engineers, and Navajo Nation requested that the U.S. Geological Survey update estimates of peak discharge magnitude for gaging stations in the region and update regional equations for estimation of peak discharge and frequency at ungaged sites. Equations were developed for estimating the magnitude of peak discharges for recurrence intervals of 2, 5, 10, 25, 50, 100, and 500 years at ungaged sites using data collected through 1999 at 146 gaging stations, an additional 13 years of peak-discharge data since a 1997 investigation, which used gaging-station data through 1986. The equations for estimation of peak discharges at ungaged sites were developed for flood regions 8, 11, high elevation, and 6 and are delineated on the basis of the hydrologic codes from the 1997 investigation. Peak discharges for selected recurrence intervals were determined at gaging stations by fitting observed data to a log-Pearson Type III distribution with adjustments for a low-discharge threshold and a zero skew coefficient. A low-discharge threshold was applied to frequency analysis of 82 of the 146 gaging stations. This application provides an improved fit of the log-Pearson Type III frequency distribution. Use of the low-discharge threshold generally eliminated the peak discharge having a recurrence interval of less than 1.4 years in the probability-density function. Within each region, logarithms of the peak discharges for selected recurrence intervals were related to logarithms of basin and climatic characteristics using stepwise ordinary least-squares regression techniques for exploratory data analysis. Generalized least-squares regression techniques, an improved regression procedure that accounts for time and spatial sampling errors, then was applied to the same data used in the ordinary least-squares regression analyses. The average standard error of prediction for a peak discharge have a recurrence interval of 100-years for region 8 was 53 percent (average) for the 100-year flood. The average standard of prediction, which includes average sampling error and average standard error of regression, ranged from 45 to 83 percent for the 100-year flood. Estimated standard error of prediction for a hybrid method for region 11 was large in the 1997 investigation. No distinction of floods produced from a high-elevation region was presented in the 1997 investigation. Overall, the equations based on generalized least-squares regression techniques are considered to be more reliable than those in the 1997 report because of the increased length of record and improved GIS method. Techniques for transferring flood-frequency relations to ungaged sites on the same stream can be estimated at an ungaged site by a direct application of the regional regression equation or at an ungaged site on a stream that has a gaging station upstream or downstream by using the drainage-area ratio and the drainage-area exponent from the regional regression equation of the respective region.

Scientific Investigations Report↗
USGS2012

Preliminary analysis of the hydrologic effects of temporary shutdowns of the Rondout-West Branch Water Tunnel on the groundwater-flow system in Wawarsing, New York

Flooding of streets and residential basements, and bacterial contamination of private-supply wells with Escherichia coli (E. coli) are recurring problems in the Rondout Valley near the Town of Wawarsing, Ulster County, New York. Leakage from the Rondout-West Branch (RWB) Water Tunnel and above-normal precipitation have been suspected of causing elevated groundwater levels and basement flooding. The hydrology of a 12-square-mile study area within the Town of Wawarsing was studied during 2008-10. A network of 41 wells (23 unconsolidated-aquifer and 18 bedrock wells) and 2 surface-water sites was used to monitor the hydrologic effects of four RWB Water Tunnel shutdowns. The study area is underlain by a sequence of northeast-trending sedimentary rocks that include limestone, shale, and sandstone. The bedrock contains dissolution features, fractures, and faults. Inflows that ranged from less than 1 to more than 9,000 gallons per minute from the fractured bedrock were documented during construction of the 13.5-foot-diameter RWB Water Tunnel through the sedimentary-rock sequence 710 feet (ft) beneath the study-area valley. Glacial sediments infill the valley above the bedrock sequence and consist of clay, silt, sand, and gravel. The groundwater-flow system in the valley consists of both fractured-rock and unconsolidated aquifers. Water levels in both the bedrock and unconsolidated aquifers respond to variations in seasonal precipitation. During the past 9 years (2002-10), annual precipitation at Central Park, N.Y., has exceeded the 141-year mean. Potentiometric-surface maps indicate that groundwater in the bedrock flows from the surrounding hills on the east and west sides of the valley toward the center of the valley, and ultimately toward the northeast. On average, water levels in the bedrock aquifer had seasonal differences of 5.3 ft. Analysis of hydrographs of bedrock wells indicates that many of these wells are affected by the RWB Tunnel leakage. Tunnel-leakage influences (water level and temperature changes) in the bedrock aquifer were measured at distances up to 7,000 ft from the RWB Tunnel. Water levels in the bedrock changed as much as 12 ft within 0.5 hour during tunnel shutdowns. Nine of the 10 wells that responded to the shutdowns showed a water-level response of 5 ft or greater. Changes in water levels ranged from 1.5 to 12 ft, with tunnel-leakage influence delay times ranging from 0.5 to 60 hours. Many of the longest tunnel-influence delay times and smallest water-level changes were in wells located closest to the tunnel in shale. Tunnel-influence response of the bedrock aquifer is consistent with its preliminary characterization as an anisotropic aquifer with greater transmissivity along bedding strike than across bedding strike. This tunnel-influence response is also consistent with the likely presence of discrete high-transmissivity networks along fractured limestone beds that have undergone dissolution. A lack of bedrock observation wells in half of the study area hampered a more thorough analysis of the extent of leakage from the RWB Tunnel in the study area. On average, water levels in the unconsolidated aquifer had a seasonal difference of 5.0 ft. Some unconsolidated-aquifer wells indicated water-level changes due to tunnel leakage. The locations of unconsolidated-aquifer wells with measurable water-level changes due to tunnel leakage correlated with those in the bedrock. Water levels in the unconsolidated aquifer changed as much as 2.5 ft within 18 hours of tunnel shutdowns, but water-level changes in some unconsolidated-aquifer wells were smaller or nonexistent.

New York↗
USGS2014

Status and trends of Caribbean coral reefs: 1970-2012

This it the 9th status report since the Global Coral Reef Monitoring Network (GCRMN) was founded in 1995 was the data arm of the International Coral Reef Initiative (ICRI) to document the ecological condition or corral reefs, strengthen monitoring efforts, and link existing organizations and people working on reefs worldwide. The US Government provided the initial funding to help set up a global network of coral reef workers and has continued to provide core support. Since then, the series of reports have aimed to present the current status of coral reefs of the world or particular regions, the major threats to reefs and their consequences, and any initiative undertaken under the auspices of ICRI or other bodies to arrest or reverse the decline of coral reefs. IUCN assumed responsibility for hosting the global coordination of the GCRMN in 2010 under the scientific direction of Jeremy Jackson with the following objectives: 1. Document quantitatively the global status and trends for corals, macroalgae, sea urchins, and fishes based on available data from individual scientists as well as the peer reviewed scientific literature, monitoring programs, and report. 2. Bring together regional experts in a series of workshops to involve them in data compilation, analysis, and synthesis. 3. Integrate coral reef status and trends with independent environmental, management, and socioeconomic data to better understand the primary factors responsible for coral reef decline, the possible synergies among factors that may further magnify their impacts, and how these stresses may be more effectively alleviated. Work with GCRMN partners to establish simple and practical standardized protocols for future monitoring and assessment. Disseminate information and results to help guide member state policy and actions. The overarching objective is to understand why some reefs are much healthier than others, to identify what kinds of actions have been particularly beneficial or harmful, and to vigorously communicate results in simple and straightforward terms to foster more effective conservation and management. This and subsequent reports will focus on separate biogeographic regions in a stepwise fashion and combine all of the results for a global synthesis in the coming years. We began in the wide Caribbean region because the historical data are so extensive and to refine methods of analysis before moving on to other regions. This report documents quantitative trends for Caribbean reef corals, macroalgae, sea urchins, and fishes based on data from 90 reef locations over the past 43 tears. This is the first report to combine all these disparate kinds of data in a single place to explore how the different major components of coral reef ecosystems interact on a broadly regional oceanic scale. We obtained data from more than 35,000 ecological surveys carried out by 78 principal investigators (PIs) and some 200 colleagues working in 34 countries, states, and territories throughout the wide Caribbean region. We conducted two workshops in Panama and Brisbane, Australia to bring together people who provided the data to assist in data quality control, analysis, and synthesis. The first workshop at the Smithsonian Tropical Research Institute (STRI) in the Republic of Panama 29 April to 5 May, 2012 included scientists from 18 countries and territories to verify and expand the database and to conduct exploratory analyses of status and trends. Preliminary results based on the Panama workshop were presented to the DC Marine Community and Smithsonian Institution Senate of Scientists in May 2012 and at the International Coral Reef Symposium (ICRS) and annual ICRI meeting in Cairns, Australia in July 2012. The second workshop in Brisbane, Australia in December 2012 brought together eight coral reef scientists for more detailed data analysis and organization of results for this report and subsequent publications. Subsequent presentations to solicit comments while the report was being finalized were made in 2013-2014 at the ICRI General meeting in Belize, the biennial meeting of the Association of Island Marine Laboratories in Jamaica, the Panamerican Coral Reef Congress in Merida, Mexico, the annual meeting of he Western Society of Naturalists, and numerous universities in Costa Rica, the USA and Europe. The main body of the report is in two sections. Part I provides an overview of overall status and trends and detailed analyses of the multiple factors responsible for the decline of reef corals throughout the entire wider Caribbean region. The editors are grateful to all the people who have so generously provided data and expertise, but we assume responsibility for the many statements, conclusions and recommendations and final wording of the text. Part II provides a more detailed analysis of the status and trends of coral reef ecosystems in the 32 countries, states, and territories for which we have data. The format includes maps indicating all locations sampled, a detailed table of data sources and sites surveyed, timelines of ecologically important evens, and relevant references. Each of these reports was compiled in consultation with local experts and all those who provided data and advice are listed as authors of each country report.

Report↗
USGS2013

Results of repeat bathymetric and velocimetric surveys at the Amelia Earhart Bridge on U.S. Highway 59 over the Missouri River at Atchison, Kansas, 2009-2013

Bathymetric and velocimetric data were collected six times by the U.S. Geological Survey, in cooperation with the Kansas Department of Transportation, in the vicinity of Amelia Earhart Bridge on U.S. Highway 59 over the Missouri River at Atchison, Kansas. A multibeam echosounder mapping system and an acoustic Doppler current meter were used to obtain channel-bed elevations and depth-averaged velocities for a river reach approximately 2,300 feet long and extending across the active channel of the Missouri River. The bathymetric and velocimetric surveys provide a “snapshot” of the channel conditions at the time of each survey, and document changes to the channel-bed elevations and velocities during the course of construction of a new bridge for U.S. Highway 59 downstream from the Amelia Earhart Bridge. The baseline survey in June 2009 revealed substantial scour holes existed at the railroad bridge piers upstream from and at pier 10 of the Amelia Earhart Bridge, with mostly uniform flow and velocities throughout the study reach. After the construction of a trestle and cofferdam on the left (eastern) bank downstream from the Amelia Earhart Bridge, a survey on June 2, 2010, revealed scour holes with similar size and shape as the baseline for similar flow conditions, with slightly higher velocities and a more substantial contraction of flow near the bridges than the baseline. Subsequent surveys during flooding conditions in June 2010 and July 2011 revealed substantial scour near the bridges compared to the baseline survey caused by the contraction of flow; however, the larger flood in July 2011 resulted in less scour than in June 2010, partly because the removal of the cofferdam for pier 5 of the new bridge in March 2011 diminished the contraction near the bridges. Generally, the downstream part of the study reach exhibited varying amounts of scour in all of the surveys except the last when compared to the baseline. During the final survey, velocities throughout the study area were the lowest of all the surveys, resulting in overall deposition throughout the reach compared to the baseline survey—despite the presence of the trestle in the final survey. The multiple surveys at the Amelia Earhart Bridge document the effects of moderate- to high-flow conditions on scour, compounded by the effects of adding and removing a constriction in the channel. Additional factors such as pier shape and angle of approach flow also were documented.

Kansas;Missouri↗
USGS2008

Mobility statistics and automated hazard mapping for debris flows and rock avalanches

Power-law equations that are physically motivated and statistically tested and calibrated provide a basis for forecasting areas likely to be inundated by debris flows, rock avalanches, and lahars with diverse volumes. The equations A=α 1 V 2/3 and B=α 2 V 2/3 are based on the postulate that the maximum valley cross-sectional area (A) and total valley planimetric area (B) likely to be inundated by a flow depend only on its volume (V) and the topography of the flow path. Testing of these equations involves determining whether or not they fit data for documented flows satisfactorily, and calibration entails determining best-fit values of the coefficients α 1 and α 2 for debris flows, rock avalanches, and lahars. This report describes statistical testing and calibration of the equations by using field data compiled from many sources, and it describes application of the equations to delineation of debris-flow hazard zones. Statistical results show that for each type of flow (debris flows, rock avalanches, and lahars), the dependence of A and B on V is described well by power laws with exponents equal to 2/3. This value of the exponent produces fits that are effectively indistinguishable from the best fits obtained by using adjustable power-law exponents. Statistically calibrated values of the coefficients α 1 and α 2 provide scale-invariant indices of the relative mobilities of rock avalanches (α 1 = 0.2, α 2 = 20), nonvolcanic debris flows (α 1 = 0.1, α 2 = 20), and lahars (α 1 = 0.05, α 2 = 200). These values show, for example, that a lahar of specified volume can be expected to inundate a planimetric area ten times larger than that inundated by a rock avalanche or nonvolcanic debris flow of the same volume. The utility of the calibrated debris-flow inundation equations A=0.1V 2/3 and B=20V 2/3 is demonstrated by using them within the GIS program LAHARZ to delineate nested hazard zones for future debris flows in an area bordering the Umpqua River in the south-central Oregon Coast Range. This application requires use of high-resolution topographic data derived form LIDAR surveys, knowledge of local geology to specify a suitable range of prospective debris-flow volumes, and development and use of a new algorithm for identification of prospective debris-flow source areas in finely dissected terrain.

Oregon↗
USGS2020

Optimization assessment of a groundwater-level observation network in the Middle Rio Grande Basin, New Mexico

The U.S. Geological Survey, in cooperation with the Albuquerque Bernalillo County Water Utility Authority (ABCWUA), measures groundwater levels continuously (hourly) and discretely (semiannually and annually) at a network of wells and piezometers (hereafter called the observation network) within the Middle Rio Grande Basin in central New Mexico. Groundwater levels that are measured in this observation network provide a long-term hydrologic dataset that is heavily relied upon to make water management decisions. The desire to upgrade and perform maintenance on this observation network initiated this study, which assesses the spatial and temporal importance of measurements towards optimization of monitoring the observation network to reduce or redirect monitoring costs. This report describes the methods and results of the optimization assessment of this observation network, which included separate spatial and temporal methodologies and an evaluation using principal component analysis (PCA). Results from the spatial optimization assessment can be used to help identify observation network sites that do not significantly affect the generation of winter groundwater-elevation contour maps of the production zone. Results from the temporal optimization assessment and PCA can also be consulted when deciding which sites to remove from the network, especially for sites that are monitored more frequently than annually. Results from the temporal optimization assessment can be used to inform the minimum monitoring frequency at the observation network required to capture the trends shown in higher frequency monitoring. The PCA results distinguish spatially distributed characteristic water-level trends that can inform the management decisions that are made when using the spatial and temporal optimization assessment results. Reducing the temporal frequency or spatial density of monitoring is ultimately a management decision that depends on the amount of data loss or degradation that is deemed acceptable while still meeting the network objectives of the ABCWUA. This study can also serve as a starting point to a data gap analysis of local aquifer characteristics and help guide enhanced observation network design as needs arise or in advance of future water management decisions. The results of the spatial optimization assessment indicate that as many as about 20 specified sites can be removed from the observation network with a relatively small loss in the ability to represent the kriged groundwater-elevation surfaces of the production zone that were generated by using median groundwater elevations for two periods: the 2001 time interval and 2015 time interval. This analysis also demonstrated the importance of wells at the margin of the study area and in areas where there are large hydrologic gradients. At many of the 47 hourly monitored sites analyzed in the temporal optimization assessment, temporal trends were well represented for at least one of the reduced monitoring frequencies tested, indicating that a reduced frequency may be sufficient to adequately characterize seasonal and long-term trends. PCA and k-means clustering analysis of the 15 hourly monitored sites that are screened within the production zone indicate that the sites can be categorized into four groups, or clusters, of differing groundwater-level hydrograph characteristics. Except for one cluster, all of the clusters have the potential to be well represented by fewer index monitoring sites.

New Mexico↗
USGS2018

Effects of surface-water and groundwater inflows and outflows on the hydrology of the Tsala Apopka Lake Basin in Citrus County, Florida

The U.S. Geological Survey, in cooperation with the Southwest Florida Water Management District, initiated a study to quantify the inflows and outflows in the Floral City, Inverness, and Hernando pools of the Tsala Apopka Lake Basin in Citrus County, Florida. This study assesses hydrologic changes in pool stages, groundwater levels, spring flows, and streamflows caused by the diversion of streamflow from the Withlacoochee River to the Tsala Apopka Lake Basin through water-control structures. A surface-water/groundwater flow model was developed using hydraulic parameters for lakes, streams, the unsaturated zone, and the underlying surficial and Upper Floridan aquifers estimated using an inverse modeling calibration technique. After calibration, the model was used to assess the relation between inflows and outflows in the Tsala Apopka Lake Basin and changes in pool stages. Simulation results using the calibrated surface-water/groundwater flow model showed that leakage rates from the pools to the Upper Floridan aquifer were largest at the deep lake cells and that these leakage rates to the Upper Floridan aquifer were the highest in the model area. Downward leakage to the Upper Floridan aquifer occurred beneath most of the extent of the Floral City, Inverness, and Hernando pools. These leakage rates depended on the lakebed leakance and the difference between lake stages and heads in the Upper Floridan aquifer. Leakage rates were higher for the Floral City pool than for the Inverness pool, and higher for the Inverness pool than for the Hernando pool. Lakebed leakance was higher for the Floral City pool than for the Hernando pool, and higher for the Hernando pool than for the Inverness pool. Simulation results showed that the average recharge rate to the surficial aquifer was 10.3 inches per year for the 2004 to 2012 simulation period. Areas that recharge the surficial aquifer covered about 86 percent of the model area. Simulations identified areas along segments of the Withlacoochee River and within land-surface depressions that receive water from the surficial aquifer. Recharge rates were largest in physiographic regions having a deep water table. Simulated heads in the Upper Floridan aquifer indicated the general flow directions in the active flow model area were from the northeast toward the southwest and then westward toward the coast, and from the southeast toward the northwest and then westward toward the coast, consistent with flow directions inferred from the estimated potentiometric surface map for May 2010. The largest inflow in the water budget of the Upper Floridan aquifer was downward leakage from the overlying hydrogeologic unit. The largest outflow in the water budget of the Upper Floridan aquifer was spring flow. The calibrated surface-water and groundwater flow model was used to simulate hydrologic scenarios that included changes in rainfall rates, projected increases in groundwater pumping rates for 2025 and 2035, no flow for the 2004–12 period through the eight water-control structures in the Tsala Apopka Lake Basin, and the removal of the Inglis Dam and the Inglis Bypass Spillway on Lake Rousseau. Scenario simulation results were compared to annual average calibrated water levels and flows from 2004 to 2012. Simulated declines in the Tsala Apopka Lake pool stages under the 10-percent lower rainfall scenario were about 0.8, 0.3, and 1.3 feet (ft) for the Floral City, Inverness, and Hernando pools, respectively. Simulated groundwater levels under the same scenario declined up to 5.4 ft in the surficial aquifer and up to 2.9 ft in the Upper Floridan aquifer. Under the projected increases in groundwater pumping rates for 2035 that represented an increase of 36 percent from average 2004 to 2012 pumping rates, the simulated declines in the Floral City, Inverness, and Hernando pool stages were, in downstream order, 0.02, 0.06, and 0.04 ft. The largest drawdown under the projected increases in groundwater pumping rates for 2035 was 2.1 ft in the surficial aquifer and about 1.8 ft in the Upper Floridan aquifer. A scenario of decreased rainfall by 10 percent caused greater declines in water levels and pool stages than projected increases in groundwater pumping rates. The simulation with no flow through the eight Tsala Apopka Lake water-control structures resulted in simulated declines in average pool stage of 1.8, 1.9, and 0.5 ft in the Floral City, Inverness, and Hernando pools, respectively. The simulated removal of the two water-control structures in Lake Rousseau caused flow to increase at Rainbow Springs by 28 cubic feet per second, an increase of 4.7 percent from the average calibrated flow for 2004 to 2012.

Florida↗
USGS2021

Sediment budget for watersheds of West Maui, Hawaii

Episodic runoff brings suspended sediment to West Maui’s nearshore waters, turning them from blue to brown. This pollution degrades the ecological, cultural, and recreational value of these iconic nearshore waters. We used mapping, monitoring, and modeling to identify and quantify the watershed sources for fine sediment that pollutes the nearshore each year. These results focus strategies to reduce pollution on the outstanding sources for this sediment. Terrestrial runoff causing coastal plumes now occurs when two or more hours of rain falls at rates greater than 10–20 millimeter (mm) per hour in source watersheds. Analysis of recent and historical rainfall indicates that West Maui communities can expect rainfalls to bring coastal plumes at least 3–5 times per year. Former agricultural fields and some unimproved roads are possible sources for the fine sediment of these plumes. We found, however, that these obvious sources do not produce plumes during small annual storms, because they drain water at rates that far exceed most annual rainfalls and because there is no evidence for runoff from rains that caused recent plumes. Streambanks now eroding into historic fill terraces of sands, silts, and clays are a more plausible source. These terraces are found only downstream of historical agricultural fields and are composed of silt and fine sand. Surveys show that the fill terraces occupy ~40 percent of streambank length, making them extensive. During 2015–2016, these deposits eroded at median rates of 5–24 mm per year. Summed over West Maui’s watersheds, these rates imply sediment loads carried to the coast that can be ten times or more than modeled pre-human values. A sediment budget indicates that bank erosion of fill terraces from a few watersheds likely dominates the current annual fine-sediment load to the nearshore, with Kahana Stream watershed producing the largest annual input of 285 metric tons, the equivalent of 29 dump-truck loads every year. Although past large storms have contributed to sediment loading, annual plume generation is now caused by smaller rainfalls eroding these near-stream terrace deposits, a legacy of historic agriculture. Treatments of former agricultural fields, roads, and reserve forests are consequently not likely to measurably effect sediment pollution from smaller, more frequent storms. Increased runoff from residential and commercial development of West Maui has the potential to exacerbate sediment plumes from such storms.

Hawaii↗
USGS2009

Estimation of Leakage Potential of Selected Sites in Interstate and Tri-State Canals Using Geostatistical Analysis of Selected Capacitively Coupled Resistivity Profiles, Western Nebraska, 2004

With increasing demands for reliable water supplies and availability estimates, groundwater flow models often are developed to enhance understanding of surface-water and groundwater systems. Specific hydraulic variables must be known or calibrated for the groundwater-flow model to accurately simulate current or future conditions. Surface geophysical surveys, along with selected test-hole information, can provide an integrated framework for quantifying hydrogeologic conditions within a defined area. In 2004, the U.S. Geological Survey, in cooperation with the North Platte Natural Resources District, performed a surface geophysical survey using a capacitively coupled resistivity technique to map the lithology within the top 8 meters of the near-surface for 110 kilometers of the Interstate and Tri-State Canals in western Nebraska and eastern Wyoming. Assuming that leakage between the surface-water and groundwater systems is affected primarily by the sediment directly underlying the canal bed, leakage potential was estimated from the simple vertical mean of inverse-model resistivity values for depth levels with geometrically increasing layer thickness with depth which resulted in mean-resistivity values biased towards the surface. This method generally produced reliable results, but an improved analysis method was needed to account for situations where confining units, composed of less permeable material, underlie units with greater permeability. In this report, prepared by the U.S. Geological Survey in cooperation with the North Platte Natural Resources District, the authors use geostatistical analysis to develop the minimum-unadjusted method to compute a relative leakage potential based on the minimum resistivity value in a vertical column of the resistivity model. The minimum-unadjusted method considers the effects of homogeneous confining units. The minimum-adjusted method also is developed to incorporate the effect of local lithologic heterogeneity on water transmission. Seven sites with differing geologic contexts were selected following review of the capacitively coupled resistivity data collected in 2004. A reevaluation of these sites using the mean, minimum-unadjusted, and minimum-adjusted methods was performed to compare the different approaches for estimating leakage potential. Five of the seven sites contained underlying confining units, for which the minimum-unadjusted and minimum-adjusted methods accounted for the confining-unit effect. Estimates of overall leakage potential were lower for the minimum-unadjusted and minimum-adjusted methods than those estimated by the mean method. For most sites, the local heterogeneity adjustment procedure of the minimum-adjusted method resulted in slightly larger overall leakage-potential estimates. In contrast to the mean method, the two minimum-based methods allowed the least permeable areas to control the overall vertical permeability of the subsurface. The minimum-adjusted method refined leakage-potential estimation by additionally including local lithologic heterogeneity effects.

Scientific Investigations Report↗
USGS2020

Effects of huisache removal on rangeland evapotranspiration in Victoria County, south-central Texas, 2015–18

The U.S. Geological Survey and Desert Research Institute, in cooperation with the Natural Resources Conservation Service, Texas State Soil and Water Conservation Board, Victoria County Groundwater Conservation District, Victoria Soil and Water Conservation District, and the San Antonio River Authority, evaluated the hydrologic effects of Vachellia farnesiana var. farnesiana (huisache) removal on rangeland evapotranspiration in Victoria County, Texas. Measurements of evapotranspiration, rainfall, and related properties were made at two sites during March 2015 through August 2018. One site was predominantly grassland. The other site was dominated by dense huisache vegetation that was removed about halfway through the study period. The resulting evapotranspiration data were examined for differences between the locations and differences between the pre-removal (2015–16) and post-removal (2017–18) periods to assess the effects of huisache removal on evapotranspiration. Evapotranspiration measurements were made using the eddy-covariance technique and were supplemented by remote-sensing estimates of evapotranspiration derived from thermal and optical satellite images. A map of remotely sensed evapotranspiration was generated for the area surrounding the study sites for 2015 and demonstrates the capability of remote sensing to evaluate land-management effects on evapotranspiration for larger scale areas, such as a county or stream-basin area. During the pre-removal period (March 2015–December 2016), evapotranspiration was greater at the huisache site than at the grassland site. Evapotranspiration at the grassland site (average of the eddy-covariance evapotranspiration and average remotely sensed evapotranspiration) was 87.6 millimeters per month (mm/mo) and at the huisache site was 100.8 mm/mo, with the differences in evapotranspiration rates being attributed to the difference in site vegetation. After huisache was removed in January 2017, evapotranspiration at the huisache site was substantially lower than at the grassland site, the changes in evapotranspiration rates being attributed not only to removal of huisache vegetation but also to possible disruption of soil runoff and infiltration characteristics. During the post-removal period (February 2017–August 2018), evapotranspiration was 88.5 mm/mo at the grassland site and 72.9 mm/mo at the huisache site (average of the eddy-covariance and average remotely sensed evapotranspiration). The monthly differences in evapotranspiration between the grassland and huisache sites, determined by eddy-covariance and remote-sensing methods, were statistically significant between the pre-removal and post-removal periods. Also, the pre-removal period provided the best conditions to evaluate the differences between huisache site and grassland site evapotranspiration. During the pre-removal period, evapotranspiration from the huisache site as measured by the eddy-covariance method was, on average, 10.7 mm/mo greater than evapotranspiration measured at the grassland site. As determined by the average of the remotely sensed methods, huisache site evapotranspiration was 15.8 mm/mo greater than grassland site evapotranspiration. These average differences in evapotranspiration rates by the two methods indicate that evapotranspiration at the grassland site was, on average, 13.2 mm/mo less than that at the huisache site during the pre-removal period. This average difference in evapotranspiration rates also indicates potential increased groundwater recharge and (or) surface-water runoff at the grassland site.

Texas↗
USGS2011

Water availability and use pilot: Methods development for a regional assessment of groundwater availability, southwest alluvial basins, Arizona

Executive Summary: Arizona is located in an arid to semiarid region in the southwestern United States and is one of the fastest growing States in the country. Population in Arizona surpassed 6.5 million people in 2008, an increase of 140 percent since 1980, when the last regional U.S. Geological Survey (USGS) groundwater study was done as part of the Regional Aquifer System Analysis (RASA) program. The alluvial basins of Arizona are part of the Basin and Range Physiographic Province and cover more than 73,000 mi 2 , 65 percent of the State's total land area. More than 85 percent of the State's population resides within this area, accounting for more than 95 percent of the State's groundwater use. Groundwater supplies in the area are expected to undergo further stress as an increasing population vies with the State's important agricultural sector for access to these limited resources. To provide updated information to stakeholders addressing issues surrounding limited groundwater supplies and projected increases in groundwater use, the USGS Groundwater Resources Program instituted the Southwest Alluvial Basins Groundwater Availability and Use Pilot Program to evaluate the availability of groundwater resources in the alluvial basins of Arizona. The principal products of this evaluation of groundwater resources are updated groundwater budget information for the study area and a proof-of-concept groundwater-flow model incorporating several interconnected groundwater basins. This effort builds on previous research on the assessment and mapping of groundwater conditions in the alluvial basins of Arizona, also supported by the USGS Groundwater Resources Program. Regional Groundwater Budget: The Southwest Alluvial Basins-Regional Aquifer System Analysis (SWAB-RASA) study produced semiquantitative groundwater budgets for each of the alluvial basins in the SWAB-RASA study area. The pilot program documented in this report developed new quantitative estimates of groundwater budget components using recent (2000-2007) data and methods of data analysis. Estimates of inflow components, including mountain-front recharge, incidental recharge from irrigation of agriculture, managed recharge from recharge facilities, interbasin underflow from upgradient basins, and streamflow losses, are quantified for recent time periods. Mountain-front recharge is the greatest inflow component to the groundwater system and was estimated using two methods: a basin characteristic model and new precipitation information used in a previously developed regression equation. Annual mountain-front recharge for the study area for 1940-2007 estimated by the two methods is 730,000 acre-ft for the basin characteristic model and 643,000 acre-ft for the regression equation, representing 1.5 percent and 1.3 percent of precipitation, respectively. Outflow components, including groundwater withdrawals, evapotranspiration, and interbasin flow to downgradient basins, are also presented for recent time periods. Groundwater withdrawals accounted for the largest share of the water budget, with nearly 2.4 million acre-ft per year withdrawn from the study area in recent years. Evapotranspiration from groundwater was estimated at nearly 1.3 million acre-ft per year for the study area using a newly developed method incorporating vegetation indices from satellite images and land cover information. For water-budget components with temporal variation that could be assessed from available data, estimates for intervening time periods since before development were also developed. An estimate of aquifer storage change, representing both gains to and losses from the groundwater system since before development, was derived for the most developed basins in the study area using available estimates of groundwater-level changes and storage coefficients. An overall storage loss of 74.5 million acre-ft was estimated for these basins within the study area. Demonstration

Arizona↗
USGS2007

Hydrogeologic characteristics of the St. Croix River basin, Minnesota and Wisconsin: Implications for the susceptibility of ground water to potential contamination

Population growth in the St. Croix River Basin in Minnesota and Wisconsin has intensified concerns of county resource managers and the National Park Service, which is charged with protecting the St. Croix National Scenic Riverway, about the potential for ground-water contamination in the basin. This report describes a previously developed method that was adapted to illustrate potential ground-water-contamination susceptibility in the St. Croix River Basin. The report also gives an estimate of ground-water-residence time and surface-water/ground-water interaction as related to natural attenuation and movement of contaminants in five tributary basins. A ground-water-contamination-susceptibility map was adapted from a state-wide map of Wisconsin to the St. Croix River Basin by use of well-driller construction records and regional maps of aquifer properties in Minnesota and Wisconsin. Measures of various subsurface properties were combined to generate a spatial index of susceptibility. The subjective index method developed for the State of Wisconsin by Schmidt (1987) 1 was not derived from analyses of water-quality data or physical processes. Nonetheless, it was adapted for this report to furnish a seamless map across state boundaries that would be familiar to many resource managers. Following this method, areas most susceptible to contamination appear to have coarse-grained sediments (sands or gravels) and shallow water tables or are underlain by carbonate-bedrock aquifers. The least susceptible areas appear to have fine-grained sediments and deep water tables. If an aquifer becomes contaminated, the ground-water-residence time can affect potential natural attenuation along the ground-water-flow path. Mean basin ground-water-residence times were computed for the Apple, Kettle, Kinnickinnic, Snake and Sunrise River Basins, which are tributary basins to the St. Croix Basin, by use of average aquifer properties of saturated thickness, porosity, and recharge rates. The Apple River Basin had the shortest mean ground-water-residence times (20–120 years), owing largely to the moderate saturated thickness and high recharge rate in the basin. The Kinnickinnic and Sunrise River Basins had the longest mean residence times (60–350 and 70–390 years, respectively) chiefly because of the relatively large saturated thickness of the basins. Owing to limitations of the residence-time calculations, actual ground-water-residence times will vary around the mean values within each basin and may range from days or weeks in karst carbonate aquifers to millennia in deep confined sandstone aquifers. Areas of relatively short residence time (less than the median residence time in each basin) were identified by use of ground-water-flow models for each of the five tributary basins. Results of simulations show that these areas, in which contaminants may have relatively less time for natural attenuation along the short flow paths, generally occur near streams and rivers where ground water discharges to the surface. Finally, the ground-water-flow models were used to simulate ground-water/surface-water interaction in the five tributary basins. Results of simulations show that some lakes and reservoirs leak surface water into the ground-water-flow system on their downgradient side, where the surface-water outflow has been restricted by a dam or a naturally constricted outlet. These locations are noteworthy because contaminated surface waters could potentially enter the ground-water-flow system at these locations.

Minnesota, Wisconsin↗