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 685 records · Page 38Linked to original sources

Methods for Processing and Summarizing Time-Series Temperature Data Collected as Part of the National Water-Quality Assessment Program Studies on the Effects of Urbanization on Stream Ecosystems

Temperature data and summary statistics are presented for 256 sites in 9 metropolitan areas as part of the U.S. Geological Survey National Water-Quality Assessment Program studies of the effects of urbanization on stream ecosystems. The computer program (GRAN) that was developed to derive uniform data granularity and calculate temperature statistics (means, standard deviations, rates of change, degree days) is described, as are the methods used to estimate missing daily mean temperatures, degree days (annual and summer periods), and 7-day running averages of daily mean temperatures.

Data Series↗

Williston Basin groundwater availability, United States and Canada

The Williston Basin contains important oil and gas resources for the Nation. Freshwater supplies are limited in this semiarid area, and oil and gas development can require large volumes of freshwater. Groundwater is the primary source of water for many water users in the Williston Basin, so to better understand these resources, the U.S. Geological Survey (USGS) assessed the groundwater availability in this area. The final phase of this assessment included a computer model that simulates how groundwater flows in the aquifer systems and simulates how changes in water use and natural conditions may affect the water resources. These results provide a tool for land and water-resource managers to determine how water can be used for multiple purposes in the Williston Basin. For additional information about this assessment and more in-depth descriptions and results, see Long and others (2018).

Williston Basin↗

A compendium of P- and S-wave velocities from surface-to-borehole logging; summary and reanalysis of previously published data and analysis of unpublished data

For over 28 years, the U.S. Geological Survey (USGS) has been acquiring seismic velocity and geologic data at a number of locations in California, many of which were chosen because strong ground motions from earthquakes were recorded at the sites. The method for all measurements involves picking first arrivals of P- and S-waves from a surface source recorded at various depths in a borehole (as opposed to noninvasive methods, such as the SASW method [e.g., Brown et al., 2002]). The results from most of the sites are contained in a series of U.S. Geological Survey Open-File Reports (see References). Until now, none of the results have been available as computer files, and before 1992 the interpretation of the arrival times was in terms of piecemeal interval velocities, with no attempt to derive a layered model that would fit the travel times in an overall sense (the one exception is Porcella, 1984). In this report I reanalyze all of the arrival times in terms of layered models for P- and for S-wave velocities at each site, and I provide the results as computer files. In addition to the measurements reported in the open-file reports, I also include some borehole results from other reports, as well as some results never before published. I include data for 277 boreholes (at the time of this writing; more will be added to the web site as they are obtained), all in California (I have data from boreholes in Washington and Utah, but these will be published separately). I am also in the process of interpreting travel time data obtained using a seismic cone penetrometer at hundreds of sites; these data can be interpreted in the same way of those obtained from surface-to-borehole logging. When available, the data will be added to the web site (see below for information on obtaining data from the World Wide Web (WWW)). In addition to the basic borehole data and results, I provide information concerning strong-motion stations that I judge to be close enough to the boreholes that the borehole velocity models can be used as the velocity models beneath the stations.

California↗

Volcanogenic Massive Sulfide Deposits of the World - Database and Grade and Tonnage Models

Grade and tonnage models are useful in quantitative mineral-resource assessments. The models and database presented in this report are an update of earlier publications about volcanogenic massive sulfide (VMS) deposits. These VMS deposits include what were formerly classified as kuroko, Cyprus, and Besshi deposits. The update was necessary because of new information about some deposits, changes in information in some deposits, such as grades, tonnages, or ages, revised locations of some deposits, and reclassification of subtypes. In this report we have added new VMS deposits and removed a few incorrectly classified deposits. This global compilation of VMS deposits contains 1,090 deposits; however, it was not our intent to include every known deposit in the world. The data was recently used for mineral-deposit density models (Mosier and others, 2007; Singer, 2008). In this paper, 867 deposits were used to construct revised grade and tonnage models. Our new models are based on a reclassification of deposits based on host lithologies: Felsic, Bimodal-Mafic, and Mafic volcanogenic massive sulfide deposits. Mineral-deposit models are important in exploration planning and quantitative resource assessments for two reasons: (1) grades and tonnages among deposit types vary significantly, and (2) deposits of different types occur in distinct geologic settings that can be identified from geologic maps. Mineral-deposit models combine the diverse geoscience information on geology, mineral occurrences, geophysics, and geochemistry used in resource assessments and mineral exploration. Globally based deposit models allow recognition of important features and demonstrate how common different features are. Well-designed deposit models allow geologists to deduce possible mineral-deposit types in a given geologic environment and economists to determine the possible economic viability of these resources. Thus, mineral-deposit models play a central role in presenting geoscience information in a useful form to policy makers. The foundation of mineral-deposit models is information about known deposits. The purpose of this publication is to present the latest geologic information and newly developed grade and tonnage models for VMS deposits in digital form. This publication contains computer files with information on VMS deposits from around the world. It also presents new grade and tonnage models for three subtypes of VMS deposits and a text file allowing locations of all deposits to be plotted in geographic information system (GIS) programs. The data are presented in FileMaker Pro and text files to make the information available to a wider audience. The value of this information and any derived analyses depends critically on the consistent manner of data gathering. For this reason, we first discuss the rules used in this compilation. Next, we provide new grade and tonnage models and analysis of the information in the file. Finally, the fields of the data file are explained. Appendix A gives the summary statistics for the new grade-tonnage models and Appendix B displays the country codes used in the database.

Open-File Report↗

Geographic information system programs for use in the water-supply-allocation permitting process

Computer programs designed for use in a geographic information system as an aid in the water-supply- allocation permitting process are described. These programs were developed by the U.S. Geological Survey during a project conducted in cooperation with the New Jersey Department of Environmental Protection. The programs enable a user to display proposed water-supply-allocation sites in a defined area together with present sites and important hydrologic and geographic features on a computer screen or on hardcopy plots. The programs are menu-driven and do not require familiarity with geographic information systems. Source codes for the programs are included in appendixes.

Open-File Report↗

Comparison of peak-flow estimation methods for small drainage basins in Maine

Understanding the accuracy of commonly used methods for estimating peak streamflows is important because the designs of bridges, culverts, and other river structures are based on these flows. Different methods for estimating peak streamflows were analyzed for small drainage basins in Maine. For the smallest basins, with drainage areas of 0.2 to 1.0 square mile, nine peak streamflows from actual rainfall events at four crest-stage gaging stations were modeled by the Rational Method and the Natural Resource Conservation Service TR-20 method and compared to observed peak flows. The Rational Method had a root mean square error (RMSE) of -69.7 to 230 percent (which means that approximately two thirds of the modeled flows were within -69.7 to 230 percent of the observed flows). The TR-20 method had an RMSE of -98.0 to 5,010 percent. Both the Rational Method and TR-20 underestimated the observed flows in most cases. For small basins, with drainage areas of 1.0 to 10 square miles, modeled peak flows were compared to observed statistical peak flows with return periods of 2, 50, and 100 years for 17 streams in Maine and adjoining parts of New Hampshire. Peak flows were modeled by the Rational Method, the Natural Resources Conservation Service TR-20 method, U.S. Geological Survey regression equations, and the Probabilistic Rational Method. The regression equations were the most accurate method of computing peak flows in Maine for streams with drainage areas of 1.0 to 10 square miles with an RMSE of -34.3 to 52.2 percent for 50-year peak flows. The Probabilistic Rational Method was the next most accurate method (-38.5 to 62.6 percent). The Rational Method (-56.1 to 128 percent) and particularly the TR-20 method (-76.4 to 323 percent) had much larger errors. Both the TR-20 and regression methods had similar numbers of underpredictions and overpredictions. The Rational Method overpredicted most peak flows and the Probabilistic Rational Method tended to overpredict peak flows from the smaller (less than 5 square miles) drainage basins and underpredict peak flows from larger drainage basins. The results of this study are consistent with the most comprehensive analysis of observed and modeled peak streamflows in the United States, which analyzed statistical peak flows from 70 drainage basins in the Midwest and the Northwest.

Maine↗

PRMS-IV, the precipitation-runoff modeling system, version 4

Computer models that simulate the hydrologic cycle at a watershed scale facilitate assessment of variability in climate, biota, geology, and human activities on water availability and flow. This report describes an updated version of the Precipitation-Runoff Modeling System. The Precipitation-Runoff Modeling System is a deterministic, distributed-parameter, physical-process-based modeling system developed to evaluate the response of various combinations of climate and land use on streamflow and general watershed hydrology. Several new model components were developed, and all existing components were updated, to enhance performance and supportability. This report describes the history, application, concepts, organization, and mathematical formulation of the Precipitation-Runoff Modeling System and its model components. This updated version provides improvements in (1) system flexibility for integrated science, (2) verification of conservation of water during simulation, (3) methods for spatial distribution of climate boundary conditions, and (4) methods for simulation of soil-water flow and storage.

Techniques and Methods↗

Development and Application of Watershed Regressions for Pesticides (WARP) for Estimating Atrazine Concentration Distributions in Streams

Regression models were developed for predicting atrazine concentration distributions in rivers and streams, using the Watershed Regressions for Pesticides (WARP) methodology. Separate regression equations were derived for each of nine percentiles of the annual distribution of atrazine concentrations and for the annual time-weighted mean atrazine concentration. In addition, seasonal models were developed for two specific periods of the year--the high season, when the highest atrazine concentrations are expected in streams, and the low season, when concentrations are expected to be low or undetectable. Various nationally available watershed parameters were used as explanatory variables, including atrazine use intensity, soil characteristics, hydrologic parameters, climate and weather variables, land use, and agricultural management practices. Concentration data from 112 river and stream stations sampled as part of the U.S. Geological Survey's National Water-Quality Assessment and National Stream Quality Accounting Network Programs were used for computing the concentration percentiles and mean concentrations used as the response variables in regression models. Tobit regression methods, using maximum likelihood estimation, were used for developing the models because some of the concentration values used for the response variables were censored (reported as less than a detection threshold). Data from 26 stations not used for model development were used for model validation. The annual models accounted for 62 to 77 percent of the variability in concentrations among the 112 model development stations. Atrazine use intensity (the amount of atrazine used in the watershed divided by watershed area) was the most important explanatory variable in all models, but additional watershed parameters significantly increased the amount of variability explained by the models. Predicted concentrations from all 10 models were within a factor of 10 of the observed concentrations at most model development and model validation stations. Results for the two sets of seasonal models were similar. Concentration distributions derived from the seasonal-model predictions provided additional information compared to distributions derived from the annual models.

Water-Resources Investigations Report↗

DIGITAL GEOLOGIC MAP OF THE UNITED STATES.

The geologic map of the United States was published in 1974 by the U. S. Geological Survey. This major publication contains an enormous amount of information on the surficial geology of the United States. Many geologists have used this map as a research tool. Most have needed information from only specific parts of the map, and have manually extracted data from these areas. These data have then been combined with other geological information, much of which - especially that concerning minerals and energy - either is already in machine-readable computer files or is rapidly being converted to that form.

Conference Paper↗

Mineral resource of the month: iron oxide pigments

The article discusses iron oxide pigments, which have been used as colorants since human began painting as they resist color change due to sunlight exposure, have good chemical resistance and are stable under normal ambient conditions. Cyprus, Italy and Spain are among the countries that are known for the production of iron oxide pigments. Granular forms of iron oxides and nano-sized materials are cited as developments in the synthetic iron oxide pigment industry which are being used in computer disk drives and nuclear magnetic resonance imaging.

Earth↗

Two-dimensional steady-state dispersion in a saturated porous medium

A previously developed analytical solution for two-dimensional dispersion is computed for various conditions. These results were then compared with solution of previously developed approximate models of transverse dispersion which were used to analyze experimentally derived concentration distribution. Comparison established that, whenever steady state was reached, the values of dispersion coefficient computed using the approximate expression agreed with values derived from exact expression and in addition allowed a quick computation of the parameter.

Journal of Research of the U.S. Geological Survey↗

Effective site coefficients for the 2024 International Building Code (IBC)

The U.S. National Seismic Hazard Models (NSHMs), developed by the U.S. Geological Survey (USGS), have long been the scientific foundation of seismic design guidelines and have been used to compute design ground motions for construction of new buildings and retrofit of existing buildings in the United States and its territories. The 2018 NSHM is adopted by the 2024 International Building Code (IBC). Prior to the 2018 NSHM update, hazard calculations were required only at one reference site condition defined by V S30 =760 m/s (where V S30 is the time-averaged shear wave velocity from the surface to a depth of 30 m) and three periods (peak ground acceleration, PGA, and pseudo spectral accelerations at periods of 0.2 and 1 s, S s and S 1 ). Site coefficients, F PGA , F a , and F v , were then defined by the Building Seismic Safety Council (BSSC) Provisions Update Committee (PUC) in the site-specific procedures of National Earthquake Hazard Reduction Program (NEHRP) Recommended Seismic Provisions to calculate ground motions for other site classes with different V S30 values at the given periods. The design ground motions at other periods were then estimated using a generic spectral shape that was also defined by the BSSC PUC in NEHRP provisions. In recent years, the engineering community has realized there were deficiencies with the F a and F v site coefficients and the generic spectral shape. To avoid potentially dangerous underestimations of design ground motions for long period structures on soft site conditions in high seismicity regions, the BSSC PUC recommended the use of multi-period response spectra (MPRS) in 2017. As a result, the USGS produced multi-period and multi-V S30 response spectral values in the 2018 NSHM for calculations of design ground motions and the site coefficients F a and F v were eliminated from the 2020 NEHRP Provisions. As these site coefficients were widely used inside and outside of the United States, in this study we back-calculate the “effective” site coefficients F a,eff , and F v,eff by comparing MPRS for various site classes with the MPRS for the reference site condition, and discuss the changes that are observed in the 2024 IBC compared to its previous version in 2021. The effective site coefficients are presented for test site locations and their dependence on various factors including period, ground motion intensity, and regional models are discussed. Ratio maps between the new effective site coefficients and the old ones are then presented for soft site classes and for short and long periods. For soft site classes at short periods, the new effective site coefficients are lower than the old site coefficients for high seismicity regions and higher for low seismicity regions. As it was expected, for soft site classes at long periods and high seismicity regions, the new effective site coefficients are much larger than the old site coefficients without imposing the 50% increase as a penalty that was suggested in the 2021 IBC, whereas they could be much smaller if the 50% increase would have been imposed particularly around New Madrid and Charleston high seismicity regions. For low seismicity regions, the long period effective site coefficients can be smaller or larger by 20% compared to the 2021 IBC coefficients.

Conference Paper↗

Groundwater and surface-water interaction, water quality, and processes affecting loads of dissolved solids, selenium, and uranium in Fountain Creek, near Pueblo, Colorado, 2012–2014

In 2012, the U.S. Geological Survey, in cooperation with the Arkansas River Basin Regional Resource Planning Group, initiated a study of groundwater and surface-water interaction, water quality, and loading of dissolved solids, selenium, and uranium to Fountain Creek near Pueblo, Colorado, to improve understanding of sources and processes affecting loading of these constituents to streams in the Arkansas River Basin. Fourteen monitoring wells were installed in a series of three transects across Fountain Creek near Pueblo, and temporary streamgages were established at each transect to facilitate data collection for the study. Groundwater and surface-water interaction was characterized by using hydrogeologic mapping, groundwater and stream-surface levels, groundwater and stream temperatures, vertical hydraulic-head gradients and ratios of oxygen and hydrogen isotopes in the hyporheic zone, and streamflow mass-balance measurements. Water quality was characterized by collecting periodic samples from groundwater, surface water, and the hyporheic zone for analysis of dissolved solids, selenium, uranium, and other selected constituents and by evaluating the oxidation-reduction condition for each groundwater sample under different hydrologic conditions throughout the study period. Groundwater loads to Fountain Creek and in-stream loads were computed for the study area, and processes affecting loads of dissolved solids, selenium, and uranium were evaluated on the basis of geology, geochemical conditions, land and water use, and evapoconcentration. During the study period, the groundwater-flow system generally contributed flow to Fountain Creek and its hyporheic zone (as a single system) except for the reach between the north and middle transects. However, the direction of flow between the stream, the hyporheic zone, and the near-stream aquifer was variable in response to streamflow and stage. During periods of low streamflow, Fountain Creek generally gained flow from groundwater. However, during periods of high streamflow, the hydraulic gradient between groundwater and the stream temporarily reversed, causing the stream to lose flow to groundwater. Concentrations of dissolved solids, selenium, and uranium in groundwater generally had greater spatial variability than surface water or hyporheic-zone samples, and constituent concentrations in groundwater generally were greater than in surface water. Constituent concentrations in the hyporheic zone typically were similar to or intermediate between concentrations in groundwater and surface water. Concentrations of dissolved solids, selenium, uranium, and other constituents in groundwater samples collected from wells located on the east side of the north monitoring well transect were substantially greater than for other groundwater, surface-water, and hyporheic-zone samples. With one exception, groundwater samples collected from wells on the east side of the north transect exhibited oxic to mixed (oxic-anoxic) conditions, whereas most other groundwater samples exhibited anoxic to suboxic conditions. Concentrations of dissolved solids, selenium, and uranium in surface water generally increased in a downstream direction along Fountain Creek from the north transect to the south transect and exhibited an inverse relation to streamflow with highest concentration occurring during periods of low streamflow and lowest concentrations occurring during periods of high streamflow. Groundwater loads of dissolved solids, selenium, and uranium to Fountain Creek were small because of the small amount of groundwater flowing to the stream under typical low-streamflow conditions. In-stream loads of dissolved solids, selenium, and uranium in Fountain Creek varied by date, primarily in relation to streamflow at each transect and were much larger than computed constituent loads from groundwater. In-stream loads generally decreased with decreases in streamflow and increased as streamflow increased. In-stream loads of dissolved solids and selenium increased between the north and middle transects but generally decreased between the middle and south transects. By contrast, uranium loads generally decreased between the north and middle transects but increased between the middle and south transects. In-stream load differences between transects appear primarily to be related to differences in streamflow. However, because groundwater typically flows to Fountain Creek under low-flow conditions, and groundwater has greater concentrations of dissolved solids, selenium, and uranium than surface water in Fountain Creek, increases in loads between transects likely are affected by inflow of groundwater to the stream, which can account for a substantial proportion of the in-stream load difference between transects. When loads decreased between transects, the primary cause likely was decreased streamflow as a result of losses to groundwater and flow through the hyporheic zone. However, localized groundwater inflow likely attenuated the magnitude by which the in-stream loads decreased. The combination of localized soluble geologic sources and oxic conditions likely is the primary reason for the occurrence of high concentrations of dissolved solids, selenium, and uranium in groundwater on the east side of the north monitoring well transect. To evaluate conditions potentially responsible for differences in water quality and redox conditions, physical characteristics such as depth to water, saturated thickness, screen depth below the water table, screen height above bedrock, and aquifer hydraulic conductivity were compared by using Wilcoxon rank-sum tests. Results indicated no significant difference between depth to water, screen height above bedrock, and hydraulic conductivity for groundwater samples collected from wells on the east side of the north transect and groundwater samples from all other wells. However, saturated thickness and screen depth below the water table both were significantly smaller for groundwater samples collected from wells on the east side of the north transect than for groundwater samples from other wells, indicating that these characteristics might be related to the elevated constituent concentrations found at that location. Similarly, saturated thickness and screen depth below the water table were significantly smaller for groundwater samples under oxic or mixed (oxic-anoxic) conditions than for those under anoxic to suboxic conditions. The greater constituent concentrations at wells on the east side of the north transect also could, in part, be related to groundwater discharge from an unnamed alluvial drainage located directly upgradient from that location. Although the quantity and quality of water discharging from the drainage is not known, the drainage appears to collect water from a residential area located upgradient to the east of the wells, and groundwater could become concentrated in nitrate and other dissolved constituents before flowing through the drainage. High levels of nitrate, whether from anthropogenic or natural geologic sources, could promote more soluble forms of selenium and other constituents by affecting the redox condition of groundwater. Whether oxic conditions at wells on the east side of the north transect are the result of physical characteristics or of groundwater inflow from the alluvial drainage, the oxic conditions appear to cause increased dissolution of minerals from the shallow shale bedrock at that location. Because ratios of hydrogen and oxygen isotopes indicate evaporation likely has not had a substantial effect on groundwater, constituent concentrations at that location likely are not the result of evapoconcentration.

Colorado↗

Updated computations and estimates of streamflows tributary to Carson Valley, Douglas County, Nevada, and Alpine County, California, 1990-2002

Rapid population growth in Carson Valley has caused concern over the continued availability of water resources to sustain future growth. The U.S. Geological Survey, in cooperation with Douglas County, began a study to update estimates of water-budget components in Carson Valley for current climatic conditions. Data collected at 19 sites included 9 continuous records of tributary streamflows, 1 continuous record of outflow from the valley, and 408 measurements of 10 perennially flowing but ungaged drainages. These data were compiled and analyzed to provide updated computations and estimates of streamflows tributary to Carson Valley, 1990-2002. Mean monthly and annual flows were computed from continuous records for the period 1990-2002 for five streams, and for the period available, 1990-97, for four streams. Daily mean flow from ungaged drainages was estimated using multi-variate regressions of individual discharge measurements against measured flow at selected continuous gages. From the estimated daily mean flows, monthly and annual mean flows were calculated from 1990 to 2002. These values were used to compute estimates of mean monthly and annual flows for the ungaged perennial drainages. Using the computed and estimated mean annual flows, annual unit-area runoff was computed for the perennial drainages, which ranged from 0.30 to 2.02 feet. For the period 1990-2002, estimated inflow of perennial streams tributary to Carson Valley totaled about 25,900 acre-feet per year. Inflow computed from gaged perennial drainages totaled 10,300 acre-feet per year, and estimated inflow from ungaged perennial drainages totaled 15,600 acre-feet per year. The annual flow of perennial streams ranges from 4,210 acre-feet at Clear Creek to 450 acre-feet at Stutler Canyon Creek. Differences in unit-area runoff and in the seasonal timing of flow likely are caused by differences in geologic setting, altitude, slope, or aspect of the individual drainages. The remaining drainages are ephemeral and supply inflow to the valley floor only during spring runoff in wet years or during large precipitation events. Annual unit-area runoff for the perennial drainages was used to estimate inflow from ephemeral drainages totaling 11,700 acre-feet per year. The totaled estimate of perennial and ephemeral tributary inflows to Carson Valley is 37,600 acre-feet per year. Gaged perennial inflow is 27 percent of the total, ungaged perennial inflow is 42 percent, and ephemeral inflow is 31 percent. The estimate is from 50 to 60 percent greater than three previous estimates, one made for a larger area and similar to two other estimates made for larger areas. The combined uncertainty of the estimates totaled about 33 percent of the total inflow or about 12,000 acre-feet per year.

Scientific Investigations Report↗

FEMFLOW3D; a finite-element program for the simulation of three-dimensional aquifers; version 1.0

This document describes a computer program that simulates three-dimensional ground- water systems using the finite-element method. The program was developed to simulate regional ground-water systems, but it can be applied to small-scale problems as well. This program can be used to simulate both confined and water-table aquifers. The program simulates a linearized three-dimensional free-surface ground-water system with a fixed grid. FEMFLOW3D is applicable to the simulation of various free-surface ground-water systems for which the change in aquifer thickness is small relative to the overall aquifer thickness. The finite-element method provides flexibility in the design of a geometric grid that represents the physical dimensions of an aquifer system. For example, features that can be well represented with a finite-element grid include irregular, random geographic and geologic features; irregular boundaries; and increased detail within localized areas of particular interest within the study area. The structure of the computer program consists of a main program, which serves as a simple driver, and a set of subroutines in which the calculations are performed. The background, mathematical basis, structure, and inputs for each of the subroutines are described in the document where applicable. Each subroutine generally handles (1) a part of the mathematical calculations related to the finite-element method, (2) a specific feature of the hydrologic system, or (3) special features related to the management input or output data. Hydrologic features that can be represented with the program include stream-aquifer interactions, phreatophytic evapotranspiration, highly permeable fault zones, land subsidence, and land-aquifer interactions associated with land-use activities. The program can also represent the primary features associated with complex irrigation systems, such as irrigated agriculture, and can calculate the ground-water recharge that results from these activities. Three boundary conditions, including specified-head boundaries, specified-flux boundaries, and variable-flux boundaries, can be represented with the program. The program also provides a method for identifying aquifer and river-bed parameters that can be used in the calibration of models. This document also includes model validation, source code, and example input and output files. Model validation was performed using four test problems. For each test problem, the results of a model simulation with FEMFLOW3D were compared with either an analytic solution or the results of an independent numerical approach. The source code, written in the ANSI x3.9-1978 FORTRAN standard, and the complete input and output of an example problem are listed in the appendixes.

Open-File Report↗

Oil spills, 1971-75, Gulf of Mexico outer continental shelf

Oil spillage connected with federally supervised drilling and production activities has been a matter of wide public concern. In its supervision of mineral-resource development on the Outer Continental Shelf (OCS), the U.S. Geological Survey is responsible for the day-to-day inspection and monitoring of OCS oil and gas operations. During these activities, the U.S. Geological Survey records and investigates hydrocarbon discharges resulting from such operations. Beginning in 1971, all spills have been recorded, and a computer file has been maintained on all spills of 1 barrel or more. The total Gulf of Mexico OCS oil spillage recorded during January 1, 1971-December 31, 1975, amounted to 51,421 barrels. Production during that period amounted to 35,219 barrels per barrel spilled. In all, 5,857 spills were recorded, but 85.5 percent of the total spill volume was contributed by just 5 incidents. The environmental effect of these incidents apparently was minimal and of short duration. No spills of more than 50 barrels resulted from drilling operations during the period. The only spillage resulting from blowouts was caused by nondrilling incidents, including completion, production, and workover. The amount of oil discharged from spills of less than 50 barrels decreased by more than half between 1971 and 1975. The improvement reflects changes in the operating philosophy of the offshore industry, tightening of U.S. Geological Survey operating orders, and substantial increases in the inspection force. Most production-platform spills involve failures in the sump system, the separator system, or other hydrocarbon-handling equipment; improved sump-system designs and better high-low-level controls have reduced both the number and the volume of spills. Pipeline and pump spills also declined significantly, although the decline appears less attributable to revisions in OCS operating requirements. No operator consistently contributed a disproportionate amount of spillage. Most of the platforms having production spills had only one spill during the period. Of the 304 platforms having incidents, 255 had 2 spills or less and 188 had only 1 spill. Comparison of the volume spilled with the volume produced for Gulf of Mexico OCS operations during 1971-75 results in a spillage rate of 0.0028 percent. Spillage from tankers is estimated to be greater. The OCS drilling and production spillage rate is not expected to increase during exploration and development of the frontier areas.

Circular↗

Scaling variables and interpretation of eigenvalues in principal component analysis of geologic data

The dominant feature distinguishing one method of principal components analysis from another is the manner in which the original data are transformed prior to the other computations. The only other distinguishing feature of any importance is whether the eigenvectors of the inner product-moment of the transformed data matrix are taken directly as the Q-mode scores or scaled by the square roots of their associated eigenvalues and called the R-mode loadings. If the eigenvectors are extracted from the product-moment correlation matrix, the variables, in effect, were transformed by column standardization (zero means and unit variances), and the sum of the p-largest eigenvalues divided by the sum of all the eigenvalues indicates the degree to which a model containing p components will account for the total variance in the original data. However, if the data were transformed in any manner other than column standardization, the eigenvalues cannot be used in this manner, but can only be used to determine the degree to which the model will account for the transformed data. Regardless of the type of principal components analysis that is performed-even whether it is R or Q-mode-the goodness-of-fit of the model to the original data is given better by the eigenvalues of the correlation matrix than by those of the matrix that was actually factored. ?? 1980 Plenum Publishing Corporation.

Journal of the International Association for Mathe↗

Comparison of methods for computing streamflow statistics for Pennsylvania streams

Methods for computing streamflow statistics intended for use on ungaged locations on Pennsylvania streams are presented and compared to frequency distributions of gaged streamflow data. The streamflow statistics used in the comparisons include the 7-day 10-year low flow, 50-year flood flow, and the 100-year flood flow; additional statistics are presented. Streamflow statistics for gaged locations on streams in Pennsylvania were computed using three methods for the comparisons: 1) Log-Pearson type III frequency distribution (Log-Pearson) of continuous-record streamflow data, 2) regional regression equations developed by the U.S. Geological Survey in 1982 (WRI 82-21), and 3) regional regression equations developed by the Pennsylvania State University in 1981 (PSU-IV). Log-Pearson distribution was considered the reference method for evaluation of the regional regression equations. Low-flow statistics were computed using the Log-Pearson distribution and WRI 82-21, whereas flood-flow statistics were computed using all three methods. The urban adjustment for PSU-IV was modified from the recommended computation to exclude Philadelphia and the surrounding areas (region 1) from the adjustment. Adjustments for storage area for PSU-IV were also slightly modified. A comparison of the 7-day 10-year low flow computed from Log-Pearson distribution and WRI-82- 21 showed that the methods produced significantly different values for about 7 percent of the state. The same methods produced 50-year and 100-year flood flows that were significantly different for about 24 percent of the state. Flood-flow statistics computed using Log-Pearson distribution and PSU-IV were not significantly different in any regions of the state. These findings are based on a statistical comparison using the t-test on signed ranks and graphical methods.

Water-Resources Investigations Report↗