Geology ReportsSearch

SEARCH · Geology Reports

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

Research on interactive genetic-geological models to evaluate favourability for undiscovered uranium resources

Current methods of evaluating favourability for undiscovered uranium resources are unduly subjective, quite possibly inconsistent and, as a consequence, of questionable reliability. This research is aimed at reducing the subjectivity and increasing the reliability by designing an improved method that depends largely on geological data and their statistical frequency of occurrence. This progress report outlines a genetic approach to modelling the geological factors that controlled uranium mineralization in order to evaluate the favourability for the occurrence of undiscovered uranium deposits of the type modeled. A genetic model is constructed from all the factors that describe the processes, in chronological sequence, that formed uranium deposits thought to have a common origin. The field and laboratory evidence for the processes constitute a geologic occurrence base that parallels the chronological sequence of events. The genetic model and the geologic-occurrence base are portrayed as two columns of an interactive matrix called the "genetic-geologic model". For each column, eight chronological stages are used to describe the overall formation of the uranium deposits. These stages consist of (1) precursor processes; (2) host-rock formation; (3) preparation of host-rock; (4) uranium-source development; (5) transport of uranium; (6) primary uranium deposition; (7) post-deposition modification; and (8) preservation. To apply the genetic-geological model to evaluate favourability, a question is posed that determines the presence or absence of each attribute listed under the geologic-occurrence base. By building a logic circuit of the attributes according to either their essential or non-essential nature, the resultant match between a well-documented control area and the test area may be determined. The degree of match is a measure of favourability for uranium occurrence as hypothesized in the genetic model. This process of geological decision analysis results in a series of favourability maps that can be combined into a final composite favourability map.

Conference Paper

Progress report on daily flow-routing simulation for the Carson River, California and Nevada

A physically based flow-routing model using Hydrological Simulation Program-FORTRAN (HSPF) was constructed for modeling streamflow in the Carson River at daily time intervals as part of the Truckee-Carson Program of the U.S. Geological Survey (USGS). Daily streamflow data for water years 1978-92 for the mainstem river, tributaries, and irrigation ditches from the East Fork Carson River near Markleeville and West Fork Carson River at Woodfords down to the mainstem Carson River at Fort Churchill upstream from Lahontan Reservoir were obtained from several agencies and were compiled into a comprehensive data base. No previous physically based flow-routing model of the Carson River has incorporated multi-agency streamflow data into a single data base and simulated flow at a daily time interval. Where streamflow data were unavailable or incomplete, hydrologic techniques were used to estimate some flows. For modeling purposes, the Carson River was divided into six segments, which correspond to those used in the Alpine Decree that governs water rights along the river. Hydraulic characteristics were defined for 48 individual stream reaches based on cross-sectional survey data obtained from field surveys and previous studies. Simulation results from the model were compared with available observed and estimated streamflow data. Model testing demonstrated that hydraulic characteristics of the Carson River are adequately represented in the models for a range of flow regimes. Differences between simulated and observed streamflow result mostly from inadequate data characterizing inflow and outflow from the river. Because irrigation return flows are largely unknown, irrigation return flow percentages were used as a calibration parameter to minimize differences between observed and simulated streamflows. Observed and simulated streamflow were compared for daily periods for the full modeled length of the Carson River and for two major subreaches modeled with more detailed input data. Hydrographs and statistics presented in this report describe these differences. A sensitivity analysis of four estimated components of the hydrologic system evaluated which components were significant in the model. Estimated ungaged tributary streamflow is not a significant component of the model during low runoff, but is significant during high runoff. The sensitivity analysis indicates that changes in the estimated irrigation diversion and estimated return flow creates a noticeable change in the statistics. The modeling for this study is preliminary. Results of the model are constrained by current availability and accuracy of observed hydrologic data. Several inflows and outflows of the Carson River are not described by time-series data and therefore are not represented in the model.

Open-File Report

Random forest regression models for estimating low-streamflow statistics at ungaged locations in New York, excluding Long Island

This report is currently being updated with re‑calibrated random‑forest models for StreamStats implementation. The basin‑characteristic inputs were modified to align with StreamStats requirements. Abstract Models to estimate low-streamflow statistics at ungaged locations in New York, excluding Long Island and including hydrologically connected basins from bordering States, were developed for the first time by the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation. A total of 224 basin characteristics were developed for 213 unaltered streamgages (locations where the human effects on streamflow were limited), across the following categories: basin geometry, climate, land cover, soils, surficial geology, and other characteristics. The basins with unaltered streamgages were evaluated for potential redundancy, and streamgages in close proximity and with similar drainage areas were flagged and removed from the testing and cross-validation datasets to prevent data leaking from the training dataset to the testing dataset. Random forest regression models were created by using basin characteristics as predictor variables and by developing a workflow to train, tune, and test the model. Models were developed to estimate the ungaged lowest annual 7-day and 30-day average streamflow that occurs (on average) once every 10 years (7Q10 and 30Q10). The top four basin characteristics used for the 7Q10 and 30Q10 models were drainage area, total stream length, perimeter of the basin, and length of the longest flow path. Results for the 7Q10 and 30Q10 models had coefficients of determination ( R 2 ) of 0.796 and 0.853, respectively. The output model results were bias-corrected for ungaged locations across New York and are available within the interactive StreamStats tool.

New York

Effects of Land-Use Changes and Ground-Water Withdrawals on Stream Base Flow, Pocono Creek Watershed, Monroe County, Pennsylvania

The Pocono Creek watershed drains 46.5 square miles in eastern Monroe County, Pa. Between 2000 and 2020, the population of Monroe County is expected to increase by 70 percent, which will result in substantial changes in land-use patterns. An evaluation of the effect of reduced recharge from land-use changes and additional ground-water withdrawals on stream base flow was done by the U.S. Geological Survey (USGS) in cooperation with the U.S. Environmental Protection Agency (USEPA) and the Delaware River Basin Commission as part of the USEPA?s Framework for Sustainable Watershed Management Initiative. Two models were used. A Soil and Water Assessment Tool (SWAT) model developed by the USEPA provided areal recharge values for 2000 land use and projected full buildout land use. The USGS MODFLOW-2000 ground-water-flow model was used to estimate the effect of reduced recharge from changes in land use and additional ground-water withdrawals on stream base flow. This report describes the ground-water-flow-model simulations. The Pocono Creek watershed is underlain by sedimentary rock of Devonian age, which is overlain by a veneer of glacial deposits. All water-supply wells are cased into and derive water from the bedrock. In the ground-water-flow model, the surficial geologic units were grouped into six categories: (1) moraine deposits, (2) stratified drift, (3) lake deposits, (4) outwash, (5) swamp deposits, and (6) undifferentiated deposits. The unconsolidated surficial deposits are not used as a source of water. The ground-water and surface-water systems are well connected in the Pocono Creek watershed. Base flow measured on October 13, 2004, at 27 sites for model calibration showed that streams gained water between all sites measured except in the lower reach of Pocono Creek. The ground-water-flow model included the entire Pocono Creek watershed. Horizontally, the modeled area was divided into a 53 by 155 cell grid with 6,060 active cells. Vertically, the modeled area was discretized into four layers. Layers 1 and 2 represented the unconsolidated surficial deposits where they are present and bedrock where the surficial deposits are absent. Layer 3 represented shallow bedrock and was 200 ft (feet) thick. Layer 4 represented deep bedrock and was 300 ft thick. A total of 873 cells representing streams were assigned to layer 1. Recharge rates for model calibration were provided by the USEPA SWAT model for 2000 land-use conditions. Recharge rates for 2000 for the 29 subwatersheds in the SWAT model ranged from 6.11 to 22.66 inches per year. Because the ground-water-flow model was calibrated to base-flow data collected on October 13, 2004, the 2000 recharge rates were multiplied by 1.18 so the volume of recharge was equal to the volume of streamflow measured at the mouth of Pocono Creek. During model calibration, adjustments were made to aquifer hydraulic conductivity and streambed conductance. Simulated base flows and hydraulic heads were compared to measured base flows and hydraulic heads using the root mean squared error (RMSE) between measured and simulated values. The RMSE of the calibrated model for base flow was 4.7 cubic feet per second for 27 locations, and the RMSE for hydraulic heads for 15 locations was 35 ft. The USEPA SWAT model was used to provide areal recharge values for 2000 and full buildout land-use conditions. The change in recharge ranged from an increase of 37.8 percent to a decrease of 60.8 percent. The ground-water-flow model was used to simulate base flow for 2000 and full buildout land-use conditions using steady-state simulations. The decrease in simulated base flow ranged from 3.8 to 63 percent at the streamflow-measurement sites. Simulated base flow at streamflow-gaging station Pocono Creek above Wigwam Run near Stroudsburg, Pa. (01441495), decreased 25 percent. This is in general agreement with the SWAT model, which estimated a 30.6-percent loss in base flow at the streamflow-gaging station.

Scientific Investigations Report

A theory for modeling ground-water flow in heterogeneous media

Construction of a ground-water model for a field area is not a straightforward process. Data are virtually never complete or detailed enough to allow substitution into the model equations and direct computation of the results of interest. Formal model calibration through optimization, statistical, and geostatistical methods is being applied to an increasing extent to deal with this problem and provide for quantitative evaluation and uncertainty analysis of the model. However, these approaches are hampered by two pervasive problems: 1) nonlinearity of the solution of the model equations with respect to some of the model (or hydrogeologic) input variables (termed in this report system characteristics) and 2) detailed and generally unknown spatial variability (heterogeneity) of some of the system characteristics such as log hydraulic conductivity, specific storage, recharge and discharge, and boundary conditions. A theory is developed in this report to address these problems. The theory allows construction and analysis of a ground-water model of flow (and, by extension, transport) in heterogeneous media using a small number of lumped or smoothed system characteristics (termed parameters). The theory fully addresses both nonlinearity and heterogeneity in such a way that the parameters are not assumed to be effective values. The ground-water flow system is assumed to be adequately characterized by a set of spatially and temporally distributed discrete values, ?, of the system characteristics. This set contains both small-scale variability that cannot be described in a model and large-scale variability that can. The spatial and temporal variability in ? are accounted for by imagining ? to be generated by a stochastic process wherein ? is normally distributed, although normality is not essential. Because ? has too large a dimension to be estimated using the data normally available, for modeling purposes ? is replaced by a smoothed or lumped approximation y?. (where y is a spatial and temporal interpolation matrix). Set y?. has the same form as the expected value of ?, y 'line' ? , where 'line' ? is the set of drift parameters of the stochastic process; ?. is a best-fit vector to ?. A model function f(?), such as a computed hydraulic head or flux, is assumed to accurately represent an actual field quantity, but the same function written using y?., f(y?.), contains error from lumping or smoothing of ? using y?.. Thus, the replacement of ? by y?. yields nonzero mean model errors of the form E(f(?)-f(y?.)) throughout the model and covariances between model errors at points throughout the model. These nonzero means and covariances are evaluated through third and fifth-order accuracy, respectively, using Taylor series expansions. They can have a significant effect on construction and interpretation of a model that is calibrated by estimating ?.. Vector ?.. is estimated as 'hat' ? using weighted nonlinear least squares techniques to fit a set of model functions f(y'hat' ?) to a. corresponding set of observations of f(?), Y. These observations are assumed to be corrupted by zero-mean, normally distributed observation errors, although, as for ?, normality is not essential. An analytical approximation of the nonlinear least squares solution is obtained using Taylor series expansions and perturbation techniques that assume model and observation errors to be small. This solution is used to evaluate biases and other results to second-order accuracy in the errors. The correct weight matrix to use in the analysis is shown to be the inverse of the second-moment matrix E(Y-f(y?.))(Y-f(y?.))', but the weight matrix is assumed to be arbitrary in most developments. The best diagonal approximation is the inverse of the matrix of diagonal elements of E(Y-f(y?.))(Y-f(y?.))', and a method of estimating this diagonal matrix when it is unknown is developed using a special objective function to compute 'hat' ?. When considered to be an estimate of f

Professional Paper

Estimation of groundwater use for a groundwater-flow model of the Lake Michigan Basin and adjacent areas, 1864-2005

The U.S. Geological Survey, at the request of Congress, is assessing the availability and use of the Nation's water resources to help characterize how much water is available now, how water availability is changing, and how much water can be expected to be available in the future. The Great Lakes Basin Pilot project of the U.S. Geological Survey national assessment of water availability and use focused on the Great Lakes Basin and included detailed studies of the processes governing water availability in the Great Lakes Basin. One of these studies included the development of a groundwater-flow model of the Lake Michigan Basin. This report describes the compilation and estimation of the groundwater withdrawals in those areas in Wisconsin, Michigan, Indiana, and Illinois that were needed for the Lake Michigan Basin study groundwater-flow model. These data were aggregated for 12 model time intervals spanning 1864 to 2005 and were summarized by model area, model subregion, category of water use, aquifer system, aquifer type, and hydrogeologic unit model layer. The types and availability of information on groundwater withdrawals vary considerably among states because water-use programs often differ in the types of data collected and in the methods and frequency of data collection. As a consequence, the methods used to estimate and verify the data also vary. Additionally, because of the different sources of data and different terminologies applied for the purposes of this report, the water-use data published in this report may differ from water-use data presented in other reports. These data represent only a partial estimate of groundwater use in each state because estimates were compiled only for areas in Wisconsin, Michigan, Indiana, and Illinois within the Lake Michigan Basin model area. Groundwater-withdrawal data were compiled for both nearfield and farfield model areas in Wisconsin and Illinois, whereas these data were compiled primarily for the nearfield model area in Michigan and Indiana. Overall water use for the selected areas in Wisconsin, Michigan, Indiana, and Illinois was less during early time intervals than during more recent intervals, with large increases beginning around the 1960s. Total estimated groundwater withdrawals for model input range from 18.01 million gallons per day (Mgal/d) for interval 1 (1864-1900) to 1,280.25 Mgal/d for interval 12 (2001-5). Withdrawals for the public-supply category make up the majority of the withdrawals in each of the four states. In Wisconsin and Michigan, the second largest withdrawals are for the irrigation category; in Indiana and Illinois, industrial withdrawals account for the second largest withdrawal amounts. The smallest withdrawals are for miscellaneous uses in Wisconsin and irrigation uses in Indiana and Illinois. Estimated groundwater withdrawals in the Southern Lower Peninsula of Michigan, Northeastern Illinois, and the farfield model area are generally larger than in the other model subregions. Withdrawals in Michigan and Indiana are predominantly from the Quaternary aquifer system, whereas withdrawals in Illinois are predominantly from the Cambrian-Ordovician aquifer systems. Withdrawals in Wisconsin are about equal from the Quaternary and Cambrian-Ordovician aquifer systems. Estimated groundwater withdrawals in Michigan and Indiana are predominantly from the unconfined unconsolidated aquifer type. Withdrawals in Illinois are largely from the deep confined bedrock aquifer type, although they decreased considerably in more recent time intervals. Wisconsin withdrawals are about equal from unconfined unconsolidated and deep confined bedrock aquifer types. Groundwater-withdrawal estimates in Wisconsin were compiled for the 47 easternmost counties within the boundary of the Lake Michigan Basin model, of which 32 counties, though not entirely contained, are at least partly within the Lake Michigan Basin. Overall, 6,457 withdrawal locations were estima

Scientific Investigations Report

A model for evaluating stream temperature response to climate change in Wisconsin

Expected climatic changes in air temperature and precipitation patterns across the State of Wisconsin may alter future stream temperature and flow regimes. As a consequence of flow and temperature changes, the composition and distribution of fish species assemblages are expected to change. In an effort to gain a better understanding of how climatic changes may affect stream temperature, an approach was developed to predict and project daily summertime stream temperature under current and future climate conditions for 94,341 stream kilometers across Wisconsin. The approach uses a combination of static landscape characteristics and dynamic time-series climatic variables as input for an Artificial Neural Network (ANN) Model integrated with a Soil-Water-Balance (SWB) Model. Future climate scenarios are based on output from downscaled General Circulation Models (GCMs). The SWB model provided a means to estimate the temporal variability in groundwater recharge and provided a mechanism to evaluate the effect of changing air temperature and precipitation on groundwater recharge and soil moisture. The Integrated Soil-Water-Balance and Artificial Neural Network version 1 (SWB-ANNv1) Model was used to simulate daily summertime stream temperature under current (1990–2008) climate and explained 76 percent of the variation in the daily mean based on validation at 67 independent sites. Results were summarized as July mean water temperature, and individual stream segments were classified by thermal class (cold, cold transition, warm transition, and warm) for comparison of current (1990–2008) with future climate conditions. Integrating the SWB Model with the ANN Model provided a mechanism by which downscaled global or regional climate model results could be used to estimate the potential effects of climate change on future stream temperature on a daily time step. To address future climate scenarios, statistically downscaled air temperature and precipitation projections from 10 GCMs and 2 time periods were used with the SWB-ANNv1 Model to project future stream temperature. Projections of future stream temperatures at mid- (2046–65) and late- (2081–2100) 21st century showed the July mean water temperature increasing for all stream segments with about 80 percent of stream kilometers increasing by 1 to 2 degrees Celsius (°C) by mid-century and about 99 percent increasing by 1 to 3 °C by late-century. Projected changes in stream temperatures also affected changes in thermal classes with a loss in the total amount of cold-water, cold-transition, and warm-transition thermal habitat and a gain in warm-water and very warm thermal habitat for both mid- and late-21st century time periods. The greatest losses occurred for cold-water streams and the greatest gains for warm-water streams, with a contraction of cold-water streams in the Driftless Area of western and southern Wisconsin and an expansion of warm-water streams across northern Wisconsin. Results of this study suggest that such changes will affect the composition of fish assemblages, with a loss of suitable habitat for cold-water fishes and gain in suitable habitat for warm-water fishes. In the end, these projected changes in thermal habitat attributable to climate may result in a net loss of fisheries, because many warm-water species may be unable to colonize habitats formerly occupied by cold-water species because of other habitat limitations (e.g., stream size, gradient). Although projected stream temperatures may vary greatly, depending on the emissions scenario and models used, the results presented in this report represent one possibility. The relative change in stream temperature can provide useful information for planning for potential climate impacts to aquatic ecosystems. Model results can be used to help identify vulnerabilities of streams to climate change, guide stream surveys and thermal classifications, prioritize the allocation of scarce financial resources, identify approaches to climate adaptation to best protect and enhance resiliency in stream thermal habitat, and provide information to make quantitative assessments of statewide stream resources.

Wisconsin

Audiomagnetotelluric data and preliminary two-dimensional models from Spring, Dry Lake, and Delamar Valleys, Nevada

This report presents audiomagnetotelluric (AMT) data along fourteen profiles in Spring, Delamar, and Dry Lake Valleys, and the corresponding preliminary two-dimensional (2-D) inverse models. The AMT method is a valuable tool for estimating the electrical resistivity of the Earth over depth ranges from a few meters to less than one kilometer, and it is important for revealing subsurface structure and stratigraphy within the Basin and Range province of eastern Nevada, which can be used to define the geohydrologic framework of the region. We collected AMT data by using the Geometrics StrataGem EH4 system. Profiles were 0.7 - 3.2 km in length with station spacing of 50-400 m. Data were recorded in a coordinate system parallel to and perpendicular to the regional geologic-strike direction with Z positive down. We show AMT station locations, sounding curves of apparent resistivity, phase, and coherency, and 2-D models of subsurface resistivity along the profiles. The 2-D inverse models are computed from the transverse electric (TE), transverse magnetic (TM), and TE+TM mode data by using a conjugate gradient, finite-difference method. Preliminary interpretation of the 2-D models defines the structural framework of the basins and the resistivity contrasts between alluvial basin-fill, volcanic units, and carbonate basement rocks.

Nevada

Hydrogeology and simulation of ground-water flow at Dover Air Force Base, Delaware

Dover Air Force Base in Kent County, Delaware, has many contaminated sites that are in active remediation. To assist in this remediation, a steady-state model of ground-water flow was developed to aid in understanding the hydrology of the system, and for use as a ground-watermanagement tool. This report describes the hydrology on which the model is based, a description of the model itself, and some applications of the model.Dover Air Force Base is underlain by unconsolidated sediments of the Atlantic Coastal Plain. The primary units that were investigated include the upper Calvert Formation and the overlying Columbia Formation. The uppermost sand unit in the Calvert Formation at Dover Air Force Base is the Frederica aquifer, which is the deepest unit investigated in this report. A confining unit of clayey silt in the upper Calvert Formation separates the Frederica aquifer from the lower surficial aquifer, which is the basal Columbia Formation. North and northwest of Dover Air Force Base, the Frederica aquifer subcrops beneath the Columbia Formation and the upper Calvert Formation confining unit is absent. The Calvert Formation dips to the southeast. The Columbia Formation consists predominately of sands, silts, and gravels, although in places there are clay layers that separate the surficial aquifer into an upper and lower surficial aquifer. The areal extent of these clay layers has been mapped by use of gamma logs. Long-term hydrographs reveal substantial changes in both seasonal and annual ground-water recharge. These variations in recharge are related to temporal changes in evaporation, transpiration, and precipitation. The hydrographs show areas where extensive silts and clays are present in the surficial aquifer. In these areas, the vertical gradient between water levels in wells screened above and below the clays can be as large as several feet, and local ground-water highs typically form during normal recharge conditions. When drought conditions persist, water drains off these highs and the vertical gradients decrease. At the south end of Dover Air Force Base, hydrographs of water levels in the Frederica aquifer show that off-Base pumping can cause the water levels to decline below sea level during part of the year.A 4-layer, steady-state numerical model of ground-water flow was developed for Dover Air Force Base and the surrounding area. The upper two layers represent the upper and lower surficial aquifers, which are in the Columbia Formation. In some areas of the model, a semi-confining unit is used to represent an intermittent clay layer between the upper and lower surficial aquifer. This semi-confining unit causes the local groundwater highs in the surficial aquifer. The third model layer represents the upper part of the Calvert Formation, a confining unit. The fourth model layer represents the Frederica aquifer. The model was calibrated to hydraulic heads and to ground-water discharge in Pipe Elm Branch, both of which were measured in September 1997. For the calibrated model, the root-mean-squared errors for the hydraulic heads and the ground-water discharge in the Pipe Elm Branch were 9 percent of the range of head and 3 percent of discharge, respectively. Heads simulated by use of the model were consistent with a map showing average water levels in the region. The U.S. Geological Survey?s MODPATH program was used to simulate ground-water-flow directions for several areas on the Base. This analysis showed the effects of the local groundwater highs. In these areas, ground water can flow from the highs and then dramatically change flow direction as it enters the lower surficial aquifer. The steady-state model has several limitations. The entire ground-water system is under transient hydraulic conditions, due mainly to seasonal and yearly changes in recharge and to withdrawal from irrigation wells. Yet this steady-state model is still considered to be an effective tool for understanding the ground-water-flow system u

Water-Resources Investigations Report

Documentation of a groundwater flow model (SJRRPGW) for the San Joaquin River Restoration Program study area, California

To better understand the potential effects of restoration flows on existing drainage problems, anticipated as a result of the San Joaquin River Restoration Program (SJRRP), the U.S. Geological Survey (USGS), in cooperation with the U.S. Bureau of Reclamation (Reclamation), developed a groundwater flow model (SJRRPGW) of the SJRRP study area that is within 5 miles of the San Joaquin River and adjacent bypass system from Friant Dam to the Merced River. The primary goal of the SJRRP is to reestablish the natural ecology of the river to a degree that restores salmon and other fish populations. Increased flows in the river, particularly during the spring salmon run, are a key component of the restoration effort. A potential consequence of these increased river flows is the exacerbation of existing irrigation drainage problems along a section of the river between Mendota and the confluence with the Merced River. Historically, this reach typically was underlain by a water table within 10 feet of the land surface, thus requiring careful irrigation management and (or) artificial drainage to maintain crop health. The SJRRPGW is designed to meet the short-term needs of the SJRRP; future versions of the model may incorporate potential enhancements, several of which are identified in this report. The SJRRPGW was constructed using the USGS groundwater flow model MODFLOW and was built on the framework of the USGS Central Valley Hydrologic Model (CVHM) within which the SJRRPGW model domain is embedded. The Farm Process (FMP2) was used to simulate the supply and demand components of irrigated agriculture. The Streamflow-Routing Package (SFR2) was used to simulate the streams and bypasses and their interaction with the aquifer system. The 1,300-square mile study area was subdivided into 0.25-mile by 0.25-mile cells. The sediment texture of the aquifer system, which was used to distribute hydraulic properties by model cell, was refined from that used in the CVHM to better represent the natural heterogeneity of aquifer-system materials within the model domain. In addition, the stream properties were updated from the CVHM to better simulate stream-aquifer interactions, and water-budget subregions were refined to better simulate agricultural water supply and demand. External boundary conditions were derived from the CVHM. The SJRRPGW was calibrated for April 1961 to September 2003 by using groundwater-level observations from 133 wells and streamflow observations from 19 streamgages. The model was calibrated using public-domain parameter estimation software (PEST) in a semi-automated manner. The simulated groundwater-level elevations and trends (including seasonal fluctuations) and surface-water flow magnitudes and trends reasonably matched observed data. The calibrated model is planned to be used to assess the potential effects of restoration flows on agricultural lands and the relative capabilities of proposed SJRRP actions to reduce these effects.

California

Mount Rainier National Park and Olympic National Park elk monitoring program annual report 2011

Fiscal year 2011 was the first year of implementing an approved elk monitoring protocol in Mount Rainier (MORA) and Olympic (OLYM) National Parks in the North Coast and Cascades Network (NCCN) (Griffin et al. 2012). However, it was the fourth and second year of gathering data according to protocol in MORA and OLYM respectively; data gathered during the protocol development phase followed procedures that are laid out in the protocol. Elk monitoring in these large wilderness parks relies on aerial surveys from a helicopter. Summer surveys are intended to provide quantitative estimates of abundance, sex and age composition, and distribution of migratory elk in high elevation trend count areas. An unknown number of elk is not detected during surveys; however the protocol estimates the number of missed elk by applying a model that accounts for detection bias. Detection bias in elk surveys in MORA is estimated using a double-observer sightability model that was developed using survey data from 2008-2010 (Griffin et al. 2012). That model was developed using elk that were previously equipped with radio collars by cooperating tribes. At the onset of protocol development in OLYM there were no existing radio-collars on elk. Consequently the majority of the effort in OLYM in the past 4 years has been focused on capturing and radio-collaring elk and conducting sightability trials needed to develop a double-observer sightability model in OLYM. In this annual report we provide estimates of abundance and composition for MORA elk, raw counts of elk made in OLYM, and describe sightability trials conducted in OLYM. At MORA the North trend count area was surveyed twice and the South once (North Rainier herd, and South Rainier herd). We counted 373 and 267 elk during two replicate surveys of the North Rainier herd, and 535 elk in the South Rainier herd. Using the model, we estimated that 413 and 320 elk were in the North and 652 elk were in the South trend count areas during the time of the respective surveys. At OLYM, the Core and Northwest trend count areas were completely surveyed, as were portions of the Quinault. In addition, we surveyed 10 survey units specifically to get resight data. Two-hundred and forty eight elk were counted in the Core, 19 in the Northwest, and 169 in the Quinault. We conducted double-observer sightability trials associated with 14 collared elk groups for use in developing the double-observer sightability model for OLYM.

Washington

Assessment of Effectiveness and Limitations of Habitat Suitability Models for Wetland Restoration

Habitat suitability index (HSI) models developed for wildlife in the Louisiana Coastal Area Comprehensive Ecosystem Restoration Plan (LCA study) have been assessed for parameter and overall model quality. The success of the suitability models from the South Florida Water Management District for The Everglades restoration project and from the Spatially Explicit Species Index Models (SESI) of the Across Trophic Level System Simulation (ATLSS) Program of Florida warranted investigation with possible application of modeling theory to the current LCA study. General HSI models developed by the U.S. Fish and Wildlife Service were also investigated. This report presents examinations of theoretical formulae and comparisons of the models, performed by using diverse hypothetical settings of hydrological/biological ecosystems to highlight weaknesses as well as strengths among the models, limited to the American alligator and selected wading bird species (great blue heron, great egret, and white ibis). Recommendations were made for the LCA study based on these assessments. An enhanced HSI model for the LCA study is proposed for the American alligator, and a new HSI model for wading birds is introduced for the LCA study. Performance comparisons of the proposed models with the other suitability models are made by using the aforementioned hypothetical settings.

Open-File Report

Peak streamflow trends in Illinois and their relation to changes in climate, water years 1921–2020

This report characterizes changes in peak streamflow in Illinois and the relation of these changes to climatic variability, and provides a foundation for future studies that can address nonstationarity in peak-flow frequency analysis in Illinois. Records of annual peak and daily streamflow at streamgages and gridded monthly climatic data (observed and modeled) were examined across four trend periods (100 years, water years 1921–2020; 75 years, 1946–2020; 50 years, 1971–2020; 30 years 1991–2020) for trends, change points, and other statistical properties indicative of changing conditions. Median peak streamflows generally exhibit upward trends across the State for the 100- and 75-year trend periods and in northern and southern Illinois for the 50- and 30-year trend periods. The medians of the trend magnitudes (normalized by median peak streamflow) range from a 23-percent increase during the 30-year trend period to a 41-percent increase during the 100-year trend period. Streamgages with trends in peak streamflow often also have change points, or abrupt changes, in streamflow magnitude. More than two-thirds of streamgages at the 100- and 75-year trend periods exhibit a trend and change point in median peak streamflow in the same direction. Temporally, clusters of change points are observed in the late 1960s through early 1980s for the 100- and 75-year trend periods and around 2006 for the 50- and 30-year trend periods. Trends in the 90-percent quantile of peak streamflow, which correspond to the 10-percent exceedance probability often used for the design of drainage structures, increased about the same amount as the 50-percent quantile peak streamflows, except at the 100-year trend period, where the 50-percent quantile peak flow increased more for almost all streamgages. The frequency of high flows has also increased in Illinois, with increases in peaks-over-threshold observed across much of the State for the 100- and 75-year trend periods and in northern and southern Illinois for the 50- and 30-year trend periods. Upward trends in observed temperature and observed annual precipitation dominate in all trend periods, with clusters of likely upward trends observed in northern and southern Illinois at the 50- and 30-year trend periods. As expected in response to increasing temperature, the modeled proportion of precipitation falling as snow has largely decreased in the study basins across the State, and modeled potential evapotranspiration has increased. Upward trends in modeled annual runoff, which in this report incorporates only the effects of climatic variation, are observed in the same geographic areas where there are increases in observed annual precipitation. The widespread upward trends in the magnitude of median peak streamflows and the frequency with which high flows occur across the State at the 100- and 75-year trend periods and in northern and southern Illinois at the 50- and 30-year trend periods appear to be driven largely by increases in precipitation based on spatial patterns of these changes and statistical relations between streamflow and climate metrics. Other effects not considered in this report, like urbanization, may be important drivers for certain streamgages in the State. The prevalence of nonstationarity in peak streamflow in Illinois has important implications for peak-flow frequency analysis. Average annual precipitation and the occurrence of extreme precipitation events are expected to increase across the State. If precipitation continues to increase as expected, peak-flow frequency estimates based on older records may no longer represent the hydrologic regime of today, and methods for nonstationary peak-flow frequency analysis may be needed.

Illinois

Empirical estimation of recreational exploitation of burbot, Lota lota, in the Wind River drainage of Wyoming using a multistate capture–recapture model

Burbot, Lota lota (Linnaeus), is a regionally popular sportfish in the Wind River drainage of Wyoming, USA, at the southern boundary of the range of the species. Recent declines in burbot abundances were hypothesised to be caused by overexploitation, entrainment in irrigation canals and habitat loss. This study addressed the overexploitation hypothesis using tagging data to generate reliable exploitation, abundance and density estimates from a multistate capture–recapture model that accounted for incomplete angler reporting and tag loss. Exploitation rate μ was variable among the study lakes and inversely correlated with density. Exploitation thresholds μ 40 associated with population densities remaining above 40% of carrying capacity were generated to characterise risk of overharvest using exploitation and density estimates from tagging data and a logistic surplus-production model parameterised with data from other burbot populations. Bull Lake (μ = 0.06, 95% CI: 0.03–0.11; μ 40 = 0.18) and Torrey Lake (μ = 0.02, 95% CI: 0.00–0.11; μ 40 = 0.18) had a low risk of overfishing, Upper Dinwoody Lake had intermediate risk (μ = 0.08, 95% CI: 0.02–0.32; μ 40 = 0.18) and Lower Dinwoody Lake had high risk (μ = 0.32, 95% CI: 0.10–0.67; μ 40 = 0.08). These exploitation and density estimates can be used to guide sustainable management of the Wind River drainage recreational burbot fishery and inform management of other burbot fisheries elsewhere.

Wyoming

Habitat Suitability Index Models: Field sparrow

Habitat preferences of the field sparrow (Spizella pusilla) are described in this report, which is one of a series of Habitat Suitability Index (HSI) models. A review and synthesis of the literature is followed by the development of a habitat model for the field sparrow throughout its breeding range in the United States. HSI models are designed to be used in conjunction with Habitat Evaluation Procedures previously developed by the U.S. Fish and Wildlife Service.

FWS/OBS

Improved earthquake monitoring in the central and eastern United States in support of seismic assessments for critical facilities

Evaluation of seismic monitoring capabilities in the central and eastern United States for critical facilities - including nuclear powerplants - focused on specific improvements to understand better the seismic hazards in the region. The report is not an assessment of seismic safety at nuclear plants. To accomplish the evaluation and to provide suggestions for improvements using funding from the American Recovery and Reinvestment Act of 2009, the U.S. Geological Survey examined addition of new strong-motion seismic stations in areas of seismic activity and addition of new seismic stations near nuclear power-plant locations, along with integration of data from the Transportable Array of some 400 mobile seismic stations. Some 38 and 68 stations, respectively, were suggested for addition in active seismic zones and near-power-plant locations. Expansion of databases for strong-motion and other earthquake source-characterization data also was evaluated. Recognizing pragmatic limitations of station deployment, augmentation of existing deployments provides improvements in source characterization by quantification of near-source attenuation in regions where larger earthquakes are expected. That augmentation also supports systematic data collection from existing networks. The report further utilizes the application of modeling procedures and processing algorithms, with the additional stations and the improved seismic databases, to leverage the capabilities of existing and expanded seismic arrays.

Open-File Report

Simulation of dissolved organic carbon flux in the Penobscot Watershed, Maine

Dissolved organic carbon (DOC) is an important component of the carbon cycle as a measure of the hydrological transport of carbon between terrestrial carbon pools into soil pools and eventually into streams. As a result, changes in DOC in rivers and streams may indicate alterations in the storage of terrestrial carbon. Exploring the complex interactions between biogeochemical cycling and hydrologic processes, as well as the micro-climate variabilities that impact the rate of DOC fluxes, are challenging because the information is not readily available from in-situ measurements or from empirical models alone. This is particularly true of large-scale watersheds. The Penobscot Watershed is the largest watershed of the Gulf of Maine and the second largest in New England. Its typical soils, with high organic matter and a large forested and wetland landscape, result in higher DOC fluxes than what has been observed previously for most rivers in the northern temperate or boreal zones (Hope et al., 1994; Mulholland, 1997; Aitkenhead and McDowell, 2000). In this study, we emphasized the simulation of streamflow and DOC fluxes from the Penobscot Watershed (and several tributaries within the Penobscot Watershed) using the spatially distributed process-based Regional Hydro-Ecological Simulation System (RHESSys) model. Simulated results were evaluated using field measurements (streamflow, DOC fluxes) and remotely sensed products (Net Primary Production (NPP) and Leaf Area Index (LAI) from Moderate Resolution Imaging Spectroradiometer (MODIS). The average DOC flux for the Penobscot Watershed during 2004-2012 using the RHESSys model was 69 kg C/ha/year. The RHESSys simulated DOC flux is shown to correlate well with observed values, as well as with results previously reported from the empirical Load Estimator (LOADEST) model (71 kg C/ha/year) for 2004-2007 ( Huntington and Aiken, 2013). Our simulated results also show a temporal variation in the amount of DOC flux, indicating that the antecedent DOC concentration from one year can impact the DOC export in following years. Thus, DOC concentration is positively correlated with streamflow and antecedent precipitation, in agreement with previous studies (Ågren et al., 2010; Huntington and Aiken, 2013; Tian et al., 2013 ). The successful application of the rigorous RHESSys model in the Penobscot Watershed makes it a reasonable platform to test future scenarios impacting the hydrology and biogeochemistry within similar large complex watersheds.

Maine

Model methodology for estimating pesticide concentration extremes based on sparse monitoring data

This report describes a new methodology for using sparse (weekly or less frequent observations) and potentially highly censored pesticide monitoring data to simulate daily pesticide concentrations and associated quantities used for acute and chronic exposure assessments, such as the annual maximum daily concentration. The new methodology is based on a statistical model that expresses log-transformed daily pesticide concentration in terms of a seasonal wave, flow-related variability, long-term trend, and serially correlated errors. Methods are described for estimating the model parameters, generating conditional simulations of daily pesticide concentration given sparse (weekly or less frequent) and potentially highly censored observations, and estimating concentration extremes based on the conditional simulations. The model can be applied to datasets with as few as 3 years of record, as few as 30 total observations, and as few as 10 uncensored observations. The model was applied to atrazine, carbaryl, chlorpyrifos, and fipronil data for U.S. Geological Survey pesticide sampling sites with sufficient data for applying the model. A total of 112 sites were analyzed for atrazine, 38 for carbaryl, 34 for chlorpyrifos, and 33 for fipronil. The results are summarized in this report; and, R functions, described in this report and provided in an accompanying model archive, can be used to fit the model parameters and generate conditional simulations of daily concentrations for use in investigations involving pesticide exposure risk and uncertainty.

Scientific Investigations Report