Search USGS⌕ Search

Geology topics

John Doherty

Publications and source records attributed to John Doherty.

9 recordsLinked to original sources
USGS2012

Approaches in highly parameterized inversion: TSPROC, a general time-series processor to assist in model calibration and result summarization

The TSPROC ( T ime S eries PROC essor) computer software uses a simple scripting language to process and analyze time series. It was developed primarily to assist in the calibration of environmental models. The software is designed to perform calculations on time-series data commonly associated with surface-water models, including calculation of flow volumes, transformation by means of basic arithmetic operations, and generation of seasonal and annual statistics and hydrologic indices. TSPROC can also be used to generate some of the key input files required to perform parameter optimization by means of the PEST ( P arameter EST imation) computer software. Through the use of TSPROC, the objective function for use in the model-calibration process can be focused on specific components of a hydrograph.

Techniques and Methods↗
USGS2011

Predictive uncertainty analysis of a saltwater intrusion model using null-space Monte Carlo

Because of the extensive computational burden and perhaps a lack of awareness of existing methods, rigorous uncertainty analyses are rarely conducted for variable-density flow and transport models. For this reason, a recently developed null-space Monte Carlo (NSMC) method for quantifying prediction uncertainty was tested for a synthetic saltwater intrusion model patterned after the Henry problem. Saltwater intrusion caused by a reduction in fresh groundwater discharge was simulated for 1000 randomly generated hydraulic conductivity distributions, representing a mildly heterogeneous aquifer. From these 1000 simulations, the hydraulic conductivity distribution giving rise to the most extreme case of saltwater intrusion was selected and was assumed to represent the "true" system. Head and salinity values from this true model were then extracted and used as observations for subsequent model calibration. Random noise was added to the observations to approximate realistic field conditions. The NSMC method was used to calculate 1000 calibration-constrained parameter fields. If the dimensionality of the solution space was set appropriately, the estimated uncertainty range from the NSMC analysis encompassed the truth. Several variants of the method were implemented to investigate their effect on the efficiency of the NSMC method. Reducing the dimensionality of the null-space for the processing of the random parameter sets did not result in any significant gains in efficiency and compromised the ability of the NSMC method to encompass the true prediction value. The addition of intrapilot point heterogeneity to the NSMC process was also tested. According to a variogram comparison, this provided the same scale of heterogeneity that was used to generate the truth. However, incorporation of intrapilot point variability did not make a noticeable difference to the uncertainty of the prediction. With this higher level of heterogeneity, however, the computational burden of generating calibration-constrained parameter fields approximately doubled. Predictive uncertainty variance computed through the NSMC method was compared with that computed through linear analysis. The results were in good agreement, with the NSMC method estimate showing a slightly smaller range of prediction uncertainty than was calculated by the linear method. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2010

Using a cloud to replenish parched groundwater modeling efforts

Groundwater models can be improved by introduction of additional parameter flexibility and simultaneous use of soft-knowledge. However, these sophisticated approaches have high computational requirements. Cloud computing provides unprecedented access to computing power via the Internet to facilitate the use of these techniques. A modeler can create, launch, and terminate “virtual” computers as needed, paying by the hour, and save machine images for future use. Such cost-effective and flexible computing power empowers groundwater modelers to routinely perform model calibration and uncertainty analysis in ways not previously possible.

Ground Water↗
USGS2010

Creative use of pilot points to address site and regional scale heterogeneity in a variable-density model

Pilot points for parameter estimation were creatively used to address heterogeneity at both the well field and regional scales in a variable-density groundwater flow and solute transport model designed to test multiple hypotheses for upward migration of fresh effluent injected into a highly transmissive saline carbonate aquifer. Two sets of pilot points were used within in multiple model layers, with one set of inner pilot points (totaling 158) having high spatial density to represent hydraulic conductivity at the site, while a second set of outer points (totaling 36) of lower spatial density was used to represent hydraulic conductivity further from the site. Use of a lower spatial density outside the site allowed (1) the total number of pilot points to be reduced while maintaining flexibility to accommodate heterogeneity at different scales, and (2) development of a model with greater areal extent in order to simulate proper boundary conditions that have a limited effect on the area of interest. The parameters associated with the inner pilot points were log transformed hydraulic conductivity multipliers of the conductivity field obtained by interpolation from outer pilot points. The use of this dual inner-outer scale parameterization (with inner parameters constituting multipliers for outer parameters) allowed smooth transition of hydraulic conductivity from the site scale, where greater spatial variability of hydraulic properties exists, to the regional scale where less spatial variability was necessary for model calibration. While the model is highly parameterized to accommodate potential aquifer heterogeneity, the total number of pilot points is kept at a minimum to enable reasonable calibration run times.

Book↗
USGS2009

Using a coupled groundwater/surface-water model to predict climate-change impacts to lakes in the Trout Lake Watershed, northern Wisconsin

A major focus of the U.S. Geological Survey’s Trout Lake Water, Energy, and Biogeochemical Budgets (WEBB) project is the development of a watershed model to allow predictions of hydrologic response to future conditions including land-use and climate change. The coupled groundwater/surface-water model GSFLOW was chosen for this purpose because it could easily incorporate an existing groundwater flow model and it provides for simulation of surface-water processes. The Trout Lake watershed in northern Wisconsin is underlain by a highly conductive outwash sand aquifer. In this area, streamflow is dominated by groudwater contributions, however, surface runoff occurs during intense rainfall periods and spring snowmelt. Surface runoff also occurs locally near stream/lake areas where the unsaturated zone is thin. A diverse data set, collected from 1992 to 2007 for the Trout Lake WEBB project and the co-located and NSF-funded North Temperate Lake LTER project, includes snowpack, solar radiation, potential evapotranspiration, lake levels, groundwater levels, and streamflow. The time-series processing software TSPROC (Doherty 2001)was used to distill the large time series data set to a smaller set of observations and summary statistics that captured the salient hydrologic information. The time-series processing reduced hundreds of thousands of observations to less than 5,000. Model calibration included specific predictions for several lakes in the study area using the PEST parameter estimation suit of software (Doherty 2007). The calibrated model was used to simulate the hydrologic response in the study lakes to a variety of climate change scenarios culled from the IPCC Fourth Assessment Report of the Intergovernmental Panel of Climate Change (Solomon et al. 2007). Results from the simulations indicate climate change could result in substantial changes to the lake levels and components of the hydrologic budget of a seepage lake in the flow system. For a drainage lake lower in the flow system, the impacts of climate change are diminished.

Wisconsin↗
USGS2009

Using a coupled groundwater/surfacewater model to predict climate-change impacts to lakes in the Trout Lake watershed, Northern Wisconsin

A major focus of the U.S. Geological Survey’s Trout Lake Water, Energy, and Biogeochemical Budgets (WEBB) project is the development of a watershed model to allow predictions of hydrologic response to future conditions including land-use and climate change. The coupled groundwater/surface-water model GSFLOW was chosen for this purpose because it could easily incorporate an existing groundwater flow model and it provides for simulation of surface-water processes. The Trout Lake watershed in northern Wisconsin is underlain by a highly conductive outwash sand aquifer. In this area, streamflow is dominated by groundwater contributions; however, surface runoff occurs during intense rainfall periods and spring snowmelt. Surface runoff also occurs locally near stream/lake areas where the unsaturated zone is thin. A diverse data set, collected from 1992 to 2007 for the Trout Lake WEBB project and the co-located and NSF-funded North Temperate Lakes LTER project, includes snowpack, solar radiation, potential evapotranspiration, lake levels, groundwater levels, and streamflow. The timeseries processing software TSPROC (Doherty 2003) was used to distill the large time series data set to a smaller set of observations and summary statistics that captured the salient hydrologic information. The timeseries processing reduced hundreds of thousands of observations to less than 5,000. Model calibration included specific predictions for several lakes in the study area using the PEST parameter estimation suite of software (Doherty 2007). The calibrated model was used to simulate the hydrologic response in the study lakes to a variety of climate change scenarios culled from the IPCC Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Solomon et al. 2007). Results from the simulations indicate climate change could result in substantial changes to the lake levels and components of the hydrologic budget of a seepage lake in the flow system. For a drainage lake lower in the flow system, the impacts of climate change are diminished.

Conference Paper↗
USGS2003

MODFLOW-2000, the U.S. Geological Survey modular ground-water model -- Documentation of MOD-PREDICT for predictions, prediction sensitivity analysis, and evaluation of uncertainty

This document describes the MOD-PREDICT program, which helps evaluate userdefined sets of observations, prior information, and predictions, using the ground-water model MODFLOW-2000. MOD-PREDICT takes advantage of the existing Observation and Sensitivity Processes (Hill and others, 2000) by initiating runs of MODFLOW-2000 and using the output files produced. The names and formats of the MODFLOW-2000 input files are unchanged, such that full backward compatibility is maintained. A new name file and input files are required for MOD-PREDICT. The performance of MOD-PREDICT has been tested in a variety of applications. Future applications, however, might reveal errors that were not detected in the test simulations. Users are requested to notify the U.S. Geological Survey of any errors found in this document or the computer program using the email address available at the web address below. Updates might occasionally be made to this document, to the MOD-PREDICT program, and to MODFLOW- 2000. Users can check for updates on the Internet at URL http://water.usgs.gov/software/ground water.html/.

Open-File Report↗