Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Computers and Geology”

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

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

At least 1,225 records · Page 68Linked to original sources

Cimarron River alluvial aquifer hydrogeologic framework, water budget, and implications for future water availability in the Pawnee Nation Tribal jurisdictional area, Payne County, Oklahoma, 2016–18

The Cimarron River is a free-flowing river and is a major source of water as it flows across Oklahoma. Increased demand for water resources within the Cimarron River alluvial aquifer in north-central Oklahoma (primarily in Payne County) has led to increases in groundwater withdrawals for agriculture, public, irrigation, industrial, and domestic supply purposes. The Pawnee Nation of Oklahoma (Pawnee Nation) is particularly concerned about the sustainability of the Cimarron River alluvial aquifer and whether the aquifer will continue to be a viable water resource for future generations of Tribal members and residents. To better understand current (2021) water resources and possible future water availability in the Pawnee Nation Tribal jurisdictional area, the U.S. Geological Survey, in cooperation with the Bureau of Indian Affairs and the Pawnee Nation of Oklahoma, compiled available hydrogeologic data and developed conceptual and numerical groundwater-flow models for the Cimarron River alluvial aquifer in Payne County, north-central Oklahoma, including a focus area in the Pawnee Nation Tribal jurisdictional area for the 2016–18 study period. A conceptual water budget was created to establish estimates of groundwater fluxes into and out of the aquifer through hydrologic boundaries and groundwater withdrawals for use in the numerical groundwater-flow model. The conceptual water budget focuses on the alluvial aquifer, meaning that inflows include sources of water to the aquifer and that outflows include sources of water out of the aquifer, such as base-flow contributions to the Cimarron River. The conceptual water budget was constructed by using data from 2017 (the most complete year of record for each data type included in the model) for the Pawnee Nation subdomain of the Cimarron River alluvial aquifer model extent (Pawnee Nation subdomain). Groundwater withdrawals were estimated from groundwater-withdrawal rate information for permanent and temporary permitted wells that was obtained from the Oklahoma Water Resources Board. One-half of each annual permitted groundwater-withdrawal rate allotted was used as the estimated annual groundwater-withdrawal amount. Halving the permitted groundwater-withdrawal rate was done because permitted withdrawal rates are the maximum permitted rate and actual groundwater withdrawals are generally appreciably lower than the maximum permitted rate. Total groundwater withdrawals were estimated as 1,300 acre-feet per year for the Pawnee Nation subdomain. Various hydrogeologic data were measured to assist with model development, including depth to bedrock and water-table altitude data. In support of the model development, analyses pertaining to groundwater flow, groundwater/surface-water interactions, base flows in the Cimarron River, and lithological interpretations in the Pawnee Nation Tribal jurisdictional area were used to compute a conceptual water budget applicable to the 2016–18 study period. A numerical groundwater-flow model was developed using the hydrogeologic framework of the Cimarron River alluvial aquifer and the conceptual water budget. The numerical model consists of a single layer representing alluvium and terrace deposits within the alluvial aquifer model area. Hydraulic conductivities were estimated and modeled for the alluvium and terrace deposits in the alluvial aquifer. Base-flow values were estimated using the base-flow index from streamflow data collected at U.S. Geological Survey streamgages. Stream seepage values were derived from the mean 2017 base-flow index between certain streamgages. Hydraulic conductivities were specified an initial (before calibration) value of 120 feet per day for the alluvium deposits and 16 feet per day for the terrace deposits. The simulated inflows in the numerical groundwater-flow model of the Pawnee Nation subdomain were higher than the inflows of conceptual water budget, and the simulated outflows were lower than the outflows of the conceptual water budget. Overall, simulated base flows matched closely to observed base flows for the 2016 and 2017 stress periods. Simulated streamflow tended to match better with the observed streamflow for 2017, which was the period with the most data for the Cimarron River alluvial aquifer model. Streamflow capture analysis was applied to the steady-state simulation to identify areas of the aquifer where base flows in the Cimarron River were most sensitive to groundwater withdrawals. The initial base-flow value was assigned the value obtained from streamflow-routing software used to simulate stream outflow for the calibrated steady-state base model. Subsequent simulations were run in each active cell in the Pawnee Nation subdomain for a specified groundwater-withdrawal rate of 180,000 cubic feet per day. The study area that includes the Pawnee Nation subdomain is in the upper Arkansas River Basin. A groundwater-withdrawal rate of 180,000 cubic feet per second per day represents a 34 percent increase compared to the highest permitted groundwater-withdrawal rate for the study area, which corresponds to the estimated 34 percent increase in groundwater withdrawals predicted by 2060 for the upper Arkansas River Basin. Simulated streamflow capture was highest in the alluvium deposits adjacent to the Cimarron River; that is, base flow in the Cimarron River decreased the most for simulated groundwater withdrawals in the alluvium deposits adjacent to the Cimarron River. Streamflow capture increased as the distance of a well from the Cimarron River decreased in the simulation. The northeastern part of the Pawnee Nation subdomain showed greater streamflow capture in a broader area; streamflow in that part of the Pawnee Nation subdomain is likely more sensitive to groundwater withdrawals compared to other parts of the Pawnee Nation subdomain.

Oklahoma↗

Analysis of vertical flow during ambient and pumped conditions in four monitoring wells at the Pantex Plant, Carson County, Texas, July-September 2008

The Pantex Plant is a U.S. Department of Energy/National Nuclear Security Administration (USDOE/NNSA)-owned, contractor-operated facility managed by Babcock & Wilcox Technical Services Pantex, LLC (B&W Pantex) in Carson County, Texas, approximately 17 miles northeast of Amarillo. The U.S. Geological Survey, in cooperation with B&W Pantex through the USDOE/NNSA, made a series of flowmeter measurements and collected other borehole geophysical logs during July–September 2008 to analyze vertical flow in screened intervals of four selected monitoring wells (PTX01–1012, PTX06–1044, PTX06–1056, and PTX06–1068) at the Pantex Plant. Hydraulic properties (transmissivity values) of the section of High Plains (Ogallala) aquifer penetrated by the wells also were computed. Geophysical data were collected under ambient and pumped flow conditions in the four monitoring wells. Unusually large drawdowns occurred at two monitoring wells (PTX06–1044 and PTX06–1056) while the wells were pumped at relatively low rates. A decision was made to redevelop those wells, and logs were run again after redevelopment in the two monitoring wells. Logs collected in monitoring well PTX01–1012 during ambient conditions indicate a dynamic environment that probably was affected by pumping of nearby irrigation or public-supply wells. During pumping, downward vertical flow of 0.2 to 2.1 gallons per minute that occurred during ambient conditions was either reversed or reduced. During pumping, a gradual trend of more positive flowmeter values (upward flow) with distance up the well was observed. Estimated total transmissivity for four production zones identified from Flow–B numerical model results taken together was calculated to be about 3,100 feet squared per day. Logs collected in monitoring well PTX06–1044 during ambient conditions before redevelopment indicate a static environment with no flow. During pumping there was upward vertical flow at rates ranging from 0.1 to about 1.5 gallons per minute. During pumping, a gradual trend of more positive flowmeter values (upward flow) with distance up the well was observed. Estimated total transmissivity before redevelopment for five production zones identified from Flow–B numerical model results, and transmissivity values for each zone, are considered to be in error because of the lack of communication between the well and the aquifer before redevelopment. After redevelopment, logs for well PTX06–1044 during ambient conditions indicate a near-static environment with minimal downward flow. During pumping there was upward vertical flow at rates ranging from 0.5 to about 4.8 gallons per minute. During pumping, a gradual trend of more positive flowmeter values with distance up the well was observed. Estimated total transmissivity after redevelopment for the same five identified production zones taken together was calculated to be about 520 feet squared per day. Logs collected in monitoring well PTX06–1056 during ambient conditions before redevelopment indicate a static environment with no flow. During pumping there was upward vertical flow at rates ranging from 0.3 to about 1.5 gallons per minute. During pumping, a gradual trend of more positive flowmeter values (upward flow) with distance up the well was observed. Estimated total transmissivity before redevelopment for four production zones identified from Flow–B numerical model results taken together was calculated to be about 450 feet squared per day. After redevelopment, logs collected in monitoring well PTX06–1056 during ambient conditions indicate a near-static environment with no flow except for a very small amount of downward flow near the bottom of the well. During pumping there was upward vertical flow at rates ranging from 0.7 to about 2.9 gallons per minute. Estimated total transmissivity after redevelopment for five production zones identified from Flow–B numerical model results taken together was calculated to be about 330 feet squared per day. Logs collected in monitoring well PTX06–1068 during ambient conditions indicate a static environment with no flow. During pumping there was upward vertical flow at rates ranging from 0.4 to 4.8 gallons per minute. During pumping, a gradual trend of more positive flowmeter values (upward flow) with distance up the well was observed. Estimated total transmissivity for four production zones identified from Flow–B numerical model results taken together was calculated to be about 200 feet squared per day.

Open-File Report↗

Ungulate migrations of the western United States, Volume 1

Across the western United States, many ungulate herds must migrate seasonally to access resources and avoid harsh winter conditions. Because these migration paths cover vast landscapes (in other words migration distances up to 150 miles [241 kilometers]), they are increasingly threatened by roads, fencing, subdivisions, and other development. Over the last decade, many new tracking studies have been conducted on migratory herds, and analytical methods have been developed that allow for population-level corridors and stopovers to be mapped and prioritized. In 2018, the U.S. Geological Survey assembled a Corridor Mapping Team to provide technical assistance to western states working to map bison, elk, moose, mule deer, and pronghorn migrations using existing Global Positioning System data. Led by the Wyoming Cooperative Fish and Wildlife Research Unit, the team consists of federal scientists, university researchers, and biologists and analysts from participating state agencies. In its first year, the team has worked to develop standardized analytical and computational methods and a workflow applicable to datasets typically collected by state agencies. In 2019, the team completed analyses necessary to map corridors, stopovers, routes and winter ranges in Arizona, Idaho, Nevada, Utah, and Wyoming. A total of 26 corridors, 16 migration routes, 25 stopovers, and 9 winter ranges were mapped across these states and are included in this report. This report and associated data release provide the means for the habitats required for migration to be taken into account by state and federal transportation officials, land and wildlife managers, planners, and other conservationists working to maintain big-game migration in the western states.

Scientific Investigations Report↗

Instrumentation of urban hydrology monitoring sites in southeast Florida

An instrumentation system developed and built in laboratories of the U.S. Geological Survey, Reston, Va., has been used since 1974 to collect synchronized rainfall, runoff, and water-quality data from urban basins. A number of field modifications were made to adapt the system to local hydrologic conditions and for collection of data on about 350 rainstorms at four sites in south Florida. The instrumentation system measures rainfall at three sites, records two pressure or water-level readings for flow computations, and collects and refrigerates up to 24 water samples and rainfall-quality samples. Rainfall and runoff data are recorded every 36 seconds on a six-channel analog recorder so that all variables are time-synchronous. Flow in the storm sewer is computed from pressure measurements in a U-shaped venturi-type constriction or, when the constriction is not used, from water levels. (Kosco-USGS)

Florida↗

Research to More Effectively Manage Critical Ground-Water Basins

As the regional management agency for two of the most heavily used ground-water basins in California, the Water Replenishment District of Southern California (WRD) plays a vital role in sheparding the water resources of southern Los Angeles County. WRD is using the results of the U.S. Geological Survey (USGS) studies to help more effectively manage the Central and West Coast basins in the most efficient, cost-effective way. In partnership with WRD, the USGS is using the latest research tools to study the geohydrology and geochemistry of the two basins. USGS scientists are: *Drilling and collecting detailed data from over 40 multiple-well monitoring sites, *Conducting regional geohydrologic and geochemical analyses, *Developing and applying a computer simulation model of regional ground-water flow. USGS science is providing a more detailed understanding of ground-water flow and quality. This research has enabled WRD to more effectively manage the basins. It has helped the District improve the efficiency of its spreading ponds and barrier injection wells, which replenish the aquifers and control seawater intrusion into the ground-water system.

Fact Sheet↗

Test and calibration of the Digital World-Wide Standardized Seismograph

During the past decade there has been steady progress in the modernization of the global seismograph network operated by the U.S. Geological Survey (USGS). The World-Wide Standardized Seismograph Network (WWSSN) has been augmented by new stations with advanced instrumentation, including the Seismic Research Observatories (SRO) and the modified High-Gain Long-Period (ASRO) stations. One goal in the modernization effort has been to improve signal resolution in the long-period band. A second goal has been to generate a global digital data base to support contemporary computer-based analysis and research. In 1976, a Panel on Seismograph Networks was established by the Committee on Seismology of the National Academy of Sciences to review progress in network seismology and recommend actions that would lead to an improved global data base for seismology. One recommendation in the Panel report (Engdahl, 1977) called for upgrading selected WWSSN stations by the installation of digital recorders. This was viewed as an economical way of expanding the digital network, which had proven itself to be a very promising new tool for earthquake and explosion research. Funds for the development and assembly of 15 digital recorders were provided to the USGS by the Defense Advanced Research Projects Agency and an ad Inoc panel of scientists was convened by the Committee on Seismology to advise the USGS on the selection of stations to be upgraded and on data recording requirements, A total of 19 digital World-Wide Standardized Seismograph (DWWSS) systems will be operational when all are installed. The additional systems were made available through purchase by the USGS and other organizations; for example, the University of Bergen purchased and installed a DWWSS-type recorder and agreed to furnish the USGS with the data. A list of operational and planned DWWSS network stations is given in Table 1.1. As one might expect, the digital recorder turned out to be somewhat more sophisticated than the original concept. It was decided to record three components of long-period data continuously, three components of intermediateperiod data in an event mode, and the vertical-component short-period data in and event mode (with the capability of adding short-period horizontal channels in the future). Special amplifiers were developed for use with the WWSS seismometers, and a 16-bit fixed-point analog-to-digital converter was chosen to provide increased resolution (as opposed to a 16-bit gain-ranged encoder). The microprocessor-based digital recording systems were developed and assembled at the USGS Albuquerque Seismological Laboratory (ASL) and ASL-based technicians began installation at WWSSN stations in 1980. The current and proposed locations of the DWWSSN stations, together with other stations in the Global Digital Seismograph Network (GDSN), are shown on the map in Figure 1.1. A system was operated at Albuquerque for about 18 months, serving as a test bed for evaluation studies. Although the network hardware has been available for some time, the installation of the DWWSSN has proceeded slowly. The National Science Foundation supported installation of six stations and the USGS is funding installation of most of the others; however, the network completion date is conjectural because of funding uncertainties. The DWWSSN stations are supported with supplies and technical assistance from ASL (subject to availability of funds). Data recorded on magnetic tapes are mailed to ASL where they are reviewed for quality, then merged with other GDSN station data on the network-day tapes. Hoffman (1980) provides a description of the network-day tape format. Zirbes and Buland (1981) have developed and published user software for reading and interpreting the day tapes. This report will serve several purposes. One is to provide nominal system transfer functions and calibration information that are needed in the analysis of DWWSSN data. A second purpose is to report on an evaluation of operating characteristics (calibration stability, noise levels, and linearity) that may limit the usefulness of the data and to determine if modifications may be needed to improve the data. It is not an exhaustive study in this respect. We continue to depend mostly on data user feedback to point out deficiencies and we solicit comments whenever anomalies are observed in the data.

Open-File Report↗

Vertical profiles of streambed hydraulic conductivity determined using slug tests in central and western Nebraska

Many issues of water-resources management rely on modeling of ground-water/surfacewater interactions, and streambed hydraulic conductivity is a key parameter controlling the water fluxes across the stream/aquifer interface. However, in central and western Nebraska, this parameter is generally undefined. The U.S. Geological Survey, in cooperation with the Nebraska Platte River Cooperative Hydrology Study Group, performed slug tests at 15 stream sites in the Platte, Republican, and Little Blue River watersheds to determine the hydraulic conductivity of streambeds in central and western Nebraska. Slug tests were completed at several discrete depth intervals using pneumatic or mechanical methods, and the water-level response was monitored on site using a pressure transducer and laptop computer. Responses were analyzed using either the Bouwer and Rice or Springer and Gelhar methods. Vertical profiles of hydraulic conductivity with depth were developed and were compared to available information on lithology. The profiles and corresponding lithology showed that different types of streambeds were tested and suggested that some streambeds display a large variability in hydraulic conductivity with depth. In some cases, hydraulic conductivity values associated with nonstreambed materials could be identified from nearby lithologic descriptions. Seven of 15 sites had streambed values that ranged over more than 3 orders of magnitude, and that variability increased significantly when the measurements considered to be from nonstreambed materials were included. Streambed profiles from the Platte and South Platte River sites generally were more homogeneous and of larger hydraulic conductivity than the other sites. No restrictive layers were detected at any of the streambed sites on the main stems or the flood plains of the main stems of their respective watersheds. Alternatively, the profiles characterized by a restrictive streambed layer at some depth below the streambed surface were all from tributary sites out of the main-stem flood plain. These profiles can be used to represent the streambed hydraulic conductivity in central and western Nebraska in various applications, including modeling ground-water/surfacewater interactions.

Nebraska↗

Littoral transport rates in the Santa Barbara Littoral Cell: a process-based model analysis

Identification of the sediment transport patterns and pathways is essential for sustainable coastal zone management of the heavily modified coastline of Santa Barbara and Ventura County (California, USA). A process-based model application, based on Delft3D Online Morphology, is used to investigate the littoral transport potential along the Santa Barbara Littoral Cell (between Point Conception and Mugu Canyon). An advanced optimalization procedure is applied to enable annual sediment transport computations by reducing the ocean wave climate in 10 wave height - direction classes. Modeled littoral transport rates compare well with observed dredging volumes, and erosion or sedimentation hotspots coincide with the modeled divergence and convergence of the transport gradients. Sediment transport rates are strongly dependent on the alongshore variation in wave height due to wave sheltering, diffraction and focusing by the Northern Channel Islands, and the local orientation of the geologically-controlled coastline. Local transport gradients exceed the net eastward littoral transport, and are considered a primary driver for hot-spot erosion.

California↗

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↗

Rayleigh wave ellipticity measurement uncertainty across the IRIS/USGS and New China Digital Seismograph Networks

Long-period Rayleigh wave horizontal to vertical amplitude (H/V) ratios at a station provide information about local earth structure that is complementary to phase velocity. However, a number of studies have observed that significant scatter appears in these measurements making it difficult to use H/V ratio measurements to resolve earth structure. Some of the scatter in these measurements has been attributed to local geological structure while some has remained unaccounted for. Most Global Seismographic Network (GSN) stations contain two nearby high-quality broad-band seismometers (e.g. in the same vault, but on different piers or in different boreholes). For each broad-band sensor in the IRIS/USGS component of the GSN, we estimate H/V ratios of fundamental mode Rayleigh waves using M > 6.5 earthquakes from 2001 to 2018 (around 19 000 measurements). We compute these ratios at a number of discrete periods (25, 50, 75, 100 and 150 s) and find that for well-isolated Rayleigh waves (windows where the correlation coefficients between radial and the phase-shifted vertical components are greater than 0.9) significant scatter in H/V ratios occurs between colocated sensors (greater than 25 per cent at 100 s period). This suggests the scatter in H/V ratio measurements can be at least partially attributed to extremely local phenomena such as sensor emplacement in the vault. We also find that H/V ratios can vary as a function of event backazimuth, indicating that care must be taken when computing average ratios for a station, as a large number of events from a given region could bias H/V ratio measurements at a station.

Geophysical Journal International↗

Elastic-wave propagation and site amplification in the Salt Lake Valley, Utah, from simulated normal faulting earthquakes

The two-dimensional seismic response of the Salt Lake valley to near- and far-field earthquakes has been investigated from simulations of vertically incident plane waves and from normal-faulting earthquakes generated on the basin-bounding Wasatch fault. The response to normal faulting earthquakes was simulated using a two-dimensional finite-element method and the plane-wave response was calculated from two-dimensional finite-difference simulations. The plane-wave simulations were then compared with observed site amplifications in the Salt Lake valley, based on seismic recordings from nuclear explosions in southern Nevada, that show 10 times greater amplification within the basin than measured values on hard-rock sites. While previous studies attribute this increased site amplification to the near-surface unconsolidated/consolidated alluvial fill contact, our synthetic seismograms suggest that in the frequency band 0.3 to 1.5 Hz at least one-half the site amplification can be attributed to the impedance contrast between the basin sediments and higher velocity basement rocks. Synthetic seismograms from vertically incident plane-wave sources and buried double-couple sources predict large amplitude Rayleigh-wave propagation from the edges of the basin and, in general, uniform site amplification. In contrast, near-field simulations of basin-bounding, normal-faulting earthquakes predict large-amplitude Rayleigh waves propagating westward from the fault across the basin. Spectra of synthetic accelerograms computed from the normal-faulting earthquakes shows that spectral amplification within the basin is primarily due to source directivity with a maxima near the surface projection of the fault that decays rapidly away from the fault. Importantly, the synthetic modeling of near-field earthquake sources show that near-field directivity effects are important and should be considered in an earthquake hazard assessment of the Salt Lake valley and similar geologic settings along the Wasatch Front.

Utah↗

Use of remote-sensing techniques to survey the physical habitat of large rivers

Remote-sensing techniques that can be used to quantitatively characterize the physical habitat in large rivers in the United States where traditional survey approaches typically used in small- and medium-sized streams and rivers would be ineffective or impossible to apply. The state-of-the-art remote-sensing technologies that we discuss here include side-scan sonar, RoxAnn, acoustic Doppler current profiler, remotely operated vehicles and camera systems, global positioning systems, and laser level survey systems. The use of these technologies will permit the collection of information needed to create computer visualizations and hard copy maps and generate quantitative databases that can be used in real-time mode in the field to characterize the physical habitat at a study location of interest and to guide the distribution of sampling effort needed to address other habitat-related study objectives. This report augments habitat sampling and characterization guidance provided by Meador et al. (1993) and is intended for use primarily by U.S. Geological Survey National Water Quality Assessment program managers and scientists who are documenting water quality in streams and rivers of the United States.

Report↗

Rationale and operational plan to upgrade the U.S. gravity database

A concerted effort is underway to prepare a substantially upgraded digital gravity anomaly database for the United States and to make this data set and associated usage tools available on the internet. This joint effort, spearheaded by the geophysics groups at the National Imagery and Mapping Agency (NIMA), University of Texas at El Paso (UTEP), U.S. Geological Survey (USGS), and National Oceanic and Atmospheric Administration (NOAA), is an outgrowth of the new geoscientific community initiative called Geoinformatics (www.geoinformaticsnetwork.org). This dominantly geospatial initiative reflects the realization by Earth scientists that existing information systems and techniques are inadequate to address the many complex scientific and societal issues. Currently, inadequate standardization and chaotic distribution of geoscience data, inadequate accompanying documentation, and the lack of easy-to-use access tools and computer codes for analysis are major obstacles for scientists, government agencies, and educators. An example of the type of activities envisioned, within the context of Geoinformatics, is the construction, maintenance, and growth of a public domain gravity database and development of the software tools needed to access, implement, and expand it. This product is far more than a high quality database; it is a complete data system for a specific type of geophysical measurement that includes, for example, tools to manipulate the data and tutorials to understand and properly utilize the data. On August 9, 2002, twenty-one scientists from the federal, private and academic sectors met at a workshop to discuss the rationale for upgrading both the United States and North American gravity databases (including offshore regions) and, more importantly, to begin developing an operational plan to effectively create a new gravity data system. We encourage anyone interested in contributing data or participating in this effort to contact G.R. Keller or T.G. Hildenbrand. This workshop was the first step in building a web-based data system for sharing quality gravity data and methodology, and it builds on existing collaborative efforts. This compilation effort will result in significant additions to and major refinement of the U.S. database that is currently released publicly by NOAA’s National Geophysical Data Center and will also include an additional objective to substantially upgrade the North American database, released over 15 years ago (Committee for the Gravity Anomaly Map of North America, 1987).

Open-File Report↗

Preliminary slope-stability analysis of Augustine Volcano: Chapter 14 in The 2006 eruption of Augustine Volcano, Alaska

Augustine Volcano has been a prolific producer of large debris avalanches during the Holocene. Originating as landslides from the steep upper edifice, these avalanches typically slide into the surrounding ocean. At least one debris avalanche that occurred in 1883 during an eruption initiated a far-traveled tsunami. The possible occurrence of another edifice collapse and ensuing tsunami was a concern during the 2006 eruption of Augustine. To aid in hazard assessments, we have evaluated the slope stability of Augustine's edifice, using a quasi-three-dimensional, geotechnically based slope-stability model implemented in the computer program SCOOPS. We analyzed the effects of topography, variations in rock strength, and earthquake-induced strong ground motion on the relative stability of millions of potential large (>0.1 km 3 volume) slope failures throughout the edifice. Preliminary results from pre-2006 topography provide three insights. First, the predicted stability of all parts of the upper edifice is approximately the same, suggesting an equal likelihood of slope failure, in agreement with geologic observations that debris avalanches have swept all sectors of the volcano. Second, the least stable (by a small amount) sector is on the east flank where a debris avalanche would flow into deeper ocean water and a resulting tsunami would be directed toward the southwestern part of the Kenai Peninsula. Third, most model scenarios predict stable edifice slopes, and only scenarios assuming extensive weak rocks and moderate to strong ground shaking predict potential large collapses. Because other transient triggering mechanisms, such as shallow magma intrusion, may be needed to instigate slope instability, monitoring ground deformation and seismicity could

Alaska↗

Surface-geophysical investigation of a formerly used defense site, Machiasport, Maine, February 2003

Surface-geophysical surveys were conducted in February 2003 at a formerly used defense site in Maine, where residual chlorinated solvents are affecting off-site domestic water-supply wells. The U.S. Geological Survey and Argonne National Laboratory used surface-geophysical methods, including ground-penetrating radar and seismic-refraction tomography, to characterize the lithology and structure of the bedrock at the site and to identify highly fractured areas that may provide pathways for ground-water flow and contaminant transport. Multifrequency electromagnetic and inductive terrain-conductivity methods also were evaluated, but these techniques were adversely affected by a nearby naval computer and telecommunications station. Interpretation of the data from ground-penetrating radar indicates that depth to the weathered bedrock surface is approximately 0.5 to 3 meters. Reflections from within the bedrock are visible throughout all ground-penetrating radar profiles, and zones of scattered electromagnetic energy may correlate to zones of highly fractured bedrock. Interpretation of the data from seismic-refraction tomography inversion indicates that zones of relatively low seismic velocity and topographic lows may correlate with fractured and water-producing intervals within the bedrock. Integrated interpretation of the results from ground-penetrating radar and seismic-refraction tomography was used to locate boreholes along the surface-geophysical profiles. An integrated analysis of information obtained from the surface- and borehole-geophysical surveys and test drilling will be used by the U.S. Army Corps of Engineers to develop a conceptual model of ground-water flow and solute transport at the site.

Maine↗

Flood-inundation maps for the Amite and Comite Rivers from State Highway 64 to U.S. Highway 190 at Central, Louisiana

Flood-inundation maps for a 14.5-mile reach of the Amite River and a 20.2-mile reach of the Comite River from State Highway 64 to U.S. Highway 190 were created by the U.S. Geological Survey (USGS) in cooperation with the City of Central, Louisiana. These maps, which can be accessed through an interactive mapper at the USGS Flood Inundation Mapping Program website and from a companion USGS data release, depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at the USGS streamgages Amite River at Magnolia, La., (07377300) and Comite River near Comite, La. (07378000). Flood profiles were computed for the Amite and Comite River reaches by using the two-dimensional (2D), finite-volume numerical modeling options in the U.S. Army Corps of Engineers Hydrologic Engineering Center’s River Analysis System (USACE HEC-RAS) software version 5.0.3. Models were calibrated to the current (2018) stage-discharge relations at the Amite River at Magnolia, La., and Comite River near Comite, La., streamgages, water-surface profiles from the March and August 2016 floods, and documented high-water marks from the flood of August 2016. The hydraulic models were used to compute 37 individual water-surface profiles (21 for the Amite River and 16 for the Comite River) at 1.0-foot intervals ranging from the National Weather Service flood stage to the highest peak on record at the two streamgages. The 37 simulated water-surface profiles were used with a light detection and ranging-derived digital elevation model to delineate the flood extent and associated depth at each water level. The delineated areas (inundation maps) were merged into 127 combinations or possible flooding scenarios based on annual peak stage information from the two streamgaging stations. The availability of these maps, along with real-time data delivered via the internet, will provide emergency management personnel and residents with information that is critical for flood-response activities such as evacuations and road closures, as well as for recovery efforts after floods.

Louisiana↗

Hurricane Mitch: Peak Discharge for Selected River Reachesin Honduras

Hurricane Mitch began as a tropical depression in the Caribbean Sea on 22 October 1998. By 26 October, Mitch had strengthened to a Category 5 storm as defined by the Saffir-Simpson Hurricane Scale (National Climate Data Center, 1999a), and on 27 October was threatening the northern coast of Honduras (fig. 1). After making landfall 2 days later (29 October), the storm drifted south and west across Honduras, wreaking destruction throughout the country before reaching the Guatemalan border on 31 October. According to the National Climate Data Center of the National Oceanic and Atmospheric Administration (National Climate Data Center, 1999b), Hurricane Mitch ranks among the five strongest storms on record in the Atlantic Basin in terms of its sustained winds, barometric pressure, and duration. Hurricane Mitch also was one of the worst Atlantic storms in terms of loss of life and property. The regionwide death toll was estimated to be more than 9,000; thousands of people were reported missing. Economic losses in the region were more than $7.5 billion (U.S. Agency for International Development, 1999). Honduras suffered the most widespread devastation during the storm. More than 5,000 deaths, and economic losses of more than $4 billion, were reported by the Government of Honduras. Honduran officials estimated that Hurricane Mitch destroyed 50 years of economic development. In addition to the human and economic losses, intense flooding and landslides scarred the Honduran landscape - hydrologic and geomorphologic processes throughout the country likely will be affected for many years. As part of the U.S. Government's response to the disaster, the U.S. Geological Survey (USGS) conducted post-flood measurements of peak discharge at 16 river sites throughout Honduras (fig. 2). Such measurements, termed 'indirect' measurements, are used to determine peak flows when direct measurements (using current meters or dye studies, for example) cannot be made. Indirect measurements of peak discharge are based on post-flood surveys of the river channel (observed high-water marks, cross sections, and hydraulic properties) and model computation of peak discharge. Determination of the flood peaks associated with Hurricane Mitch will help scientists understand the magnitude of this devastating hurricane. Peak-discharge information also is critical for the proper design of hydraulic structures (such as bridges and levees), delineation of theoretical flood boundaries, and development of stage-discharge relations at streamflow-monitoring sites.

Water-Resources Investigations Report↗

Documentation of a graphical display program for the saturated- unsaturated transport (SUTRA) finite-element simulation model

This report documents a graphical display program for the U. S. Geological Survey finite-element groundwater flow and solute transport model. Graphic features of the program, SUTRA-PLOT (SUTRA-PLOT = saturated/unsaturated transport), include: (1) plots of the finite-element mesh, (2) velocity vector plots, (3) contour plots of pressure, solute concentration, temperature, or saturation, and (4) a finite-element interpolator for gridding data prior to contouring. SUTRA-PLOT is written in FORTRAN 77 on a PRIME 750 computer system, and requires Version 9.0 or higher of the DISSPLA graphics library. The program requires two input files: the SUTRA input data list and the SUTRA simulation output listing. The program is menu driven and specifications for individual types of plots are entered and may be edited interactively. Installation instruction, a source code listing, and a description of the computer code are given. Six examples of plotting applications are used to demonstrate various features of the plotting program. (Author 's abstract)

Water-Resources Investigations Report↗