USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modeling Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21Linked to original sources

Helicopter electromagnetic survey of the Model Land Area, Southeastern Miami-Dade County, Florida

This report describes a helicopter electromagnetic survey flown over the Model Land Area in southeastern Miami-Dade County, Florida, to map saltwater intrusion in the Biscayne aquifer. The survey, which is located south and east of Florida City, Florida, covers an area of 115 square kilometers with a flight-line spacing of 400 meters. A five-frequency, horizontal, coplanar bird with frequencies ranging from 400 to 100,000 Hertz was used. The data were interpreted using differential resistivity analysis and inversion to produce cross sections and resistivity depth-slice maps. The depth of investigation is as deep as 100 meters in freshwater-saturated portions of the Biscayne aquifer and the depth diminishes to about 50 meters in areas that are intruded by saltwater. The results compare favorably with ground-based, time-domain electromagnetic soundings and induction logs from observation wells in the area. The base of a high-resistivity, freshwater-saturated zone mapped in the northern 2 kilometers of the survey area corresponds quite well with the base of the surficial aquifer that has been determined by drilling. In general, saltwater in the survey area extends 9 to 12 kilometers inland from the coast; however, there is a long nose of saltwater centered along the Card Sound Road Canal that extends 15 kilometers inland. The cause of this preferential intrusion is likely due to uncontrolled surface flow along the canal and subsequent leakage of saltwater into the aquifer. Saltwater also extends farther inland in the area between U.S. Highway 1 and Card Sound Road than it does to the west of this area. Until 1944, a railroad grade occupied the current location of U.S. Highway 1. Borrow ditches associated with the railroad grade connected to Barnes Sound and allowed saltwater to flow during droughts and storm surges to within a few kilometers of Florida City. Relicts of this saltwater that settled to the bottom of the Biscayne aquifer can be seen in the helicopter electromagnetic data. The area to the west of U.S. Highway 1 is more resistive in the upper 10 meters than the area to the east of the road; this reflects the influence of surface-water flows that are blocked by U.S. Highway 1. Between Card Sound Road and U.S. Highway 1, resistivities are slightly lower compared to adjacent areas. In the southern portion of the survey area, the surficial aquifer underlying the Biscayne aquifer is more resistive; this indicates that it contains fresher water than that found at the base of the Biscayne aquifer.

Florida↗

An individual-based model for predicting dynamics of a newly established Mexican wolf (Canis lupus baileyi) population—Final report

Project Summary The Mexican wolf recovery team proposed to establish other populations of Mexican wolves ( Canis lupus baileyi ) in the Southwest (U.S. Fish and Wildlife Service, 1982). We were tasked to conduct an extensive simulation modeling exercise to determine release strategies (in conjunction with management actions) that best predict establishment of a new Mexican wolf population. Our objectives were to determine optimal release and management strategies for population establishment and growth. This is a retrospective analysis utilizing data from 1998 to 2014, and during this period, we divided management strategies into two phases; (1) 1998–2008, where nuisance wolves (i.e., wolves that exhibit nuisance behavior or depredate livestock) were managed primarily through lethal removals and removals to captivity, and (2) 2009–2014, when lethal removals ceased and diversionary feeding was provided to denning packs to dissuade wolves from conflict with humans. Management strategies from the second phase are being used for management of the current Mexican wolf population, and demographic rates derived from alternate population modeling in Vortex incorporating post-2008 wolf data are being used to guide future recovery efforts. Therefore, demographic rates estimated from our retrospective analysis will differ (i.e., due to our unique approach to the analyses and the demographic rates being derived from a different dataset), and are intended solely to address the objectives of this report, and are not intended as basis for the development of management recommendations for the current Mexican wolf population. Using individual-based models, we tested dozens of scenarios and derived an optimal release strategy that had the highest probability of establishing a new population and which maximized subsequent post-release growth, and in this report, we present these model results. Findings from this research will improve our understanding of release strategies that yield growing populations, advance our understanding of the demands of reintroducing large carnivores, and provide insight into beneficial strategies that could aid other species reintroduction programs.

Arizona, New Mexico↗

Factors controlling storm impacts on coastal barriers and beaches - A preliminary basis for near real-time forecasting

Analysis of ground conditions and meteorological and oceanographic parameters for some of the most severe Atlantic and Gulf Coast storms in the U.S. reveals the primary factors affecting morphological storm responses of beaches and barrier islands. The principal controlling factors are storm characteristics, geographic position relative to storm path, timing of storm events, duration of wave exposure, wind stress, degree of flow confinement, antecedent topography and geologic framework, sediment textures, vegetative cover, and type and density of coastal development. A classification of commonly observed storm responses demonstrates the sequential interrelations among (1) land elevations, (2) water elevations in the ocean and adjacent lagoon (if present), and (3) stages of rising water during the storm. The predictable coastal responses, in relative order from high frequency beach erosion to low frequency barrier inundation, include: beach erosion, berm migration, dune erosion, washover terrace construction, perched fan deposition, sheetwash, washover channel incision, washout formation, and forced and unforced ebb flow. Near real-time forecasting of expected storm impacts is possible if the following information is available for the coast: a detailed morphological and topographic characterization, accurate storm-surge and wave-runup models, the real-time reporting of storm parameters, accurate forecasts of the storm position relative to a particular coastal segment, and a conceptual model of geological processes that encompasses observed morphological changes caused by extreme storms.

Journal of Coastal Research↗

Dynamic computer model for the metallogenesis and tectonics of the Circum-North Pacific

The digital files on this report consist of a dynamic computer model of the metallogenesis and tectonics of the Circum-North Pacific, and background articles, figures, and maps. The tectonic part of the dynamic computer model is derived from a major analysis of the tectonic evolution of the Circum-North Pacific which is also contained in directory tectevol. The dynamic computer model and associated materials on this CD-ROM are part of a project on the major mineral deposits, metallogenesis, and tectonics of the Russian Far East, Alaska, and the Canadian Cordillera. The project provides critical information on bedrock geology and geophysics, tectonics, major metalliferous mineral resources, metallogenic patterns, and crustal origin and evolution of mineralizing systems for this region. The major scientific goals and benefits of the project are to: (1) provide a comprehensive international data base on the mineral resources of the region that is the first, extensive knowledge available in English; (2) provide major new interpretations of the origin and crustal evolution of mineralizing systems and their host rocks, thereby enabling enhanced, broad-scale tectonic reconstructions and interpretations; and (3) promote trade and scientific and technical exchanges between North America and Eastern Asia.

Open-File Report↗

Three-dimensional hydrogeologic framework model of the Russian River watershed, California

The Russian River watershed is in northern Sonoma County and southern Mendocino County, California, in the northern part of the California Coast Ranges. The Russian River serves as a supply for agricultural irrigation and for municipal, domestic, and commercial uses. Through a cooperative agreement with the California State Water Resources Control Board and Sonoma County Water Agency, the U.S. Geological Survey has completed studies to better understand the hydrogeologic system and develop numerical hydrologic modeling tools to evaluate and aid in managing groundwater resources. This report focuses on the development of a digital three-dimensional hydrogeologic framework model of the Russian River watershed for use in groundwater resource assessment and numerical models. The digital three-dimensional hydrogeologic framework model of the Russian River watershed portrays the altitude, thickness, and extent of five hydrogeologic units. These five hydrogeologic units include (1) a basement unit, (2) the Sonoma Volcanics, (3) a consolidated sedimentary rock unit, (4) an unconsolidated sediment unit, and (5) channel alluvium. Model input data were compiled from published geologic maps, interpreted well data, and a model of the top of basement derived from gravity data. These data were used to construct surfaces that represent the upper and lower subsurface boundaries of each hydrogeologic unit. Top surfaces were created for the five hydrogeologic units and then stacked in three dimensions to create a solid-volume digital model. The digital three-dimensional hydrogeologic framework model described in this report and the corresponding data represent the generalized geometry of the subsurface geologic units; the model reproduces the input geologic data with reasonable accuracy and is consistent with previously published subsurface conceptualizations of the region. The model indicates the overall geometry of the basement within the watershed and the spatial extent, altitude, and thickness of the basin-filling units. The hydrogeologic framework model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

California↗

Late Quaternary sedimentological and climate changes at Lake Bosumtwi Ghana: new constraints from laminae analysis and radiocarbon age modeling

The Lake Bosumtwi sediment record represents one of the longest and highest-resolution terrestrial records of paleoclimate change available from sub-Saharan Africa. Here we report a new sediment age model framework for the last ~ 45 cal kyr of sedimentation using a combination of high-resolution radiocarbon dating, Bayesian age-depth modeling and lamination counting. Our results highlight the practical limits of these methods for reducing age model uncertainties and suggest that even with very high sampling densities, radiocarbon uncertainties of at least a few hundred years are unavoidable. Age model uncertainties are smallest during the Holocene (205 yr) and the glacial (360 yr) but are large at the base of the record (1660 yr), due to a combination of decreasing sample density, larger calibration uncertainties and increases in radiocarbon age scatter. For portions of the chronology older than ~ 35 cal kyr, additional considerations, such as the use of a low-blank graphitization system and more rigorous sample pretreatment were necessary to generate a reliable age depth model because of the incorporation of small amounts of younger carbon. A comparison of radiocarbon age model results and lamination counts over the time interval ~ 15–30 cal kyr agree with an overall discrepancy of ~ 10% and display similar changes in sedimentation rate, supporting the annual nature of sediment laminations in the early part of the record. Changes in sedimentation rates reconstructed from the age-depth model indicate that intervals of enhanced sediment delivery occurred at 16–19, 24 and 29–31 cal kyr, broadly synchronous with reconstructed drought episodes elsewhere in northern West Africa and potentially, with changes in Atlantic meridional heat transport during North Atlantic Heinrich events. These data suggest that millennial-scale drought events in the West African monsoon region were latitudinally extensive, reaching within several hundred kilometers of the Guinea coast. This is inconsistent with a simple southward shift in the mean position of the monsoon rainbelt, and requires changes in moisture convergence as a result of either a reduction in the moisture content of the tropical rainbelt, decreased convection, or both.

Lake Bosumtwi↗

Water-quality data from 2002 to 2003 and analysis of data gaps for development of total maximum daily loads in the Lower Klamath River Basin, California

The U.S. Geological Survey (USGS) collected water-quality data during 2002 and 2003 in the Lower Klamath River Basin, in northern California, to support studies of river conditions as they pertain to the viability of Chinook and Coho salmon and endangered suckers. To address the data needs of the North Coast Regional Water Quality Control Board for the development of Total Maximum Daily Loads (TMDLs), water temperature, dissolved oxygen, specific conductance, and pH were continuously monitored at sites on the Klamath, Trinity, Shasta, and Lost Rivers. Water-quality samples were collected and analyzed for selected nutrients, organic carbon, chlorophyll-a, pheophytin-a, and trace elements. Sediment oxygen demand was measured on the Shasta River. Results of analysis of the data collected were used to identify locations in the Lower Klamath River Basin and periods of time during 2002 and 2003 when river conditions were more likely to be detrimental to salmonid or sucker health because of occasional high water temperatures, low dissolved oxygen, and conditions that supported abundant populations of algae and aquatic plants. The results were also used to assess gaps in data by furthering the development of the conceptual model of water flow and quality in the Lower Klamath River Basin using available data and the current understanding of processes that affect water quality and by assessing needs for the develoment of mathematical models of the system. The most notable gap in information for the study area is in sufficient knowledge about the occurrence and productivity of algal communities. Other gaps in data include vertical water-quality profiles for the reservoirs in the study area, and in an adequate understanding of the chemical oxygen demands and the sediment oxygen demands in the rivers and of the influence of riparian shading on the rivers. Several mathematical models are discussed in this report for use in characterizing the river systems in the study area; also discussed are the specific data needed for the models, and the spatial and temporal data available as boundary conditions. The models will be useful for the future development of TMDLs for temperature, nutrients, and dissolved oxygen and for assessing the role of natural and anthropogenic sources of heat, oxygen-producing and -consuming substances, and nutrients in the Klamath, Shasta, and Lost Rivers.

Scientific Investigations Report↗

Fault segmentation: New concepts from the Wasatch Fault Zone, Utah, USA

The question of whether structural segment boundaries along multisegment normal faults such as the Wasatch fault zone (WFZ) act as persistent barriers to rupture is critical to seismic hazard analyses. We synthesized late Holocene paleoseismic data from 20 trench sites along the central WFZ to evaluate earthquake rupture length and fault segmentation. For the youngest (<3 ka) and best-constrained earthquakes, differences in earthquake timing across prominent primary segment boundaries, especially for the most recent earthquakes on the north-central WFZ, are consistent with segment-controlled ruptures. However, broadly constrained earthquake times, dissimilar event times along the segments, the presence of smaller-scale (subsegment) boundaries, and areas of complex faulting permit partial-segment and multisegment (e.g., spillover) ruptures that are shorter (~20–40 km) or longer (~60–100 km) than the primary segment lengths (35–59 km). We report a segmented WFZ model that includes 24 earthquakes since ~7 ka and yields mean estimates of recurrence (1.1–1.3 kyr) and vertical slip rate (1.3–2.0 mm/yr) for the segments. However, additional rupture scenarios that include segment boundary spatial uncertainties, floating earthquakes, and multisegment ruptures are necessary to fully address epistemic uncertainties in rupture length. We compare the central WFZ to paleoseismic and historical surface ruptures in the Basin and Range Province and central Italian Apennines and conclude that displacement profiles have limited value for assessing the persistence of segment boundaries but can aid in interpreting prehistoric spillover ruptures. Our comparison also suggests that the probabilities of shorter and longer ruptures on the WFZ need to be investigated.

Utah↗

Time-averaged fluxes of lead and fallout radionuclides to sediments in Florida Bay

Recent, unmixed sediments from mud banks of central Florida Bay were dated using 210 Pb/ 226 Ra, and chronologies were verified by comparing sediment lead temporal records with Pb/Ca ratios in annual layers of coral ( Montastrea annularis ) located on the ocean side of the Florida Keys. Dates of sediment lead peaks (1978±2) accord with prior observations of a 6 year lag between the occurrence of maximum atmospheric lead in 1972 and peak coral lead in 1978. Smaller lags of 1–2 years occur between the maximum atmospheric radionuclide fallout and peaks in sediment temporal records of 137 Cs and Pu. Such lags are consequences of system time averaging (STA) in which atmospherically delivered particle-associated constituents accumulate and mix in a (sedimentary?) reservoir before transferring to permanent sediments and coral. STA model calculations, using time-dependent atmospheric inputs, produced optimized profiles in excellent accord with measured sediment 137 Cs, Pu, lead, and coral lead distributions. Derived residence times of these particle tracers (16±1, 15.7±0.7, 19±3, and 16±2 years, respectively) are comparable despite differences in sampling locations, in accumulating media, and in element loading histories and geochemical properties. For a 16 year weighted mean residence time, STA generates the observed 6 year lead peak lag. Evidently, significant levels of nondegradable, particle-associated contaminants can persist in Florida Bay for many decades following elimination of external inputs. Present results, in combination with STA model analysis of previously reported radionuclide profiles, suggest that decade-scale time averaging may occur widely in recent coastal marine sedimentary environments.

Florida↗

Case study of a full-scale evapotranspiration cover

The design, construction, and performance analyses of a 6.1ha evapotranspiration (ET) landfill cover at the semiarid U.S. Army Fort Carson site, near Colorado Springs, Colo. are presented. Initial water-balance model simulations, using literature reported soil hydraulic data, aided selection of borrow-source soil type(s) that resulted in predictions of negligible annual drainage (⩽1mm∕year). Final construction design was based on refined water-balance simulations using laboratory determined soil hydraulic values from borrow area natural soil horizons that were described with USDA soil classification methods. Cover design components included a 122cm thick clay loam (USDA), compaction ⩽80% of the standard Proctor maximum dry density (dry bulk density ∼1.3Mg/m3), erosion control measures, top soil amended with biosolids, and seeding with native grasses. Favorable hydrologic performance for a 5year period was documented by lysimeter-measured and Richards’-based calculations of annual drainage that were all <0.4mm∕year. Water potential data suggest that ET removed water that infiltrated the cover and contributed to a persistent driving force for upward flow and removal of water from below the base of the cover.

Journal of Geotechnical and Geoenvironmental Engin↗

'Is Ash Falling?', an online ashfall reporting tool in support of improved ashfall warnings and investigations of ashfall processes

The primary volcano hazard in Alaska is airborne ash, which endangers aircraft flying the busy North Pacific air routes and consequently affects global commerce. Downwind ashfall is also a significant threat to commerce, transportation and day-to-day activities in nearby Alaska communities. A web-enabled database, "Is Ash Falling?" has been developed to collect ashfall observations and encourage sample collections from the public during eruptions, enabling volcano observatory staff to concentrate on eruption response. Knowing the locations of filed ashfall reports improves public ashfall warnings and forecasts by providing on-the-ground checks for ash dispersion and fallout computer models and satellite imagery interpretation. Reports of ashfall are shared with emergency management agencies and the wider public. These reports also give scientists a more complete record of the amount, duration and other conditions of ashfall.

Journal of Applied Volcanology↗

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, situated in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the impact of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study we constructed the Strawberry Creek integrated hydrological model (SCIHM) to increase this understanding. The SCIHM is an integrated surface runoff and groundwater model that uses the coupled groundwater and surface-water flow model (GSFLOW), which is based on the integration of the precipitation-runoff modeling system (PRMS) and the modular groundwater flow model commonly called MODFLOW, version MODFLOW-2005 software to simulate surface runoff and infiltration and groundwater flow. The model has three layers, 263 rows, and 176 columns. The model area includes the Strawberry Creek and four adjacent watersheds. The PRMS model was calibrated using two streamflow gaging stations and the GSFLOW model was calibrated to reported spring diversion discharge and a sparse number of groundwater-level measurements. The SCIHM was run with and without diversions active and simulated streamflow was compared, finding that in the headwaters of Strawberry Creek about 35 percent of the diversion flow was captured from base flow.

EarthArXiv↗

Water-level database update for the Death Valley regional groundwater flow system, Nevada and California, 1907-2007

The water-level database for the Death Valley regional groundwater flow system in Nevada and California was updated. The database includes more than 54,000 water levels collected from 1907 to 2007, from more than 1,800 wells. Water levels were assigned a primary flag and multiple secondary flags that describe hydrologic conditions and trends at the time of the measurement and identify pertinent information about the well or water-level measurement. The flags provide a subjective measure of the relative accuracy of the measurements and are used to identify which water levels are appropriate for calculating head observations in a regional transient groundwater flow model. Included in the report appendix are all water-level data and their flags, selected well data, and an interactive spreadsheet for viewing hydrographs and well locations.

Data Series↗

Logistic and linear regression model documentation for statistical relations between continuous real-time and discrete water-quality constituents in the Kansas River, Kansas, July 2012 through June 2015

The Kansas River is a primary source of drinking water for about 800,000 people in northeastern Kansas. Source-water supplies are treated by a combination of chemical and physical processes to remove contaminants before distribution. Advanced notification of changing water-quality conditions and cyanobacteria and associated toxin and taste-and-odor compounds provides drinking-water treatment facilities time to develop and implement adequate treatment strategies. The U.S. Geological Survey (USGS), in cooperation with the Kansas Water Office (funded in part through the Kansas State Water Plan Fund), and the City of Lawrence, the City of Topeka, the City of Olathe, and Johnson County Water One, began a study in July 2012 to develop statistical models at two Kansas River sites located upstream from drinking-water intakes. Continuous water-quality monitors have been operated and discrete-water quality samples have been collected on the Kansas River at Wamego (USGS site number 06887500) and De Soto (USGS site number 06892350) since July 2012. Continuous and discrete water-quality data collected during July 2012 through June 2015 were used to develop statistical models for constituents of interest at the Wamego and De Soto sites. Logistic models to continuously estimate the probability of occurrence above selected thresholds were developed for cyanobacteria, microcystin, and geosmin. Linear regression models to continuously estimate constituent concentrations were developed for major ions, dissolved solids, alkalinity, nutrients (nitrogen and phosphorus species), suspended sediment, indicator bacteria ( Escherichia coli , fecal coliform, and enterococci), and actinomycetes bacteria. These models will be used to provide real-time estimates of the probability that cyanobacteria and associated compounds exceed thresholds and of the concentrations of other water-quality constituents in the Kansas River. The models documented in this report are useful for characterizing changes in water-quality conditions through time, characterizing potentially harmful cyanobacterial events, and indicating changes in water-quality conditions that may affect drinking-water treatment processes.

Kansas↗

Integrated watershed-scale response to climate change for selected basins across the United States

A study by the U.S. Geological Survey (USGS) evaluated the hydrologic response to different projected carbon emission scenarios of the 21st century using a hydrologic simulation model. This study involved five major steps: (1) setup, calibrate and evaluated the Precipitation Runoff Modeling System (PRMS) model in 14 basins across the United States by local USGS personnel; (2) acquire selected simulated carbon emission scenarios from the World Climate Research Programme's Coupled Model Intercomparison Project; (3) statistical downscaling of these scenarios to create PRMS input files which reflect the future climatic conditions of these scenarios; (4) generate PRMS projections for the carbon emission scenarios for the 14 basins; and (5) analyze the modeled hydrologic response. This report presents an overview of this study, details of the methodology, results from the 14 basin simulations, and interpretation of these results. A key finding is that the hydrological response of the different geographical regions of the United States to potential climate change may be different, depending on the dominant physical processes of that particular region. Also considered is the tremendous amount of uncertainty present in the carbon emission scenarios and how this uncertainty propagates through the hydrologic simulations.

Scientific Investigations Report↗

Linear regression model documentation for computing water-quality constituent concentrations or densities using continuous real-time water-quality data for the Kansas River above Topeka Weir at Topeka, Kansas, November 2018 through June 2021

The Kansas River and its associated alluvial aquifer provide drinking water to more than 950,000 people in northeastern Kansas. Water suppliers that rely on the Kansas River as a water-supply source use physical and chemical processes to treat and remove contaminants before public distribution. An early-notification system of changing water-quality conditions allows water suppliers to proactively make decisions that affect water treatment. The U.S. Geological Survey (USGS), in cooperation with the Kansas Water Office (funded in part through the Kansas Water Plan), the Kansas Department of Health and Environment, The Nature Conservancy, the City of Lawrence, the City of Manhattan, the City of Olathe, the City of Topeka, WaterOne, and Evergy, began collecting water-quality data at the Kansas River above Topeka Weir at Topeka, Kansas (USGS site 06888990, hereafter referred to as the “Topeka site”), during November 2018 to develop linear regression models that relate continuous in situ water-quality sensor measurements to discretely sampled water-quality constituent concentrations or densities. The addition of the Topeka site expanded an existing water-quality monitoring network, which included the upstream Kansas River at Wamego, Kans., and downstream Kansas River at De Soto, Kans., sites. Linear regression analysis was used to develop models that compute real-time concentrations or densities for total dissolved solids, major ions, hardness as calcium carbonate, nutrients (nitrogen and phosphorus species), chlorophyll a , total suspended solids, suspended sediment, and Escherichia coli at the Topeka site using data collected during November 2018 through June 2021. Water-quality constituent concentrations or densities computed from the models documented in this report are available at the USGS National Real-Time Water-Quality website (https://nrtwq.usgs.gov), are useful to the public for cultural and recreational purposes, and can be used to guide water-treatment processes, compare conditions with Federal and State water-quality criteria, and characterize changes in Kansas River water-quality conditions through time.

Kansas↗

One-Water Hydrologic Flow Model (MODFLOW-OWHM)

The One-Water Hydrologic Flow Model (MF-OWHM) is a MODFLOW-based integrated hydrologic flow model (IHM) that is the most complete version, to date, of the MODFLOW family of hydrologic simulators needed for the analysis of a broad range of conjunctive-use issues. Conjunctive use is the combined use of groundwater and surface water. MF-OWHM allows the simulation, analysis, and management of nearly all components of human and natural water movement and use in a physically-based supply-and-demand framework. MF-OWHM is based on the Farm Process for MODFLOW-2005 (MF-FMP2) combined with Local Grid Refinement (LGR) for embedded models to allow use of the Farm Process (FMP) and Streamflow Routing (SFR) within embedded grids. MF-OWHM also includes new features such as the Surface-water Routing Process (SWR), Seawater Intrusion (SWI), and Riparian Evapotrasnpiration (RIP-ET), and new solvers such as Newton-Raphson (NWT) and nonlinear preconditioned conjugate gradient (PCGN). This IHM also includes new connectivities to expand the linkages for deformation-, flow-, and head-dependent flows. Deformation-dependent flows are simulated through the optional linkage to simulated land subsidence with a vertically deforming mesh. Flow-dependent flows now include linkages between the new SWR with SFR and FMP, as well as connectivity with embedded models for SFR and FMP through LGR. Head-dependent flows now include a modified Hydrologic Flow Barrier Package (HFB) that allows optional transient HFB capabilities, and the flow between any two layers that are adjacent along a depositional or erosional boundary or displaced along a fault. MF-OWHM represents a complete operational hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for consumption by irrigated agriculture, but also of water used in urban areas and by natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, MF-OWHM has the unique set of capabilities to simulate and analyze historical, present, and future conjunctive-use conditions. MF-OWHM is especially useful for the analysis of agricultural water use where few data are available for pumpage, land use, or agricultural information. The features presented in this IHM include additional linkages with SFR, SWR, Drain-Return (DRT), Multi-Node Wells (MNW1 and MNW2), and Unsaturated-Zone Flow (UZF). Thus, MF-OWHM helps to reduce the loss of water during simulation of the hydrosphere and helps to account for “all of the water everywhere and all of the time.” In addition to groundwater, surface-water, and landscape budgets, MF-OWHM provides more options for observations of land subsidence, hydraulic properties, and evapotranspiration (ET) than previous models. Detailed landscape budgets combined with output of estimates of actual evapotranspiration facilitates linkage to remotely sensed observations as input or as additional observations for parameter estimation or water-use analysis. The features of FMP have been extended to allow for temporally variable water-accounting units (farms) that can be linked to land-use models and the specification of both surface-water and groundwater allotments to facilitate sustainability analysis and connectivity to the Groundwater Management Process (GWM). An example model described in this report demonstrates the application of MF-OWHM with the addition of land subsidence and a vertically deforming mesh, delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in demand caused by deficiency in supply, and changes in multi-aquifer pumpage caused by constraints imposed through the Farm Process and the MNW2 Package, and changes in surface water such as runoff, streamflow, and canal flows through SFR and SWR linkages.

Techniques and Methods↗

Power analysis and trend detection for water quality monitoring data. An application for the Greater Yellowstone Inventory and Monitoring Network

An important consideration for long term monitoring programs is determining the required sampling effort to detect trends in specific ecological indicators of interest. To enhance the Greater Yellowstone Inventory and Monitoring Network’s water resources protocol(s) (O’Ney 2006 and O’Ney et al. 2009 [under review]), we developed a set of tools to: (1) determine the statistical power for detecting trends of varying magnitude in a specified water quality parameter over different lengths of sampling (years) and different within-year collection frequencies (monthly or seasonal sampling) at particular locations using historical data, and (2) perform periodic trend analyses for water quality parameters while addressing seasonality and flow weighting. A power analysis for trend detection is a statistical procedure used to estimate the probability of rejecting the hypothesis of no trend when in fact there is a trend, within a specific modeling framework. In this report, we base our power estimates on using the seasonal Kendall test (Helsel and Hirsch 2002) for detecting trend in water quality parameters measured at fixed locations over multiple years. We also present procedures (R-scripts) for conducting a periodic trend analysis using the seasonal Kendall test with and without flow adjustment. This report provides the R-scripts developed for power and trend analysis, tutorials, and the associated tables and graphs. The purpose of this report is to provide practical information for monitoring network staff on how to use these statistical tools for water quality monitoring data sets.

Natural Resource Report↗