Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Environmental Modeling & Assessment”

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
USGS2014

Comparing species distribution models constructed with different subsets of environmental predictors

Aim To assess the usefulness of combining climate predictors with additional types of environmental predictors in species distribution models for range-restricted species, using common correlative species distribution modelling approaches. Location Florida, USA Methods We used five different algorithms to create distribution models for 14 vertebrate species, using seven different predictor sets: two with bioclimate predictors only, and five ‘combination’ models using bioclimate predictors plus ‘additional’ predictors from groups representing: human influence, land cover, extreme weather or noise (spatially random data).We use a linear mixed-model approach to analyse the effects of predictor set and algorithm on model accuracy, variable importance scores and spatial predictions. Results Regardless of modelling algorithm, no one predictor set produced significantly more accurate models than all others, though models including human influence predictors were the only ones with significantly higher accuracy than climate-only models. Climate predictors had consistently higher variable importance scores than additional predictors in combination models, though there was variation related to predictor type and algorithm. While spatial predictions varied moderately between predictor sets, discrepancies were significantly greater between modelling algorithms than between predictor sets. Furthermore, there were no differences in the level of agreement between binary ‘presence–absence’ maps and independent species range maps related to the predictor set used. Main conclusions Our results indicate that additional predictors have relatively minor effects on the accuracy of climate-based species distribution models and minor to moderate effects on spatial predictions. We suggest that implementing species distribution models with only climate predictors may provide an effective and efficient approach for initial assessments of environmental suitability.

Florida↗
USGS2009

Predictive Models of the Hydrological Regime of Unregulated Streams in Arizona

Three statistical models were developed by the U.S. Geological Survey in cooperation with the Arizona Department of Environmental Quality to improve the predictability of flow occurrence in unregulated streams throughout Arizona. The models can be used to predict the probabilities of the hydrological regime being one of four categories developed by this investigation: perennial, which has streamflow year-round; nearly perennial, which has streamflow 90 to 99.9 percent of the year; weakly perennial, which has streamflow 80 to 90 percent of the year; or nonperennial, which has streamflow less than 80 percent of the year. The models were developed to assist the Arizona Department of Environmental Quality in selecting sites for participation in the U.S. Environmental Protection Agency's Environmental Monitoring and Assessment Program. One model was developed for each of the three hydrologic provinces in Arizona - the Plateau Uplands, the Central Highlands, and the Basin and Range Lowlands. The models for predicting the hydrological regime were calibrated using statistical methods and explanatory variables of discharge, drainage-area, altitude, and location data for selected U.S. Geological Survey streamflow-gaging stations and a climate index derived from annual precipitation data. Models were calibrated on the basis of streamflow data from 46 stations for the Plateau Uplands province, 82 stations for the Central Highlands province, and 90 stations for the Basin and Range Lowlands province. The models were developed using classification trees that facilitated the analysis of mixed numeric and factor variables. In all three models, a threshold stream discharge was the initial variable to be considered within the classification tree and was the single most important explanatory variable. If a stream discharge value at a station was below the threshold, then the station record was determined as being nonperennial. If, however, the stream discharge was above the threshold, subsequent decisions were made according to the classification tree and explanatory variables to determine the hydrological regime of the reach as being perennial, nearly perennial, weakly perennial, or nonperennial. Using model calibration data, misclassification rates for each model were 17 percent for the Plateau Uplands, 15 percent for the Central Highlands, and 14 percent for the Basin and Range Lowlands models. The actual misclassification rate may be higher; however, the model has not been field verified for a full error assessment. The calibrated models were used to classify stream reaches for which the Arizona Department of Environmental Quality had collected miscellaneous discharge measurements. A total of 5,080 measurements at 696 sites were routed through the appropriate classification tree to predict the hydrological regime of the reaches in which the measurements were made. The predictions resulted in classification of all stream reaches as perennial or nonperennial; no reaches were predicted as nearly perennial or weakly perennial. The percentages of sites predicted as being perennial and nonperennial, respectively, were 77 and 23 for the Plateau Uplands, 87 and 13 for the Central Highlands, and 76 and 24 for the Basin and Range Lowlands.

Open-File Report↗
USGS2019

Designing flows to enhance ecosystem functioning in heavily altered rivers

More than a century of dam construction and water development in the western United States has led to extensive ecological alteration of rivers. Growing interest in improving river function is compelling practitioners to consider ecological restoration when managing dams and water extraction. We developed an Ecological Response Model (ERM) for the Cache la Poudre River, northern Colorado, USA, to illuminate effects of current and possible future water management and climate change. We used empirical data and modeled interactions among multiple ecosystem components to capture system-wide insights not possible with the unintegrated models commonly used in environmental assessments. The ERM results showed additional flow regime modification would further alter the structure and function of Poudre River aquatic and riparian ecosystems due to multiple and interacting stressors. Model predictions illustrated that specific peak flow magnitudes in spring and early summer are critical for substrate mobilization, dynamic channel morphology, and overbank flows, with strong subsequent effects on instream and riparian biota that varied seasonally and spatially, allowing exploration of nuanced management scenarios. Instream biological indicators benefitted from higher and more stable base flows and high peak flows, but stable base flows with low peak flows were only half as effective to increase indicators. Improving base flows while reducing peak flows, as currently proposed for the Cache la Poudre River, would further reduce ecosystem function. Modeling showed that even presently depleted annual flow volumes can achieve substantially different ecological outcomes in designed flow scenarios, while still supporting social demands. Model predictions demonstrated that implementing designed flows in a natural pattern, with attention to base and peak flows, may be needed to preserve or improve ecosystem function of the Poudre River. Improved regulatory policies would include preservation of ecosystem-level, flow-related processes and adaptive management when water development projects are considered.

Colorado↗
USGS2013

Forecasting the effects of coastal protection and restoration projects on wetland morphology in coastal Louisiana under multiple environmental uncertainty scenarios

Few landscape scale models have assessed the effects of coastal protection and restoration projects on wetland morphology while taking into account important uncertainties in environmental factors such as sea-level rise (SLR) and subsidence. In support of Louisiana's 2012 Coastal Master Plan, we developed a spatially explicit wetland morphology model and coupled it with other predictive models. The model is capable of predicting effects of protection and restoration projects on wetland area, landscape configuration, surface elevation, and soil organic carbon (SOC) storage under multiple environmental uncertainty scenarios. These uncertainty scenarios included variability in parameters such as eustatic SLR (ESLR), subsidence rate, and Mississippi River discharge. Models were run for a 2010–2060 simulation period. Model results suggest that under a “future-without-action” condition (FWOA), coastal Louisiana is at risk of losing between 2118 and 4677 km 2 of land over the next 50 years, but with protection and restoration projects proposed in the Master Plan, between 40% and 75% of that loss could be mitigated. Moreover, model results indicate that under a FWOA condition, SOC storage (to a depth of 1 m) could decrease by between 108 and 250 million metric tons, a loss of 12% to 30% of the total coastwide SOC, but with the Master Plan implemented, between 35% and 74% of the SOC loss could be offset. Long-term maintenance of project effects was best attained in areas of low SLR and subsidence, with a sediment source to support marsh accretion. Our findings suggest that despite the efficacy of restoration projects in mitigating losses in certain areas, net loss of wetlands in coastal Louisiana is likely to continue. Model results suggest certain areas may eventually be lost regardless of proposed restoration investment, and, as such, other techniques and strategies of adaptation may have to be utilized in these areas.

Louisiana↗
USGS2016

Influence of demography and environment on persistence in toad populations

Effective conservation of rare species requires an understanding of how potential threats affect population dynamics. Unfortunately, information about population demographics prior to threats (i.e., baseline data) is lacking for many species. Perturbations, caused by climate change, disease, or other stressors can lead to population declines and heightened conservation concerns. Boreal toads ( Anaxyrus boreas boreas ) have undergone rangewide declines due mostly to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd), with only a few sizable populations remaining in the southern Rocky Mountains, USA, that are disease-free. Despite the apparent region-wide occurrence of Bd, our focal populations in central Colorado were disease free over a 14-year capture-mark-recapture study until the recent discovery of Bd at one of the sites. We used recapture data and the Pradel reverse-time model to assess the influence of environmental and site-specific conditions on survival and recruitment. We then forecast changes in the toad populations with 2 growth models; one using an average lambda value to initiate the projection, and one using the most recent value to capture potential effects of the incursion of disease into the system. Adult survival was consistently high at the 3 sites, whereas recruitment was more variable and markedly low at 1 site. We found that active season moisture, active season length, and breeding shallows were important factors in estimating recruitment. Population growth models indicated a slight increase at 1 site but decreasing trends at the 2 other sites, possibly influenced by low recruitment. Insight into declining species management can be gained from information on survival and recruitment and how site-specific environmental factors influence these demographic parameters. Our data are particularly useful because they provide baseline data on demographics in populations before a disease outbreak and enhance our ability to detect changes in population parameters potentially caused by the disease.

Journal of Wildlife Management↗
USGS2008

Sediment ingestion rates in waterfowl (Anatidae) and their use in environmental risk assessment

When waterfowl (Anatidae) ingest sediment as they feed, they are exposed to the environmental contaminants in those sediments. The rate of ingestion may be key to assessing environmental risk. Rates of sediment ingestion were estimated as from <2% to 22% in 16 species of waterfowl collected in the northeastern United States. The piscivorous red-breasted merganser ( Mergus serrator ) ingested sediment at the lowest rate and the benthos-feeding canvasback ( Aythya valisineria ) at the highest rate. Sediment ingestion rates were related to diet and to the sediments where waterfowl fed. Waterfowl ingested the least sediment from hard-bottomed habitats with fast-moving water and ingested the most sediment from soft-bottomed areas with slow-moving water. Understanding the greater hazards from contaminants associated with low-flow habitats may help in prioritizing sites to be remediated. The tundra swan ( Cygnus columbianus ), which ingests sediment at an estimated 8.4% of its diet, dry weight, is suggested as a potential generic model for use in environmental risk assessments designed to protect waterfowl.

Integrated Environmental Assessment and Management↗
USGS2021

Modelling presence versus abundance for invasive species risk assessment

Aim Invasive species prevention and management can be guided by comparisons of invasion risk across space and among species. Species distribution models are widely used to assess invasion risk and typically estimate suitability for species presence. However, suitability for presence may not capture patterns of abundance and impact. We asked how models estimating suitability for presence versus suitability for abundance aligned in their implications for risk assessment. Location Western United States. Methods We developed ensembles of species distribution models for presence and for abundance for four invasive plants. We visualized the distribution of presence and abundance in environmental and geographic space and compared model outputs using criteria relevant for decision-making: a comparison of risk across management units for each species, and a ranking of risk among species for each management unit. Results We found good overall agreement between models of presence versus abundance in the relative risk across management units and among species. However, the area predicted to be suitable for invasive species presence was often substantially higher than the area predicted to be suitable for abundance, especially within uninvaded management units. Main conclusions Models of suitability for invasive species presence and abundance yielded similar assessments of relative risk in comparisons across space and species. In addition, we found patterns of presence and abundance in environmental space can guide modelling decisions and model interpretation. Suitability for abundance can improve relative risk assessment when abundance locations occupy a well-defined subset of the environmental space corresponding to presence. Where abundance locations occur throughout this environmental space, as was particularly striking for Taeniatherum caput-medusae, suitability for presence may better reflect risk of ongoing population increases and spread. This species is at risk of becoming abundant across a substantial portion of the western United States.

Diversity and Distributions↗
USGS2013

Selenium in ecosystems within the mountaintop coal mining and valley-fill region of southern West Virginia-assessment and ecosystem-scale modeling

Coal and associated waste rock are among environmental selenium (Se) sources that have the potential to affect reproduction in fish and aquatic birds. Ecosystems of southern West Virginia that are affected by drainage from mountaintop coal mines and valleys filled with waste rock in the Coal, Gauley, and Lower Guyandotte watersheds were assessed during 2010 and 2011. Sampling data from earlier studies in these watersheds (for example, Upper Mud River Reservoir) and other mining-affected watersheds also are included to assess additional hydrologic settings and food webs for comparison. Basin schematics give a comprehensive view of sampled species and Se concentration data specific to location and date. Food-web diagrams document the progression of Se trophic transfer across suspended particulate material, invertebrates, and fish for each site to serve as the basis for developing an ecosystem-scale model to predict Se exposure within the hydrologic conditions and food webs of southern West Virginia. This approach integrates a site-specific predator’s dietary exposure pathway into modeling to ensure an adequate link to Se toxicity and, thus, to species vulnerability. Site-specific fish abundance and richness data in streams documented various species of chub, shiner, dace, darters, bass, minnow, sunfish, sucker, catfish, and central stoneroller (Campostoma anomalum), mottled sculpin (Cottus bairdii), and least brook lamprey (Lampetra aepyptera). However, Se assessment species for streams, and hence, model species for streams, were limited to creek chub (Semotilus atromaculatus) and central stoneroller. Both of these species of fish are generally considered to have a high tolerance for environmental stress based on traditional comparative fish community assessment, with creek chub being present at all sites. Aquatic insects (mayfly, caddisfly, stonefly, dobsonfly, chironomid) were the main invertebrates sampled in streams. Collection of suspended particulate material acted as an integrator of organic-rich, fine-grained biomass present in streams. The base-case food web modeled for streams was suspended particulate material to aquatic insect to creek chub, with comparative modeling of a direct particulate-to-stoneroller food web. Model species for a reservoir setting were based on an earlier study of bluegill sunfish (Lepomis macrochirus), green sunfish (Lepomis cyanellus), and largemouth bass (Micropterus salmoides). Several reservoir food webs were considered based on a variety of invertebrates (insect, snail, clam). For stream and reservoir settings, predicted Se concentrations in exposure scenarios showed a high degree of correlation (r 2 = 0.91 for invertebrates and 0.75 for fish) with field observations of Se concentrations when modeling was initiated from suspended-particulate-material Se concentrations and model transfer parameters defined previously in the literature were used. These strong correlations validate the derived site-specific model and establish sufficient confidence that the predictions from the developed model can be quantitatively applied to the ecosystems in southern West Virginia. An application of modeling used a metric describing the partitioning of Se between particulate material and dissolved phases (K d ) to allow determination of a dissolved Se concentration that would be necessary to attain a site-specific Se fish body burden. The operationally defined K d quantifies the complex process of transformation at the base of a food web on a site-specific basis. The magnitude of this metric is known to vary with such factors as Se speciation, particulate-material type, and hydrology. This application (1) ties dissolved Se concentrations to fish tissue concentrations; (2) allows consideration of different choices for intervening site-specific exposure steps that set Se bioaccumulation, partitioning, and bioavailability; and (3) generates implications for management decisions that define protection through different regulatory pathways and guidelines. The range of model outcomes accounts for critical sources of variability and establishes whether site and food-web characterization were adequate to represent the dynamics of the system with certainty. This is especially true in terms of particulate-material phases at the base of the food web and utilization of K d in different hydrologic settings. For streams, a range of field-derived K d ds were applied to food-web exposure scenarios within a framework of locational and hydrologic variables (area of stream basin; stream gradient and discharge) that may affect the magnitude of K d . Overlaying even a coarse temporal scale that acknowledges variability in stream dissolved Se and Se speciation, such as through seasonal derivation of K d , can substantially narrow model uncertainty. Modeling that constrains the place and time of greatest ecosystem Se sensitivity within a specified food web gives insight into Se risk and identifies controlling management alternatives within a watershed or stream basin. If there is a range of hydrologic settings, specificity is needed to establish a hierarchy of in-stream and off-stream habitats for a watershed approach that takes into account Se-enriched water moving through different K d and food web environments. If there is a range of predator vulnerabilities (measured as a combination of food-web Se biodynamics and response in Se toxicity tests) within the site-specific community of fish species to be protected, then choice of fish species is critical to protection because it determines the food web and, hence, the magnitude of biotransfer through which Se is modeled. Whether creek chub is representative of the vulnerability to Se of all fish species encountered within the study-site ecosystems will require additional species-specific data and analysis. A range of site-specific scenarios illustrated here set model outcomes, but the final quantitative evaluation of alternatives and their implications will be those generated through choices and guidance formulated by state and other agencies in their decisionmaking processes. Proposed additions and refinements to the ecosystem-scale site-specific approach developed here include consideration of: measurement of temporally matched pairs of dissolved and suspended-particulate-material Se concentrations across a broader range of stream sites to expand the stream K d database and to test the representativeness of a suspended-particulate-material sample within a stream; characterization of different phases of particulate material across seasons to better define the base of the food web and connect to invertebrate feeding; refinement of model assumptions concerning dietary preferences and composition for fish to develop additional trophic transfer factors (TTFs) (for example, calculation of TTFinvertebrate composite for mixed diets); expansion of modeling of fish species and their food webs to include Se-vulnerable species; temporal characterization of a predator’s life cycle and habitat use as additional model layers to integrate with Se biodynamics in streams; investigation of the effect of stream gradient on K d based on a finer scale than presented here in terms of such variables as residence time, watershed dilution, and physical habitat attributes (for example, amount of ponding versus run or riffle within a stream); and linkage to discharge through use of stream gaging to record variability and enable model organization within water-year types and discharge seasons. Investigating the presence and variability of prey and predator species in demographically open systems such as streams also is key to model outcomes given the overall environmental stressors (for example, general landscape change, food-web disruption, recolonization potential) imposed on the composition of biological communities in coal mining and valley-fill affected watersheds

West Virginia↗
USGS2017

Virginia flow-ecology modeling results—An initial assessment of flow reduction effects on aquatic biota

Background The U.S. Geological Survey (USGS), in cooperation with the Virginia Department of Environmental Quality (DEQ), reviewed a previously compiled set of linear regression models to assess their utility in defining the response of the aquatic biological community to streamflow depletion. As part of the 2012 Virginia Healthy Watersheds Initiative (HWI) study conducted by Tetra Tech, Inc., for the U.S. Environmental Protection Agency (EPA) and Virginia DEQ, a database with computed values of 72 hydrologic metrics, or indicators of hydrologic alteration (IHA), 37 fish metrics, and 64 benthic invertebrate metrics was compiled and quality assured. Hydrologic alteration was represented by simulation of streamflow record for a pre-water-withdrawal condition (baseline) without dams or developed land, compared to the simulated recent-flow condition (2008 withdrawal simulation) including dams and altered landscape to calculate a percent alteration of flow. Biological samples representing the existing populations represent a range of alteration in the biological community today. For this study, all 72 IHA metrics, which included more than 7,272 linear regression models, were considered. This extensive dataset provided the opportunity for hypothesis testing and prioritization of flow-ecology relations that have the potential to explain the effect(s) of hydrologic alteration on biological metrics in Virginia streams.

Virginia↗
USGS2025

Habitat and predator influences on the spatial ecology of nine-banded armadillos

Mesopredator suppression has implications for community structure, biodiversity, and ecosystem function, but mesopredators with physical defenses may not avoid apex predators. We investigated nine-banded armadillos ( Dasypus novemcinctus ) in southwestern Oklahoma (USA) to evaluate if a species with physical defenses was influenced by a dominant predator, the coyote ( Canis latrans ). We sampled nine-banded armadillos and coyotes with motion-activated cameras. We used single-species and conditional two-species occupancy models to assess the influences of environmental factors and coyotes on nine-banded armadillo occurrence and site-use intensity (i.e., detection). We used camera-based detections to characterize the diel activity of each species and their overlap. Nine-banded armadillo occupancy was greater at sites closer to cover, with lower slopes, and further from water, whereas coyote space use was greater at higher elevations; both species were positively associated with recent burns. Nine-banded armadillo occurrence was not influenced by coyotes, but site-use intensity was suppressed by the presence of coyotes. Nine-banded armadillos (strictly nocturnal) and coyotes (predominantly nocturnal) had a high overlap in summer diel activity. Nine-banded armadillos are ecosystem engineers but are often considered a threat to species of concern and/or a nuisance. Thus, understanding the role of interspecific interactions on nine-banded armadillos has important implications for conservation and management.

Oklahoma↗
USGS2025

Human perturbations to mercury in global rivers

Mercury compounds are potent neurotoxins that pose threats to human health, primarily through fish consumption. Rivers, critical for drinking water and food supply, have seen rapid increases in mercury concentrations and export to coastal margins since the Industrial Revolution (~1850). However, patterns of these changes remain understudied, limiting assessments of environmental policies. Here, we develop a global model to simulate preindustrial riverine total mercury and assess human perturbations by comparing it to present-day conditions. We find that global rivers transported ~390 megagrams annually of mercury to the oceans in the preindustrial era, with spatial variability. Human activities have elevated riverine mercury budgets by two to three times in the present day. Establishing a baseline riverine mercury level, our findings reveal rapid responses of riverine mercury to human perturbations and could be used to inform targets for global riverine mercury restoration. Total riverine mercury concentrations could also be used as indicators to comprehensively understand the effectiveness of mercury pollution governance.

Science Advances↗
USGS2014

Temporal patterns in the foraging behavior of sea otters in Alaska

Activity time budgets in apex predators have been proposed as indicators of population status relative to resource limitation or carrying capacity. We used archival time-depth recorders implanted in 15 adult female and 4 male sea otters (Enhydra lutris) from the northernmost population of the species, Prince William Sound, Alaska, USA, to examine temporal patterns in their foraging behavior. Sea otters that we sampled spent less time foraging during summer (females 8.8 hr/day, males 7.9 hr/day) than other seasons (females 10.1&ndash;10.5 hr/day, males 9.2&ndash;9.5 hr/day). Both sexes showed strong preferences for diurnal foraging and adjusted their foraging effort in response to the amount of available daylight. One exception to this diurnal foraging mode occurred after females gave birth. For approximately 3 weeks post-partum, females switched to nocturnal foraging, possibly in an effort to reduce the risk of predation by eagles on newborn pups. We used multilevel mixed regression models to assess the contribution of several biological and environmental covariates to variation in the daily foraging effort of parous females. In the random effects only model, 87% of the total variation in foraging effort was within-otter variation. The relatively small among-otter variance component (13%) indicates substantial consistency in the foraging effort of sea otters in this northern population. In the top 3 models, 17% of the within-otter variation was explained by reproductive stage, day length, wind speed, air temperature and a wind speed &times; air temperature interaction. This study demonstrates the potential importance of environmental and reproductive effects when using activity budgets to assess population status relative to carrying capacity.

Alaska↗
USGS2015

An evaluation of fish behavior upstream of the water temperature control tower at Cougar Dam, Oregon, using acoustic cameras, 2013

This report describes the initial year of a 2-year study to determine the feasibility of using acoustic cameras to monitor fish movements to help inform decisions about fish passage at Cougar Dam near Springfield, Oregon. Specifically, we used acoustic cameras to measure fish presence, travel speed, and direction adjacent to the water temperature control tower in the forebay of Cougar Dam during the spring (May, June, and July) and fall (September, October, and November) of 2013. Cougar Dam is a high-head flood-control dam, and the water temperature control tower enables depth-specific water withdrawals to facilitate adjustment of water temperatures released downstream of the dam. The acoustic cameras were positioned at the upstream entrance of the tower to monitor free-ranging subyearling and yearling-size juvenile Chinook salmon ( Oncorhynchus tshawytscha ). Because of the large size discrepancy, we could distinguish juvenile Chinook salmon from their predators, which enabled us to measure predators and prey in areas adjacent to the entrance of the tower. We used linear models to quantify and assess operational and environmental factors&mdash;such as time of day, discharge, and water temperature&mdash;that may influence juvenile Chinook salmon movements within the beam of the acoustic cameras. Although extensive milling behavior of fish near the structure may have masked directed movement of fish and added unpredictability to fish movement models, the acoustic-camera technology enabled us to ascertain the general behavior of discrete size classes of fish. Fish travel speed, direction of travel, and counts of fish moving toward the water temperature control tower primarily were influenced by the amount of water being discharged through the dam.

Oregon↗
USGS1983

Evaluating environmental and economic consequences of alternative pest management strategies: results of modeling workshops

The U.S. Environmental Protection Agency (EPA) needs a comprehensive method to evaluate the human health and environmental effects of alternative agricultural pest management strategies. This project explored the utility of Adaptive Environmental Assessment (AEA) techniques for meeting this need. The project objectives were to produce models for environmental impact analysis, improve communications, identify research needs and data requirements, and demonstrate a process for resolving conflicts. The project was structured around the construction (in an initial 2 1/2-day workshop) and examination (in a second 2 1/2-day workshop) of a simulation model of a corn agroecosystem. The model conceptualized at the first workshop simulates the effect of corn agrecosystem decisions on crop production, economic returns, and environmental indicators. The model is composed of five interacting submodels: 1) a Production Strategies submodel which makes decisions concerning tillage, planting, fertilizer and pesticide applications, and harvest; 2) a Hydrology/Chemical Transport submodel which represents soil hydrology, erosion, and concentrations of fertilizers and pesticides in the soil, runoff, surface waters, and percolation; 3) a Vegetation submodel which simulates growth of agricultural crops (corns and soybeans) and weeds; 4) a Pests submodel which calculates pest population levels and resulting crop damage; and 5) an Environmental Effects submodel which calculates indicators of potential fish kills, human health effects, and wildlife habitat. The most persistent data gaps encountered in quantifying the model were coefficients to relate environmental consequences to alternative pest management strategies. While the model developed in the project is not yet accurate enough to be used for real-world decisions about the use of pesticides on corn, it does contain the basic structure upon which such a model could be built. More importantly at this stage of development, the project has shown that very complex systems can be modeled in short periods of time and that the process of building such models increases understanding among disciplinary specialists and between diverse institutional interests. This process can be useful to EPA as the agency cooperates with other institutions to meet its responsibilities in less costly ways. Activities at the second 2 1/2-day workshop included a review of the model, incorporation of necessary corrections, simulation of policy scenarios, and examination of techniques to address remaining institutional conflicts. Participants were divided into three groups representing environmental, production or industry, and regulatory interests. Each group developed scenarios that would be most appealing to their particular interest and the scenarios were simulated by the agroecosystem computer model. Negotiators from each of the interest groups decided whether a hypothetical herbicide should be relabeled and if certain restrictions should be imposed on its use. Other participants functioned as experts and consultants on caucus teams. A solution to the hypothetical problem was successfully negotiated. Workshop participants and project staff agreed that the model and processes developed during the project should be used in training students, extension specialists, farmers, researchers, and chemical producers in collaborative problem solving methods. More productive research can be planned, and more realistic models of complex systems can be built in this way. More importantly, greater trust of decisionmakers in computer models, better understanding by technical experts about disciplines other than their own, and improved cooperation between institutional interests can be achieved. This trust, understanding, and cooperation are critical ingredients in solving problems that are too complex to be resolved by independent disciplinary activity and unilateral decision authority.

Report↗
USGS2020

Understanding and documenting the scientific basis of selenium ecological protection in support of site-specific guidelines development for Lake Koocanusa, Montana, U.S.A., and British Columbia, Canada

Modeling of ecosystems is a part of the U.S. Environmental Protection Agency’s protocol for developing site-specific selenium guidelines for protection of aquatic life. Selenium as an environmental contaminant is known to bioaccumulate and cause reproductive effects in fish and wildlife. Here we apply a modeling methodology—ecosystem-scale selenium modeling—to understand and document the scientific basis for predicting and validating ecological protection for Lake Koocanusa, a transboundary reservoir between Montana and British Columbia. A comprehensive set of site-specific data compiled from public databases (Federal, State, and Provincial) and reports by Teck Coal Ltd., is available in a companion U.S. Geological Survey data release. The tissue guideline used within modeling here to assess protection is the U.S. Environmental Protection Agency’s national selenium guideline for whole-body fish (dry weight); however, other numeric values for a whole-body guideline or other tissue types may be assumed if applicable tissue-to-tissue conversion factors are available. We consider the report assembled here as a working document that presents a model that can effectively address and structure the needs of (1) scientific understanding in representing the lake’s ecosystem and selenium biodynamics and (2) policy and management development during a decision-making process, but it is open to modification and updating as more ecologically detailed data become available. The approach brings together the main concerns involved in selenium toxicity: likelihood of high exposure, inherent species sensitivity, and close connectivity of ecosystem characteristics and behavioral ecology of predators. Detailed site-specific modeling equations are provided to document the linked factors that determine the responses of ecosystems to selenium. A series of scenarios quantifies the implications of choices of site-specific variables including food-web species, bioavailability of particulate material, and partitioning between the dissolved and particulate phases at the base of food webs. A gradient mapping tool applied to Lake Koocanusa provides a precedent for ecosystem-scale modeling of lakes by recognizing the importance of lake strata and hydrodynamics as components of modeling. Data requirements for ecosystem modeling, including ecological and hydrological process information fundamental to the dietary biodynamics of selenium in site-specific food webs, were assessed as a precursor to model validation for Lake Koocanusa. Understanding these relationships is necessary to connect modeling outcomes to reproductive effects and establish boundaries, in the case of Lake Koocanusa, for the influences of dam operation, fish-community viability, and its Clean Water Act impaired 303(d)-listing status on ecosystem function. We find that an assemblage of conditions affects the representation of Lake Koocanusa’s ecosystem within modeling scenarios but that the constructed gradient maps, mechanistic model, and associated bioaccumulation potentials portray and quantify the variables that are determinative to protection of predator species. Ecological and hydrological sorting of compiled individual data points on a site- and species-specific basis helps identify and address model uncertainties. Sources of uncertainty include (1) the scarcity of data for some environmental media compartments across time and locations, (2) the complexity of hydrodynamic conditions that can lead to seasonal ecological disconnects such as in selenium partitioning from water into particulates, and (3) the functional status of Lake Koocanusa’s ecosystem because of cumulative effects of various environmental stresses (for example, fish-community changes, flow regime changes, parasites, gonadal dysfunction, and increasing mining input-selenium concentrations since 1984). To this last point, it is important to determine where Lake Koocanusa is in an impairment-restoration cycle so as not to base protection on survivor bias, the maintenance of a currently degraded ecosystem, or normalized toxicity. In a broader context, one of the overall consequences of revised selenium regulations is that their derivation is now dependent on being able to define and understand the status of the ecosystem on which protection is based.

Montana, British Columbia↗
USGS2020

Effect of stream permanence on predation risk of lotic crayfish by riparian predators

Given the importance of crayfish in stream ecosystems, gaining insight into the role of stream permanence in maintaining predator–prey interactions is critical. Our objectives were to determine the influence of stream permanence and season on crayfish predation and assess the role of stream permanence and crayfish density on the presence of predators, while accounting for imperfect detection. We conducted surveys of crayfish density, mammalian scat, and environmental variables within 10 intermittent and 10 permanent streams in the Ozark Highlands. We used occupancy modeling to assess the relationship between predator presence, crayfish density, and environmental variables. Stream permanence did not play a role in determining relative frequency of occurrence or percent volume of crayfish prey in mammalian diets. However, percent volume and relative frequency of crayfish prey found in scats differed by season, with both highest in spring and summer. The relative frequency and percent volume of fish prey showed a significant interaction of season by stream permanence, which may be the first instance of this observation. Procyon lotor (Raccoon) had the highest detection probability ( p = 0.39), whereas Neovison vison (American Mink; p = 0.15) and Lontra canadensis (River Otter; p = 0.03) had low detection probabilities. Further study into predator–prey interactions in the context of hydrology, particularly when related to imperiled groups like freshwater crayfishes, is needed since climate change is expected to alter hydrologic patterns.

Arkansas, Missouri↗
USGS2012

Automating calibration, sensitivity and uncertainty analysis of complex models using the R package Flexible Modeling Environment (FME): SWAT as an example

Parameter optimization and uncertainty issues are a great challenge for the application of large environmental models like the Soil and Water Assessment Tool (SWAT), which is a physically-based hydrological model for simulating water and nutrient cycles at the watershed scale. In this study, we present a comprehensive modeling environment for SWAT, including automated calibration, and sensitivity and uncertainty analysis capabilities through integration with the R package Flexible Modeling Environment (FME). To address challenges (e.g., calling the model in R and transferring variables between Fortran and R) in developing such a two-language coupling framework, 1) we converted the Fortran-based SWAT model to an R function (R-SWAT) using the RFortran platform, and alternatively 2) we compiled SWAT as a Dynamic Link Library (DLL). We then wrapped SWAT (via R-SWAT) with FME to perform complex applications including parameter identifiability, inverse modeling, and sensitivity and uncertainty analysis in the R environment. The final R-SWAT-FME framework has the following key functionalities: automatic initialization of R, running Fortran-based SWAT and R commands in parallel, transferring parameters and model output between SWAT and R, and inverse modeling with visualization. To examine this framework and demonstrate how it works, a case study simulating streamflow in the Cedar River Basin in Iowa in the United Sates was used, and we compared it with the built-in auto-calibration tool of SWAT in parameter optimization. Results indicate that both methods performed well and similarly in searching a set of optimal parameters. Nonetheless, the R-SWAT-FME is more attractive due to its instant visualization, and potential to take advantage of other R packages (e.g., inverse modeling and statistical graphics). The methods presented in the paper are readily adaptable to other model applications that require capability for automated calibration, and sensitivity and uncertainty analysis.

Environmental Modelling and Software↗