Search USGS⌕ Search

SEARCH · Search USGS

Results for “Nature, Scientific Data”

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 937 records · Page 52Linked to original sources

Nutrient concentrations in surface water and groundwater, and nitrate source identification using stable isotope analysis, in the Barnegat Bay-Little Egg Harbor watershed, New Jersey, 2010–11

Five streams in the Barnegat Bay-Little Egg Harbor (BB-LEH) watershed in southern New Jersey were sampled for nutrient concentrations and stable isotope composition under base-flow and stormflow conditions, and during the growing and nongrowing seasons, to help quantify and identify sources of nutrient loading. Samples were analyzed for concentrations of total nitrogen, ammonia, nitrate plus nitrite, organic nitrogen, total phosphorus, and orthophosphate, and for nitrogen and oxygen stable isotope ratios. Concentrations of total nitrogen in the five streams appear to be related to land use, such that streams in subbasins characterized by extensive urban development (and historical agricultural land use)—North Branch Metedeconk and Toms Rivers—exhibited the highest total nitrogen concentrations (0.84–1.36 milligrams per liter (mg/L) in base flow). Base-flow total nitrogen concentrations in these two streams were dominated by nitrate; nitrate concentrations decreased during storm events as a result of dilution by storm runoff. The two streams in subbasins with the least development—Cedar Creek and Westecunk Creek—exhibited the lowest total nitrogen concentrations (0.16–0.26 mg/L in base flow), with organic nitrogen as the dominant species in both base flow and stormflow. A large proportion of these subbasins lies within forested parts of the Pinelands Area, indicating the likelihood of natural inputs of organic nitrogen to the streams that increase during periods of storm runoff. Base-flow total nitrogen concentrations in Mill Creek, in a moderately developed basin, were 0.43 to 0.62 mg/L and were dominated by ammonia, likely associated with leachate from a landfill located upstream. Total phosphorus and orthophosphate were not found at detectable concentrations in most of the surface-water samples, with the exception of samples collected from the North Branch Metedeconk River, where concentrations ranged from 0.02 to 0.09 mg/L for total phosphorus and 0.008 to 0.011 mg/L for orthophosphate. Measurements of nitrogen and oxygen stable isotope ratios of nitrate in surface-water samples revealed that a mixture of multiple subsurface sources, which may include some combination of animal and septic waste, soil nitrogen, and commercial fertilizers, likely contribute to the base-flow nitrogen load. The results also indicate that atmospheric deposition is not a predominant source of nitrogen transported to the BB-LEH estuary from the watershed, although the contribution of nitrate from the atmosphere increases during stormflow. Atmospheric deposition of nitrate has a greater influence in the less developed subbasins within the BB-LEH watershed, likely because few other major sources of nitrogen (animal and septic waste, fertilizers) are present in the less developed subbasins. Atmospheric sources appear to contribute proportionally less of the overall nitrate as development increases within the BB-LEH watershed. Groundwater samples collected from five wells located within the BB-LEH watershed and screened in the unconfined Kirkwood-Cohansey aquifer system were analyzed for nutrient and stable isotope composition. Concentrations of nitrate ranged from not detected to 3.63 mg/L, with the higher concentrations occurring in the highly developed northern portion of the watershed, indicating the likelihood of anthropogenic sources of nitrogen. Isotope data for the two wells with the highest nitrate concentrations are more consistent with fertilizer sources than with animal or septic waste. Total phosphorus was not detected in any of the wells sampled, and orthophosphate was either not detected or measured at very low concentrations (0.005–0.009 mg/L) in each of the wells sampled.

New Jersey↗

Methods for Estimating Water Withdrawals for Mining in the United States, 2005

The mining water-use category includes groundwater and surface water that is withdrawn and used for nonfuels and fuels mining. Nonfuels mining includes the extraction of ores, stone, sand, and gravel. Fuels mining includes the extraction of coal, petroleum, and natural gas. Water is used for mineral extraction, quarrying, milling, and other operations directly associated with mining activities. For petroleum and natural gas extraction, water often is injected for secondary oil or gas recovery. Estimates of water withdrawals for mining are needed for water planning and management. This report documents methods used to estimate withdrawals of fresh and saline groundwater and surface water for mining during 2005 for each county and county equivalent in the United States, Puerto Rico, and the U.S. Virgin Islands. Fresh and saline groundwater and surface-water withdrawals during 2005 for nonfuels- and coal-mining operations in each county or county equivalent in the United States, Puerto Rico, and the U.S. Virgin Islands were estimated. Fresh and saline groundwater withdrawals for oil and gas operations in counties of six states also were estimated. Water withdrawals for nonfuels and coal mining were estimated by using mine-production data and water-use coefficients. Production data for nonfuels mining included the mine location and weight (in metric tons) of crude ore, rock, or mineral produced at each mine in the United States, Puerto Rico, and the U.S. Virgin Islands during 2004. Production data for coal mining included the weight, in metric tons, of coal produced in each county or county equivalent during 2004. Water-use coefficients for mined commodities were compiled from various sources including published reports and written communications from U.S. Geological Survey National Water-use Information Program (NWUIP) personnel in several states. Water withdrawals for oil and gas extraction were estimated for six States including California, Colorado, Louisiana, New Mexico, Texas, and Wyoming, by using data from State agencies that regulate oil and gas extraction. Total water withdrawals for mining in a county were estimated by summing estimated water withdrawals for nonfuels mining, coal mining, and oil and gas extraction. The results of this study were distributed to NWUIP personnel in each State during 2007. NWUIP personnel were required to submit estimated withdrawals for numerous categories of use in their States to a national compilation team for inclusion in a national report describing water use in the United States during 2005. NWUIP personnel had the option of submitting the estimates determined by using the methods described in this report, a modified version of these estimates, or their own set of estimates or reported data. Estimated withdrawals resulting from the methods described in this report may not be included in the national report; therefore the estimates are not presented herein in order to avoid potential inconsistencies with the national report. Water-use coefficients for specific minerals also are not presented to avoid potential disclosure of confidential production data provided by mining operations to the U.S. Geological Survey.

Scientific Investigations Report↗

Hydrogeology and simulation of groundwater flow in the Lucerne Valley groundwater basin, California

The Lucerne Valley is in the southwestern part of the Mojave Desert and is about 75 miles northeast of Los Angeles, California. The Lucerne Valley groundwater basin encompasses about 230 square miles and is separated from the Upper Mojave Valley groundwater basin by splays of the Helendale Fault. Since its settlement, groundwater has been the primary source of water for agricultural, industrial, municipal, and domestic uses. Groundwater withdrawal from pumping has exceeded the amount of water recharged to the basin, causing groundwater declines of more than 100 feet between 1917 and 2016 in the center of the basin. The continued withdrawal has resulted in an increase in pumping costs, reduced well efficiency, and land subsidence near Lucerne Lake. Although the volume of pumping has declined in recent years, there is concern that new agricultural growth and limits on imported water will continue to strain the sustainability of the groundwater system. To address these concerns, the U.S. Geological Survey entered into a cooperative agreement with the Mojave Water Agency to develop a better understanding of the Lucerne Valley hydrogeologic system and provide tools to help evaluate and manage the effects of future development in the Lucerne Valley. The objectives of this study were to (1) improve the understanding of the aquifer system, (2) improve the understanding of subsidence in the basin, and (3) incorporate the understanding into a groundwater-flow model that can be used to help manage the groundwater resources in the Lucerne Valley. The model developed for this study covers the period of 1942–2016 and can help evaluate various proposed water-management scenarios during different climatic and hydrologic conditions. The aquifer system consists of a shallow aquifer, a confining unit, and middle and lower aquifers. These layered water-bearing units were identified based on geologic units of the mostly unconsolidated sediments and hydrologic properties. These alluvial deposits consist of clay, silt, sand, and gravel; some places also contain clay and silty clay lacustrine deposits. Several faults act, at least in part, as barriers to groundwater flow on the eastern, southern, and western edges of the basin. Present-day natural recharge is primarily from the infiltration of runoff from the San Bernardino Mountains to the south; however, stable and radioactive isotopes show that groundwater from the middle of the Lucerne Valley was older than about 10,000 years and probably was recharged as infiltration from streams draining the mountains in the Mojave Desert to the north, which probably does not occur under present-day climatic conditions. The annual average natural recharge for 1942–2016, estimated by a Basin Characterization Model, was about 635 acre-feet per year; the average amount of treated wastewater effluent transferred to the Lucerne Valley for artificial recharge annually ranged from about 1,500 to 4,000 acre-feet per year during 1980–2016. Pumpage estimates for 1942–2016 ranged from about 3,000 acre-feet in 1942 to about 18,300 acre-feet in 1984. The total cumulative amount of groundwater removed from the basin by pumping between 1942 and 2016 was estimated to be about 700,000 acre-feet, which was about 10 times greater than the cumulative amount of recharge to the entire Lucerne Valley groundwater basin. Before groundwater development, the direction of groundwater flow was from the southern part of the basin northward to discharge areas near Lucerne Lake, where it discharged through springs along the Helendale Fault and by evapotranspiration. Since the early 1900s, groundwater-level declines have mostly eliminated the areas where natural discharge occurred and exceeded 100 feet in the middle of the basin between the early 1950s and mid-1990s, and as much as 25 feet near the margins from about the mid-1950s to 2000s. A decrease in the rate of pumping after the mid-1990s lessened the hydraulic stress on the middle and lower aquifers and enabled hydraulic heads in the middle of the basin to recover slightly as groundwater near the margins of the basin moved toward the pumping depression. Although trends in groundwater levels in the center of the basin have reversed since the mid-1990s, levels at the basin margins continue to decline as the movement of groundwater from the margins fills the pumping depression and gradually flattens the groundwater table throughout the basin. The long-term extraction of groundwater and associated dewatering of the fine-grained sediments present within the aquifer system has resulted in aquifer compaction and consequently land subsidence, primarily near Lucerne Lake. Analysis of interferometric synthetic aperture radar data shows that almost 11 inches of land subsidence has occurred south of Lucerne Lake between April 1992 and November 2009; less subsidence occurred elsewhere in the basin during this period. This differential land subsidence has caused fissures and cracks in the ground surface, which have buckled the pavement and undercut roads in several locations. The Lucerne Valley Hydrologic Model was developed using the finite-difference groundwater modeling software One Water Hydrologic Model to represent the hydrologic conditions and stresses during 1942–2016. The model has a uniform grid of approximately 92 acres per cell (2,000 feet by 2,000 feet) and has four layers representing the water-bearing units. The results from the calibrated model simulations indicated that groundwater pumpage exceeded recharge, resulting in an estimated net cumulative depletion of groundwater storage (discharge minus recharge) of about 465,000 acre-feet from 1942 to 2016. The model simulated as much as 7.5 feet (90 inches; 2,286 millimeters) of aquifer compaction, which indicates the extensive fine-grained deposits and measured subsidence near Lucerne Lake.

California↗

Correlation analysis of groundwater and hydrologic data, Kaloko-Honokōhau National Historical Park, Hawai‘i

Designated in 1978, Kaloko-Honokōhau National Historical Park is located on the west coast of the Island of Hawaiʻi. The Kaloko-Honokōhau National Historical Park encompasses about 1,200 acres of coastal land and nearshore ecosystems, which include wetlands, anchialine pools (landlocked bodies of brackish water with hydrologic connections to the ocean), fishponds, a fishtrap, and coral reefs. These nearshore ecosystems are dependent on groundwater discharge with a freshwater component and provide habitat for threatened and endangered, endemic species, such as the orangeblack Hawaiian damselfly ( Megalagrion xanthomelas ) and the Hawaiian coot (ʻAlae keʻokeʻo, Fulica alai ). The populations of these native species, however, are threatened because of habitat loss related to urban development and environmental changes. Kaloko-Honokōhau National Historical Park is within the Keauhou aquifer system and the North Kona District, which experienced a 52 percent resident-population increase between 2000 and 2020 and a 41 percent visitor increase between 2008 and 2019. To support the current water demand associated with this growing population, groundwater is the primary source of freshwater used in the North Kona District, with about 15 million gallons of groundwater withdrawn from the Keauhou aquifer system per day since 2009. With anticipated development, future (2015–35) groundwater withdrawal from the Keauhou aquifer system is projected to be about 55 percent greater than recent (2012–14) withdrawal. Because Kaloko-Honokōhau National Historical Park is located within a coastal aquifer, natural and human-induced changes can affect the quality and quantity of groundwater, which can threaten groundwater-dependent ecosystems. To improve understanding of recent groundwater conditions, the U.S. Geological Survey, in cooperation with the National Park Service, undertook this study to document correlations between hydrologic time-series datasets from sites in and near Kaloko-Honokōhau National Historical Park using the nonparametric (distribution-free) Kendall’s tau statistical test. For the statistical analyses, dependent variables representing the groundwater system include groundwater level, the groundwater-level difference between pairs of sites, and specific conductance, and independent variables include datasets of sea level, rainfall, and groundwater withdrawal. About 34 percent of the 140 non-time-lagged Kendall’s tau statistical tests evaluated in this report are statistically significant ( p -value ≤ 0.050) with generally weak (0.1 ≤ tau ≤ 0.2) to moderate (0.2 ≤ tau ≤ 0.3) correlations. Groundwater levels measured at monitoring sites have the strongest correlation with the multivariate El Niño–Southern Oscillation index and withdrawal from production wells at the nearby Kohanaiki Private Club Community. Specific conductance is not consistently and significantly correlated with the independent hydrologic variables investigated in this report. Because the relations between hydrologic variables are commonly not instantaneous, a second set of correlations was evaluated after applying a range of time lags to the independent variable datasets. Relative to the non-time-lagged case (the set of correlations that did not use time-lagged independent variables), some of the time-lagged independent variables improved correlations with some of the dependent variables. For a particular independent variable, similar time lags were expected between the independent variable and dependent variable at all four monitoring sites. However, different time lags among the four sites sometimes produced the strongest correlations. This study identified several correlations that are statistically significant and hydrologically plausible, but the correlations could indicate that multiple concurrent factors are controlling the observed groundwater-system response, which might be better addressed using multivariate analyses. This study only investigates bivariate correlations, which may not explain all the variance in the data. The correlations analyzed in this report are limited by the quantity of available hydrologic data in the area near Kaloko-Honokōhau National Historical Park and are based on 14 years of time-series data, which were aggregated to a relatively coarse monthly temporal resolution that represents the minimum resolution common to all datasets.

Hawaii↗

Analyses of meteorological and hydrological records support Tribal members’ accounts of changing climate on the Fort Apache Reservation, east–central Arizona

The Fort Apache Reservation in east–central Arizona, home to the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, contains several climate zones because of the large variation in surface elevation within the reservation. This study was carried out in cooperation with the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, to raise awareness of how the changing climate affects the Fort Apache Reservation. This report documents the evaluation of existing multidecadal meteorological and hydrological datasets for the Fort Apache Reservation, used to evaluate the effects of a changing climate on the reservation. In this evaluation, near-surface air temperature, snow depth, snow water equivalent, precipitation, and streamflow datasets were analyzed for monotonic trends indicative of changing climatic conditions during specified periods of time. The results of these trend analyses were then compared with the Tribal community's memories of the changing climate. Trend analysis of near-surface air temperatures from a U.S. Historical Climatological Network station on the Fort Apache Reservation at Whiteriver, Arizona, indicated that mean annual air temperatures have increased by an average of 2.48 degrees Fahrenheit from 1980 to 2023. Records from the same station also indicated that average monthly maximum temperatures recorded for March increased by 5.39 degrees Fahrenheit for the same time period. Annual precipitation at the five precipitation stations used in this study decreased greatly from the 1980s to 2023. The largest total decrease was 10.07 inches, or 34.7 percent. However, only one of the two precipitation stations with longer term data available prior to 1980 had a significant negative trend when data from the entire period of record, from 1901 to 2023, were analyzed. Trend analyses show a decrease in the annual maximum snow water equivalent and an earlier disappearance of the snowpack at two Natural Resources Conservation Service snow telemetry stations in the mountainous region just east of the Fort Apache Reservation from 1981 to 2023. Based on the trend analyses, the average annual maximum snow water equivalent has decreased by more than 40 percent at both stations, and the average date when the snowpack was fully melted at the stations in the spring has moved earlier in time from late April to early April or late March. However, a statistically significant trend was not determined for the early April snow water equivalent measured at a nearby Natural Resources Conservation Service snow course across its period of record, indicating that the history of mountain snowpack in this area is not fully understood. Analysis of snowfall data from a National Oceanic and Atmospheric Administration Cooperative Observer Program network station on the Fort Apache Reservation at McNary 2N, AZ (station 025412) indicated that, on average, the measured total annual snowfall at the station decreased 42.4 percent from 1935 to 2023. Streamflow data from six U.S. Geological Survey streamgages on the Fort Apache Reservation were analyzed for trends. For most streamflow gages, statistically significant trends were not determined for tested parameters when the entire streamflow period of record was used for stations with records going back to at least the 1960s. However, when the data from 1980 to 2023 was tested, most of the streamflow parameters had statistically significant negative trends. All six streamgages showed a decrease in average annual runoff of at least 50 percent from 1980 to 2023; one streamgage showed an 81.8 percent decrease. A similar statistical finding was observed in the analysis of the annual spring snowmelt peak from one of the six streamgages used in the study and located in an area receiving measurable amounts of snowmelt runoff. When data from the entire period of record (1958–2023) was used, no trend in streamflow was determined; however, a significant negative trend was determined from 1980 to 2023, indicating a decrease in average annual springtime runoff of 62.6 percent. Statistical analysis on the timing of the annual spring snowmelt peak at the same streamgage indicated the snowmelt peak is happening on average about 12 days earlier now (2023) than it did in the past. The trend results for the timing of the annual spring snowmelt peak were the same and statistically significant for both periods tested (1958–2023 and 1980–2023). Two of the streamflow records from the Fort Apache Reservation were compared to the Palmer Hydrological Drought Index computed for Arizona Climate Division 4 (East Central) by the National Centers for Environmental Information. The comparison showed that the streamflow records generally tracked the Palmer Hydrological Drought Index. In interviews, Tribal community members living on the Fort Apache Reservation described the changes in climate that they observed during their lifetimes. Common themes reported were that air temperatures have become warmer, and the weather is less predictable with changes in seasonal patterns. Drier conditions, lower snowfall, shorter winters, and lower river levels were also reported. These community member observations align with the results of this study.

Arizona↗

Hydrochemical Regions of the Glacial Aquifer System, Northern United States, and Their Environmental and Water-Quality Characteristics

The glacial aquifer system in the United States is a large (953,000 square miles) regional aquifer system of heterogeneous composition. As described in this report, the glacial aquifer system includes all unconsolidated geologic material above bedrock that lies on or north of the line of maximum glacial advance within the United States. Examining ground-water quality on a regional scale indicates that variations in the concentrations of major and minor ions and some trace elements most likely are the result of natural variations in the geologic and physical environment. Study of the glacial aquifer system was designed around a regional framework based on the assumption that two primary characteristics of the aquifer system can affect water quality: intrinsic susceptibility (hydraulic properties) and vulnerability (geochemical properties). The hydrochemical regions described in this report were developed to identify and explain regional spatial variations in ground-water quality in the glacial aquifer system within the hypothetical framework context. Data analyzed for this study were collected from 1991 to 2003 at 1,716 wells open to the glacial aquifer system. Cluster analysis was used to group wells with similar ground-water concentrations of calcium, chloride, fluoride, magnesium, potassium, sodium, sulfate, and bicarbonate into five unique groups. Maximum Likelihood Classification was used to make the extrapolation from clustered groups of wells, defined by points, to areas of similar water quality (hydrochemical regions) defined in a geospatial model. Spatial data that represented average annual precipitation, average annual temperature, land use, land-surface slope, vertical soil permeability, average soil clay content, texture of surficial deposits, type of surficial deposit, and potential for ground-water recharge were used in the Maximum Likelihood Classification to classify the areas so the characteristics of the hydrochemical regions would resemble the characteristics of the clusters. The result of the Maximum Likelihood Classification is a map showing five hydrochemical regions of the glacial aquifer system. Statistical analysis of ion concentrations (calcium, chloride, fluoride, magnesium, sodium, potassium, sulfate, and bicarbonate) in samples collected from wells completed in the glacial aquifer system illustrates that variations in water quality can be explained, in part, by related environmental characteristics that control the movement of ground water through the aquifer system. A comparison of median concentrations of chemical constituents in ground water among the five hydrochemical regions indicates that ground water in the Midwestern Agricultural Region, the Urban-Influenced Region, and the Western Agriculture and Grassland Region has the highest concentrations of major and minor ions, whereas ground water in the Northern and Great Lakes Forested Region and the Mountain and Coastal Forested Region has the lowest concentrations of these ions. Median concentrations of barium, arsenic, lithium, boron, strontium, and nitrite plus nitrate as nitrogen also are significantly different among the hydrochemical regions.

Scientific Investigations Report↗

Hydrogeologic and Hydraulic Characterization of the Surficial Aquifer System, and Origin of High Salinity Groundwater, Palm Beach County, Florida

Previous studies of the hydrogeology of the surficial aquifer system in Palm Beach County, Florida, have focused mostly on the eastern one-half to one-third of the county in the more densely populated coastal areas. These studies have not placed the hydrogeology in a framework in which stratigraphic units in this complex aquifer system are defined and correlated between wells. Interest in the surficial aquifer system has increased because of population growth, westward expansion of urbanized areas, and increased utilization of surface-water resources in the central and western areas of the county. In 2004, the U.S. Geological Survey, in cooperation with the South Florida Water Management District, initiated an investigation to delineate the hydrogeologic framework of the surficial aquifer system in Palm Beach County, based on a lithostratigraphic framework, and to evaluate hydraulic properties and characteristics of units and permeable zones within this framework. A lithostratigraphic framework was delineated by correlating markers between all wells with data available based primarily on borehole natural gamma-ray geophysical log signatures and secondarily, lithologic characteristics. These correlation markers approximately correspond to important lithostratigraphic unit boundaries. Using the markers as guides to their boundaries, the surficial aquifer system was divided into three main permeable zones or subaquifers, which are designated, from shallowest to deepest, zones 1, 2, and 3. Zone 1 is above the Tamiami Formation in the Anastasia and Fort Thompson Formations. Zone 2 primarily is in the upper part or Pinecrest Sand Member of the Tamiami Formation, and zone 3 is in the Ochopee Limestone Member of the Tamiami Formation or its correlative equivalent. Differences in the lithologic character exist between these three zones, and these differences commonly include differences in the nature of the pore space. Zone 1 attains its greatest thickness (50 feet or more) and highest transmissivity in coastal areas. Zone 2, the most transmissive and extensive zone, is thickest (80 feet or more) and most transmissive in the inland eastern areas near Florida's Turnpike. In this area, zone 1 is absent, and the semiconfining unit above zone 2 extends to the land surface with a thickness commonly ranging from 50 to 100 feet. The thickness of zone 2 decreases to zero in most wells near the coast. Zone 3 attains its greatest thickness (100 feet or more) in the southwestern and south-central areas; zone 3 is equivalent to the gray limestone aquifer. The distribution of transmissivity was mapped by zone; however, zones 2 and 3 were commonly combined in aquifer tests. Maximum transmissivities for zone 1, zones 2 and 3, and zone 3 were 90,000, 180,000, and 70,000 ft2/d (feet-squared per day), respectively. The northern extent of the area with transmissivity greater than 50,000 ft2/d for zones 2 and 3 in the inland northeastern area along Florida's Turnpike has not been defined based on available data and could extend 5 to 10 miles farther north than mapped. Based on the thickness of zone 2 and a limited number of aquifer tests, a large area of zone 2 with transmissivity greater than 10,000 ft2/d, and possibly as much as 30,000 ft2/d, extends to the west across Water Conservation Area 1 from the inland southeastern area into the south-central area and some of the southwestern area. In contrast to the Biscayne aquifer present to the south of Palm Beach County, zones 2 and 3 are interpreted to be present principally in the Tamiami Formation and are commonly overlain by a thick semiconfining unit of moderate permeability. These zones have been referred to as the 'Turnpike' aquifer in the inland eastern areas of Palm Beach County, and the extent of greatest thickness and transmissivity follows, or is adjacent to, Florida's Turnpike. Where it is thick and transmissive, zone 1 may be considered equivalent to the Biscayne aquifer. Areas

Scientific Investigations Report↗

Lunar remote sensing and measurements

Remote sensing and measurements of the Moon from Apollo orbiting spacecraft and Earth form a basis for extrapolation of Apollo surface data to regions of the Moon where manned and unmanned spacecraft have not been and may be used to discover target regions for future lunar exploration which will produce the highest scientific yields. Orbital remote sensing and measurements discussed include (1) relative ages and inferred absolute ages, (2) gravity, (3) magnetism, (4) chemical composition, and (5) reflection of radar waves (bistatic). Earth-based remote sensing and measurements discussed include (1) reflection of sunlight, (2) reflection and scattering of radar waves, and (3) infrared eclipse temperatures. Photographs from the Apollo missions, Lunar Orbiters, and other sources provide a fundamental source of data on the geology and topography of the Moon and a basis for comparing, correlating, and testing the remote sensing and measurements. Relative ages obtained from crater statistics and then empirically correlated with absolute ages indicate that significant lunar volcanism continued to 2.5 b.y. (billion years) ago-some 600 m.y. (million years) after the youngest volcanic rocks sampled by Apollo-and that intensive bombardment of the Moon occurred in the interval of 3.84 to 3.9 b.y. ago. Estimated fluxes of crater-producing objects during the last 50 m.y. agree fairly well with fluxes measured by the Apollo passive seismic stations. Gravity measurements obtained by observing orbiting spacecraft reveal that mare basins have mass concentrations and that the volume of material ejected from the Orientale basin is near 2 to 5 million km 3 depending on whether there has or has not been isostatic compensation, little or none of which has occurred since 3.84 b.y. ago. Isostatic compensation may have occurred in some of the old large lunar basins, but more data are needed to prove it. Steady fields of remanent magnetism were detected by the Apollo 15 and 16 subsatellites, and the lunar dipole field was revised to no more than 6x 10 19 gauss. High-resolution mapping of fields of weak remanent magnetism (to 0.1 gamma) was made possible by the Apollo plasma and energetic-particle experiment. Although the causes of remanent magnetism are poorly understood, correlations with geologic units suggest the results may ultimately have farreaching significance to lunar history. Maria are much less structured by strong surface magnetic anomalies than the highlands. The strongest anomalies are associated with ejecta of farside basins, plains materials filling pre-Imbrian craters, and other old Imbrian to pre-Imbrian units. The high remanent fields could be due to cooling of ejecta units in an ancient magnetic field, lunar regolith maturity, extensive reworking and disruption of a magnetized layer, or simply surface roughness. Orbital geochemical experiments have shown that lunar high lands have larger Al: Si ratios and smaller Mg: Si ratios than maria. These two ratios are inversely related on a regional basis. With the exception of fresh craters, albedo and Al : Si ratios vary directly, showing that compositional differences as well as exposure of fresh materials are responsible for high albedos. Statistically treated data show that geologic contacts and compositional boundaries are concentric and can be roughly matched. Some craters on mare material have penetrated the mare fill, bringing highland-type materials to the surface. Natural radioactivity from thorium, potassium, and uranium is inversely correlated with elevation. Mare regions are enriched in iron, titanium, and magnesium relative to the highlands. Orbital bistatic-radar results provide estimates of surface roughness at two scale lengths (about 30 m and 250 m), which agree with visual estimates of roughness. The dielectric constant of the lunar surface, where sampled, is uniform to 13-cm radar and near 3. Slope frequency distributions measured by the radar vary and

Professional Paper↗

Characterization of and relations among precipitation, streamflow, suspended-sediment, and water-quality data at the U.S. Army Garrison Fort Carson and Piñon Canyon Maneuver Site, Colorado, water years 2016–18

Frequent and prolonged military training maneuvers are an intensive type of land use that may disturb land cover, compact soils, and have lasting effects on adjacent stream hydrology and ecosystems. To better understand the potential effect of military training on hydrologic and environmental processes, the U.S. Geological Survey in cooperation with the U.S. Army established hydrologic and water-quality data-collection networks at the U.S. Army Garrison Fort Carson (AGFC) in 1978 and at the Piñon Canyon Maneuver Site (PCMS) in 1982. The purpose of this report is to present precipitation, streamflow, suspended-sediment, and water-quality data collected by the U.S. Geological Survey at the AGFC and PCMS for water years (WYs) 2016–18 and to evaluate those data in relation to long-term data from the AGFC and PCMS. In WYs 2016–18, the U.S. Geological Survey monitored 26 sites on the AGFC and 17 sites on the PCMS for precipitation amount, streamflow, suspended sediment, and (or) water quality. On the AGFC, total annual precipitation in WYs 2016–18 was larger than the long-term mean for all 3 years at Rod and Gun Meteorologic Station at Fort Carson, CO (Rod and Gun). There were statistically significant upward trends in annual precipitation at Rod and Gun and Young Hollow Meteorologic Station at Fort Carson, CO (Young Hollow) with slopes of 1.25 and 0.66 inches per year (in/yr), respectively. The precipitation totals for WY 2017 were either the largest on record or in the top three for both sites and at Sullivan Park Meteorologic Station at Fort Carson, CO. On the PCMS, total annual precipitation was larger than the long-term mean in WYs 2016–18 at Brown Sheep Camp Meteorologic Station near Tyrone, CO; CIG Pipeline South Meteorologic Station near Simpson, CO; Bear Springs Hills Meteorologic Station near Houghton, CO (Bear Springs); and Upper Red Rock Canyon Meteorologic Station near Houghton, CO (Upper Red Rock). There were statistically significant upward trends in precipitation at Bear Springs and Upper Red Rock with slopes of 0.16 and 0.19 in/yr, respectively. The precipitation totals for WY 2017 were the largest on record for all sites except for Upper Bent Canyon Meteorological Station near Delhi, CO. Streamflow was calculated at 18 sites on the AGFC and 7 sites on the PCMS in at least 1 of WYs 2016–18. At AGFC, mean annual (or seasonal) streamflow in WYs 2016–18 was less than the long-term mean at 7 sites and greater than the long-term mean at 3 sites. There were statistically significant downward trends in mean annual or seasonal streamflow at Womack Ditch from Little Fountain Creek near Fort Carson, CO, and Ripley Ditch from Little Fountain Creek at Fort Carson, CO, with slopes of −0.036 and −0.028 cubic feet per second per year (ft 3 /s/y), respectively; and a significant upward trend in streamflow at Turkey Creek West Seepage below Teller Reservoir near Stone City, CO, with a slope of less than 0.001 ft 3 /s/y. Unlike for precipitation, the mean annual or seasonal streamflow for WY 2017 was not in the top 3 for any of the 12 sites with measured flow. At the PCMS, mean annual (or seasonal) streamflow was less than the long-term mean streamflow in WYs 2016–18 at the Taylor Arroyo below Rock Crossing near Thatcher, CO, and Bent Canyon Creek at Mouth near Timpas, CO, sites; and in WYs 2016 and 2018 at the Purgatoire River near Thatcher, CO (Purgatoire Thatcher), and Purgatoire River at Rock Crossing near Timpas, CO (Purgatoire Rock Crossing). There were no statistically significant trends in mean annual (or seasonal) streamflow at sites on the PCMS, and unlike for precipitation, the mean streamflow for WY 2017 was not in the top three for any sites except Purgatoire Rock Crossing. In WYs 2016–18, streamflow from sites on the AGFC and PCMS represented only a small fraction of streamflow in Fountain Creek or the Purgatoire River, and changes in streamflow that resulted from military maneuvers on the AGFC and PCMS were not likely to be detected in the downstream receiving waters. Suspended-sediment concentrations, loads, and yields for WYs 2016–18, were analyzed at two sites on the AGFC and five sites on the PCMS. On the AGFC, mean seasonal suspended-sediment concentrations ranged from 3.10 to 155 milligrams per liter (mg/L), mean seasonal suspended-sediment loads ranged from 0.04 to 27.1 tons per day (t/d), and seasonal suspended-sediment yields ranged from 0.28 to 216 tons per season per square mile (t/s/mi 2 ). Suspended-sediment yields at the two AGFC sites in WYs 2016–18 were all less than the long-term means. On the PCMS, mean seasonal suspended-sediment concentrations (at sites with some streamflow during a WY) ranged from 1.12 to 41.8 mg/L, mean suspended-sediment loads ranged from 0.01 to 13.1 t/d, and seasonal suspended-sediment yields ranged from 0.06 to 57.4 t/s/mi 2 . Suspended-sediment yields at the five PCMS sites in WYs 2016–18 were all less than the long-term means. In WYs 2016–18, mean daily suspended-sediment loads at Little Fountain were 1.3, 2.5, and 7.6 percent, respectively, of the mean daily suspended-sediment load at Fountain Creek at Security, Colorado. Likewise, the total of mean daily suspended-sediment loads from the five tributary sites to the Purgatoire River in WYs 2016–18 were about 0.25, 0.17, and 3.2 percent, respectively, of the historical mean daily suspended-sediment load at Purgatoire Thatcher. Spearman’s rank correlation coefficient was used to evaluate the strength and form of the relations between daily total precipitation and daily mean streamflow and between daily mean streamflow and suspended-sediment concentration and load for WYs 2016–18. For the sites on the AGFC and PCMS, there were weak or statistically insignificant positive correlations between precipitation and streamflow at nearby streamgauges, but strong statistically significant positive correlations between streamflow and suspended-sediment concentration and load. The ephemeral nature of the streams, quantity and timing of precipitation, air temperature, seasonal soil-moisture deficits, and effective runoff detention in erosion-control ponds could all contribute to inconsistent conversion of precipitation to streamflow. Water-quality data were analyzed for as many as 43 parameters from 9 samples collected from 3 sites on the AGFC and from 37 samples collected from 4 sites on PCMS during WYs 2016–18. The concentrations of selected water-quality parameters were compared to regulatory standards for aquatic life from the Colorado Department of Public Health and Environment (CDPHE) or aquatic-life criteria from the U.S. Environmental Protection Agency (EPA). There is at least 1 CDPHE standard or EPA criterion for 30 of the 43 water-quality parameters. For all samples from both the AGFC and the PCMS in WYs 2016–18, the concentrations of most water-quality parameters were compliant with the associated standards or criteria. However, there were some exceedances of standards or criteria: 11 samples exceeded the CDPHE recreational class standard for Escherichia coli concentration, 9 samples exceeded the CDPHE chronic unfiltered phosphorus aquatic-life standard, 36 samples exceeded the CDPHE chronic sulfate aquatic-life standard, 5 samples exceeded the EPA criterion for selenium, 7 samples exceeded the EPA criterion for aluminum, 2 samples exceeded the CDPHE chronic standard for iron, and 15 samples exceeded the CDPHE chronic standard for manganese. Identifying potential effects of military training on water quality in adjacent streams on the AGFC and PCMS is difficult due to the ephemeral nature of streamflow, limited number of sampling locations and samples, and limited access to the study areas. At the PCMS, pairs of water-quality samples were collected in March and May 2017 before and after an April–May 2017 military training event. At the Purgatoire Rock Crossing site, streamflow at the time of the May sample was approximately 35 times larger than streamflow for the March sample. The absolute percent differences of concentrations for 27 parameters ranged from −71.7 to 183 percent, and 7 parameters had increases in concentration whereas 22 parameters had no change or decreases in concentrations. The absolute percent differences of loads for 24 parameters ranged from 141 to 198 percent. The generally lower concentrations and higher loads were expected given the higher streamflows at the time of collection of the May compared to the March samples.

Colorado↗

Changes in streamflow and the flux of nutrients in the Mississippi-Atchafalaya River Basin, USA, 1980-2007

Nutrients and freshwater delivered by the Mississippi and Atchafalaya Rivers drive algal production in the northern Gulf of Mexico, which eventually results in the widespread occurrence of hypoxic bottom waters along the Louisiana and Texas coast. Researchers have demonstrated a relation between the extent of the hypoxic zone and the magnitude of streamflow, nutrient fluxes, and nutrient concentrations in the Mississippi River, with springtime streamflows and fluxes being the most predictive. In 1999 the U.S. Geological Survey (USGS) estimated the flux of nitrogen, phosphorus, and silica at selected sites in the Mississippi Basin and to the Gulf of Mexico for 1980-1996. These flux estimates provided the baseline information used by the Mississippi River/Gulf of Mexico Watershed Nutrient Task Force to develop an Action Plan for reducing hypoxia in the northern Gulf of Mexico. The primary goal of the Action Plan was to achieve a reduction in the size (areal extent) of the hypoxic zone from an average of approximately 14,000 square kilometers in 1996-2000 to a 5-year moving average of less than 5,000 square kilometers by 2015. Improved statistical models and adjusted maximum likelihood estimation using USGS Load Estimator (LOADEST) software were used to estimate annual and seasonal nutrient fluxes for 1980-2007 at selected sites on the Mississippi River and its tributaries. These data provide a means to evaluate the influence of natural and anthropogenic effects on delivery of water and nutrients to the Gulf of Mexico; to define subbasins that are the most important contributors of nutrients to the gulf; and to investigate the relations among streamflow, nutrient fluxes, and the size and duration of the Gulf of Mexico hypoxic zone. A comparative analysis between the baseline period of 1980-1996 and 5-year moving averages thereafter indicate that the average annual streamflow and fluxes of total nitrogen, nitrate, orthophosphate, and silica to the Gulf of Mexico have decreased. However, the flux of total phosphorus between the baseline period and subsequent 5-year periods has increased. The average spring (April, May, and June) streamflow and fluxes of silica, total nitrogen, nitrate, and orthophosphate to the Gulf of Mexico also decreased, whereas the spring flux of total phosphorus has increased. Similar changes in streamflow and nutrient flux were observed at many sites Buxtonwithin the basin. The inputs of water, total nitrogen, and total phosphorus from the major subbasins of the Mississippi-Atchafalaya River Basin as a percentage of the to-the-gulf totals have increased from the Ohio River Basin, decreased from the Missouri River Basin, and remained relatively unchanged from the Upper Mississippi, Red, and Arkansas River Basins. Changes in streamflow and nutrient fluxes are related, but short-term variations in sources of streamflow and nutrients complicate the interpretation of factors that affect nutrient delivery to the Gulf of Mexico. Parametric time-series models are used to try and separate natural variability in nutrient flux from changes due to other causes. Results indicate that the decrease in annual nutrient fluxes that has occurred between the 1980-1996 baseline period and more recent years can be largely attributed to natural causes (climate and streamflow) and not management actions or other human controlled activities in the Mississippi-Atchafalaya River Basin. The downward trends in total nitrogen, nitrate, ammonium, and orthophosphate that were detected at either the Mississippi River near St. Francisville, La., or the Atchafalaya River at Melville, La., occurred prior to 1995. In spite of the general decrease in nutrient flux, the average size of the Gulf of Mexico hypoxic zone has increased between 1997 and 2007. The reasons for this are not clear but could be due to the type or nature of nutrient delivery. Whereas the annual flux of total nitrogen to the Gulf of Mexico has decreased, the proporti

Scientific Investigations Report↗

Bathymetric and sediment facies maps for China Bend and Marcus Flats, Franklin D. Roosevelt Lake, Washington, 2008 and 2009

The U.S. Geological Survey (USGS) created bathymetric and sediment facies maps for portions of two reaches of Lake Roosevelt in support of an interdisciplinary study of white sturgeon (Acipenser transmontanus) and their habitat areas within Franklin D. Roosevelt Lake, Washington. In October 2008, scientists from the USGS used a boat-mounted multibeam echo sounder (MBES) to describe bathymetric data to characterize surface relief at China Bend and Marcus Flats, between Northport and Kettle Falls, Washington. In March 2009, an underwater video camera was used to view and record sediment facies that were then characterized by sediment type, grain size, and areas of sand deposition. Smelter slag has been identified as having the characteristics of sand-sized black particles; the two non-invasive surveys attempted to identify areas containing black-colored particulate matter that may be elements and minerals, organic material, or slag. The white sturgeon population in Lake Roosevelt is threatened by the failure of natural recruitment, resulting in a native population that consists primarily of aging fish and that is gradually declining as fish die and are not replaced by nonhatchery reared juvenile fish. These fish spawn and rear in the riverine and upper reservoir reaches where smelter slag is present in the sediment of the river lake bed. Effects of slag on the white sturgeon population in Lake Roosevelt are largely unknown. Two recent studies demonstrated that copper and other metals are mobilized from slag in aqueous environments with concentrations of copper and zinc in bed sediments reaching levels of 10,000 and 30,000 mg/kg due to the presence of smelter slag. Copper was found to be highly toxic to 30-day-old white sturgeon with 96-h LC50 concentrations ranging from 3 to 5 (u or mu)g copper per liter. Older juvenile and adult sturgeons commonly ingest substantial amounts of sediment while foraging. Future study efforts in Lake Roosevelt should include sampling of bottom material to confirm the presence or absence of slag material indicated on the map. In addition, follow-up acoustic work to determine stream velocities at varying discharges, in conjunction with sediment mapping, would be helpful to more accurately identify areas of scour and areas of sediment deposition where slag deposits may accumulate within the Marcus Flats and China Bend reaches. MBES mapping could also be used to determine changes in bed elevation and sedimentation in the study reaches and could help evaluate annual deposition and provide estimates on fine-sediment thickness.

Washington↗

Fluvial sediment dynamics in the Shoshone River and tributaries around Willwood Dam, Park County, Wyoming

Sedimentation affects many of the aging reservoirs in the United States. Dams and water diversions from rivers have been central elements of infrastructure supporting agricultural irrigation in the arid and semiarid regions of the Western United States for more than a century. The Willwood Irrigation District diversion dam (hereafter referred to as “Willwood Dam”) in Park County, Wyoming, is approximately 12 miles northeast of Cody, Wyo.; has a structural height of 70 feet; and impounds the Shoshone River for diversion into the Willwood Canal. Willwood Dam is part of a larger irrigation scheme supported by water storage in the much larger Buffalo Bill Dam, which is approximately 20 miles upstream. In October 2016, renovation construction activities at Willwood Dam and the Willwood Canal caused an unplanned evacuation of nearly 96,000 cubic yards of fine sediment. The fine sediment release in 2016 raised concerns that ongoing sediment management at Willwood Dam could impose limits on the long-term health of the aquatic ecosystem and fish populations. The U.S. Geological Survey, in cooperation with Wyoming Department of Environmental Quality and Willwood Work Groups 2 and 3, initiated an investigation of the dynamics of sediment transport in the Shoshone River and selected tributaries between Buffalo Bill Dam and Willwood Dam. The goal of the study was to quantify sediment transport into and out of Willwood Dam on an annual, seasonal, and event basis to better understand the relative quantities of sediment coming from natural sources and human activities on the landscape. The study ran from March 2019 through October 2021 and used observations of streamflow, turbidity, and acoustic backscatter collected at streamgages upstream and downstream from Willwood Dam to quantify suspended-sediment loads into and out of the dam during irrigation and fallow seasons, precipitation-runoff events, and deliberate sediment releases. Each tributary’s relative contribution to the sediment load upstream from Willwood Dam was examined using discrete measurements of suspended-sediment concentration and bedload during irrigation and fallow seasons, precipitation events, and stable conditions. Analysis of daily precipitation and temperature data indicated that conditions in the study area during the 2019 agricultural year were wetter and colder than period of record normal, and drier and near normal temperatures for the 2020 and 2021 agricultural years. Not all sediment load records between 2019 and 2021 are complete because of rejected observations (outliers), instrument failures or fouling, and instrument removal for calibrations. Statistical modeling of suspended-sediment concentration using paired values of turbidity and acoustic backscatter produced four models that, after refinement, had coefficients of determination indicating that more than 84 percent of the variance was explained by either turbidity or acoustic backscatter. A system of rules was developed to select the model predictions based on the seasonal operations of Willwood Dam, assumptions about the grain sizes mobilized during these operations, and assumed accuracy of the models at the downstream streamgage (Shoshone River below Willwood Dam, near Ralston, Wyo. [streamgage 06284010]) under different operational conditions. The sediment budget between upstream and downstream estimates of loads was interpreted using the mean predicted values bound by their respective model prediction intervals. When mean predicted loads of one streamgage were contained in the prediction intervals of the other streamgage, and vice-versa, difference in the sediment budget were interpreted as “indeterminate.” Modeled sediment load balances demonstrated the depositional and erosional behaviors expected from the conceptual model of dam operations whereby sediment tends to accumulate during irrigation seasons when the dam is spilling over the top, and sediment tends to evacuate during the fallow seasons when it is flowing through the sluice gates at the base of the dam. The sediment load calculations using the rules-based model criteria indicated that between 14,200 and 380,000 tons of suspended sediment moved through the Shoshone River around Willwood Dam during the irrigation seasons of 2019, 2020, and 2021; 380,000 tons of suspended sediment were transported during the cool, wet year of 2019, and 14,200 tons of suspended sediment were transported in 2020, which was relatively dry. During fallow seasons 2019, 2020, and 2021, which had fewer complete records, between 1,140 and 106,000 tons of suspended sediment was estimated to have moved through the river. For all seasons except fallow season 2022, the models estimated that more sediment was released from the dam than entered the dam, but the modeled mean loads at each streamgage were nearly always within the prediction intervals of each other, making the sediment balance indeterminant. Examination of suspended-sediment loads during irrigation seasons indicated that between 65 and 85 percent of fine sediment was transported during annual high flows and storm events, with the remainder transported during steady, lower streamflows. Examination of suspended loads during fallow seasons indicated that deliberate sediment releases through Willwood Dam accounted for between 39 and 67 percent of the total sediment moved during the fallow seasons. Deliberate sediment releases from Willwood Dam had estimated net exports of between 1,360 and 22,400 tons. Between August 2017 and July 2023, suspended-sediment concentration and bedload sediment samples were collected from 9 tributaries to the Shoshone River during 137 sampling events, including stable and precipitation-runoff conditions. During irrigation season precipitation events, the mean total sediment yields ranged from 0.33 to 9.51 tons per day per square mile; during fallow season precipitation events, the mean total yields ranged from 0.04 to 0.95 ton per day per square mile. The mean total sediment yield per unit area across all samples at each tributary site ranged from 0.26 to 3.08 tons per day per square mile. Bedload was a minor fraction of the total load, constituting a mean of 4 percent across all samples; 3 and 6 percent for events and nonevents, respectively, during irrigation season; and 3 and 1 percent for events and nonevents, respectively, during the fallow season. With the exception of one tributary, Dry Creek, these mean yield values were within the range of watershed-scale background sediment yield values estimated from reservoir surveys and previous suspended-sediment studies. Imagery from irrigation seasons 2012, 2015, 2017, 2019, and 2022 was used to determine the planimetric backwater extent of the pool area in the Shoshone River behind Willwood Dam to identify any changes in sediment storage. Active river channel widths in the Shoshone River upstream from Willwood Dam were all similar between years except 2015, which was determined to be statistically different from all other years. Bathymetric data taken in the pool behind Willwood Dam during three different surveys between November 2017 and April 2022 indicated no statistically significant differences in bed elevations between the years. Results from the planimetric and bathymetric survey data provide multiple lines of evidence indicating that sediment did not accumulate behind the dam within the error of the methods used. Examination of how precipitation affects sediment transport in the Shoshone River upstream from Willwood Dam indicated that accumulated rainfall from the natural runoff events captured during the study period varied from a trace to as much as 4.26 inches, with associated predicted suspended-sediment loads varying from 112 to 232,000 tons of suspended sediment. The behavior of the sediment loads relative to accumulated precipitation did not appear to change depending on irrigation or fallow season. A model of suspended-sediment concentrations relative to the 2-day accumulated precipitation indicated that suspended-sediment concentrations in the Shoshone River upstream from Willwood Dam increased exponentially for accumulations of 0.3 inch or more; such storms accounted for 10 percent or less of precipitation events observed during the 1981 to 2018 period of record. The gaps in records, precision of the instrumentation, and large variation in grain sizes in suspended-sediment mixtures downstream from the dam made closing the sediment budgets for most seasons unattainable. The biggest recent change in sediment storage measured using the planimetric area of deposits behind Willwood Dam took place between 2015 and 2017. The main event between these two measurements was the installation of new Willwood Canal gates in October 2016, which resulted in the large unplanned sediment release. Because the sediment budgets were nearly always indeterminate and the planimetric and bathymetric data indicated little change in the bed and bank material, it is likely that the change in sediment storage behind the dam during the study period was small relative to the precision of the statistical models and other uncertainties. This body of evidence suggests that, averaged during the 3-year study period, no major changes in storage took place, and that the current operations may be keeping storage at near-equilibrium. This condition could have been initiated because the middle sluice gate has now been operational since 2014, and the sediment release in October 2016 evacuated a large amount of legacy sediment from storage. Although the uncertainties are large, sluicing events allow for controlled releases of sediment that contributed to the near equilibrium conditions observed over an annual basis during this study.

Wyoming↗

Hydrogeology, hydrologic effects of development, and simulation of groundwater flow in the Borrego Valley, San Diego County, California

Executive Summary The Borrego Valley is a small valley (110 square miles) in the northeastern part of San Diego County, California. Although the valley is about 60 miles northeast of city of San Diego, it is separated from the Pacific Ocean coast by the mountains to the west and is mostly within the boundaries of Anza-Borrego Desert State Park. From the time the basin was first settled, groundwater has been the only source of water to the valley. Groundwater is used for agricultural, recreational, and municipal purposes. Over time, groundwater withdrawal through pumping has exceeded the amount of water that has been replenished, causing groundwater-level declines of more than 100 feet in some parts of the basin. Continued pumping has resulted in an increase in pumping lifts, reduced well efficiency, dry wells, changes in water quality, and loss of natural groundwater discharge. As a result, the U.S. Geological Survey began a cooperative study of the Borrego Valley with the Borrego Water District (BWD) in 2009. The purpose of the study was to develop a greater understanding of the hydrogeology of the Borrego Valley Groundwater Basin (BVGB) and to provide tools to help evaluate the potential hydrologic effects of future development. The objectives of the study were to (1) improve the understanding of groundwater conditions and land subsidence, (2) incorporate this improved understanding into a model that would assist in the management of the groundwater resources in the Borrego Valley, and (3) use this model to test several management scenarios. This model provides the capability for the BWD and regional stakeholders to quantify the relative benefits of various options for increasing groundwater storage. The study focuses on the period 1945–2010, with scenarios 50 years into the future. This report documents and presents (1) an analysis of the conceptual model, (2) a description of the hydrologic features, (3) a compilation and analysis of water-quality data, (4) the measurement and analysis of land subsidence by using geophysical and remote sensing techniques, (5) the development and calibration of a two-dimensional borehole-groundwater-flow model to estimate aquifer hydraulic conductivities, (6) the development and calibration of a three-dimensional (3-D) integrated hydrologic flow model, (7) a water-availability analysis with respect to current climate variability and land use, and (8) potential future management scenarios. The integrated hydrologic model, referred to here as the “Borrego Valley Hydrologic Model” (BVHM), is a tool that can provide results with the accuracy needed for making water-management decisions, although potential future refinements and enhancements could further improve the level of spatial and temporal resolution and model accuracy. Because the model incorporates time-varying inflows and outflows, this tool can be used to evaluate the effects of temporal changes in recharge and pumping and to compare the relative effects of different water-management scenarios on the aquifer system. Overall, the development of the hydrogeologic and hydrologic models, data networks, and hydrologic analysis provides a basis for assessing surface and groundwater availability and potential water-resource management guidelines.

California↗

Development of a flood-warning system and flood-inundation mapping in Licking County, Ohio

Digital flood-inundation maps for selected reaches of South Fork Licking River, Raccoon Creek, North Fork Licking River, and the Licking River in Licking County, Ohio, were created by the U.S. Geological Survey (USGS), in cooperation with the Ohio Department of Transportation; U.S. Department of Transportation, Federal Highway Administration; Muskingum Watershed Conservancy District; U.S. Department of Agriculture, Natural Resources Conservation Service; and the City of Newark and Village of Granville, Ohio. The inundation maps depict estimates of the areal extent of flooding corresponding to water levels (stages) at the following USGS streamgages: South Fork Licking River at Heath, Ohio (03145173); Raccoon Creek below Wilson Street at Newark, Ohio (03145534); North Fork Licking River at East Main Street at Newark, Ohio (03146402); and Licking River near Newark, Ohio (03146500). The maps were provided to the National Weather Service (NWS) for incorporation into a Web-based flood-warning system that can be used in conjunction with NWS flood-forecast data to show areas of predicted flood inundation associated with forecasted flood-peak stages. As part of the flood-warning streamflow network, the USGS re-installed one streamgage on North Fork Licking River, and added three new streamgages, one each on North Fork Licking River, South Fork Licking River, and Raccoon Creek. Additionally, the USGS upgraded a lake-level gage on Buckeye Lake. Data from the streamgages and lake-level gage can be used by emergency-management personnel, in conjunction with the flood-inundation maps, to help determine a course of action when flooding is imminent. Flood profiles for selected reaches were prepared by calibrating steady-state step-backwater models to selected, established streamgage rating curves. The step-backwater models then were used to determine water-surface-elevation profiles for up to 10 flood stages at a streamgage with corresponding streamflows ranging from approximately the 50 to 0.2-percent chance annual-exceedance probabilities for each of the 4 streamgages that correspond to the flood-inundation maps. The computed flood profiles were used in combination with digital elevation data to delineate flood-inundation areas. Maps of Licking County showing flood-inundation areas overlain on digital orthophotographs are presented for the selected floods. The USGS also developed an unsteady-flow model for a reach of South Fork Licking River for use by the NWS to enhance their ability to provide advanced flood warning in the region north of Buckeye Lake, Ohio. The unsteady-flow model was calibrated based on data from four flooding events that occurred from June 2008 to December 2011. Model calibration was approximate due to the fact that there were unmeasured inflows to the river that were not able to be considered during the calibration. Information on unmeasured inflow derived from NWS hydrologic models and additional flood-event data could enable the NWS to further refine the unsteady-flow model.

Ohio↗

Measured and Estimated Sodium-Adsorption Ratios for Tongue River and its Tributaries, Montana and Wyoming, 2004-06

The Tongue River drains an area of about 5,400 square miles and flows northward from its headwaters in the Bighorn National Forest of northeastern Wyoming to join the Yellowstone River at Miles City, Montana. Water from the Tongue River and its tributaries is extensively used for irrigation in both Wyoming and Montana. The Tongue River watershed contains vast coal deposits that are extracted at several surface mines. In some areas of the watershed, the coal beds also contain methane gas (coal-bed methane or natural gas), which has become the focus of intense exploration and development. Production of coal-bed methane requires the pumping of large volumes of ground water from the coal beds to reduce water pressure within the formation and release the stored gas. Water from the coal beds typically is high in sodium and low in calcium and magnesium, resulting in a high sodium-adsorption ratio (SAR). Disposal of ground water with high sodium concentrations into the Tongue River has the potential to increase salinity and SAR of water in the river, and potentially reduce the quality of water for irrigation purposes. This report documents SAR values measured in water samples collected at 12 monitoring sites in the Tongue River watershed and presents regression relations between specific conductance (SC) and SAR at each site for the years 2004-06. SAR in water samples was determined from laboratory-measured concentrations of sodium, calcium, and magnesium. The results of regression analysis indicated that SC and SAR were significantly related (p-values < 0.05) at most sites. The regression relations developed for most monitoring sites in the Tongue River watershed were used with continuous SC data to estimate daily SAR during the 2004 and 2005 irrigation seasons and to estimate 2006 provisional SAR values, which were displayed on the Web in real-time. Water samples were collected and analyzed from seven sites on the main stem of the Tongue River located at: (1) Monarch, Wyoming, station 06299980, (2) State line near Decker, Montana, station 06306300, (3) Tongue River Dam near Decker, Montana, station 06307500, (4) Birney Day School Bridge near Birney, Montana, station 06307616, (5) below Brandenberg Bridge near Ashland, Montana, station 06307830, (6) above T&Y Diversion Dam near Miles City, Montana, station 06307990, and (7) Miles City, Montana, station 06308500. Water samples were collected and analyzed from five sites on tributaries located at: (1) Goose Creek near Acme, Wyoming, station 06305700, (2) Prairie Dog Creek near Acme, Wyoming, station 06306250, (3) Hanging Woman Creek near Birney, Montana, station 06307600, (4) Otter Creek at Ashland, Montana, station 06307740, and (5) Pumpkin Creek near Miles City, Montana, station 06308400. All water-quality data for samples collected at these 12 sites can be accessed at Web sites http://waterdata.usgs.gov/mt/nwis or http://waterdata.usgs.gov/wy/nwis.

Scientific Investigations Report↗

Groundwater/surface-water interactions in the Bad River Watershed, Wisconsin

A groundwater-flow model was developed for the Bad River Watershed and surrounding area by using the U.S. Geological Survey (USGS) finite-difference code MODFLOW-NWT. The model simulates steady-state groundwater-flow and base flow in streams by using the streamflow routing (SFR) package. The objectives of this study were to: (1) develop an improved understanding of the groundwater-flow system in the Bad River Watershed at the regional scale, including the sources of water to the Bad River Band of Lake Superior Chippewa Reservation (Reservation) and groundwater/surface-water interactions; (2) provide a quantitative platform for evaluating future impacts to the watershed, which can be used as a starting point for more detailed investigations at the local scale; and (3) identify areas where more data are needed. This report describes the construction and calibration of the groundwater-flow model that was subsequently used for analyzing potential locations for the collection of additional field data, including new observations of water-table elevation for refining the conceptualization and corresponding numerical model of the hydrogeologic system. The study area can be conceptually divided into three primary hydrogeologic environments. The first encompasses the southern uplands with relatively low topographic relief, where groundwater-flow is unconfined and occurs primarily in sandy till and glacial outwash overlying Archean-aged crystalline bedrock. The second includes a transitional area of higher topographic relief and shallow depth to bedrock, in the vicinity of ridges formed by steeply dipping, early-Proterozoic aged metasedimentary units of the Marquette Range Supergroup (including the Ironwood Formation), and late-Proterozoic igneous units associated with the Midcontinent Rift System (MRS). Groundwater-flow in this area likely occurs primarily through connected networks of bedrock fractures that are not well characterized, and also in isolated pockets of Quaternary deposits. The third and last hydrogeologic environment includes lowlands along Lake Superior where a deep sandstone aquifer is confined by thick deposits of clay-rich till. Model input was compiled by using both published and unpublished data. Constant flux boundary conditions for the model perimeter were developed from a regional analytic element model described in appendix 1 of this report. Pumping from 26 high-capacity wells within the model area was included. The SFR stream network was developed from the National Hydrography Dataset (NHDPlus Version 2) and hydrography from the Wisconsin Department of Natural Resources (WDNR). Hydraulic conductivity values were determined for each model cell by interpolation from a network of pilot points, within zones representing major hydrogeologic units. Recharge to the groundwater system was estimated on a cell-by-cell basis by using the Soil Water Balance code (SWB), with gridded daily temperature and precipitation data for the period 1980&ndash;2011, and GIS coverages of soil and land-surface conditions. Estimated recharge varies considerably, following spatial patterns in the precipitation and soil hydrologic group inputs. The lowest recharge values occur in the Superior lowlands, whereas the highest values occur in the upland areas, especially those underlain by sandy soils, and in the vicinity of bedrock hills. The model was calibrated to groundwater-levels and base flows obtained from the USGS National Water Information System (NWIS) database, and groundwater-levels obtained from the WDNR and Band River Band well-construction databases. Calibration was performed via nonlinear regression by using the parameter-estimation software suite PEST. Groundwater levels and base-flow observations in the calibration dataset were well simulated by the calibrated model, with reasonable values of hydraulic conductivity. The pilot-point parameters that were most constrained by observations during model calibration coincided with the locations containing the most wells (head observations)&mdash;especially the population centers of Ashland, Mellen, and other communities along the major highway corridors. Results from the calibrated model illustrate differences in the nature of groundwater-flow within the watershed. In the southern part of the watershed, where bedrock is shallow, groundwater flow paths are relatively short, extending from local recharge areas to adjacent first and second-order streams. In contrast, laterally continuous deposits of clay-rich till covering the Superior Lowlands isolate most smaller streams from the sandstone aquifer, allowing for longer flow paths toward larger streams such as the Bad, Marengo, and White Rivers. Approximately three-quarters of all first-order stream cells were dry in the Superior Lowlands, compared to only half of first-order stream cells in the southern bedrock uplands. The model was used to delineate the groundwatershed for the Bad and Kakagon Rivers. &ldquo;Groundwatershed&rdquo; is defined as the area contributing groundwater discharge to one of these streams and their tributaries. The groundwatershed was found to align closely with the surface-watershed, with the most notable exception occurring along the southwestern half of Birch Hill, where surface water drains southwest towards the Potato River, and groundwater flows north and east towards Lake Superior. Similarly, the contributing area of groundwater-flow to the Reservation was delineated. Results indicate the off-Reservation groundwater contributing area to be limited in comparison to the extent of the watershed, extending southward into the highlands underlain by MRS igneous rock units, but not further into the area underlain by the Marquette Range Supergroup. Stable isotope samples were collected from 54 wells within the watershed, to investigate sources of groundwater. Oxygen-18 (&delta; 18 O) values lower than -13.0 per mil were documented in the sampling, and likely indicate the presence of recharge water from the last glacial period (>9,500 years old) beneath the northern portion of the Reservation, in the vicinity of Odanah, Wisconsin. Finally, a new data-worth analysis of potential new monitoring-well locations was performed by using the model. The relative worth of new measurements was evaluated based on their ability to increase confidence in model predictions of groundwater levels and base flows at 35 locations, under the condition of a proposed open-pit iron mine. Results of the new data-worth analysis, and other inputs and outputs from the Bad River model, are available through an online dynamic web mapping service at ( http://wim.usgs.gov/badriver/ ).

Wisconsin↗

Hypolimnetic dissolved-oxygen dynamics within selected White River reservoirs, northern Arkansas-southern Missouri, 1974-2008

Dissolved oxygen is a critical constituent in reservoirs and lakes because it is essential for metabolism by all aerobic aquatic organisms. In general, hypolimnetic temperature and dissolved-oxygen concentrations vary from summer to summer in reservoirs, more so than in natural lakes, largely in response to the magnitude of flow into and release out of the water body. Because eutrophication is often defined as the acceleration of biological productivity resulting from increased nutrient and organic loading, hypolimnetic oxygen consumption rates or deficits often provide a useful tool in analyzing temporal changes in water quality. This report updates a previous report that evaluated hypolimnetic dissolved-oxygen dynamics for a 21-year record (1974-94) in Beaver, Table Rock, Bull Shoals, and Norfork Lakes, as well as analyzed the record for Greers Ferry Lake. Beginning in 1974, vertical profiles of temperature and dissolved-oxygen concentrations generally were collected monthly from March through December at sites near the dam of each reservoir. The rate of change in the amount of dissolved oxygen present below a given depth at the beginning and end of the thermal stratification period is referred to as the areal hypolimnetic oxygen deficit. Areal hypolimnetic oxygen deficit was normalized for each reservoir based on seasonal flushing rate between April 15 and October 31 to adjust for wet year and dry year variability. Annual cycles in thermal stratification within Beaver, Table Rock, Bull Shoals, Norfork, and Greers Ferry Lakes exhibited typical monomictic (one extended turnover period per year) characteristics. Flow dynamics drive reservoir processes and need to be considered when analyzing areal hypolimnetic oxygen deficit rates. A nonparametric, locally weighted scatter plot smooth line describes the relation between areal hypolimnetic oxygen deficit and seasonal flushing rates, without assuming linearity or normality of the residuals. The results in this report are consistent with earlier findings that oxygen deficit rates and flushing-rate adjusted areal hypolimnetic oxygen deficit in Beaver and Table Rock Lakes were decreasing between 1974 and 1994. The additional data (1995-2008) demonstrate that the decline in flushing-rate adjusted areal hypolimnetic oxygen deficit in Beaver Lake has continued, whereas that in Table Rock Lake has flattened out in recent years. The additional data demonstrate the flushing-rate adjusted areal hypolimnetic oxygen deficit in Bull Shoals and Norfork Lakes have declined since 1995 (improved water quality), which was not indicated in earlier studies, while Greers Ferry Lake showed little net change over the period of record. Given the amount of data (35 years) for these reservoirs, developing an equation or model to predict areal hypolimnetic oxygen deficit and, therefore, areal hypolimnetic oxygen content, on any given day during future stratification seasons may be useful for reservoir managers.

Missouri;Arkansas↗

Groundwater recharge estimates for Maine using a Soil-Water-Balance model—25-year average, range, and uncertainty, 1991 to 2015

To address the lack of information on the spatial and temporal variability of recharge to groundwater systems in Maine, a study was initiated in cooperation with the Maine Geological Survey to use the U.S. Geological Survey Soil-Water-Balance model to evaluate annual average potential recharge across the State over a 25-year period from 1991 to 2015. The Maine Soil-Water-Balance model was calibrated using annual observations of recharge, runoff, and evapotranspiration for 32 calibration watersheds in the State during 2001–12 (902 total observations). Observations of recharge, runoff, and evapotranspiration were developed for each watershed to reduce the possibility of nonunique combinations of model parameters during the calibration. The Maine Soil-Water-Balance model was run using an optional evapotranspiration calculation method that provides more control for calibration than the standard method. The model was calibrated using the Parameter ESTimation software suite. The overall mean model error (average of all annual residuals for recharge, runoff, and precipitation) was 0.39 inch. The mean of the absolute value of the residuals, or the mean absolute error, was 2.32 inches. The root mean squared error for the calibrated model overall was 3.14 inches. Statistical tests indicated that the model residuals are normally distributed. To determine the potential uncertainty in the median annual potential recharge that results from uncertainty in the parameters as they relate to information contained in the observations, 300 alternate model realizations were run, and the standard deviation of the median potential recharge value at every pixel was calculated. Simulated 25-year median potential recharge across the State is widely variable; this variability closely follows patterns of precipitation, with additional variability contributed by the patchwork nature of the combinations of land-use class and hydrologic soil group inputs, and distribution of available water capacity in the soil across the State. Overall, the 25-year median annual potential recharge across the State is 7.5 inches, ranging from a low of about 5 inches to over 30 inches. The statewide range in the 25-year minimum values is from just over 2 inches to just over 20 inches. The statewide range in the 25-year maximum potential recharge is between 15 and 48 inches per year. The model areas with the highest simulated median potential recharge include areas underlain by type A soils (sandy and well drained), particularly those that also have land uses with low or little vegetation (blueberry barrens, developed, open space, scrub/shrub, and cropland, for example). The potential recharge values for these areas are similar to previously published values for comparable soil types. The 25-year average potential recharge grids were compared to recharge evaluated through groundwater-flow models or other methods in four hydrogeologic settings at six study areas in the State. A key factor in the ability of the Soil-Water-Balance model to reproduce the earlier study results was whether the available water-capacity data were an appropriate match for the hydrologic soil groups. The Maine Soil-Water-Balance model does a good job in representing an accurate potential recharge under circumstances where the surficial mapped soils extend below the surface to the water-table aquifer and where the available water-capacity data are in an appropriate range for the hydrologic soil group. One hydrogeologic setting that was challenging for the model was where a silt and clay layer was below a shallow soil unit that did not have available water-capacity data that were appropriate for the hydrologic soil group. In these cases, typically the available water-capacity data were very low, not accounting for the impedance of water flow provided by the underlying soil. The model also does not simulate well areas where bedrock surfaces are above the water table but below the plant rooting zone. The data products accompanying this report are intended to be used to provide first-cut estimates of recharge for geographic areas no smaller than the smallest watersheds used in the calibration of the model—or about 1.5 square miles. It is recommended that the grids are used to calculate an area-wide average potential recharge for any given area of study, and an uncertainty around the mean should be calculated from the standard deviation grid at the same time.

Maine↗