USGS Science⌕ Search

SEARCH · USGS Science

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 523 records · Page 29Linked to original sources

Two-dimensional hydraulic analyses of Joachim Creek, De Soto, Missouri

A two-dimensional hydraulic model; water-surface profiles; and digital maps of water-surface elevation, velocities, and water depths were developed for a 6.7-mile reach of Joachim Creek within and near the city of De Soto, Missouri. Water-surface profiles were generated for the 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probability (10-, 25-, 50-, 100-, and 500-year recurrence interval) flows. Digital maps of water-surface elevation, water depth, and velocity were generated for the 1- and 0.2-percent annual exceedance probability flows. Water-surface elevations and inundation extents of generated profiles and maps were substantially lower than similar products produced for the 2019 flood-insurance study that included the study reach. The differences in water-surface elevations can be attributed to differences in input streamflows and hydraulic simulation techniques. The water-surface elevations generated for the 1- and 0.2-percent annual exceedance probability flows were used to assess the vulnerability and inundation depths of 231 selected structures within the city of De Soto. Results indicate that 157 to 177 of the 231 structures were affected at the 1-percent annual exceedance probability flow, depending on the adjacent grade elevation used for reference. Between 185 and 198 structures were affected at the 0.2-percent annual exceedance probability flow, depending on grade elevation. Inundation depths at the affected structures were 0.02 to 9.28 feet (ft), depending on the flow and adjacent grade reference. Flood elevations were computed for Joachim Creek using a two-dimensional, finite-volume numerical modeling application for river hydraulics. The hydraulic model was calibrated using high-water marks from the April 18, 2013, flood and the maximum measured streamflow at the U.S. Geological Survey streamgage Joachim Creek at De Soto, Mo. (station 07019500), on September 8, 2018. The calibrated model was then used to compute the hydraulic conditions associated with the 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probability flows. The simulated water-surface elevations and digital elevation model (derived from light detection and ranging data having a 0.60-ft vertical accuracy and a 1.97-ft horizontal resolution) were used to generate products including water-surface profiles and maps of inundated area, water depth, and velocities using model postprocessing software.

Missouri↗

Stochastic simulation of spatially correlated geo-processes

In this study, developments in the theory of stochastic simulation are discussed. The unifying element is the notion of Radon projection in Euclidean spaces. This notion provides a natural way of reconstructing the real process from a corresponding process observable on a reduced dimensionality space, where analysis is theoretically easier and computationally tractable. Within this framework, the concept of space transformation is defined and several of its properties, which are of significant importance within the context of spatially correlated processes, are explored. The turning bands operator is shown to follow from this. This strengthens considerably the theoretical background of the geostatistical method of simulation, and some new results are obtained in both the space and frequency domains. The inverse problem is solved generally and the applicability of the method is extended to anisotropic as well as integrated processes. Some ill-posed problems of the inverse operator are discussed. Effects of the measurement error and impulses at origin are examined. Important features of the simulated process as described by geomechanical laws, the morphology of the deposit, etc., may be incorporated in the analysis. The simulation may become a model-dependent procedure and this, in turn, may provide numerical solutions to spatial-temporal geologic models. Because the spatial simu??lation may be technically reduced to unidimensional simulations, various techniques of generating one-dimensional realizations are reviewed. To link theory and practice, an example is computed in detail. ?? 1987 International Association for Mathematical Geology.

Mathematical Geology↗

Relations between continuous real-time physical properties and discrete water-quality constituents in the Little Arkansas River, south-central Kansas, 1998-2014

Water from the Little Arkansas River is used as source water for artificial recharge of the Equus Beds aquifer, one of the primary water-supply sources for the city of Wichita, Kansas. The U.S. Geological Survey has operated two continuous real-time water-quality monitoring stations since 1995 on the Little Arkansas River in Kansas. Regression models were developed to establish relations between discretely sampled constituent concentrations and continuously measured physical properties to compute concentrations of those constituents of interest. Site-specific regression models were originally published in 2000 for the near Halstead and near Sedgwick U.S. Geological Survey streamgaging stations and the site-specific regression models were then updated in 2003. This report updates those regression models using discrete and continuous data collected during May 1998 through August 2014. In addition to the constituents listed in the 2003 update, new regression models were developed for total organic carbon. The real-time computations of water-quality concentrations and loads are available at http://nrtwq.usgs.gov . The water-quality information in this report is important to the city of Wichita because water-quality information allows for real-time quantification and characterization of chemicals of concern (including chloride), in addition to nutrients, sediment, bacteria, and atrazine transported in the Little Arkansas River. The water-quality information in this report aids in the decision making for water treatment before artificial recharge.

Kansas↗

Floods in Central Texas, September 7-14, 2010

Severe flooding occurred near the Austin metropolitan area in central Texas September 7–14, 2010, because of heavy rainfall associated with Tropical Storm Hermine. The U.S. Geological Survey, in cooperation with the Upper Brushy Creek Water Control and Improvement District, determined rainfall amounts and annual exceedance probabilities for rainfall resulting in flooding in Bell, Williamson, and Travis counties in central Texas during September 2010. We documented peak streamflows and the annual exceedance probabilities for peak streamflows recorded at several streamflow-gaging stations in the study area. The 24-hour rainfall total exceeded 12 inches at some locations, with one report of 14.57 inches at Lake Georgetown. Rainfall probabilities were estimated using previously published depth-duration frequency maps for Texas. At 4 sites in Williamson County, the 24-hour rainfall had an annual exceedance probability of 0.002. Streamflow measurement data and flood-peak data from U.S. Geological Survey surface-water monitoring stations (streamflow and reservoir gaging stations) are presented, along with a comparison of September 2010 flood peaks to previous known maximums in the periods of record. Annual exceedance probabilities for peak streamflow were computed for 20 streamflow-gaging stations based on an analysis of streamflow-gaging station records. The annual exceedance probability was 0.03 for the September 2010 peak streamflow at the Geological Survey's streamflow-gaging stations 08104700 North Fork San Gabriel River near Georgetown, Texas, and 08154700 Bull Creek at Loop 360 near Austin, Texas. The annual exceedance probability was 0.02 for the peak streamflow for Geological Survey's streamflow-gaging station 08104500 Little River near Little River, Texas. The lack of similarity in the annual exceedance probabilities computed for precipitation and streamflow might be attributed to the small areal extent of the heaviest rainfall over these and the other gaged watersheds.

Texas↗

Computing daily mean streamflow at ungaged locations in Iowa by using the Flow Anywhere and Flow Duration Curve Transfer statistical methods

The U.S. Geological Survey (USGS) maintains approximately 148 real-time streamgages in Iowa for which daily mean streamflow information is available, but daily mean streamflow data commonly are needed at locations where no streamgages are present. Therefore, the USGS conducted a study as part of a larger project in cooperation with the Iowa Department of Natural Resources to develop methods to estimate daily mean streamflow at locations in ungaged watersheds in Iowa by using two regression-based statistical methods. The regression equations for the statistical methods were developed from historical daily mean streamflow and basin characteristics from streamgages within the study area, which includes the entire State of Iowa and adjacent areas within a 50-mile buffer of Iowa in neighboring states. Results of this study can be used with other techniques to determine the best method for application in Iowa and can be used to produce a Web-based geographic information system tool to compute streamflow estimates automatically. The Flow Anywhere statistical method is a variation of the drainage-area-ratio method, which transfers same-day streamflow information from a reference streamgage to another location by using the daily mean streamflow at the reference streamgage and the drainage-area ratio of the two locations. The Flow Anywhere method modifies the drainage-area-ratio method in order to regionalize the equations for Iowa and determine the best reference streamgage from which to transfer same-day streamflow information to an ungaged location. Data used for the Flow Anywhere method were retrieved for 123 continuous-record streamgages located in Iowa and within a 50-mile buffer of Iowa. The final regression equations were computed by using either left-censored regression techniques with a low limit threshold set at 0.1 cubic feet per second (ft3/s) and the daily mean streamflow for the 15th day of every other month, or by using an ordinary-least-squares multiple linear regression method and the daily mean streamflow for the 15th day of every other month. The Flow Duration Curve Transfer method was used to estimate unregulated daily mean streamflow from the physical and climatic characteristics of gaged basins. For the Flow Duration Curve Transfer method, daily mean streamflow quantiles at the ungaged site were estimated with the parameter-based regression model, which results in a continuous daily flow-duration curve (the relation between exceedance probability and streamflow for each day of observed streamflow) at the ungaged site. By the use of a reference streamgage, the Flow Duration Curve Transfer is converted to a time series. Data used in the Flow Duration Curve Transfer method were retrieved for 113 continuous-record streamgages in Iowa and within a 50-mile buffer of Iowa. The final statewide regression equations for Iowa were computed by using a weighted-least-squares multiple linear regression method and were computed for the 0.01-, 0.05-, 0.10-, 0.15-, 0.20-, 0.30-, 0.40-, 0.50-, 0.60-, 0.70-, 0.80-, 0.85-, 0.90-, and 0.95-exceedance probability statistics determined from the daily mean streamflow with a reporting limit set at 0.1 ft 3 /s. The final statewide regression equation for Iowa computed by using left-censored regression techniques was computed for the 0.99-exceedance probability statistic determined from the daily mean streamflow with a low limit threshold and a reporting limit set at 0.1 ft 3 /s. For the Flow Anywhere method, results of the validation study conducted by using six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 1,016 to 138 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 1,690 to 237 ft 3 /s. Values of the percent root-mean-square error ranged from 115 percent to 26.2 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 13.0 to 5.3 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.80 to 0.40. Percent-bias values ranged from 25.4 to 4.0 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.35. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.86 to 0.56. For the streamgage with the best agreement between observed and estimated streamflow, higher streamflows appear to be underestimated. For the streamgage with the worst agreement between observed and estimated streamflow, low flows appear to be overestimated whereas higher flows seem to be underestimated. Estimated cumulative streamflows for the period October 1, 2004, to September 30, 2009, are underestimated by -25.8 and -7.4 percent for the closest and poorest comparisons, respectively. For the Flow Duration Curve Transfer method, results of the validation study conducted by using the same six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 437 to 93.9 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 906 to 169 ft 3 /s. Values of the percent root-mean-square-error ranged from 67.0 to 25.6 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 12.5 to 4.4 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.79 to 0.40. Percent-bias values ranged from 22.7 to 0.94 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.38. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.89 to 0.48. For the streamgage with the closest agreement between observed and estimated streamflow, there is relatively good agreement between observed and estimated streamflows. For the streamgage with the poorest agreement between observed and estimated streamflow, streamflows appear to be substantially underestimated for much of the time period. Estimated cumulative streamflow for the period October 1, 2004, to September 30, 2009, are underestimated by -9.3 and -22.7 percent for the closest and poorest comparisons, respectively.

Illinois;Iowa;Minnesota;Missouri;Nebraska;Wisconsi↗

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana↗

Probabilistic estimates of surface coseismic slip and afterslip for Hayward fault earthquakes

We examine the partition of long‐term geologic slip on the Hayward fault into interseismic creep, coseismic slip, and afterslip. Using Monte Carlo simulations, we compute expected coseismic slip and afterslip at three alinement array sites for Hayward fault earthquakes with nominal moment magnitudes ranging from about 6.5 to 7.1. We consider how interseismic creep might affect the coseismic slip distribution as well as the variability in locations of large and small slip patches and the magnitude of an earthquake for a given rupture area. We calibrate the estimates to be consistent with the ratio of interseismic creep rate at the alinement array sites to the geologic slip rate for the Hayward fault. We find that the coseismic slip at the surface is expected to comprise only a small fraction of the long‐term geologic slip. The median values of coseismic slip are less than 0.2 m in nearly all cases as a result of the influence of interseismic creep and afterslip. However, afterslip makes a substantial contribution to the long‐term geologic slip and may be responsible for up to 0.5–1.5 m (median plus one standard deviation [S.D.]) of additional slip following an earthquake rupture. Thus, utility and transportation infrastructure could be severely impacted by afterslip in the hours and days following a large earthquake on the Hayward fault that generated little coseismic slip. Inherent spatial variability in earthquake slip combined with the uncertainty in how interseismic creep affects coseismic slip results in large uncertainties in these slip estimates.

California↗

Developing a state water plan: Ground-water conditions in Utah, spring of 1964

This report is the first in a series of annual reports which will describe ground-water conditions in Utah. It was prepared cooperatively by the U.S. Geological Survey and the Utah Water and Power Board and was designed to provide the data for interested parties, such as legislators, administrators, and planners to keep abreast of changing ground-water conditions in the state. Because this report is the first of the series, it necessarily includes certain background and descriptive information which gives a broad general picture of ground-water conditions. Subsequent reports will discuss only changes that have taken place during the previous year. Many of the data used in the preparation of the report were collected by the Geological Survey in cooperation with the Utah State Engineer during past and continuing programs. The well-location map and some statistical information about numbers of wells in the State were prepared by digital computer from the Utah Resources Information System, University of Utah, utilizing records which were compiled largely from the files of the Utah State Engineer. R.E. Marsell, geological consultant to the Utah Water and Power Board, first suggested that this report be prepared.

Utah↗

Partnerships for progress at the U.S. Geological Survey

This is about opportunity for the private sector. It is about combining the research capabilities of Government scientists with the commercial development potential of private companies. It is, consequently, about partnerships leading to products and services to enhance the quality of life and strengthen the American economy. The image at the left shows a computer screen image of Washington, DC using RevPG software for map revision.

General Information Product↗

Undiscovered lode tin resources of the Seward Peninsula, Alaska

The United States is a net importer of many important minerals, including tin. Consumption of primary tin in the United States is about 36,000 metric tons per year. Identified U.S. tin resources consist of about 40,000 metric tons. Although such figures provide insight about vulnerability to supply disruptions in the short term, they do not provide information about the potential of a country to meet its needs for a commodity from undiscovered domestic sources. Recent developments, including the preparation of models of mineral deposits and their grades and tonnages, and the application of computer simulation techniques to the estimation of metallic mineral resources, make it possible to estimate the magnitude of undiscovered resources, by deposit type, for relatively small areas such as the Seward Peninsula. This paper uses these developments and geophysical data to estimate undiscovered lode tin resources on the Seward Peninsula of Alaska. The assessment is based on a three-step methodology that (1) identifies the types of tin deposits that may be present in the region, (2) identifies the geophysical characteristics of unroofed granites and shallow granitoids, and (3) estimates, on the basis of various combinations of geologic and geophysical conditions, the number of undiscovered deposits present within or near exposed or concealed granite plutons. Computer simulation was used to combine the estimates of the number of deposits with available grade and tonnage models. Simulation experiments were designed to estimate the amount of tin in undiscovered deposits under a variety of limiting conditions. Results of simulation experiments indicate there is a 90 percent chance that the Seward Peninsula contains at least 51,000 metric tons, a 50 percent chance that it contains at least 390,000 metric tons, and a 10 percent chance it contains at least 1,100,000 metric tons of tin in undiscovered greisen, vein, and replacement deposits that have average grades of at least 0.5 percent tin.

Economic Geology↗

Sources of shaking and flooding during the Tohoku-Oki earthquake: a mixture of rupture styles

Modeling strong ground motions from great subduction zone earthquakes is one of the great challenges of computational seismology. To separate the rupture characteristics from complexities caused by 3D sub-surface geology requires an extraordinary data set such as provided by the recent Mw9.0 Tohoku-Oki earthquake. Here we combine deterministic inversion and dynamically guided forward simulation methods to model over one thousand high-rate GPS and strong motion observations from 0 to 0.25 Hz across the entire Honshu Island. Our results display distinct styles of rupture with a deeper generic interplate event (~Mw8.5) transitioning to a shallow tsunamigenic earthquake (~Mw9.0) at about 25 km depth in a process driven by a strong dynamic weakening mechanism, possibly thermal pressurization. This source model predicts many important features of the broad set of seismic, geodetic and seafloor observations providing a major advance in our understanding of such great natural hazards.

Earth and Planetary Science Letters↗

Estimates of ground-water recharge based on streamflow-hydrograph methods: Pennsylvania

This study, completed by the U.S. Geological Survey (USGS) in cooperation with the Pennsylvania Department of Conservation and Natural Resources, Bureau of Topographic and Geologic Survey (T&GS), provides estimates of ground-water recharge for watersheds throughout Pennsylvania computed by use of two automated streamflow-hydrograph-analysis methods--PART and RORA. The PART computer program uses a hydrograph-separation technique to divide the streamflow hydrograph into components of direct runoff and base flow. Base flow can be a useful approximation of recharge if losses and interbasin transfers of ground water are minimal. The RORA computer program uses a recession-curve displacement technique to estimate ground-water recharge from each storm period indicated on the streamflow hydrograph. Recharge estimates were made using streamflow records collected during 1885-2001 from 197 active and inactive streamflow-gaging stations in Pennsylvania where streamflow is relatively unaffected by regulation. Estimates of mean-annual recharge in Pennsylvania computed by the use of PART ranged from 5.8 to 26.6 inches; estimates from RORA ranged from 7.7 to 29.3 inches. Estimates from the RORA program were about 2 inches greater than those derived from the PART program. Mean-monthly recharge was computed from the RORA program and was reported as a percentage of mean-annual recharge. On the basis of this analysis, the major ground-water recharge period in Pennsylvania typically is November through May; the greatest monthly recharge typically occurs in March.

Pennsylvania↗

Water-quality data from oil and gas wells in part of the Permian Basin, southeastern New Mexico and western Texas

Approximately 8,000 chemical analyses of water produced from formations of several geologic ages in south eastern New Mexico and western Texas are tabulated by both geographic location and geologic formation. Empirical and mathematical relationships between (1) the dissolved solids and the measured and computed resistivity of water, the chloride-ion concentration and the density, and (2) the concentration of chloride-ion and the density and the resistivity of water, were established. Resistivities of water computed following a method developed by Moore, Szasz, and Whitney (1966) using the concentration of the six most common constituents dissolved in natural water correlate favorably with the measured resistivities of natural water. The water analyses were obtained chiefly from oil companies and oil-service companies operating in the western part of the Permian basin. The data were associated with the proper well by searching location, operator, and reference number cross indexes prepared from the Permian Basin Well Data System files. The cross indexes were determined to be essential to efficient data processing operations. After preliminary screening for nonrepresentative data, all data were indexed to the Permian Basin Well Data System magnetic tape file scout records by the unique reference number. The data were then coded, punched into tabulating cards, and stored on magnetic tape for computer operations. Data verification steps not normally required when using data received from a laboratory following standardized analytical techniques on samples obtained under standardized and controlled procedures were performed in order to upgrade reliability of the information. Data errors detected included poorly located or identified samples; mistakes including reproducing, or key punching; laboratory errors; inconsistent reporting; and nonrepresentative or contaminated samples.

New Mexico, Texas↗

The Water Resources Division water level recorder rental program; history and operation

The U.S. Geological Survey introduced a new water level recorder in the early 1960 's to automate computation of streamflow records. At the same time it developed a rental program as a means of providing the recorders to field offices. This report documents the experience gained over the past 20 years and offers suggestions for operation of a rental program in the future, as preparation is made to introduce a third generation of water level recorders. (USGS)

Open-File Report↗

Earthquake locations determined by the Southern Alaska seismograph network for October 1971 through May 1989

This report describes the instrumentation and evolution of the U.S. Geological Survey’s regional seismograph network in southern Alaska, provides phase and hypocenter data for seismic events from October 1971 through May 1989, reviews the location methods used, and discusses the completeness of the catalog and the accuracy of the computed hypocenters. Included are arrival time data for explosions detonated under the Trans-Alaska Crustal Transect (TACT) in 1984 and 1985. The U.S. Geological Survey (USGS) operated a regional network of seismographs in southern Alaska from 1971 to the mid 1990s. The principal purpose of this network was to record seismic data to be used to precisely locate earthquakes in the seismic zones of southern Alaska, delineate seismically active faults, assess seismic risks, document potential premonitory earthquake phenomena, investigate current tectonic deformation, and study the structure and physical properties of the crust and upper mantle. A task fundamental to all of these goals was the routine cataloging of parameters for earthquakes located within and adjacent to the seismograph network. The initial network of 10 stations, 7 around Cook Inlet and 3 near Valdez, was installed in 1971. In subsequent summers additions or modifications to the network were made. By the fall of 1973, 26 stations extended from western Cook Inlet to eastern Prince William Sound, and 4 stations were located to the east between Cordova and Yakutat. A year later 20 additional stations were installed. Thirteen of these were placed along the eastern Gulf of Alaska with support from the National Oceanic and Atmospheric Administration (NOAA) under the Outer Continental Shelf Environmental Assessment Program to investigate the seismicity of the outer continental shelf, a region of interest for oil exploration. Since then the region covered by the network remained relatively fixed while efforts were made to make the stations more reliable through improved electronic instrumentation and strengthened antenna systems. The majority of the stations installed since 1980 were operated only temporarily (from one to several years) for special studies in various areas within the network. Due to reduced funding, the network was trimmed substantially in the summer of 1985 with the closure of 15 stations, 13 of which were located in and around the Yakataga seismic gap. To further reduce costs, two telephone circuits were dropped and multiple radio relays were installed in their place. This economy reduced the reliability of these telemetry links. In addition, data collection from the areas around Cordova and Yakutat was compromised by the necessity of relying on triggered event recording using PC-based systems (Rogers, 1993) that were not fully developed and which proved to be less reliable than anticipated. The principal means of recording throughout the time period of this catalog was 20-channel oscillographs on 16-mm film (Teledyne Geotech Develocorder, Model RF400 and 4000D). Initially one Develocorder was operated at the USGS Alaskan headquarters in Anchorage, but in 1972 recording was shifted to the National Oceanic and Atmospheric Administration (NOAA) Palmer Observatory (currently the West Coast and Alaska Tsunami Warning Center). The Develocorders were turned off at the end of May 1989, and after that time recording was done in digital format at the Geophysical Institute of the University of Alaska in Fairbanks (GIUA). Thus, this catalog covers the entire period of film recording.

Alaska↗

Documentation of computer program (FHB1) for assignment of transient specified-flow and specified-head boundaries in applications of the modular finite-diference ground-water flow model (MODFLOW)

A computer program called the Flow and Head Boundary Package (FHB1) was developed for the U.S. Geological Survey three-dimensional finite-difference modular ground-water flow model, commonly referred to as MODFLOW. FHB1 allows MODFLOW users to specify flow or head boundary conditions that vary at times other than starting and ending times of stress periods and associated time steps. Values of flow and (or) head at each time step are calculated by linear interpolation of user- specified values. The ability to assign variable flow and head conditions defined at times not corresponding with the model stress periods allows greater flexibility in simulating natural geohydrologic systems and, at the same time, improves the efficiency of the methods used to represent these systems. The package also provides a way to apply specified-flow and specified-head boundaries in embedded, or nested, smaller-scale models using flow and (or) head values from larger-scale models. Using FHB1, the two models can have different simulation stress periods and time steps. Specification of variable-flow pumped wells in ground-water models is another example application.

Open-File Report↗