Search USGS⌕ Search

SEARCH · Search USGS

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

Estimation of unaltered daily mean streamflow at ungaged streams of New York, excluding Long Island, water years 1961-2010

The lakes, rivers, and streams of New York State provide an essential water resource for the State. The information provided by time series hydrologic data is essential to understanding ways to promote healthy instream ecology and to strengthen the scientific basis for sound water management decision making in New York. The U.S. Geological Survey, in cooperation with The Nature Conservancy and the New York State Energy Research and Development Authority, has developed the New York Streamflow Estimation Tool to estimate a daily mean hydrograph for the period from October 1, 1960, to September 30, 2010, at ungaged locations across the State. The New York Streamflow Estimation Tool produces a complete estimated daily mean time series from which daily flow statistics can be estimated. In addition, the New York Streamflow Estimation Tool provides a means for quantitative flow assessments at ungaged locations that can be used to address the objectives of the Clean Water Act—to restore and maintain the chemical, physical, and biological integrity of the Nation’s waters. The New York Streamflow Estimation Tool uses data from the U.S. Geological Survey streamflow network for selected streamgages in New York (excluding Long Island) and surrounding States with shared hydrologic boundaries, and physical and climate basin characteristics to estimate the natural unaltered streamflow at ungaged stream locations. The unaltered streamflow is representative of flows that are minimally altered by regulation, diversion, or mining, and other anthropogenic activities. With the streamflow network data, flow-duration exceedance probability equations were developed to estimate unaltered streamflow exceedance probabilities at an ungaged location using a methodology that equates streamflow as a percentile from a flow-duration curve for a particular day at a hydrologically similar reference streamgage with streamflow as a percentile from the flow-duration curve for the same day at an ungaged location. The reference streamgage is selected using map correlation, a geostatistical method in which variogram models are developed that correlate streamflow at one streamgage with streamflows at all other locations in the study area. Regression equations used to predict 17 flow-duration exceedance probabilities were developed to estimate the flow-duration curves at ungaged locations for New York using geographic information system-derived basin characteristics. A graphical user interface, with an integrated spreadsheet summary report, has been developed to estimate and display the daily mean streamflows and statistics and to evaluate different water management or water withdrawal scenarios with the estimated monthly data. This package of regression equations, U.S. Geological Survey streamgage data, and spreadsheet application produces an interactive tool to estimate an unaltered daily streamflow hydrograph and streamflow statistics at ungaged sites in New York. Among other uses, the New York Streamflow Estimation Tool can assist water managers with permitting water withdrawals, implementing habitat protection, estimating contaminant loads, or determining the potential affect from chemical spills.

New York↗

Scientific results of the Hydrate-01 Stratigraphic Test Well Program, Western Prudhoe Bay Unit, Alaska North Slope

The United States Department of Energy, the MH21-S Research Consortium of Japan, and the United States Geological Survey are collaborating to enable gas hydrate scientific drilling and extended-duration reservoir response testing on the Alaska North Slope. To feasibly execute such a test, a location is required that is accessible from existing roads and gravel pads and that can be occupied without disrupting ongoing industry operations. A review of potential locations meeting these criteria determined the likely occurrence of gas hydrate in two fine-grained marginal-marine sands of Tertiary age in the vicinity of the inactive “Kuparuk State 7-11-12” exploration pad in the western Prudhoe Bay Unit (PBU). Existing well and seismic data for that site were insufficient to preclude the potential for free gas occurrence within the deeper (and most prospective) target sand. Therefore, with support from the PBU Working Interest Owners, Alaska Department of Natural Resources, and Petrotechnical Resources Alaska, the Hydrate-01 Stratigraphic Test Well (STW) was drilled in December 2018 to confirm the suitability of the site for future gas hydrate scientific testing. The Hydrate-01 well was successfully drilled to −3290 ft (1003 m) subsea vertical depth at a bottom hole location of approximately 900 ft (∼275 m) east of the surface location. The drilling program featured acquisition of a full suite of logging while drilling data, the collection of side-wall pressure cores, and the installation of distributed temperature and distributed acoustic sensor fiber-optic cables. The log data acquired confirmed the occurrence of gas hydrate at high saturation in two target sands. Integrated evaluation of log and sidewall core data provide petrophysical and geomechanical property information that allow for potential reservoir response to depressurization to be simulated. The deeper “B1 sand” is deemed to be most favorable for reservoir response testing as a result of confirmed gas hydrate occurrence in sediments of high intrinsic permeability, location within 100 ft (30 m) of the base of gas hydrate stability, and minimal risk for direct communication with permeable water-bearing (hydrate-free) zones. The shallower “D1 sand” provides a secondary target that is differentiated by colder in situ temperatures and the interpreted direct hydraulic communication to a lower section of non-hydrate-bearing, water-saturated sand. The Hydrate-01 log data also confirm the occurrence of at least one sub-seismic fault in close proximity to the B1 sand reservoir. To better image the distribution of the gas-hydrate-bearing reservoir sections and associated faults, a three-dimensional (3D) vertical seismic profile was conducted in early 2019 using the distributed acoustic sensors installed as part of the Hydrate-01 STW completion. Detailed two-dimensional (2D) and 3D geologic models have been constructed to enable numerical simulations to inform the planning for potential future scientific tests of reservoir response to depressurization at the site.

Alaska↗

Suwannee river basin and estuary integrated science workshop: September 22-24, 2004 Cedar Key, Florida

In response to the growing number of environmental concerns in the mostly pristine Suwannee River Basin and the Suwannee River Estuary system, the States of Florida and Georgia, the Federal government, and other local organizations have identified the Suwannee River as an ecosystem in need of protection because of its unique biota and important water resources. Organizations with vested interests in the region formed a coalition, the Suwannee Basin Interagency Alliance (SBIA), whose goals are to promote coordination in the identification, management, and scientific knowledge of the natural resources in the basin and estuary. To date, an integrated assessment of the physical, biological, and water resources has not been completed. A holistic, multi-disciplinary approach is being pursued to address the research needs in the basin and estuary and to provide supportive data for meeting management objectives of the entire ecosystem. The USGS is well situated to focus on the larger concerns of the basin and estuary by addressing specific research questions linking water supply and quality to ecosystem function and health across county and state boundaries. A strategic plan is being prepared in cooperation with Federal, State, and local agencies to identify and implement studies to address the most compelling research issues and management questions, and to conduct fundamental environmental monitoring studies. The USGS, Suwannee River Water Management District and the Florida Marine Research Institute are co-sponsoring this scientific workshop on the Suwannee River Basin and Estuary to: Discuss current and past research findings, Identify information gaps and research priorities, and Develop an action plan for coordinated and relevant research activities in the future. This workshop builds on the highly successful basin-wide conference sponsored by the Suwannee Basin Interagency Alliance that was held three years ago in Live Oak, Florida. This years workshop will focus on identifying information needs and priorities and developing partnerships. The USGS is seeking to define the role of the USGS Florida Integrated Science Center (FISC) in conducting integrated research in the Suwannee River Basin, and to establish a cooperative program with other agencies. Participants interested in river, floodplain, springs, estuary, or basin-wide issues are encouraged to attend. Topics for this years workshop include: Water quality and geochemistry: nutrient enrichment, reduction of nutrient loading to ground water, contaminants, and land use, Hydrogeology: interactions among ground water, surface water and ecosystem, modeling, and baseline mapping, Ecosystem dynamics: structure, process, species, and habitats (estuarine, riverine, floodplain, and wetland), and Information management: data sharing, database development, geographic information system (GIS), and basin-wide models.

Open-File Report↗

Assessment experimental semivariogram uncertainty in the presence of a polynomial drift

The semivariogram, which measures the spatial variability between experimental data, is generally used as a structural input in all two-point geostatistical procedures. However, in most geoscience applications, experimental semivariograms are usually computed from a limited number of sparsely spaced measurements, which results in uncertainty associated with the semivariance values estimated for a specified number of lags. More importantly, considering a spatial variable modelled by a nonstationary random field, uncertainty is not only in the experimental semivariogram of the residuals, but also in the coefficients of the drift model estimated from the available experimental data. Therefore, when assessing the reliability of an experimental semivariogram (or estimated semivariances) in the nonstationary case, both aforementioned uncertainties should be taken into account. The aim of this paper is to extend the “Generalised Bootstrap” procedure to the nonstationary model by propagating the uncertainty associated with the estimated drift coefficients into the uncertainty in the experimental semivariogram of the residuals. The proposed methodology is demonstrated in a case study using abundant geophysical measurements characterised by a nonstationary random function. Two scenarios are evaluated in the case study: (1) it is assumed that the drift coefficients can be estimated without any uncertainty, and (2) uncertainty of the drift coefficients is taken into account. We have explained the methodology that allows to assess the uncertainty of the semivariogram lag estimates in the presence of the drift in the mean. Considering the second scenario, uncertainty is obviously larger than the case where uncertainty of the drift in the mean is ignored. This evaluation should be considered in applications where the data is often rather limited, such as subsurface hydrology (i.e. porosity, transmissivity), soil science (i.e. heavy metal content, soil moisture) and mining (i.e. scoping or pre-feasibility stage of the project). In fact, modern geostatistics should provide not only the semivariogram estimates but also estimation of its uncertainty.

Natural Resources Research↗

Prefabricated Zero Ascend Omnispecies (ZAO) modular fish passage modules using advanced manufacturing techniques

The Zero Ascend Omnispecies (ZAO) fish passage attraction module “ZAO Attractor” is a novel system designed to facilitate upstream passage for a broad spectrum of fish species, with a particular focus on American shad and river herring (alewife and/or blueback herring). Combining the use of prefabricated, modular components with advanced manufacturing, the ZAO aims to provide a cost-effective and scalable solution for fish passage attractions across various hydropower sites. The system features a surface level entry, so fish do not need to climb. In addition, it is designed to be adjusted and repositioned to create the most effective attraction flows by managing the outflows from hydropower turbines. Budget Period 1 (BP1) focused on the design and modeling of the ZAO Attractor and demonstrated the engineering feasibility of a modular, prefabricated, hydrodynamically efficient structure with surface level entry and passage. During BP1, the hypothesis was that certain structural adaptations, such as the addition of long fins, could direct turbine outflows to create a desirable downstream flow. Simulations, however, were inconclusive, revealing the complex nature of fish preferences for hydraulic conditions. It was discovered that there were limited published studies on the direct correlation between quantified flow patterns and fish attraction, emphasizing the need for controlled studies to observe and measure these effects. In Budget Period 2 (BP2), controlled studies with live fish were prioritized to address these complexities. A partially scaled ZAO Attractor was evaluated in laboratory conditions using native river water to assess its capability to attract and guide fish. Despite challenges with flume flow conditions and the impact of flume dimensions on flow directionality and velocity, the studies provided insights into fish behavior and preferences, particularly the requirement for a consistent flow stream for migration cues. Adjustments to the fin structure demonstrated potential effectiveness in guiding fish toward the module. However, the hydraulic conditions within the flume environment did not accurately reflect typical in-river scenarios, therefore the results were inconclusive with respect to the optimal system configuration for in-river use. In the future, conducting comprehensive field studies can enhance understanding of the interactions between fish behavior and flow conditions created by the ZAO Attractor. Additionally, a design review is advised, particularly regarding the load scenarios throughout the seasons and positioning of the ZAO Attractor. These studies should continue to refine the design to achieve its intended benefits: economical and effective fish attraction and passage that can be adapted to specific site conditions while allowing for continued energy generation at sites up for re-licensure and adding hydropower generation to previously untapped resources in non-powered dams and new stream reaches. In conclusion, the ZAO fish passage attraction module presents a promising approach to fish passage at hydropower stations, with the potential to balance environmental concerns with renewable energy development. However, further field validation and optimization are necessary to ensure its effectiveness and practicality in real-world applications.

Final Technical Report↗

Modeling, mapping, and measuring the risk of freshwater invasive species across Alaska

Freshwater ecosystems of the Alaskan Arctic and Subarctic provide resources that are culturally, ecologically, and economically invaluable. Presently, these regions are relatively free of the impacts from invasive species compared to southern latitudes. To date, there have been relatively few verified introductions of aquatic invasive species (AIS) to freshwater ecosystems in Alaska. The expanding list and distribution of AIS has led to significant negative ecological and economic impacts (e.g., waterweed Elodea nuttalli ; E. canadensis and northern pike Esox Lucius introduced outside its native range in Alaska). Escalating human activity across Alaskan lands and waters, coupled with rapidly shifting environmental conditions, increases the potential for new species introductions and subsequent establishment. Creating a proactive framework for well-informed decision-making and action can improve the effectiveness of prevention efforts and bolster decision support tools that help resource managers direct limited resources. Prioritizing AIS that may be introduced and become established, as well as the locations at highest risk of invasion, is foundational to building a proactive invasive species management framework in Alaska. This project sought to identify and prioritize AIS known to be invasive in the contiguous United States, evaluate current and future habitat suitability for AIS in Alaska, and assess potential for AIS to be transported to habitats across Alaska, utilizing similar assessment methods as implemented for Bering Sea marine invasive species and non-native plants in Alaska. To accomplish this goal, the objectives of the project were to: 1) develop a formal ranked list of potential AIS to freshwater systems of Alaska; 2) assess the level of establishment risk for potential AIS by developing habitat suitability models for waterbodies across Alaska; and 3), identify potential pathways and specific vectors for high-risk AIS to invade Alaska and develop a framework for how vector analysis will be completed to understand transport risk. Overall, our goal is horizon scanning which is defined by Roy et al. (2019) as “a systematic examination of potential threats and opportunities, within a given context, and likely future developments, which are at the margin of current thinking and planning.” The scans include pathway analyses and risk screening of species present at pathway origin points, with a focus on identifying species at high risk of being introduced, becoming established, spreading, and causing harm. We refined a list of 28 AIS from a list of hundreds based on characterizations of species’ invasiveness and species’ proximity to Alaska (USGS 2020; GBIF 2022). Next, we evaluated the relative invasiveness of individual species to create an initial AIS ranking. We sought to characterize habitat suitability of AIS by selecting variables that were continental in scale, covering North America to include Alaska as well as the lower 48 states comparing natural discharge, sub-basin average terrain slope (degrees), average silt fraction, average organic carbon, lithological class, and human footprint in sub-basin in 2009. We estimated AIS habitat suitability across the entire state of Alaska using the physiological tolerances of the AIS (Appendix 2). We also evaluated pathways and vectors for the introduction of AIS (Appendix 2). Many pathways and vectors considered did not meet the criteria for Alaska or freshwater systems. Of the 28 ranked species that we categorized as very high, high, and moderate levels of invasiveness; all three risk groups included fish and mollusks (Appendix 2). One commonality of the very high-invasiveness-ranked species was the availability of Ecological Risk Screening Summary documents (USFWS, 2022) produced by U.S. Fish and Wildlife Service (USFWS), except for the goldfish ( Carassius auratus ) and the New Zealand mudsnail ( Potamopyrgus antipodarum ). The Ecological Risk Screening Summary is now available for New Zealand mudsnails. In general, fish species often ranked very high or high in invasiveness and included sportfish and aquarium fish, suggesting the importance of pathways such as aquarium trade, fishing industry, intentional (but illegal) introductions of sportfishes and aquarium fishes for establishment. The technique we used for habitat suitability models necessitated aquatic environmental datasets that were continental in scale, which was often interpolated from very coarse resolution source data layers, particularly in Alaska. Better spatial data representing aquatic environments would likely improve this approach. While the lack of introductions in Alaska and nearby provinces and states is encouraging, the lack of occurrence data for the focal species also created complications for habitat suitability modeling. Despite the challenges, the habitat suitability models indicated limited suitability for warmwater species while some species, such as Brook trout (S alvelinus fontinalis ), have high habitat suitability across Alaska no matter what threshold approach is taken. Some environmental predictors were more important than others. Specifically, the most important predictor variable, ‘frost free days,’ was critical for 15 out of 28 species as expected due to harsh winter conditions in Arctic and Subarctic regions. The second most important predictor was ‘subbasin land surface runoff’, a variable that indicates the amount of discharge and runoff, while the third most important predictor was ‘snow cover’ another indication of winter conditions. Overall, the ability to understand the effect of future climate scenarios on the establishment of AIS was challenging. A detailed dataset of freshwater temperatures and water chemistry (e.g., pH, calcium) would greatly improve the ability to predict invasiveness of freshwater species to Alaska’s ecosystems on a regional basis. Future studies may benefit from a more focused geographic scope examining a group of subbasins or a regional basin rather than the entire state. These drainages could be selected based upon the mostly likely locations of introduction pathways. The two most prevalent pathway risks for AIS are in-state transfer and stowaways/contaminants. Although there are examples of introductions from other pathways, the risk is somewhat mitigated by Alaska’s climate and regulations. However, variable application of protocols for inspection and cleaning of fishing gear, watercraft, and other similar items while traveling into Alaska as well as transferring from waterbody to waterbody within the state creates a substantial risk in introducing invasive species. We plot cumulative invasive vulnerability for all subbasins and for the top 10% of subbasins (Appendix 3).

Alaska↗

Water quality, bed-sediment quality, and simulation of potential contaminant transport in Foster Creek, Berkeley County, South Carolina, 1991-93

Foster Creek, a freshwater tidal creek in Berkeley County, South Carolina, is located in an area of potential contaminant sources from residential, commercial, light industrial, and military activities. The creek is used as a secondary source of drinking water for the surrounding Charleston area. Foster Creek meets most of the freshwater- quality requirements of State and Federal regulatory agencies, but often contains low concentrations of dissolved oxygen and has been characterized as eutrophic. Investigations of water- and bed-sediment quality were made between 1991 and 1993 to assess the effects of anthropogenic sources of contamination on Foster Creek. Low-flow surface-water samples were generally free of toxic compounds with the exception of laboratory artifacts and naturally occurring trace metals. Storm-runoff samples generally contained very low concentrations (near detection limits) of a small number of volatile and semivolatile organics and naturally occurring trace metals. Concentrations of toxic compounds in excess of current (1995) South Carolina Department of Health and Environmental Control and U.S. Environmental Protection Agency regulations were not detected in surface-water samples collected from Foster Creek. Chemical analyses of streambed sediments indicated minimal anthropogenic effects on sediment quality. The particle-tracking option of the U.S. Geological Survey one-dimensional unsteady-flow model (BRANCH) indicated that as the simulated volume of rainfall runoff increased in the Foster Creek Basin, simulated particles in Foster Creek were transported greater distances. Simulating flow through the Bushy Park Dam (also known as Back River Dam) had little effect on particle movement in Foster Creek. Simulating typical withdrawal rates at a water-supply intake resulted in a slight attraction of particles toward the intake during conditions of relatively low runoff. These withdrawals had a greater influence on particles downstream of the intake than on those upstream of the intake. Simulations confirmed earlier findings which suggested that the creek would not flush during baseflow conditions, with the exception of the lower 1-mile reach, where flushing results from tidal movements. According to the simulations, Foster Creek will fully flush if a 2-year, 7-day storm occurs. Flushing appears to be affected more by the total volume of storm runoff than by typical municipal withdrawals or tidal effects.

South Carolina↗

Assessing denitrification in groundwater using natural gradient tracer tests with 15N: In situ measurement of a sequential multistep reaction

Denitrification was measured within a nitrate‐contaminated aquifer on Cape Cod, Massachusetts, using natural gradient tracer tests with 15 N nitrate. The aquifer contained zones of relatively high concentrations of nitrite (up to 77 μM) and nitrous oxide (up to 143 μM) and has been the site of previous studies examining ground water denitrification using the acetylene block technique. Small‐scale (15–24 m travel distance) tracer tests were conducted by injecting 15 N nitrate and bromide as tracers into a depth interval that contained nitrate, nitrite, nitrous oxide, and excess nitrogen gas. The timing of the bromide breakthrough curves at down‐gradient wells matched peaks in 15 N abundance above background for nitrate, nitrite, nitrous oxide, and nitrogen gas after more than 40 days of travel. Results were simulated with a one‐dimensional transport model using linked reaction kinetics for the individual steps of the denitrification reaction pathway. It was necessary to include within the model spatial variations in background concentrations of all nitrogen oxide species. The model indicated that nitrite production (0.036–0.047 μmol N (L aquifer) −1 d −1 ) was faster than the subsequent denitrification steps (0.013–0.016 μmol N (L aquifer) −1 d −1 for nitrous oxide and 0.013–0.020 μmol N (L aquifer) −1 d −1 for nitrogen gas) and that the total rate of reaction was slower than indicated by both acetylene block tracer tests and laboratory incubations. The rate of nitrate removal by denitrification was much slower than the rate of transport, indicating that nitrate would migrate several kilometers down‐gradient before being completely consumed.

Water Resources Research↗

Potential for, and possible effects of, artificial recharge in Carson Valley, Douglas County, Nevada

Rapid population growth in Carson Valley, west- central Nevada, requires a dependable municipal water source. Artificial recharge of aquifers using available flow of the Carson River is one way to increase the amount of water in underground storage and maintain a dependable ground-water supply. Ground water can be artificially recharged by routing excess surface water or, after proper treatment, routing wastewater to infiltration basins or injection wells. Withdrawal wells would remove stored water when needed. As a first step, maps showing areas in Carson Valley with high, low, moderate and unknown potential for artificial recharge were developed on the basis of the distribution of geologic units, depth to water, specific yield, infiltration rate, and location of natural recharge and discharge. For recharge by means of infiltration, areas totaling 5,700 acres have high potential, 23,900 acres have moderate potential, and 6,200 acres have low potential. For recharge through injection, areas totaling 7,800 acres have high potential and 43,500 acres have moderate potential; 23,000 acres have unknown potential because data are lacking on subsurface conditions. A ground-water-flow model was used to assess the possible results of artificial recharge. Simulations with no accompanying ground-water withdrawal show that, when recharge by injection is simulated near the valley floor, heads in the semiconfined aquifer increase over much of the valley, floor; only about 20 percent of the recharged water is stored in the aquifer after 5 years and as much as 80 percent is lost to streamflow and evapotranspiration. When recharge is simulated on the eastern side of the valley, 80 percent of the recharged water remains in storage after 5 years. When recharge is simulated near the valley floor, more water is lost to discharge than when recharge is on the eastern side of the valley. When recharge is applied for long periods without accompanying withdrawal, recharged water moves downgradient to discharge areas. The recharge water that discharges to the surface-water system could in turn replenish base flow of the Carson River and benefit downstream users.

Water-Resources Investigations Report↗

Deep well injection of brine from Paradox Valley, Colorado: Potential major precipitation problems remediated by nanofiltration

Groundwater brine seepage into the Dolores River in Paradox Valley, Colorado, increases the dissolved solids load of the Colorado River annually by ∼2.0 × 10 8 kg. To abate this natural contamination, the Bureau of Reclamation plans to pump ∼3540 m 3 /d of brine from 12 shallow wells located along the Dolores River. The brine, with a salinity of 250,000 mg/L, will be piped to the deepest (4.9 km) disposal well in the world and injected mainly into the Mississippian Leadville Limestone. Geochemical modeling indicates, and water-rock experiments confirm, that a huge mass of anhydrite (∼1.0 × 10 4 kg/d) likely will precipitate from the injected brine at downhole conditions of 120°C and 500 bars. Anhydrite precipitation could increase by up to 3 times if the injected brine is allowed to mix with the highly incompatible formation water of the Leadville Limestone and if the Mg in this brine dolomitizes the calcite of the aquifer. Laboratory experiments demonstrate that nanofiltration membranes, which are selective to divalent anions, provide a new technology that remediates the precipitation problem by removing ∼98% of dissolved SO 4 from the hypersaline brine. The fluid pressure used (50 bars) is much lower than would be required for traditional reverse osmosis membranes because nanofiltration membranes have a low rejection efficiency (5–10%) for monovalent anions. Our results indicate that the proportion of treatable brine increases from ∼60% to >85% with the addition of trace concentrations of a precipitation inhibitor and by blending the raw brine with the effluent stream.

Water Resources Research↗

Flood characteristics of Oklahoma streams techniques for calculating magnitude and frequency of floods in Oklahoma, with compilations of flood data through 1971

The 2-, 5-, 10-, 25-, 50-, and 100-year recurrence interval floods are related to basin and climatic parameters for natural streams in Oklahoma by multiple regression techniques through the mathematical model, Q x =aA b S c P d , where Q x is peak discharge for recurrence interval x, A is contributing drainage area, S is main channel slope, P is mean annual precipitation, and a, b, c, and d are regression constants and coefficients. One equation for each recurrence interval applies statewide for all natural streams of less than 2,500 mil (6,500 km 2 ), except where manmade works, such as dams, flood-detention structures, levees, channelization, and urban development, appreciably affect flood runoff. The equations can be used to estimate flood frequency of a stream at an ungaged site if drainage area size, main channel slope, and mean annual precipitation are known. At or near gaged sites, a weighted average of the regression results and the gaging station data is recommended. Individual relations of flood magnitude to contributing drainage area are given for all or parts of the main stems of the Arkansas, Salt Fork Arkansas, Cimarron, North Canadian, Canadian, Washita, North Fork Red, and Red Rivers. Parts of some of these streams, and all of the Neosho and Verdigris Rivers are not included because the effects of. major regulation from large reservoirs cannot be evaluated within the scope of the report. Graphical relations of maximum floods of record for eastern and western Oklahoma provide a guide to maximum probable floods. A random sampling of the seasonal occurrence of floods indicated about two-thirds of all annual floods in Oklahoma occur during. April through July. Less than one-half of one percent of annual floods occur in December. A compilation of flood records at all gaging sites in Oklahoma and some selected sites in adjacent States is given in an appendix. Basin and climatic parameters and log-Pearson Type III frequency data and statistics are given for most station records. A second appendix gives a reprint of the U.S. Water Resources Council Bulletin 15 which describes procedures for fitting a log-Pearson Type III distribution to gaging station data.

Oklahoma↗

Maps of estimated nitrate and arsenic concentrations in basin-fill aquifers of the southwestern United States

Human-health concerns and economic considerations associated with meeting drinking-water standards motivated a study of the vulnerability of basin-fill aquifers to nitrate contamination and arsenic enrichment in the southwestern United States. Statistical models were developed by using the random forest classifier algorithm to predict concentrations of nitrate and arsenic across a model grid representing about 190,600 square miles of basin-fill aquifers in parts of Arizona, California, Colorado, Nevada, New Mexico, and Utah. The statistical models, referred to as classifiers, reflect natural and human-related factors that affect aquifer vulnerability to contamination and relate nitrate and arsenic concentrations to explanatory variables representing local- and basin-scale measures of source and aquifer susceptibility conditions. Geochemical variables were not used in concentration predictions because they were not available for the entire study area. The models were calibrated to assess model accuracy on the basis of measured values. Only 2 percent of the area underlain by basin-fill aquifers in the study area was predicted to equal or exceed the U.S. Environmental Protection Agency drinking-water standard for nitrate as N (10 milligrams per liter), whereas 43 percent of the area was predicted to equal or exceed the standard for arsenic (10 micrograms per liter). Areas predicted to equal or exceed the drinking-water standard for nitrate include basins in central Arizona near Phoenix; the San Joaquin Valley, the Santa Ana Inland, and San Jacinto Basins of California; and the San Luis Valley of Colorado. Much of the area predicted to equal or exceed the drinking-water standard for arsenic is within a belt of basins along the western portion of the Basin and Range Physiographic Province that includes almost all of Nevada and parts of California and Arizona. Predicted nitrate and arsenic concentrations are substantially lower than the drinking-water standards in much of the study area-about 93 percent of the area underlain by basin-fill aquifers was less than one-half the standard for nitrate as N (5.0 milligrams per liter), and 50 percent was less than one-half the standard for arsenic (5.0 micrograms per liter). The predicted concentrations and the improved understanding of the susceptibility and vulnerability of southwestern basin-fill aquifers to nitrate contamination and arsenic enrichment can be used by water managers as a qualitative tool to assess and protect the quality of groundwater resources in the Southwest.

Arizona, California, Colorado, Idaho, Nevada, New ↗

Ground-water resources and simulated effects of withdrawals in the East Shore area of Great Salt Lake, Utah

The ground-water resources in the East Shore area of Great Salt Lake, Utah, were studied to better define the ground-water system; to document changes in ground-water levels, quality, and storage; and to simulate effects of an increase in ground-water withdrawals. The East Shore aquifer system is in basin-fill deposits, and is primarily a confined system with unconfined parts near the mountain front. Recharge to and discharge from the East Shore aquifer system were estimated to average about 160,000 acre-feet per year during 1969-84, with minor amounts of water being removed from storage during that period. Major sources of ground-water recharge are seepage from surface water in natural channels and irrigation canals, and subsurface inflow from consolidated rock to the basin-fill deposits. Discharge of ground water is primarily to wells, water courses, springs, and as diffuse seepage to Great Salt Lake. Average annual surface-water inflow to the study area was estimated to be 860,000 acre-feet for the period 1969-84. Annual withdrawal of ground water for municipal and industrial use increased from about 10,000 acre-feet in 1960 to more than 30,000 acre-feet in 1980 to supply a population that increased from 175,000 in 1960 to 290,000 in 1980. Long-term trends of ground-water levels indicate a steady decline at most observation wells since 1952, despite near normal or increased precipitation since the late 1960's. Water levels declined as much as 50 feet near the principal pumping center in the east-central part of the study area. They declined as much as 35 feet more than five miles from the pumping center. The increase in withdrawals and subsequent water-level declines have caused about 700 wells within 30 square miles to cease flowing since 1954. A numerical model of the East Shore aquifer system in the Weber Delta area was constructed and calibrated using water-level data and changes in ground-water withdrawals for 1955-85. Predictive simulations were made based on doubling the 1980-84 rate of municipal and industrial withdrawals for 20 years, and using both average and below-average recharge rates. The simulations indicated water-level declines of an additional 35 to 50 feet near the principal pumping center; a decrease in natural discharge to drains, evapotranspiration, and Great Salt Lake; and a decrease in ground-water storage of 80,000 to 115,000 acre-feet after 20 years.

Utah↗

Intra-community implications of implementing multiple tsunami-evacuation zones in Alameda, California

Tsunami-evacuation planning in coastal communities is typically based on maximum evacuation zones for a single scenario or a composite of sources; however, this approach may over-evacuate a community and overly disrupt the local economy and strain emergency-service resources. To minimize the potential for future over-evacuations, multiple evacuation zones based on arrival time and inundation extent are being developed for California coastal communities. We use the coastal city of Alameda, California (USA), as a case study to explore population and evacuation implications associated with multiple tsunami-evacuation zones. We use geospatial analyses to estimate the number and type of people in each tsunami-evacuation zone and anisotropic pedestrian evacuation models to estimate pedestrian travel time out of each zone. Results demonstrate that there are tens of thousands of individuals in tsunami-evacuation zones on the two main islands of Alameda, but they will likely have sufficient time to evacuate before wave arrival. Quality of life could be impacted by the high number of government offices, schools, day-care centers, and medical offices in certain evacuation zones and by potentially high population density at one identified safe area after an evacuation. Multi-jurisdictional evacuation planning may be warranted, given that many at-risk individuals may need to evacuate to neighboring jurisdictions. The use of maximum evacuation zones for local tsunami sources may be warranted given the limited amount of available time to confidently recommend smaller zones which would result in fewer evacuees; however, this approach may also result in over-evacuation and the incorrect perception that successful evacuations are unlikely.

California↗

Vegetation classification and mapping, Vicksburg National Military Park, Mississippi

Executive Summary The National Park Service (NPS) Gulf Coast Inventory and Monitoring Network, with the support of the National Park Service Vegetation Inventory completed vegetation classification and mapping for Vicksburg National Military Park (VICK), in Warren County, Mississippi, from 2004 to 2010. VICK is one of more than 250 NPS units to be covered by the Vegetation Inventory. Methods and procedures follow those of the NPS Vegetation Inventory as of August, 2012 ( http://science.nature.nps.gov/im/inventory/veg/index.cfm ). Ecologists collected floristic and environmental data from 45 vegetation field plots and classified and described 13 plant community types corresponding to US National Vegetation Classification (USNVC) associations from these data. This classification was used to map 712 hectares (1,698 acres) of Vicksburg National Military Park, in ten map classes corresponding to vegetation associations and two non-vegetated land cover classes, from a digital orthomosaic created from 1:24,000 color infrared aerial photographs and digital elevation model (DEM) data. Methods used to delineate stands varied according to vegetation type and were a combination of incorporating stand boundaries from an existing map, new “heads-up” interpretation and modeling of digital environmental spatial data. A thematic accuracy assessment collected vegetation data from 118 field observations stratified across the vegetation classes. The overall accuracy of the vegetation map was estimated to be 60.5% at the thematic level of USNVC association and 92.4% at the thematic level of USNVC group. Estimates of areas of vegetation types, as adjusted for mapping error, were calculated for all types mapped and/or observed during the thematic accuracy assessment. An important finding of this study is that the composition of the forests of VICK is substantially different from that described by earlier, more qualitative, accounts. Implications of these findings and possible applications of the data for resource monitoring and management are presented. This report summarizes the methods and general results of the Vegetation Inventory for Vicksburg National Military Park. Appendices include descriptions of vegetation types, a field key to the types, a list of plant species, and accuracy assessment contingency tables. This report is supplemented by the project data, which include: spatial database containing map polygon attributes. Geospatial products are in Universal Transverse Mercator (UTM) projection, Zone 15, using the North American Datum of 1983. all field data (plot, observation point, and accuracy assessment point) stored in a Microsoft Access database metadata for all digital products The data and reports produced by this investigation reside at: https://irma.nps.gov/App/ and/or at: http://biology.usgs.gov/npsveg/

Mississippi↗

Evaluation of salinity and nutrient conditions in the Heart River Basin, North Dakota, 1970–2020

The Heart River Basin is predominantly an agricultural basin in western North Dakota and is approximately 3,350 square miles. The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service and the Grant County Soil Conservation District, completed a study to assess spatial and temporal patterns of water quality in the Heart River Basin. The purpose of this report is to describe the methods and results of a study to evaluate salinity and nutrients in the Heart River Basin in western North Dakota. Water-quality and streamflow data used in the study were compiled from 1970 to 2020 using the National Water Quality Monitoring Council Water Quality Portal and National Water Information System. Changes in streamflow characteristics were investigated at three sites from 1970 to 2020, and changes in water quality were investigated at four sites from 1974 to 2019. Streamflow analysis indicated decreasing streamflow from 1970 until the late 1990s followed by increasing streamflow through 2020, with the largest increase in the 7-day minimum streamflow or base flow. For the historical water-quality trend period (1974–2019), total dissolved solids, sulfate, sodium, chloride, and sodium adsorption ratio concentrations have increased since the mid-1970s through 2019. Potassium concentrations during the historical period remained mostly constant with some small fluctuations. Calcium and magnesium concentrations increased since the mid-1970s at all sites, except for a decrease at one site between 1974 and 1999. During the recent trend period (1999–2019), increasing concentrations in total dissolved solids, sulfate, sodium, chloride, calcium, magnesium, and sodium adsorption ratios were observed across the Heart River Basin. The magnitude of the increases was smaller at tributary sites compared to main-stem sites. During the recent period, potassium was mostly constant, although small (−0.9 milligram per liter or less) decreases on tributaries and minor (1.3 milligrams per liter) increases on the main-stem sites were detected. Unlike dissolved ion concentrations, significant increases in nutrient concentrations were not detected from 1999 to 2019, but nitrate plus nitrite concentrations most likely decreased upstream from Lake Tschida. Inverse modeling for period 1 (1974–99) in model zone 1 (Heart River reach from site 5 to site 6) had eight reasonable models that indicated the clay mineral-water interactions and dissolution of evaporites control the geochemistry. Results of the inverse modeling for period 2 (1999–2019) in model zone 1 also had eight reasonable models that indicated that the dissolution of evaporites was the major geochemical control. Results of the geochemical modeling for period 1 (1974–99) in model zone 2 (Heart River and Sweetbriar Creek reach from sites 20 and 21 to site 22) produced seven reasonable models, and the geochemical control of the system was the dissolution of sulfate evaporite minerals. Geochemical modeling results for period 2 (1999–2019) in model zone 2 produced 11 reasonable models and was also controlled by the dissolution of sulfate evaporite minerals. Differences between the two model zones indicated that geology controls some of the water-quality changes in the Heart River Basin. Loads were estimated for total dissolved solids, sulfate, sodium, and chloride and total phosphorus. Annual loads estimated for the Heart River from 2013 through 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) were generally greatest in 2014 and least in 2016 for total dissolved solids, sulfate, sodium, and chloride. Most of the annual loads of total dissolved solids, sulfate, sodium, and chloride are delivered in March through July in the Heart River at these sites and are likely from snowmelt and spring and summer rains. The mean annual yields of total dissolved solids and sodium from 2013 to 2020 generally were largest in Big Muddy Creek (site 18), whereas yields of sulfate and chloride were largest at Sweetbriar Creek (site 21) compared to the other selected sites in the Heart River Basin. Larger yields of total dissolved solids, sulfate, sodium, and chloride at sites located on Big Muddy Creek and Sweet Briar Creek in the lower Heart River Basin were likely a result of differences in geology and soils upstream from the selected sites. A mass balance of total dissolved solids, sulfate, sodium, and chloride was estimated for the lower Heart River Basin, specifically the reach below Lake Tschida to Mandan (site 7 to site 22). Intervening flow was the largest contributor to the dissolved ion loads in the lower Heart River Basin and is an important part of understanding the transport of dissolved ions in the basin. The intervening load can include groundwater discharge, irrigation return flow, local runoff, and input from smaller ephemeral tributaries. Tributaries in the lower Heart River Basin contributed portions of the total dissolved solids, sulfate, sodium, and chloride loads at the Heart River near Mandan (site 22) that generally were proportional to the streamflow contributions. Annual loads for total phosphorus between 2013 and 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) generally were largest in 2019 and smallest in 2016. Most of the total phosphorus loads for main-stem sites 5 and 22 were transported in March, April, and June, likely from snowmelt and early summer rains. The mean annual yields of total phosphorus for 2013–20 were largest on the main-stem site upstream from Lake Tschida (site 5) and Sweetbriar Creek (site 21), whereas the smallest yields were in Big Muddy Creek (site 18). Much of the phosphorus that enters Lake Tschida from the upper basin does not get transported downstream to the lower basin, and much of the phosphorus in the lower basin was attributed to intervening flow.

North Dakota↗

Simulation of Fish, Mud, and Crystal Lakes and the shallow ground-water system, Dane County, Wisconsin

A new MODFLOW lake package (LAK3) that simulates ground-water/lake interaction was used in simulation of Fish, Mud and Crystal Lakes?three shallow seepage lakes located in northwestern Dane County, Wis. The simulations were done to help determine the cause of increasing lake stages and provide a tool to estimate the effect of pumping water from Fish lake on future lake stages. The ground-water-flow model was developed using a telescopic-mesh refinement of the Dane and southwestern Columbia Counties regional model previously developed by the U.S. Geological Survey and the Wisconsin Geological and Natural History Survey. The parameter estimation model, UCODE, was coupled to the steadystate ground-water model to automate and optimize the calibration procedure. The steady-state model was calibrated to measured ground-water levels, Spring Creek streamflow measured at Lodi, and Fish and Crystal Lake stages. The results of the steady-state model were used as initial conditions in a transient simulation beginning in 1966 and ending in 1998. Recharge based on annual baseflow in Black Earth Creek, runoff based on measured coefficients, and precipitation and evaporation from the lake surfaces, were varied during the transient simulation. Measured Fish Lake stage was matched to simulated stage to calibrate the transient model. Model results suggest that the increase in regional ground-water recharge resulted in increased ground-water flow to the lake, which in turn resulted in increased lake stages. Simulation results of withdrawal of water from Fish Lake at 500 gallons per minute, assuming 1990?98 climatic conditions, indicate that after 1 year of pumping the stage of Fish and Mud Lakes would be reduced more than 1 foot and the stage of Crystal Lake would be reduced by less than 0.2 foot. When pumping is stopped, the lake stages would recover to near pre-pumping levels within about 3 years. When pumping is extended to 5 years, Fish and Mud Lake stage would be reduced by a maximum of 3.8 feet and Crystal Lake stage is reduced a maximum of 0.8 feet. After 4 years of recovery, Fish and Mud Lake stages are within 0.9 foot of prepumping levels and Crystal Lake stage is within 0.7 foot.

Wisconsin↗

Geochemical and isotopic composition of ground water, with emphasis on sources of sulfate, in the upper Floridan aquifer and intermediate aquifer system in southwest Florida

In southwest Florida, sulfate concentrations in water from the Upper Floridan aquifer and overlying intermediate aquifer system are commonly above 250 milligrams per liter (the drinking water standard), particularly in coastal areas. Possible sources of sulfate include dissolution of gypsum from the deeper part of the Upper Floridan aquifer or the middle confining unit, saltwater in the aquifer, and saline waters from the middle confining unit and Lower Floridan aquifer. The sources of sulfate and geochemical processes controlling ground-water composition were evaluated for the Peace and Myakka River Basins and adjacent coastal areas of southwest Florida. Samples were collected from 63 wells and a saline spring, including wells finished at different depth intervals of the Upper Floridan aquifer and intermediate aquifer system at about 25 locations. Sampling focused along three ground-water flow paths (selected based on a predevelopment potentiometric-surface map). Ground water was analyzed for major ions, selected trace constituents, dissolved organic carbon, and stable isotopes (delta deuterium, oxygen-18, carbon-13 of inorganic carbon, and sulfur-34 of sulfate and sulfide); the ratio of strontium-87 to strontium-86 was analyzed for waters along one of the flow paths. Chemical and isotopic data indicate that dedolomitization reactions (gypsum and dolomite dissolution and calcite precipitation) control the chemical composition of water in the Upper Floridan aquifer in inland areas. This is confirmed by mass-balance modeling between wells in the shallowest interval in the aquifer along the flow paths. However, gypsum occurs deeper in the aquifer than these wells. Upwelling of sulfate-rich water that previously dissolved gypsum in deeper parts of the aquifer is a more likely source of sulfate than gypsum dissolution in shallow parts of the aquifer. This deep ground water moves to shallower zones in the aquifer discharge area. Saltwater from the Upper Floridan aquifer has not dissolved significant amounts of gypsum compared to fresher water in the aquifer. This is consistent with a shallow seawater source for the saltwater, rather than a deeper source from the underlying middle confining unit or Lower Floridan aquifer, which would have elevated sulfate concentrations. Ion exchange and dolomitization may be important reactions for saltwater in the aquifer. According to geochemical modeling, the freshwater end member for water in the saltwater mixing zone in the southwestern part of the study area is not upgradient water from the Upper Floridan aquifer that dissolved gypsum. Instead, this water appears to be isolated from the regional freshwater flow system and may be part of a more localized flow system. The chemical and isotopic composition of water in the intermediate aquifer system is controlled by differences in extent of reactions with aquifer minerals, upward leakage from the Upper Floridan aquifer, and saltwater mixing. In inland areas, water generally is characterized by relatively low sulfate concentrations (less than 250 milligrams per liter) and differences in extent of carbonate mineral dissolution. Some inland waters have elevated chloride concentrations, which may be related to evaporation prior to recharge. In coastal Sarasota County and in isolated inland areas, water from the intermediate aquifer system has high sulfate concentrations characteristic of dedolomitization waters from the Upper Floridan aquifer. The chemical and isotopic composition of these waters is controlled by upward leakage from the Upper Floridan aquifer, which naturally occurs in the discharge area but may be locally enhanced by pumping or interconnection of wells open to both aquifer systems. In western Charlotte County, the waters are dominated by sodium and chloride, and their compositions are consistent with mixing between saltwater and inland intermediate aquifer system water that has not been influenced by discharge from the

Florida↗