USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modeling and Using Context”

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 37 records · Page 2Linked to original sources
USGS2019

Using social-context matching to improve spatial function-transfer performance for cultural ecosystem service models

Recreational and aesthetic enjoyment of public lands is increasing across a wide range of activities, highlighting the need to assess and adapt management to accommodate these uses. Despite a growing number of studies on mapping cultural ecosystem services, most are local-scale assessments that rely on costly and time-consuming primary data collection. As a result, the availability of spatial information on non-market values associated with cultural ecosystem services (social values) remains limited. Spatial function transfer, if it could be justified for social-value models, would expedite the development of social-value information and promote its more regular inclusion in ecosystem service assessments. We used survey data from six national forests in Colorado and Wyoming to explore the potential for transferring cultural ecosystem service models between forests and specifically to test the hypothesis that transfer performance increases with social-context similarity between transferring and receiving areas. Results confirm this relationship but fall just short of being able to predict with certainty when transferred models will meet the minimum performance criterion needed for defensible use by managers. Social values are highly variable and can be difficult to predict, but our results suggest that with the right combination of indicators spatial function transfer can become a defensible means of generating social-value information when primary data collection is not feasible.

Colorado, Wyoming↗
USGS2025

Simulated soundscapes and transfer learning boost the performance of acoustic classifiers under data scarcity

1. The biodiversity crisis necessitates spatially extensive methods to monitor multiple taxonomic groups for evidence of change in response to evolving environmental conditions. Programs that combine passive acoustic monitoring and machine learning are increasingly used to meet this need. These methods require large, annotated datasets, which are time-consuming and expensive to produce, creating potential barriers to adoption in data- and funding-poor regions. Recently released pre-trained avian acoustic classification models provide opportunities to reduce the need for manual labelling and accelerate the development of new acoustic classification algorithms through transfer learning. Transfer learning is a strategy for developing algorithms under data scarcity that uses pre-trained models from related tasks to adapt to new tasks. 2. Our primary objective was to develop a transfer learning strategy using the feature embeddings of a pre-trained avian classification model to train custom acoustic classification models in data-scarce contexts. We used three annotated avian acoustic datasets to test whether transfer learning and soundscape simulation-based data augmentation could substantially reduce the annotated training data necessary to develop performant custom acoustic classifiers. We also conducted a sensitivity analysis for hyperparameter choice and model architecture. We then assessed the generalizability of our strategy to increasingly novel non-avian classification tasks. 3. With as few as two training examples per class, our soundscape simulation data augmentation approach consistently yielded new classifiers with improved performance relative to the pre-trained classification model and transfer learning classifiers trained with other augmentation approaches. Performance increases were evident for three avian test datasets, including single-class and multi-label contexts. We observed that the relative performance among our data augmentation approaches varied for the avian datasets and nearly converged for one dataset when we included more training examples. 4. We demonstrate an efficient approach to developing new acoustic classifiers leveraging open-source sound repositories and pre-trained networks to reduce manual labelling. With very few examples, our soundscape simulation approach to data augmentation yielded classifiers with performance equivalent to those trained with many more examples, showing it is possible to reduce manual label-ling while still achieving high-performance classifiers and, in turn, expanding the potential for passive acoustic monitoring to address rising biodiversity monitoring needs.

Methods in Ecology and Evolution↗
USGS2016

Sediment accumulation in prairie wetlands under a changing climate: The relative roles of landscape and precipitation

Sediment accumulation threatens the viability and hydrologic functioning of many naturally formed depressional wetlands across the interior regions of North America. These wetlands provide many ecosystem services and vital habitats for diverse plant and animal communities. Climate change may further impact sediment accumulation rates in the context of current land use patterns. We estimated sediment accretion in wetlands within a region renowned for its large populations of breeding waterfowl and migrant shorebirds and examined the relative roles of precipitation and land use context in the sedimentation process. We modeled rates of sediment accumulation from 1971 through 2100 using the Revised Universal Soil Loss Equation (RUSLE) with a sediment delivery ratio and the Unit Stream Power Erosion Deposition model (USPED). These models predicted that by 2100, 21–33 % of wetlands filled completely with sediment and 27–46 % filled by half with sediments; estimates are consistent with measured sediment accumulation rates in the region reported by empirical studies. Sediment accumulation rates were strongly influenced by size of the catchment, greater coverage of tilled landscape within the catchment, and steeper slopes. Conservation efforts that incorporate the relative risk of infilling of wetlands with sediments, thus emphasizing areas of high topographic relief and large watersheds, may benefit wetland-dependent biota.

Prairie Pothole region↗
USGS1998

Estimate of aquifer properties by numerically simulating ground-water/surface-water interactions, Fort Wainwright, Alaska

MODFLOW, a finite-difference model of ground-water flow, was used to simulate the flow of water between the aquifer and the Chena River at Fort Wainwright, Alaska. The model was calibrated by comparing simulated ground-water hydrographs to those recorded in wells during periods of fluctuating river levels. The best fit between simulated and observed hydrographs occurred for the following: 20 feet per day for vertical hydraulic conductivity, 400 feet per day for horizontal hydraulic conductivity, 1:20 for anisotropy (vertical to horizontal hydraulic conductivity), and 350 per feet for riverbed conductance. These values include a 30 percent adjustment for geometry effects. The estimated values for hydraulic conductivities of the alluvium are based on assumed values of 0.25 for specific yield and 0.000001 per foot for specific storage of the alluvium; the values assumed for bedrock are 0.1 foot per day horizontal hydraulic conductivity, 0.005 foot per day vertical hydraulic conductivity, and 0.0000001 per foot for specific storage. The resulting diffusivity for the alluvial aquifer is 1,600 feet per day. The estimated values of these hydraulic properties are nearly proportional to the assumed value of specific yield. These values were not found to be sensitive to the assumed values for bedrock. The hydrologic parameters estimated using the cross-sectional model are only valid when taken in context with the other values (both estimated and assumed) used in this study. The model simulates horizontal and vertical flow directions near the river during periods of varying river stage. This information is useful for interpreting bank-storage effects, including the flow of contaminants in the aquifer near the river.

Alaska↗
USGS2018

Combining conflicting Bayesian models to develop paleoseismic records—An example from the Wasatch Fault Zone, Utah

Bayesian statistical analyses of paleoseismic data result in the probabilistic determination of earthquake times using geochronological data evaluated in the context of a stratigraphic model. However, a fundamental problem in paleoseismology is how to use the Bayesian approach to model sparse and/or conflicting geochronological datasets, such as those derived from sites exhibiting episodic sedimentary and pedogenic processes in moderate‐ to high‐energy environments (e.g., a normal‐faulted alluvial fan). Using paleoseismic data for the Corner Canyon site on the Salt Lake City segment of the Wasatch fault zone (Utah), we develop an approach by which multiple Bayesian models are combined to generate an earthquake history at a site. This approach accommodates mutually exclusive interpretations of the geochronological data and thereby limits the influence of sparse data, stratigraphically inconsistent ages, or a single, subjective model interpretation. For the Corner Canyon site, we integrate four OxCal Bayesian models to generate a chronology of six events between ∼ 4.8 "> ∼ 4.8 and ∼ 0.5    ka "> ∼ 0.5 ka ⁠ . Late Holocene (post‐5 ka) mean recurrence and vertical slip‐rate estimates are ∼ 0.9    ky "> ∼ 0.9 ky (0.7–1.0 ky; 95% confidence) and 1.1    mm / yr "> 1.1 mm / yr ( ⁠ 0.8 – 1.7    mm / yr "> 0.8 – 1.7 mm / yr range), respectively. Although our method increases the uncertainty in the timing of individual earthquakes, it more objectively accounts for potential geochronological errors and different interpretations of stratigraphic age control. By relaxing the need to select a single age model, our approach yields more accurate earthquake‐timing results that will better facilitate evaluations of along‐fault event correlation and earthquake rupture length.

Utah↗
USGS1999

A GIS modeling method applied to predicting forest songbird habitat

We have developed an approach for using “presence” data to construct habitat models. Presence data are those that indicate locations where the target organism is observed to occur, but that cannot be used to define locations where the organism does not occur. Surveys of highly mobile vertebrates often yield these kinds of data. Models developed through our approach yield predictions of the amount and the spatial distribution of good-quality habitat for the target species. This approach was developed primarily for use in a GIS context; thus, the models are spatially explicit and have the potential to be applied over large areas. Our method consists of two primary steps. In the first step, we identify an optimal range of values for each habitat variable to be used as a predictor in the model. To find these ranges, we employ the concept of maximizing the difference between cumulative distribution functions of (1) the values of a habitat variable at the observed presence locations of the target organism, and (2) the values of that habitat variable for all locations across a study area. In the second step, multivariate models of good habitat are constructed by combining these ranges of values, using the Boolean operators “and” and “or.” We use an approach similar to forward stepwise regression to select the best overall model. We demonstrate the use of this method by developing species-specific habitat models for nine forest-breeding songbirds (e.g., Cerulean Warbler, Scarlet Tanager, Wood Thrush) studied in southern Ohio. These models are based on species’ microhabitat preferences for moisture and vegetation characteristics that can be predicted primarily through the use of abiotic variables. We use slope, land surface morphology, land surface curvature, water flow accumulation downhill, and an integrated moisture index, in conjunction with a land-cover classification that identifies forest/nonforest, to develop these models. The performance of these models was evaluated with an independent data set. Our tests showed that the models performed better than random at identifying where the birds occurred and provided useful information for predicting the amount and spatial distribution of good habitat for the birds we studied. In addition, we generally found positive correlations between the amount of habitat, as predicted by the models, and the number of territories within a given area. This added component provides the possibility, ultimately, of being able to estimate population sizes. Our models represent useful tools for resource managers who are interested in assessing the impacts of alternative management plans that could alter or remove habitat for these birds.

Ecological Applications↗
USGS2001

. Ecological conceptual models: a framework and case study on ecosystem management for South Florida sustainability

The Everglades and South Florida ecosystems are the focus of national and international attention because of their current degraded and threatened state. Ecological risk assessment, sustainability and ecosystem and adaptive management principles and processes are being used nationally as a decision and policy framework for a variety of types of ecological assessments. The intent of this study is to demonstrate the application of these paradigms and principles at a regional scale. The effects-directed assessment approach used in this study consists of a retrospective, eco-epidemiological phase to determine the causes for the current conditions and a prospective predictive risk-based assessment using scenario analysis to evaluate future options. Embedded in these assessment phases is a process that begins with the identification of goals and societal preferences which are used to develop an integrated suite of risk-based and policy relevant conceptual models. Conceptual models are used to illustrate the linkages among management (societal) actions, environmental stressors, and societal/ecological effects, and provide the basis for developing and testing causal hypotheses. These models, developed for a variety of landscape units and their drivers, stressors, and endpoints, are used to formulate hypotheses to explain the current conditions. They are also used as the basis for structuring management scenarios and analyses to project the temporal and spatial magnitude of risk reduction and system recovery. Within the context of recovery, the conceptual models are used in the initial development of performance criteria for those stressors that are determined to be most important in shaping the landscape, and to guide the use of numerical models used to develop quantitative performance criteria in the scenario analysis. The results will be discussed within an ecosystem and adaptive management framework that provides the foundation for decision making.

Science of the Total Environment↗
USGS2007

OPR-PPR, a computer program for assessing data importance to model predictions using linear statistics

The OPR-PPR program calculates the Observation-Prediction (OPR) and Parameter-Prediction (PPR) statistics that can be used to evaluate the relative importance of various kinds of data to simulated predictions. The data considered fall into three categories: (1) existing observations, (2) potential observations, and (3) potential information about parameters. The first two are addressed by the OPR statistic; the third is addressed by the PPR statistic. The statistics are based on linear theory and measure the leverage of the data, which depends on the location, the type, and possibly the time of the data being considered. For example, in a ground-water system the type of data might be a head measurement at a particular location and time. As a measure of leverage, the statistics do not take into account the value of the measurement. As linear measures, the OPR and PPR statistics require minimal computational effort once sensitivities have been calculated. Sensitivities need to be calculated for only one set of parameter values; commonly these are the values estimated through model calibration. OPR-PPR can calculate the OPR and PPR statistics for any mathematical model that produces the necessary OPR-PPR input files. In this report, OPR-PPR capabilities are presented in the context of using the ground-water model MODFLOW-2000 and the universal inverse program UCODE_2005. The method used to calculate the OPR and PPR statistics is based on the linear equation for prediction standard deviation. Using sensitivities and other information, OPR-PPR calculates (a) the percent increase in the prediction standard deviation that results when one or more existing observations are omitted from the calibration data set; (b) the percent decrease in the prediction standard deviation that results when one or more potential observations are added to the calibration data set; or (c) the percent decrease in the prediction standard deviation that results when potential information on one or more parameters is added.

Techniques and Methods↗
USGS2012

Explaining differences between bioaccumulation measurements in laboratory and field data through use of a probabilistic modeling approach

In the regulatory context, bioaccumulation assessment is often hampered by substantial data uncertainty as well as by the poorly understood differences often observed between results from laboratory and field bioaccumulation studies. Bioaccumulation is a complex, multifaceted process, which calls for accurate error analysis. Yet, attempts to quantify and compare propagation of error in bioaccumulation metrics across species and chemicals are rare. Here, we quantitatively assessed the combined influence of physicochemical, physiological, ecological, and environmental parameters known to affect bioaccumulation for 4 species and 2 chemicals, to assess whether uncertainty in these factors can explain the observed differences among laboratory and field studies. The organisms evaluated in simulations including mayfly larvae, deposit-feeding polychaetes, yellow perch, and little owl represented a range of ecological conditions and biotransformation capacity. The chemicals, pyrene and the polychlorinated biphenyl congener PCB-153, represented medium and highly hydrophobic chemicals with different susceptibilities to biotransformation. An existing state of the art probabilistic bioaccumulation model was improved by accounting for bioavailability and absorption efficiency limitations, due to the presence of black carbon in sediment, and was used for probabilistic modeling of variability and propagation of error. Results showed that at lower trophic levels (mayfly and polychaete), variability in bioaccumulation was mainly driven by sediment exposure, sediment composition and chemical partitioning to sediment components, which was in turn dominated by the influence of black carbon. At higher trophic levels (yellow perch and the little owl), food web structure (i.e., diet composition and abundance) and chemical concentration in the diet became more important particularly for the most persistent compound, PCB-153. These results suggest that variation in bioaccumulation assessment is reduced most by improved identification of food sources as well as by accounting for the chemical bioavailability in food components. Improvements in the accuracy of aqueous exposure appear to be less relevant when applied to moderate to highly hydrophobic compounds, because this route contributes only marginally to total uptake. The determination of chemical bioavailability and the increase in understanding and qualifying the role of sediment components (black carbon, labile organic matter, and the like) on chemical absorption efficiencies has been identified as a key next steps.

Integrated Environmental Assessment and Management↗
USGS2003

Diagnostic tools for mixing models of stream water chemistry

Mixing models provide a useful null hypothesis against which to evaluate processes controlling stream water chemical data. Because conservative mixing of end‐members with constant concentration is a linear process, a number of simple mathematical and multivariate statistical methods can be applied to this problem. Although mixing models have been most typically used in the context of mixing soil and groundwater end‐members, an extension of the mathematics of mixing models is presented that assesses the “fit” of a multivariate data set to a lower dimensional mixing subspace without the need for explicitly identified end‐members. Diagnostic tools are developed to determine the approximate rank of the data set and to assess lack of fit of the data. This permits identification of processes that violate the assumptions of the mixing model and can suggest the dominant processes controlling stream water chemical variation. These same diagnostic tools can be used to assess the fit of the chemistry of one site into the mixing subspace of a different site, thereby permitting an assessment of the consistency of controlling end‐members across sites. This technique is applied to a number of sites at the Panola Mountain Research Watershed located near Atlanta, Georgia.

Water Resources Research↗
USGS2002

Development and evaluation of a local grid refinement method for block-centered finite-difference groundwater models using shared nodes

A new method of local grid refinement for two-dimensional block-centered finite-difference meshes is presented in the context of steady-state groundwater-flow modeling. The method uses an iteration-based feedback with shared nodes to couple two separate grids. The new method is evaluated by comparison with results using a uniform fine mesh, a variably spaced mesh, and a traditional method of local grid refinement without a feedback. Results indicate: (1) The new method exhibits quadratic convergence for homogeneous systems and convergence equivalent to uniform-grid refinement for heterogeneous systems. (2) Coupling the coarse grid with the refined grid in a numerically rigorous way allowed for improvement in the coarse-grid results. (3) For heterogeneous systems, commonly used linear interpolation of heads from the large model onto the boundary of the refined model produced heads that are inconsistent with the physics of the flow field. (4) The traditional method works well in situations where the better resolution of the locally refined grid has little influence on the overall flow-system dynamics, but if this is not true, lack of a feedback mechanism produced errors in head up to 3.6% and errors in cell-to-cell flows up to 25%. ?? 2002 Elsevier Science Ltd. All rights reserved.

Advances in Water Resources↗
USGS2019

Evaluating and using existing models to map probable suitable habitat for rare plants to inform management of multiple-use public lands in the California desert

Multiple-use public lands require balancing diverse resource uses and values across landscapes. In the California desert, there is strong interest in renewable energy development and important conservation concerns. The Bureau of Land Management recently completed a land-use plan for the area that provides protection for modeled suitable habitat for multiple rare plants. Three sets of habitat models were commissioned for plants of conservation concern as part of the planning effort. The Bureau of Land Management then needed to determine which model or combination of models to use to implement plan requirements. Our goals were to: 1) develop a process for evaluating the existing habitat models and 2) use the evaluation results to map probable and potential suitable habitat. We developed a method for evaluating the construction (input data and methods) and performance of existing models and applied it to 88 habitat models for 43 rare plant species. We also developed a process for mapping probable and potential suitable habitat based on the existing models; potential habitat maps are intended only to guide future field surveys. We were able to map probable suitable habitat for 26 of the 43 species and potential suitable habitat for 41 species. Forty percent of the project area contains probable suitable habitat for at least one species (43,338 km2), with much of that habitat (43%) occurring on lands managed by the Bureau of Land Management. Lands prioritized for renewable energy development contain 3% of the habitat modeled as suitable for at least one species. Our products can be used by agencies to review proposed projects and plan future plant surveys and by developers to target sites likely to minimize conflicts with rare plant conservation goals. Our methods can be broadly applied to understand and quantify the defensibility of models used in conservation and regulatory contexts.

California↗
USGS2021

Multiple feedbacks due to biotic interactions across trophic levels can lead to persisten novel conditions that hinder restoration

Unlike traditional successional theory, Alternate Stable Equilibrium (ASE) theory posits that more than one community state is possible in a single environment, depending on the order that species arrive. ASE theory is often invoked in management situations where initial stressors have been removed, but native-dominated communities are not returning to degraded areas. Fundamental to this theory is the assumption that equilibria are maintained by positive feedbacks between colonizers and their environment. While ASE has been relatively well studied in aquatic ecosystems, more complex terrestrial systems offer multiple challenges, including species interactions across trophic levels that can lead to multiple feedbacks. Here, we discuss ASE theory as it applies to terrestrial, invaded ecosystems, and detail a case study from Hawaii that exemplifies how species interactions can favour the persistence of invaders, and how an understanding of interactions and feedbacks can be used to guide management. Our system includes intact native-dominated mesic forest and areas cleared for pasture, planted with non-native grasses, and later planted with a monoculture of a native nitrogen-fixing tree in an effort to restore forests. We discuss interactions between birds, understorey fruiting native species, understorey non-native grasses, soils and bryophytes in separate feedback mechanisms, and explain our efforts to identify which of these feedbacks is most important to address in a management context. Finally, we suggest that using models can help overcome some of the challenges that terrestrial ecosystems pose when studying ASE.

Book chapter↗
USGS2006

Combining multistate capture-recapture data with tag recoveries to estimate demographic parameters

Matrix population models that allow an animal to occupy more than one state over time are important tools for population and evolutionary ecologists. Definition of state can vary, including location for metapopulation models and breeding state for life history models. For populations whose members can be marked and subsequently re-encountered, multistate mark-recapture models are available to estimate the survival and transition probabilities needed to construct population models. Multistate models have proved extremely useful in this context, but they often require a substantial amount of data and restrict estimation of transition probabilities to those areas or states subjected to formal sampling effort. At the same time, for many species, there are considerable tag recovery data provided by the public that could be modeled in order to increase precision and to extend inference to a greater number of areas or states. Here we present a statistical model for combining multistate capture-recapture data (e.g., from a breeding ground study) with multistate tag recovery data (e.g., from wintering grounds). We use this method to analyze data from a study of Canada Geese (Branta canadensis) in the Atlantic Flyway of North America. Our analysis produced marginal improvement in precision, due to relatively few recoveries, but we demonstrate how precision could be further improved with increases in the probability that a retrieved tag is reported.

Ecology↗
USGS2019

Selecting a landscape model for natural resource management applications

Purpose of Review: Climate change and associated ecological impacts have challenged many conventional, observation-based approaches for predicting ecosystem and landscape responses to natural resource management. Complex spatial ecological models provide powerful, flexible tools which managers and others can use to make inferences about management impacts on future, no-analog landscape conditions. However, land managers who wish to use ecosystem and landscape models for natural resource applications are faced with the difficult task of deciding among many models that differ in important ways. Here, we summarize a process to aid managers in the selection of an appropriate model for natural resource management. Recent Findings: To guide management planning, scientifically credible information on how landscapes will respond to management actions under changing climate is required. Landscape models are increasingly used in a management context to evaluate of impacts of changing climate and interacting stressors on ecosystems and to test effects of alternative management options on desired conditions. However, the wide range of available models makes selection of appropriate and viable models a complex process. Summary: We present a series of preliminary steps to define critical scales of time, space, and ecological organization to guide an experimental design for a modeling project and then list a set of criteria for selecting a landscape or ecological model. Material presented includes the preliminary steps (crafting modeling objective, designing modeling project), organizational concerns (resources available, expertise on hand, timelines), and modeling details (complexity, design, documentation) of model selection.

Current Landscape Ecology Reports↗
USGS2010

Flood hazard awareness and hydrologic modelling at Ambos Nogales, United States–Mexico border

Appropriate land‐use, watershed‐management, and flood‐attenuation plans are critical in the cross‐border urban environment known collectively as Ambos Nogales. This paper summarizes methodologies for predicting the watershed response associated with land‐use change within a spatial and temporal context through the use of a hydrological model in a cross‐border setting. The KINEROS2 model is implemented via the Automated Geospatial Watershed Assessment 2.0 geographic information system interface to evaluate the watershed of Nogales, Arizona, and Nogales, Sonora, Mexico, to assess flood vulnerability by quantifying volumes of runoff and peak flow, based on alternative land‐use scenarios. Cross‐border geospatial data acquisition and input to models are described. Discussions about the KINEROS2 model results identify flood‐prone areas, simulate the impact of land‐use change, and evaluate the impact of potential flood‐control interventions in the Ambos Nogales watershed. Products from this research are being used in a comprehensive plan for sustainable development in Ambos Nogales.

Ambos Nogales watershed↗
USGS2015

Updated numerical model with uncertainty assessment of 1950-56 drought conditions on brackish-water movement within the Edwards aquifer, San Antonio, Texas

In 2010, the U.S. Geological Survey, in cooperation with the San Antonio Water System, began a study to assess the brackish-water movement within the Edwards aquifer (more specifically the potential for brackish-water encroachment into wells near the interface between the freshwater and brackish-water transition zones, referred to in this report as the transition-zone interface) and effects on spring discharge at Comal and San Marcos Springs under drought conditions using a numerical model. The quantitative targets of this study are to predict the effects of higher-than-average groundwater withdrawals from wells and drought-of-record rainfall conditions of 1950–56 on (1) dissolved-solids concentration changes at production wells near the transition-zone interface, (2) total spring discharge at Comal and San Marcos Springs, and (3) the groundwater head (head) at Bexar County index well J-17. The predictions of interest, and the parameters implemented into the model, were evaluated to quantify their uncertainty so the results of the predictions could be presented in terms of a 95-percent credible interval. The model area covers the San Antonio and Barton Springs segments of the Edwards aquifer; the history-matching effort was focused on the San Antonio segment. A previously developed diffuse-flow model of the Edwards aquifer, which forms the basis for the model in this assessment, is primarily based on a conceptualization in which flow in the aquifer is predominately through a network of numerous small fractures and openings. Primary updates to this model include an extension of the active area downdip, a conversion to an 8-layer SEAWAT variable-density flow and transport model to simulate dissolved-solids concentration effects on water density, history matching to 1999–2009 conditions, and parameter estimation in a highly parameterized context using automated methods in PEST (a model-independent Parameter ESTimation code). In addition to the best-fit parameter values derived from history matching, the uncertainty of model parameters was also estimated by using linear uncertainty analysis. Comparison of “prior” (before history matching) and “posterior” (after history matching) variances of parameters indicate that the information within the observation dataset used for history matching informs many parameters. The concentration threshold parameters were well-informed by the observation dataset as their posterior distributions were much narrower than their prior distributions. The transition-zone scaling parameters of hydraulic conductivity, effective porosity, and specific storage were all informed by the observation dataset, as evidenced by the difference between the prior and posterior variances. Saline-zone scaling parameters, alternatively, were not informed by the observation dataset for effective porosity and specific storage. Resulting posterior drier-month, wetter-month, and annual recharge multiplier parameter variances are important to understanding how well recharge is estimated and implemented within the model. The shifts of the posterior distributions left and right indicate that there were zones where less or more water was needed in the model. The widths of the distributions were not decreased substantially, indicating that many of the best-fit recharge parameters are not statistically different from the initial values specified in the history-matching effort. Recharge from rainfall is the driving force behind groundwater flow and heads in the aquifer; therefore, an increase in understanding of this process would benefit model development by potentially decreasing the uncertainty of this parameter. The history-matching effort was most helpful in informing the parameters in the model that control discharge at springs, namely, the spring orifice (drain) altitude and drain conductance parameters for each spring. The uncertainty assessment of the predictive model (a hypothetical recurrence of 1950–56 drought conditions and higher-than-average groundwater withdrawals from wells) provided insights into the potential effects of these conditions on dissolved-solids concentration changes at production wells near the transition-zone interface, discharges at Comal and San Marcos Springs, and heads at Bexar County index well J-17. Results at the 25 production wells near the transition-zone interface indicate that the uncertainty of model input parameters based on expert knowledge yielded an upper bound of the 95-percent credible interval of dissolved-solids concentrations that exceeds the secondary drinking water standards of 1,000 milligrams per liter (mg/L) of the Texas Commission on Environmental Quality (TCEQ) for many wells. However, the history-matching process provided key information to inform prediction-sensitive model parameters and therefore, contributed to a substantial decrease of the upper bound of the 95-percent credible interval to below the secondary drinking water standards. Reductions in dissolved-solids concentration changes were on the order of 400 mg/L to 1,300 mg/L. The reduction in uncertainty in regards to this prediction implies that this prediction of dissolved-solids concentration change can be made with some certainty using this current model and that those parameters that control this prediction are informed by the observation dataset. Even though predictive uncertainty was reduced for this prediction, dissolved-solids concentration changes were still greater than zero, indicating a minimal increase in concentration at these 25 production wells during the 7-year simulation period is likely. However, this minimal concentration increase indicates a small potential for movement of the brackish-water transition zone near these wells during the 7-year simulation period of drought-ofrecord (1950–56) rainfall conditions with higher-than-average groundwater withdrawals by wells. Predictive results of total spring discharge during the 7-year period, as well as head predictions at Bexar County index well J-17, were much different than the dissolved-solids concentration change results at the production wells. These upper bounds are an order of magnitude larger than the actual prediction which implies that (1) the predictions of total spring discharge at Comal and San Marcos Springs and head at Bexar County index well J-17 made with this model are not reliable, and (2) parameters that control these predictions are not informed well by the observation dataset during historymatching, even though the history-matching process yielded parameters to reproduce spring discharges and heads at these locations during the history-matching period. Furthermore, because spring discharges at these two springs and heads at Bexar County index well J-17 represent more of a cumulative effect of upstream conditions over a larger distance (and longer time), many more parameters (with their own uncertainties) are potentially controlling these predictions than the prediction of dissolved-solids concentration change at the prediction wells, and therefore contributing to a large posterior uncertainty.

Texas↗
USGS2022

Modeling functional flows in California rivers

Environmental flows are critical to the recovery and conservation of freshwater ecosystems worldwide. However, estimating desired ranges of environmental flows across large, diverse landscapes is challenging. To advance protections of environmental flows for streams in California, USA, we developed a statewide modeling approach focused on functional components of the natural flow regime. Functional flow components in California streams—fall pulse flows, wet season peak flows and base flows, the spring flow recession, and dry season baseflows—support essential physical and ecological processes in riverine ecosystems. These functional flow components can be represented by functional flow metrics (FFMs) and quantified by their magnitude, timing, frequency, duration, and rate-of-change from daily streamflow records. After quantifying FFMs at reference-quality streamflow gages in California, we used machine-learning methods to estimate their natural range of values for all stream reaches in the state based on physical watershed characteristics and climatic factors. We found that the models performed well in predicting FFMs in streams across a diversity of landscape and climate contexts, according to several model performance criteria. Using the predicted FFM values, we established initial estimates of ecological flows that are expected to support critical functions and are broadly protective of ecosystem health. Modeling functional flows statewide offers a pathway for increasing the pace and scale of environmental flow protections in California and beyond.

California↗