Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Natural Resource Modeling”

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

Estimating hydraulic properties of the Floridan Aquifer System by analysis of earth-tide, ocean-tide, and barometric effects, Collier and Hendry Counties, Florida

Aquifers are subjected to mechanical stresses from natural, non-anthropogenic, processes such as pressure loading or mechanical forcing of the aquifer by ocean tides, earth tides, and pressure fluctuations in the atmosphere. The resulting head fluctuations are evident even in deep confined aquifers. The present study was conducted for the purpose of reviewing the research that has been done on the use of these phenomena for estimating the values of aquifer properties, and determining which of the analytical techniques might be useful for estimating hydraulic properties in the dissolved-carbonate hydrologic environment of southern Florida. Fifteen techniques are discussed in this report, of which four were applied. An analytical solution for head oscillations in a well near enough to the ocean to be influenced by ocean tides was applied to data from monitor zones in a well near Naples, Florida. The solution assumes a completely non-leaky confining unit of infinite extent. Resulting values of transmissivity are in general agreement with the results of aquifer performance tests performed by the South Florida Water Management District. There seems to be an inconsistency between results of the amplitude ratio analysis and independent estimates of loading efficiency. A more general analytical solution that takes leakage through the confining layer into account yielded estimates that were lower than those obtained using the non-leaky method, and closer to the South Florida Water Management District estimates. A numerical model with a cross-sectional grid design was applied to explore additional aspects of the problem. A relation between specific storage and the head oscillation observed in a well provided estimates of specific storage that were considered reasonable. Porosity estimates based on the specific storage estimates were consistent with values obtained from measurements on core samples. Methods are described for determining aquifer diffusivity by comparing the time-varying drawdown in an open well with periodic pressure-head oscillations in the aquifer, but the applicability of such methods might be limited in studies of the Floridan aquifer system.

Florida↗

Research, monitoring, and evaluation of emerging issues and measures to recover the Snake River Fall Chinook salmon ESU

The portion of the Snake River fall Chinook salmon Oncorhynchus tshawytscha evolutionary significant unit (ESU) that spawns upstream of Lower Granite Dam transitioned from low to high abundance during 1992–2020 in response to U.S. Endangered Species Act recovery efforts and other federally mandated actions. This annual report focuses on changes in population abundance and habitat use by natural- and hatchery-origin spawners. Typically, we also report on population attributes of natural-origin juveniles, but data on juveniles were not collected in 2020 due to Covid-19. Spawners have located and used most of the available spawning habitat and that habitat is gradually approaching the point that no more redds can be supported. Timing of spawning and fry emergence have been relatively stable, but effects of density dependence are evident in juvenile life stages. Apparent abundance of juvenile fall Chinook salmon has increased and we noted the following changes: parr dispersal from riverine rearing habitat into Lower Granite Reservoir has become earlier; growth rate (g/d) and dispersal size of parr has declined; and passage timing of smolts from the two Snake River reaches has become earlier and downstream movement rate has increased. These findings coupled with stock-recruitment analyses presented in this report provide evidence for density-dependence in the Snake River reaches and in Lower Granite Reservoir resulting from the expansion of the recovery program. The long-term goal is to use this information in a comprehensive modeling effort to conduct action-effectiveness and uncertainty research and to inform Fish Population, Hydrosystem, Harvest, Hatchery, and Predation and Invasive Species Management Research, Monitoring, and Evaluation (RM&E) programs. In 2020, the U.S. Geological Survey (USGS) focused survey efforts in the Snake River on deepwater redd searches and fish collection for parentage-based tagging (PBT) analyses. We use a boat-mounted underwater video camera to count 170 deepwater redds at 19 of the 28 sites surveyed. Redd depths averaged 4.2 m. We collected genetic samples from 297 live fall Chinook salmon and 16 carcasses at 44 unique geographic locations that spanned 89 river kilometers. Seventy-two fish were recovered at Eureka Bar (rkm 307.1) and Corral Creek (rkm 349.7), which accounted for 23% of all collected fish in 2020. Most (238 fish) post-spawned salmon were collected from early to mid-November just after peak spawning. A summary of 2019 PBT results can be found in Appendix A.1. In 2020, we PIT tagged subyearling fall Chinook salmon in the Clearwater River to obtain population and growth data. In the Clearwater River, we tagged 2,192 subyearlings and recaptured 79 (3.6%) fish in the river and 187 fish (78 tagged by the U.S. Geological Survey, 109 tagged by the Nez Perce Tribe) at Lower Granite Dam during October which provided information for growth estimation. Within riverine habitats, growth in both length and mass were higher for fish tagged with 8-mm tags than with 9- and 12-mm tags. Estimated growth in length and mass of subyearlings was generally lower in Lower Granite Reservoir than in riverine habitats. Information on prey resources and juvenile fall Chinook salmon prey consumption was collected to better understand the growth opportunity of late-migrating fish in Lower Granite Reservoir. Zooplankton and surface drifting prey were collected from three reservoir locations from July through October during 2019 and 2020. Fall Chinook salmon diet data were collected from angled fish using gastric lavage. Cladocera and Copepoda were the most abundant zooplankton taxa collected while Diptera was the most common invertebrate taxon collected in surface drift samples. Totals of 49 and 94 juvenile fall Chinook salmon were captured in 2019 and 2020, respectively, and most fish were caught in the lower reach of the reservoir in October of each year. Juvenile fall Chinook salmon consumed mainly dipterans in July and October 2019 and mainly Daphnia during August–October in 2020. Fish showed selection mainly for dipterans in 2019 but strong selection for Daphnia during October 2020. Stomach fullness values were relatively low (<1.1%) during both years. Results show that prey resources are adequate in Lower Granite Reservoir to support positive fish growth during late summer and early autumn.

Idaho, Oregon, Washington↗

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↗

Water quality in the upper Shoal Creek basin, southwestern Missouri, 1999-2000

Results of a water-quality investigation of the upper Shoal Creek Basin in southwestern Missouri indicate that concentrations of total nitrite plus nitrate as nitrogen (NO2t+NO3t) in water samples from Shoal Creek were unusually large [mean of 2.90 mg/L (milligrams per liter), n (sample size)=60] compared to other Missouri streams (mean of 1.02 mg/L, n=1,340). A comparison of instantaneous base-flow loads of NO2t+NO3t indicates that at base-flow conditions, most NO2t+NO3t discharged by Shoal Creek is from nonpoint sources. Nearly all the base-flow instantaneous load of total phosphorus as P (Pt) discharged by Shoal Creek can be attributed to effluent from a municipal wastewater treatment plant. Samples collected from a single runoff event indicate that substantial quantities of Pt can be transported during runoff events compared to base-flow transport. Only minor quantities of NO2t+NO3t are transported during runoff events compared to base-flow transport. Fecal coliform bacteria densities at several locations exceed the Missouri Department of Natural Resources (MDNR) standard of 200 col/100 mL (colonies per 100 milliliters) for whole-body contact recreation. During 13 months of monitoring at 13 stream sites, fecal coliform densities (median of 277 and 400 col/100 mL) at two sites (sites 2 and 3) on Shoal Creek exceeded the MDNR standard at base-flow conditions. The maximum fecal coliform density of 120,000 col/100 mL was detected at site 3 (MDNR monitoring site) during a runoff event in April 1999 at a peak discharge of 1,150 ft3/s (cubic feet per second). Fecal coliform densities also exceeded the MDNR standard in three tributaries with the largest densities (median of 580 col/100 mL) detected in Pogue Creek. Results of ribopattern analyses indicate that most Escherichia coli (E. coli) bacteria in water samples from the study area probably are from nonhuman sources. The study area contains about 25,000 cattle, and has an estimated annual production of 33 million broilers and 300,000 turkeys. Probable nonhuman sources included turkeys, horses, chickens, and cattle; however, wildlife sources such as deer, raccoon, muskrat, and opossum were not evaluated. Human waste was an important source of E. coli in water samples collected at the MDNR monitoring site (site 3) on Shoal Creek and at two tributary sites (Joyce Creek and Woodward Creek). In general, the detection of human ribopatterns was consistent with the detection of organic compounds commonly associated with human wastewater such as caffeine, triclosan, or phenol, and the fecal indicators cholesterol and 3B-coprostanol. Ribopattern analysis indicate that horses were an important source of E. coli in Woodward Creek, which was consistent with horses being pastured immediately upstream from the sampling site on this creek. Pogue Creek contains a large density of turkey barns and five of eight E. coli isolates from one sample from Pogue Creek were matched to turkeys. Water samples from Pogue Creek generally did not contain detectable concentrations of human wastewater compounds, but one sample did contain detectable quantities of the antibiotics tylosin and lincomycin (widely used in the animal industry), and sulfamethoxazole (human use only). Although promising, the ability of ribopattern analyses to positively identify the source of a particular isolate is uncertain because of the small sample size, possible differences between animal source patterns in the study area and database used, lack of native wildlife source patterns, and variation in results depending on the number of possible animal host considered. Results of this study indicate that a trend of increasing fecal coliform densities with increasing time detected by the MDNR is, in part, caused by trends in annual precipitation and stream discharge, and not necessarily changes in land use or densities of animal operations. A multiple linear regression (MLR) model using specific conductance and wate

Water-Resources Investigations Report↗

Sea lamprey mark type, wounding rate, and parasite-host preference and abundance relationships for lake trout and other species in Lake Ontario

We examined how the frequency of attacks by Sea Lamprey on fishes in Lake Ontario varied in response to Sea Lamprey abundance and preferred host abundance (Lake Trout >432mm). For this analysis we assembled seven data sets. Two fishery independent surveys for Lake Trout: US Geological Survey (USGS)/New York State Department of Environmental Conservation (NYSDEC) south shore September gillnet assessment of adult Lake Trout (USGS/NYSDEC SGNS)( Lantry and Lantry 2011); and Ontario Ministry of Natural Resources (OMNR) monthly June-November community index gillnetting in northeastern Lake Ontario (OMNR CIS) (Ontario Ministry of Natural Resources 2011). One angler survey: NYSDEC April-September Fishing Boat Survey data collected along the south shore for Chinook and Coho salmon, and Rainbow and Brown trout (NYSDEC FBS) ( Lantry and Eckert 2012). Two spawning run datasets: OMNR north shore data including spring spawning runs of Rainbow Trout in the Ganaraska River and electroshocking data for fall spawning runs of Chinook and Coho salmon in the Credit River (Ontario Ministry of Natural Resources 2011); and NYSDEC data from the Salmon River on the southeastern shore including October spawning runs of Chinook and Coho salmon. One Sea Lamprey spawning survey: Department of Fisheries and Oceans Canada (DFO)/US Fish and Wildlife Service (USFWS) data for spawning-phase Sea Lamprey abundance collected from known spawning streams distributed throughout the Lake Ontario drainage basin (Mullet et al. 2003). One assessment of the abundance of dead Lake Trout: USGS/NYSDEC October-November bottom trawl collection of Lake Trout carcasses (Schneider et al. 1996). Annual patterns in A1, A2 and A3 wound stages did not track well in plots of wounding rates for the USGS/NYSDEC SGNS and correlations between A1 and later stages did not exist. A1 rates were not correlated to either Lake Trout abundance or Sea Lamprey numbers when considered alone, but were strongly correlated to the ratio between Sea Lamprey numbers and Lake Trout abundance (parasite/host ratio). While A2 and A3 rates were correlated to each other, neither was consistently correlated to any of the Lake Trout abundance or Sea Lamprey abundance parameters and sums of A1 to A3 rates did not improve correlations over those for A1 rates considered alone. Our analysis of the strain-specific susceptibility of Lake Trout to attack by Sea Lampreys extended the previous Schneider et al. (1996) analysis of three strains (SUP, CWL, and SEN) and 11 years of data 1982-1992 to an analysis of seven strains (SUP, CWL, SEN, JEN, LEW, ONT, and OXS) and two groups of unmarked fish (1983-1995 and 1996-2010) and included 18 more years of data through 2010. The susceptibility to attack for CWLs and SENs were below SUPs and nearly identical to the earlier values, new values for LEWs were greater than SUPS and values of unmarked Lake Trout prior to 1996 were unexpectedly greater than SUPs. By reexamining the Schneider et al. (1996) regression relationship between A1 wounding on Lake Trout and the incidence of Lake Trout carcasses recovered in fall bottom trawls (including three additional years of data), and substituting A1 wounding rate for total numbers of A1s observed which was used as the independent variable in the previous version, we were able to increase the variance explained by the relationship from an r2 of 0.60 to 0.88. Healing rate of wounds was explored by examining the monthly incidence of A1 and A2 wounds on Lake Trout from the OMNR CIS. Because wounding intensity varied between years and monthly sample size was frequently low, the ratio of A2 to A1s wounds was used to index how wounds accumulated or disappeared from the Lake Trout populations across seasons. The A2/A1 ratio decreased between June/July and October. A simple wounding model for Lake Trout was constructed to simulate the pattern of ratios by inflicting wounds on the population each month according to a predetermined distribution and including parameters for lethality (fixed distribution) and healing rate of A1 and A2 wounds (discretely varied between simulations). The best simulated representation (ratio size and monthly pattern) of the OMNR CIS data came from an A1 healing rate of 0.5 months and an A2 rate of 2 months. To examine whether alternate hosts provided reliable data to index damage caused by Sea Lampreys we compared September values of abundance and A1 wounding rates on Lake Trout, Sea Lamprey abundance, and the parasite/host ratio with NYSDEC Creel data for observations of Sea lampreys attached to sport-caught salmonids and to wounding observations for salmonids sampled in OMNR and NYSDEC spawning run assessments. Attachment frequency on NYSDEC Creel Chinook Salmon, Brown Trout, and Rainbow Trout were strongly correlated to all measures of Lake Trout abundance and wounding and to the parasite/host ratio. Chinook Salmon and Coho Salmon wounding observations for the NYSDEC Salmon River spawning run assessments and Rainbow Trout from the OMNR Ganaraska River spawning run assessment were strongly correlated to wounding measures for nearly all salmonids and to the parasite/host ratio.

Lake Ontario↗

Incorporating climate data into emergency planning and exercises: A primer for emergency management practioners and data developers

Climate change has and will continue to sharpen climate-related risks to communities and natural resources in California and elsewhere, as the probabilities of more extreme weather, floods, and fires continue to increase. This poses a problem of novel situations for emergency management. Progress has been made in terms of formally incorporating climate projections, data, and research on expected changes in climate-driven hazards into long-term hazard mitigation and climate adaptation strategies at both state and national levels. However, there are fewer examples of how climate change considerations have, or could be, incorporated into shorter-term emergency preparedness and response strategies. This is an important gap to fill, as climate resilience depends not only on mitigation and prevention measures, but also on the ability of agencies to coordinate and effectively minimize impacts when prevention measures fall short. The goal of this primer is to provide guidance on how to incorporate the best available information on climate variability and change into emergency management planning, with a focus on the development and use of extreme weather event scenarios for use in exercises. The first section is aimed toward a broad audience, including emergency management practitioners who use extreme weather event scenarios. It provides an overview of available data and tools that can inform scenario design as well as techniques for scenario design based on the hazard of interest, the audience and application, and the technical skills and resources required to develop, summarize, and/or visualize the data. This section concludes with an overview of approaches and lessons learned related to extreme event response planning and exercise design. Overall, this section highlights the advantages of developing quantitative scenarios based on spatial data, which allows visualizations and interactive data explorations that can provide greater specificity in discussions related to preparedness and response strategies. It further highlights the advantages of developing a core expert working group to guide planning, holding pre-exercise workshops to engage diverse communities outside of the emergency management sector, and engaging decisionmakers post-exercise to communicate key issues and outcomes as well as potential approaches for mitigating consequences that were identified by participants. The second section is aimed toward the scientific community and data developers involved in the creation of extreme weather event scenarios. This section provides technical guidance and detailed descriptions of four types of data resources and five analytical approaches that can be used to create extreme weather event scenarios based on the design considerations highlighted in section one. The computational resources and expertise required varies substantially across the options presented and is a primary consideration. These requirements, in addition to considerations related to audience and application, may determine the novelty and detail of the event, the detail of weather forecast information that can be provided, and the spatial extent across which the event can reasonably be modeled. The importance of, and approaches for, delivering information in a form that is accessible to emergency management practitioners is also discussed.

Report↗

Water-table contours, directions of ground-water movement, and measurements of inflow to American Falls Reservoir, southeastern Idaho, April 1984

In 1978 the U.S. Geological Survey began a 5-year study of the High Plains regional aquifer system to provide hydrologic information for evaluating the effects of long-term development of the aquifer and to develop a capability for predicting aquifer response to alternative changes in ground-water management. By use of a digital model, this report presents a quantitative description of the High Plains aquifer in Oklahoma. The High Plains aquifer consists predominantly of the Tertiary Ogallala Formation and overlying Quaternary alluvium and terrace deposits which are hydraulically connected to the High Plains aquifer. Much of the aquifer is underlain by formations of Permian through Cretaceous age, which generally have very small hydraulic conductivities. In some areas parts of underlying Triassic, Jurassic, or Cretaceous rocks are hydraulically connected with the aquifer. The High Plains aquifer is a water-table aquifer in which water moves generally to the east-southeast. Before the beginning of extensive irrigation of the 1960's, the aquifer was essentially in dynamic equilibrium with recharge from precipitation balanced by natural discharge from the aquifer. Ground-water discharge appeared in streams leaving the area or was returned to the atmosphere through evapotranspiration. Accurate records of irrigation pumpage are not available from the High Plains. In order to estimate irrigation pumpage, published records of crop distribution were used and a consumptive use was assigned to each principal irrigated crop. This method gave an estimated irrigation demand. Pumpage was taken as a percentage of the total irrigation demand. Irrigation has decreased ground-water discharge from the High Plains aquifer. Ground-water discharge was estimated as approximately 118 cubic feet per second in 1980. A finite-difference digital model was used to simulate flow in the High Plains aquifer. The recharge was adjusted so that 1980 ground-water discharge was 118 cubic feet per second, the estimated ground-water discharge for 1980. Recharge in the eastern half of the modeled area was 0.45 inch per year; one-half this value was used in the western half of the modeled area. Hydraulic conductivity was divided into three zones: 19.3 feet per day in the eastern zone; 16.2 feet per day in the central zone; and 8.28 feet per day in the western zone. A specific yield of 14.7 percent was used in the model. Using all these parameters, the model was calibrated so that the mean difference between predevelopment modeled and measured head was -0.044 foot. Following the calibration procedure, the model was used to predict the volumes of water in storage and distribution of saturated thickness in 1993 and 2020 using the 1980 pumping rates. The calculated quantity of water in storage in the aquifer in 1941 (predevelopment) was approximately 135.2 million acre-feet; in 1980, approximately 121.9 million acre-feet; in 1993 approximately 112.7 million acre-feet; and in 2020, approximately 96.2 million acre-feet. The High Plains aquifer in Oklahoma will continue to be an important source of water past the year 2000. As withdrawals continue from the aquifer at the present rate, the water table will continue to decline and when the water table drops below the streambed in any part of the area, ground-water discharge to streams will cease in that area. Based on the calculated volumes of water in storage, the volume of water remaining in storage as compared to the predevelopment volume is as follows: 90 percent in 1980, 83 percent in 1993, and 71 percent in 2020.

Idaho↗

Preliminary assessment of injection, storage, and recovery of freshwater in the lower Hawthorn aquifer, Cape Coral, Florida

A preliminary assessment of subsurface injection, storage and recovery of fresh canal water was made in the naturally brackish lower Hawthorn aquifer in Cape Coral, southwestern Florida. A digital modeling approach was used for this preliminary assessment, incorporating available data on hydrologic conditions, aquifer properties, and water quality to simulate density-dependent ground-water flow and advective-dispersive transport of a conservative ground-water solute (chloride ion). A baseline simulation was used as reference to compare the effects of changing various operational factors on the recovery efficiency. A recovery efficiency of 64 percent was estimated for the baseline simulation. Based on the model, the recovery efficiency increases if the injection rate and recovery rates are increased and if the ratio of recovery rate to injection rate is increased. Recovery efficiency decreases if the amount of water injected is increased; slightly decreases if the storage time is increased; is not changed significantly if the water is injected to a specific flow zone; increases with successive cycles of injection, storage, and recovery; and decreases if the chloride concentrations in either the injection water or native aquifer water are increased. In everal hypothetical tests, the recovery efficiency fluctuated between 22 and about 100 percent. Two successive cycles could bring the recovery efficiency from 60 to about 80 percent. Interlayer solute mass movement across the upper and lower boundaries seems to be the most important factor affecting the recovery efficiency. A sensitivity analysis was performed applying a technique in which the change in the various factors and the corresponding model responses are normalized so that meaningful comparisons among the responses could be made. The general results from the sensitivity analysis indicated that the permeabilities of the upper and lower flow zones were the most important factors that produced the greatest changes in the relative sensitivity of the recovery efficiency. Almost equally significant changes occurred in the relative sensitivity of the recovery efficiency when all porosity values of the upper and lower flow zones and the leaky confining units and the vertical anisotropy ratio were changed. The advective factors are the most important in the Cape Coral area according to the sensitivity analysis. However, the dispersivity values used in the model were extrapolated from studies conducted at the nearby Lee County Water Treatment Plant, and these values might not be representative of the actual dispersive characteristics of the lower Hawthorn aquifer in the Cape Coral area.

Florida↗

Simulation of ground-water flow and potential land subsidence, upper Santa Cruz Basin, Arizona

A numerical ground-water flow model of the upper Santa Cruz basin in Pinal, Pima, and Santa Cruz Counties was developed to evaluate predevelopment conditions in 1940, ground-water withdrawals for 1940-86, and potential water-level declines and land subsidence for 1987-2024. Simulations of steady-state ground-water conditions indicate 12,900 acre-feet of ground-water inflow, 15,260 acre-feet of outflow, 53,000 acre-feet of pre- development pumpage, 29,840 acre-feet of mountain- front recharge, and 34,020 acre-feet of streamflow infiltration in 1940. Simulations of transient ground-water conditions indicate a total of 6.6 million acre-feet of net pumpage and 3.4 million acre-feet of water removed from aquifer storage for 1941-86. A difference of 1.2 million acre-feet between estimated and net pumpage is attributed to increased recharge from irrigation return flow, mine return flow, and infiltration of sewage effluent. Estimated natural recharge represents 40 percent of pumpage for 1966-86 and averaged 63,860 acre-feet per year for 1940-57 and 76,250 acre-feet per year for 1958-86. The increase in recharge after 1958 was coincident with above- average winter streamflow in the Santa Cruz River for 1959-86. Increased recharge after 1958 and decreased pumpage after 1975 contributed to decreased water-level declines or to recoveries after 1977 in wells near the Santa Cruz River and its tributaries. The results of projection simu- lations indicate that a maximum potential subsi- dence for 1987-2024 ranges from 1.2 feet for an inelastic specific storage of .0001 ft to 12 feet for an inelastic specific storage of .0015 ft. The simulations were made on the basis of pumpage and recharge rates from 1986 and by using a preconso- lidation-stress threshold of 100 feet. A permanent reduction in acquitard storage can range from 1 to 12 percent of the potential loss of 3.9 million acre-feet in aquifer-system storage for 1987-2024.

Water-Resources Investigations Report↗

Historical changes in bed elevation and water depth within the Nehalem Bay, Oregon, 1891–2019

Estuaries, at the nexus of rivers and the ocean, are depositional areas that respond to changes in streamflow, tides, sea level, and inputs of sediment from marine and watershed sources. Understanding changes in bed elevations, deposited and eroded sediment, and water depth throughout estuaries is relevant for understanding their present-day status and long-term evolution, identifying potential hazards to human communities, and informing estuarine conservation. In response to observations of sedimentation in the Nehalem Bay, northwestern Oregon, by the Port of Nehalem, the magnitudes and patterns of bathymetric change in the Bay were documented and described by two approaches. The first approach compared changes in bed elevation with estimated volumes of erosion and deposition from overlapping survey data acquired in 1957 and 2019 for the area of the Nehalem Bay from upstream of the Highway 101 bridge to downstream of Fishery Point. The second approach examined changes in water depth for seven zones from the confluence of the North Fork and Nehalem Rivers to the mouth of the Nehalem River using nautical charts (1891, 1947, 1970, 1990, and 2004). These two approaches were used because the bathymetric surveys from 1957 and 2019 could be tied to a common vertical datum, allowing for a direct comparison of changes in bed elevations, whereas the nautical charts could not be tied to a common vertical datum, which limited the analyses to a comparison of changes in water depths over a broader time frame. Bed elevation changes from 1957 to 2019 were assessed from upstream of the Highway 101 bridge to downstream of Fishery Point where the two surveys overlapped (2 square kilometers) using thalweg longitudinal profiles, channel cross sections, and digital elevation models (DEMs) showing the elevation differences between the two surveys (or DEMs of difference). The most prominent change between 1957 and 2019 was the migration of the thalweg (or deepest part of the channel) between the downstream end of Lazarus Island and downstream of Fishery Point; this migration resulted in sediment deposition in the former thalweg and sediment erosion in formerly shallow areas to form the new thalweg. Bed elevation changes in the thalweg also varied longitudinally between 1957 and 2019. The bed elevation of the thalweg in both surveys, however, was generally less than 1 meter (m). The thalweg in the area of overlapping surveys shortened from about 7.0 to 6.7 kilometers in length over that same period. The bed elevation changes between the DEMs showed that maximum erosion and deposition was 4.3 and 4.5 m, respectively. In this same time period, the net change in sediment volume was 230,000 cubic meters (m 3 ), indicating net deposition. However, the error estimated for the 95 percent confidence interval analyses is ±315,000 m 3 , and therefore does not preclude the possibility that net erosion may have occurred. Historical changes in water depth from soundings depicted on nautical charts from 1891, 1947, 1970, 1990, and 2004 were evaluated by assessing spatial and temporal changes for seven zones of the Nehalem Bay. Across all years and zones, water depths ranged from about 0.2 to 9.4 m, whereas median water depths ranged from 0.3 to 6.4 m. Median depths and the range of water depths did not systematically increase or decrease throughout all zones during the same periods. In all nautical charts, the zone at the mouth of the Nehalem River consistently had the deepest soundings (7.9 to 9.4 m) and the greatest range of water depths (7.3 to 8.8 m). Qualitative evaluation of the nautical charts showed minimal changes in the overall shape of the Nehalem Bay. The exception to this observation was at the mouth of the Bay, where two historical outlets to the Pacific Ocean depicted in the 1891 nautical chart were reduced to one outlet following the construction of jetties (1916 and 1918). The results of this study emphasize that bed elevations and water depths within the Nehalem Bay have varied between 1891 and 2019, as illustrated by the lateral and vertical changes in the thalweg and changes in water depths over time. Changes in thalweg position and related patterns of sediment erosion and deposition are expected in the future as the Nehalem Bay continues to respond to changes in tides, sea level, streamflow, and sediment inputs from watershed and marine sources. The results of this study and the surveys from 1957 and 2019 provide a foundation for documenting and evaluating future changes in the Nehalem Bay and prioritizing actions to manage and protect natural resources and recreational access to the Nehalem Bay.

Oregon↗

Projected effects of proposed chloride-control projects on shallow ground water– Preliminary results for the Wichita River basin, Texas

The U.S. Army Corps of Engineers' plan to control the natural chloride pollution in the Wichita River basin includes the construction of Truscott Brine Lake on a tributary of the North Wichita River. In connection with the proposed brine lake, the U.S. Geological Survey was requested to: (1) Define the existing ground-water conditions in the shallow fresh-water system of the project area; and (2) project the post-construction effects of the proposed lake on the fresh-water aquifer, especially in relation to hydraulic-head changes but also with respect to possible changes in the chemical quality of the ground water. The fresh-water aquifer in the project area is a shallow water-table system with relatively fresh water (generally a calcium sulfate type) that contains approximately 500 to 5,000 milligrams per liter of dissolved solids. The aquifer consists of Permian rocks with very small values of hydraulic conductivity and overlies a brine system (sodium chloride type) that is even less permeable. A thin transition zone separates the fresh-water aquifer and the brine system. Small quantities of infiltration from precipitation throughout the area's watershed constitute the recharge to the aquifer. Discharge from the aquifer consists of the small base flow along creeks; well discharge is negligible. Two-dimensional mathematical computer models were developed for aquifer simulation of: (1) Steady-state conditions in a fresh-water system and (2) transient conditions in a brine- fresh-water system where the density effects of the brine are considered. The main results 'of projecting the effects of the proposed Truscott Brine Lake on the fresh-water aquifer are: (1) Hydraulic head rises of 5 to 40 feet would be confined to areas near the proposed dam and along the lake shoreline, and (2) migration of salt water downstream from the dam generally would be limited to less than 1 mile and apparently would not reach equilibrium during the 100-year duration of the project. The modeling efforts did not include possible effects related to hydrodynamic dispersion in the brine- fresh-water system. Possible changes in the hydraulic conductivity of the aquifer, due to physical and chemical interactions in the brine and fresh-water environments, also were not considered.

Texas↗

Ground-water quality in the Appalachian Plateaus, Kanawha River basin, West Virginia

Water samples collected from 30 privately-owned and small public-supply wells in the Appalachian Plateaus of the Kanawha River Basin were analyzed for a wide range of constituents, including bacteria, major ions, nutrients, trace elements, radon, pesticides, and volatile organic compounds. Concentrations of most constituents from samples analyzed did not exceed U.S. Environmental Protection Agency (USEPA) standards. Constituents that exceeded drinking-water standards in at least one sample were total coliform bacteria, Escherichia coli (E. coli) , iron, manganese, and sulfate. Total coliform bacteria were present in samples from five sites, and E. coli were present at only one site. USEPA secondary maximum contaminant levels (SMCLs) were exceeded for three constituents -- sulfate exceeded the SMCL of 250 mg/L (milligrams per liter) in samples from 2 of 30 wells; iron exceeded the SMCL of 300 µg/L (micrograms per liter) in samples from 12 of the wells, and manganese exceeded the SMCL of 50 µg/L in samples from 17 of the wells sampled. None of the samples contained concentrations of nutrients that exceeded the USEPA maximum contaminant levels (MCLs) for these constituents. The maximum concentration of nitrate detected was only 4.1 mg/L, which is below the MCL of 10 mg/L. Concentrations of nitrate in precipitation and shallow ground water are similar, potentially indicating that precipitation may be a source of nitrate in shallow ground water in the study area. Radon concentrations exceeded the recently proposed maximum contaminant level of 300 pCi/L at 50 percent of the sites sampled. The median concentration of radon was only 290 pCi/L. Radon-222 is a naturally occurring, carcinogenic, radioactive decay product of uranium. Concentrations, however, did not exceed the alternate maximum contaminant level (AMCL) for radon of 4,000 pCi/L in any of the 30 samples. Arsenic concentrations exceeded the proposed MCL of 5µg/L at 4 of the 30 sites. No samples exceeded the current MCL of 50 µg/L. Neither pesticides nor volatile organic compounds (VOCs) were prevalent in the study area, and the concentrations of the compounds that were detected did not exceed any USEPA MCLs. Pesticides were detected in only two of the 30 wells sampled, but four pesticides -- atrazine, carbofuran, DCPA, and deethylatrazine -- were detected in one well; molinate was detected in the other well. All of the pesticides detected were at estimated concentrations of only 0.002 µg/L. Of the VOCs detected, trihalomethane compounds (THMs), which can result from chlorination of a well, were the most common. THMs were detected in 13 of the 30 wells sampled. Gasoline by-products, such as benzene, toluene, ethylbenzene and xylene (BTEX compounds) were detected in 10 of the 30 wells sampled. The maximum concentration of any of the VOCs detected in this study, however, was only 1.040 µg/L, for the THM dichlorofluoromethane. Water samples from 25 of the wells were analyzed for chlorofluorocarbons (CFCs) to estimate the apparent age of ground water. The analyses indicated that age of water ranged from 10 to greater than 57 years, and that the age of ground water could be correlated with the topographic setting of the wells sampled. Thus the apparent age of water in wells on hilltops was youngest (median of 13 years) and that of water in wells in valleys was oldest (median of 42 years). Water from wells on hillsides was intermediate in age (median of 29 years). These data can be used to define contributing areas to wells, corroborate or revise conceptual ground-water flow models, estimate contaminant travel times from spills to other sources such as nearby domestic or public supply wells, and to manage point and nonpoint source activities that may affect critical aquifers.

West Virginia↗

Status of sea otter populations in southcentral and southeast Alaska, 2002-2003

During the years 2002-2004 estimated sea otter population sizes were calculated for Southeast Alaska, Prince William Sound, and the Kenai Peninsula and Cook Inlet regions of Alaska. Aerial surveys were conducted by a single observer from a float-equipped Bellanca Scout fixed-wing aircraft flying at 91 m altitude and 65 mph. The surveys followed protocols written by Bodkin and Udevitz (1999). The survey design consisted of systematic sampling of 400 m wide transects that were uniformly placed throughout the survey area. Selection and sampling of transects was proportional to expected sea otter abundance, with most effort taking place in transects over waters 0-40 m in depth. Intensive searches were periodically conducted within transects to estimate the proportion of sea otters not detected on strips. To obtain an adjusted population size estimate, strip counts are adjusted for the area not surveyed and by a correction factor.

Alaska↗

Hydrology and simulation of ground-water flow in Cedar Valley, Iron County, Utah

Cedar Valley, located in the eastern part of Iron County in southwestern Utah, is experiencing rapid population growth. Cedar Valley traditionally has supported agriculture, but the growing population needs a larger share of the available water resources. Water withdrawn from the unconsolidated basin fill is the primary source for public supply and is a major source of water for irrigation. Water managers are concerned about increasing demands on the water supply and need hydrologic information to manage this limited water resource and minimize flow of water unsuitable for domestic use toward present and future public-supply sources. Surface water in the study area is derived primarily from snowmelt at higher altitudes east of the study area or from occasional large thunderstorms during the summer. Coal Creek, a perennial stream with an average annual discharge of 24,200 acre-feet per year, is the largest stream in Cedar Valley. Typically, all of the water in Coal Creek is diverted for irrigation during the summer months. All surface water is consumed within the basin by irrigated crops, evapotranspiration, or recharge to the ground-water system. Ground water in Cedar Valley generally moves from primary recharge areas along the eastern margin of the basin where Coal Creek enters, to areas of discharge or subsurface outflow. Recharge to the unconsolidated basin-fill aquifer is by seepage of unconsumed irrigation water, streams, direct precipitation on the unconsolidated basin fill, and subsurface inflow from consolidated rock and Parowan Valley and is estimated to be about 42,000 acre-feet per year. Stable-isotope data indicate that recharge is primarily from winter precipitation. The chloride mass-balance method indicates that recharge may be less than 42,000 acre-feet per year, but is considered a rough approximation because of limited chloride concentration data for precipitation and Coal Creek. Continued declining water levels indicate that recharge is not sufficient to meet demand. Water levels in many areas are at or close to historic lows. In 2000, withdrawal from wells was estimated to be 36,000 acre-feet per year. About 4,000 acre-feet per year are estimated to discharge to evapotranspiration or as subsurface outflow. Prior to large-scale ground-water development, ground-water discharge by evapotranspiration and discharge to springs was much larger. Ground water along the eastern margin of the valley between Cedar City and Enoch is unsuitable for domestic use because of high dissolved-solids and nitrate concentrations. The predominant ions of Ca and SO4 in this area indicate dissolution of gypsum in the Markagunt Plateau to the east. Data collected during this study were compared to historic data; there is no evidence to indicate deterioration in ground-water quality. The spatial distribution of ground water with high nitrate concentration does not appear to be migrating beyond its previously known extent. No single source can be identified as the cause for elevated nitrate concentrations in ground water. Low nitrogen-15 values north of Cedar City indicate a natural geologic source. Higher nitrogen-15 values toward the center of the basin and associated hydrologic data indicate probable recharge from waste-water effluent. Excess dissolved nitrogen gas and low nitrate concentrations in shallow ground water indicate that denitrification is occurring in some areas. A computer ground-water flow model was developed to simulate flow in the unconsolidated basin fill. The method of determining recharge from irrigation was changed during the calibration process to incorporate more areal and temporal variability. In general, the model accurately simulates water levels and water-level fluctuations and can be considered an adequate tool to help determine the valley-wide effects on water levels of additional ground-water withdrawals and changes in water use. The model was used to simulated water-level changes caused by projecting current withdrawal rates, increased withdrawal rates, and a 10-year drought. Water levels declined 20 to 275 feet in the southern and central parts of the valley and less than 20 feet north of Enoch

Utah↗

Simulated transport and biodegradation of chlorinated ethenes in a fractured dolomite aquifer near Niagara Falls, New York

Leakage of trichloroethene (TCE) from a neutralization pond at a former manufacturing facility near Niagara Falls, N.Y. during 1950-87 into the Guelph Formation of the Lockport Group, a fractured dolomite aquifer, created a plume of TCE and its metabolites that, by 1990, extended about 4,300 feet south of the facility. A smaller plume of dense, nonaqueous-phase liquids (DNAPL) probably serves as a continuing source of TCE. The presence of the TCE metabolites cis-1,2-dichloroethene (DCE), vinyl chloride (VC), and ethene in the plume, and the results of previous laboratory microcosm studies, indicate that the TCE is being degraded by naturally occurring microorganisms. Biodegradation rates of TCE and its metabolites were estimated through simulation with BIOMOC, a solute-transport model that represents multispecies reactions through Monod kinetics. A fracture zone in the Guelph Formation was represented as a porous medium containing an extensive, 3-foot thick layer with several interconnected fractures; this layer is bounded above and below by subhorizontal stratigraphic contacts. The Monod reaction constants were estimated through nonlinear regression to minimize the difference between computed concentrations of TCE and its metabolites, and the concentrations measured before and during 5 years of pump-and-treat remediation. Transport simulations indicated that, by April 1998, the chlorinated ethene plume had reached a dynamic equilibrium between the rate of TCE dissolution and the rate of removal through pumping and biodegradation. Biodegradation of chlorinated ethenes at the site can be simulated as first-order reactions because the concentrations are generally less than the half-saturation constants estimated for Monod kinetics (320 mg/L for TCE, 10 mg/L for DCE, and 1 mg/L for VC). Computed degradation rates are proportional to the estimated ground-water velocity, which could vary by more than an order magnitude at the site, as indicated by the estimated range of fracture porosity--3 to 0.3 percent. Half-lives corresponding to first-order rate constants estimated for the lower velocity (5 to 15 feet per day) ranged from 21 to 25 days for TCE, 170 to 230 days for DCE, and 18 to 23 days for VC. Chlorinated ethene concentrations of April 1998 were better reproduced when the TCE source was represented as a constant concentration than as a constant flux, because the latter predicted that the plume would dissipate after 5 years of pump-and-treat remediation. This result indicates that the rate of TCE dissolution is not limited by the mass of TCE in the DNAPL plume. Simulation of diffusion by the transport model MOC3D indicated that concentrations of these contaminants within the rock matrix surrounding the fracture zone were relatively unchanged after 5 years of pump-and-treat remediation. The principal sources of uncertainty in the prediction of biodegradation rates and of the fate of chlorinated ethenes at the site are the fracture porosity and DNAPL mass in the aquifer.

New York↗

Detection of underground voids in Ohio by use of geophysical methods

Geophysical methods are generally classified as electrical, potential field, and seismic methods. Each method type relies on contrasts of physical properties in the subsurface. Forward models based on the physical properties of air- and water-filled voids within common geologic materials indicate that several geophysical methods are technically feasible for detection of subsurface voids in Ohio, but ease of use and interpretation varies widely between the methods. Ground-penetrating radar is the most rapid and cost-effective method for collection of subsurface data in areas associated with voids under roadways. Electrical resistivity, gravity, or seismic reflection methods have applications for direct delineation of voids, but data-collection and analytical procedures are more time consuming. Electrical resistivity, electromagnetic, or magnetic methods may be useful in locating areas where conductive material, such as rail lines, are present in abandoned underground coal mines. Other electrical methods include spontaneous potential and very low frequency (VLF); these latter two methods are considered unlikely candidates for locating underground voids in Ohio. Results of ground-penetrating radar surveys at three highway sites indicate that subsurface penetration varies widely with geologic material type and amount of cultural interference. Two highway sites were chosen over abandoned underground coal mines in eastern Ohio. A third site in western Ohio was chosen in an area known to be underlain by naturally occurring voids in lime stone. Ground-penetrating radar surveys at Interstate 470, in Belmont County, Ohio, indicate subsurface penetration of less than 15 feet over a mined coal seam that was known to vary in depth from 0 to 40 feet. Although no direct observations of voids were made, anomalous areas that may be related to collapse structures above voids were indicated. Cultural interference dominated the radar records at Interstate 70, Guernsey County, Ohio, where coal was mined under the site at a depth of about 50 feet. Interference from overhead powerlines, the field vehicle, and guardrails complicated an interpretation of the radar records where the depth of penetration was estimated to be less than 5 feet. Along State Route 33, in Logan County, Ohio, bedding planes and structures possibly associated with dissolution of limestone were profiled with ground-penetrating radar. Depth of penetration was estimated to be greater than 50 feet.

Water-Resources Investigations Report↗

Use of environmental tracers to evaluate ground-water age and water-quality trends in a buried-valley aquifer, Dayton area, southwestern Ohio

Chlorofluorocarbons (CFC method) and tritium and helium isotopes (3H-3He method) were used as environmental tracers to estimate ground-water age in conjunction with efforts to develop a regional ground-water flow model of the buried-valley aquifer in the Dayton area, southwestern Ohio. This report describes results of CFC and water-quality sampling, summarizes relevant aspects of previously published work, and describes the use of 3H-3He ages to characterize temporal trends in ground-water quality of the buried-valley aquifer near Dayton, Ohio. Results of CFC sampling indicate that approximately 25 percent of the 137 sampled wells were contaminated with excess CFC's that rendered the ground water unsuitable for age dating. Evaluation of CFC ages obtained for the remaining samples indicated that the CFC compounds used for dating were being affected by microbial degradation. The degradation occurred under anoxic conditions that are found in most parts of the buried-valley aquifer. As a result, ground-water ages derived by the CFC method were too old and were inconsistent with measured tritium concentrations and independently derived 3H-3He ages. Limited data indicate that dissolved methane may play an important role in the degradation of the CFC's. In contrast, the 3H-3He technique was found to yield ground-water ages that were chemically and hydrologically reasonable. Ground-water ages derived by the 3H-3He technique were compared to values for selected water- quality characteristics to evaluate temporal trends in ground-water quality in the buried- valley aquifer. Distinct temporal trends were not identified for pH, alkalinity, or calcium and magnesium because of rapid equilibration of ground-water with calcite and dolomite in aquifer sediments. Temporal trends in which the amount of scatter and the number of outlier concentrations increased as ground-water age decreased were noted for sodium, potassium, boron, bromide, chloride, ammonia, nitrate, phosphate, sulfate, and organic carbon. Elevated concentrations of these constituents in shallow ground water are probably related to human activities. Temporal trends in which concentrations declined as ground-water age increased may reflect natural processes that reduce constituent concentrations to low levels. For example, the absence of nitrate detections in ground water recharged before 1980 may indicate natural removal of nitrate by bacterially mediated denitrification. Temporal trends observed for dissolved oxygen, iron, nitrate and silica indicate that these constituents may help identify recently (post-1990) recharged ground water.

Ohio↗

Community estimate of global glacier mass changes from 2000 to 2023

Glaciers are indicators of ongoing anthropogenic climate change 1 . Their melting leads to increased local geohazards 2 , and impacts marine 3 and terrestrial 4 , 5 ecosystems, regional freshwater resources 6 , and both global water and energy cycles 7 , 8 . Together with the Greenland and Antarctic ice sheets, glaciers are essential drivers of present 9 , 10 and future 11 , 12 , 13 sea-level rise. Previous assessments of global glacier mass changes have been hampered by spatial and temporal limitations and the heterogeneity of existing data series 14 , 15 , 16 . Here we show in an intercomparison exercise that glaciers worldwide lost 273 ± 16 gigatonnes in mass annually from 2000 to 2023, with an increase of 36 ± 10% from the first (2000–2011) to the second (2012–2023) half of the period. Since 2000, glaciers have lost between 2% and 39% of their ice regionally and about 5% globally. Glacier mass loss is about 18% larger than the loss from the Greenland Ice Sheet and more than twice that from the Antarctic Ice Sheet 17 . Our results arise from a scientific community effort to collect, homogenize, combine and analyse glacier mass changes from in situ and remote-sensing observations. Although our estimates are in agreement with findings from previous assessments 14 , 15 , 16 at a global scale, we found some large regional deviations owing to systematic differences among observation methods. Our results provide a refined baseline for better understanding observational differences and for calibrating model ensembles 12 , 16 , 18 , which will help to narrow projection uncertainty for the twenty-first century 11 , 12 , 18 .

Nature↗