Search USGS⌕ Search

SEARCH · Search USGS

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 649 records · Page 36Linked to original sources
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2013

A scenario study of seismically induced landsliding in Seattle using broadband synthetic seismograms

We demonstrate the value of utilizing broadband synthetic seismograms to assess regional seismically induced landslide hazard. Focusing on a case study of an M w 7.0 Seattle fault earthquake in Seattle, Washington, we computed broadband synthetic seismograms that account for rupture directivity and 3D basin amplification. We then adjusted the computed motions on a fine grid for 1D amplifications based on the site response of typical geologic profiles in Seattle and used these time‐series ground motions to trigger shallow landsliding using the Newmark method. The inclusion of these effects was critical in determining the extent of landsliding triggered. We found that for inertially triggered slope failures modeled by the Newmark method, the ground motions used to simulate landsliding must have broadband frequency content in order to capture the full slope displacement. We applied commonly used simpler methods based on ground‐motion prediction equations for the same scenario and found that they predicted far fewer landslides if only the mean values were used, but far more at the maximum range of the uncertainties, highlighting the danger of using just the mean values for such methods. Our results indicate that landsliding triggered by a large Seattle fault earthquake will be extensive and potentially devastating, causing direct losses and impeding recovery. The high impact of landsliding predicted by this simulation shows that this secondary effect of earthquakes should be studied with as much vigor as other earthquake effects.

Washington↗
USGS1986

Preliminary map showing the thickness of glacial deposits in Ohio

This map was compiled as part of a U.S. Geological Survey project to portray the character and thickness of glacial deposits east of the Rocky Mountains. It is a preliminary map, compiled at a scale of 1:250,000 from sources of differing type and quality (fig. I, table 1). County drift-thickness maps produced by the Ohio Geological Survey are the most reliable source of data used. These county maps, either published or in manuscript form, were available for 36 of the 72 counties that contain glacial or glacially related deposits more than 50 ft thick. An earlier, reconnaissance drift-thickness map of northwestern Ohio (Ohio Department of Natural Resources, 1960) provided data for parts of 12 additional counties. Published maps of bedrock topography were utilized for 9 counties and for parts of 6 counties. The computation of drift thickness in these counties is less reliable then the drift-thickness maps published by the Ohio Geological Survey because the bedrock topography maps, most of which are at least 20 years old, were compiled from a more limited set of data. In the remaining counties, only well-log data were available. Drift-thickness data in these counties are the least reliable; the determination of thickness trends was commonly guided by the configuration of preglacial drainage channels (as drawn by Stout and others, 1943) and by topography. Throughout the glaciated areas of Ohio, and especially near the glacial margin, a network of drainage channels is buried by glacial drift. Some of these channels arc of a preglacial age while others were formed during interglacial periods; some were cut into the bedrock by rivers— either preglacial rivers or meltwater streams flowing along the glacial margin—while others, such as the lower Cuyahoga River valley, may have been overdeepened by the erosive power of glacial ice. The thickest deposits of drift in the state occur in the lower Cuyahoga River valley near Lake Erie, and above the buried valley in west-central Ohio that has been historically considered a part of the Teays Valley. The thickest drift reliably measured in Ohio to date is 602 ft in a drill hole near Cleveland. In contrast to the extreme variations in drift thickness encountered in the vicinity of buried channels, drift on the upland arcus is generally thinner and the variations in thickness are much less pronounced. Worthy of note, however, are three large areas where the drift sheet is relatively thick. In northwestern Ohio, a large volume of drift was deposited along the flanks of the Erie ice lobe (fig. 2) near the interlobate position with the Saginaw lobe to the northwest; drift thickness there exceeds 200 ft. Thick drift was also deposited in a roughly cast-west band across the Miami lobe. The mechanism that produced this band of thick drift is not obvious, but it may have been influenced in part by bedrock topography. Bedrock control of drift thickness is more clearly indicated to the cast of Columbus, along the eastern flank of the Scioto lobe, where ice slow was resisted by rocks of the Allegheny plateau. The edge of the plateau, or the Allegheny escarpment, is obscured by glacial deposits but its likely position (Fenneman, 1938; Stout and others, 1913; Dove, 1960; and Root and others, 1961) is shown on the map. Southward from the ice margin's reentrant position in southern Richland County, ice flowing eastward from the Scioto lobe encountered the topographically higher plateau, which constrained the ice and caused drift to accumulate in significant thicknesses just to the west of the escarpment.

Ohio↗
USGS2005

Estimation of agricultural pesticide use in drainage basins using land cover maps and county pesticide data

A geographic information system (GIS) was used to estimate agricultural pesticide use in the drainage basins of streams that are studied as part of the U.S. Geological Survey?s National Water-Quality Assessment (NAWQA) Program. Drainage basin pesticide use estimates were computed by intersecting digital maps of drainage basin boundaries with an enhanced version of the National Land Cover Data 1992 combined with estimates of 1992 agricultural pesticide use in each United States county. This report presents the methods used to quantify agricultural pesticide use in drainage basins using a GIS and includes the estimates of atrazine use applied to row crops, small-grain crops, and fallow lands in 150 watersheds in the conterminous United States. Basin atrazine use estimates are presented to compare and analyze the results that were derived from 30-meter and 1-kilometer resolution land cover and county pesticide use data, and drainage basin boundaries at various grid cell resolutions. Comparisons of the basin atrazine use estimates derived from watershed boundaries, county pesticide use, and land cover data sets at different resolutions, indicated that overall differences were minor. The largest potential for differences in basin pesticide use estimates between those derived from the 30-meter and 1-kilometer resolution enhanced National Land Cover Data 1992 exists wherever there are abrupt agricultural land cover changes along the basin divide. Despite the limitations of the drainage basin pesticide use data described in this report, the basin estimates provide consistent and comparable indicators of agricultural pesticide application in surface-water drainage basins studied in the NAWQA Program.

Open-File Report↗
USGS2017

National assessment of shoreline change—Summary statistics for updated vector shorelines and associated shoreline change data for the Gulf of Mexico and Southeast Atlantic coasts

Long-term rates of shoreline change for the Gulf of Mexico and Southeast Atlantic regions of the United States have been updated as part of the U.S. Geological Survey’s National Assessment of Shoreline Change project. Additional shoreline position data were used to compute rates where the previous rate-of-change assessment only included four shoreline positions at a given location. The long-term shoreline change rates also incorporate the proxy-datum bias correction to account for the unidirectional onshore bias of the proxy-based high water line shorelines relative to the datum-based mean high water shorelines. The calculation of uncertainty associated with the long-term average rates has also been updated to match refined methods used in other study regions of the National Assessment project. The average rates reported here have a reduced amount of uncertainty relative to those presented in the previous assessments for these two regions.

Alabama, Florida, Georgia, Louisiana, Mississippi,↗
USGS1983

MAGPOLY; a modification of a three-dimensional magnetic modelling program

The three-dimensional magnetic modelling program of Donald Plouff (Plouff, 1975) has been made operational on a Honeywell Multics 68/80 computer. Modifications were made to input and output procedures to make the program compatible with U.S. Geological Survey (USGS) standard geophysical input and output formats.

Open-File Report↗
USGS1983

GRAVPOLY; a modification of a three-dimensional gravity modeling program

The three-dimensional magnetic modelling program of Donald Plouff (Plouff, 1975) has been made operational on a Honeywell Multics 68/80 computer. Modifications were made to input and output procedures to make the program compatible with U.S. Geological Survey (USGS) standard geophysical input and output formats.

Open-File Report↗
USGS2004

Estimating the Magnitude and Frequency of Floods in Small Urban Streams in South Carolina, 2001

The magnitude and frequency of floods at 20 streamflowgaging stations on small, unregulated urban streams in or near South Carolina were estimated by fitting the measured wateryear peak flows to a log-Pearson Type-III distribution. The period of record (through September 30, 2001) for the measured water-year peak flows ranged from 11 to 25 years with a mean and median length of 16 years. The drainage areas of the streamflow-gaging stations ranged from 0.18 to 41 square miles. Based on the flood-frequency estimates from the 20 streamflow-gaging stations (13 in South Carolina; 4 in North Carolina; and 3 in Georgia), generalized least-squares regression was used to develop regional regression equations. These equations can be used to estimate the 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence-interval flows for small urban streams in the Piedmont, upper Coastal Plain, and lower Coastal Plain physiographic provinces of South Carolina. The most significant explanatory variables from this analysis were mainchannel length, percent impervious area, and basin development factor. Mean standard errors of prediction for the regression equations ranged from -25 to 33 percent for the 10-year recurrence-interval flows and from -35 to 54 percent for the 100-year recurrence-interval flows. The U.S. Geological Survey has developed a Geographic Information System application called StreamStats that makes the process of computing streamflow statistics at ungaged sites faster and more consistent than manual methods. This application was developed in the Massachusetts District and ongoing work is being done in other districts to develop a similar application using streamflow statistics relative to those respective States. Considering the future possibility of implementing StreamStats in South Carolina, an alternative set of regional regression equations was developed using only main channel length and impervious area. This was done because no digital coverages are currently available for basin development factor and, therefore, it could not be included in the StreamStats application. The average mean standard error of prediction for the alternative equations was 2 to 5 percent larger than the standard errors for the equations that contained basin development factor. For the urban streamflow-gaging stations in South Carolina, measured water-year peak flows were compared with those from an earlier urban flood-frequency investigation. The peak flows from the earlier investigation were computed using a rainfall-runoff model. At many of the sites, graphical comparisons indicated that the variance of the measured data was much less than the variance of the simulated data. Several statistical tests were applied to compare the variances and the means of the measured and simulated data for each site. The results indicated that the variances were significantly different for 11 of the 13 South Carolina streamflow-gaging stations. For one streamflow-gaging station, the test for normality, which is one of the assumptions of the data when comparing variances, indicated that neither the measured data nor the simulated data were distributed normally; therefore, the test for differences in the variances was not used for that streamflow-gaging station. Another statistical test was used to test for statistically significant differences in the means of the measured and simulated data. The results indicated that for 5 of the 13 urban streamflowgaging stations in South Carolina there was a statistically significant difference in the means of the two data sets. For comparison purposes and to test the hypothesis that there may have been climatic differences between the period in which the measured peak-flow data were measured and the period for which historic rainfall data were used to compute the simulated peak flows, 16 rural streamflow-gaging stations with long-term records were reviewed using similar techniques as those used for the measured an

South Carolina↗
USGS2005

The Modular Modeling System (MMS): A toolbox for water- and environmental-resources management

The increasing complexity of water- and environmental-resource problems require modeling approaches that incorporate knowledge from a broad range of scientific and software disciplines. To address this need, the U.S. Geological Survey (USGS) has developed the Modular Modeling System (MMS). MMS is an integrated system of computer software for model development, integration, and application. Its modular design allows a high level of flexibility and adaptability to enable modelers to incorporate their own software into a rich array of built-in models and modeling tools. These include individual process models, tightly coupled models, loosely coupled models, and fully- integrated decision support systems. A geographic information system (GIS) interface, the USGS GIS Weasel, has been integrated with MMS to enable spatial delineation and characterization of basin and ecosystem features, and to provide objective parameter-estimation methods for models using available digital data. MMS provides optimization and sensitivity-analysis tools to analyze model parameters and evaluate the extent to which uncertainty in model parameters affects uncertainty in simulation results. MMS has been coupled with the Bureau of Reclamation object-oriented reservoir and river-system modeling framework, RiverWare, to develop models to evaluate and apply optimal resource-allocation and management strategies to complex, operational decisions on multipurpose reservoir systems and watersheds. This decision support system approach has been developed, tested, and implemented in the Gunnison, Yakima, San Joaquin, Rio Grande, and Truckee River basins of the western United States. MMS is currently being coupled with the U.S. Forest Service model SIMulating Patterns and Processes at Landscape Scales (SIMPPLLE) to assess the effects of alternative vegetation-management strategies on a variety of hydrological and ecological responses. Initial development and testing of the MMS-SIMPPLLE integration is being conducted on the Colorado Plateau region of the western United Sates.

Conference Paper↗
USGS2007

A method to estimate groundwater depletion from confining layers

Although depletion of storage in low‐permeability confining layers is the source of much of the groundwater produced from many confined aquifer systems, it is all too frequently overlooked or ignored. This makes effective management of groundwater resources difficult by masking how much water has been derived from storage and, in some cases, the total amount of water that has been extracted from an aquifer system. Analyzing confining layer storage is viewed as troublesome because of the additional computational burden and because the hydraulic properties of confining layers are poorly known. In this paper we propose a simplified method for computing estimates of confining layer depletion, as well as procedures for approximating confining layer hydraulic conductivity ( K ) and specific storage ( S s ) using geologic information. The latter makes the technique useful in developing countries and other settings where minimal data are available or when scoping calculations are needed. As such, our approach may be helpful for estimating the global transfer of groundwater to surface water. A test of the method on a synthetic system suggests that the computational errors will generally be small. Larger errors will probably result from inaccuracy in confining layer property estimates, but these may be no greater than errors in more sophisticated analyses. The technique is demonstrated by application to two aquifer systems: the Dakota artesian aquifer system in South Dakota and the coastal plain aquifer system in Virginia. In both cases, depletion from confining layers was substantially larger than depletion from the aquifers.

Water Resources Research↗
USGS1996

Review of Seismic Hazard Issues Associated with Auburn Dam Project, Sierra Nevada Foothills, California

Summary -- The U.S. Geological Survey was requested by the U.S. Department of the Interior to review the design values and the issue of reservoir-induced seismicity for a concrete gravity dam near the site of the previously-proposed Auburn Dam in the western foothills of the Sierra Nevada, central California. The dam is being planned as a flood-control-only dam with the possibility of conversion to a permanent water-storage facility. As a basis for planning studies the U.S. Army Corps of Engineers is using the same design values approved by the Secretary of the Interior in 1979 for the original Auburn Dam. These values were a maximum displacement of 9 inches on a fault intersecting the dam foundation, a maximum earthquake at the site of magnitude 6.5, a peak horizontal acceleration of 0.64 g, and a peak vertical acceleration of 0.39 g. In light of geological and seismological investigations conducted in the western Sierran foothills since 1979 and advances in the understanding of how earthquakes are caused and how faults behave, we have developed the following conclusions and recommendations: Maximum Displacement. Neither the pre-1979 nor the recent observations of faults in the Sierran foothills precisely define the maximum displacement per event on a fault intersecting the dam foundation. Available field data and our current understanding of surface faulting indicate a range of values for the maximum displacement. This may require the consideration of a design value larger than 9 inches. We recommend reevaluation of the design displacement using current seismic hazard methods that incorporate uncertainty into the estimate of this design value. Maximum Earthquake Magnitude. There are no data to indicate that a significant change is necessary in the use of an M 6.5 maximum earthquake to estimate design ground motions at the dam site. However, there is a basis for estimating a range of maximum magnitudes using recent field information and new statistical fault relations. We recommend reevaluating the maximum earthquake magnitude using current seismic hazard methodology. Design Ground Motions. A large number of strong-motion records have been acquired and significant advances in understanding of ground motion have been achieved since the original evaluations. The design value for peak horizontal acceleration (0.64 g) is larger than the median of one recent study and smaller than the median value of another. The value for peak vertical acceleration (0.39 g) is somewhat smaller than median values of two recent studies. We recommend a reevaluation of the design ground motions that takes into account new ground motion data with particular attention to rock sites at small source distances. Reservoir-Induced Seismicity. The potential for reservoir-induced seismicity must be considered for the Auburn Darn project. A reservoir-induced earthquake is not expected to be larger than the maximum naturally occurring earthquake. However, the probability of an earthquake may be enhanced by reservoir impoundment. A flood-control-only project may involve a lower probability of significant induced seismicity than a multipurpose water-storage dam. There is a need to better understand and quantify the likelihood of this hazard. A methodology should be developed to quantify the potential for reservoir induced seismicity using seismicity data from the Sierran foothills, new worldwide observations of induced and triggered seismicity, and current understanding of the earthquake process. Reevaluation of Design Parameters. The reevaluation of the maximum displacement, maximum magnitude earthquake, and design ground motions can be made using available field observations from the Sierran foothills, updated statistical relations for faulting and ground motions, and current computational seismic hazard methodologies that incorporate uncertainty into the analysis. The reevaluation does not require significant new geological field studies.

Open-File Report↗
USGS1984

Relationship of suspended sediment to streamflow in the Green River Basin, Wyoming

The relation of suspended-sediment concentration to water discharge was computed using regression analysis for 33 stations and 2 miscellaneous sites in the Green River basin, Wyoming. These relations show that suspended-sediment concentrations vary in direct proportion to streamflow. The number of coincident instantaneous suspended-sediment concentration and water-discharge values available for each regression analysis ranged from 6 to 98. Standard error of estimates for 35 relationships were usually large and correlation were relatively small. Large standard errors and small correlations are expected because of the complexity of the fluvial system. The effects of other factors, such as geology, topography, climate, soils, vegetation, and land use were not evaluated. Daily suspended-sediment discharges for four stations were computed by the ' sediment-transport-curve ' method described in this report , and compared with values computed by the traditional ' temporal-concentration-graph ' method. Annual sediment yields in tons per square mile drainage area, were compared in like manner. For those stations examined, values computed by the ' sediment-transport-curve ' method were within 58 percent of those computed by the temporal-concentration-graph ' method when only one year of record was used, but were within 12 to 21 percent when 2 to 4 years of record were averaged. (USGS)

Water-Resources Investigations Report↗
USGS1986

Management of the life and death of an earth-science database: some examples from geotherm

Productive earth-science databases require managers who are familiar with and skilled at using available software developed specifically for database management. There also should be a primary user with a clearly understood mission. The geologic phenomenon addressed by the database must be sufficiently understood, and adequate appropriate data must be available to construct a useful database. The database manager, in concert with the primary user, must ensure that data of adequate quality are available in the database, as well as prepare for mechanisms of releasing the data when the database is terminated. The primary user needs to be held accountable along with the database manager to ensure that a useful database will be created. Quality of data and maintenance of database relevancy to the user's mission are important issues during the database's lifetime. Products prepared at termination may be used more than the operational database and thus are of critical importance. These concepts are based, in part, on both the shortcomings and successes of GEOTHERM, a comprehensive system of databases and software used to store, locate, and evaluate the geology, geochemistry, and hydrology of geothermal systems. ?? 1986.

Computers & Geosciences↗
USGS1998

Chlorofluorocarbons, sulfur hexafluoride, and dissolved permanent gases in ground water from selected sites in and near the Idaho National Engineering and Environmental Laboratory, Idaho, 1994-97

From July 1994 through May 1997, the U.S. Geological Survey in cooperation with the Department of Energy, sampled 86 wells completed in the Snake River Plain aquifer at and near the Idaho N ationa1 Engineering and Environmental Laboratory (INEEL). The wells were sampled for a variety of constituents including one- and two-carbon halocarbons. Concentrations of dichlorodifluoromethane (CFC-12), trichlorofluoromethane (CFC-11) and trichlorotrifluororoethane (CFC-113) were determined. The samples for halocarbon analysis were collected in 62-milliliter flame sealed borosilicate glass ampoules in the field. The data will be used to evaluate the ages of ground waters at INEEL. The ages of the ground water will be used to determine recharge rates, residence time, and travel time of water in the Snake River Plain aquifer in and near INEEL. The chromatograms of 139 ground waters are presented showing a large number of halomethanes, haloethanes, and haloethenes present in the ground waters underlying the INEEL. The chromatograms can be used to qualitatively evaluate a large number of contaminants at parts per trillion to parts per billion concentrations. The data can be used to study temporal and spatial distribution of contaminants in the Snake River Plain aquifer. Representative compressed chromatograms for all ground waters sampled in this study are available on two 3.5-inch high density computer disks. The data and the program required to decompress the data can be obtained from the U.S. Geological Survey office at Idaho Falls, Idaho. Sulfur hexafluoride (SF6) concentrations were measured in selected wells to determine the feasibility of using this environmental tracer as an age dating tool of ground water. Concentrations of dissolved nitrogen, argon, carbon dioxide, oxygen, and methane were measured in 79 ground waters. Concentrations of dissolved permanent gases are tabulated and will be used to evaluate the temperature of recharge of ground water in and near the INEEL.

Idaho↗
USGS1993

Regression modeling of ground-water flow; Supplement 1 : modifications to the computer code for nonlinear regression solution of steady-state ground-water flow problems

Modifications to the computer code for nonlinear regression solution of steady-state ground-water flow problems in "Regression Modeling of Ground-Water Flow", by R.L. Cooley and R.L. Naff (U.S. Geological Survey, Techniques of Water Resources Investigations, book 3, chap. B4, 1990), improve the performance of the code. The changes 1) allow any model parameters to be transformed to natural logarithms, and 2) improve the procedure for computing the parameter that damps changes in the values of model parameters for each iteration of the solution method. In addition, the user is shown how to read prior information separately from the initial set of parameter values.

Techniques of Water-Resources Investigations↗
USGS1982

Maps showing ground-water conditions in the Bill Williams area, Mohave, Yavapai, and Yuma Counties, Arizona - 1980

The Bill Williams area includes about 3,200 mi 2 in Mohave, Yavapai, and Yuma Counties in west-central Arizona. The west half of the area is in the Basin and Range lowlands water province, and the east half is in the Central high-lands water province (see index map). The Basin and Range lowlands province generally is characterized by high mountains separated by broad valleys filled with deposits that commonly store large amounts of ground water. The Central highlands province consists mostly of rugged mountain masses made up of igneous, metamorphic, and well-consolidated sedimentary rocks that contain little space for the storage of ground water except where highly fractured or faulted. A few small valleys between the mountains contain varying thicknesses of water.-bearing deposits. The area is drained by the Bill Williams River and its major tributaries-the Big Sandy River and the Santa Maria River. Many reaches of the Big Sandy and Santa Maria Rivers and their major tributaries are perennial; the flow is sustained by ground-water discharge (Brown and others, 1978, sheet 2). In the Bill Williams area most of the water used is from ground water, although a small amount of surface water also may be diverted. About 18,000 acre-ft of ground water was withdrawn in 1979 (U.S. Geological Survey, 1981). About 17,000 acre-ft was used for the irrigation of 5,200 acres, and the rest was used for domestic, stock, and public supplies. Most of the irrigated land is in Skull Valley and along lower Kirkland Creek and the Bill Williams River. Only selected wells are shown on the maps in areas of high well density. The hydrologic data on which these maps are based are available, for the most part, in computer-printout form and may be consulted at the Arizona Department of Water Resources, 99 East Virginia, Phoenix, and at U.S. Geological Survey offices in: Federal Building, 301 West Congress Street, Tucson, and Valley Center, Suite 1880, Phoenix. Material from which copies can be made at private expense is available at the Tucson and Phoenix offices of the U.S. Geological Survey.

Arizona↗
USGS2009

Estimation of Leakage Potential of Selected Sites in Interstate and Tri-State Canals Using Geostatistical Analysis of Selected Capacitively Coupled Resistivity Profiles, Western Nebraska, 2004

With increasing demands for reliable water supplies and availability estimates, groundwater flow models often are developed to enhance understanding of surface-water and groundwater systems. Specific hydraulic variables must be known or calibrated for the groundwater-flow model to accurately simulate current or future conditions. Surface geophysical surveys, along with selected test-hole information, can provide an integrated framework for quantifying hydrogeologic conditions within a defined area. In 2004, the U.S. Geological Survey, in cooperation with the North Platte Natural Resources District, performed a surface geophysical survey using a capacitively coupled resistivity technique to map the lithology within the top 8 meters of the near-surface for 110 kilometers of the Interstate and Tri-State Canals in western Nebraska and eastern Wyoming. Assuming that leakage between the surface-water and groundwater systems is affected primarily by the sediment directly underlying the canal bed, leakage potential was estimated from the simple vertical mean of inverse-model resistivity values for depth levels with geometrically increasing layer thickness with depth which resulted in mean-resistivity values biased towards the surface. This method generally produced reliable results, but an improved analysis method was needed to account for situations where confining units, composed of less permeable material, underlie units with greater permeability. In this report, prepared by the U.S. Geological Survey in cooperation with the North Platte Natural Resources District, the authors use geostatistical analysis to develop the minimum-unadjusted method to compute a relative leakage potential based on the minimum resistivity value in a vertical column of the resistivity model. The minimum-unadjusted method considers the effects of homogeneous confining units. The minimum-adjusted method also is developed to incorporate the effect of local lithologic heterogeneity on water transmission. Seven sites with differing geologic contexts were selected following review of the capacitively coupled resistivity data collected in 2004. A reevaluation of these sites using the mean, minimum-unadjusted, and minimum-adjusted methods was performed to compare the different approaches for estimating leakage potential. Five of the seven sites contained underlying confining units, for which the minimum-unadjusted and minimum-adjusted methods accounted for the confining-unit effect. Estimates of overall leakage potential were lower for the minimum-unadjusted and minimum-adjusted methods than those estimated by the mean method. For most sites, the local heterogeneity adjustment procedure of the minimum-adjusted method resulted in slightly larger overall leakage-potential estimates. In contrast to the mean method, the two minimum-based methods allowed the least permeable areas to control the overall vertical permeability of the subsurface. The minimum-adjusted method refined leakage-potential estimation by additionally including local lithologic heterogeneity effects.

Scientific Investigations Report↗
USGS2001

Assessment method for epithermal gold deposits in northeast Washington State using weights-of-evidence GIS modeling

The weights-of-evidence analysis, a quantitative mineral resource mapping tool, is used to delineate favorable areas for epithermal gold deposits and to predict future exploration activity of the mineral industry for similar deposits in a four-county area (222 x 277 km), including the Okanogan and Colville National Forests of northeastern Washington. Modeling is applied in six steps: (1) building a spatial digital database, (2) extracting predictive evidence for a particular deposit, based on an exploration model, (3) calculating relative weights for each predictive map, (4) combining the geologic evidence maps to predict the location of undiscovered mineral resources and (5) measuring the intensity of recent exploration activity by use of mining claims on federal lands, and (6) combining mineral resource and exploration activity into an assessment model of future mining activity. The analysis is accomplished on a personal computer using ArcView GIS platform with Spatial Analyst and Weights-of-Evidence software. In accord with the descriptive model for epithermal gold deposits, digital geologic evidential themes assembled include lithologic map units, thrust faults, normal faults, and igneous dikes. Similarly, geochemical evidential themes include placer gold deposits and gold and silver analyses from stream sediment (silt) samples from National Forest lands. Fifty mines, prospects, or occurrences of epithermal gold deposits, the training set, define the appropriate a really-associated terrane. The areal (or spatial) correlation of each evidential theme with the training set yield predictor theme maps for lithology, placer sites and normal faults. The weights-of-evidence analysis disqualified the thrust fault, dike, and gold and silver silt analyses evidential themes because they lacked spatial correlation with the training set. The decision to accept or reject evidential themes as predictors is assisted by considering probabilistic data consisting of weights and contrast values calculated for themes according to areal correlation with the training sites. Predictor themes having acceptable weights and contrast values are combined into a preliminary model to predict the locations of undiscovered epithermal gold deposits. This model facilitates ranking of tracts as non-permissive, permissive or favorable categories based on exclusionary, passive, and active criteria through evaluation of probabilistic data provided by interaction of predictor themes. The method is very similar to the visual inspection method of drawing conclusions from anomalies on a manually overlain system of maps. This method serves as a model for future mineral assessment procedures because of its objective nature. To develop a model to predict future exploration activity, the locations of lode mining claims were summarized for 1980, 1985, 1990, and 1996. Land parcels containing historic claims were identified either as those with mining claims present in 1980 or valid claims present in 1985. Current claim parcels were identified as those containing valid lode claims in either 1990 or 1996. A consistent parcel contains both historic and current claims. The epithermal gold and mining claim activity models were combined into an assessment (or mineral resource-activity) model to assist in land use decisions by providing a prediction of mineral exploration activity on federal land in the next decade. Ranks in the assessment model are: (1) no activity, (2) low activity, (3) low to moderate activity, (4) moderate activity and (5) high activity.

Washington↗