Geology Reports⌕ Search

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.

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

At least 1,135 records · Page 63Linked to original sources

A filter circuit board for the Earthworm Seismic Data Acquisition System

The Earthworm system is a seismic network data acquisition and processing system used by the Northern California Seismic Network as well as many other seismic networks. The input to the system is comprised of many realtime electronic waveforms fed to a multi-channel digitizer on a PC platform. The digitizer consists of one or more National Instruments Corp. AMUX–64T multiplexer boards attached to an A/D converter board located in the computer. Originally, passive filters were installed on the multiplexers to eliminate electronic noise picked up in cabling. It was later discovered that a small amount of crosstalk occurred between successive channels in the digitizing sequence. Though small, this crosstalk will cause what appear to be small earthquake arrivals at the wrong time on some channels. This can result in erroneous calculation of earthquake arrival times, particularly by automated algorithms. To deal with this problem, an Earthworm filter board was developed to provide the needed filtering while eliminating crosstalk. This report describes the tests performed to find a suitable solution, and the design of the circuit board. Also included are all the details needed to build and install this board in an Earthworm system or any other system using the AMUX–64T board. Available below is the report in PDF format as well as an archive file containing the circuit board manufacturing information.

Open-File Report↗

Estimation of hydraulic parameters from an unconfined aquifer test conducted in a glacial outwash deposit, Cape Cod, Massachusetts

An aquifer test conducted in a sand and gravel, glacial outwash deposit on Cape Cod, Massachusetts was analyzed by means of a model for flow to a partially penetrating well in a homogeneous, anisotropic unconfined aquifer. The model is designed to account for all significant mechanisms expected to influence drawdown in observation piezometers and in the pumped well. In addition to the usual fluid-flow and storage processes, additional processes include effects of storage in the pumped well, storage in observation piezometers, effects of skin at the pumped-well screen, and effects of drainage from the zone above the water table. The aquifer was pumped at a rate of 320 gallons per minute for 72-hours and drawdown measurements were made in the pumped well and in 20 piezometers located at various distances from the pumped well and depths below the land surface. To facilitate the analysis, an automatic parameter estimation algorithm was used to obtain relevant unconfined aquifer parameters, including the saturated thickness and a set of empirical parameters that relate to gradual drainage from the unsaturated zone. Drainage from the unsaturated zone is treated in this paper as a finite series of exponential terms, each of which contains one empirical parameter that is to be determined. It was necessary to account for effects of gradual drainage from the unsaturated zone to obtain satisfactory agreement between measured and simulated drawdown, particularly in piezometers located near the water table. The commonly used assumption of instantaneous drainage from the unsaturated zone gives rise to large discrepancies between measured and predicted drawdown in the intermediate-time range and can result in inaccurate estimates of aquifer parameters when automatic parameter estimation procedures are used. The values of the estimated hydraulic parameters are consistent with estimates from prior studies and from what is known about the aquifer at the site. Effects of heterogeneity at the site were small as measured drawdowns in all piezometers and wells were very close to the simulated values for a homogeneous porous medium. The estimated values are: specific yield, 0.26; saturated thickness, 170 feet; horizontal hydraulic conductivity, 0.23 feet per minute; vertical hydraulic conductivity, 0.14 feet per minute; and specific storage, 1.3x10-5 per foot. It was found that drawdown in only a few piezometers strategically located at depth near the pumped well yielded parameter estimates close to the estimates obtained for the entire data set analyzed simultaneously. If the influence of gradual drainage from the unsaturated zone is not taken into account, specific yield is significantly underestimated even in these deep-seated piezometers. This helps to explain the low values of specific yield often reported for granular aquifers in the literature. If either the entire data set or only the drawdown in selected deep-seated piezometers was used, it was found unnecessary to conduct the test for the full 72-hours to obtain accurate estimates of the hydraulic parameters. For some piezometer groups, practically identical results would be obtained for an aquifer test conducted for only 8-hours. Drawdowns measured in the pumped well and piezometers at distant locations were diagnostic only of aquifer transmissivity.

Massachusetts↗

A parallel-processing approach to computing for the geographic sciences

The overarching goal of this project is to build a spatially distributed infrastructure for information science research by forming a team of information science researchers and providing them with similar hardware and software tools to perform collaborative research. Four geographically distributed Centers of the U.S. Geological Survey (USGS) are developing their own clusters of low-cost personal computers into parallel computing environments that provide a costeffective way for the USGS to increase participation in the high-performance computing community. Referred to as Beowulf clusters, these hybrid systems provide the robust computing power required for conducting research into various areas, such as advanced computer architecture, algorithms to meet the processing needs for real-time image and data processing, the creation of custom datasets from seamless source data, rapid turn-around of products for emergency response, and support for computationally intense spatial and temporal modeling.

Open-File Report↗

Digital recovery, modification, and analysis of Tetra Tech seismic horizon mapping, National Petroleum Reserve Alaska (NPRA), northern Alaska

We have digitized, modified, and analyzed seismic interpretation maps of 12 subsurface stratigraphic horizons spanning portions of the National Petroleum Reserve in Alaska (NPRA). These original maps were prepared by Tetra Tech, Inc., based on about 15,000 miles of seismic data collected from 1974 to 1981. We have also digitized interpreted faults and seismic velocities from Tetra Tech maps. The seismic surfaces were digitized as two-way travel time horizons and converted to depth using Tetra Tech seismic velocities. The depth surfaces were then modified by long-wavelength corrections based on recent USGS seismic re-interpretation along regional seismic lines. We have developed and executed an algorithm to identify and calculate statistics on the area, volume, height, and depth of closed structures based on these seismic horizons. These closure statistics are tabulated and have been used as input to oil and gas assessment calculations for the region. Directories accompanying this report contain basic digitized data, processed data, maps, tabulations of closure statistics, and software relating to this project.

Open-File Report↗

Analysis of magnetotelluric profile data from the Ruby Mountains metamorphic core complex and southern Carlin Trend region, Nevada

We have collected about 150 magnetotelluric (MT) soundings in northeastern Nevada in the region of the Ruby Mountains metamorphic core complex uplift and southern Carlin mineral trend, in an effort to illuminate controls on core complex evolution and deposition of world-class gold deposits. The region has experienced a broad range of tectonic events including several periods of compressional and extensional deformation, which have contributed to the total expression of electrical resistivity. Most of the soundings are in three east-west profiles across increasing degrees of core uplift to the north (Bald Mountain, Harrison Pass and Secret Pass latitudes). Two shorter lines cross a prominent east-west structure to the north of the northern profile. MT impedance tensor and vertical magnetic field rotations imply a N-NNE average regional geoelectric strike, similar to surface geologic trends. Model resistivity cross sections were derived using a 2-D inversion algorithm, which damps departures of model parameters from an a priori structure, emphasizing the transverse magnetic (TM) mode and vertical magnetic field data. Geological interpretation of the resistivity combines previous seismic, potential field and isotope models, structural and petrological models for regional compression and extension, and detailed structural/stratigraphic interpretations incorporating drilling for petroleum and mineral exploration. To first order, the resistivity structure is one of a moderately conductive, Phanerozoic sedimentary section fundamentally disrupted by intrusion and uplift of resistive crystalline rocks. Late Devonian and early Mississippian shales of the Pilot and Chainman Formations together form an important conductive marker sequence in the stratigraphy and show pronounced increases in conductance (conductivity-thickness product) from east to west. These increases in conductance are attributed to graphitization caused by Elko-Sevier era compressional shear deformation and possibly by intrusive heating. The resistive crystalline central massifs adjoin the host stratigraphy across crustal-scale, subvertical fault zones. These zones provide electric current pathways to the lower crust for heterogeneous, upper crustal induced current flow. Resistive core complex crust may be steeply bounded under the middle of the neighboring grabens and not deepen at a shallow angle to arbitrary distances to the west. The numerous crustal breaks imaged with MT may contribute to the low effective elastic thickness estimated regionally for the Great Basin and exemplify the mid-crustal, steeply dipping slip zones in which major earthquakes nucleate. An east-west oriented conductor in the crystalline upper crust spans the East Humboldt Range and northern Ruby Mountains. The conductor may be related to an inferred ArcheanProterozoic suture or nearby graphitic metasediments, with possible alteration by middle Tertiary magmatic activity. Lower crustal resistivity everywhere under the profiles is low and appears quasi one-dimensional. It is consistent with a low rock porosity (

Open-File Report↗

Regional sediment budget of the Columbia River littoral cell, USA: Analysis of bathymetric- and topographic-volume change

In this Open-File Report we present calculations of changes in bathymetric and topographic volumes for the Grays Harbor, Willapa Bay, and Columbia River entrances and the adjacent coasts of North Beach, Grayland Plains, Long Beach, and Clatsop Plains for four intervals: pre-jetty - 1920s (Interval 1), 1920s - 1950s (Interval 2), 1950s - 1990s (Interval 3), and 1920s 1990s (Interval 4). This analysis is part of the Southwest Washington Coastal Erosion Study (SWCES), the goals of which are to understand and predict the morphologic behavior of the Columbia River littoral cell on a management scale of tens of kilometers and decades. We obtain topographic Light Detection and Ranging (LIDAR) data from a joint project by the U.S. Geological Survey (USGS), National Oceanic and Atmospheric Administration (NOAA), National Aeronautic and Space Administration (NASA), and the Washington State Department of Ecology (DOE) and bathymetric data from the U.S. Coast and Geodetic Survey (USC&GS), U.S. Army Corps of Engineers (USACE), USGS, and the DOE. Shoreline data are digitized from T-Sheets and aerial photographs from the USC&GS and National Ocean Service (NOS). Instead of uncritically adjusting each survey to NAVD88, a common vertical land-based datum, we adjust some surveys to produce optimal results according to the following criteria. First, we minimize offsets in overlapping surveys within the same era, and second, we minimize bathymetric changes (relative to the 1990s) in deep water, where we assume minimal change has taken place. We grid bathymetric and topographic datasets using kriging and triangulation algorithms, calculate bathymetric-change surfaces for each interval, and calculate volume changes within polygons that are overlaid on the bathymetric-change surfaces. We find similar morphologic changes near the entrances to Grays Harbor and the Columbia River following jetty construction between 1898 and 1916 at the Grays Harbor entrance and between 1885 and 1913 at the Columbia River entrance. The inlets and inner deltas eroded and the outer deltas moved offshore and accreted. The adjacent coasts experienced accretion over alongshore distances of tens of kilometers. North of the Grays Harbor entrance along North Beach and north of the Columbia River entrance along Long Beach the shoreface and the beach-dune complex mainly prograded, whereas south of the Grays Harbor entrance along Grayland Plains and south of the Columbia River entrance along Clatsop Plains the beach-dune complex above -10 m NAVD88 prograded and the shoreface between approximately -30 m and -10 m NAVD88 eroded. In the decades following jetty construction, the rates of erosion and accretion at the entrances decreased and the centers of deposition along the adjacent coasts moved away from the entrances. The rates of change have decreased, suggesting the systems are approaching dynamic equilibrium. Exceptions to this behaviour are the accretion of the beach-dune complex of Long Beach, the erosion of Cape Shoalwater, and the northward migration of the Willapa Bay ebb-tidal delta during all intervals. The net shoreline advance of Long Beach increases from 0.28 m/yr in pre-jetty conditions to 3.78 m/yr during Interval 4. The erosion of Cape Shoalwater and the northward migration of the Willapa Bay ebb-tidal delta are related to the northern migration of the Willapa Bay North Channel. Volume changes at the Grays Harbor, Willapa Bay, and Columbia River entrances and the Columbia River estuary are balanced against losses and gains due to littoral transport and sand supply from the Columbia River. Based on these sediment balances, we infer the following pathways: sand that eroded from the inlets and inner deltas at the Grays Harbor and Columbia River entrances moved offshore and northward to accrete the outer deltas and the beaches to the north; sand from the south flank of the Grays Harbor delta and shelf along Grayland Plains moved onshore to accrete the beach dune complex of Grayland Plains and moved northward to maintain accretion of the outer delta and the beach-dune complex of North Beach; sand that eroded from the south flank of the Columbia River delta and shelf along Clatsop Plains contributed to the accretion of the beach-dune complex of Clatsop Plains and the Columbia River outer delta. The net volume change for Interval 1 and 3 at the Grays Harbor entrance and for Interval 1 at the Columbia River entrance is erosion, whereas the net change for the other intervals is accretion. For the entire CRLC, there is a net loss of 185 Mm 3 for Interval 1, a net gain of 357 Mm 3 for Interval 2, and a net gain of 187 Mm 3 for Interval 3. These imbalances can be the result of incomplete bathymetric coverage of the bays and shoreface, uncertainties in the adjustments of vertical tidal datums, inconsistencies in the bathymetric data, and uncertainties in the sediment supply of the Columbia River.

Oregon, Washington↗

Distributed Processing of Projections of Large Datasets: A Preliminary Study

Modern information needs have resulted in very large amounts of data being used in geographic information systems. Problems arise when trying to project these data in a reasonable amount of time and accuracy, however. Current single-threaded methods can suffer from two problems: fast projection with poor accuracy, or accurate projection with long processing time. A possible solution may be to combine accurate interpolation methods and distributed processing algorithms to quickly and accurately convert digital geospatial data between coordinate systems. Modern technology has made it possible to construct systems, such as Beowulf clusters, for a low cost and provide access to supercomputer-class technology. Combining these techniques may result in the ability to use large amounts of geographic data in time-critical situations.

Open-File Report↗

Methodology and Results for the Assessment of Oil and Gas resources, National Petroleum Reserve, Alaska

Oil and gas resources in each of the 24 plays within the National Petroleum Reserve in Alaska (NPRA) were estimated using a play analysis. Assessors specified geologic attributes, risks, and number of prospects for each play. Some specifications established distributions, while others were given as single values. From this information, sizes of oil and gas accumulations were simulated using a Monte Carlo algorithm. The number of such accumulations considered in a given simulation run was obtained from the distribution of the number of prospects. Each prospect in each successful simulation run was risked. This process yielded size-frequency distributions and summary statistics for the various petroleum categories. Estimates of remaining resources from individual plays were then aggregated, and measures of uncertainty computed. Technically recoverable, undiscovered oil beneath the Federal part of NPRA likely ranges between 5.9 and 13.2 billion barrels, with a mean (expected) value of 9.3 billion barrels. Technically recoverable, undiscovered nonassociated natural gas for the same area likely ranges between 39.1 and 83.2 trillion cubic feet, with a mean (expected) value of 59.7 trillion cubic feet. Mean values of the corresponding associated dissolved gas and natural gas liquid are 10.3 trillion cubic feet and 1.4 billion barrels respectively.

Open-File Report↗

Extending the Beowulf Cluster to the Desktop

Existing Beowulf clusters are normally limited to the actual number of compute nodes physically connected to the cluster's network switch. At some point, it may become necessary to extend the size of the cluster beyond the capacity of the switch. By using existing computers on a campus network, one can extend the computing power of the cluster by including these machines during times they would normally be idle. We have done a survey of methods for extending clusters and devised a parallel solution to a computationally complex all-possible-regressions algorithm and tested it by using idle computing power on the campus network.

Open-File Report↗

Direct current resistivity profiling to study distribution of water in the unsaturated zone near the Amargosa Desert Research Site, Nevada

In order to study the distribution of water in the unsaturated zone and potential for ground-water recharge near the Amargosa Desert Research Site south of Beatty, Nevada, the U.S. Geological Survey collected direct-current resistivity measurements along three profiles in May 2003 using an eight-channel resistivity imaging system. Resistivity data were collected along profiles across the ADRS, across a poorly incised (distributary) channel system of the Amargosa River southwest of the ADRS, and across a well-incised flood plain of the Amargosa River northwest of the ADRS. This report describes results of an initial investigation to estimate the distribution of water in the unsaturated zone and to evaluate the shallow subsurface stratigraphy near the ADRS. The geophysical method of dc resistivity was employed by using automated data collection with numerous electrodes. "Cross sections" of resistivity, produced by using an inversion algorithm on the field data, at the three field sites are presented and interpreted.

Nevada↗

Digital shaded-relief map of Venezuela

The digital shaded relief map of Venezuela is a composite of 3-arc-second (approximately 90-meter) elevation data, captured during the Shuttle Radar Topography Mission (SRTM) in February 2000. The SRTM, a joint project between the National Geospatial-Intelligence Agency (NGA) and the National Aeronautics and Space Administration (NASA), provides the most comprehensive international digital elevation dataset ever produced. The 11-day flight mission aboard the U.S. Space Shuttle Endeavour obtained elevation data for about 80% of the world's landmass through the use of single-pass Interferometric Synthetic Aperture Radar (InSAR) technology. This publication includes a colorized relief image of Venezuela (GeoTIFF format) and the elevation data (ESRI Grid format) used for rendering the image. These datasets are intended to accompany U.S. Geological Survey Open-File Report 2005–1038 ( http://pubs.usgs.gov/of/2005/1038/ ) and U.S. Geological Survey Data Series 199 ( http://pubs.usgs.gov/ds/2006/199/ ). SRTM data in version 2.0 have been processed using the latest data-void filling algorithms and represent a significant improvement to the data used in version 1.0.

Open-File Report↗

A new method of edge detection for object recognition

Traditional edge detection systems function by returning every edge in an input image. This can result in a large amount of clutter and make certain vectorization algorithms less accurate. Accuracy problems can then have a large impact on automated object recognition systems that depend on edge information. A new method of directed edge detection can be used to limit the number of edges returned based on a particular feature. This results in a cleaner image that is easier for vectorization. Vectorized edges from this process could then feed an object recognition system where the edge data would also contain information as to what type of feature it bordered.

Open-File Report↗

Version 3.0 of EMINERS - Economic Mineral Resource Simulator

Quantitative mineral resource assessment, as developed by the U.S. Geological Survey (USGS), consists of three parts: (1) development of grade and tonnage mineral deposit models; (2) delineation of tracts permissive for each deposit type; and (3) probabilistic estimation of the numbers of undiscovered deposits for each deposit type. The estimate of the number of undiscovered deposits at different levels of probability is the input to the EMINERS (Economic Mineral Resource Simulator) program. EMINERS uses a Monte Carlo statistical process to combine probabilistic estimates of undiscovered mineral deposits with models of mineral deposit grade and tonnage to estimate mineral resources. Version 3.0 of the EMINERS program is available as this USGS Open-File Report 2004-1344. Changes from version 2.0 include updating 87 grade and tonnage models, designing new templates to produce graphs showing cumulative distribution and summary tables, and disabling economic filters. The economic filters were disabled because embedded data for costs of labor and materials, mining techniques, and beneficiation methods are out of date. However, the cost algorithms used in the disabled economic filters are still in the program and available for reference for mining methods and milling techniques. The release notes included with this report give more details on changes in EMINERS over the years. EMINERS is written in C++ and depends upon the Microsoft Visual C++ 6.0 programming environment. The code depends heavily on the use of Microsoft Foundation Classes (MFC) for implementation of the Windows interface. The program works only on Microsoft Windows XP or newer personal computers. It does not work on Macintosh computers. For help in using the program in this report, see the "Quick-Start Guide for Version 3.0 of EMINERS-Economic Mineral Resource Simulator" (W.J. Bawiec and G.T. Spanski, 2012, USGS Open-File Report 2009-1057, linked at right). It demonstrates how to execute EMINERS software using default settings and existing deposit models.

Open-File Report↗

Real-time forecasts of tomorrow's earthquakes in California: A new mapping tool

We have derived a multi-model approach to calculate time-dependent earthquake hazard resulting from earthquake clustering. This file report explains the theoretical background behind the approach, the specific details that are used in applying the method to California, as well as the statistical testing to validate the technique. We have implemented our algorithm as a real-time tool that has been automatically generating short-term hazard maps for California since May of 2002, at http://step.wr.usgs.gov

California↗

Incorporating uncertainty into mercury-offset decisions with a probabilistic network for National Pollutant Discharge Elimination System permit holders: An interim report

This interim report describes an alternative approach for evaluating the efficacy of using mercury (Hg) offsets to improve water quality. Hg-offset programs may allow dischargers facing higher-pollution control costs to meet their regulatory obligations by making more cost effective pollutant-reduction decisions. Efficient Hg management requires methods to translate that science and economics into a regulatory decision framework. This report documents the work in progress by the U.S. Geological Survey’s Western Geographic Science Center in collaboration with Stanford University toward developing this decision framework to help managers, regulators, and other stakeholders decide whether offsets can cost effectively meet the Hg total maximum daily load (TMDL) requirements in the Sacramento River watershed. Two key approaches being considered are: (1) a probabilistic approach that explicitly incorporates scientific uncertainty, cost information, and value judgments; and (2) a quantitative approach that captures uncertainty in testing the feasibility of Hg offsets. Current fate and transport-process models commonly attempt to predict chemical transformations and transport pathways deterministically. However, the physical, chemical, and biologic processes controlling the fate and transport of Hg in aquatic environments are complex and poorly understood. Deterministic models of Hg environmental behavior contain large uncertainties, reflecting this lack of understanding. The uncertainty in these underlying physical processes may produce similarly large uncertainties in the decisionmaking process. However, decisions about control strategies are still being made despite the large uncertainties in current Hg loadings, the relations between total Hg (HgT) loading and methylmercury (MeHg) formation, and the relations between control efforts and Hg content in fish. The research presented here focuses on an alternative analytical approach to the current use of safety factors and deterministic methods for Hg TMDL decision support, one that is fully compatible with an adaptive management approach. This alternative approach uses empirical data and informed judgment to provide a scientific and technical basis for helping National Pollutant Discharge Elimination System (NPDES) permit holders make management decisions. An Hg-offset system would be an option if a wastewater-treatment plant could not achieve NPDES permit requirements for HgT reduction. We develop a probabilistic decision-analytical model consisting of three submodels for HgT loading, MeHg, and cost mitigation within a Bayesian network that integrates information of varying rigor and detail into a simple model of a complex system. Hg processes are identified and quantified by using a combination of historical data, statistical models, and expert judgment. Such an integrated approach to uncertainty analysis allows easy updating of prediction and inference when observations of model variables are made. We demonstrate our approach with data from the Cache Creek watershed (a subbasin of the Sacramento River watershed). The empirical models used to generate the needed probability distributions are based on the same empirical models currently being used by the Central Valley Regional Water Quality Control Cache Creek Hg TMDL working group. The significant difference is that input uncertainty and error are explicitly included in the model and propagated throughout its algorithms. This work demonstrates how to integrate uncertainty into the complex and highly uncertain Hg TMDL decisionmaking process. The various sources of uncertainty are propagated as decision risk that allows decisionmakers to simultaneously consider uncertainties in remediation/implementation costs while attempting to meet environmental/ecologic targets. We must note that this research is on going. As more data are collected, the HgT and cost-mitigation submodels are updated and the uncertainties may be reduced. Subsequently, the value of using a probabilistic framework for estimating and explicitly stating these uncertainties within a decisionmaking process can be estimated when new data are collected. Future work includes the design and implementation of a Bayesian network decision support system (BN-DSS) to produce mitigation scenarios for offset-project evaluation in the Cache Creek watershed. The decisionmaker, a wastewater-treatment plant, is expected to evaluate potential Hg-offset programs in terms of changes in HgT load changes, MeHg-production potential, project cost, and other suitability criteria. Subsequently, scenarios can be analyzed by performing sensitivity analyses and ranking environmental and economic uncertainties in terms of the decisionmaker’s preferences and risk choices. Such an analysis allows decisionmakers and stakeholders to explore various scenarios and predict the consequences of different stated preferences over outcomes.

Open-File Report↗

Electrical resistivity surveys in Prospect Gulch, San Juan County, Colorado

Prospect Gulch is a major source of naturally occurring and mining related metals to Cement Creek, a tributary of the upper Animas River in southwestern Colorado. Efforts to improve water quality in the watershed have focused on Prospect Gulch because many of its abandoned mines and are located on federal lands. Information on sources and pathways of metals, and related ground-water flow, will be useful to help prioritize and develop remediation strategies. It has been shown that the occurrence of sulfate, aluminum, iron, zinc and other metals associated with historical mining and the natural weathering of pyritic rock is substantial. In this study, direct current resistivity surveys were conducted to determine the subsurface resistivity distribution and to identify faults and fractures that may act as ground-water conduits or barriers to flow. Five lines of resistivity data were collected in the vicinity of Prospect Gulch, and cross-section profiles were constructed from the field data using a two-dimensional inversion algorithm. The conductive anomalies in the profiles are most likely caused by wet or saturated rocks and sediments, clay rich deposits, or high TDS ground water. Resistive anomalies are likely bedrock, dry surficial and sub-surface deposits, or deposits of ferricrete.

Open-File Report↗

Multiple-Objective Stepwise Calibration Using Luca

This report documents Luca (Let us calibrate), a multiple-objective, stepwise, automated procedure for hydrologic model calibration and the associated graphical user interface (GUI). Luca is a wizard-style user-friendly GUI that provides an easy systematic way of building and executing a calibration procedure. The calibration procedure uses the Shuffled Complex Evolution global search algorithm to calibrate any model compiled with the U.S. Geological Survey's Modular Modeling System. This process assures that intermediate and final states of the model are simulated consistently with measured values.

Open-File Report↗

GP Workbench Manual: Technical Manual, User's Guide, and Software Guide

GP Workbench is an open-source general-purpose geophysical data processing software package written primarily for ground penetrating radar (GPR) data. It also includes support for several USGS prototype electromagnetic instruments such as the VETEM and ALLTEM. The two main programs in the package are GP Workbench and GP Wave Utilities. GP Workbench has routines for filtering, gridding, and migrating GPR data; as well as an inversion routine for characterizing UXO (unexploded ordinance) using ALLTEM data. GP Workbench provides two-dimensional (section view) and three-dimensional (plan view or time slice view) processing for GPR data. GP Workbench can produce high-quality graphics for reports when Surfer 8 or higher (Golden Software) is installed. GP Wave Utilities provides a wide range of processing algorithms for single waveforms, such as filtering, correlation, deconvolution, and calculating GPR waveforms. GP Wave Utilities is used primarily for calibrating radar systems and processing individual traces. Both programs also contain research features related to the calibration of GPR systems and calculating subsurface waveforms. The software is written to run on the Windows operating systems. GP Workbench can import GPR data file formats used by major commercial instrument manufacturers including Sensors and Software, GSSI, and Mala. The GP Workbench native file format is SU (Seismic Unix), and subsequently, files generated by GP Workbench can be read by Seismic Unix as well as many other data processing packages.

Open-File Report↗