SEARCH · Geology Reports

Results for “Algorithms”

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

SourceResearch connectionFrom yearThrough year

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

USGS2007

Hawaiian Volcano Observatory 1985 Annual Administrative Report

INTRODUCTORY NOTE The Hawaiian Volcano Observatory Summaries have been published in the current format since 1956. The Quarterly Summaries (1956 through 1973) and the Annual Summaries (1974 through 1985) were originally published as Administrative Reports. These reports have been compiled and published as U.S. Geological Survey Open-File Reports. The quarterly reports have been combined and published as one annual summary. All the summaries from 1956 to the present are now available as .pdf files at http://www.usgs.gov/pubprod. The earthquake summary data are presented as a listing of origin time, depth, magnitude, and other location parameters. Network instrumentation, field station sites, and location algorithms are described. Tilt and other deformation data are included until Summary 77, January to December 1977. From 1978, the seismic and deformation data are published separately, due to differing schedules of data reduction. There are eight quarters - from the fourth quarter of 1959 to the third quarter of 1961 - that were never published. Two of these (4th quarter 1959, 1st quarter 1960) have now been published, using handwritten notes of Jerry Eaton (HVO seismologist at the time) and his colleagues. The seismic records for the remaining six summaries went back to California in 1961 with Jerry Eaton. Other responsibilities intervened, and the seismic summaries were never prepared.

Open-File Report↗

USGS2007

Geosoft eXecutables (GX's) Developed by the U.S. Geological Survey, Version 2.0, with Notes on GX Development from Fortran Code

Introduction Geosoft executables (GX's) are custom software modules for use with the Geosoft Oasis montaj geophysical data processing system, which currently runs under the Microsoft Windows 2000 or XP operating systems. The U.S. Geological Survey (USGS) uses Oasis montaj primarily for the processing and display of airborne geophysical data. The ability to add custom software modules to the Oasis montaj system is a feature employed by the USGS in order to take advantage of the large number of geophysical algorithms developed by the USGS during the past half century. This main part of this report, along with Appendix 1, describes Version 2.0 GX's developed by the USGS or specifically for the USGS by contractors. These GX's perform both basic and advanced operations. Version 1.0 GX's developed by the USGS were described by Phillips and others (2003), and are included in Version 2.0. Appendix 1 contains the help files for the individual GX's. Appendix 2 describes the new method that was used to create the compiled GX files, starting from legacy Fortran source code. Although the new method shares many steps with the approach presented in the Geosoft GX Developer manual, it differs from that approach in that it uses free, open-source Fortran and C compilers and avoids all Fortran-to-C conversion.

Open-File Report↗

USGS2008

Shear wave structure of Umbria and Marche, Italy, strong motion seismometer sites Affected by the 1997-98 Umbria-Marche, Italy, earthquake sequence

A long sequence of earthquakes, eight with magnitudes between 5 and 6, struck the Umbria and Marche regions of central Italy between September 26, 1997 and July 1998. The earthquake swarm caused severe structural damage, particularly to masonry buildings, and resulted in the loss of twelve lives and about 150 injuries. The source of the events was a single seismogenic structure that consists of several faults with a prevailing northwest-southeast strike and crosses the Umbria-Marche border. The focal mechanism of the largest shocks indicates that the events were the product of shallow extensional normal faulting along a NE-SW extension perpendicular to the trend of the Apennines. The network of analog seismometer stations in the Umbria and Marche regions recorded motions of the main September and October 1997 events and a dense array of mobile digital stations, installed since September 29, recorded most of the swarm. The permanent national network Rete Accelerometrica Nazionale (RAN) is administered and maintained by Dipartimento delle Protezione Civile (DPC: Civil Protection Department); the temporary array was managed by Servizio Sismico Nazionale (SSN) in cooperation with small agencies and Universities. ENEA, the operator of many seismometer stations in Umbria, is the public Italian National Agency for New Technologies, Energy and the Environment. Many of the temporary and permanent stations in the Italian seismic network have little or no characterization of seismic velocities. In this study, we investigated 17 Italian sites using an active-source approach that employs low frequency harmonic waves to measure the dispersive nature of surface waves in the ground. We used the Spectral Analysis of Surface Wave (SASW) approach, coupled with an array of harmonic-wave electro-mechanical sources that are driven in-phase to excite the ground. An inversion algorithm using a non-linear least-squares best-fit method is used to compute shear wave velocities for up to 100 meters of the soil column. A draft report was published in the summer of 2008, followed by a comment period, lengthy discussions with Italian colleagues, and improved knowledge of the subsurface at the sites from soil logs. Four of the sites were reprocessed in order to correct issues with phase unwrapping of the field dispersion curves that complicated the velocity profile calculations at the lowest velocity sites. This report presents the final results from the reprocessing effort.

Open-File Report↗

USGS2009

Modeling habitat of the desert tortoise (Gopherus agassizii) in the Mojave and parts of the Sonoran Deserts of California, Nevada, Utah, and Arizona

Habitat modeling is an important tool used to simulate the potential distribution of a species for a variety of basic and applied questions. The desert tortoise ( Gopherus agassizii ) is a federally listed threatened species in the Mojave Desert and parts of the Sonoran Desert of California, Nevada, Utah, and Arizona. Land managers in this region require reliable information about the potential distribution of desert tortoise habitat to plan conservation efforts, guide monitoring activities, monitor changes in the amount and quality of habitat available, minimize and mitigate disturbances, and ultimately to assess the status of the tortoise and its habitat toward recovery of the species. By applying information from the literature and our knowledge or assumptions of environmental variables that could potentially explain variability in the quality of desert tortoise habitat, we developed a quantitative habitat model for the desert tortoise using an extensive set of field-collected presence data. Sixteen environmental data layers were converted into a grid covering the study area and merged with the desert tortoise presence data that we gathered for input into the Maxent habitat-modeling algorithm. This model provides output of the statistical probability of habitat potential that can be used to map potential areas of desert tortoise habitat. This type of analysis, while robust in its predictions of habitat, does not account for anthropogenic changes that may have altered habitat with relatively high potential into areas with lower potential.

California, Nevada, Utah, Arizona↗

USGS2009

Paper birch decline in the Niobrara River Valley, Nebraska: Weather, microclimate, and birch stand conditions

The Niobrara River Valley in north-central Nebraska supports scattered stands of paper birch (Betula papyrifera Marsh), a species more typical of boreal forests. These birch stands are considered to be relictual populations that have persisted since the end of the Wisconsin glaciation, when regional flora was more boreal in nature (Wright 1970, Kaul and others, 1988). Dieback of canopy-sized birch has been observed throughout the Niobrara Valley in recent years, although no onset dates are documented. The current dieback event probably started around or after the early 1980’s. The study objectives were to understand microclimatic conditions in birch stands relative to nearby weather stations and historic weather conditions, and to assess current health conditions of individual birch trees. Temperature was measured every half-hour from June 2005 through October 2007 in 12 birch stands and individual birch tree health was measured as expressed by percent living canopy in these and 13 additional stands in spring 2006 and 2007. Birch site microclimate was compared to data from a National Weather Service station in Valentine, Nebraska, and to an automated weather station at The Nature Conservancy Niobrara Valley Preserve 24 kilometers north of Johnstown, Nebraska. Historic weather data from the Valentine station and another National Weather Service Station at Ainsworth, Nebraska, were used to reconstruct minimum and maximum temperature at The Nature Conservancy and one microclimate monitoring station using Kalman filtering and smoothing algorithms. Birch stand microclimate differed from local weather stations as well as among stands. Birch health was associated with annual minimum temperature regimes; those stands whose annual daily minimum temperature regimes were most like The Nature Conservancy station contained smaller proportions of living trees. Frequency of freeze/thaw conditions capable of inducing rootlet injury and subsequent crown dieback significantly have increased in the second one-half of the period of record (1978–2007) as compared to the first one-half (1948–1977). River location was associated with birch health; upper river sites had significantly healthier trees than north bank sites. Localized microclimates in the birch stands have likely facilitated the persistence of the birch populations in a region otherwise unsuitable for the species. These microclimate differences may reduce frequency of thaw/freeze conditions that can induce root injury and potential crown dieback. A large population decline in the context of increased frequency of potentially injurious climatic events would make population recovery much more difficult now than from 1948 to 1977, when thaw/freeze conditions were less frequent. These conditions, combined with little evidence of recruitment of young birch and great geographic distances from potential immigrant sources, make the future persistence of birch in the Niobrara River Valley stands uncertain.

Open-File Report↗

USGS2010

Simulation of Runoff and Reservoir Inflow for Use in a Flood-Analysis Model for the Delaware River, Pennsylvania, New Jersey, and New York, 2004-2006

A model was developed to simulate inflow to reservoirs and watershed runoff to streams during three high-flow events between September 2004 and June 2006 for the main-stem subbasin of the Delaware River draining to Trenton, N.J. The model software is a modified version of the U.S. Geological Survey (USGS) Precipitation-Runoff Modeling System (PRMS), a modular, physically based, distributed-parameter modeling system developed to evaluate the impacts of various combinations of precipitation, climate, and land use on surface-water runoff and general basin hydrology. The PRMS model simulates time periods associated with main-stem flooding that occurred in September 2004, April 2005, and June 2006 and uses both daily and hourly time steps. Output from the PRMS model was formatted for use as inflows to a separately documented reservoir and riverrouting model, the HEC-ResSim model, developed by the U.S. Army Corps of Engineers Hydrologic Engineering Center to evaluate flooding. The models were integrated through a graphical user interface. The study area is the 6,780 square-mile watershed of the Delaware River in the states of Pennsylvania, New Jersey, and New York that drains to Trenton, N.J. A geospatial database was created for use with a geographic information system to assist model discretization, determine land-surface characterization, and estimate model parameters. The USGS National Elevation Dataset at 100-meter resolution, a Digital Elevation Model (DEM), was used for model discretization into streams and hydrologic response units. In addition, geospatial processing was used to estimate initial model parameters from the DEM and other data layers, including land use. The model discretization represents the study area using 869 hydrologic response units and 452 stream segments. The model climate data for point stations were obtained from multiple sources. These sources included daily data for 22 National Weather Service (NWS) Cooperative Climate Station network stations, hourly data for 15 stations from the National Climatic Data Center, hourly data for 1 station from the NWS Middle Atlantic River Forecast Center records, and daily and hourly data for 7 stations operated by the New York City Department of Environmental Protection. The NWS Multisensor Precipitation Estimate data set for 2001-2007 was used for computing daily precipitation for the model and for computing hourly precipitation for storm simulation periods. Calibration of the PRMS model included regression and optimization algorithms, as well as manual adjustments of model parameters. The general goal of the calibration procedure was to minimize the difference between discharge measured at USGS streamgages and the corresponding discharge simulated by the model. Daily streamflow data from 35 USGS streamgages were used in model calibration. The streamflow data represent areas draining from 20.2 to 6,780 square miles. The PRMS model simulates reservoir inflow and watershed runoff for use as input into HECResSim for the purpose of evaluating and comparing the effects of different watershed conditions on main-stem flooding in the Delaware River watershed draining to Trenton, N.J. The PRMS model is useful as a planning tool to simulate the effects of land-use changes and different antecedent conditions on local runoff and reservoir inflow and, as input to the HEC-ResSim model, on flood flows in the main stem of the Delaware River.

Open-File Report↗

USGS2011

Improved earthquake monitoring in the central and eastern United States in support of seismic assessments for critical facilities

Evaluation of seismic monitoring capabilities in the central and eastern United States for critical facilities - including nuclear powerplants - focused on specific improvements to understand better the seismic hazards in the region. The report is not an assessment of seismic safety at nuclear plants. To accomplish the evaluation and to provide suggestions for improvements using funding from the American Recovery and Reinvestment Act of 2009, the U.S. Geological Survey examined addition of new strong-motion seismic stations in areas of seismic activity and addition of new seismic stations near nuclear power-plant locations, along with integration of data from the Transportable Array of some 400 mobile seismic stations. Some 38 and 68 stations, respectively, were suggested for addition in active seismic zones and near-power-plant locations. Expansion of databases for strong-motion and other earthquake source-characterization data also was evaluated. Recognizing pragmatic limitations of station deployment, augmentation of existing deployments provides improvements in source characterization by quantification of near-source attenuation in regions where larger earthquakes are expected. That augmentation also supports systematic data collection from existing networks. The report further utilizes the application of modeling procedures and processing algorithms, with the additional stations and the improved seismic databases, to leverage the capabilities of existing and expanded seismic arrays.

Open-File Report↗

USGS2011

Environmental flow allocation and statistics calculator

The Environmental Flow Allocation and Statistics Calculator (EFASC) is a computer program that calculates hydrologic statistics based on a time series of daily streamflow values. EFASC will calculate statistics for daily streamflow in an input file or will generate synthetic daily flow series from an input file based on rules for allocating and protecting streamflow and then calculate statistics for the synthetic time series. The program reads dates and daily streamflow values from input files. The program writes statistics out to a series of worksheets and text files. Multiple sites can be processed in series as one run. EFASC is written in MicrosoftRegistered Visual BasicCopyright for Applications and implemented as a macro in MicrosoftOffice Excel 2007Registered. EFASC is intended as a research tool for users familiar with computer programming. The code for EFASC is provided so that it can be modified for specific applications. All users should review how output statistics are calculated and recognize that the algorithms may not comply with conventions used to calculate streamflow statistics published by the U.S. Geological Survey.

Open-File Report↗

USGS2012

A multiple-point geostatistical method for characterizing uncertainty of subsurface alluvial units and its effects on flow and transport

This report provides a proof-of-concept to demonstrate the potential application of multiple-point geostatistics for characterizing geologic heterogeneity and its effect on flow and transport simulation. The study presented in this report is the result of collaboration between the U.S. Geological Survey (USGS) and Stanford University. This collaboration focused on improving the characterization of alluvial deposits by incorporating prior knowledge of geologic structure and estimating the uncertainty of the modeled geologic units. In this study, geologic heterogeneity of alluvial units is characterized as a set of stochastic realizations, and uncertainty is indicated by variability in the results of flow and transport simulations for this set of realizations. This approach is tested on a hypothetical geologic scenario developed using data from the alluvial deposits in Yucca Flat, Nevada. Yucca Flat was chosen as a data source for this test case because it includes both complex geologic and hydrologic characteristics and also contains a substantial amount of both surface and subsurface geologic data. Multiple-point geostatistics is used to model geologic heterogeneity in the subsurface. A three-dimensional (3D) model of spatial variability is developed by integrating alluvial units mapped at the surface with vertical drill-hole data. The SNESIM (Single Normal Equation Simulation) algorithm is used to represent geologic heterogeneity stochastically by generating 20 realizations, each of which represents an equally probable geologic scenario. A 3D numerical model is used to simulate groundwater flow and contaminant transport for each realization, producing a distribution of flow and transport responses to the geologic heterogeneity. From this distribution of flow and transport responses, the frequency of exceeding a given contaminant concentration threshold can be used as an indicator of uncertainty about the location of the contaminant plume boundary.

Open-File Report↗

USGS2012

Dissolved oxygen analysis, TMDL model comparison, and particulate matter shunting—Preliminary results from three model scenarios for the Klamath River upstream of Keno Dam, Oregon

Efforts are underway to identify actions that would improve water quality in the Link River to Keno Dam reach of the Upper Klamath River in south-central Oregon. To provide further insight into water-quality improvement options, three scenarios were developed, run, and analyzed using previously calibrated CE-QUAL-W2 hydrodynamic and water-quality models. Additional scenarios are under development as part of this ongoing study. Most of these scenarios evaluate changes relative to a "current conditions" model, but in some cases a "natural conditions" model was used that simulated the reach without the effect of point and nonpoint sources and set Upper Klamath Lake at its Total Maximum Daily Load (TMDL) targets. These scenarios were simulated using a model developed by the U.S. Geological Survey (USGS) and Watercourse Engineering, Inc. for the years 2006–09, referred to here as the "USGS model." Another model of the reach was developed by Tetra Tech, Inc. for years 2000 and 2002 to support the Klamath River TMDL process; that model is referred to here as the "TMDL model." The three scenarios described in this report included (1) an analysis of whether this reach of the Upper Klamath River would be in compliance with dissolved oxygen standards if sources met TMDL allocations, (2) an application of more recent datasets to the TMDL model with comparison to results from the USGS model, and (3) an examination of the effect on dissolved oxygen in the Klamath River if particulate material were stopped from entering Klamath Project diversion canals. Updates and modifications to the USGS model are in progress, so in the future these scenarios will be reanalyzed with the updated model and the interim results presented here will be superseded. Significant findings from this phase of the investigation include: * The TMDL analysis used depth-averaged dissolved oxygen concentrations from model output for comparison with dissolved oxygen standards. The Oregon dissolved oxygen standards do not specify whether the numeric criteria are based on depth-averaged dissolved oxygen concentration; this was an interpretation of the standards rule by the Oregon Department of Environmental Quality (ODEQ). In this study, both depth-averaged and volume-averaged dissolved oxygen concentrations were calculated from model output. Results showed that modeled depth-averaged concentrations typically were lower than volume-averaged dissolved oxygen concentrations because depth-averaging gives a higher weight to small volume areas near the channel bottom that often have lower dissolved oxygen concentrations. Results from model scenarios in this study are reported using volume-averaged dissolved oxygen concentrations. * Under all scenarios analyzed, violations of the dissolved oxygen standard occurred most often in summer. Of the three dissolved oxygen criteria that must be met, the 30-day standard was violated most frequently. Under the base case (current conditions), fewer violations occurred in the upstream part of the reach. More violations occurred in the down-stream direction, due in part to oxygen demand from the decay of algae and organic matter from Link River and other inflows. * A condition in which Upper Klamath Lake and its Link River outflow achieved Upper Klamath Lake TMDL water-quality targets was most effective in reducing the number of violations of the dissolved oxygen standard in the Link River to Keno Dam reach of the Klamath River. The condition in which point and nonpoint sources within the Link River to Keno Dam reach met Klamath River TMDL allocations had no effect on dissolved oxygen compliance in some locations and a small effect in others under current conditions. On the other hand, meeting TMDL allocations for nonpoint and point sources was predicted to be important in meeting dissolved oxygen criteria when Upper Klamath Lake and Link River also met Upper Klamath TMDL water-quality targets. * The location of greatest dissolved oxygen improvement from nutrient and organic matter reductions was downstream from point and nonpoint source inflows because time and distance are required for decay to occur and for oxygen demand to be exerted. * After assessing compliance with dissolved oxygen standards at all 102 model segments in the Link River to Keno Dam reach, it was determined that the seven locations used by ODEQ appear to be a representative subset of the reach for dissolved oxygen analysis. * The USGS and TMDL models were qualitatively compared by running both models for the 2006–09 period but preserving the essential characteristics of each, such as organic matter partitioning, bathymetric representation, and parameter rates. The analysis revealed that some constituents were not greatly affected by the differing algorithms, rates, and assumptions in the two models. Conversely, other constituents, especially organic matter, were simulated differently by the two models. Organic matter in this river system is best represented by a mixture of relatively labile particulate material and a substantial concentration of refractory dissolved material. In addition, the use of a first-order sediment oxygen demand, as in the USGS model, helps to capture the seasonal and dynamic effect of settled organic and algal material. * Simulation of shunting (diverting) particulate material away from the intake of four Klamath Project diversion canals, so that the material stayed in the river and out of the Project area, caused higher concentrations of particulate material to occur in the river. In all cases modeled, the increase in in-river particulate material also produced decreased dissolved oxygen concentrations and an increase in the number of days when dissolved oxygen standards were violated. * If particulate material were shunted back into the river at the Klamath Project diversion canals, less organic matter and nutrients would be taken into the Klamath Project area and the Lost River basin, resulting in return flows to the Klamath River via Lost River Diversion Channel that may have reduced nutrient concentrations. Model scenarios bracketing potential end-member nutrient concentrations showed that the composition of the return flows had little to no effect on dissolved oxygen compliance under simulated conditions.

Oregon↗

USGS2012

Thermodynamic method for generating random stress distributions on an earthquake fault

This report presents a new method for generating random stress distributions on an earthquake fault, suitable for use as initial conditions in a dynamic rupture simulation. The method employs concepts from thermodynamics and statistical mechanics. A pattern of fault slip is considered to be analogous to a micro-state of a thermodynamic system. The energy of the micro-state is taken to be the elastic energy stored in the surrounding medium. Then, the Boltzmann distribution gives the probability of a given pattern of fault slip and stress. We show how to decompose the system into independent degrees of freedom, which makes it computationally feasible to select a random state. However, due to the equipartition theorem, straightforward application of the Boltzmann distribution leads to a divergence which predicts infinite stress. To avoid equipartition, we show that the finite strength of the fault acts to restrict the possible states of the system. By analyzing a set of earthquake scaling relations, we derive a new formula for the expected power spectral density of the stress distribution, which allows us to construct a computer algorithm free of infinities. We then present a new technique for controlling the extent of the rupture by generating a random stress distribution thousands of times larger than the fault surface, and selecting a portion which, by chance, has a positive stress perturbation of the desired size. Finally, we present a new two-stage nucleation method that combines a small zone of forced rupture with a larger zone of reduced fracture energy.

Open-File Report↗

USGS2012

Producing fractional rangeland component predictions in a sagebrush ecosystem, a Wyoming sensitivity analysis

Remote sensing information has been widely used to monitor vegetation condition and variations in a variety of ecosystems, including shrublands. Careful application of remotely sensed imagery can provide additional spatially explicit, continuous, and extensive data on the composition and condition of shrubland ecosystems. Historically, the most widely available remote sensing information has been collected by Landsat, which has offered large spatial coverage and moderate spatial resolution data globally for nearly three decades. Such medium-resolution satellite remote sensing information can quantify the distribution and variation of terrestrial ecosystems. Landsat imagery has been frequently used with other high-resolution remote sensing data to classify sagebrush components and quantify their spatial distributions (Ramsey and others, 2004; Seefeldt and Booth, 2004; Stow and others, 2008; Underwood and others, 2007). Modeling algorithms have been developed to use field measurements and satellite remote sensing data to quantify the extent and evaluate the quality of shrub ecosystem components in large geographic areas (Homer and others, 2009). The percent cover of sagebrush ecosystem components, including bare-ground, herbaceous, litter, sagebrush, and shrub, have been quantified for entire western states (Homer and others, 2012). Furthermore, research has demonstrated the use of current measurements with historical archives of Landsat imagery to quantify the variations of these components for the last two decades (Xian and others, 2012). The modeling method used to quantify the extent and spatial distribution of sagebrush components over a large area also has required considerable amounts of training data to meet targeted accuracy requirements. These training data have maintained product accuracy by ensuring that they are derived from good quality field measurements collected during appropriate ecosystem phenology and subsequently maximized by extrapolation on high-resolution remote sensing data (Homer and others, 2012). This method has proven its utility; however, to develop these products across even larger areas will require additional cost efficiencies to ensure that an adequate product can be developed for the lowest cost possible. Given the vast geographic extent of shrubland ecosystems in the western United States, identifying cost efficiencies with optimal training data development and subsequent application to medium resolution satellite imagery provide the most likely areas for methodological efficiency gains. The primary objective of this research was to conduct a series of sensitivity tests to evaluate the most optimal and practical way to develop Landsat scale information for estimating the extent and distribution of sagebrush ecosystem components over large areas in the conterminous United States. An existing dataset of sagebrush components developed from extensive field measurements, high-resolution satellite imagery, and medium resolution Landsat imagery in Wyoming was used as the reference database (Homer and others, 2012). Statistical analysis was performed to analyze the relation between the accuracy of sagebrush components and the amount and distribution of training data on Landsat scenes needed to obtain accurate predictions.

Wyoming↗

USGS2013

Hydrothermal alteration maps of the central and southern Basin and Range Province of the United States compiled from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data

Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data and Interactive Data Language (IDL) logical operator algorithms were used to map hydrothermally altered rocks in the central and southern parts of the Basin and Range province of the United States. The hydrothermally altered rocks mapped in this study include (1) hydrothermal silica-rich rocks (hydrous quartz, chalcedony, opal, and amorphous silica), (2) propylitic rocks (calcite-dolomite and epidote-chlorite mapped as separate mineral groups), (3) argillic rocks (alunite-pyrophyllite-kaolinite), and (4) phyllic rocks (sericite-muscovite). A series of hydrothermal alteration maps, which identify the potential locations of hydrothermal silica-rich, propylitic, argillic, and phyllic rocks on Landsat Thematic Mapper (TM) band 7 orthorectified images, and geographic information systems shape files of hydrothermal alteration units are provided in this study.

central and southern Basin and Range Province↗

USGS2014

Hydrologic Drought Decision Support System (HyDroDSS)

The hydrologic drought decision support system (HyDroDSS) was developed by the U.S. Geological Survey (USGS) in cooperation with the Rhode Island Water Resources Board (RIWRB) for use in the analysis of hydrologic variables that may indicate the risk for streamflows to be below user-defined flow targets at a designated site of interest, which is defined herein as data-collection site on a stream that may be adversely affected by pumping. Hydrologic drought is defined for this study as a period of lower than normal streamflows caused by precipitation deficits and (or) water withdrawals. The HyDroDSS is designed to provide water managers with risk-based information for balancing water-supply needs and aquatic-habitat protection goals to mitigate potential effects of hydrologic drought. This report describes the theory and methods for retrospective streamflow-depletion analysis, rank correlation analysis, and drought-projection analysis. All three methods are designed to inform decisions made by drought steering committees and decisionmakers on the basis of quantitative risk assessment. All three methods use estimates of unaltered streamflow, which is the measured or modeled flow without major withdrawals or discharges, to approximate a natural low-flow regime. Retrospective streamflow-depletion analysis can be used by water-resource managers to evaluate relations between withdrawal plans and the potential effects of withdrawal plans on streams at one or more sites of interest in an area. Retrospective streamflow-depletion analysis indicates the historical risk of being below user-defined flow targets if different pumping plans were implemented for the period of record. Retrospective streamflow-depletion analysis also indicates the risk for creating hydrologic drought conditions caused by use of a pumping plan. Retrospective streamflow-depletion analysis is done by calculating the net streamflow depletions from withdrawals and discharges and applying these depletions to a simulated record of unaltered streamflow. Rank correlation analysis in the HyDroDSS indicates the persistence of hydrologic measurements from month to month for the prediction of developing hydrologic drought conditions and quantitatively indicates which hydrologic variables may be used to indicate the onset of hydrologic drought conditions. Rank correlation analysis also indicates the potential use of each variable for estimating the monthly minimum unaltered flow at a site of interest for use in the drought-projection analysis. Rank correlation analysis in the HyDroDSS is done by calculating Spearman’s rho for paired samples and the 95-percent confidence limits of this rho value. Rank correlation analysis can be done by using precipitation, groundwater levels, measured streamflows, and estimated unaltered streamflows. Serial correlation analysis, which indicates relations between current and future values, can be done for a single site. Cross correlation analysis, which indicates relations among current values at one site and current and future values at a second site, also can be done. Drought-projection analysis in the HyDroDSS indicates the risk for being in a hydrologic drought condition during the current month and the five following months with and without pumping. Drought-projection analysis also indicates the potential effectiveness of water-conservation methods for mitigating the effect of withdrawals in the coming months on the basis of the amount of depletion caused by different pumping plans and on the risk of unaltered flows being below streamflow targets. Drought-projection analysis in the HyDroDSS is done with Monte Carlo methods by using the position analysis method. In this method the initial value of estimated unaltered streamflows is calculated by correlation to a measured hydrologic variable (monthly precipitation, groundwater levels, or streamflows from an index station identified with the rank correlation analysis). Then a pseudorandom number generator is used to create 251 six-month-long flow traces by using a bootstrap method. Serial correlation of the estimated unaltered monthly minimum streamflows determined from the rank correlation analysis is preserved within each flow trace. The sample of unaltered streamflows indicates the risk of being below flow targets in the coming months under simulated natural conditions (without historic withdrawals). The streamflow-depletion algorithms are then used to estimate risks of flow being below targets if selected pumping plans are used. This report also describes the implementation of the HyDroDSS. The HyDroDSS was developed as a Microsoft Access® database application to facilitate storage, handling, and use of hydrologic datasets with a simple graphical user interface. The program is implemented in the database by using the Visual Basic for Applications® (VBA) programming language. Program source code for the analytical techniques is provided in the HyDroDSS and in electronic text files accompanying this report. Program source code for the graphical user interface and for data-handling code, which is specific to Microsoft Access® and the HyDroDSS, is provided in the database. An installation package with a run-time version of the software is available with this report for potential users who do not have a compatible copy of Microsoft Access®. Administrative rights are needed to install this version of the HyDroDSS. A case study, to demonstrate the use of HyDroDSS and interpretation of results for a site of interest, is detailed for the USGS streamgage on the Hunt River (station 01117000) near East Greenwich in central Rhode Island. The Hunt River streamgage was used because it has a long record of streamflow and is in a well-studied basin with a substantial amount of hydrologic and water-use data including groundwater pumping for municipal water supply.

Rhode Island↗

USGS2014

Laharz_py: GIS tools for automated mapping of lahar inundation hazard zones

Laharz_py is written in the Python programming language as a suite of tools for use in ArcMap Geographic Information System (GIS). Primarily, Laharz_py is a computational model that uses statistical descriptions of areas inundated by past mass-flow events to forecast areas likely to be inundated by hypothetical future events. The forecasts use physically motivated and statistically calibrated power-law equations that each has a form A = cV 2/3 , relating mass-flow volume (V) to planimetric or cross-sectional areas (A) inundated by an average flow as it descends a given drainage. Calibration of the equations utilizes logarithmic transformation and linear regression to determine the best-fit values of c. The software uses values of V, an algorithm for idenitifying mass-flow source locations, and digital elevation models of topography to portray forecast hazard zones for lahars, debris flows, or rock avalanches on maps. Laharz_py offers two methods to construct areas of potential inundation for lahars: (1) Selection of a range of plausible V values results in a set of nested hazard zones showing areas likely to be inundated by a range of hypothetical flows; and (2) The user selects a single volume and a confidence interval for the prediction. In either case, Laharz_py calculates the mean expected A and B value from each user-selected value of V. However, for the second case, a single value of V yields two additional results representing the upper and lower values of the confidence interval of prediction. Calculation of these two bounding predictions require the statistically calibrated prediction equations, a user-specified level of confidence, and t-distribution statistics to calculate the standard error of regression, standard error of the mean, and standard error of prediction. The portrayal of results from these two methods on maps compares the range of inundation areas due to prediction uncertainties with uncertainties in selection of V values. The Open-File Report document contains an explanation of how to install and use the software. The Laharz_py software includes an example data set for Mount Rainier, Washington. The second part of the documentation describes how to use all of the Laharz_py tools in an example dataset at Mount Rainier, Washington.

Washington↗

USGS2014

Technical implementation plan for the ShakeAlert production system: an Earthquake Early Warning system for the West Coast of the United States

Earthquake Early Warning (EEW) systems can provide as much as tens of seconds of warning to people and automated systems before strong shaking arrives. The United States Geological Survey (USGS) and its partners are developing such an EEW system, called ShakeAlert, for the West Coast of the United States. This document describes the technical implementation of that system, which leverages existing stations and infrastructure of the Advanced National Seismic System (ANSS) regional networks to achieve this new capability. While significant progress has been made in developing the ShakeAlert early warning system, improved robustness of each component of the system and additional testing and certification are needed for the system to be reliable enough to issue public alerts. Major components of the system include dense networks of ground motion sensors, telecommunications from those sensors to central processing systems, algorithms for event detection and alert creation, and distribution systems to alert users. Capital investment costs for a West Coast EEW system are projected to be $38.3M, with additional annual maintenance and operations totaling $16.1M—in addition to current ANSS expenditures for earthquake monitoring. An EEW system is complementary to, but does not replace, other strategies to mitigate earthquake losses. The system has limitations: false and missed alerts are possible, and the area very near to an earthquake epicenter may receive little or no warning. However, such an EEW system would save lives, reduce injuries and damage, and improve community resilience by reducing longer-term economic losses for both public and private entities.

California;Oregon;Washington↗

USGS2014

Stochastic modeling of a lava-flow aquifer system

This report describes preliminary three-dimensional geostatistical modeling of a lava-flow aquifer system using a multiple-point geostatistical model. The purpose of this study is to provide a proof-of-concept for this modeling approach. An example of the method is demonstrated using a subset of borehole geologic data and aquifer test data from a portion of the Calico Hills Formation, a lava-flow aquifer system that partially underlies Pahute Mesa, Nevada. Groundwater movement in this aquifer system is assumed to be controlled by the spatial distribution of two geologic units—rhyolite lava flows and zeolitized tuffs. The configuration of subsurface lava flows and tuffs is largely unknown because of limited data. The spatial configuration of the lava flows and tuffs is modeled by using a multiple-point geostatistical simulation algorithm that generates a large number of alternative realizations, each honoring the available geologic data and drawn from a geologic conceptual model of the lava-flow aquifer system as represented by a training image. In order to demonstrate how results from the geostatistical model could be analyzed in terms of available hydrologic data, a numerical simulation of part of an aquifer test was applied to the realizations of the geostatistical model.

Nevada↗

USGS2015

National Earthquake Information Center systems overview and integration

Overview The primary mission of the U.S. Geological Survey (USGS) National Earthquake Information Center (NEIC) is comprehensive global earthquake monitoring (M4.5 or larger) and complete seismic monitoring of the United States for all significant earthquakes (M3.0 or larger or felt). In recent years, the NEIC has assumed a more prominent role in local and regional seismic monitoring, backup capabilities, and special studies. The NEIC uses numerous scientific algorithms to produce information products for local, national, and international disaster responders; infrastructure corporations (for example, utility providers and insurers); U.S. embassies worldwide; news agencies; national and international earth science research communities; and millions of public users. To support this mission, the NEIC operates a constant service to derive numerous information products as rapidly and accurately as possible. Derived information products and data sets provide a wealth of related earthquake information and include products such as location, magnitude, maps of shaking intensity, moment tensors, phase arrivals, and loss and damage reports. These information products are automatically derived in near real-time and refined over time. Earthquake information products vary from well formatted text files, to binary files, to directories full of many files. A team of analysts, who provide constant coverage, 24 hours a day, seven days a week, 365 days a year, review the acquired data and automatically derived products (1) in near real-time for rapid release and (2) in the weeks following events to produce historical catalogs of global earthquakes in support of scientific research. Data and products are reviewed both in near real-time for rapid release and in the weeks following events to produce historical catalogs of global earthquakes in support of scientific research. At the present time, the NEIC publishes data for approximately 25,000 earthquakes annually. The NEIC software system that supports this service and mission comprises many subsystems. The internals of each subsystem are often complex, but the integration between all of the subsystems is well defined, standardized, and straightforward. This Open-File Report describes the systems’ integration infrastructure and how all the subsystems fit together. It provides summary views of the system as a whole, along with detailed views of the inputs, outputs, and dependencies for each subsystem. It also outlines the many shared data services that support these subsystems and applications as well as external customers. The NEIC systems-integration effort is ongoing. Additional integration steps, needed to reach our goals, are discussed. It is important to note that this document provides a brief introduction to the work of dozens of software developers and IT specialists, spanning in many cases more than a decade. References to significant amounts of supporting documentation, code, and information are supplied within.

Open-File Report↗