Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Inverse Problems”

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 145 records · Page 8Linked to original sources
USGS1995

Analysis of nutrient and ancillary water-quality data for surface and ground water of the Willamette Basin, Oregon, 1980-90

An analysis of historical water-quality data for surface and ground water collected in the Willamette and Sandy River Basins during the 1980-90 water years was performed. For surface water, most data were concentrated at sites on the main stem Willamette River or near the mouths of major tributaries. All seasons were represented. Data for nitrogen and phosphorus species were readily available, but simultaneously collected discharge measurements frequently were not. Seven primary sites were used for a quantitative analysis of nutrient data in surface water. At six of the seven primary sites, median concentrations of nitrite plus nitrate and total phosphorus were less than 0.5 milligrams per liter (mg/L) as nitrogen (N) and 0.1 mg/L as phosphorus (P), respectively. These concen- trations were lower than national median concen- trations for basins with similar land uses. At the Pudding River site, which received significant point-source and agricultural nonpoint-source inputs, median values for nitrite plus nitrate and total phosphorus were 1.4 mg/L as N and 0.2 mg/L as P, respectively. Nitrite-plus-nitrate concentrations at this site were seasonally dependent, with the highest concentrations occurring during winter. Forested sites had significantly lower water temperatures and nutrient concentrations than urban or agricultural sites. Evidence of diel variations in dissolved oxygen concentrations and pH values suggest that, at some sites, low dissolved oxygen concentrations may be a problem during early morning hours in the summer. Historical nutrient data for ground water were limited primarily to nitrate and nitrite-plus-nitrate determinations for wells completed in the basin-fill and alluvial aquifer. Nitrate and nitrite-plus-nitrate concentrations in this aquifer showed a weak inverse relationship with depth.

Water-Resources Investigations Report↗
USGS2002

Groundwater contamination downstream of a contaminant penetration site. II. Horizontal penetration of the contaminant plume

Part I of this study (Rubin, H.; Buddemeier, R.W. Groundwater Contamination Downstream of a Contaminant Penetration Site Part 1: Extension-Expansion of the Contaminant Plume. J. of Environmental Science and Health Part A (in press).) addressed cases, in which a comparatively thin contaminated region represented by boundary layers (BLs) developed within the freshwater aquifer close to contaminant penetration site. However, at some distance downstream from the penetration site, the top of the contaminant plume reaches the top or bottom of the aquifer. This is the location of the "attachment point," which comprises the entrance cross section of the domain evaluated by the present part of the study. It is shown that downstream from the entrance cross section, a set of two BLs develop in the aquifer, termed inner and outer BLs. It is assumed that the evaluated domain, in which the contaminant distribution gradually becomes uniform, can be divided into two sections, designated: (a) the restructuring section, and (b) the establishment section. In the restructuring section, the vertical concentration gradient leads to expansion of the inner BL at the expense of the outer BL, and there is almost no transfer of contaminant mass between the two layers. In the establishment section, each of the BLs occupies half of the aquifer thickness, and the vertical concentration gradient leads to transfer of contaminant mass from the inner to the outer BL. By use of BL approximations, changes of salinity distribution in the aquifer are calculated and evaluated. The establishment section ends at the uniformity point, downstream from which the contaminant concentration profile is practically uniform. The length of the restructuring section, as well as that of the establishment section, is approximately proportional to the aquifer thickness squared, and is inversely proportional to the transverse dispersivity. The study provides a convenient set of definitions and terminology that are helpful in visualizing the gradual development of uniform contaminant concentration distribution in an aquifer subject to contaminant plume penetration. The method developed in this study can be applied to a variety of problems associated with groundwater quality, such as initial evaluation of field data, design of field data collection, the identification of appropriate boundary conditions for numerical models, selection of appropriate numerical modeling approaches, interpretation and evaluation of field monitoring results, etc.

Journal of Environmental Science and Health - Part↗
USGS2015

Point spread functions for earthquake source imaging: An interpretation based on seismic interferometry

Recently, various methods have been proposed and applied for earthquake source imaging, and theoretical relationships among the methods have been studied. In this study, we make a follow-up theoretical study to better understand the meanings of earthquake source imaging. For imaging problems, the point spread function (PSF) is used to describe the degree of blurring and degradation in an obtained image of a target object as a response of an imaging system. In this study, we formulate PSFs for earthquake source imaging. By calculating the PSFs, we find that waveform source inversion methods remove the effect of the PSF and are free from artifacts. However, the other source imaging methods are affected by the PSF and suffer from the effect of blurring and degradation due to the restricted distribution of receivers. Consequently, careful treatment of the effect is necessary when using the source imaging methods other than waveform inversions. Moreover, the PSF for source imaging is found to have a link with seismic interferometry with the help of the source-receiver reciprocity of Green’s functions. In particular, the PSF can be related to Green’s function for cases in which receivers are distributed so as to completely surround the sources. Furthermore, the PSF acts as a low-pass filter. Given these considerations, the PSF is quite useful for understanding the physical meaning of earthquake source imaging.

Geophysical Journal International↗
USGS2010

Estimating 14C groundwater ages in a methanogenic aquifer

This paper addresses the problem of 14 C age dating of groundwaters in a confined regional aquifer affected by methanogenesis. Increasing CH 4 concentrations along the groundwater flow system and 13 C and 14 C isotopic data for dissolved inorganic carbon, dissolved organic carbon, and CH 4 clearly show the effect of methanogenesis on groundwater chemistry. Inverse reaction path modeling using NETPATH indicates the predominant geochemical reactions controlling the chemical evolution of groundwater in the aquifer are incongruent dissolution of dolomite, ion exchange, methanogenesis, and oxidation of sedimentary organic matter. Modeling of groundwater 14 C ages using NETPATH indicates that a significant part of groundwater in the Alliston aquifer is less than 13,000 years old; however, older groundwater in the range of 15,000–23,000 years is also present in the aquifer. This paper demonstrates that 14C ages calculated using NETPATH, incorporating the effects of methanogenesis on the carbon pools, provide reasonable groundwater ages that were not possible by other isotopic methods.

Water Resources Research↗
USGS2016

The Earthquake‐Source Inversion Validation (SIV) Project

Finite‐fault earthquake source inversions infer the (time‐dependent) displacement on the rupture surface from geophysical data. The resulting earthquake source models document the complexity of the rupture process. However, multiple source models for the same earthquake, obtained by different research teams, often exhibit remarkable dissimilarities. To address the uncertainties in earthquake‐source inversion methods and to understand strengths and weaknesses of the various approaches used, the Source Inversion Validation (SIV) project conducts a set of forward‐modeling exercises and inversion benchmarks. In this article, we describe the SIV strategy, the initial benchmarks, and current SIV results. Furthermore, we apply statistical tools for quantitative waveform comparison and for investigating source‐model (dis)similarities that enable us to rank the solutions, and to identify particularly promising source inversion approaches. All SIV exercises (with related data and descriptions) and statistical comparison tools are available via an online collaboration platform, and we encourage source modelers to use the SIV benchmarks for developing and testing new methods. We envision that the SIV efforts will lead to new developments for tackling the earthquake‐source imaging problem.

Seismological Research Letters↗
USGS2015

A random-walk algorithm for modeling lithospheric density and the role of body forces in the evolution of the Midcontinent Rift

This paper develops a Monte Carlo algorithm for extracting three-dimensional lithospheric density models from geophysical data. Empirical scaling relationships between velocity and density create a 3D starting density model, which is then iteratively refined until it reproduces observed gravity and topography. This approach permits deviations from uniform crustal velocity-density scaling, which provide insight into crustal lithology and prevent spurious mapping of crustal anomalies into the mantle. We test this algorithm on the Proterozoic Midcontinent Rift (MCR), north-central U.S. The MCR provides a challenge because it hosts a gravity high overlying low shear-wave velocity crust in a generally flat region. Our initial density estimates are derived from a seismic velocity/crustal thickness model based on joint inversion of surface-wave dispersion and receiver functions. By adjusting these estimates to reproduce gravity and topography, we generate a lithospheric-scale model that reveals dense middle crust and eclogitized lowermost crust within the rift. Mantle lithospheric density beneath the MCR is not anomalous, consistent with geochemical evidence that lithospheric mantle was not the primary source of rift-related magmas and suggesting that extension occurred in response to far-field stress rather than a hot mantle plume. Similarly, the subsequent inversion of normal faults resulted from changing far-field stress that exploited not only warm, recently faulted crust but also a gravitational potential energy low in the MCR. The success of this density modeling algorithm in the face of such apparently contradictory geophysical properties suggests that it may be applicable to a variety of tectonic and geodynamic problems.

Iowa, Kansas, Michigan, Minnesota, Nebraska, Wisco↗
USGS2017

Improved efficiency of maximum likelihood analysis of time series with temporally correlated errors

Most time series of geophysical phenomena have temporally correlated errors. From these measurements, various parameters are estimated. For instance, from geodetic measurements of positions, the rates and changes in rates are often estimated and are used to model tectonic processes. Along with the estimates of the size of the parameters, the error in these parameters needs to be assessed. If temporal correlations are not taken into account, or each observation is assumed to be independent, it is likely that any estimate of the error of these parameters will be too low and the estimated value of the parameter will be biased. Inclusion of better estimates of uncertainties is limited by several factors, including selection of the correct model for the background noise and the computational requirements to estimate the parameters of the selected noise model for cases where there are numerous observations. Here, I address the second problem of computational efficiency using maximum likelihood estimates (MLE). Most geophysical time series have background noise processes that can be represented as a combination of white and power-law noise, 1 / f α "> 1 / f α 1/fα with frequency, f . With missing data, standard spectral techniques involving FFTs are not appropriate. Instead, time domain techniques involving construction and inversion of large data covariance matrices are employed. Bos et al. (J Geod, 2013 . doi: 10.1007/s00190-012-0605-0 ) demonstrate one technique that substantially increases the efficiency of the MLE methods, yet is only an approximate solution for power-law indices >1.0 since they require the data covariance matrix to be Toeplitz. That restriction can be removed by simply forming a data filter that adds noise processes rather than combining them in quadrature. Consequently, the inversion of the data covariance matrix is simplified yet provides robust results for a wider range of power-law indices.

Journal of Geodesy↗
USGS2015

Hydrogeochemical effects of a bulkhead in the Dinero mine tunnel, Sugar Loaf mining district, near Leadville, Colorado

The Dinero mine drainage tunnel is an abandoned, draining mine adit near Leadville, Colorado, that has an adverse effect on downstream water quality and aquatic life. In 2009, a bulkhead was constructed (creating a mine pool and increasing water-table elevations behind the tunnel) to limit drainage from the tunnel and improve downstream water quality. The goal of this study was to document changes to hydrology and water quality resulting from bulkhead emplacement, and to understand post-bulkhead changes in source water and geochemical processes that control mine-tunnel discharge and water quality. Comparison of pre-and post-bulkhead hydrology and water quality indicated that tunnel discharge and zinc and manganese loads decreased by up to 97 percent at the portal of Dinero tunnel and at two downstream sites (LF-537 and LF-580). However, some water-quality problems persisted at LF-537 and LF-580 during high-flow events and years, indicating the effects of the remaining mine waste in the area. In contrast, post-bulkhead water quality degraded at three upstream stream sites and a draining mine tunnel (Nelson tunnel). Water-quality degradation in the streams likely occurred from increased contributions of mine-pool groundwater to the streams. In contrast, water-quality degradation in the Nelson tunnel was likely from flow of mine-pool water along a vein that connects the Nelson tunnel to mine workings behind the Dinero tunnel bulkhead. Principal components analysis, mixing analysis, and inverse geochemical modeling using PHREEQC indicated that mixing and geochemical reactions (carbonate dissolution during acid weathering, precipitation of goethite and birnessite, and sorption of zinc) between three end-member water types generally explain the pre-and post-bulkhead water composition at the Dinero and Nelson tunnels. The three end members were (1) a relatively dilute groundwater having low sulfate and trace element concentrations; (2) mine pool water, and (3) water that flowed from a structure in front of the bulkhead after bulkhead emplacement. Both (2) and (3) had high sulfate and trace element concentrations. These results indicate how analysis of monitoring information can be used to understand hydrogeochemical changes resulting from bulkhead emplacement. This understanding, in turn, can help inform future decisions on the disposition of the remaining mine waste and water-quality problems in the area.

Colorado↗
USGS1967

The ternary system K2SO4 MgSO4 CaSO4

Melting and subsolidus relations in the system K 2 SO 4 MgSO 4 CaSO 4 were studied using heating-cooling curves, differential thermal analysis, optics, X-ray diffraction at room and high temperatures and by quenching techniques. Previous investigators were unable to study the binary MgSO 4 CaSO 4 system and the adjacent area in the ternary system because of the decomposition of MgSO 4 and CaSO 4 at high temperatures. This problem was partly overcome by a novel sealed-tube quenching method, by hydrothermal synthesis, and by long-time heating in the solidus. As a result of this study, we found: (1) a new compound, CaSO 4 ·3MgSO 4 (m.p. 1201°C) with a field extending into the ternary system; (2) a high temperature form of MgSO 4 with a sluggishly reversible inversion. An X-ray diffraction pattern for this polymorphic form is given; (3) the inversion of β-CaSO 4 (anhydrite) to α-CaSO 4 at 1195°C, in agreement with grahmann ; (1) (4) the melting point of MgSO 4 is 1136°C and that of CaSO 4 is 1462°C (using sealed tube methods to prevent decomposition of the sulphates); (5) calcium langbeinite (K 2 SO 4 ·2CaSO 4 ) is the only compound in the K 2 SO 4 CaSO 4 binary system. This resolved discrepancies in the results of previous investigators; (6) a continuous solid solution series between congruently melting K 2 SOP 4 ·2MgSO 4 (langbeinite) and incongruently melting K 2 SO 4 ·2CaSO 4 (calcium langbeinite); (7) the liquidus in the ternary system consists of primary phase fields of K 2 SO 4 , MgSO 4 , CaSO 4 , langbeinite-calcium langbeinite solid solution, and CaSO 4 ·3MgSO 4 . The CaSO 4 field extends over a large portion of the system. Previously reported fields for the compounds (K 2 SO 4 ·MgSO 4 ·nCaSO 4 ), K 2 SO 4 ·3CaSO 4 and K 2 SO 4 ·CaSO 4 were not found; (8) a minimum in the ternary system at: 740°C, 25% MgSO 4 , 6% CaSO 4 , 69% K 2 SO 4 ; and ternary eutectics at 882°C, 49% MgSO 4 , 19% CaSO 4 , 32% K 2 SO 4 ; and 880°, 67·5% MgSO 4 , 5% CaSO 4 , 27·5% K 2 SO 4 .

Journal of Inorganic and Nuclear Chemistry↗
USGS2006

Thickness distribution of a cooling pyroclastic flow deposit on Augustine Volcano, Alaska: Optimization using InSAR, FEMs, and an adaptive mesh algorithm

Interferometric synthetic aperture radar (InSAR) imagery documents the consistent subsidence, during the interval 1992–1999, of a pyroclastic flow deposit (PFD) emplaced during the 1986 eruption of Augustine Volcano, Alaska. We construct finite element models (FEMs) that simulate thermoelastic contraction of the PFD to account for the observed subsidence. Three-dimensional problem domains of the FEMs include a thermoelastic PFD embedded in an elastic substrate. The thickness of the PFD is initially determined from the difference between post- and pre-eruption digital elevation models (DEMs). The initial excess temperature of the PFD at the time of deposition, 640 °C, is estimated from FEM predictions and an InSAR image via standard least-squares inverse methods. Although the FEM predicts the major features of the observed transient deformation, systematic prediction errors (RMSE = 2.2 cm) are most likely associated with errors in the a priori PFD thickness distribution estimated from the DEM differences. We combine an InSAR image, FEMs, and an adaptive mesh algorithm to iteratively optimize the geometry of the PFD with respect to a minimized misfit between the predicted thermoelastic deformation and observed deformation. Prediction errors from an FEM, which includes an optimized PFD geometry and the initial excess PFD temperature estimated from the least-squares analysis, are sub-millimeter (RMSE = 0.3 mm). The average thickness (9.3 m), maximum thickness (126 m), and volume (2.1 × 10 7 m 3 ) of the PFD, estimated using the adaptive mesh algorithm, are about twice as large as the respective estimations for the a priori PFD geometry. Sensitivity analyses suggest unrealistic PFD thickness distributions are required for initial excess PFD temperatures outside of the range 500–800 °C.

Alaska↗
USGS2026

Computational electromagnetic geophysics for groundwater system studies: A review on established practices and recent advances

Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies. As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.

Journal of Hydrology↗
USGS1989

Short note; Detectability levels for central induction transient soundings

A question usually asked during the planning of a transient sounding survey is, "Can information useful to solving the geologic problem at hand be obtained from the measurements?" This question is usually answered by constructing a geologic model for the survey area based upon the best available information to determine which model parameters can be resolved. Specifically, this determination can be accomplished by computing a forward model, varying a parameter of interest, and observing whether the responses of the original and perturbed models are different. Alternatively, inversion-based methods can be used to estimate parameter resolution (Inman, 1975; Glenn and Ward, 1976).

Geophysics↗
USGS1999

A new multistage groundwater transport inverse method: presentation, evaluation, and implications

More computationally efficient methods of using concentration data are needed to estimate groundwater flow and transport parameters. This work introduces and evaluates a three‐stage nonlinear‐regression‐based iterative procedure in which trial advective‐front locations link decoupled flow and transport models. Method accuracy and efficiency are evaluated by comparing results to those obtained when flow‐ and transport‐model parameters are estimated simultaneously. The new method is evaluated as conclusively as possible by using a simple test case that includes distinct flow and transport parameters, but does not include any approximations that are problem dependent. The test case is analytical; the only flow parameter is a constant velocity, and the transport parameters are longitudinal and transverse dispersivity. Any difficulties detected using the new method in this ideal situation are likely to be exacerbated in practical problems. Monte‐Carlo analysis of observation error ensures that no specific error realization obscures the results. Results indicate that, while this, and probably other, multistage methods do not always produce optimal parameter estimates, the computational advantage may make them useful in some circumstances, perhaps as a precursor to using a simultaneous method.

Water Resources Research↗
USGS2020

Generalizing the inversion‐based PSHA source model for an interconnected fault system

This article represents a step toward generalizing and simplifying the procedure for constructing an inversion‐based seismic hazard source model for an interconnected fault system, including the specification of adjustable segmentation constraints. A very simple example is used to maximize understandability and to counter the notion that an inversion approach is only applicable when an abundance of data is available. Also exemplified is how to construct a range of models to adequately represent epistemic uncertainties (which should be a high priority in any hazard assessment). Opportunity is also taken to address common concerns and misunderstandings associated with the third Uniform California Earthquake Rupture Forecast, including the seemingly disproportionate number of large‐magnitude events, and how well hazard is resolved given the overall problem is very underdetermined. However, the main aim of this article is to provide a general protocol for constructing such models.

Bulletin of the Seismological Society of America↗
USGS1996

Earthquake mechanisms from linear-programming inversion of seismic-wave amplitude ratios

The amplitudes of radiated seismic waves contain far more information about earthquake source mechanisms than do first-motion polarities, but amplitudes are severely distorted by the effects of heterogeneity in the Earth. This distortion can be reduced greatly by using the ratios of amplitudes of appropriately chosen seismic phases, rather than simple amplitudes, but existing methods for inverting amplitude ratios are severely nonlinear and require computationally intensive searching methods to ensure that solutions are globally optimal. Searching methods are particularly costly if general (moment tensor) mechanisms are allowed. Efficient linear-programming methods, which do not suffer from these problems, have previously been applied to inverting polarities and wave amplitudes. We extend these methods to amplitude ratios, in which formulation on inequality constraint for an amplitude ratio takes the same mathematical form as a polarity observation. Three-component digital data for an earthquake at the Hengill-Grensdalur geothermal area in southwestern Iceland illustrate the power of the method. Polarities of P , SH , and SV waves, unusually well distributed on the focal sphere, cannot distinguish between diverse mechanisms, including a double couple. Amplitude ratios, on the other hand, clearly rule out the double-couple solution and require a large explosive isotropic component.

Bulletin of the Seismological Society of America↗
USGS2023

Reconstructing missing data by comparing interpolation techniques: Applications for long-term water quality data

Missing data are typical yet must be addressed for proper inferences or expanding datasets to guide our limnological understanding and management of aquatic systems. Interpolation methods (i.e., estimating missing values using known values within the dataset) can alleviate data gaps and common problems. We compared seven popular interpolation methods for predicting substantial missingness in a long-term water quality dataset from the Upper Mississippi River, U.S.A. The dataset included 80,000 sampling sites collected over 30 yr that had substantial missingness for total nitrogen (TN), total phosphorus (TP), and water velocity. For all three interpolated water quality variables, random forests had very high prediction accuracy and outperformed the methods of ordinary kriging, polynomial regressions, regression trees, and inverse distance weighting. TP had a mean absolute error (MAE) of 0.03 mg (L-TP) −1 , TN had a MAE of 0.39 mg (L-TN) −1 , and water velocity had a MAE of 0.10 m s −1 . The random forests' error rates were mapped and showed low spatiotemporal variability across the riverscape, indicating high model performance across many habitat types and large spatial scales. In the current era of “big data,” interpolation becomes an imperative step prior to ecological analyses yet remains unfamiliar and underutilized. Our research briefly describes the importance of addressing missingness and provides a roadmap to conduct model intercomparisons of other big datasets. We also share adaptable data analysis scripts, which allows others to readily conduct interpolation comparisons for many limnology applications and contexts.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗
USGS1991

Accumulation of trace elements and organochlorines by surf scoters wintering in the Pacific northwest

Selenium, cadmium, mercury, copper, manganese, zinc, aluminum, lead, PCBs and DDE were accumulated by segments of the surf scoter (Melanitta perspicillata) population that winters in the Pacific Northwest, but whether the uptake occurred on breeding and/or wintering grounds was uncertain for some contaminants. Surf scoters collected in Puget Sound and San Francisco Bay (in another study) during the same period (January 1985) contained similar concentrations of cadmium, but Alsea Bay scoters contained more. Cadmium was inversely related to both liver and body weights of Northwest scoters in January; similar weight losses were reported in experimental laboratory studies. Northwest and north San Francisco Bay scoters contained similar mercury concentrations, but those in south San Francisco Bay contained higher concentrations. San Francisco Bay scoters contained higher arsenic and selenium concentrations than those in the Northwest; however, the 43.4 ppm (geometric mean, dry wt) selenium in livers at Commencement Bay in January was above levels associated with the reproductive problems in aquatic birds at Kesterson National Wildlife Refuge. Even higher concentrations of some elements may be found in surf scoters in March, because a later collection (March) at San Francisco Bay yielded higher concentrations than found there in January. Trace element concentrations in birds at a given wintering location are variable among species and may be influenced by diet, breeding grounds, and physiology (e.g., at Commencement Bay surf scoters with a sediment-associated diet contained 50X more cadmium in their kidneys than did fish-eating western grebes [Aechmophorus occidentalis]). The numerous wildlife species that live on estuaries require further attention.

Northwestern Naturalist↗
USGS1977

SIPB: a seismic refraction inverse modeling program for batch computer systems

SIPB is an interactive Fortran computer program that was developed for use with a timeshare computer system with program control information submitted from a remote terminal, and output data displayed on the terminal or printed on a line printer. The program is an upgraded version of FSIPI (Scott, Tibbetts, and Burdick, 1972) with several major improvements in addition to .its adaptation to timeshare operation. The most significant improvement was made in the procedure for handling data from in-line offset shotpoints beyond the end shotpoints of the geophone spread. The changes and improvements are described, user's instructions are outlined, examples of input and output data for a test problem are presented, and the Fortran program is listed in this report. An upgraded batch-mode program, SIPB, is available for users who do not have a timeshare computer system available (Scott, 1977).

Open-File Report↗