USGS Science⌕ Search

SEARCH · USGS Science

Results for “Investigator Annual Report”

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,387 records · Page 77Linked to original sources

Hydrologic assessment of the Edwin B. Forsythe National Wildlife Refuge

The Edwin B. Forsythe National Wildlife Refuge (hereafter Forsythe refuge or the refuge) is situated along the central New Jersey coast and provides a mixture of freshwater and saltwater habitats for numerous bird, wildlife, and plant species. Little data and information were previously available regarding the freshwater dynamics that support the refuge’s ecosystems. In cooperation with the U.S. Fish and Wildlife Service, the U.S. Geological Survey conducted an assessment of the hydrologic resources and processes in the refuge and surrounding areas to provide baseline information for evaluating restoration projects and future changes in the hydrologic system associated with climate change and other anthropogenic stressors. During spring 2015, water levels were measured at groundwater and surface-water sites in and near the Forsythe refuge. These water-level measurements, along with surface-water elevations obtained from digital elevation models, were used to construct water-table-elevation and depth-to-water maps of the refuge and surrounding areas. Water-table elevations in the refuge ranged from sea level to approximately 65 feet above sea level; in most of the refuge, the water-table elevation was within 3 feet of sea level. The water-table-elevation map indicates that the direction of shallow groundwater flow at the regional scale is generally from west to east (much of it from the northwest to the southeast), and groundwater moves downgradient from the uplands toward major groundwater discharge areas consisting of coastal streams and wetlands. The depth to water is estimated to be less than 2 feet for approximately 86 percent of the refuge, which coincides closely with the percentage of wetland area in the refuge. Depth to water in excess of 20 feet below land surface is limited to higher elevation areas of the refuge. Streamflow data collected at continuous-record streamgages and partial-record stations within the Mullica-Toms Basin were summarized. Hydrograph separation of streamflow data for eight streamgages (2004–13) reveals that base flow accounts for 68–94 percent of streamflow in basins upstream from the refuge. The high base-flow inputs underscore the importance of groundwater as a source of freshwater that supports both the streams that flow into the refuge and the hydroecology of the contributing basins. Mean annual flow typically ranged from 1.7 to 2.1 cubic feet per second per square mile at the streamgages (2004–13) and between 1.2 and 2.3 cubic feet per second per square mile at the partial-record stations (1965–2015) but was notably greater or lower than these ranges at several stations. Mean annual water budgets were estimated for multiple regions of the refuge for 2004–13 using data compiled from nearby meteorological stations and groundwater flows derived from previously calibrated groundwater-flow models. Precipitation, groundwater recharge, and evapotranspiration were estimated from available data; direct runoff was calculated as the residual component of the water balance. Groundwater recharge rates were greatest in the upland-dominated areas of the refuge with estimates of 14.4 to 18.9 inches per year, which are equivalent to 30 to 40 percent of precipitation. Groundwater recharge rates were nearly zero in the central coastal areas because these areas are major groundwater discharge zones, the water table is near land surface, the subsurface is close to saturation and cannot accept much recharge, and much of the area is underlain by thick marsh deposits likely with low permeability. Estimates of evapotranspiration varied from about 26 inches per year in the upland-dominated areas to more than 35 inches per year in the coastal wetlands, equivalent to 55–79 percent of mean annual precipitation, indicating that it is a major component of the hydrodynamics of the Forsythe refuge. On the basis of output from previously calibrated groundwater-flow models, nearly all of the groundwater exiting the surficial aquifer system in the central coastal areas of the refuge is discharged to wetlands, which highlights the importance of groundwater discharge in supporting the ecosystems of the Forsythe refuge. In the central coastal areas, horizontal flow contributes more than 90 percent of the groundwater flow to the surficial system, indicating that the upbasin areas are a substantial source of water that ultimately discharges to streams and wetlands in the refuge.

New Jersey↗

Results of a Two-Dimensional Hydrodynamic and Sediment-Transport Model to Predict the Effects of the Phased Construction and Operation of the Olmsted Locks and Dam on the Ohio River near Olmsted, Illinois

The Olmsted two-dimensional hydrodynamic and sediment-transport model was developed in cooperation with the U.S. Army Corps of Engineers, Louisville District. The model was used to estimate the effects that the phased-construction sequence and operation of the Olmsted Locks and Dam had on sediment-transport patterns in the 11.9-mile study reach (Ohio River miles 962.6 to 974.5), particularly over an area of endangered orange-footed pearly mussel (Plethobasus cooperianus) beds beginning approximately 2 miles downstream of the dam construction. A Resource Management Associates?2 (RMA-2) two-dimensional hydrodynamic model for the reach was calibrated to a middle-flow hydraulic survey (350,000 cubic feet per second) and verified with data collected during low- and high-flow hydraulic surveys (72,500 and 770,000 cubic feet per second, respectively). The calibration and validation process included matching water-surface elevations at the construction site and velocity profiles at 15 cross sections throughout the study reach. The sediment-transport aspect of the project was simulated with the Waterways Experiment Station's Sed2D model for a 6-year planned-construction period (construction-phase modeling) and a subsequent 3-year operational period (operational-phase modeling). The sediment-transport results from the construction and operational models both were compared to results of concurrent baseline simulations to determine the changes in erosional and depositional patterns induced by the dam construction and operation throughout the study reach and more importantly over the area of the endangered mussel beds. Simulation of the phased-in-the-wet Olmsted Locks and Dam construction and subsequent operation period resulted in a maximum additional deposition of approximately 2 feet over a localized region of the mussel beds when compared to the bed change simulated with baseline conditions (river conditions that included only the completed locks section). Most areas on the mussel beds experienced less than 0.5 feet of cumulative bed change between the baseline and construction phases during the nine annual hydrographs. The bed change over the 9 year Olmsted Locks and Dam simulation reveals a continuous downstream progression and deepening of the main channel and deposition along the right bank with limited lateral migration toward the more densely populated mussel-bed areas. The sensitivity of the mussels to sediment deposition is difficult to quantify; therefore, the effect of simulated deposition on the welfare of the mussels is uncertain. The model also will provide the U.S. Army Corps of Engineers a tool to predict the locations of high deposition in navigable sections, which can save engineers time and resources when monitoring the need for dredging operations.

Water-Resources Investigations Report↗

Digital simulation of the effects of urbanization on runoff in the upper Santa Ana Valley, California

The Stanford Watershed Model was used to simulate the effects of urbanization on the discharge from five drainage basins in the upper Santa Ana Valley, an area with an average annual precipitation of 15 inches. The drainage basins ranged in size from 3.72 to 83.4 square miles. Using the model, synthetic records of streamflow for each basin were generated to represent various degrees of urban development. Examination of the synthetic records indicated that urbanization has the following effects on streamflow in the area: Average annual runoff from a drainage basin with an effective impervious area of 10 percent of the drainage area is approximately 2 inches, and increases by 1 inch for each increase in effective impervious cover equal to 10 percent of the drainage area. About 30 percent of a fully urbanized area is effectively impervious. Urbanization can increase the magnitude of peak discharge and daily mean discharge with a recurrence interval of 2 years by a factor of three to six. Peak discharges and daily mean discharges that have recurrence intervals greater than a limiting value ranging from 50 to 200 years or more are little affected by urbanization.

California↗

The effect of discharge and water quality of the Alafia River, Hillsborough River, and the Tampa Bypass Canal on nutrient loading to Hillsborough Bay, Florida

Techniques to measure discharge and nutrient loads in the tidally affected portions of two major rivers tributary to Tampa Bay, the Alafia River and the Hillsborough River, were developed and tested. Discharge, water quality, and total phosphorus and total nitrogen loads for the period April 1, 1991, through March 31, 1992, were evaluated and compared with discharge,water quality, and loads at long-term, nontidal gages in the basins. Long-term discharge and water-quality characteristics at selected sites in the Alafia river and Hillsborough River basins were evaluated. A long-term, decreasing trend in annual-mean discharge was observed for discharges at the Alafia River, Sulphur Springs, and Hillsborough River. Low-flow and high-flow characteristics in the Alafia River and Hillsborough River have changed as well. The decreasing trend in the Alafia River discharges is not due to deficient rainfall but probably is due to decreased ground-water inflow to the river because of long-term declines in the potentiometric surface of the Upper Floridan aquifer. Daily-mean discharges at the mouth of the Alafia River were more variable than discharges at the long-term gage upstream. Daily-mean discharge near the mouth of the river was negative at times, indicating a net loss of water from the river. Daily-mean discharge from the Hillsborough River was minimal from Apil to May 1991, and from late September 1991 to March 1992. During these periods, discharge from Sulphur Springs was a major source of freshwater to the tidally affected reach of the river. Concentrations of total phosphorus and orthophosphorus in the Alafia River above Lithia Springs were the greatest in the 1960's and have generally declined since then. Total nitrogen concentrations have been declining since about 1981. However, increases in nitrate plus nitrite nitrogen concentrations are occurring in Lithia Springs, a second-magnitude spring that flows into the Alafia River. Specific conductance of water discharging from Sulphur Springs to the Hillsborough River has increased from about 124 to more than 2,000 microsiemens per centimeter since 1945. Water quality at the mouth of the Alafia River and Hillsborough River is the result of mixing of freshwater and estuarine water from Hillsborough Bay. Large daily variations in water quality occur at these site because of tidal currents, and vertical stratification of specific conductance is a common feature. Concentrations of phosphorus, nitrate plus nitrite nitrogen, organic carbon, and silica are inversely related to specific conductance at the mouth of the Alafia River. Constituent concentration and discharge data were used to compute loads during the study period. Average daily phosphorus loads were 2.4 tons per day at the mouth of the Alafia River; 0.35 ton per day at the mouth of the Hillsborough River; and 0.06 ton per day at the Tampa Bypass Canal. Average daily nitrogen loads were 1.7 tons per day at the mouth of the Alafia River; 0.86 ton per day at the mouth of the Hillsborough River; and 0.26 ton per day at the Tampa Bypass Canal. The greatest annual loads of phosphorus and nitrogen from the major tributaries to Hillsborough Bay are from the Alafia River, with the greatest loads at the river mouth. Total phosphorus load from the Alafia River was about 894 tons during April 1991 through March 1992, more than six times greater than phosphorus loads form the Hillsborough River. Annual total nitrogen load at the mouth of the Alafia River was about 630 tons, two times greater than at the mouth of the Hillsborough River and more than six times greater than loads from the Tampa Bypass Canal. Basinwide yields of total phosphorus during April 1991 through March 1992 were about 2 tons per square mile at the mouth of the Alafia river basin and were about 0.2 ton per square mile at the mouth of the Hillsborough River. Total nitrogen yield was about 1.5 tons per square mile at the mouth of the Alafia River and about 0.5 ton per s

Water-Resources Investigations Report↗

Estimation of annual agricultural pesticide use for counties of the conterminous United States, 1992-2009

A method was developed to calculate annual county level pesticide use for selected herbicides, insecticides, and fungicides applied to agricultural crops grown in the conterminous United States from 1992 through 2009. Pesticide-use data compiled by proprietary surveys of farm operations located within Crop Reporting Districts were used in conjunction with annual harvested-crop acreage reported by the U.S. Department of Agriculture National Agricultural Statistics Service (NASS) to calculate use rates per harvested crop acre, or an 'estimated pesticide use' (EPest) rate, for each crop by year. Pesticide-use data were not available for all Crop Reporting Districts and years. When data were unavailable for a Crop Reporting District in a particular year, EPest extrapolated rates were calculated from adjoining or nearby Crop Reporting Districts to ensure that pesticide use was estimated for all counties that reported harvested-crop acreage. EPest rates were applied to county harvested-crop acreage differently to obtain EPest-low and EPest-high estimates of pesticide-use for counties and states, with the exception of use estimates for California, which were taken from annual Department of Pesticide Regulation Pesticide Use Reports. Annual EPest-low and EPest-high use totals were compared with other published pesticide-use reports for selected pesticides, crops, and years. EPest-low and EPest-high national totals for five of seven herbicides were in close agreement with U.S. Environmental Protection Agency and National Pesticide Use Data estimates, but greater than most NASS national totals. A second set of analyses compared EPest and NASS annual state totals and state-by-crop totals for selected crops. Overall, EPest and NASS use totals were not significantly different for the majority of crop-stateyear combinations evaluated. Furthermore, comparisons of EPest and NASS use estimates for most pesticides had rank correlation coefficients greater than 0.75 and median relative errors of less than 15 percent. Of the 48 pesticide-by-crop combinations with 10 or more state-year combinations, 12 of the EPest-low and 17 of the EPest-high totals showed significant differences (p < 0.05) from NASS use estimates. The differences between EPest and NASS estimates did not follow consistent patterns related to particular crops, years, or states, and most correlation coefficients were greater than 0.75. EPest values from this study are suitable for making national, regional, and watershed assessments of annual pesticide use from 1992 to 2009. Although estimates are provided by county to facilitate estimation of watershed pesticide use for a wide variety of watersheds, there is a greater degree of uncertainty in individual county-level estimates when compared to Crop Reporting District or state-level estimates because (1) EPest crop-use rates were developed on the basis of pesticide use on harvested acres in multi-county areas (Crop Reporting Districts) and then allocated to county harvested cropland; (2) pesticide-by-crop use rates were not available for all Crop Reporting Districts in the conterminous United States, and extrapolation methods were used to estimate pesticide use for some counties; and (3) it is possible that surveyed pesticide-by-crop use rates do not reflect all agricultural use on all crops grown. The methods developed in this study also are applicable to other agricultural pesticides and years.

Scientific Investigations Report↗

Simulation of ground-water flow in the Albuquerque Basin, central New Mexico, 1901-1994, with projections to 2020

This report describes a three-dimensional finite-difference ground-water-flow model of the Santa Fe Group aquifer system in the Albuquerque Basin, which comprises the Santa Fe Group (late Oligocene to middle Pleistocene age) and overlying valley and basin-fill deposits (Pleistocene to Holocene age). The model is designed to be flexible and adaptive to new geologic and hydrologic information as it becomes available, by using a geographic information system as a data-base manager to interface with the model. The aquifer system was defined and quantified in the model consistent with the current (July 1994) understanding of the structural and geohydrologic framework of the basin. Rather than putting the model through a rigorous calibration process, discrepancies between simulated and measured responses in hydraulic head were taken to indicate that the understanding of a local part of the aquifer system was incomplete or incorrect. The model simulates ground-water flow over an area of about 2,400 square miles to a depth of 1,730 to about 2,020 feet below the water table with 244 rows, 178 columns, and 11 layers. Of the 477,752 cells in the model, 310,376 are active. The top four model layers approximate the 80- foot thickness of alluvium in the incised and refilled valley of the Rio Grande to provide detail of the effect of ground-water withdrawals on the surface-water system. Away from the valley, these four layers represent the interval within the Santa Fe Group aquifer system between the computed predevelopment water table and a level 80 feet below the grade of the Rio Grande. The simulations include initial conditions (steady-state), the 1901-1994 historical period, and four possible ground-water withdrawal scenarios from 1994 to 2020. The model indicates that for the year ending in March 1994, net surface-water loss in the basin resulting from the City of Albuquerque's ground-water withdrawal totaled about 53,000 acre-feet. The balance of the about 123,000 acre-feet of withdrawal came from aquifer storage depletion (about 67,800 acre-feet) and captured or salvaged evapotranspiration (about 2,500 acrefeet). In the four scenarios projected from 1994 to 2020, City of Albuquerque annual withdrawals ranged from about 98,700 to about 177,000 acre-feet by the year 2020. The range of resulting surface-water loss was from about 62,000 to about 77,000 acre-feet. The range of aquifer storage depletion was from about 33,400 to about 95,900 acre-feet. Captured evapotranspiration and drain-return flow remained nearly constant for all scenarios. From 1994 to 2020, maximum projected declines in hydraulic head in the primary water-production zone of the aquifer (model layer 9) for the four scenarios ranged from 55 to 164 feet east of the Rio Grande and from 91 to 258 feet west of the river. Average declines in a 383.7-square-mile area around Albuquerque ranged from 28 to 65 feet in the production zone for the same period.

New Mexico↗

Simulation of net infiltration and potential recharge using a distributed-parameter watershed model of the Death Valley region, Nevada and California

This report presents the development and application of the distributed-parameter watershed model, INFILv3, for estimating the temporal and spatial distribution of net infiltration and potential recharge in the Death Valley region, Nevada and California. The estimates of net infiltration quantify the downward drainage of water across the lower boundary of the root zone and are used to indicate potential recharge under variable climate conditions and drainage basin characteristics. Spatial variability in recharge in the Death Valley region likely is high owing to large differences in precipitation, potential evapotranspiration, bedrock permeability, soil thickness, vegetation characteristics, and contributions to recharge along active stream channels. The quantity and spatial distribution of recharge representing the effects of variable climatic conditions and drainage basin characteristics on recharge are needed to reduce uncertainty in modeling ground-water flow. The U.S. Geological Survey, in cooperation with the Department of Energy, developed a regional saturated-zone ground-water flow model of the Death Valley regional ground-water flow system to help evaluate the current hydrogeologic system and the potential effects of natural or human-induced changes. Although previous estimates of recharge have been made for most areas of the Death Valley region, including the area defined by the boundary of the Death Valley regional ground-water flow system, the uncertainty of these estimates is high, and the spatial and temporal variability of the recharge in these basins has not been quantified. To estimate the magnitude and distribution of potential recharge in response to variable climate and spatially varying drainage basin characteristics, the INFILv3 model uses a daily water-balance model of the root zone with a primarily deterministic representation of the processes controlling net infiltration and potential recharge. The daily water balance includes precipitation (as either rain or snow), snow accumulation, sublimation, snowmelt, infiltration into the root zone, evapotranspiration, drainage, water content change throughout the root-zone profile (represented as a 6-layered system), runoff (defined as excess rainfall and snowmelt) and surface water run-on (defined as runoff that is routed downstream), and net infiltration (simulated as drainage from the bottom root-zone layer). Potential evapotranspiration is simulated using an hourly solar radiation model to simulate daily net radiation, and daily evapotranspiration is simulated as an empirical function of root zone water content and potential evapotranspiration. The model uses daily climate records of precipitation and air temperature from a regionally distributed network of 132 climate stations and a spatially distributed representation of drainage basin characteristics defined by topography, geology, soils, and vegetation to simulate daily net infiltration at all locations, including stream channels with intermittent streamflow in response to runoff from rain and snowmelt. The temporal distribution of daily, monthly, and annual net infiltration can be used to evaluate the potential effect of future climatic conditions on potential recharge. The INFILv3 model inputs representing drainage basin characteristics were developed using a geographic information system (GIS) to define a set of spatially distributed input parameters uniquely assigned to each grid cell of the INFILv3 model grid. The model grid, which was defined by a digital elevation model (DEM) of the Death Valley region, consists of 1,252,418 model grid cells with a uniform grid cell dimension of 278.5 meters in the north-south and east-west directions. The elevation values from the DEM were used with monthly regression models developed from the daily climate data to estimate the spatial distribution of daily precipitation and air temperature. The elevation values were also used to simulate atmosp

Water-Resources Investigations Report↗

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗

Nutrient Loading and Algal Response in West Thompson Lake, Thompson, Connecticut, 2003-2005

Water quality and nutrient loads were characterized for parts of the Quinebaug River and West Thompson Lake in northeastern Connecticut during 2003 to 2005. The West Thompson Lake watershed is a mainly forested watershed that receives treated municipal wastewater from several point sources in Massachusetts. The lake is a flood-control reservoir formed in 1966 by impoundment of the Quinebaug River. Median concentrations of total phosphorus in two inflow (upstream) and one outflow (downstream) sampling stations on the Quinebaug River were higher than the nutrient criteria recommended by the U.S. Environmental Protection Agency (USEPA) for rivers and streams in aggregate Ecoregion XIV. In general, concentrations of total phosphorus in West Thompson Lake also were above the nutrient criteria recommended by USEPA for lakes and impoundments in aggregate Ecoregion XIV. The trophic status of West Thompson Lake has changed since 1995 from a hypereutrophic lake to a eutrophic lake; however, the lake still has large algal blooms. These blooms are predominated by blue-green algae, with chlorophyll-a concentrations of more than 30 micrograms per liter and algal cell counts as high as 73,000 cells/mL. Water samples collected during the summer of 2005 identified phosphorus as the primary limiting nutrient early in the season, but algal growth is probably co-limited by phosphorus and nitrogen later in the season. Lake-bottom sediments were collected from several areas throughout the lake and ranged in thickness from less than 1 foot (ft) to more than 3 ft. Concentrations of phosphorus in sediments differed throughout the lake; the highest values were found in the middle of the lake. Concentrations of total phosphorus also increased from an average 1,800 milligrams per kilogram (mg/kg) in the upper layers of sediment to more than 6,000 mg/kg at depth in the sediment. Annual, seasonal, and monthly loads and yields of nutrients were calculated for the three sampling locations on the Quinebaug River to develop a nutrient mass-balance model (budget) for West Thompson Lake. The average annual yields of total phosphorus during 2000 to 2005 were 115 pounds per square mile per year (lb/mi2/yr) at Quinebaug (inflow station), 116 lb/mi2/yr at Red Bridge Road (inflow station), and 97.9 lb/mi2/yr at West Thompson (outflow station). The 18-percent decrease in the average annual yield of total phosphorus between the inflow station at Red Bridge Road and the outlet of West Thompson Lake at West Thompson indicates that a significant part of the phosphorus load is retained in the lake. Annual yields of total phosphorus at Quinebaug have decreased significantly since the 1980s, from 362 lb/mi2/yr (for 1981-1990) to 115 lb/mi2/yr (1996-2005). The annual net export of phosphorus in West Thompson Lake during water years 2000 to 2005 ranged from -36 percent (2005) to 1 percent (2002) of the incoming load. Seasonal mass-balance data for total phosphorus during the summers of 2000 to 2003, when streamflow was at or lower than normal, indicated a net export of phosphorus that ranged from 3.4 percent (2003) to 30.7 percent (2002) of the incoming load. During the summer of 2004, however, streamflows were much higher than normal, and there was a negative export of phosphorus in West Thompson Lake of -3.9 percent. The annual net export of nitrogen in West Thompson Lake during water years 2000 to 2005 ranged from -5 percent (2002) to 4 percent (2001) of the incoming load. No clear pattern was evident to relate total nitrogen export to seasonal variables or runoff. Removal of phosphorus during the summer by wastewater-treatment plants (WWTPs) in Massachusetts reduces the concentration and load of total phosphorus entering West Thompson Lake in the summer; however, the large amount of phosphorus retained in the lake during the other seasons, in addition to the phosphorus stored in the lake-bottom sediments, may become available to fuel algal blooms in the lake

Scientific Investigations Report↗

Simulation of regional ground-water flow in the Upper Deschutes Basin, Oregon

This report describes a numerical model that simulates regional ground-water flow in the upper Deschutes Basin of central Oregon. Ground water and surface water are intimately connected in the upper Deschutes Basin and most of the flow of the Deschutes River is supplied by ground water. Because of this connection, ground-water pumping and reduction of artificial recharge by lining leaking irrigation canals can reduce the amount of ground water discharging to streams and, consequently, streamflow. The model described in this report is intended to help water-management agencies and the public evaluate how the regional ground-water system and streamflow will respond to ground-water pumping, canal lining, drought, and other stresses. Ground-water flow is simulated in the model by the finite-difference method using MODFLOW and MODFLOWP. The finite-difference grid consists of 8 layers, 127 rows, and 87 columns. All major streams and most principal tributaries in the upper Deschutes Basin are included. Ground-water recharge from precipitation was estimated using a daily water-balance approach. Artificial recharge from leaking irrigation canals and on-farm losses was estimated from diversion and delivery records, seepage studies, and crop data. Ground-water pumpage for irrigation and public water supplies, and evapotranspiration are also included in the model. The model was calibrated to mean annual (1993-95) steady-state conditions using parameter-estimation techniques employing nonlinear regression. Fourteen hydraulic-conductivity parameters and two vertical conductance parameters were determined using nonlinear regression. Final parameter values are all within expected ranges. The general shape and slope of the simulated water-table surface and overall hydraulic-head distribution match the geometry determined from field measurements. The fitted standard deviation for hydraulic head is about 76 feet. The general magnitude and distribution of ground-water discharge to streams is also well simulated throughout the model. Ground-water discharge to streams in the area of the confluence of the Deschutes, Crooked, and Metolius Rivers is closely matched. The model was also calibrated to transient conditions from 1978 to 1997 using traditional trial-and-error methods. Climatic cycles during this period provided an excellent regional hydrologic signal for calibration. Climate-driven water-level fluctuations are simulated with reasonable accuracy over most of the model area. The timing and magnitude of simulated water-level fluctuations caused by annual pulses of recharge from precipitation match those observed reasonably well, given the limitations of the time discretization in the model. Water-level fluctuations caused by annual canal leakage are simulated very well over most of the area where such fluctuations occur. The transient model also simulates the volumetric distribution and temporal variations in ground-water discharge reasonably well. The match between simulated and measured volume of and variations in ground-water discharge is, however, somewhat dependent on geographic scale. The rates of and variations in ground-water discharge are matched best at regional scales. Example simulations were made to demonstrate the utility of the model for evaluating the effects of ground-water pumping or canal lining. Pumping simulations show that pumped water comes largely from aquifer storage when pumping begins, but as the water table stabilizes, the pumping increasingly diminishes the discharge to streams and, hence, streamflow. The time it takes for pumping to affect streamflow varies spatially depending, in general, on the location of pumping relative to the discharge areas. Canal-lining simulations show similar effects.

Water-Resources Investigations Report↗

Estimated effects of pumping on groundwater storage and Walker River stream efficiencies in Smith and Mason Valleys, west-central Nevada

The Walker River originates in the Sierra Nevada Mountains and flows nearly 160 miles to its terminus at Walker Lake in west-central Nevada. The river provides a source of irrigation water for tens of thousands of acres of agricultural lands in California and Nevada and is the principal source of inflow to Walker Lake. Extraction of groundwater for agricultural use became prevalent in the late 1950s and early 1960s to supplement irrigation demands not met by surface-water diversions during times of drought. There is growing concern that continued groundwater withdrawals within the Walker River Basin are likely contributing to depleted streamflow of the Walker River and the long-term depletion of groundwater storage in the basin. This report documents changes in groundwater storage-volume and trends in Walker River stream efficiency, a measure of change in flow due to gaining or losing conditions, in the two largest agricultural valleys in the Walker River Basin, Smith and Mason Valleys, for a multi-decade period. Groundwater-level maps from previous studies were used for the beginning (1970) and middle (2006) points of this study. Groundwater levels measured from 1991–95 and 2016–20 were used to construct median groundwater-level maps that represented conditions in 1995 and 2020. Valley wide groundwater-level change was calculated by comparing groundwater-level maps for the periods 1970–95, 1996–2006, and 2007–20 and by observing the overall change from 1970 to 2020. Groundwater storage-volume change was calculated using groundwater-level change and previously defined specific yield values. Between 1970 and 2020, groundwater storage-volume declined 287,600 acre-feet in Smith Valley and 269,000 acre-feet in Mason Valley. Using groundwater storage-volume decline and annual groundwater pumpage rates, a maximum groundwater pumpage rate can be computed to support management of water resources. In Smith Valley, groundwater pumping in excess of 22,300 acre-feet per year would likely result in groundwater storage decline. In Mason Valley, groundwater pumping in excess of 75,200 acre-feet per year would likely result in groundwater storage decline. Stream efficiency was calculated using continuous streamflow data and monthly diversion volumes on two reaches: (1) the West Walker River in Smith Valley, from 1948 to 2020 and (2) the Walker River in Mason Valley, from 1958 to 2020. Stream efficiency during non-irrigation season in Smith and Mason Valleys declined at a statistically significant rate of 1.1 and 0.6 percent per year, respectively. Trends in stream efficiency corresponded to occurrence of prolonged drought, deviation from average annual streamflows, and total groundwater pumpage. Long-term declines in groundwater storage-volume and stream efficiency demonstrate that the alluvial aquifer system is becoming increasingly depleted, such that the river can no longer replenish groundwater storage while simultaneously balancing groundwater and surface-water withdrawals. The introduction of supplemental groundwater pumpage was intended to offset surface-water deficits during dry years; however, pumpage occurs even in years when average or above average streamflows meet surface-water demands. Reliance on supplemental groundwater pumpage has resulted in widespread groundwater storage-volume decline and decreased stream efficiency. With each successive drought cycle, the ability of Walker River to sustain streamflows and convey water downstream has diminished. Above average wet periods have a marginal and short-lived effect on rebounding the groundwater levels outside of the river corridor. Moreover, if the trend continues, each future drought cycle may further reduce groundwater supplies and that may further decrease streamflow reliability.

Nevada↗

Evapotranspiration rate measurements of vegetation typical of ground-water discharge areas in the Basin and Range carbonate-rock aquifer system, Nevada and Utah, September 2005-August 2006

Evapotranspiration was measured at six eddy-correlation sites for a 1-year period between September 1, 2005, and August 31, 2006. Five sites were in phreatophytic shrubland dominated by greasewood, and one site was in a grassland meadow. The measured annual evapotranspiration ranged from 10.02 to 12.77 inches at the shrubland sites and 26.94 inches at the grassland site. Evapotranspiration rates correlated to measured vegetation densities and to satellite-derived vegetation indexes. Evapotranspiration rates were greater at sites with denser vegetation. The primary water source supporting evapotranspiration was water derived from local precipitation at the shrubland sites, and ground water at the grassland site. Measured precipitation, ranging from 6.21 to 11.41 inches, was within 20 percent of the computed long-term annual mean. The amount of ground water consumed by phreatophytes depends primarily on local precipitation and vegetation density. The ground-water contribution to local evapotranspiration ranged from 6 to 38 percent of total evapotranspiration at the shrubland sites, and 70 percent of total evapotranspiration at the grassland site. Average depth to water ranged from 7.2 to 32.4 feet below land surface at the shrubland sites, and 3.9 feet at the grassland site. Water levels declined throughout the growing season and recovered during the non-growing season. Diurnal water-level fluctuations associated with evapotranspiration were evident at some sites but not at others.

Nevada, Utah↗

Physical and Vegetative Characteristics of a Newly Constructed Wetland and Modified Stream Reach, Tredyffrin Township, Chester County, Pennsylvania, 2000-2006

To compensate for authorized disturbance of naturally occurring wetlands and streams during roadway improvements to U.S. Highway 202 in Chester and Montgomery Counties, Pa., the Pennsylvania Department of Transportation (PennDOT) constructed 0.42 acre of emergent wetland and 0.94 acre of scrub-shrub/forested wetland and modified sections of a 1,600-foot reach of Valley Creek with woody riparian plantings and streambank-stabilization structures (including rock deflectors). In accordance with project permits and additional guidance issued by the U.S. Army Corps of Engineers, the U.S. Geological Survey (USGS), in cooperation with PennDOT, collected data from 2000 through 2006 to quantify changes in 1) the vegetation, soils, and extent of emergent and scrub-shrub/forested parts of the constructed wetland, 2) the profile, dimension, and substrate in the vicinity of rock deflectors placed at two locations within the modified stream reach, and 3) the woody vegetation within the planted riparian buffer. The data for this investigation were collected using an approach adapted from previous investigations so that technology and findings may be more easily transferred among projects with similar objectives. Areal cover by planted and non-planted vegetation growing within the emergent and scrub-shrub/forested parts of the constructed wetland exceeded 85 percent at the end of each growing season, a criterion in special condition 25c in the U.S. Army Corps of Engineers project permit. Areal cover of vegetation in emergent and scrub-shrub/forested parts of the constructed wetland exceeded 100 percent in all but one growing season. Frequent and long-lasting soil saturation favored obligate-wetland species like Typha latifolia (broadleaf cattail) and Scirpus validus (great bulrush), both of which maintained dominance in the emergent wetland throughout the study (percent cover was 20 and 78 percent, respectively, in 2006). Echinocloa crusgalli (barnyard grass), an annual invasive from Eurasia, initially established in the newly disturbed soils of the scrub-shrub/forested wetland (areal cover was 56 percent in 2000), but by 2002, E. crusgalli was not growing in any sample plots and other species including Agrostis stolonifera (creeping bent grass), Festuca rubra (red fescue), Cornus spp. (dogwood species), and Salix nigra (black willow) were becoming more common. Sal. nigra contributed 30-percent cover in the scrub-shrub/forested wetland part by fall 2003. Rapid colonization of this species in subsequent years increased annual cover through 2006, when 15- to 25-foot tall Sal. nigra trees dominated the tree/shrub stratum (48 percent of the areal cover in 2006). The understory of the scrub-shrub/forested wetland was mostly shaded because of the canopy of Sal. nigra trees. Herbaceous species growing under and near the margins of the canopy included Ag. stolonifera and Ty. latifolia (29- and 23-percent areal cover, respectively). Flows in Valley Creek are responsible for transporting sediment and shaping the channel. Annual mean streamflow during the period the modified stream reach was monitored ranged from 15.2 cubic feet per second (ft3/s) in the 2002 water year to 53.0 ft3/s in the 2004 water year. This is a range of about 55 percent lower to 58 percent higher than the annual mean streamflow for the period of record. Despite the variability in streamflow, longitudinal profiles surveyed near rock deflectors in two short (100-foot) reaches within the modified stream reach maintained a constant slope throughout the monitoring period, most likely because of the presence of bedrock control. Cross-section geometry in the upstream reach was virtually unchanged during the monitoring period but 10 feet of bank migration was measured downstream, leaving the rock deflectors in mid-stream. As indicated by the change in channel morphology at the downstream reach, it is apparent that the rock deflectors were ineffective at adequately protecting the bank

Scientific Investigations Report↗

Flood-inundation maps for Indian Creek and Tomahawk Creek, Johnson County, Kansas, 2014

Digital flood-inundation maps for a 6.4-mile upper reach of Indian Creek from College Boulevard to the confluence with Tomahawk Creek, a 3.9-mile reach of Tomahawk Creek from 127th Street to the confluence with Indian Creek, and a 1.9-mile lower reach of Indian Creek from the confluence with Tomahawk Creek to just beyond the Kansas/Missouri border at State Line Road in Johnson County, Kansas, were created by the U.S. Geological Survey in cooperation with the city of Overland Park, Kansas. The flood-inundation maps, which can be accessed through the U.S. Geological Survey Flood Inundation Mapping Science Web site at http://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at the U.S. Geological Survey streamgages on Indian Creek at Overland Park, Kansas; Indian Creek at State Line Road, Leawood, Kansas; and Tomahawk Creek near Overland Park, Kansas. Near real time stages at these streamgages may be obtained on the Web from the U.S. Geological Survey National Water Information System at http://waterdata.usgs.gov/nwis or the National Weather Service Advanced Hydrologic Prediction Service at http://water.weather.gov/ahps/ , which also forecasts flood hydrographs at these sites. Flood profiles were computed for the stream reaches by means of a one-dimensional step-backwater model. The model was calibrated for each reach by using the most current stage-discharge relations at the streamgages. The hydraulic models were then used to determine 15 water-surface profiles for Indian Creek at Overland Park, Kansas; 17 water-surface profiles for Indian Creek at State Line Road, Leawood, Kansas; and 14 water-surface profiles for Tomahawk Creek near Overland Park, Kansas, for flood stages at 1-foot intervals referenced to the streamgage datum and ranging from bankfull to the next interval above the 0.2-percent annual exceedance probability flood level (500-year recurrence interval). The simulated water-surface profiles were then combined in a geographic information system with a digital elevation model derived from light detection and ranging data (having a 0.429-foot vertical and 0.228-foot horizontal accuracy) to delineate the area flooded at each water level. The availability of these maps, along with Web information regarding current stage from the U.S. Geological Survey streamgages and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel and residents with information that is critical for flood response activities such as evacuations, road closures, and postflood recovery efforts.

Kansas↗

Hydrologic appraisal of the Pine Barrens, Suffolk County, New York

The groundwater resources of the Pine Barrens, in southeastern Suffolk County, were appraised during 1982-83 for their potential for use as a supply of potable water. The study area encompasses approximately 50 square miles. The Precambrian bedrock is overlain by Cretaceous, Pleistocene, and Holocene deposits. The surficial material consists of morainal and outwash deposits overlain by recent beach and marsh deposits. The groundwater reservoir consists of three aquifers--the Lloyd and the Magothy of Cretaceous age and the upper glacial of Pleistocene age. The water is generally of suitable quality for drinking and most other uses except in several localized spots where point-source contamination has been documented. Precipitation is the sole source of recharge. Average annual precipitation recorded in the area during 1943-82 was 43 inches, about half of which reaches the groundwater reservoir. Overland runoff is estimated to be 0.5 inches per year, and evapotranspiration 23 inches per year. Total public water-supply withdrawal in the area in 1982 is estimated to have been 9.09 million gallons per day (Mgal/d). The upper glacial aquifer is the major source of public water; pumpage in 1982 was 7.88 Mgal/d; pumpage from the Magothy was 1.02 Mgal/d, and pumpage from the Lloyd aquifer was 0.19 Mgal/d. Stream discharge is measured periodically at 17 streamflow measurement sites on the 11 streams in the Pine Barrens. (USGS)

Water-Resources Investigations Report↗

Tributary-stream infiltration in Marsh Creek Valley, north-central Pennsylvania

The geohydrology of infiltration from five tributary streams along a 3.6-mile reach of Marsh Creek valley in north-central Pennsylvania was investigated during 1983-85. Marsh Creek valley is underlain by up to 100 feet of stratified drift that overlies Devonian bedrock. The stratified drift is overlain by up to 30 feet of alluvial-fan deposits near the tributary streams. Four of the five tributary streams lose large amounts of water to the stratified-drift aquifer in Marsh Creek valley. Along reaches away from the valley wall, infiltration losses from the streams averaged about 2 cubic feet per second per 1,000 feet of wetted channel length. Estimated hydraulic conductivity of the deposits near these streams ranges from 31 to 100 feet per day and averages 61 feet per day. Silty beds of lower permeability near the streams may significantly affect infiltration. The low permeability of the sediments near the fifth stream, which probably consist largely of fine-grained alluvium and swamp deposits, may account for the lack of infiltration losses along this stream. Tributary-stream infiltration accounted for more than 70 percent of the estimated recharge to the stratified-drift aquifer along the reach investigated during water year 1985, in which annual precipitation was below average. The sources of recharge and their estimated rates were: (1) direct infiltration of precipitation on the valley, 1.7 cubic feet per second; (2) unchanneled runoff and ground-water inflow from the uplands, 2.7 cubic feet per second; and (3) tributary-stream infiltration from Asaph Run, 3.7 cubic feet per second, Straight Run, 3.7 cubic feet per second, Dantz Run, 1.2 cubic feet per second, and Canada Run, 1.9 cubic feet per second. The temporal variation in recharge from tributary-stream infiltration greatly affects drawdowns caused by pumping from the wellfield at the National Fisheries Research and Development Laboratory near Straight Run.

Pennsylvania↗

Potential effects of climate change on streamflow, eastern and western slopes of the Sierra Nevada, California and Nevada

Precipitation-runoff models of the East Fork Carson and North Fork American Rivers were developed and calibrated for use in evaluating the sensitivity of streamflow in the north-central Sierra Nevada to climate change. The East Fork Carson River drains part of the rain-shadowed, eastern slope of the Sierra Nevada and is generally higher than the North Fork American River, which drains the wetter, western slope. First, a geographic information system was developed to describe the spatial variability of basin characteristics and to help estimate model parameters. The result was a partitioning of each basin into noncontiguous, but hydrologically uniform, land units. Hydrologic descriptions of these units were developed and the Precipitation- Runoff Modeling System (PRMS) was used to simulate water and energy balances for each unit in response to daily weather conditions. The models were calibrated and verified using historical streamflows over 22-year (Carson River) and 42-year (American River) periods. Simulated annual streamflow errors average plus 10 percent of the observed flow for the East Fork Carson River basin and plus 15 percent for the North Fork American River basin. Interannual variability is well simulated overall, but, at daily scales, wet periods are simulated more accurately than drier periods. The simulated water budgets for the two basins are significantly different in seasonality of streamflow, sublimation, evapotranspiration, and snowmelt. The simulations indicate that differences in snowpack and snowmelt timing can play pervasive roles in determining the sensitivity of water resources to climate change, in terms of both resource availability and amount. The calibrated models were driven by more than 25 hypothetical climate-change scenarios, each 100 years long. The scenarios were synthesized and spatially disaggregated by methods designed to preserve realistic daily, monthly, annual, and spatial statistics. Simulated streamflow timing was not very sensitive to changes in mean precipitation, but was sensitive to changes in mean temperatures. Changes in annual streamflow amounts were amplified reflections of imposed mean precipitation changes, with especially large responses to wetter climates. In contrast, streamflow amount was surprisingly insensitive to mean temperature changes as a result of temporal links between peak snowmelt and the beginning of warm-season evapotranspiration. Comparisons of simulations driven by temporally detailed climate-model changes in which mean temperature changes vary from month to month and simulations in which uniform climate changes were imposed throughout the year indicate that the snowpack accumulates the influences of short-term conditions so that season average climate changes were more important than shorter term changes.

Water-Resources Investigations Report↗

Regional analysis of annual precipitation maxima in Montana

Dimensionless precipitation-frequency curves for estimating precipitation depths having large recurrence intervals were developed for 2-, 6-, and 24-hour storm durations for three homogeneous regions in Montana. Within each homogeneous region, at-site annual precipitation maxima were made dimensionless by dividing by the at-site mean and grouped so that a single frequency curve would be applicable for each duration. L-moment statistics were used to help define the homogeneous regions and to develop the dimensionless precipitation- frequency curves. Data from 459 precipitation stations were used after application of statistical tests to ensure that the data were not serially correlated and were stationary over the general period of data collection (1900-92). The data were found to have a small, but significant, degree of interstation correlation. The GEV distribution was used to construct dimensionless frequency curves of annual precipitation maxima for each duration within each region. Each dimensionless frequency curve was considered to be reliable for recurrence intervals up to the effective record length. Because of significant, though small, interstation correlation in all regions for all durations, and because the selected regions exhibited some heterogeneity, the effective record length was considered to be less than the total number of station-years of data. The effective record length for each duration in each region was estimated using a graphical method and found to range from 500 years for 6-hour duration data in Region 2 to 5,100 years for 24-hour duration data in Region 3.

Water-Resources Investigations Report↗