Search USGS⌕ Search

SEARCH · Search USGS

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 217 records · Page 12Linked to original sources

An inventory of three-dimensional geologic models—U.S. Geological Survey, 2004–22

A database of spatial footprints and characteristics of three-dimensional geological models that were constructed by the U.S. Geological Survey between 2004 and 2022 was compiled as part of ongoing development of subsurface geologic information by the USGS National Cooperative Geologic Mapping Program. This initial inventory resulted in the compilation of 38 three-dimensional geological models that vary widely in their spatial extent, the type and purpose of the model, the number of subsurface units characterized by the model, and the software platforms used to create the model. This Data Report provides the scientific rationale and explanation of the contents of a companion USGS digital data release of spatial data and attributes associated with each three-dimensional model.

Data Report↗

Persistent slip rate discrepancies in the eastern California (USA) shear zone

Understanding fault slip rates in the eastern California shear zone (ECSZ) using GPS geodesy is complicated by potentially overlapping strain signals due to many sub-parallel strike-slip faults and by inconsistencies with geologic slip rates. The role of fault system geometry in describing ECSZ deformation may be investigated with total variation regularization, which algorithmically determines a best-fitting geometry from an initial model with numerous faults, constrained by a western United States GPS velocity field. The initial dense model (1) enables construction of the first geodetically constrained block model to include all ECSZ faults with geologic slip rates, allowing direct geologic-geodetic slip rate comparisons, and (2) permits fault system geometries with many active faults that are analogous to distributed interseismic deformation. Beginning with 58 ECSZ blocks, a model containing 10 ECSZ blocks is most consistent with geologic slip rates, reproducing five of 11 within their reported uncertainties. The model fits GPS observations with a mean residual velocity of 1.5 mm/yr. Persistent geologic-geodetic slip rate discrepancies occur on the Calico and Garlock faults, on which we estimate slip rates of 7.6 mm/yr and <2 mm/yr, respectively, indicating that inconsistencies between geology and geodesy may be concentrated on or near these faults and are not due to pervasive distributed deformation in the region. Discrepancies may in part be due to postseismic relaxation following the A.D. 1992 M w 7.3 Landers and 1999 M w 7.1 Hector Mine earthquakes. Otherwise, resolving geologic-geodetic discrepancies would require as much as 11.4 mm/yr of off-fault deformation within <10 km of the main ECSZ faults, with ∼5 mm/yr concentrated near the Calico fault.

California↗

Simulation of the Regional Ground-Water-Flow System and Ground-Water/Surface-Water Interaction in the Rock River Basin, Wisconsin

A regional, two-dimensional, areal ground-water-flow model was developed to simulate the ground-water-flow system and ground-water/surface-water interaction in the Rock River Basin. The model was developed by the U.S. Geological Survey (USGS), in cooperation with the Rock River Coalition. The objectives of the regional model were to improve understanding of the ground-water-flow system and to develop a tool suitable for evaluating the effects of potential regional water-management programs. The computer code GFLOW was used because of the ease with which the model can simulate ground-water/surface-water interactions, provide a framework for simulating regional ground-water-flow systems, and be refined in a stepwise fashion to incorporate new data and simulate ground-water-flow patterns at multiple scales. The ground-water-flow model described in this report simulates the major hydrogeologic features of the modeled area, including bedrock and surficial aquifers, ground-water/surface-water interactions, and ground-water withdrawals from high-capacity wells. The steady-state model treats the ground-water-flow system as a single layer with hydraulic conductivity and base elevation zones that reflect the distribution of lithologic groups above the Precambrian bedrock and a regionally significant confining unit, the Maquoketa Formation. In the eastern part of the Basin where the shale-rich Maquoketa Formation is present, deep ground-water flow in the sandstone aquifer below the Maquoketa Formation was not simulated directly, but flow into this aquifer was incorporated into the GFLOW model from previous work in southeastern Wisconsin. Recharge was constrained primarily by stream base-flow estimates and was applied uniformly within zones guided by regional infiltration estimates for soils. The model includes average ground-water withdrawals from 1997 to 2006 for municipal wells and from 1997 to 2005 for high-capacity irrigation, industrial, and commercial wells. In addition, the model routes tributary base flow through the river network to the Rock River. The parameter-estimation code PEST was linked to the GFLOW model to select the combination of parameter values best able to match more than 8,000 water-level measurements and base-flow estimates at 9 streamgages. Results from the calibrated GFLOW model show simulated (1) ground-water-flow directions, (2) ground-water/surface-water interactions, as depicted in a map of gaining and losing river and lake sections, (3) ground-water contributing areas for selected tributary rivers, and (4) areas of relatively local ground water captured by rivers. Ground-water flow patterns are controlled primarily by river geometries, with most river sections gaining water from the ground-water-flow system; losing sections are most common on the downgradient shore of lakes and reservoirs or near major pumping centers. Ground-water contributing areas to tributary rivers generally coincide with surface watersheds; however the locations of ground-water divides are controlled by the water table, whereas surface-water divides are controlled by surface topography. Finally, areas of relatively local ground water captured by rivers generally extend upgradient from rivers but are modified by the regional flow pattern, such that these areas tend to shift toward regional ground-water divides for relatively small rivers. It is important to recognize the limitations of this regional-scale model. Heterogeneities in subsurface properties and in recharge rates are considered only at a very broad scale (miles to tens of miles). No account is taken of vertical variations in properties or pumping rates, and no provision is made to account for stacked ground-water-flow systems that have different flow patterns at different depths. Small-scale flow systems (hundreds to thousands of feet) associated with minor water bodies are not considered; as a result, the model is not currently designed for simulating site-specifi

Scientific Investigations Report↗

Description and evaluation of numerical groundwater flow models for the Edwards Aquifer, south-central Texas

A substantial number of public water system wells in south-central Texas withdraw groundwater from the karstic, highly productive Edwards aquifer. However, the use of numerical groundwater flow models to aid in the delineation of contributing areas for public water system wells in the Edwards aquifer is problematic because of the complex hydrogeologic framework and the presence of conduit-dominated flow paths in the aquifer. The U.S. Geological Survey, in cooperation with the Texas Commission on Environmental Quality, evaluated six published numerical groundwater flow models (all deterministic) that have been developed for the Edwards aquifer San Antonio segment or Barton Springs segment, or both. This report describes the models developed and evaluates each with respect to accessibility and ease of use, range of conditions simulated, accuracy of simulations, agreement with dye-tracer tests, and limitations of the models. These models are (1) GWSIM model of the San Antonio segment, a FORTRAN computer-model code that pre-dates the development of MODFLOW; (2) MODFLOW conduit-flow model of San Antonio and Barton Springs segments; (3) MODFLOW diffuse-flow model of San Antonio and Barton Springs segments; (4) MODFLOW Groundwater Availability Modeling [GAM] model of the Barton Springs segment; (5) MODFLOW recalibrated GAM model of the Barton Springs segment; and (6) MODFLOW–DCM (dual conductivity model) conduit model of the Barton Springs segment. The GWSIM model code is not commercially available, is limited in its application to the San Antonio segment of the Edwards aquifer, and lacks the ability of MODFLOW to easily incorporate newly developed processes and packages to better simulate hydrologic processes. MODFLOW is a widely used and tested code for numerical modeling of groundwater flow, is well documented, and is in the public domain. These attributes make MODFLOW a preferred code with regard to accessibility and ease of use. The MODFLOW conduit-flow model incorporates improvements over previous models by using (1) a user-friendly interface, (2) updated computer codes (MODFLOW–96 and MODFLOW–2000), (3) a finer grid resolution, (4) less-restrictive boundary conditions, (5) an improved discretization of hydraulic conductivity, (6) more accurate estimates of pumping stresses, (7) a long transient simulation period (54 years, 1947–2000), and (8) a refined representation of high-permeability zones or conduits. All of the models except the MODFLOW–DCM conduit model have limitations resulting from the use of Darcy’s law to simulate groundwater flow in a karst aquifer system where non-Darcian, turbulent flow might actually dominate. The MODFLOW–DCM conduit model is an improvement in the ability to simulate karst-like flow conditions in conjunction with porous-media-type matrix flow. However, the MODFLOW–DCM conduit model has had limited application and testing and currently (2008) lacks commercially available pre- and post-processors. The MODFLOW conduit-flow and diffuse-flow Edwards aquifer models are limited by the lack of calibration for the northern part of the Barton Springs segment (Travis County) and their reliance on the use of the calibrated hydraulic conductivity and storativity values from the calibrated Barton Springs segment GAM model. The major limitation of the Barton Springs segment GAM and recalibrated GAM models is that they were calibrated to match measured water levels and springflows for a restrictive range of hydrologic conditions, with each model having different hydraulic conductivity and storativity values appropriate to the hydrologic conditions that were simulated. The need for two different sets of hydraulic conductivity and storativity values increases the uncertainty associated with the accuracy of either set of values, illustrates the non-uniqueness of the model solution, and probably most importantly demonstrates the limitations of using a one-layer model to represent the heterogeneous hydrostratigraphic units composing the Edwards aquifer. In general, the best matches or agreement between groundwater flow directions inferred by numerical model simulation, and by dye-tracer tests, are observed where model outputs accurately reproduce the configuration of the potentiometric surface with regard to the positions of major and minor groundwater troughs and divides. None of the models, with the possible exception of the MODFLOW–DCM conduit model, has a documented capability to accurately simulate travel times for conduit-dominated velocities in the Edwards aquifer. Public water system assessments of wells in the Barton Springs segment of the Edwards aquifer, and elsewhere where conduit-flow conditions are thought to dominate aquifer hydraulic behavior, might be enhanced by use of either the MODFLOW–DCM model or the newly developed U.S. Geological Survey MODFLOW Conduit-Flow Process module for MODFLOW–2005, because each incorporates a type of dual or triple hydraulic conductivity approach and has the capability to explicitly simulate turbulent flow and conduit hydraulic characteristics.

Texas↗

Digital data sets that describe aquifer characteristics of the Elk City Aquifer in western Oklahoma

ARC/INFO export and nonproprietary format files This diskette contains digitized aquifer boundaries and maps of hydraulic conductivity, recharge, and ground-water level elevation contours for the Elk City aquifer in western Oklahoma. The aquifer covers an area of approximately 193,000 acres and supplies ground water for irrigation, domestic, and industrial purposes in Beckham, Custer, Roger Mills, and Washita Counties along the divide between the Washita and Red River basins. The Elk City aquifer consists of the Elk City Sandstone and overlying terrace deposits, made up of clay, silt, sand and gravel, and dune sands in the eastern part and sand and gravel of the Ogallala Formation (or High Plains aquifer) in the western part of the aquifer. The Elk City aquifer is unconfined and composed of very friable sandstone, lightly cemented with clay, calcite, gypsum, or iron oxide. Most of the grains are fine-sized quartz but the grain size ranges from clay to cobble in the aquifer. The Doxey Shale underlies the Elk City aquifer and acts as a confining unit, restricting the downward movement of ground water. All of the data sets were digitized and created from information and maps in a ground-water modeling thesis and report of the Elk City aquifer. The maps digitized were published at a scale of 1:63,360. Ground-water flow models are numerical representations that simplify and aggregate natural systems. Models are not unique; different combinations of aquifer characteristics may produce similar results. Therefore, values of hydraulic conductivity and recharge used in the model and presented in this data set are not precise, but are within a reasonable range when compared to independently collected data.

Oklahoma↗

Application of the Stream Salmonid Simulator (S3) model to assess fall Chinook salmon (Oncorhynchus tshawytscha) production in the American River, California

Executive Summary Anadromous fish returning to the lower American River are restricted to 36 kilometers of free-flowing river between Nimbus Dam and American River’s confluence with the Sacramento River, California. Salmon in the American River provide an important freshwater recreational fishery. However, annual salmon production in the American River in recent years has been low relative to the mid-1990s (Surface Water Resources, Inc., 2001). To investigate the low production of fall-run Chinook salmon ( Oncorhynchus tshawytscha ), the Bureau of Reclamation requested that the U.S. Geological Survey apply the Stream Salmonid Simulator (S3) model to the population of fall-run Chinook salmon on the American River. The American River was chosen among seven candidate Sacramento Basin rivers for S3 application. The American River was selected because of its management and public interest, recently low anadromous fish production, and rich time series of key demographic data needed for S3 application. Data that were not available, however, were empirical estimates on juvenile salmon habitat suitability in the American River. Therefore, a large component of applying S3 to the American River was devoted to the estimation of juvenile salmon habitat suitability and capacity. This entailed snorkeling the lower American River for 3 weeks in March 2021 during the early out-migration period for juvenile Chinook salmon. These efforts were fruitful and showed that the typically small fish (<55 millimeters) in the American River preferred much shallower depths than predicted by habitat suitability criteria derived from the literature for this population. Having empirical estimates on juvenile salmon in the American River provided a solid foundation from which to simulate the population using the S3 model. The S3 model is a spatially explicit population model that runs on a daily time step to simulate redd superimposition, egg maturation, fry emergence and the subsequent growth, survival, and emigration of juvenile Chinook salmon from the river. The key features of this model relevant to this report include (1) a temperature-dependent bioenergetics model driving daily growth rates; (2) density-dependent dynamics that are influenced by the effect of flow on suitable habitat area; and (3) within-year habitat, river flow, and water temperature effects specific to spawning, egg incubation, and fry, parr, and smolt life stages. We used estimates of spawning escapement and geo-referenced redd locations to quantify the spatial and temporal distribution of female spawners for brood years 2014–19. These estimates of female spawners initiate the simulation of each year’s juvenile salmon emergence and emigration over a spatial domain extending from Nimbus Dam to the river’s confluence with the Sacramento River. Using weekly estimates of juvenile salmon abundance and size (fork length) that passed the Watt Avenue fish trap (river kilometer 14.7), we calibrated the S3 model by estimating three key demographic parameters for each year, y : (1) S y , the average daily survival probability, (2) M 0y , the intercept for density-dependence in movement, representing the average daily probability of remaining in a habitat at zero abundance, and (3) C y , the average daily proportion of maximum consumption. These parameters were obtained by minimizing the Mallow’s distance (Lupu and others, 2017) between distributions of weekly abundances and sizes of fish at the traps and weekly simulated abundances and sizes (by S3). Investigation of model fit showed excellent agreement between simulated annual abundances and the abundance of fish passing the fish trap. However, when we compared weekly abundances at the fish trap, S3 under-predicted peaks and over-predicted troughs in the time series of weekly abundances at the fish trap. Thus, some unknown within-year effects have yet to be identified and incorporated in the S3 model. Identifying these important effects and incorporating them in the S3 model would help explain the lack of fit between estimated and simulated weekly abundances. We estimated parameters for 6 years that included a wide range of female spawner abundances (3,057–10,753) and water year types (Critical–Wet). We contrast our estimated parameters to the corresponding number of female spawners and the water year type for the Sacramento Valley. By happenstance, years having higher annual spawner abundances concurred with Critical to Dry water year types. Estimates of survival trended lower with higher spawner abundances and Critical to Dry conditions. In contrast, the extremely wet water year of 2017 had the lowest M 0y , suggesting less density-dependence in fish movement, and the lowest C y , suggesting lower average consumption in this year. When this high-flow year was excluded, a trend towards higher probabilities of fish remaining in a habitat at low abundance and lower proportions of maximum consumption was apparent from Critical to Wet conditions, but only 5 years of data were included. Except for 2017, daily proportions of maximum consumption were relatively high ( C y > 0.83), suggesting that fish were feeding at reasonably high proportions relative to the expected maximum consumption as defined by the “Wisconsin” bioenergetics model (Stewart and Ibarra, 1991). Survival estimates from fry emergence to outmigration at the Sacramento River confluence were generally low when integrated over time. The highest daily survival probability was S y = 0.93 in 2019, or 50 percent total mortality after 10 days. In contrast, our lowest daily survival probability was S y = 0.74 in 2015, or 95 percent total mortality after 10 days. Consequently, even our highest estimated daily survival probability might be considered low. This is especially true given that S y was estimated over a relatively short distance (<14.7 kilometers) from emergence to the Watt Avenue fish trap. Several factors, including our assumed and relatively high daily egg survival rate of 0.9975, could influence juvenile survival estimates. For example, an egg survival rate of 0.9975 results in 3-percent total mortality after 10 days. Egg mortality estimates used in S3 calibration were approximated from egg survivorship studies in the Yakima River, Washington (Johnson and others, 2012), and remains one of the greater uncertainties in S3 when estimating survival across life stages. By including bona fide estimates of egg survival in S3 simulations, the validity of the S3’s current daily egg survival rate could be assessed specifically for the American River. Tagging studies also could provide S3 with direct estimates of juvenile survival and movement; survival during egg incubation then could be estimated indirectly via model fitting.

California↗

U. S. G. S. MODULAR GROUND-WATER FLOW MODEL: DESIGNED TO BE UNDERSTOOD AND ADAPTED.

The paper discusses a carefully designed model program and its complementary complete description of all of the physical and mathematical concepts used in the model. The model program consists of a series of independent subroutines called modules. Modules are grouped by hydrologic function into 'packages. ' A report describing the model was written in three levels. At one level the hydrologic concepts on which the model is based are explained. The second level describes how the concepts were incorporated into the program. The third level describes specific modules in sufficient depth to permit modifications to the model.

Conference Paper↗

Simulation of flow and sediment mobility using a multidimensional flow model for the White Sturgeon critical-habitat reach, Kootenai River near Bonners Ferry, Idaho

In 1994, the Kootenai River white sturgeon (Acipenser transmontanus) was listed as an Endangered Species as a direct result of two related observations. First, biologists observed that the white sturgeon population in the Kootenai River was declining. Second, they observed a decline in recruitment of juvenile sturgeon beginning in the 1950s with an almost total absence of recruitment since 1974, following the closure of Libby Dam in 1972. This second observation was attributed to changes in spawning and (or) rearing habitat resulting from alterations in the physical habitat, including flow regime, sediment-transport regime, and bed morphology of the river. The Kootenai River White Sturgeon Recovery Team was established to find and implement ways to improve spawning and rearing habitat used by white sturgeon. They identified the need to develop and apply a multidimensional flow model to certain reaches of the river to quantify physical habitat in a spatially distributed manner. The U.S. Geological Survey has addressed these needs by developing, calibrating, and validating a multidimensional flow model used to simulate streamflow and sediment mobility in the white sturgeon critical-habitat reach of the Kootenai River. This report describes the model and limitations, presents the results of a few simple simulations, and demonstrates how the model can be used to link physical characteristics of streamflow to biological or other habitat data. This study was conducted in cooperation with the Kootenai Tribe of Idaho along a 23-kilometer reach of the Kootenai River, including the white sturgeon spawning reach near Bonners Ferry, Idaho that is about 108 to 131 kilometers below Libby Dam. U.S. Geological Survey's MultiDimensional Surface-Water Modeling System was used to construct a flow model for the critical-habitat reach of the Kootenai River white sturgeon, between river kilometers 228.4 and 245.9. Given streamflow, bed roughness, and downstream water-surface elevation, the model computes the velocity field, water-surface elevations, and boundary shear stress throughout the modeled reach. The 17.5 kilometer model reach was subdivided into two segments on the basis of predominant grain size: a straight reach with a sand, gravel, and cobble substrate located between the upstream model boundary at river kilometer 245.9 and the upstream end of Ambush Rock at river kilometer 244.6, and a meandering reach with a predominately sand substrate located between upstream end of Ambush Rock and the downstream model boundary at river kilometer 228.4. Model cell size in the x and y (horizontal) dimensions is 5 meters by 5 meters along the computational grid centerline with 15 nodes in the z (vertical) dimension. The model was calibrated to historical streamflows evenly distributed between 141.6 and 2,548.9 cubic meters per second. The model was validated by comparing simulated velocities with velocities measured at 15 cross sections during steady streamflow. These 15 cross sections were each measured multiple (7-13) times to obtain velocities suitable for comparison to the model results. Comparison of modeled and measured velocities suggests that the model does a good job of reproducing flow patterns in the river, although some discrepancies were noted. The model was used to simulate water-surface elevation, depth, velocity, bed shear stress, and sediment mobility for Kootenai River streamflows of 170, 566, 1,130, 1,700, and 2,270 cubic meters per second (6,000, 20,000, 40,000, 60,000, and 80,000 cubic feet per second). The three lowest streamflow simulations represent a range of typical river conditions before and since the construction of Libby Dam, and the highest streamflow simulation (2,270 cubic meters per second) is approximately equal to the annual median peak streamflow prior to emplacement of Libby Dam in 1972. Streamflow greater than 566 cubic meters per second were incrementally increased by 570 cubic meters per second. For each

Scientific Investigations Report↗

LakeVOC; A Deterministic Model to Estimate Volatile Organic Compound Concentrations in Reservoirs and Lakes

This report documents LakeVOC, a model to estimate volatile organic compound (VOC) concentrations in lakes and reservoirs. LakeVOC represents the lake or reservoir as a two-layer system and estimates VOC concentrations in both the epilimnion and hypolimnion. The air-water flux of a VOC is characterized in LakeVOC in terms of the two-film model of air-water exchange. LakeVOC solves the system of coupled differential equations for the VOC concentration in the epilimnion, the VOC concentration in the hypolimnion, the total mass of the VOC in the lake, the volume of the epilimnion, and the volume of the hypolimnion. A series of nine simulations were conducted to verify LakeVOC representation of mixing, dilution, and gas exchange characteristics in a hypothetical lake, and two additional estimates of lake volume and MTBE concentrations were done in an actual reservoir under environmental conditions. These 11 simulations showed that LakeVOC correctly handled mixing, dilution, and gas exchange. The model also adequately estimated VOC concentrations within the epilimnion in an actual reservoir with daily input parameters. As the parameter-input time scale increased (from daily to weekly to monthly, for example), the differences between the measured-averaged concentrations and the model-estimated concentrations generally increased, especially for the hypolimnion. This may be because as the time scale is increased from daily to weekly to monthly, the averaging of model inputs may cause a loss of detail in the model estimates.

Open-File Report↗

MODFLOW–LGR—Documentation of ghost node local grid refinement (LGR2) for multiple areas and the boundary flow and head (BFH2) package

This report documents the addition of ghost node Local Grid Refinement (LGR2) to MODFLOW-2005, the U.S. Geological Survey modular, transient, three-dimensional, finite-difference groundwater flow model. LGR2 provides the capability to simulate groundwater flow using multiple block-shaped higher-resolution local grids (a child model) within a coarser-grid parent model. LGR2 accomplishes this by iteratively coupling separate MODFLOW-2005 models such that heads and fluxes are balanced across the grid-refinement interface boundary. LGR2 can be used in two-and three-dimensional, steady-state and transient simulations and for simulations of confined and unconfined groundwater systems. Traditional one-way coupled telescopic mesh refinement methods can have large, often undetected, inconsistencies in heads and fluxes across the interface between two model grids. The iteratively coupled ghost-node method of LGR2 provides a more rigorous coupling in which the solution accuracy is controlled by convergence criteria defined by the user. In realistic problems, this can result in substantially more accurate solutions and require an increase in computer processing time. The rigorous coupling enables sensitivity analysis, parameter estimation, and uncertainty analysis that reflects conditions in both model grids. This report describes the method used by LGR2, evaluates accuracy and performance for two-and three-dimensional test cases, provides input instructions, and lists selected input and output files for an example problem. It also presents the Boundary Flow and Head (BFH2) Package, which allows the child and parent models to be simulated independently using the boundary conditions obtained through the iterative process of LGR2.

Techniques and Methods↗

Geologic and geophysical models for Osage County, Oklahoma, with implications for groundwater resources

This report summarizes a three-dimensional (3-D) geologic model that was constructed to provide a framework to investigate groundwater resources of the Osage Nation in northeastern Oklahoma. This report also presents an analysis of an airborne electromagnetic (AEM) survey that assessed the spatial variation of electrical resistivity to depths as great as 300 meters in the subsurface. The report and model provide support for a countywide assessment of groundwater resources, emphasizing the Upper Pennsylvanian rock units in the shallow subsurface of central and eastern Osage County having electrical resistivity properties that may indicate aquifers. Surface outcrops and subsurface stratigraphic picks on wire-line geophysical logs of Upper Pennsylvanian&ndash;Lower Permian sedimentary rock were used to construct a 3-D model of the geologic subsurface as an aid for evaluating groundwater resources in Osage County. Quaternary alluvium and terraces along major streams and the Arkansas River are included in the geologic framework model. Data from the AEM survey were subjected to quality-control procedures, truncated at depth of investigation (DOI), and then used to build a 3-D electrical resistivity model making use of secondary and tertiary interpolation profiles between primary data profiles. The AEM data highlight westward-inclined resistivity gradients that parallel the shallow dip of bedrock strata; bodies have resistivity >30 ohm-meters, and extend as much as 10 kilometers (km) down the dip of host geologic units. Volume analysis and internal imaging of an integrated 3-D geology and electrical resistivity model give a proxy for likely aquifer units with large relative volumes of high resistivity: Quaternary alluvium, Elgin Sandstone Lentil in the upper part of the Vamoosa Group, Tallant Formation, and parts of a combined Wann-Iola-Chanute Formation. Less voluminous, high-resistivity bodies correspond to intervals in the lower part of the Vamoosa Group in the east-central part of the county and probable limestone intervals in the upper part of the Vanoss Group in the northwest part of the county. Northwestern and eastern troughs of potable water previously defined for central Osage County generally correspond to down-dip projections of high-resistivity bodies associated with the Elgin Sandstone Lentil of the Vamoosa Group and Tallant Formation, respectively.

Oklahoma↗

Transient calibration of a groundwater-flow model of Chimacum Creek Basin and vicinity, Jefferson County, Washington: a supplement to Scientific Investigations Report 2013-5160

A steady-state groundwater-flow model described in Scientific Investigations Report 2013-5160, ”Numerical Simulation of the Groundwater-Flow System in Chimacum Creek Basin and Vicinity, Jefferson County, Washington” was developed to evaluate potential future impacts of growth and of water-management strategies on water resources in the Chimacum Creek Basin. This supplement to that report describes the unsuccessful attempt to perform a calibration to transient conditions on the model. The modeled area is about 64 square miles on the Olympic Peninsula in northeastern Jefferson County, Washington. The geologic setting for the model area is that of unconsolidated deposits of glacial and interglacial origin typical of the Puget Sound Lowlands. The hydrogeologic units representing aquifers are Upper Aquifer (UA, roughly corresponding to recessional outwash) and Lower Aquifer (LA, roughly corresponding to advance outwash). Recharge from precipitation is the dominant source of water to the aquifer system; discharge is primarily to marine waters below sea level and to Chimacum Creek and its tributaries. The model is comprised of a grid of 245 columns and 313 rows; cells are a uniform 200 feet per side. There are six model layers, each representing one hydrogeologic unit: (1) Upper Confining unit (UC); (2) Upper Aquifer unit (UA); (3) Middle Confining unit (MC); (4) Lower Aquifer unit (LA); (5) Lower Confining unit (LC); and (6) Bedrock unit (OE). The transient simulation period (October 1994–September 2009) was divided into 180 monthly stress periods to represent temporal variations in recharge, discharge, and storage. An attempt to calibrate the model to transient conditions was unsuccessful due to instabilities stemming from oscillations in groundwater discharge to and recharge from streamflow in Chimacum Creek. The model as calibrated to transient conditions has mean residuals and standard errors of 0.06 ft ±0.45 feet for groundwater levels and 0.48 ± 0.06 cubic feet per second for flows. Although the expected seasonal trends were observed in model results, the typical observed annual variation of groundwater levels of about 2 feet was not. Streamflow at the most downstream observation point was about three times larger than simulated streamflow. Because the transient version of the model proved inherently unstable, it was not used to simulate forecast conditions for alternate hydrologic or anthropogenic changes. Adaptation of alternate stream simulation packages, such as RIV, or newer versions of MODFLOW, such as MODFLOW-NWT, could possibly assist with achieving calibration to transient conditions.

Washington↗

Sensitivity of Earthquake Loss Estimates to Source Modeling Assumptions and Uncertainty

Introduction: This report explores how uncertainty in an earthquake source model may affect estimates of earthquake economic loss. Specifically, it focuses on the earthquake source model for the San Francisco Bay region (SFBR) created by the Working Group on California Earthquake Probabilities. The loss calculations are made using HAZUS-MH, a publicly available computer program developed by the Federal Emergency Management Agency (FEMA) for calculating future losses from earthquakes, floods and hurricanes within the United States. The database built into HAZUS-MH includes a detailed building inventory, population data, data on transportation corridors, bridges, utility lifelines, etc. Earthquake hazard in the loss calculations is based upon expected (median value) ground motion maps called ShakeMaps calculated for the scenario earthquake sources defined in WGCEP. The study considers the effect of relaxing certain assumptions in the WG02 model, and explores the effect of hypothetical reductions in epistemic uncertainty in parts of the model. For example, it addresses questions such as what would happen to the calculated loss distribution if the uncertainty in slip rate in the WG02 model were reduced (say, by obtaining additional geologic data)? What would happen if the geometry or amount of aseismic slip (creep) on the region's faults were better known? And what would be the effect on the calculated loss distribution if the time-dependent earthquake probability were better constrained, either by eliminating certain probability models or by better constraining the inherent randomness in earthquake recurrence? The study does not consider the effect of reducing uncertainty in the hazard introduced through models of attenuation and local site characteristics, although these may have a comparable or greater effect than does source-related uncertainty. Nor does it consider sources of uncertainty in the building inventory, building fragility curves, and other assumptions adopted in the loss calculations. This is a sensitivity study aimed at future regional earthquake source modelers, so that they may be informed of the effects on loss introduced by modeling assumptions and epistemic uncertainty in the WG02 earthquake source model.

Open-File Report↗

A latest Pleistocene and Holocene composite tephrostratigraphic framework for northeastern North America

Lakes and bogs in northeastern North America preserve tephra deposits sourced from multiple volcanic systems in the Northern Hemisphere. However, most studies of these deposits focus on specific Holocene intervals and the latest Pleistocene, providing snapshots rather than a full picture. We combine new data with previous work, supplemented by a broad review of the characteristics and ages of potential source regions and volcanoes, to develop the first composite tephrostratigraphic framework covering the last ~14,000 years for this region. We report new cryptotephra records from three ombrotrophic peat bogs—Irwin Smith (Michigan), Bloomingdale (New York), and Sidney Bog (Maine)—as well as new analyses and age models from previously reported sites, Nordan’s Pond Bog (Newfoundland) and Thin-Ice Pond (Nova Scotia). A new tephra (Iliinsky) from the NGRIP and GRIP ice cores is also presented as it can be correlated to new data from these terrestrial records and helps validate radiocarbon age models. We identify 21 new tephra in addition to the 15 already known, several of which cover the entire region – the White River Ash east, Newberry Pumice, Ruppert (NDN230), and Mazama. For the first time we find Mount St. Helens Yn (ca. 3660 cal yr BP) and a set P tephra (~3000–2550 cal yr BP), and confirm the presence of Jala Pumice from Volcan Ceboruco, Mexico, and KS1 from Ksudach volcano, Kamchatka. We describe new “ultra-distal” tephra, including the early Holocene KS 2 eruption, and propose correlations to volcanoes Iliinsky and Shiveluch of Kamchatka, and Ushishir of the Kurile Islands. Not all of these tephra represent large eruptions, with several plausible correlations to sub-Plinian events. Using Bayesian age-modeling, we present new age estimates for the newly described tephra, for tephra with previously poor age control, and for several proximal correlatives. Overall, we demonstrate northeastern North America’s importance for providing transcontinental linkages between paleoenvironmental records and providing insights into ash distribution from different styles and sizes of eruptions.

northeastern North America↗

Simulation of the Groundwater-Flow System in Pierce, Polk, and St. Croix Counties, Wisconsin

Groundwater is the sole source of residential water supply in Pierce, Polk, and St. Croix Counties, Wisconsin. A regional three-dimensional groundwater-flow model and three associated demonstration inset models were developed to simulate the groundwater-flow systems in the three-county area. The models were developed by the U.S. Geological Survey in cooperation with the three county governments. The objectives of the regional model of Pierce, Polk, and St. Croix Counties were to improve understanding of the groundwaterflow system and to develop a tool suitable for evaluating the effects of potential water-management programs. The regional groundwater-flow model described in this report simulates the major hydrogeologic features of the modeled area, including bedrock and surficial aquifers, groundwater/surface-water interactions, and groundwater withdrawals from high-capacity wells. Results from the regional model indicate that about 82 percent of groundwater in the three counties is from recharge within the counties; 15 percent is from surface-water sources, consisting primarily of recirculated groundwater seepage in areas with abrupt surface-water-level changes, such as near waterfalls, dams, and the downgradient side of reservoirs and lakes; and 4 percent is from inflow across the county boundaries. Groundwater flow out of the counties is to streams (85 percent), outflow across county boundaries (14 percent), and pumping wells (1 percent). These results demonstrate that the primary source of groundwater withdrawn by pumping wells is water that recharges within the counties and would otherwise discharge to local streams and lakes. Under current conditions, the St. Croix and Mississippi Rivers are groundwater discharge locations (gaining reaches) and appear to function as 'fully penetrating' hydraulic boundaries such that groundwater does not cross between Wisconsin and Minnesota beneath them. Being hydraulic boundaries, however, they can change in response to water withdrawals. Tributary rivers act as 'partially penetrating' hydraulic boundaries such that groundwater can flow underneath them through the deep sandstone aquifers. The model also demonstrates the effects of development on groundwater in the study area. Water-level declines since predevelopment (no withdrawal wells) are most pronounced where pumping is greatest and flow between layered aquifers is impeded by confining units or faults. The maximum simulated water-level decline is about 40 feet in the deep Mount Simon aquifer below the city of Hudson, Wisconsin. Three inset models were extracted from the regional model to demonstrate the process and additional capabilities of the U.S. Geological Survey MODFLOW code. Although the inset models were designed to provide information about the groundwater-flow system, results from the inset models are presented for demonstration purposes only and are not sufficiently detailed or calibrated to be used for decisionmaking purposes without refinement. Simulation of groundwater/lake-water interaction around Twin Lakes near Roberts, in St. Croix County, Wisconsin, showed that groundwater represents approximately 5 to 20 percent of the overall lake-water budget. Groundwater-contributing areas to streams in western Pierce County are generally similar in size to the surface-water-contributing areas but do not necessarily correspond to the same land area. Transient streamflow simulations of Osceola Creek in Polk County demonstrate how stream base flow can be influenced not only by seasonal precipitation and recharge variability but also by systematic changes to the system, such as groundwater withdrawal from wells.

Scientific Investigations Report↗

Simulation of the regional groundwater-flow system of the Menominee Indian Reservation, Wisconsin

A regional, two-dimensional, steady-state groundwater-flow model was developed to simulate the groundwater-flow system and groundwater/surface-water interactions within the Menominee Indian Reservation. The model was developed by the U.S. Geological Survey (USGS), in cooperation with the Menominee Indian Tribe of Wisconsin, to contribute to the fundamental understanding of the region&rsquo;s hydrogeology. The objectives of the regional model were to improve understanding of the groundwater-flow system, including groundwater/surface-water interactions, and to develop a tool suitable for evaluating the effects of potential regional water-management programs. The computer code GFLOW was used because of the ease with which the model can simulate groundwater/surface-water interactions, provide a framework for simulating regional groundwater-flow systems, and be refined in a stepwise fashion to incorporate new data and simulate groundwater-flow patterns at multiple scales. Simulations made with the regional model reproduce groundwater levels and stream base flows representative of recent conditions (1970&ndash;2013) and illustrate groundwater-flow patterns with maps of (1) the simulated water table and groundwater-flow directions, (2) probabilistic areas contributing recharge to high-capacity pumped wells, and (3) estimation of the extent of infiltrated wastewater from treatment lagoons. The groundwater-flow model described in this report simulates the major hydrogeologic features of the modeled area, including surficial unconsolidated aquifers, groundwater/surface-water interactions, and groundwater withdrawals from existing high-capacity production wells. Areas contributing recharge to pumped high-capacity wells on the Menominee Indian Reservation were delineated by tracking simulated water particles from the water table to wells in combination with Monte Carlo techniques, and maps of the probability of capture for each well nest were produced. Groundwater-agebased areas contributing recharge to wells were simulated by using the calibrated set of parameters and porosity values adjusted to account for bias in simulated saturated thickness. Simulations were performed for current (2013) pumping rates. The simulations show a range in sensitivity of the simulated areas contributing recharge to wells given the parameters evaluated through the Monte Carlo analysis. The areas contributing recharge to supply wells for the villages of Zoar and Neopit are long and narrow, with a sharp gradation from high to low probability of capture. The areas contributing recharge to supply wells for Middle Village and the village of Keshena exhibit a sharp gradation from high to low probability over a relatively small area between the well and a local groundwater mound. The highest probability areas contributing recharge to the supply wells for the Villages of Onekewat and Redwing are in the immediate vicinity of the wells. These wells also have an extensive area with low probability for capturing water that is likely due to a locally low hydraulic gradient and the large degree of uncertainty associated with the lakebed resistance parameters that control interaction between groundwater and local lakes. Additional field investigations and associated local model refinements would facilitate further reductions in uncertainty associated with simulated areas contributing recharge to the wells. The likely extent of the Neopit wastewater plume was simulated by using the groundwater-flow model and Monte Carlo techniques to evaluate the sensitivity of predictive simulations to a range of model parameter values. Wastewater infiltrated from the currently operating lagoons flows predominantly south toward Tourtillotte Creek. Some of the infiltrated wastewater is simulated as having a low probability of flowing beneath Tourtillotte Creek to the nearby West Branch Wolf River. Results for the probable extent of the wastewater plume are considered to be qualitative because the method only considers advective flow and does not account for processes affecting contaminant transport in porous media. Therefore, results for the probable extent of the wastewater plume are sensitive to the number of particles used to represent flow from the lagoon and the resolution of a synthetic grid used for the analysis. Nonetheless, it is expected that the qualitative results may be of use for identifying potential downgradient areas of concern that can then be evaluated using the quantitative &ldquo;area contributing recharge to wells&rdquo; method or traditional contaminant-transport simulations.

Wisconsin↗

Four models used for numerical simulation of a borehole radar antenna

In this report are four different models that represent an antenna used by personnel at the U.S. Geological Survey for crosswell investigations. The four models vary in complexity and concomitantly the accuracy with which they represent the actual antenna. These models are used in numerical simulations of the antenna to determine how it radiates radar waves.

Open-File Report↗

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington↗