USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modeling Report”

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 595 records · Page 33Linked to original sources

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗

Equation-free modeling unravels the behavior of complex ecological systems

Ye et al. (1) address a critical problem confronting the management of natural ecosystems: How can we make forecasts of possible future changes in populations to help guide management actions? This problem is especially acute for marine and anadromous fisheries, where the large interannual fluctuations of populations, arising from complex nonlinear interactions among species and with varying environmental factors, have defied prediction over even short time scales. The empirical dynamic modeling (EDM) described in Ye et al.’s report, the latest in a series of papers by Sugihara and his colleagues, offers a promising quantitative approach to building models using time series to successfully project dynamics into the future. With the term “equation-free” in the article title, Ye et al. (1) are suggesting broader implications of their approach, considering the centrality of equations in modern science. From the 1700s on, nature has been increasingly described by mathematical equations, with differential or difference equations forming the basic framework for describing dynamics. The use of mathematical equations for ecological systems came much later, pioneered by Lotka and Volterra, who showed that population cycles might be described in terms of simple coupled nonlinear differential equations. It took decades for Lotka–Volterra-type models to become established, but the development of appropriate differential equations is now routine in modeling ecological dynamics. There is no question that the injection of mathematical equations, by forcing “clarity and precision into conjecture” (2), has led to increased understanding of population and community dynamics. As in science in general, in ecology equations are a key method of communication and of framing hypotheses. These equations serve as compact representations of an enormous amount of empirical data and can be analyzed by the powerful methods of mathematics.

PNAS↗

Estimating Water Storage Capacity of Existing and Potentially Restorable Wetland Depressions in a Subbasin of the Red River of the North

Executive Summary Concern over flooding along rivers in the Prairie Pothole Region has stimulated interest in developing spatially distributed hydrologic models to simulate the effects of wetland water storage on peak river flows. Such models require spatial data on the storage volume and interception area of existing and restorable wetlands in the watershed of interest. In most cases, information on these model inputs is lacking because resolution of existing topographic maps is inadequate to estimate volume and areas of existing and restorable wetlands. Consequently, most studies have relied on wetland area to volume or interception area relationships to estimate wetland basin storage characteristics by using available surface area data obtained as a product from remotely sensed data (e.g., National Wetlands Inventory). Though application of areal input data to estimate volume and interception areas is widely used, a drawback is that there is little information available to provide guidance regarding the application, limitations, and biases associated with such approaches. Another limitation of previous modeling efforts is that water stored by wetlands within a watershed is treated as a simple lump storage component that is filled prior to routing overflow to a pour point or gaging station. This approach does not account for dynamic wetland processes that influence water stored in prairie wetlands. Further, most models have not considered the influence of human-induced hydrologic changes, such as land use, that greatly influence quantity of surface water inputs and, ultimately, the rate that a wetland basin fills and spills. The goals of this study were to (1) develop and improve methodologies for estimating and spatially depicting wetland storage volumes and interceptions areas and (2) develop models and approaches for estimating/simulating the water storage capacity of potentially restorable and existing wetlands under various restoration, land use, and climatic scenarios. To address these goals, we developed models and approaches to spatially represent storage volumes and interception areas of existing and potentially restorable wetlands in the upper Mustinka subbasin within Grant County, Minn. We then developed and applied a model to simulate wetland water storage increases that would result from restoring 25 and 50 percent of the farmed and drained wetlands in the upper Mustinka subbasin. The model simulations were performed during the growing season (May-October) for relatively wet (1993; 0.79 m of precipitation) and dry (1987; 0.40 m of precipitation) years. Results from the simulations indicated that the 25 percent restoration scenario would increase water storage by 21-24 percent and that a 50 percent scenario would increase storage by 34-38 percent. Additionally, we estimated that wetlands in the subbasin have potential to store 11.57-20.98 percent of the total precipitation that fell over the entire subbasin area (52,758 ha). Our simulation results indicated that there is considerable potential to enhance water storage in the subbasin; however, evaluation and calibration of the model is necessary before simulation results can be applied to management and planning decisions. In this report we present guidance for the development and application of models (e.g., surface area-volume predictive models, hydrology simulation model) to simulate wetland water storage to provide a basis from which to understand and predict the effects of natural or human-induced hydrologic alterations. In developing these approaches, we tried to use simple and widely available input data to simulate wetland hydrology and predict wetland water storage for a specific precipitation event or a series of events. Further, the hydrology simulation model accounted for land use and soil type, which influence surface water inputs to wetlands. Although information presented in this report is specific to the Mustinka subbasin, the approaches

Open-File Report↗

Reconnaissance of land-use sources of pesticides in drinking water, McKenzie River, Oregon

The Eugene Water and Electric Board (EWEB) provides water and electricity to the City of Eugene, Oregon, from the McKenzie River. In the spring of 2002, EWEB initiated a pesticide monitoring program in cooperation with the U.S. Geological Survey as part of their Drinking Water Source Protection Plan. Approximately twice yearly pesticide samples were collected from 2002 to 2010 at a suite of sampling sites representing varying land uses in the lower McKenzie River basin. A total of 117 ambient samples were collected from 28 tributary and mainstem sites, including those dominated by forestry, urban, and agricultural activities, as well as the mouths of major tributaries characterized by a mixture of upstream land use. Constituents tested included 175 compounds in filtered water (72 herbicides, 43 insecticides, 10 fungicides, and 36 of their degradation products, as well as 14 pharmaceutical compounds). No attempt was made to sample different site types equivalently; sampling was instead designed primarily to characterize representative storm events during spring and fall runoff conditions in order to assess or confirm the perceived importance of the different site types as sources for pesticides. Sampling was especially limited for agricultural sites, which were only sampled during two spring storm surveys. A total of 43 compounds were detected at least once, with many of these detected only at low concentrations (<0.1 micrograms per liter). Nine compounds were detected at the drinking- water intake, and most of these were reported as estimates less than the laboratory reporting level. Human-health benchmark concentrations were consistently several orders of magnitude higher than detected concentrations at the intake, indicating that pesticide concentrations present a negligible threat to human health. The largest number of pesticide detections occurred during spring storm surveys and primarily were associated with urban stormwater drains. Urban sites also were associated with the highest concentrations, occasionally exceeding 1 microgram per liter. Many of the compounds detected at urban sites were relatively hydrophobic (do not mix easily with water), persistent, and suspected of endocrine disruption. In contrast, forestry compounds were rarely detectable in the McKenzie River, even though forest land predominates in the basin and forestry pesticide use was detected in small tributaries draining forested lands following application. Agricultural pesticide runoff was not well characterized by the limited data available, although a large number of compounds was estimated to be used in the basin and concentrations were moderately high in the few samples collected from small tributaries draining agricultural lands. Results from this analysis indicate that urban pesticide use is potentially an important source for pesticides of concern for drinking water, not limited exclusively to storm conditions. Forestry pesticide use is not considered a likely threat to drinking water quality at the present time (2012). A more complete understanding of agricultural chemicals in runoff in the McKenzie River basin requires further investigation. In addition to evaluating the data collected in this study, a conceptual model describing pesticide contamination in the McKenzie River basin is provided, based on current scientific understanding that is consistent with the data analysis presented in this report. This model is intended to provide a foundation for future monitoring in the basin.

Oregon↗

Integration of a numerical model and remotely sensed data to study urban/rural land surface climate processes

Simulation of urban/rural land surface climate processes using boundary layer climate models requires accurate input data with regard to surface thermal and radiative properties. The research reported here resulted in development of a procedure to integrate the satellite-derived surface biophysical parameters with a boundary layer climate model for simulating spatial surface energy exchange. The procedure was tested through spatial surface energy balance simulation of an urban/rural landscape in eastern Nebraska. The modeled surface temperature and net radiation were compared to those derived from the concurrent satellite data. The errors of the modeled surface temperature were small, and were mainly attributed to uncertainties in the estimation of surface moisture availability and satellite-derived surface radiant temperature. Modeled net radiation was also in agreement with the values calculated from satellite data. Modeled turbulent heat fluxes were in general agreement as compared to those reported in the literature, but the model tended to overestimate the latent heat flux for most rural land cover types. It was concluded that by incorporation of satellite-derived surface physical parameters into a boundary layer model, simulation of spatial land surface climate processes was much improved. The method and procedures developed from this study can be utilized in other boundary layer climate models.

Nebraska↗

Methods for computing water-quality concentrations and loads at sites operated by the U.S. Geological Survey Kansas Water Science Center

The U.S. Geological Survey (USGS) Kansas Water Science Center (KSWSC) has published time-series computations of water-quality concentrations and loads based on in situ sensor data since 1995. Water-quality constituent concentrations or densities are computed using regression models that relate in situ sensor values to laboratory analyses of periodically collected samples. These regression models currently (2024) follow no uniform published guidance and are individually documented through USGS reports. This report describes updated (2024) procedures designed to improve the consistency, quality, and timeliness of computed continuous water-quality data produced by the USGS KSWSC. Beginning in 2024, models developed by the USGS KSWSC that follow specific procedures and requirements related to sample collection, model fit, and model documentation outlined in this report are planned to be published and stored in the USGS National Real-Time Water Quality Data for the Nation Data Service. This report also describes USGS KSWSC procedures for evaluating and publishing time-series water-quality computations after initial model development and documentation. This guidance can be used to improve USGS KSWSC model development and data computation consistency and streamline the time-series water-quality computation process from model development to publication.

Open-File Report↗

Spatiotemporal methane emission from global reservoirs

Inland aquatic systems, such as reservoirs, contribute substantially to global methane (CH 4 ) emissions; yet are among the most uncertain components of the total CH 4 budget. Reservoirs have received recent attention as they may generate high CH 4 fluxes. Improved quantification of these CH 4 fluxes, particularly their spatiotemporal distribution, is key to realistically incorporating them in CH 4 modeling and budget studies. Here we report on a new global, gridded (0.25° lat × 0.25° lon) study of reservoir CH 4 emissions, accounting for new knowledge regarding reservoir areal extent and distribution, and spatiotemporal emission patterns influenced by diurnal variability, temperature-dependent seasonality, satellite-derived freeze-thaw dynamics, and eco-climatic zone. The results of this new data set comprise daily CH 4 emissions throughout the full annual cycle and show that reservoirs cover 297 × 10 3 km 2 globally and emit 10.1 Tg CH 4 yr −1 (1σ uncertainty range of 7.2–12.9 Tg CH 4 yr −1 ) from diffusive (1.2 Tg CH 4 yr −1 ) and ebullitive (8.9 Tg CH 4 yr −1 ) emission pathways. This analysis of reservoir CH 4 emission addresses multiple gaps and uncertainties in previous studies and represents an important contribution to studies of the global CH 4 budget. The new data sets and methodologies from this study provide a framework to better understand and model the current and future role of reservoirs in the global CH 4 budget and to guide efforts to mitigate reservoir-related CH 4 emissions.

Journal of Geophysical Research: Biogeosciences↗

Summary of Survival Data from Juvenile Coho Salmon in the Klamath River, Northern California, 2007

A study to estimate the effects of Iron Gate Dam discharge on ESA-listed juvenile coho salmon during their seaward migration to the ocean was begun in 2005. Estimates of survival through various reaches of river downstream from the dam were completed in 2006 and 2007 as part of this process. This report describes the estimates of survival during 2007, and is a complement to a similar report from 2006. Further analyses will be included in a separate report. In 2007, a series of models were evaluated to determine what survival and capture probabilities of radio-tagged hatchery fish were in several reaches between Iron Gate Hatchery at river kilometer 309 and a site at river kilometer 33. The results indicate trends in survival among reaches were similar to those found in 2006, but the survival in 2007 was lower than in 2006. The differences in survivals from Iron Gate Hatchery to river kilometer 33 in 2006 (0.653 SE 0.039) and 2007 (0.497 SE 0.044) were caused primarily by differences in survivals upstream from the Scott River. This document is a brief summary of 2007 survival results.

Open-File Report↗

An Overview of the Factors Involved in Evaluating the Geochemical Effects of Highway Runoff on the Environment

Materials washed by rain and snowmelt from highways into adjacent surface waters, ground waters, and ecosystems can pollute water and affect biota. To understand the chemical behavior of any one of these materials and its effects on the environment requires knowledge of the chemistry of the material and how it interacts with other components in the local geochemical system. An integrated watershed approach, therefore, would be the most effective method to assess the effects of highway runoff on local receiving waters. Analysis of one or a few specific contaminants will provide limited and incomplete information and may be misleading in terms of environmental effects. This report addresses the background geochemistry required to model highway runoff and to make realistic assessments of the potential effects of runoff on the environment.

Open-File Report↗

Mapping ground water in three dimensions: An analysis of airborne geophysical surveys of the Upper San Pedro River Basin, Cochise County, southeastern Arizona

This report summarizes the results of two airborne geophysical surveys conducted in the upper San Pedro Valley of southeastern Arizona in 1997 and 1999. The combined surveys cover about 1,000 square kilometers and extend from the Huachuca Mountains on the west to the Mule Mountains and Tombstone Hills on the east and from north of the Babocomari River to near the Mexican border on the south. The surveys included the acquisition of high-resolution magnetic data, which were used to map depth to the crystalline basement rocks underlying the sediments filling the basin. The magnetic inversion results show a complex basement morphology, with sediment thickness in the center of the valley ranging from ~237 meters beneath the city of Sierra Vista to ~1,500 meters beneath Huachuca City and the Palominas area near the Mexican border. The surveys also included acquisition of 60-channel time-domain electromagnetic (EM) data. Extensive quality analyses of these data, including inversion to conductivity vs. depth (conductivity-depth-transform or CDT) profiles and comparisons with electrical well logs, show that the electrical conductor mapped represents the subsurface water-bearing sediments throughout most of the basin. In a few places (notably the mouth of Huachuca Canyon), the reported water table lies above where the electrical conductor places it. These exceptions appear to be due to a combination of outdated water-table information, significant horizontal displacement between the wells and the CDT profiles, and a subtle calibration issue with the CDT algorithm apparent only in areas of highly resistive (very dry) overburden. These occasional disparities appear in less than 5 percent of the surveyed area. Observations show, however, that wells drilled in the thick unsaturated zone along the Huachuca Mountain front eventually intersect water, at which point the water rapidly rises high into the unsaturated zone within the wellbore. This rising of water in a wellbore implies some sort of confinement below the thick unsaturated zone, a confinement that is not identified in the available literature. Occasional disparities notwithstanding, maps of the electrical conductor derived from the airborne EM system provide a synoptic view of the presence of water underlying the upper San Pedro Valley, including its three-dimensional distribution. The EM data even show faults previously only inferred from geologic mapping. The magnetic and electromagnetic data together appear to show the thickness of the sediments, the water in the saturated sediments down to a maximum of about 400 meters depth, and even places where the main ground-water body is not in direct contact with the San Pedro River. However, the geophysical data cannot reveal anything directly about hydraulic conductivity or ground-water flow. Estimating these characteristics requires new hydraulic modeling based in part on this report. One concern to reviewers of this report is the effect that clays may have on the electrical conductor mapped with the airborne geophysical system. Although the water in the basin is unusually conductive, averaging 338 microsiemens per centimeter, reasoning cited below suggests that the contribution of clays to the overall conductivity would be relatively small. Basic principles of sedimentary geology suggest that silts and clays should dominate the center of the basin, while sands and gravels would tend to dominate the margins. Although clay content may increase the amplitude of the observed electrical conductors somewhat, it will not affect the depths to the conductor derived from depth inversions. Further, fine-grained sediments generally have higher porosity and tend to lie toward a basin center, a fact in general agreement with the observed geophysical data.

Arizona↗

Population modeling for furbearer management

The management of furbearers has become increasingly complex as greater demands are placed on their populations. Correspondingly, needs for information to use in management have increased. Inadequate information leads the manager to err on the conservative side; unless the size of the 'harvestable surplus' is known, the population cannot be fully exploited. Conversely, information beyond what is needed becomes an unaffordable luxury. Population modeling has proven useful for organizing information on numerous game animals. Modeling serves to determine if information of the right kind and proper amount is being gathered; systematizes data collection, data interpretation, and decision making; and permits more effective management and better utilization of game populations. This report briefly reviews the principles of population modeling, describes what has been learned from previous modeling efforts on furbearers, and outlines the potential role of population modeling in furbearer management.

Book chapter↗

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

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

Techniques and Methods↗

A guide for using the transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California

This report is a guide for executing numerical simulations with the transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California using the U.S. Geological Survey modular finite-difference ground-water flow model, MODFLOW-2000. Model inputs, including observations of hydraulic head, discharge, and boundary flows, are summarized. Modification of the DVRFS transient ground-water model is discussed for two common uses of the Death Valley regional ground-water flow system model: predictive pumping scenarios that extend beyond the end of the model simulation period (1998), and model simulations with only steady-state conditions.

Open-File Report↗

Simulation of Flow, Sediment Transport, and Sediment Mobility of the Lower Coeur d'Alene River, Idaho

A one-dimensional sediment-transport model and a multi-dimensional hydraulic and bed shear stress model were developed to investigate the hydraulic, sediment transport, and sediment mobility characteristics of the lower Coeur d?Alene River in northern Idaho. This report documents the development and calibration of those models, as well as the results of model simulations. The one-dimensional sediment-transport model (HEC-6) was developed, calibrated, and used to simulate flow hydraulics and erosion, deposition, and transport of sediment in the lower Coeur d?Alene River. The HEC-6 modeled reach, comprised of 234 cross sections, extends from Enaville, Idaho, on the North Fork of the Coeur d?Alene River and near Pinehurst, Idaho, on the South Fork of the river to near Harrison, Idaho, on the main stem of the river. Bed-sediment samples collected by previous investigators and samples collected for this study in 2005 were used in the model. Sediment discharge curves from a previous study were updated using suspended-sediment samples collected at three sites since April 2000. The HEC-6 was calibrated using river discharge and water-surface elevations measured at five U.S. Geological Survey gaging stations. The calibrated HEC-6 model allowed simulation of management alternatives to assess erosion and deposition from proposed dredging of contaminated streambed sediments in the Dudley reach. Four management alternatives were simulated with HEC-6. Before the start of simulation for these alternatives, seven cross sections in the reach near Dudley, Idaho, were deepened 20 feet?removing about 296,000 cubic yards of sediments?to simulate dredging. Management alternative 1 simulated stage-discharge conditions from 2000, and alternative 2 simulated conditions from 1997. Results from alternatives 1 and 2 indicated that about 6,500 and 12,300 cubic yards, respectively, were deposited in the dredged reach. These figures represent 2 and 4 percent, respectively, of the total volume of dredged sediments removed before the start of simulation. In alternatives 3 and 4, the incoming total sediment discharges from the South Fork of the river were decreased by one-half. Management alternative 3 simulated stage-discharge conditions from 2000, and alternative 4 simulated conditions from 1997. Reducing incoming sediment discharge from the South Fork did not affect the streambed and deposition in the Dudley and downstream reaches, probably because the distance between the South Fork and the Dudley reach is long enough for sediment supply, transport capacity, and channel geometry to be balanced before reaching the Dudley and downstream reaches. Development and calibration of a multi-dimensional hydraulic and bed shear stress model (FASTMECH) allowed simulation of water-surface elevation, depth, velocity, bed shear stress, and sediment mobility in the Dudley reach (5.3 miles). The computational grid incorporated bathymetric and Light Detection and Ranging (LIDAR) data, with a node spacing of about 2.5 meters. With the exception of the fourth FASTMECH calibration simulation, results from the FASTMECH calibration simulations indicated that flow depths, flow velocities, and bed shear stresses increased as river discharge increased. Water-surface elevations in the fourth calibration simulation were about 2 feet higher than those in the other simulations because high lake levels in Coeur d?Alene Lake caused backwater conditions. Average simulated velocities along the thalweg ranged from about 3 to 5.3 feet per second, and maximum simulated velocities ranged from 3.9 to 7 feet per second. In the dredged reach, average simulated velocity along the thalweg ranged from 3.5 to 6 feet per second. The model also simulated several back-eddies (flow reversal); the largest eddy encompassed about one-third of the river width. Average bed shear stresses increased more than 200 percent from the first to the last simulation. Simulated sediment mobility, asses

Scientific Investigations Report↗

Empirical models of wind conditions on Upper Klamath Lake, Oregon

Upper Klamath Lake is a large (230 square kilometers), shallow (mean depth 2.8 meters at full pool) lake in southern Oregon. Lake circulation patterns are driven largely by wind, and the resulting currents affect the water quality and ecology of the lake. To support hydrodynamic modeling of the lake and statistical investigations of the relation between wind and lake water-quality measurements, the U.S. Geological Survey has monitored wind conditions along the lakeshore and at floating raft sites in the middle of the lake since 2005. In order to make the existing wind archive more useful, this report summarizes the development of empirical wind models that serve two purposes: (1) to fill short (on the order of hours or days) wind data gaps at raft sites in the middle of the lake, and (2) to reconstruct, on a daily basis, over periods of months to years, historical wind conditions at U.S. Geological Survey sites prior to 2005. Empirical wind models based on Artificial Neural Network (ANN) and Multivariate-Adaptive Regressive Splines (MARS) algorithms were compared. ANNs were better suited to simulating the 10-minute wind data that are the dependent variables of the gap-filling models, but the simpler MARS algorithm may be adequate to accurately simulate the daily wind data that are the dependent variables of the historical wind models. To further test the accuracy of the gap-filling models, the resulting simulated winds were used to force the hydrodynamic model of the lake, and the resulting simulated currents were compared to measurements from an acoustic Doppler current profiler. The error statistics indicated that the simulation of currents was degraded as compared to when the model was forced with observed winds, but probably is adequate for short gaps in the data of a few days or less. Transport seems to be less affected by the use of the simulated winds in place of observed winds. The simulated tracer concentration was similar between model results when simulated winds were used to force the model, and when observed winds were used to force the model, and differences between the two results did not accumulate over time.

Scientific Investigations Report↗

Application of a Geographic Information System for regridding a ground-water flow model of the Columbia Plateau Regional Aquifer System, Walla Walla River basin, Oregon-Washington

Computerized Geographic Information Systems (GIS) have become viable and valuable tools for managing,analyzing, creating, and displaying data for three-dimensional finite-difference ground-water flow models. Three GIS applications demonstrated in this study are: (1) regridding of data arrays from an existing large-area, low resolution ground-water model to a smaller, high resolution grid; (2) use of GIS techniques for assembly of data-input arrays for a ground-water model; and (3) use of GIS for rapid display of data for verification, for checking of ground-water model output, and for the cre.ation of customized maps for use in reports. The Walla Walla River Basin was selected as the location for the demonstration because (1) data from a low resolution ground-water model (Columbia Plateau Regional Aquifer System Analysis [RASA]) were available and (2) concern for long-term use of water resources for irrigation in the basin. The principal advantage of regridding is that it may provide the ability to more precisely calibrate a model, assuming chat a more detailed coverage of data is available, and to evaluate the numerical errors associated with a particular grid design.Regridding gave about an 8-fold increase in grid-node density.Several FORTRAN programs were developed to load the regridded ground-water data into a finite-difference modular model as model-compatible input files for use in a steady-state model run.To facilitate the checking and validating of the GIS regridding process, maps and tabular reports were produced for each of eight ground-water parameters by model layer. Also, an automated subroutine that was developed to view the model-calculated water levels in cross-section will aid in the synthesis and interpretation of model results.

Water-Resources Investigations Report↗

Ground water in the Fresno area, California; preliminary report

The U.S. Geological Survey, in cooperation with Fresno County, is developing a computer model of the unconfined aquifer system in the Fresno area. In the area, consolidated rocks are overlain by unconsolidated deposits. The unconsolidated deposits are divided into a lower fine-grained unit and an overlying coarse-grained unit. The system is assumed to be a single unconfined unit. The top of the system is the water table and the bottom is the base of the coarse-grained deposits. The northeastern boundary is the only natural lateral hydrologic boundary. The other lateral boundaries are placed so that errors along them will not affect heads in the area of interest. Initial heads for 1943, hydraulic conductivity, and bottom elevations were used in a steady-state model that generated net flux values. It is anticipated that most of the remaining effort will center on model calibration. (USGS)

Water-Resources Investigations Report↗

Calibration of the U.S. Geological Survey National Crustal Model

The U.S. Geological Survey National Crustal Model (NCM) is being developed to include spatially varying estimates of site response in seismic hazard assessments. Primary outputs of the NCM are continuous velocity and density profiles from the Earth’s surface to the mantle transition zone at 410-kilometer (km) depth for each location on a 1-km grid across the conterminous United States. Datasets used to produce the NCM may have a resolution of better than 1 km near the Earth’s surface in some regions, but, with increasing depth, NCM resolution decreases to tens to hundreds of kilometers in the mantle. Basic subsurface information is provided by the NCM geologic framework, thermal model, and petrologic and mineral physics database. In this report, the velocities and densities that can be extracted from the NCM are calibrated through the development of a porosity model based on Biot-Gassmann theory and more than 2,000 compressional- and (or) shear-wave velocity profiles less than 10 km deep from across the conterminous United States and southwestern Canada. Sediment and rock porosities are derived from shear-wave velocity and are found to depend on effective pressure, rock type, and age (for sedimentary and extrusive volcanic deposits). Porosity-effective pressure functions are then estimated for each rock type (and age for sedimentary and extrusive volcanic deposits). Unconsolidated sediments are found to have higher porosities than consolidated units, which have higher porosities than unweathered igneous units; young sedimentary units (for example, Quaternary age units) tend to have higher porosities than older sedimentary units (for example, pre-Cenozoic age units); porosity decreases with increasing effective pressure; and porosities can decrease quickly through the weathered layer of intrusive rocks. Comparing two Los Angeles area velocity models and the U.S. Geological Survey Bay Area velocity model with the NCM, the NCM does a better job on average of reproducing observed shear-wave velocities below 1 km per second because it has less bias and uncertainty. Approaching and above 1 km per second, the NCM tends to underpredict observed shear-wave velocity. Whereas several factors could contribute to this, the primary factor is probably bias in the NCM geologic framework. For example, the NCM will predict lower velocities in places where the depth to bedrock and basement appear shallower in the measured velocity profiles than specified in the NCM geologic framework. With regard to observed compressional-wave velocity and density, the NCM has significantly less bias than California models for the former, especially below 2 km per second, and all models tend to overpredict density for densities less than about 2,200 kilograms per cubic meter.

Open-File Report↗