Search USGSSearch

SEARCH · Search USGS

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 109 records · Page 6Linked to original sources

Geologic cross sections and preliminary geologic map of the Questa Area, Taos County, New Mexico

In 2011, the senior authors were contacted by Ron Gardiner of Questa, and Village of Questa Mayor Esther Garcia, to discuss the existing and future groundwater supply for the Village of Questa. This meeting led to the development of a plan in 2013 to perform an integrated geologic, geophysical, and hydrogeologic investigation of the Questa area by the New Mexico Bureau of Geology & Mineral Resources (NMBG), the U.S. Geological Survey (USGS), and New Mexico Tech (NMT). The NMBG was responsible for the geologic map and geologic cross sections. The USGS was responsible for a detailed geophysical model to be incorporated into the NMBG products. NMT was responsible for providing a graduate student to develop a geochemical and groundwater flow model. This report represents the final products of the geologic and geophysical investigations conducted by the NMBG and USGS. The USGS final products have been incorporated directly into the geologic cross sections. The objective of the study was to characterize and interpret the shallow (to a depth of approximately 5,000 ft) three-dimensional geology and preliminary hydrogeology of the Questa area. The focus of this report is to compile existing geologic and geophysical data, integrate new geophysical data, and interpret these data to construct three, detailed geologic cross sections across the Questa area. These cross sections can be used by the Village of Questa to make decisions about municipal water-well development, and can be used in the future to help in the development of a conceptual model of groundwater flow for the Questa area. Attached to this report are a location map, a preliminary geologic map and unit descriptions, tables of water wells and springs used in the study, and three detailed hydrogeologic cross sections shown at two different vertical scales. The locations of the cross sections are shown on the index map of the cross section sheet.

Open-File Report

A systematic literature review of forecasting and predictive models of harmful algal blooms in flowing waters

Occurrences of harmful algal blooms (HABs) in rivers challenge the belief that rivers are not susceptible to HABs because of their short residence times and fluctuating hydrology. Here we present a systematic literature review of predictive and forecasting models for HABs in flowing waters, including rivers, flowing in-stream reservoirs (e.g., run-of-river reservoirs and lock-and-dam systems) and tidal or estuarine systems with riverine processes. The review aimed to understand current and historical modeling approaches for predicting and forecasting river HABs, without restricting to specific taxa, such as cyanobacteria, or modeling endpoints. The review included 162 articles published over nearly 50 years, covering more than 80 rivers worldwide. Eutrophic, non-wadable rivers with in-stream obstruction were commonly modeled, though diverse environmental characteristics were reported. Most articles used algal biomass or chlorophyll as modeling endpoints, with a quarter using novel or unique endpoints. Algal toxins motivated model development in 23% of the articles, however just 5% used algal toxins as an endpoint. Only 6% of the articles modeled benthic HABs; the rest focused on pelagic HABs. There was no standard model used for modeling river HABs. Process-based models were more common (59%) than data-driven approaches (37%), with model formulations ranging from simple to complex, which contrasts with a lake-focused literature review of HAB models that found data-driven models were more common. Models in river settings shared similar input variables as those previously identified for lakes, such as water temperature, nutrients, and light availability. However, streamflow and other transport metrics took prominence in river models compared to lake models. Algal cell physiology (such as growth, predation, and motility) was routinely included as input data or as mathematical formulations in process-based models and these processes were frequently identified as an important predictor by the articles’ authors. Conversely, data-driven models rarely included these processes, instead using predictors related to environmental conditions, such as nutrients, water quality, water temperature, and streamflow. These important proxy predictors have apparent success with modeling overall algal biomass (irrespective of taxa) whereas other factors, such as those related to algal physiology and other biological processes, are likely responsible for more subtle shifts in community composition. These differences highlight the influence of data availability, especially for processes that are difficult, time-consuming, or expensive to measure, on model development and model outcomes, raising questions about the selection of modeling inputs and endpoints. Challenges to advancing river HAB modeling include the lack of site-specific model inputs representing key processes (e.g., photosynthetic parameters and predation rates), overlooked riverine environments like the benthos and side/back-channel areas, lack of information on environmental settings, and poorly reported model performance metrics. This review emphasizes opportunities for advancing river HAB modeling by learning from well-honed estuarine models, supporting current forecasting and operationalization efforts, and developing common datasets for river HAB model development and evaluation.

BioRxiv

MODFLOW-2000 : the U.S. Geological Survey modular ground-water model--documentation of the Advective-Transport Observation (ADV2) Package

Observations of the advective component of contaminant transport in steady-state flow fields can provide important information for the calibration of ground-water flow models. This report documents the Advective-Transport Observation (ADV2) Package, version 2, which allows advective-transport observations to be used in the three-dimensional ground-water flow parameter-estimation model MODFLOW-2000. The ADV2 Package is compatible with some of the features in the Layer-Property Flow and Hydrogeologic-Unit Flow Packages, but is not compatible with the Block-Centered Flow or Generalized Finite-Difference Packages. The particle-tracking routine used in the ADV2 Package duplicates the semi-analytical method of MODPATH, as shown in a sample problem. Particles can be tracked in a forward or backward direction, and effects such as retardation can be simulated through manipulation of the effective-porosity value used to calculate velocity. Particles can be discharged at cells that are considered to be weak sinks, in which the sink applied does not capture all the water flowing into the cell, using one of two criteria: (1) if there is any outflow to a boundary condition such as a well or surface-water feature, or (2) if the outflow exceeds a user specified fraction of the cell budget. Although effective porosity could be included as a parameter in the regression, this capability is not included in this package. The weighted sum-of-squares objective function, which is minimized in the Parameter-Estimation Process, was augmented to include the square of the weighted x-, y-, and z-components of the differences between the simulated and observed advective-front locations at defined times, thereby including the direction of travel as well as the overall travel distance in the calibration process. The sensitivities of the particle movement to the parameters needed to minimize the objective function are calculated for any particle location using the exact sensitivity-equation approach; the equations are derived by taking the partial derivatives of the semi-analytical particle-tracking equation with respect to the parameters. The ADV2 Package is verified by showing that parameter estimation using advective-transport observations produces the true parameter values in a small but complicated test case when exact observations are used. To demonstrate how the ADV2 Package can be used in practice, a field application is presented. In this application, the ADV2 Package is used first in the Sensitivity-Analysis mode of MODFLOW-2000 to calculate measures of the importance of advective-transport observations relative to head-dependent flow observations when either or both are used in conjunction with hydraulic-head observations in a simulation of the sewage-discharge plume at Cape Cod, Massachusetts. The ADV2 Package is then used in the Parameter-Estimation mode of MODFLOW-2000 to determine best-fit parameter values. It is concluded that, for this problem, advective-transport observations improved the calibration of the model and the estimation of ground-water flow parameters, and the use of formal parameter-estimation methods and related techniques produced significant insight into the physical system.

Open-File Report

Methods and guidelines for effective model calibration; with application to UCODE, a computer code for universal inverse modeling, and MODFLOWP, a computer code for inverse modeling with MODFLOW

This report documents methods and guidelines for model calibration using inverse modeling. The inverse modeling and statistical methods discussed are broadly applicable, but are presented as implemented in the computer programs UCODE, a universal inverse code that can be used with any application model, and MODFLOWP, an inverse code limited to one application model. UCODE and MODFLOWP perform inverse modeling, posed as a parameter-estimation problem, by calculating parameter values that minimize a weighted least-squares objective function using nonlinear regression. Minimization is accomplished using a modified Gauss-Newton method, and prior, or direct, information on estimated parameters can be included in the regression. Inverse modeling in many fields is plagued by problems of instability and nonuniqueness, and obtaining useful results depends on (1) defining a tractable inverse problem using simplifications appropriate to the system under investigation and (2) wise use of statistics generated using calculated sensitivities and the match between observed and simulated values, and associated graphical analyses. Fourteen guidelines presented in this work suggest ways of constructing and calibrating models of complex systems such that the resulting model is as accurate and useful as possible.

Water-Resources Investigations Report

Rio Grande transboundary integrated hydrologic model and water-availability analysis, New Mexico and Texas, United States, and northern Chihuahua, Mexico

Changes in population, agricultural development and practices (including shifts to more water-intensive crops), and climate variability are increasing demands on available water resources, particularly groundwater, in one of the most productive agricultural regions in the Southwest—the Rincon and Mesilla Valley parts of Rio Grande Valley, Doña Ana and Sierra Counties, New Mexico, and El Paso County, Texas. The goal of this study was to produce an integrated hydrological simulation model to help evaluate water-management strategies, including conjunctive use of surface water and groundwater for historical conditions, and to support long-term planning for the Rio Grande Project. This report describes model construction and applications by the U.S. Geological Survey, working in cooperation and collaboration with the Bureau of Reclamation. This model, the Rio Grande Transboundary Integrated Hydrologic Model, simulates the most important natural and human components of the hydrologic system, including selected components related to variations in climate, thereby providing a reliable assessment of surface-water and groundwater conditions and processes that can inform water users and help improve planning for future conditions and sustained operations of the Rio Grande Project (RGP) by the Bureau of Reclamation. Model development included a revision of the conceptual model of the flow system, construction of a Transboundary Rio Grande Watershed Model (TRGWM) water-balance model using the Basin Characterization Model, and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model version 2 (referred to as MF-OWHM2). The hydrologic models were developed for and calibrated to historical conditions of water and land use, and parameters were adjusted so that simulated values closely matched available measurements (calibration). The calibrated model was then used to assess the use and movement of water in the Rincon Valley, Mesilla Basin, and northern part of the Conejos-Médanos Basin, with the entire region referred to as the “Transboundary Rio Grande” or TRG. These tools provide a means to understand hydrologic system response to the evolution of water use in the region, its availability, and potential operational constraints of the RGP. The conceptual model identified surface-water and groundwater inflows and outflows that included the movement and use of water both in natural and in anthropogenic systems. The groundwater-flow system is characterized by a layered geologic sedimentary sequence combined with the effects of groundwater pumping, operation of the RGP, natural runoff and recharge, and the application of irrigation water at the land surface that is captured and reused in an extensive network of canals and drains as part of the conjunctive use of water in the region. Historical groundwater-level fluctuations followed a cyclic pattern that were aligned with climate cycles, which collectively resulted in alternating periods of wet or dry years. Periods of drought that persisted for one or more years are associated with low surface-water availability that resulted in higher rates of groundwater-level decline. Rates of groundwater-level decline also increased during periods of agricultural intensification, which necessitated increasing use of groundwater as a source of irrigation water. Agriculture in the area was initially dominated by alfalfa and cotton, but since 1970 more water-intensive pecan orchards and vegetable production have become more common. Groundwater levels substantially declined in subregions where drier climate combined with increased demand, resulting in periods of reduced streamflows. Most of the groundwater was recharged in the Rio Grande Valley floor, and most of the pumpage and aquifer storage depletion was in Mesilla Basin agricultural subregions. A cyclic imbalance between inflows and outflows resulted in the modeled cyclic depletion (groundwater withdrawals in excess of natural recharge) of the groundwater basin during the 75-year simulation period of 1940–2014. Changes in groundwater storage can vary considerably from year to year, depending on land use, pumpage, and climate conditions. Climatic drivers of wet and dry years can greatly affect all inflows, outflows, and water use. Although streamflow and, to a minor extent, precipitation during inter-decadal wet-year periods replenished the groundwater historically, contemporary water use and storage depletion could have reduced the effects of these major recharge events. The average net groundwater flow-rate deficit for 1953–2014 was estimated to be about 1,090 acre-feet per year.

Chihuahua, New Mexico, Texas

Water budget and hydraulic aspects of artificial recharge, south coast of Puerto Rico

An analog model was used to evaluate ground-water conditions on the south coast of Puerto Rico. Water levels during a normal period and during an extended drought were simulated. Recharge and discharge values are reported. The model was also used to evaluate the possibilities of using treated waste water to recharge the aquifer. Three methods were considered: infiltration basins, injection, and irrigation. The tests were planned to determine what changes in water levels would result if certain rates of application were used. Because of the limited vertical hydraulic conductivity, irrigation is suggested as the most practical method of waste-water use. (Woodard-USGS)

Puerto Rico

Oyster model inventory: Identifying critical data and modeling approaches to support restoration of oyster reefs in coastal U.S. Gulf of Mexico waters

Executive Summary Along the coast of the U.S. Gulf of Mexico, the eastern oyster ( Crassostrea virginica ) plays important ecological and economic roles. Commercial landings from this region account for more than 50 percent of all U.S. landings; these oyster reefs also provide varied ecosystem services, including nursery habitat for many fish and macroinvertebrate species, shoreline protection, and water-quality maintenance. Declining trends in both total oyster production and functional reef area across this region have spurred investment in restoration of oyster resources, with specific calls for restoration projects to develop a network of reefs and identify broodstock and sanctuary reef restoration sites. Decision making related to restoration and establishment of a network of oyster reefs in the Gulf of Mexico requires information on both the environment and the effects of the environment on the oyster life cycle (including larval movement, survival, oyster recruitment, reproduction, growth, and mortality). Here, we examined the current state of data and model development in this region with the goal of providing an overview of oyster modeling approaches and an inventory of available data and existing oyster models. This report is meant to provide an overview to managers for understanding existing efforts and identify a path forward to most efficiently inform oyster resource management and restoration planning in moving from a single reef management approach to a reef network management approach. Numerous models related to some aspect of the oyster life cycle have been built, calibrated, and validated for various Gulf of Mexico estuaries over the last few decades (over 30 models identified). These models, which could inform site restoration, can be classified into four approaches: (1) oyster Habitat Suitability Index (HSI) models; (2) larval transport models; (3) on-reef oyster models that may include oyster growth, mortality and reproduction, and substrate persistence; and (4) coupled larval transport on-reef metapopulation models that simulate the entire oyster life cycle. The data requirements, model complexity and assumptions, and transferability vary by approach. Specifically, some approaches may offer greater accessibility, flexibility, and transferability spatially or temporally, with minimal data input, but only provide broad information to support site selection. In contrast, other approaches may require significant site-specific data for their construction and validation but may provide more accurate and location-specific data to support site selection for broodstock reefs. Regardless of modeling approach used, data on environmental drivers, such as salinity, water temperature, or water flow impacting oyster metabolism and movement, are required at appropriate spatial and temporal scales. While numerous data collection platforms, environmental models, and research products exist within Gulf of Mexico estuaries to provide important environmental data to use as drivers in the oyster models, significant variability in temporal and spatial coverage of the data, and variation in the availability of future condition models, exists across estuaries. This variation influences the spatial and temporal scales at which oyster models may be developed and impacts the calibration and validation of the oyster models within a given estuary, affecting its potential ability to address specific management or restoration questions. While multiple modeling approaches exist for informing site selection of broodstock or sanctuary oyster reefs, the development, calibration, and validation of a single modeling platform presents the most efficient, transferable, and useful tool for managers across the Gulf of Mexico. The development of a single modeling platform would involve using standardized input variables, governing equations, and assumptions for the modeled oyster processes and outputs, and for standardized calibration and validation procedures that could be applied within each estuary. The differences among estuary applications would require substituting only estuary-specific environmental data, and calibrating and validating the modeling approach with local oyster data. Two modeling approaches likely to be useful include (1) development of a general geospatial HSI modeling framework that could be applied consistently across estuaries and (2) a mechanistic coupled larval transport on-reef metapopulation model requiring only estuarine specific calibration and hydrodynamic models. Both approaches benefit from existing work across multiple Gulf of Mexico estuaries and could provide valuable support for oyster restoration, but may differ in their ability to address specific questions related to oyster restoration. HSI models specifically guide restoration practitioners in determining suitable habitat based on available data. The HSI approach, while currently more widely used and accessible, requires more development of larval suitability and larval input and output components in order to inform reef connectivity. A metapopulation approach considering the full oyster life cycle that simulates both on-reef oyster growth, mortality, reproduction, substrate persistence, and larval transport (ideally with larval growth and mortality) would provide the greatest detail and level of understanding but requires significant up-front investment. The larval oyster model and on-reef oyster model are usually developed independently for systems, although the two approaches can be coupled to represent the entire oyster life cycle in order to characterize and assess a reef metapopulation. This approach may be less accessible and much more data-intensive, however, and it requires some expertise to run and apply to inform oyster resource management. Ultimately, the development of single modeling platforms for each of these approaches would provide flexible tools applicable across all Gulf of Mexico oyster supporting estuaries. By using a single platform for model development, testing, calibrating and validating, and evaluation of modeled future scenarios, oyster restoration scientists and managers would not only be able to examine different scenario outcomes within a single estuary, but could also have comparable modeled results to evaluate potential outcomes, across estuaries and regions, that are not confounded by varying modeled data inputs, governing equations, assumptions, or user judgement.

Alabama, Florida, Louisiana, Mississippi, Texas

Report on the workshop 'Global modelling of biodiversity and ecosystem services'

A three-day workshop on ‘Global Modelling of Biodiversity and Ecosystem Services’, was held in the Hague, Netherlands, from 24th to 26th June 2019. The workshop, attended by 35 modelling and scenario-building experts, was organised on behalf of the former IPBES1 expert group on scenarios and models of the first IPBES work programme by its interim technical support unit, and hosted by the PBL Netherlands Environmental Assessment Agency. The workshop drew on the ‘nature futures’ participatory scenario-building exercise initiated by the IPBES expert group on scenarios and models, and other biodiversity modelling initiatives such as the ISIMIP project2 working on adding biodiversity to the Shared Socioeconomic Pathways (SSPs) scenarios framework, the 'bending the curve' initiative3 led by IIASA4 and WWF5, and GEOBON6 working on modelling Essential Biodiversity Variables. The workshop was a step towards coordinating across biodiversity modelling initiatives, to build on each other’s work, and to seek synergies for the production of innovative scenarios on biodiversity and ecosystem services to inform the post-2020 agenda of the Convention on Biological Diversity, as well as the Sustainable Development Goals. The aims of the workshop were to: 1. Compile material as input for a first draft of the fifth Global Biodiversity Outlook (GBO-5) based on recent scenario work, including the ‘bending the curve’ scenarios and the newly developed PBL scenarios (modified from the Rio+20 scenarios), and existing models (to be completed by August 2019) 2. Develop a protocol for modelling trends and near term projections on indicators relevant to the Nature Futures Framework7 using models that are readily available (to be completed by early 2020) 3. Set the agenda and define the aims for a larger meeting at the end of 2019 to discuss the long term strategy towards the development of appropriate indicators and models to produce Nature Futures scenarios (to continue beyond Jan 2020)

Report

A review of surface energy balance models for estimating actual evapotranspiration with remote sensing at high spatiotemporal resolution over large extents

Many approaches have been developed for measuring or estimating actual evapotranspiration ( ET a ), and research over many years has led to the development of remote sensing methods that are reliably reproducible and effective in estimating ET a . Several remote sensing methods can be used to estimate ET a at the high spatial resolution of agricultural fields and the large extent of river basins. More complex remote sensing methods apply an analytical approach to ET a estimation using physically based models of varied complexity that require a combination of ground-based and remote sensing data, and are grounded in the theory behind the surface energy balance model. This report, funded through cooperation with the International Joint Commission, provides an overview of selected remote sensing methods used for estimating water consumed through ET a and focuses on Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC) and Operational Simplified Surface Energy Balance (SSEBop), two energy balance models for estimating ET a that are currently applied successfully in the United States. The METRIC model can produce maps of ET a at high spatial resolution (30 meters using Landsat data) for specific areas smaller than several hundred square kilometers in extent, an improvement in practice over methods used more generally at larger scales. Many studies validating METRIC estimates of ET a against measurements from lysimeters have shown model accuracies on daily to seasonal time scales ranging from 85 to 95 percent. The METRIC model is accurate, but the greater complexity of METRIC results in greater data requirements, and the internalized calibration of METRIC leads to greater skill required for implementation. In contrast, SSEBop is a simpler model, having reduced data requirements and greater ease of implementation without a substantial loss of accuracy in estimating ET a . The SSEBop model has been used to produce maps of ET a over very large extents (the conterminous United States) using lower spatial resolution (1 kilometer) Moderate Resolution Imaging Spectroradiometer (MODIS) data. Model accuracies ranging from 80 to 95 percent on daily to annual time scales have been shown in numerous studies that validated ET a estimates from SSEBop against eddy covariance measurements. The METRIC and SSEBop models can incorporate low and high spatial resolution data from MODIS and Landsat, but the high spatiotemporal resolution of ET a estimates using Landsat data over large extents takes immense computing power. Cloud computing is providing an opportunity for processing an increasing amount of geospatial “big data” in a decreasing period of time. For example, Google Earth Engine TM has been used to implement METRIC with automated calibration for regional-scale estimates of ET a using Landsat data. The U.S. Geological Survey also is using Google Earth Engine TM to implement SSEBop for estimating ET a in the United States at a continental scale using Landsat data.

Scientific Investigations Report

An interface for simulating radiative transfer in and around volcanic plumes with the Monte Carlo radiative transfer model McArtim

This report describes two software tools that, when used as front ends for the three-dimensional backward Monte Carlo atmospheric-radiative-transfer model (RTM) McArtim, facilitate the generation of lookup tables of volcanic-plume optical-transmittance characteristics in the ultraviolet/visible-spectral region. In particular, the differential optical depth and derivatives thereof (that is, weighting functions), with regard to a change in SO 2 column density or aerosol optical thickness, can be simulated for a specific measurement geometry and a representative range of plume conditions. These tables are required for the retrieval of SO 2 column density in volcanic plumes, using the simulated radiative-transfer/differential optical-absorption spectroscopic (SRT-DOAS) approach outlined by Kern and others (2012). This report, together with the software tools published online, is intended to make this sophisticated SRT-DOAS technique available to volcanologists and gas geochemists in an operational environment, without the need for an indepth treatment of the underlying principles or the low-level interface of the RTM McArtim.

Open-File Report

Summary of hydrologic modeling for the Delaware River Basin using the Water Availability Tool for Environmental Resources (WATER)

The Water Availability Tool for Environmental Resources (WATER) is a decision support system for the nontidal part of the Delaware River Basin that provides a consistent and objective method of simulating streamflow under historical, forecasted, and managed conditions. In order to quantify the uncertainty associated with these simulations, however, streamflow and the associated hydroclimatic variables of potential evapotranspiration, actual evapotranspiration, and snow accumulation and snowmelt must be simulated and compared to long-term, daily observations from sites. This report details model development and optimization, statistical evaluation of simulations for 57 basins ranging from 2 to 930 km 2 and 11.0 to 99.5 percent forested cover, and how this statistical evaluation of daily streamflow relates to simulating environmental changes and management decisions that are best examined at monthly time steps normalized over multiple decades. The decision support system provides a database of historical spatial and climatic data for simulating streamflow for 2001–11, in addition to land-cover and general circulation model forecasts that focus on 2030 and 2060. WATER integrates geospatial sampling of landscape characteristics, including topographic and soil properties, with a regionally calibrated hillslope-hydrology model, an impervious-surface model, and hydroclimatic models that were parameterized by using three hydrologic response units: forested, agricultural, and developed land cover. This integration enables the regional hydrologic modeling approach used in WATER without requiring site-specific optimization or those stationary conditions inferred when using a statistical model.

Delaware, Maryland, New Jersey, New York, Pennsylv

Methods for estimating drought streamflow probabilities for Virginia streams

Maximum likelihood logistic regression model equations used to estimate drought flow probabilities for Virginia streams are presented for 259 hydrologic basins in Virginia. Winter streamflows were used to estimate the likelihood of streamflows during the subsequent drought-prone summer months. The maximum likelihood logistic regression models identify probable streamflows from 5 to 8 months in advance. More than 5 million streamflow daily values collected over the period of record (January 1, 1900 through May 16, 2012) were compiled and analyzed over a minimum 10-year (maximum 112-year) period of record. The analysis yielded the 46,704 equations with statistically significant fit statistics and parameter ranges published in two tables in this report. These model equations produce summer month (July, August, and September) drought flow threshold probabilities as a function of streamflows during the previous winter months (November, December, January, and February). Example calculations are provided, demonstrating how to use the equations to estimate probable streamflows as much as 8 months in advance.

Virginia

One-Dimensional Transport with Inflow and Storage (OTIS): A Solute Transport Model for Streams and Rivers

OTIS is a mathematical simulation model used to characterize the fate and transport of water-borne solutes in streams and rivers. The governing equation underlying the model is the advection-dispersion equation with additional terms to account for transient storage, lateral inflow, first-order decay, and sorption. This equation and the associated equations describing transient storage and sorption are solved using a Crank-Nicolson finite-difference solution. OTIS may be used in conjunction with data from field-scale tracer experiments to quantify the hydrologic parameters affecting solute transport. This application typically involves a trial-and-error approach wherein parameter estimates are adjusted to obtain an acceptable match between simulated and observed tracer concentrations. Additional applications include analyses of nonconservative solutes that are subject to sorption processes or first-order decay. OTIS-P, a modified version of OTIS, couples the solution of the governing equation with a nonlinear regression package. OTIS-P determines an optimal set of parameter estimates that minimize the squared differences between the simulated and observed concentrations, thereby automating the parameter estimation process. This report details the development and application of OTIS and OTIS-P. Sections of the report describe model theory, input/output specifications, sample applications, and installation instructions.

Water-Resources Investigations Report

Assessment of vulnerabilities and opportunities to restore marsh sediment supply at Nisqually River Delta, west-central Washington

A cascading set of hazards to coastal environments is intimately tied to sediment transport and includes the flooding and erosion of shorelines and habitats that support communities, industry, infrastructure, and ecosystem functions (for example, habitats critical to fisheries). This report summarizes modeling and measurement data used to evaluate the sediment budget of the Nisqually River Delta, the vulnerability of the largest estuary restoration project in Puget Sound at the Billy Frank Jr. Nisqually National Wildlife Refuge, and the role of coastal hydrodynamics and potential restoration alternatives for recovering sediment delivery to its marshes. The 2009 Brown’s Farm Restoration achieved many goals toward recovering salmon habitat, but the understanding of the delta and restoration area sediment budgets remain poorly quantified. Specifically, quantitative estimates of the amount of sediment delivered to the delta and restored marsh areas, which had subsided in response to historical diking and draining for grazing, were identified as important information needs. Forecasts of potential outcomes of proposed adaptive distributary channel restoration actions were also prioritized to inform potential solutions. These estimates can be used to evaluate whether sufficient sediment is available for marsh recovery downstream from Alder Lake, which traps about 90 percent the Nisqually River sediment load that could reach the delta . Additionally, quantitative sediment information was identified to help prioritize opportunities to recover and maintain the area marshes and guide ecosystem restoration investments across the delta to reduce the vulnerability of the system to drowning under projected sea level rise. A coupled, numerical hydrodynamic-sediment transport model and measurements of the sediment load just upstream from the delta were used to evaluate the (1) availability of sediment for marsh recovery, (2) sediment transport dynamics across the estuary, and (3) potential outcomes of distributary reconnection alternatives under existing and projected conditions of streamflow and sea level. Modeling and measurements indicated that the volume of fluvial sediment load reaching and accumulating in the restoration area ranges from 7 to 32 percent and identified that restoration alternatives could recover about an additional 10–12 percent under current and projected sea-level rise by the year 2100. At these rates of sediment delivery, 85–200+ years may be necessary to fill for marsh vegetation development and maintenance. The model also reveals the sensitivity of sediment transport and accumulation to sediment properties, hydrodynamics, and wave conditions. The low sediment accumulation results in large part because of the role of waves in directing sediment transport offshore and challenges of restoring geomorphic processes suited to maintaining habitat structure where opportunity exists or least conflicts with land use. The findings therefore have implications for siting, phasing, and implementing strategies to route and retain sediment. This study shows that opportunities to recover sediment higher in the tidal prism, where a greater hydraulic gradient and gravity could promote progradation and greater sediment retention, may be more effective than alternatives lower in the tidal prism implemented to date and assessed in this study. Furthermore, the modeling indicates that distributary channel restoration also may provide additional benefits to society by reducing flood stage, and therefore, flood hazards surrounding the delta.

Washington

Evaluation of actual evapotranspiration rates from the Operational Simplified Surface Energy Balance (SSEBop) model in Florida and parts of Alabama and Georgia, 2000–17

Evapotranspiration (ET) is the water-vapor flux transported from the surface of the Earth into the atmosphere and is the sum of surface water directly evaporated and subsurface water transpired by plants. ET rates are commonly estimated by using potential or reference ET, which might differ from actual ET rates. Actual evapotranspiration (ETa) rates can be estimated by using the Operational Simplified Surface Energy Balance (SSEBop) model. This report evaluates SSEBop ETa rates at the point and basin scales in Florida and parts of Alabama and Georgia for 2000–17. ETa rates computed by using data from 24 micrometeorological stations in Florida are referred to as mETa rates and were used to quantify biases in the SSEBop ETa rates, stratified by generalized land-use type. Bias was computed as mETa minus SSEBop ETa rates for given generalized land-use types, and bias-correction equations were computed by using least-squares regressions. In addition to mETa rates at station locations, annual average ETa rates calculated from the application of a water-balance method to 55 basins in Florida and parts of Alabama and Georgia were used to assess the accuracy of the annual SSEBop ETa rates at the basin scale. Another independent model used to simulate ETa rates was based on monthly reference ET from the statewide daily reference evapotranspiration (ETo) gridded dataset for Florida computed by using Geostationary Operational Environmental Satellite estimates of solar radiation (GOES ETo). ETa at grid points was computed as monthly GOES ETo multiplied by ratios of monthly mETa to GOES ETo, computed at micrometeorological stations and stratified by each generalized land-use type. The coefficient of determination (R 2 ) between monthly mETa and SSEBop ETa rates for all stations combined improved from 0.37 before bias correction of SSEBop ETa rates to 0.79 after the bias correction stratified by land-use type. For individual land-uses types, R 2 varied from 0.59 for the monthly mETa at a station in the land-use type forest to 0.82 for the monthly mETa at stations in the land-use type shallow-water-table pasture. Root-mean-square error (RMSE) was computed as a function of the difference between SSEBop ETa rates and mETa rates. RMSE of monthly SSEBop ETa rates was 1.27 inches per month before the bias corrections improved to 0.73 inch per month after the bias corrections. RMSE for bias-corrected annual SSEBop ETa rates based on micrometeorological stations with complete years of records ranged from 2.01 inches per year (in/yr) for the land-use type of agriculture to 5.73 in/yr for the land-use type of deep water-table pasture, or 4.96 and 21.21 percent errors relative to annual mETa rates, respectively. Bias-corrected annual SSEBop ETa rates were also compared to annual ETa rates computed by using a water-balance method (wbETa) for 55 basins in Florida. Differences in bias-corrected average annual SSEBop ETa rates and average annual wbETa rates for the 55 basins ranged from −3.67 to 5.29 in/yr (−9.24 to 17.36 percent). RMSE when computed as a function of the differences between annual SSEBop ETa rates and wbETa rates decreased, on average, from 4.13 in/yr for the uncorrected bias SSEBop ETa rates to 1.95 in/yr for the bias-corrected SSEBop rates. The average annual bias-corrected SSEBop ETa rates, from all basins, was 36.46 in/yr or 3.41 percent lower than the average annual wbETa rate of 37.79 inches. Bias in SSEBop ETa rates varies based on time step (monthly versus annual), scale (point, basin, statewide), and land-use type. Applications to hydrologic models should consider bias relative to the inherent error in models. Bias-corrected SSEBop ETa rates could be used as calibration targets in models of hydrologic processes, such as groundwater models. Annual bias in SSEBop ETa introduced to the model calibration is typically below the margin of error associated with typical residuals in model simulations, depending on scale. Surface-water and groundwater-flow models with RMSEs on the order of a few feet could benefit from bias-corrected SSEBop values of ETa.

Alabama, Florida, Georgia

Vital signs monitoring plan for the Klamath Network: Phase I report

This report chronicles the Phase 1 stage of the vital signs monitoring program for the Klamath Network. It consists of two chapters and eleven appendixes. The purposes of Chapter One are to 1) describe the network administrative structure and approach to planning; 2) introduce the Klamath Network parks, their resources, and environmental settings; 3) explain the need for monitoring changes in resources and supporting environments; 4) identify key information gaps that limit understanding of how to best achieve these monitoring goals. The purpose of Chapter Two is to develop the descriptive information provided in Chapter One into a conceptual basis for vital signs monitoring and to present the Network’s initial suite of conceptual models. The Report Appendices provide in-depth information on a variety of topics researched in preparation of the report, including: detailed natural resource profiles for each park, supporting policies and guidelines, regional fire regimes, vegetation types of the parks, exotic species threats, interagency monitoring programs, air issues, water quality (Phase 1 Report), Network vital signs (Scoping Summary Report), rare species, and rare habitats of the parks.

Report

Simulation of regional ground-water flow in the Upper Deschutes Basin, Oregon

This report describes a numerical model that simulates regional ground-water flow in the upper Deschutes Basin of central Oregon. Ground water and surface water are intimately connected in the upper Deschutes Basin and most of the flow of the Deschutes River is supplied by ground water. Because of this connection, ground-water pumping and reduction of artificial recharge by lining leaking irrigation canals can reduce the amount of ground water discharging to streams and, consequently, streamflow. The model described in this report is intended to help water-management agencies and the public evaluate how the regional ground-water system and streamflow will respond to ground-water pumping, canal lining, drought, and other stresses. Ground-water flow is simulated in the model by the finite-difference method using MODFLOW and MODFLOWP. The finite-difference grid consists of 8 layers, 127 rows, and 87 columns. All major streams and most principal tributaries in the upper Deschutes Basin are included. Ground-water recharge from precipitation was estimated using a daily water-balance approach. Artificial recharge from leaking irrigation canals and on-farm losses was estimated from diversion and delivery records, seepage studies, and crop data. Ground-water pumpage for irrigation and public water supplies, and evapotranspiration are also included in the model. The model was calibrated to mean annual (1993-95) steady-state conditions using parameter-estimation techniques employing nonlinear regression. Fourteen hydraulic-conductivity parameters and two vertical conductance parameters were determined using nonlinear regression. Final parameter values are all within expected ranges. The general shape and slope of the simulated water-table surface and overall hydraulic-head distribution match the geometry determined from field measurements. The fitted standard deviation for hydraulic head is about 76 feet. The general magnitude and distribution of ground-water discharge to streams is also well simulated throughout the model. Ground-water discharge to streams in the area of the confluence of the Deschutes, Crooked, and Metolius Rivers is closely matched. The model was also calibrated to transient conditions from 1978 to 1997 using traditional trial-and-error methods. Climatic cycles during this period provided an excellent regional hydrologic signal for calibration. Climate-driven water-level fluctuations are simulated with reasonable accuracy over most of the model area. The timing and magnitude of simulated water-level fluctuations caused by annual pulses of recharge from precipitation match those observed reasonably well, given the limitations of the time discretization in the model. Water-level fluctuations caused by annual canal leakage are simulated very well over most of the area where such fluctuations occur. The transient model also simulates the volumetric distribution and temporal variations in ground-water discharge reasonably well. The match between simulated and measured volume of and variations in ground-water discharge is, however, somewhat dependent on geographic scale. The rates of and variations in ground-water discharge are matched best at regional scales. Example simulations were made to demonstrate the utility of the model for evaluating the effects of ground-water pumping or canal lining. Pumping simulations show that pumped water comes largely from aquifer storage when pumping begins, but as the water table stabilizes, the pumping increasingly diminishes the discharge to streams and, hence, streamflow. The time it takes for pumping to affect streamflow varies spatially depending, in general, on the location of pumping relative to the discharge areas. Canal-lining simulations show similar effects.

Water-Resources Investigations Report

Forecasting storm-induced coastal flooding for 21st century sea-level rise scenarios in the Hawaiian, Mariana, and American Samoan Islands

Oceanographic, coastal engineering, ecologic, and geospatial data and tools were combined to evaluate the increased risks of storm-induced coastal flooding in the populated Hawaiian, Mariana, and American Samoan Islands as a result of climate change and sea-level rise. We followed a hybrid (dynamical and statistical) downscaling approach to map flooding due to waves and storm surge at 10-square meter resolution along all 1,870 kilometers of these islands’ coastlines for annual (1-year), 20-year, and 100-year return-interval storm events and +0.00 meter (m), +0.25 m, +0.50 m, +1.00 m, +1.50 m, +2.00 m, and +3.00 m sea-level rise scenarios. We quantified the coastal flood depths and extents using the latest climate forcing from Intergovernmental Panel for Climate Change’s Sixth Assessment Report Coupled Model Intercomparison Project. The data generated using these methods provide stakeholders and decision makers with a spatially explicit, rigorous valuation of how, where, and when climate change and sea-level rise increase coastal storm-induced flooding to help identify areas where management and (or) restoration could potentially help reduce the risk to, and increase the resiliency of, the coastal communities in the populated Hawaiian, Mariana, and American Samoan Islands.

Hawaii