USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modeling Report”

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

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

At least 685 records · Page 38Linked to original sources

MODFLOW-2000, The U.S. Geological Survey modular ground-water model: User guide to modularization concepts and the ground-water flow process

MODFLOW is a computer program that numerically solves the three-dimensional ground-water flow equation for a porous medium by using a finite-difference method. Although MODFLOW was designed to be easily enhanced, the design was oriented toward additions to the ground-water flow equation. Frequently there is a need to solve additional equations; for example, transport equations and equations for estimating parameter values that produce the closest match between model-calculated heads and flows and measured values. This report documents a new version of MODFLOW, called MODFLOW-2000, which is designed to accommodate the solution of equations in addition to the ground-water flow equation. This report is a user's manual. It contains an overview of the old and added design concepts, documents one new package, and contains input instructions for using the model to solve the ground-water flow equation.

Open-File Report↗

Mapping habitat for multiple species in the Desert Southwest

Many utility scale renewable energy projects are currently proposed across the Mojave Ecoregion. Agencies that manage biological resources throughout this region need to understand the potential impacts of these renewable energy projects and their associated infrastructure (for example, transmission corridors, substations, access roads, etc.) on species movement, genetic exchange among populations, and species’ abilities to adapt to changing environmental conditions. Understanding these factors will help managers’ select appropriate project sites and possibly mitigate for anticipated effects of management activities. We used species distribution models to map habitat for 15 species across the Mojave Ecoregion to aid regional land-use management planning. Models were developed using a common 1 × 1 kilometer resolution with maximum entropy and generalized additive models. Occurrence data were compiled from multiple sources, including VertNet ( http://vertnet.org/ ), HerpNET ( http://www.herpnet.org ), and MaNIS ( http://manisnet.org ), as well as from internal U.S. Geological Survey databases and other biologists. Background data included 20 environmental covariates representing terrain, vegetation, and climate covariates. This report summarizes these environmental covariates and species distribution models used to predict habitat for the 15 species across the Mojave Ecoregion.

Mojave Ecoregion↗

Geohydrology of Storage Unit III and a combined flow model of the Santa Barbara and foothill ground-water basins, Santa Barbara County, California

The city of Santa Barbara pumps most of its ground water from the Santa Barbara and Foothill ground-water basins. The Santa Barbara basin is subdivided into two storage units: Storage Unit I and Storage Unit III. The Foothill basin and Storage Unit I of the Santa Barbara basin have been studied extensively and ground-water flow models have been developed for them. In this report, the geohydrology of the Santa Barbara ground- water basin is described with a special emphasis on Storage Unit III in the southwestern part of the basin. The purposes of this study were to summarize and evaluate the geohydrology of Storage Unit III and to develop an areawide model of the Santa Barbara and Foothill basins that includes the previously unmodeled Storage Unit III. Storage Unit III is in the southwestern part of the city of Santa Barbara. It is approximately 3.5 miles long and varies in width from about 2,000 feet in the southeast to 4,000 feet in the north-west. Storage Unit III is composed of the Santa Barbara Formation and overlying alluvium. The Santa Barbara Formation (the principal aquifer) consists of Pleistocene and Pliocene(?) unconsolidated marine sand, silt, and clay, and it has a maximum saturated thickness of about 160 feet. The alluvium that overlies the Santa Barbara Formation has a maximum saturated thickness of about 140 feet. The storage unit is bounded areally by faults and low-permeability deposits and is underlain by rocks of Tertiary age. The main sources of recharge to Storage Unit III are seepage from Arroyo Burro and infiltration of precipitation. Most of the recharge occurs in the northwest part of the storage unit, and ground water flows toward the southeast along the unit's long axis. Lesser amounts of recharge may occur as subsurface flow from the Hope Ranch subbasin and as upwelling from the underlying Tertiary rocks. Discharge from Storage Unit III occurs as pumpage, flow to underground drains, underflow through alluvium in the vicinity of Arroyo Burro across the Lavigia Fault, evapotranspiration, and underflow to the Pacific Ocean. The faults that bound Storage Unit III generally are considered to be effective barriers to the flow of ground water. Interbasin ground-water flow occurs where deposits of younger alluvium along stream channels cross faults. Ground-water quality in Storage Unit III deposits varies with location and depth. Upward leakage of poor-quality water from the underlying Tertiary rocks occurs in the storage unit, and such leakage can be influenced by poor well construction or by heavy localized pumping. The highest dissolved-solids concentration (4,710 milligrams per liter) in ground water resulting from this upward leakage is found in the coastal part of the storage unit. The ground-water system was modeled as two horizontal layers. In the Foothill basin and Storage Unit I the layers are separated by a confining bed. The upper layer represents the upper producing zone and the shallow zone near the coast. The lower layer represents the lower producing zone. In general, the faults in the study area were assumed to be no-flow boundaries, except for the offshore fault that forms the southeast boundary; the southeast boundary was simulated as a general-head boundary. The Storage Unit III model was combined with the preexisting Storage Unit I and Foothill basin models, using horizontal flow barriers, to form an areawide model. The areawide model was calibrated by simulating steady-state predevelopment conditions and transient conditions for 1978-92. The nonpumping steady- state simulation was used to verify that the calibrated model yielded physically reasonable results for predevelopment conditions. The calibrated areawide model calculates water levels in Storage Unit III that are within 10 feet of measured water levels at all sites of comparison. In addition, the model adequately simulates water levels in the Storage Unit I and Foothill basin areas. A total of 33,430 acre-feet of water was pumped from the study area during the simulation period. Model results indicate that 2,833 acre-feet came from storage and 5, 332 acre-feet crossed the general-head boundary from the ocean, thus indicating that seawater intrusion could occur. A sensitivity analysis indicates that, in general, the model is most sensitive to changes in transmissivity and total recharge.

California↗

Listings of model values for the simulation of ground-water flow in the Cimarron River alluvium and terrace deposits from Freedom to Guthrie, Oklahoma

This report contains MODFLOW input and output listings for the simulation of ground-water flow in alluvium and terrace deposits associated with the Cimarron River from Freedom to Guthrie, Oklahoma. These values are to be used in conjuction with the report, 'Geohydrology of alluvium and terrace deposits of the Cimarron River from Freedom to Guthrie, Oklahoma,' by G.P. Adams and D.L. Bergman, published as U.S. Geological Survey Water-Resources Investigatons Report 95-4066. The simulation used a digital ground-water flow model and was evaluated by a management and statistical program.

Open-File Report↗

Forecasting selenium discharges to the San Francisco Bay-Delta Estuary: Ecological effects of a proposed San Luis Drain extension

Selenium discharges to the San Francisco Bay-Delta Estuary (Bay-Delta) could change significantly if federal and state agencies (1) approve an extension of the San Luis Drain to convey agricultural drainage from the western San Joaquin Valley to the North Bay (Suisun Bay, Carquinez Strait, and San Pablo Bay); (2) allow changes in flow patterns of the lower San Joaquin River and Bay-Delta while using an existing portion of the San Luis Drain to convey agricultural drainage to a tributary of the San Joaquin River; or (3) revise selenium criteria for the protection of aquatic life or issue criteria for the protection of wildlife. Understanding the biotransfer of selenium is essential to evaluating effects of selenium on Bay-Delta ecosystems. Confusion about selenium threats to fish and wildlife stem from (1) monitoring programs that do not address specific protocols necessary for an element that bioaccumulates; and (2) failure to consider the full complexity of the processes that result in selenium toxicity. Past studies show that predators are more at risk from selenium contamination than their prey, making it difficult to use traditional methods to predict risk from environmental concentrations alone. This report presents an approach to conceptualize and model the fate and effects of selenium under various load scenarios from the San Joaquin Valley. For each potential load, progressive forecasts show resulting (1) water-column concentration; (2) speciation; (3) transformation to particulate form; (4) particulate concentration; (5) bioaccumulation by invertebrates; (6) trophic transfer to predators; and (7) effects on those predators. Enough is known to establish a first-order understanding of relevant conditions, biological response, and ecological risks should selenium be discharged directly into the North Bay through a conveyance such as a proposed extension of the San Luis Drain. The approach presented here, the Bay-Delta selenium model, determines the mass, fate, and effects of selenium released to the Bay-Delta through use of (1) historical land-use, drainage, alluvial-fill, and runoff databases; (2) existing knowledge concerning biogeochemical reactions and physiological parameters of selenium (e.g., speciation, partitioning between dissolved and particulate forms, and bivalve assimilation efficiency); and (3) site-specific data mainly from 1986 to 1996 for clams and bottom-feeding fish and birds. Selenium load scenarios consider effluents from North Bay oil refineries and discharges of agricultural drainage from the San Joaquin Valley to enable calculation of (a) a composite freshwater endmember selenium concentration at the head of the estuary; and (b) a selenium concentration at a selected seawater location (Carquinez Strait) as a foundation for modeling. Analysis of selenium effects also takes into account the mode of conveyance for agricultural drainage (i.e., the San Luis Drain or San Joaquin River); and flows of the Sacramento River and San Joaquin River on a seasonal or monthly basis. Load scenarios for San Joaquin Valley mirror predictions made since 1955 of a worsening salt (and by inference, selenium) build-up exacerbated by an arid climate and massive irrigation. The reservoir of selenium in the San Joaquin Valley is sufficient to provide loading at an annual rate of approximately 42,500 pounds of selenium to a Bay-Delta disposal point for 63 to 304 years at the lower range of projections presented here, even if influx of selenium from the California Coast Ranges could be curtailed. Disposal of wastewaters on an annual basis outside of the San Joaquin Valley may slow the degradation of valley resources, but drainage alone cannot alleviate the salt and selenium build-up in the San Joaquin Valley, at least within a century. Load scenarios also show the different proportions of selenium loading to the Bay-Delta. Oil refinery loads from 1986 to 1992 ranged from 8.5 to 20 pounds of selenium per day; with treatment and cleanup, loads decreased to 3.0 pounds of selenium per day in 1999. In contrast, San Joaquin Valley agricultural drainage loads disposed of in a San Luis Drain extension could range from 45 to 117 pounds of selenium per day across a set of historical and future conditions. Components of this valley-wide load include five source subareas (i.e., Grassland, Westlands, Tulare, Kern, and Northern) defined by water and drainage management. Loads vary per subarea mainly because of proximity of the subarea to geologic sources of selenium and irrigation history. Loads from the Sacramento River, depending on flow conditions, range from 0.8 to 10 pounds of selenium per day. Loads from the San Joaquin River vary depending on restoration and flow conditions, which are considered. A consistent picture of ecological risk emerges under modeled selenium discharges from a proposed San Luis Drain extension. The threat to the estuary is greatest during low flow seasons and critically dry years. Where selenium undergoes reactions typical of low flow or longer residence time, highly problematic bioaccumulation in prey (food) is forecast. Surf scoter, greater and lesser scaup, and white sturgeon appear to be most at risk because these Bay-Delta predators feed on deposit and filter-feeding bivalves. Recent findings add Sacramento splittail and Dungeness crab to that list. During the low flow season of critically dry years, forecasted selenium concentrations in water, particulate matter, prey (diet), and predator tissue exceed guidelines with a high certainty of producing adverse effects under the most likely load scenario from a proposed San Luis Drain extension. High flows afford some protection under certain conditions in modeled San Joaquin River scenarios. However, meeting a combined goal of releasing a specific load during maximum flows and keeping selenium concentrations in the river below a certain objective to protect against bioaccumulation may not always be attainable. Management of the San Joaquin River on a constant concentration basis also could create problematic bioaccumulation during a wet year, especially during the low flow season, because high flows translate to high loads that are not always offset by seasonal river inflows. Prior to refinery cleanup, selenium contamination was sufficient to threaten reproduction in key species within the Bay-Delta ecosystems and human health advisories were posted based on selenium concentrations in tissues of diving ducks. During this time, selenium concentrations in the Bay-Delta were well below the most stringent recommended water quality criterion [1 microgram per liter (1 µg/L)]. Enhanced biogeochemical transformations to bioavailable particulate selenium and efficient bioaccumulation by bivalves characterized the system. If these biogeochemical conditions continue to prevail and agricultural selenium sources replace or exceed refinery sources, ecological forecasts suggest the risk of adverse effects will be difficult to eliminate under an out-of-valley resolution to the selenium problem. The Bay-Delta selenium model presented here is a systematic approach for conducting forecasts of the ecological effects from selenium on aquatic food webs. It is a new tool that links and models the major processes leading from loads through consumer organisms to predators. It also is a feasible approach for site-specific analysis and could provide a framework for developing new protective selenium foodweb guidelines and predator criteria. Model components that help ensure understanding ecosystems and the basis of environmental protection are (1) contaminant concentrations and speciation in sources, such as particulate material, that most influence bioavailability; (2) bioaccumulation models that calculate concentrations in diet, specifically in bivalves of the Bay-Delta that act as sensitive indicators of selenium contamination; (3) food-web type that determines what animals are threatened and when; and (4) multiple media concentrations (water, particulate material, and tissue of prey and predators) that, in-combination, determine risk or hazard.

California↗

Simulation of Ground-Water Flow in the Middle Rio Grande Basin Between Cochiti and San Acacia, New Mexico

This report describes a three-dimensional, finite difference, ground-water-flow model of the Santa Fe Group aquifer system within the Middle Rio Grande Basin between Cochiti and San Acacia, New Mexico. The aquifer system is composed of the Santa Fe Group of middle Tertiary to Quaternary age and post-Santa Fe Group valley and basin-fill deposits of Quaternary age. Population increases in the basin since the 1940's have caused dramatic increases in ground-water withdrawals from the aquifer system, resulting in large ground-water-level declines. Because the Rio Grande is hydraulically connected to the aquifer system, these ground-water withdrawals have also decreased flow in the Rio Grande. Concern about water resources in the basin led to the development of a research plan for the basin focused on the hydrologic interaction of ground water and surface water (McAda, D.P., 1996, Plan of study to quantify the hydrologic relation between the Rio Grande and the Santa Fe Group aquifer system near Albuquerque, central New Mexico: U.S. Geological Survey Water-Resources Investigations Report 96-4006, 58 p.). A multiyear research effort followed, funded and conducted by the U.S. Geological Survey and other agencies (Bartolino, J.R., and Cole, J.C., 2002, Ground-water resources of the Middle Rio Grande Basin, New Mexico: U.S. Geological Survey Circular 1222, 132 p.). The modeling work described in this report incorporates the results of much of this work and is the culmination of this multiyear study. The purpose of the model is (1) to integrate the components of the ground-water-flow system, including the hydrologic interaction between the surface-water systems in the basin, to better understand the geohydrology of the basin and (2) to provide a tool to help water managers plan for and administer the use of basin water resources. The aquifer system is represented by nine model layers extending from the water table to the pre-Santa Fe Group basement rocks, as much as 9,000 feet below the NGVD 29. The horizontal grid contains 156 rows and 80 columns, each spaced 3,281 feet (1 kilometer) apart. The model simulates predevelopment steady-state conditions and historical transient conditions from 1900 to March 2000 in 1 steady-state and 52 historical stress periods. Average annual conditions are simulated prior to 1990, and seasonal (winter and irrigation season) conditions are simulated from 1990 to March 2000. The model simulates mountain-front, tributary, and subsurface recharge; canal, irrigation, and septic-field seepage; and ground-water withdrawal as specified-flow boundaries. The model simulates the Rio Grande, riverside drains, Jemez River, Jemez Canyon Reservoir, Cochiti Lake, riparian evapotranspiration, and interior drains as head-dependent flow boundaries. Hydrologic properties representing the Santa Fe Group aquifer system in the ground-water-flow model are horizontal hydraulic conductivity, vertical hydraulic conductivity, specific storage, and specific yield. Variable horizontal anisotropy is applied to the model so that hydraulic conductivity in the north-south direction (along model columns) is greater than hydraulic conductivity in the east-west direction (along model rows) over much of the model. This pattern of horizontal anisotropy was simulated to reflect the generally north-south orientation of faulting over much of the modeled area. With variable horizontal anisotropy, horizontal hydraulic conductivities in the model range from 0.05 to 60 feet per day. Vertical hydraulic conductivity is specified in the model as a horizontal to vertical anisotropy ratio (calculated to be 150:1 in the model) multiplied by the horizontal hydraulic conductivity along rows. Specific storage was estimated to be 2 x 10-6 per foot in the model. Specific yield was estimated to be 0.2 (dimensionless). A ground-water-flow model is a tool that can integrate the complex interactions of hydrologic boundary conditions, aquifer materials

Water-Resources Investigations Report↗

Annual dissolved nitrite plus nitrate and total phosphorous loads for the Susquehanna, St. Lawrence, Mississippi-Atchafalaya, and Columbia River basins, 1968-2004

Annual stream-water loads were calculated near the outlet of four of the larger river basins (Susquehanna, St. Lawrence, Mississippi-Atchafalaya, and Columbia) in the United States for dissolved nitrite plus nitrate (NO2 + NO3) and total phosphorus using LOADEST load estimation software. Loads were estimated for the period 1968-2004; although loads estimated for individual river basins and chemical constituent combinations typically were for shorter time periods due to limitations in data availability. Stream discharge and water-quality data for load estimates were obtained from the U.S. Geological Survey (USGS) with additional stream discharge data for the Mississippi-Atchafalaya River Basin from the U.S. Army Corps of Engineers. The loads were estimated to support national assessments of changes in stream nutrient loads that are periodically conducted by Federal agencies (for example, U.S. Environmental Protection Agency) and other water- and land-resource organizations. Data, methods, and results of load estimates are summarized herein; including World Wide Web links to electronic ASCII text files containing the raw data. The load estimates are compared to dissolved NO2 + NO3 loads for three of the large river basins from 1971 to 1998 that the USGS provided during 2001 to The H. John Heinz III Center for Science, Economics and the Environment (The Heinz Center) for a report The Heinz Center published during 2002. Differences in the load estimates are the result of using the most up-to-date monitoring data since the 2001 analysis, differences in how concentrations less than the reporting limit were handled by the load estimation models, and some errors and exclusions in the 2001 analysis datasets (which resulted in some inaccurate load estimates).

Open-File Report↗

Proceedings of a workshop on fish habitat suitability index models

One of the habitat-based methodologies for impact assessment currently in use by the U.S. Fish and Wildlife Service is the Habitat Evaluation Procedures (HEP) (U.S. Fish and Wildlife Service 1980). HEP is based on the assumption that the quality of an area as wildlife habitat at a specified target year can be described by a single number, called a Habitat Suitability Index (HSI). An HSI of 1.0 represents optimum habitat: an HSI of 0.0 represents unsuitable habitat. The verbal or mathematical rules by which an HSI is assigned to an area are called an HSI model. A series of Habitat Suitability Index (HSI) models, described by Schamberger et al. (1982), have been published to assist users in applying HEP. HSI model building approaches are described in U.S. Fish and Wildlife Service (1981). One type of HSI model described in detail requires the development of Suitability Index (SI) graphs for habitat variables believed to be important for the growth, survival, standing crop, or other measure of well-being for a species. Suitability indices range from 0 to 1.0, with 1.0 representing optimum conditions for the variable. When HSI models based on suitability indices are used, habitat variable values are measured, or estimated, and converted to SI's through the use of a Suitability Index graph for each variable. Individual SI's are aggregated into an HSI. Standard methods for testing this type of HSI model did not exist at the time the studies reported in this document were performed. A workshop was held in Fort Collins, Colorado, February 14-15, 1983, that brought together biologists experienced in the use, development, and testing of aquatic HSI models, in an effort to address the following objectives: (1) review the needs of HSI model users; (2) discuss and document the results of aquatic HSI model tests; and (3) provide recommendations for the future development, testing, modification, and use of HSI models. Individual presentations, group discussions, and group decision techniques were used to develop and present information at the meeting. A synthesis of the resulting concepts, results, and recommendations follows this preface. Subsequent papers describe individual tests of selected HSI models. Most of the tests involved comparison of values from HSI models or Suitability index (SI) curves with standing crop, as required contractually. Time and budget constraints generally limited tests to the use of data previously collected for other purposes. These proceedings are intended to help persons responsible for the development, testing, or use of HSI models by increasing their understanding of potential uses and limitations of testing procedures and models based on aggregated Suitability Indices. Problems encountered when testing HSI models are described, model performance during tests is documents, and recommendations for future model development and testing presented by the participants are listed and interpreted.

Biological Report↗

Documentation of models describing relations between continuous real-time and discrete water-quality constituents in the Little Arkansas River, south-central Kansas, 1998–2019

Data were collected at two monitoring sites along the Little Arkansas River in south-central Kansas that bracket most of the easternmost part of the Equus Beds aquifer. The data were used as part of the city of Wichita’s aquifer storage and recovery project to evaluate source water quality. The U.S. Geological Survey, in cooperation with the City of Wichita, has continued to monitor the water quality of these sites through 2019 to update previously published regression-based models using continuously measured physicochemical properties and discretely sampled water-quality constituents of interest. The purpose of this report is to provide an update of the previously published linear regression models that have been used to continuously compute estimates of water-quality constituent concentrations or densities at these two sites. Water-quality constituent model updates include those for dissolved and suspended solids, suspended-sediment concentration, hardness, alkalinity, primary ions (bicarbonate, calcium, sodium, chloride, and sulfate), nutrients (total Kjeldahl nitrogen and total phosphorus), total organic carbon, indicator bacteria ( Escherichia coli and fecal coliform bacteria), a trace element (arsenic), and a pesticide (atrazine). Regression analyses were used to develop surrogate models that related continuously measured physicochemical properties, streamflow, and seasonal components to discretely sampled water-quality constituent concentrations or densities. Specific conductance was an explanatory variable for dissolved solids, primary ions, and atrazine. Turbidity was an explanatory variable for total suspended solids and sediment, nutrients, total organic carbon, and indicator bacteria. Streamflow and water temperature were explanatory variables for dissolved arsenic. Seasonal components were included as explanatory variables for atrazine models. The amount of variance explained by most of the updated models was within 5 percent of previously published models.

Kansas↗

Precipitation runoff modeling system (PRMS) as part of an integrated hydrologic model for the Osage Nation, northeastern Oklahoma, 1915–2014

Executive Summary The Osage Nation lacks a comprehensive tribal water plan to describe the quality and quantity of water resources in the Osage Nation, a 2,304-square-mile (mi 2 ) area of rolling pastures, tallgrass prairie, and mixed woodlands in northeastern Oklahoma. A tribal water plan can be used to help manage the sustainable development of surface and groundwater resources, thereby helping to provide a better future for the Osage Nation and their neighbors, while preserving water resources for the benefit of the surrounding environment and future generations. To help meet these goals and contribute to increased knowledge of the quantity and quality of water resources and the hydrologic processes and factors affecting those resources, the U.S. Geological Survey (USGS) in cooperation with the Osage Nation began studies to evaluate the surface-water and groundwater resources of the Osage Nation. An important component of these studies is the development and application of numerical models to improve quantification and understanding of the hydrologic system. These models are needed to estimate and quantify the effects of historical and potential future water resource development for the Osage Nation. This report describes the development and application of a precipitation-runoff model, the Osage Nation watershed model (ONWM). The ONWM is needed as a component of the Osage Nation integrated hydrologic model (ONIHM). At the time of this study, the ONIHM was being developed using the USGS computer software MODFLOW-One Water Hydrologic Flow Model (MODFLOW-OWHM). The intended use of the ONIHM is to simulate all surface-water and groundwater components of the hydrologic system for a 2,905-mi 2 study area centered on the Osage Nation. The ONWM was developed using the USGS Precipitation Runoff Modeling System, version 4 (PRMS-IV) computer software, also referred to as PRMS in this report, for an 8,343-mi 2 study area in northeastern Oklahoma and southeastern Kansas, centered on and including the areas of the Osage Nation and the ONIHM. The ONWM is to be used as part of the ONIHM to provide a direct coupling with spatially and temporally varying daily climate conditions affecting the ONIHM study area. As an integral part of the ONIHM, the ONWM (1) simulates the inflow boundary conditions from tributary basins in the region outside and surrounding the ONIHM area; (2) provides estimates of spatially and temporally distributed precipitation, air temperature, potential evapotranspiration (PET), actual evapotranspiration (ET), soil moisture, recharge, and streamflow in the ONIHM area; and (3) provides a preliminary water budget for the ONIHM area and the surrounding region, including tributary drainage basins outside of and next to the ONIHM. The specific objectives of this study were to use the ONWM to (1) provide a systematic inventory of the historical distribution of water inflows from precipitation (rain or snow) falling on the land surface and flowing through the surface-water network, (2) provide a historical context of the variability and spatial and temporal distribution of these waters, and (3) provide estimates of water inflows and potential observations to the ONIHM. The application of the ONWM as a component of the ONIHM is needed for planned simulations using the ONIHM to improve the understanding of the hydrologic system and to develop a fully comprehensive water budget, including the use and movement of water across the landscape, in the surface-water network, and in groundwater aquifers under historical and potential future conditions.

Oklahoma↗

2023 Coastal master plan: Model improvement plan, ICM-wetlands, vegetation, and soil

As part of the model improvement effort for the 2023 Coastal Master Plan, the wetland processes captured by the morphology and vegetation models used during previous master plans were reevaluated to assess how Integrated Compartment Model (ICM) subroutines could be improved. This process considered technical reviews, comments, and suggested improvements provided by model developers, advisory groups, and other experts during previous master plan cycles. The availability of new data and information that could be used to make model improvements was also considered. In many cases, the team considered and tested multiple options or approaches. As a result of this effort, recommended improvements are provided here. The improvements recommended to be included in the 2023 Coastal Master Plan include: adjusting marsh collapse thresholds, refining organic matter accretion calculations, developing an unstructured grid for modeling vegetation, improving flotant marsh and forested wetlands algorithms, creating and applying an updated map of existing vegetation, adjusting model code, and updating the submerged aquatic vegetation (SAV) module. This report describes the team’s work through a series of 7 distinct activities to identify and test options for model improvements to ensure the updated ICM used for the 2023 Coastal Master Plan appropriately captures ecological and morphological processes observed in Coastal Louisiana. As appropriate, relevant literature and data are discussed. Test runs to evaluate how changes influence model outputs are also documented. A final list of recommended updates, taking into account consideration of all options and results from test runs, is summarized at the end of the report. A later report will describe the final ICM-LAVegMod and ICM-Morph subroutines for the 2023 Coastal Master Plan, detailing the updates that have been incorporated.

Louisiana↗

Development and calibration of a ground-water flow model for the Sparta Aquifer of southeastern Arkansas and north-central Louisiana and simulated response to withdrawals, 1998-2027

The Sparta aquifer, which consists of the Sparta Sand, in southeastern Arkansas and north-central Louisiana is a major water resource and provides water for municipal, industrial, and agricultural uses. In recent years, the demand in some areas has resulted in withdrawals from the Sparta aquifer that substantially exceed replenishment of the aquifer. Considerable drawdown has occurred in the potentiometric surface forming regional cones of depression as water is removed from storage by withdrawals. These cones of depression are centered beneath the Grand Prairie area and the cities of Pine Bluff and El Dorado in Arkansas, and Monroe in Louisiana. The rate of decline for hydraulic heads in the aquifer has been greater than 1 foot per year for more than a decade in much of southern Arkansas and northern Louisiana where hydraulic heads are now below the top of the Sparta Sand. Continued hydraulic-head declines have caused water users and managers alike to question the ability of the aquifer to supply water for the long term. Concern over protecting the Sparta aquifer as a sustainable resource has resulted in a continued, cooperative effort by the Arkansas Soil and Water Conservation Commission, U.S. Army Corps of Engineers, and the U.S. Geological Survey to develop, maintain, and utilize numerical ground-water flow models to manage and further analyze the ground-water system. The work presented in this report describes the development and calibration of a ground-water flow model representing the Sparta aquifer to simulate observed hydraulic heads, documents major differences in the current Sparta model compared to the previous Sparta model calibrated in the mid-1980's, and presents the results of three hypothetical future withdrawal scenarios. The current Sparta model-a regional scale, three-dimensional numerical ground-water flow model-was constructed and calibrated using available hydrogeologic, hydraulic, and water-use data from 1898 to 1997. Significant changes from the previous model include grid rediscretization of the aquifer, extension of the active model area northward beyond the Cane River Formation facies change, and representation of model boundaries. The current model was calibrated with the aid of parameter estimation, a nonlinear regression technique, combined with trial and error parameter adjustment using a total of 795 observations from 316 wells over 4 different years-1970, 1985, 1990, and 1997. The calibration data set provides broad spatial and temporal coverage of aquifer conditions. Analysis of the residual statistics, spatial distribution of residuals, simulated compared to observed hydrographs, and simulated compared to observed potentiometric surfaces were used to analyze the ability of the calibrated model to simulate aquifer conditions within acceptable error. The calibrated model has a root mean square error of 18 feet for all observations, an improvement of more than 12 feet from the previous model. The current Sparta model was used to predict the effects of three hypothetical withdrawal scenarios on hydraulic heads over the period 1998-2027 with one of those extended indefinitely until equilibrium conditions were attained, or steady state. In scenario 1a, withdrawals representing the time period from 1990 to 1997 was held constant for 30 years from 1998 to 2027. Hydraulic heads in the middle of the cone of depression centered on El Dorado decreased by 10 feet from the 1997 simulation to 222 feet below NGVD of 1929 in 2027. Hydraulic heads in the Pine Bluff cone of depression showed a greater decline from 61 feet below NGVD of 1929 to 78 feet below NGVD of 1929 in the center of the cone. With these same withdrawals extended to steady state (scenario 1b), hydraulic heads in the Pine Bluff cone of depression center declined an 2 Development and Calibration of a Ground-Water Flow Model for the Sparta Aquifer of Southeastern Arkansas and North-Central Louisiana and Simulated Response to Withdrawa

Water-Resources Investigations Report↗

Spectral characteristics and mapping of lithium-rich playas in the Basin and Range Province, western United States

This report presents a comprehensive set of prospectivity maps for lithium in brines located within the playas of the Basin and Range Province of western United States. The prospectivity assessment is based on criteria established in the descriptive deposit model outlined in a previous U.S. Geological Survey report, with a focus on regional remote sensing methodologies. The playas within the study area have been systematically ranked according to several key factors, including size, the presence and abundance of source rocks, vegetation as an indicator of water availability, documented prospects, and data products derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor. Notably, the deposits in the playa of Clayton Valley make it stand out as the sole lithium brine-producing playa in the United States. The ASTER sensor provides optimal regional coverage, characterized by a 60-kilometer swath width and unique multispectral capabilities that are not available with other remote sensing instruments. This report and the accompanying maps are intended to serve as resources for further multidisciplinary assessments of lithium or other elements or commodities where prospectivity can be inferred from the presence of proxies such as alteration minerals.

Arizona, Baja California, California, Chihuahua, I↗

Assessing potential effects of highway runoff on receiving-water quality at selected sites in Oregon with the Stochastic Empirical Loading and Dilution Model (SELDM)

In 2012, the U.S. Geological Survey and the Oregon Department of Transportation began a cooperative study to demonstrate use of the Stochastic Empirical Loading and Dilution Model (SELDM) for runoff-quality analyses in Oregon. SELDM can be used to estimate stormflows, constituent concentrations, and loads from the area upstream of a stormflow discharge site, from the site of interest and in the receiving waters downstream of the discharge. SELDM also can be used to assess the potential effectiveness of best management practices (BMP) for mitigating potential effects of runoff in receiving waters. Nominally, SELDM is a highway-runoff model, but it is well suited for analysis of runoff from other land uses as well. This report provides case studies and examples to demonstrate stochastic-runoff modeling concepts and to demonstrate application of the model. Basin characteristics from six Oregon highway study sites were used to demonstrate various applications of the model. The highway catchment and upstream basin drainage areas of these study sites ranged from 3.85 to 11.83 acres and from 0.16 to 6.56 square miles, respectively. The upstream basins of two sites are urbanized, and the remaining four sites are less than 5 percent impervious. SELDM facilitates analysis by providing precipitation, pre-storm streamflow, and other variables by region or from hydrologically similar sites. In Oregon, there can be large variations in precipitation and streamflow among nearby sites. Therefore, spatially interpolated geographic information system data layers containing storm-event precipitation and pre-storm streamflow statistics specific to Oregon were created for the study using Kriging techniques. Concentrations and loads of cadmium, chloride, chromium, copper, iron, lead, nickel, phosphorus, and zinc were simulated at the six Oregon highway study sites by using statistics from sites in other areas of the country. Water‑quality datasets measured at hydrologically similar basins in the vicinity of the study sites in Oregon were selected and compiled to estimate stormflow-quality statistics for the upstream basins. The quality of highway runoff and some upstream stormflow constituents were simulated by using statistical moments (average, standard deviation, and skew) of the logarithms of data. Some upstream stormflow constituents were simulated by using transport curves, which are relations between stormflow and constituent concentrations. Stochastic analyses were done by using SELDM to demonstrate use of the model and to illustrate the types of information that stochastic analyses may provide: 1. An analysis was done to demonstrate use of dilution factors as an initial reconnaissance tool for comparing relative risk among sites. 2. An analysis of hardness-dependent, water-quality criteria was done to illustrate the effects of variations in hardness and flow on the application and interpretation of such criteria. This analysis shows that hardness-dependent criteria can vary by an order of magnitude among storm events because hardness is diluted by stormflows. 3. An analysis of uncertainties in input and output values was done to demonstrate that properly selected robust datasets are needed to represent conditions at a site of interest. This analysis shows that the rate of water-quality exceedances that are measured or simulated may depend on sample size and the luck of the draw. 4. An analysis was done to demonstrate that SELDM and other Monte Carlo models may generate extreme values from input statistics, which may or may not be feasible based on physicochemical or hydrological limits. 5. An analysis of BMP modeling methods was done to demonstrate use of the model for estimating treatment requirements for meeting water-quality objectives. 6. An analysis of the use of grab sampling and nonstochastic upstream modeling methods was done to evaluate the potential effects on modeling outcomes. Additional analyses using surrogate water-quality datasets for the upstream basin and highway catchment were provided for six Oregon study sites to illustrate the risk-based information that SELDM will produce. These analyses show that the potential effects of highway runoff on receiving-water quality downstream of the outfall depends on the ratio of drainage areas (dilution), the quality of the receiving water upstream of the highway, and the concentration of the criteria of the constituent of interest. These analyses also show that the probability of exceeding a water-quality criterion may depend on the input statistics used, thus careful selection of representative values is important.

Oregon↗

Near real-time updating of pager loss estimates

Initial alerts by PAGER (Prompt Assessment of Global Earthquakes for Response) within minutes following an earthquake include several uncertainties, mainly due to potential inaccuracies in location, depth, fault delineation, and shaking estimates. We enhance an updating framework by incorporating early reports of fatalities within the first 24 hours, or so, of an earthquake to update PAGER’s overall fatality estimates and its resulting alert level. Though initial loss reports by officials or the media are uncertain and often undercount the eventual reported impacts, their temporal evolution provides predictive constraints for the PAGER model. The proposed framework helps capture these in a systematic way to minimize potential large fluctuations in PAGER alerts as ShakeMap (the USGS product which estimates how an area is affected by an earthquake) gets updated in the early hours after an earthquake. The new framework also accounts for uncertainties associated with early fatality reports as well as PAGER model-related uncertainties in order to improve the overall impact forecast. This updating framework improves the loss estimate and alert level to the correct level within the first 24 hours even when the initial estimation from PAGER is assumed to be off by two levels of alert, which is plausible due to potential over- or under-estimation of the PAGER model. While test results are very encouraging, our future work aims at implementation of operational PAGER model updating, which entails additional challenges in acquiring useful data, estimating their credibility, and developing rigorously tested operational code and protocols

Conference Paper↗

SutraGUI, a graphical-user interface for SUTRA, a model for ground-water flow with solute or energy transport

This report describes SutraGUI, a flexible graphical user-interface (GUI) that supports two-dimensional (2D) and three-dimensional (3D) simulation with the U.S. Geological Survey (USGS) SUTRA ground-water-flow and transport model (Voss and Provost, 2002). SutraGUI allows the user to create SUTRA ground-water models graphically. SutraGUI provides all of the graphical functionality required for setting up and running SUTRA simulations that range from basic to sophisticated, but it is also possible for advanced users to apply programmable features within Argus ONE to meet the unique demands of particular ground-water modeling projects. SutraGUI is a public-domain computer program designed to run with the proprietary Argus ONE? package, which provides 2D Geographic Information System (GIS) and meshing support. For 3D simulation, GIS and meshing support is provided by programming contained within SutraGUI. When preparing a 3D SUTRA model, the model and all of its features are viewed within Argus 1 in 2D projection. For 2D models, SutraGUI is only slightly changed in functionality from the previous 2D-only version (Voss and others, 1997) and it provides visualization of simulation results. In 3D, only model preparation is supported by SutraGUI, and 3D simulation results may be viewed in SutraPlot (Souza, 1999) or Model Viewer (Hsieh and Winston, 2002). A comprehensive online Help system is included in SutraGUI. For 3D SUTRA models, the 3D model domain is conceptualized as bounded on the top and bottom by 2D surfaces. The 3D domain may also contain internal surfaces extending across the model that divide the domain into tabular units, which can represent hydrogeologic strata or other features intended by the user. These surfaces can be non-planar and non-horizontal. The 3D mesh is defined by one or more 2D meshes at different elevations that coincide with these surfaces. If the nodes in the 3D mesh are vertically aligned, only a single 2D mesh is needed. For nonaligned meshes, two or more 2D meshes of similar connectivity are used. Between each set of 2D meshes (and model surfaces), the vertical space in the 3D mesh is evenly divided into a user-specified number of layers of finite elements. Boundary conditions may be specified for 3D models in SutraGUI using a variety of geometric shapes that may be located freely within the 3D model domain. These shapes include points, lines, sheets, and solids. These are represented by 2D contours (within the vertically-projected Argus ONE view) with user-defined elevations. In addition, boundary conditions may be specified for 3D models as points, lines, and areas that are located exactly within the surfaces that define the model top and the bottoms of the tabular units. Aquifer properties may be specified separately for each tabular unit. If the aquifer properties vary vertically within a unit, SutraGUI provides the Sutra_Z function that can be used to specify such variation.

Open-File Report↗

Lithostratigraphy, geophysics, biostratigraphy, and strontium-isotope stratigraphy of the surficial aquifer system of eastern Collier County and northern Monroe County, Florida

In 1997, ten cores were drilled in eastern Collier County and northern Monroe County, within the limits of the Big Cypress National Preserve. These cores represent a continuation of the study of seven cores in western Collier County begun in 1996 and reported in Weedman and others (1997) and Edwards and others (1998). This joint U.S. Geological Survey and Florida Geological Survey project is designed to acquire subsurface geologic and hydrologic data in southwest Florida to extend current ground-water models, thereby expanding the utility of these models for land and water management. In this report we describe the lithostratigraphy, geophysical logging, sedimentological analysis, dinocyst biostratigraphy, and strontium-isotope stratigraphy of these ten cores. The three geophysical logs (natural gamma-ray, induction conductivity, and neutron porosity) assumed to be related to formation lithology and water quality show that a number of clay-rich zones are present in all of the boreholes, and that pore-water conductivity increases with depth. The clay-rich zones are confirmed by visual examination of core material and sedimentological analysis. The relative transmissivity calculated at 10-foot-thick intervals shows that in six of the boreholes, high values are associated with the shallow aquifer in the 0-40 ft interval. Two of the boreholes (the most northerly and the most easterly) showed relatively higher values of transmissivity in permeable zones at or somewhat below 100 ft in depth. Core geology and logs indicate that the deeper aquifers are not more permeable than similar deeper zones in the other boreholes, but rather that the shallow aquifer appears to be less permeable in these two coreholes. The Arcadia (?) Formation was only penetrated in the deepest core where it is late Miocene in age. The Peace River Formation was penetrated in all but the two westernmost cores. It yields a late Miocene age, based on both dinocysts and strontium-isotope stratigraphy. The top is an irregular surface. Age and stratigraphic relations suggest that the upper part of the Peace River and lower part of the unnamed formation are at least partially equivalent laterally. The unnamed formation was recovered in every core. It is thinnest in the northernmost core and thickest to the west. Ages calculated from strontium isotopes range from 6.9 to 4.6 million years ago (late Miocene to early Pliocene). The top of the unnamed formation is deepest to the north and it becomes shallower to the southwest. The Tamiami Formation also was recovered in every core and consistently yields early Pliocene ages; it yields late Pliocene ages near the top in two cores. The age and lateral relations strongly suggest that the lower part of the Tamiami Formation and the upper part of the unnamed formation are lateral facies of each other. The Fort Thompson (?) Formation, Miami Limestone, and undifferentiated siliciclastic sediments and limestone at the very top of the cores were not dated.

Florida↗

Temperature and solute-transport simulation in streamflow using a Lagrangian reference frame

A computer program for simulating one-dimensional, unsteady temperature and solute transport in a river has been developed and documented for general use. The solution approach to the convective-diffusion equation uses a moving reference frame (Lagrangian) which greatly simplifies the mathematics of the solution procedure and dramatically reduces errors caused by numerical dispersion. The solution procedure has the further advantages, relative to conventional Eulerian solution schemes, of being easy to understand in the physical sense, of being extremely stable numerically, and of providing an accounting system which is very useful for model calibration. The model documentation is presented as a series of four programs of increasing complexity. The conservative transport model can be used to route a single conservative substance, such as dye, through a reach of a river. The simplified temperature model is used to predict water temperature in rivers, either with or without thermal loading, when few meteorological data are available. Only equilibrium temperature and windspeed are required. It is suggested that air temperature can be used to approximate equilibrium temperature. The complete temperature model is highly accurate but requires rather complete meteorological data. Finally, the 10-parameter model can be used to route as many as 10 interacting constituents through a river reach. The mathematical description of the interaction between the constituents, which does not need to be linear, is generally up to the user to supply. An example problem is solved for a three-parameter system involving temperature, dissolved oxygen, and biochemical oxygen demand. For simplicity, all models are developed and presented assuming steady non-uniform flow. Generalization of the models to allow unsteady flow is extremely simple, involving the addition of no more than 18 cards to the program deck. The report is concluded by describing this generalization for any of the models. Before using the models with unsteady flow, a flow model must be used to calculate and store the necessary flow data at each cross section and time step. Such a flow model is available and documented.

Water-Resources Investigations Report↗