Search USGS⌕ Search

SEARCH · Search USGS

Results for “Ecological Modelling”

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 937 records · Page 52Linked to original sources

Assessing conservation relevance of organism-environment relations using predicted changes in response variables

1. Organism–environment models are used widely in conservation. The degree to which they are useful for informing conservation decisions – the conservation relevance of these relations – is important because lack of relevance may lead to misapplication of scarce conservation resources or failure to resolve important conservation dilemmas. Even when models perform well based on model fit and predictive ability, conservation relevance of associations may not be clear without also knowing the magnitude and variability of predicted changes in response variables. 2. We introduce a method for evaluating the conservation relevance of organism–environment relations that employs confidence intervals for predicted changes in response variables. The confidence intervals are compared to a preselected magnitude of change that marks a threshold (trigger) for conservation action. To demonstrate the approach, we used a case study from the Chihuahuan Desert involving relations between avian richness and broad-scale patterns of shrubland. We considered relations for three winters and two spatial extents (1- and 2-km-radius areas) and compared predicted changes in richness to three thresholds (10%, 20% and 30% change). For each threshold, we examined 48 relations. 3. The method identified seven, four and zero conservation-relevant changes in mean richness for the 10%, 20% and 30% thresholds respectively. These changes were associated with major (20%) changes in shrubland cover, mean patch size, the coefficient of variation for patch size, or edge density but not with major changes in shrubland patch density. The relative rarity of conservation-relevant changes indicated that, overall, the relations had little practical value for informing conservation decisions about avian richness. 4. The approach we illustrate is appropriate for various response and predictor variables measured at any temporal or spatial scale. The method is broadly applicable across ecological environments, conservation objectives, types of statistical predictive models and levels of biological organization. By focusing on magnitudes of change that have practical significance, and by using the span of confidence intervals to incorporate uncertainty of predicted changes, the method can be used to help improve the effectiveness of conservation efforts.

Methods in Ecology and Evolution↗

Streamflow as a stressor: Disentangling hydrology and water quality impacts to characterize flow-ecology relationships for two stream assemblages across two southeastern landscapes

Disassociating the independent effects of flow and water quality on the ecology of flowing waters is an overarching goal in water resource science needed to improve the efficacy of watershed management. However, the interrelatedness of these gradients and their subsequent alteration due to land use change has constrained progress made on this front. The objective of this study was to use benthic macroinvertebrate and fish assemblage data to characterize flow-ecology relationships that were unchanged by water quality impacts across two southeastern landscapes in the USA to help detect ecological change driven by flow alteration. General linear latent models were used to identify taxa that were responsive to high or low flow metrics and water quality gradients. Bayesian hierarchical generalized additive models were then developed using these indicator taxa and three biological metrics to identify flow-specific relationships that were unaffected by water quality impacts. Three low flow-specific relationships were identified, illustrating how potential agricultural or urban impacts to hydrology reduced stream biological health. Importantly, flow-ecology relationships developed using indicator taxa in this study effectively captured hydrology-specific impacts while biological metrics typical of state monitoring and assessment programs did not. Therefore, developing flow-specific biological metrics is a critical step when developing management strategies targeting flow alteration. Implementing standardized frameworks such as the one characterized here can limit contradictory findings and improve streamflow enhancement and restoration project efficacy. These low flow-specific relationships will enhance managers' capacity to develop environmental flow standards, monitor their success, and better understand urban and agricultural impacts on stream assemblages.

Alabama↗

Tools for quantifying isotopic niche space and dietary variation at the individual and population level.

Ecologists are increasingly using stable isotope analysis to inform questions about variation in resource and habitat use from the individual to community level. In this study we investigate data sets from 2 California sea otter ( Enhydra lutris nereis ) populations to illustrate the advantages and potential pitfalls of applying various statistical and quantitative approaches to isotopic data. We have subdivided these tools, or metrics, into 3 categories: IsoSpace metrics, stable isotope mixing models, and DietSpace metrics. IsoSpace metrics are used to quantify the spatial attributes of isotopic data that are typically presented in bivariate (e.g., δ 13 C versus δ 15 N) 2-dimensional space. We review IsoSpace metrics currently in use and present a technique by which uncertainty can be included to calculate the convex hull area of consumers or prey, or both. We then apply a Bayesian-based mixing model to quantify the proportion of potential dietary sources to the diet of each sea otter population and compare this to observational foraging data. Finally, we assess individual dietary specialization by comparing a previously published technique, variance components analysis, to 2 novel DietSpace metrics that are based on mixing model output. As the use of stable isotope analysis in ecology continues to grow, the field will need a set of quantitative tools for assessing isotopic variance at the individual to community level. Along with recent advances in Bayesian-based mixing models, we hope that the IsoSpace and DietSpace metrics described here will provide another set of interpretive tools for ecologists.

Journal of Mammalogy↗

Cohesive framework for modeling plant cover class data

The study of plant distribution and abundance is a fundamental pursuit in ecology and conservation biology. Measuring plant abundance by visually assessing percent cover and recording a cover class is a common field method that yields ordinal data. Statistical models for ordinal data exist but entail cumbersome interpretations and sometimes restrictive assumptions. We propose a Bayesian hierarchical framework for analysing cover class data that allows for linking ordinal observations to a latent beta distribution and accounts for zero inflation. Harnessing a latent beta distribution supports interpreting changes in abundance in terms of mean percent cover rather than odds ratios of cumulative cover classes as for cumulative link models. The zero augmentation allows for simultaneous inferences on both occurrence (distribution) and abundance. We show how our model can account for true and false zeros, misclassification of cover classes, multiple species and hierarchical sampling designs, using empirical examples and simulations. Simulated observation errors, when ignored, led to models overestimating abundance and underestimating occurrence. Based on simulations, we found no substantial difference between mean percent cover estimates when analyzing ordinal cover classes versus continuous percent cover as the response. Our empirical datasets displayed high probability of detection (>0.85 on average for all species), likely due to the sampling design used and training of observers. Probability of occurrence was slightly underestimated for bare ground, Artemisia tridentata , Elycap medusae , and Poa secunda using a model that ignored imperfect detection. Estimated mean percent cover was not substantially impacted by ignoring measurement error for five plant species and bare ground. Our modelling framework for cover class data allows for an explicit separation of distribution from abundance and, importantly, allows for interpreting species–environment relationships in terms of variation in mean percent cover as compared to cumulative odds ratios. The beta distribution inherently accommodates heteroscedasticity and skewness, statistical properties that are a consequence of spatially aggregated patterns common to plant survey data. Recording cover classes provides a reliable, efficient way to measure plants and our simulations suggest little loss of information compared to assuming continuous percent cover. We provide JAGS and Stan model code for implementation.

Methods in Ecology and Evolution↗

Modeling habitat dynamics accounting for possible misclassification

Land cover data are widely used in ecology as land cover change is a major component of changes affecting ecological systems. Landscape change estimates are characterized by classification errors. Researchers have used error matrices to adjust estimates of areal extent, but estimation of land cover change is more difficult and more challenging, with error in classification being confused with change. We modeled land cover dynamics for a discrete set of habitat states. The approach accounts for state uncertainty to produce unbiased estimates of habitat transition probabilities using ground information to inform error rates. We consider the case when true and observed habitat states are available for the same geographic unit (pixel) and when true and observed states are obtained at one level of resolution, but transition probabilities estimated at a different level of resolution (aggregations of pixels). Simulation results showed a strong bias when estimating transition probabilities if misclassification was not accounted for. Scaling-up does not necessarily decrease the bias and can even increase it. Analyses of land cover data in the Southeast region of the USA showed that land change patterns appeared distorted if misclassification was not accounted for: rate of habitat turnover was artificially increased and habitat composition appeared more homogeneous. Not properly accounting for land cover misclassification can produce misleading inferences about habitat state and dynamics and also misleading predictions about species distributions based on habitat. Our models that explicitly account for state uncertainty should be useful in obtaining more accurate inferences about change from data that include errors.

Landscape Ecology↗

Compilation of geologic, hydrologic, and ground-water flow modeling information for the Spokane Valley-Rathdrum Prairie aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources and Washington Department of Ecology compiled and described geologic, hydrologic, and ground-water flow modeling information about the Spokane Valley-Rathdrum Prairie (SVRP) aquifer in northern Idaho and northeastern Washington. Descriptions of the hydrogeologic framework, water-budget components, ground- and surface-water interactions, computer flow models, and further data needs are provided. The SVRP aquifer, which covers about 370 square miles including the Rathdrum Prairie, Idaho and the Spokane valley and Hillyard Trough, Washington, was designated a Sole Source Aquifer by the U.S. Environmental Protection Agency in 1978. Continued growth, water management issues, and potential effects on water availability and water quality in the aquifer and in the Spokane and Little Spokane Rivers have illustrated the need to better understand and manage the region's water resources. The SVRP aquifer is composed of sand, gravel, cobbles, and boulders primarily deposited by a series of catastrophic glacial outburst floods from ancient Glacial Lake Missoula. The material deposited in this high-energy environment is coarser-grained than is typical for most basin-fill deposits, resulting in an unusually productive aquifer with well yields as high as 40,000 gallons per minute. In most places, the aquifer is bounded laterally by bedrock composed of granite, metasedimentary rocks, or basalt. The lower boundary of the aquifer is largely unknown except along the margins or in shallower parts of the aquifer where wells have penetrated its entire thickness and reached bedrock or silt and clay deposits. Based on surface geophysics, the thickness of the aquifer is about 500 ft near the Washington-Idaho state line, but more than 600 feet within the Rathdrum Prairie and more than 700 feet in the Hillyard trough based on drilling records. Depth to water in the aquifer is greatest in the northern Rathdrum Prairie (about 500 feet) and least near the city of Spokane along the Spokane River (less than about 50 feet). Ground-water flow is south from near the southern end of Lake Pend Oreille and Hoodoo Valley, through the Rathdrum Prairie, then west toward Spokane. In Spokane, the aquifer splits and water moves north through the Hillyard Trough as well as west through the Trinity Trough. From the Trinity Trough water flows north along the western arm of the aquifer. The aquifer's discharge area is along the Little Spokane River and near Long Lake, Washington. A compilation of estimates of water-budget components, including recharge (precipitation, irrigation, canal leakage, septic tank effluent, inflow from tributary basins, and flow from the Spokane River) and discharge (withdrawals from wells, flow to the Spokane and Little Spokane Rivers, evapotranspiration, and underflow to Long Lake) illustrates that these estimated values should be compared with caution due to several variables including the area and time period of interest as well as methods employed in making the estimates. Numerous studies have documented the dynamic ground-water and surface-water interaction between the SVRP aquifer and the Spokane and Little Spokane Rivers. Gains and losses vary throughout the year, as well as the locations of gains and losses. September 2004 streamflow measurements indicated that the upper reach of the Spokane River between Post Falls and downstream at Flora Road lost 321 cubic feet per second. A gain of 736 cubic feet per second was measured between the Flora Road site and downstream at Green Street Bridge. A loss of 124 cubic feet per second was measured for the reach between the Green Street Bridge and the Spokane River at Spokane gaging station. The river gained about 87 cubic feet per second between the Spokane River at Spokane gaging station and the TJ Meenach Bridge. Overall, the Spokane River gained about 284 cubic feet per second between the Post Falls,

Scientific Investigations Report↗

Estimated probability of arsenic in groundwater from bedrock aquifers in New Hampshire, 2011

Probabilities of arsenic occurrence in groundwater from bedrock aquifers at concentrations of 1, 5, and 10 micrograms per liter (µg/L) were estimated during 2011 using multivariate logistic regression. These estimates were developed for use by the New Hampshire Environmental Public Health Tracking Program. About 39 percent of New Hampshire bedrock groundwater was identified as having at least a 50 percent chance of containing an arsenic concentration greater than or equal to 1 µg/L. This compares to about 7 percent of New Hampshire bedrock groundwater having at least a 50 percent chance of containing an arsenic concentration equaling or exceeding 5 µg/L and about 5 percent of the State having at least a 50 percent chance for its bedrock groundwater to contain concentrations at or above 10 µg/L. The southeastern counties of Merrimack, Strafford, Hillsborough, and Rockingham have the greatest potential for having arsenic concentrations above 5 and 10 µg/L in bedrock groundwater. Significant predictors of arsenic in groundwater from bedrock aquifers for all three thresholds analyzed included geologic, geochemical, land use, hydrologic, topographic, and demographic factors. Among the three thresholds evaluated, there were some differences in explanatory variables, but many variables were the same. More than 250 individual predictor variables were assembled for this study and tested as potential predictor variables for the models. More than 1,700 individual measurements of arsenic concentration from a combination of public and private water-supply wells served as the dependent (or predicted) variable in the models. The statewide maps generated by the probability models are not designed to predict arsenic concentration in any single well, but they are expected to provide useful information in areas of the State that currently contain little to no data on arsenic concentration. They also may aid in resource decision making, in determining potential risk for private wells, and in ecological-level analysis of disease outcomes. The approach for modeling arsenic in groundwater could also be applied to other environmental contaminants that have potential implications for human health, such as uranium, radon, fluoride, manganese, volatile organic compounds, nitrate, and bacteria.

New Hampshire↗

Landscape characteristics influencing the genetic structure of greater sage-grouse within the stronghold of their range: a holistic modeling approach

Given the significance of animal dispersal to population dynamics and geographic variability, understanding how dispersal is impacted by landscape patterns has major ecological and conservation importance. Speaking to the importance of dispersal, the use of linear mixed models to compare genetic differentiation with pairwise resistance derived from landscape resistance surfaces has presented new opportunities to disentangle the menagerie of factors behind effective dispersal across a given landscape. Here, we combine these approaches with novel resistance surface parameterization to determine how the distribution of high- and low-quality seasonal habitat and individual landscape components shape patterns of gene flow for the greater sage-grouse ( Centrocercus urophasianus ) across Wyoming. We found that pairwise resistance derived from the distribution of low-quality nesting and winter, but not summer, seasonal habitat had the strongest correlation with genetic differentiation. Although the patterns were not as strong as with habitat distribution, multivariate models with sagebrush cover and landscape ruggedness or forest cover and ruggedness similarly had a much stronger fit with genetic differentiation than an undifferentiated landscape. In most cases, landscape resistance surfaces transformed with 17.33-km-diameter moving windows were preferred, suggesting small-scale differences in habitat were unimportant at this large spatial extent. Despite the emergence of these overall patterns, there were differences in the selection of top models depending on the model selection criteria, suggesting research into the most appropriate criteria for landscape genetics is required. Overall, our results highlight the importance of differences in seasonal habitat preferences to patterns of gene flow and suggest the combination of habitat suitability modeling and linear mixed models with our resistance parameterization is a powerful approach to discerning the effects of landscape on gene flow.

Ecology and Evolution↗

Bayesian cross-validation for model evaluation and selection, with application to the North American Breeding Bird Survey

The analysis of ecological data has changed in two important ways over the last 15 years. The development and easy availability of Bayesian computational methods has allowed and encouraged the fitting of complex hierarchical models. At the same time, there has been increasing emphasis on acknowledging and accounting for model uncertainty. Unfortunately, the ability to fit complex models has outstripped the development of tools for model selection and model evaluation: familiar model selection tools such as Akaike's information criterion and the deviance information criterion are widely known to be inadequate for hierarchical models. In addition, little attention has been paid to the evaluation of model adequacy in context of hierarchical modeling, i.e., to the evaluation of fit for a single model. In this paper, we describe Bayesian cross-validation, which provides tools for model selection and evaluation. We describe the Bayesian predictive information criterion and a Bayesian approximation to the BPIC known as the Watanabe-Akaike information criterion. We illustrate the use of these tools for model selection, and the use of Bayesian cross-validation as a tool for model evaluation, using three large data sets from the North American Breeding Bird Survey.

Ecology↗

Mapping landslide susceptibility over large regions with limited data

Landslide susceptibility maps indicate the spatial distribution of landslide likelihood. Modeling susceptibility over large or diverse terrains remains a challenge due to the sparsity of landslide data (mapped extent of known landslides) and the variability in triggering conditions. Several different data sampling strategies of landslide locations used to train a susceptibility model are used to mitigate this challenge. However, to our knowledge, no study has systematically evaluated how different sampling strategies alter a model's predictor effects (i.e., how a predictor value influences the susceptibility output) critical to explaining differences in model outputs. Here, we introduce a statistical framework that examines the variation in predictor effects and the model accuracy (measured using receiver operator characteristics) to highlight why certain sampling strategies are more effective than others. Specifically, we apply our framework to an array of logistic regression models trained on landslide inventories collected at sub-regional scales over four terrains across the United States. Results show significant variations in predictor effects depending on the inventory used to train the models. The inconsistent predictor effects cause low accuracies when testing models on inventories outside the domain of the training data. Grouping test and training sets according to physiographic and ecological characteristics, which are thought to share similar triggering mechanisms, does not improve model accuracy. We also show that using limited landslide data distributed uniformly over the entire modeling domain is better than using dense but spatially isolated data to train a model for applications over large regions.

California, Kentucky, North Carolina, West Virgini↗

Scenarios of land use and land cover change in the conterminous United States: Utilizing the special report on emission scenarios at ecoregional scales

Global environmental change scenarios have typically provided projections of land use and land cover for a relatively small number of regions or using a relatively coarse resolution spatial grid, and for only a few major sectors. The coarseness of global projections, in both spatial and thematic dimensions, often limits their direct utility at scales useful for environmental management. This paper describes methods to downscale projections of land-use and land-cover change from the Intergovernmental Panel on Climate Change's Special Report on Emission Scenarios to ecological regions of the conterminous United States, using an integrated assessment model, land-use histories, and expert knowledge. Downscaled projections span a wide range of future potential conditions across sixteen land use/land cover sectors and 84 ecological regions, and are logically consistent with both historical measurements and SRES characteristics. Results appear to provide a credible solution for connecting regionalized projections of land use and land cover with existing downscaled climate scenarios, under a common set of scenario-based socioeconomic assumptions.

Global Environmental Change↗

Incorporating citizen science data in spatially explicit integrated population models

Information about population abundance, distribution, and demographic rates is critical for understanding a species’ ecology and for effective conservation and management. To collect data over large spatial and temporal extents for such inferences, especially for species with low densities or wide distributions, citizen science can be an efficient approach. Integrated models have also emerged as an important methodology to estimate population parameters by combining multiple types of data, including citizen science data. We developed a spatially explicit integrated model that combines opportunistically collected presence–absence (PA) data, commonly collected in citizen science efforts, with systematically collected spatial capture–recapture (SCR) data, which are often limited to small spatial and temporal extents. We conducted single and multi‐season simulations with parameters informed by North American black bear ( Ursus americanus ) populations, to evaluate the influence of varying amounts of opportunistic PA data collected at larger spatial and temporal extents on the estimation of population‐level parameters. Integrating opportunistic PA data increased the precision and accuracy of posterior estimates of abundance, and survival and recruitment rates. In some cases, adding PA locations improved abundance estimates more than increasing PA detection probability. Posterior estimates were as precise and unbiased as when higher quality, but sparse, SCR data were available. We also applied the integrated model to SCR and citizen science PA data collected on black bears in New York, with results consistent with our simulations. Our findings indicate that citizen science in integrated models can be a cost‐efficient way to improve estimates of population parameters and increase the spatiotemporal extent of inference. Continued developments with integrated models and citizen science data will offer additional ways to improve our understanding of population structure and demographics.

Ecology↗

Application of an unstructured 3D finite volume numerical model to flows and salinity dynamics in the San Francisco Bay-Delta

A linked modeling approach has been undertaken to understand the impacts of climate and infrastructure on aquatic ecology and water quality in the San Francisco Bay-Delta region. The Delft3D Flexible Mesh modeling suite is used in this effort for its 3D hydrodynamics, salinity, temperature and sediment dynamics, phytoplankton and water-quality coupling infrastructure, and linkage to a habitat suitability model. The hydrodynamic model component of the suite is D-Flow FM, a new 3D unstructured finite-volume model based on the Delft3D model. In this paper, D-Flow FM is applied to the San Francisco Bay-Delta to investigate tidal, seasonal and annual dynamics of water levels, river flows and salinity under historical environmental and infrastructural conditions. The model is driven by historical winds, tides, ocean salinity, and river flows, and includes federal, state, and local freshwater withdrawals, and regional gate and barrier operations. The model is calibrated over a 9-month period, and subsequently validated for water levels, flows, and 3D salinity dynamics over a 2 year period. Model performance was quantified using several model assessment metrics and visualized through target diagrams. These metrics indicate that the model accurately estimated water levels, flows, and salinity over wide-ranging tidal and fluvial conditions, and the model can be used to investigate detailed circulation and salinity patterns throughout the Bay-Delta. The hydrodynamics produced through this effort will be used to drive affiliated sediment, phytoplankton, and contaminant hindcast efforts and habitat suitability assessments for fish and bivalves. The modeling framework applied here will serve as a baseline to ultimately shed light on potential ecosystem change over the current century.

California↗

A predictive model to inform adaptive management of double-crested cormorants and fisheries in Michigan

The proliferation of double-crested cormorants (DCCOs; Phalacrocorax auritus ) in North America has raised concerns over their potential negative impacts on game, cultured and forage fishes, island and terrestrial resources, and other colonial water birds, leading to increased public demands to reduce their abundance. By combining fish surplus production and bird functional feeding response models, we developed a deterministic predictive model representing bird–fish interactions to inform an adaptive management process for the control of DCCOs in multiple colonies in Michigan. Comparisons of model predictions with observations of changes in DCCO numbers under management measures implemented from 2004 to 2012 suggested that our relatively simple model was able to accurately reconstruct past DCCO population dynamics. These comparisons helped discriminate among alternative parameterizations of demographic processes that were poorly known, especially site fidelity. Using sensitivity analysis, we also identified remaining critical uncertainties (mainly in the spatial distributions of fish vs. DCCO feeding areas) that can be used to prioritize future research and monitoring needs. Model forecasts suggested that continuation of existing control efforts would be sufficient to achieve long-term DCCO control targets in Michigan and that DCCO control may be necessary to achieve management goals for some DCCO-impacted fisheries in the state. Finally, our model can be extended by accounting for parametric or ecological uncertainty and including more complex assumptions on DCCO–fish interactions as part of the adaptive management process.

Natural Resource Modeling↗

Modeling effects of disturbance across life history strategies of stream fishes

A central goal of population ecology is to establish linkages between life history strategy, disturbance, and population dynamics. Globally, disturbance events such as drought and invasive species have dramatically impacted stream fish populations and contributed to sharp declines in freshwater biodiversity. Here, we used RAMAS-Metapop to construct stage-based demographic metapopulation models for stream fishes with periodic, opportunistic, and equilibrium life history strategies and assessed their responses to disturbance scenarios that approximated the effects of invasion, drought, and the additive effects of both disturbances. Our models indicated that populations respond differentially to disturbance based on life history strategy. Equilibrium strategists were best able to deal with simulated invasion. Periodic strategists did well under lower levels of drought, whereas opportunistic strategists outperformed other life histories under extreme seasonal drought. When we modeled additive effects scenarios, these disturbances interacted synergistically, dramatically increasing terminal extinction risk for all three life history strategies. Modeling exercises that examine broad life history categories can help to answer fundamental ecological questions about the relationship between disturbance resilience and life history, as well as help managers to develop generalized conservation strategies when species-specific data are lacking. Our results indicate that life history strategy is a fundamental determinant of population trajectories, and that disturbances can interact synergistically to dramatically impact extinction outcomes.

Oecologia↗

The fitting of general force-of-infection models to wildlife disease prevalence data

Researchers and wildlife managers increasingly find themselves in situations where they must deal with infectious wildlife diseases such as chronic wasting disease, brucellosis, tuberculosis, and West Nile virus. Managers are often charged with designing and implementing control strategies, and researchers often seek to determine factors that influence and control the disease process. All of these activities require the ability to measure some indication of a disease's foothold in a population and evaluate factors affecting that foothold. The most common type of data available to managers and researchers is apparent prevalence data. Apparent disease prevalence, the proportion of animals in a sample that are positive for the disease, might seem like a natural measure of disease's foothold, but several properties, in particular, its dependency on age structure and the biasing effects of disease-associated mortality, make it less than ideal. In quantitative epidemiology, the a??force of infection,a?? or infection hazard, is generally the preferred parameter for measuring a disease's foothold, and it can be viewed as the most appropriate way to a??adjusta?? apparent prevalence for age structure. The typical ecology curriculum includes little exposure to quantitative epidemiological concepts such as cumulative incidence, apparent prevalence, and the force of infection. The goal of this paper is to present these basic epidemiological concepts and resulting models in an ecological context and to illustrate how they can be applied to understand and address basic epidemiological questions. We demonstrate a practical approach to solving the heretofore intractable problem of fitting general force-of-infection models to wildlife prevalence data using a generalized regression approach. We apply the procedures to Mycobacterium bovis (bovine tuberculosis) prevalence in bison (Bison bison) in Wood Buffalo National Park, Canada, and demonstrate strong age dependency in the force of infection as well as an increased mortality hazard in positive animals.

Ecology↗

Integrated science strategy for assessing and monitoring water availability and migratory birds for terminal lakes across the Great Basin, United States

Executive Summary In 2022, the U.S. Geological Survey (USGS) established the Saline Lake Ecosystems Integrated Water Availability Assessment (IWAAs) to monitor and assess the hydrology of terminal lakes in the Great Basin and the migratory birds and other wildlife dependent on those habitats. Scientists from across the USGS (with specialties in water quantity, water quality, limnology, avian biology, data science, landscape ecology, and science communication) formed the Saline Lake Ecosystems IWAAs Team. The team has developed this regional strategic science plan to guide data collection and assessment activities at terminal lakes in the Great Basin. The U.S. Congress requested the USGS to establish the Saline Lake Ecosystems IWAAs in response to historically low water levels at terminal lakes and associated wetlands across the Great Basin. Not all Great Basin terminal lakes have high salinity; however, all terminal lakes occur in endorheic, closed, basins with no surface-water outflow. Low lake levels across the Great Basin are the result of increased water use for agriculture and municipalities, drought conditions, and a warming climate. Great Basin terminal lake water extents have decreased by as much as 90 percent over the last 150 years, and terminal lake wetlands have decreased in area by as much as 47 percent since 1984. Lake elevations and wetland areas are primarily supported by freshwater inputs from snowmelt feeding upgradient rivers, streams, and springs. These freshwater inputs have been severely reduced because of continued and increased surface-water diversions and surface-water capture through groundwater pumping for agriculture, mining, and public supply as well as unprecedented drought conditions and warming temperatures related to climate change. Water quality, specifically salinity, is highly variable for terminal lakes of the Great Basin, and this variability is a result of the balance between freshwater inflow and evaporation. Variability of salinity at each of the terminal lakes can be affected by lake morphology, hydrogeologic features of the basin, annual variability in weather patterns, and changes in upgradient water use. Hypersaline terminal lakes provide abundant food resources such as brine shrimp and brine flies that support nesting and migrating birds. The density and composition of invertebrates are closely tied to lake salinity. Increased salinity can exceed the tolerance of invertebrates, severely limiting their biomass. In contrast, decreased salinity can lead to altered invertebrate community composition, reducing the abundance of optimal avian prey resources. Great Basin terminal lake ecosystems, including open-water and adjacent aquatic and terrestrial environments, provide resources necessary to sustain many animal populations throughout the year. Although a variety of taxa use terminal lakes, these ecosystems are of acute importance for the millions of migratory waterbirds (for example, shorebirds, wading birds, and waterfowl) dependent on the network of terminal lakes and their associated wetlands. Migratory birds transiting the Pacific and Central Flyways use Great Basin terminal lake ecosystems throughout the year to feed, nest, and transit between wintering and breeding ranges. As such, successful conservation of birds and their habitats requires coordinated management of water and habitats across the Great Basin network of terminal lakes and wetlands. The linkages between water availability and ecosystem vulnerability of terminal lakes in the Great Basin are not well understood. The vulnerability of terminal lakes is related to the factors driving change and adaptive capacity of the lake ecosystem. Saline lake ecosystems are vulnerable when changes in water quantity affect ecosystem function. Water quantity affects salinity, which affects food webs and habitat; these linkages can be investigated with water-quality and food web monitoring. Water quantity also affects inundated habitat, which can be quantified through remote sensing. It is necessary to quantify hydroclimatic and water use controls on water availability to terminal lakes to assess the response of the ecosystems. Remotely sensed data can provide a broad-scale and long-term synoptic view of terminal lake hydrologic characteristics, but ground observations are required to interpret changes in water quality and ecological functions. Some terminal lake basins have ongoing monitoring and modeling efforts within the Great Basin (for example, Great Salt Lake, Carson River Basin), yet most monitoring locations are hydrologically upgradient and too far away from lake inflows to provide an accurate assessment of hydrological trends for the lake ecosystems. Other terminal lakes have no long-term hydrological monitoring in their respective watersheds (for example, Lake Abert). Ecological data collection in the Great Basin is also insufficient to understand how many birds exist on the landscape, how birds use the mosaic of terminal-lake habitats as an interconnected system, and how Great Basin terminal lakes are linked to the larger continental system of the Pacific and Central Flyways. Across agencies and organizations, tracking bird movement, abundance, and diversity is inconsistent, with some lakes having once- or twice-a-year bird survey efforts and a few locations having more intensive ecological data-gathering efforts (for example, Great Salt Lake, Lake Abert). Bridging hydrological and ecological information gaps will improve understanding of the trends in water supply and water quality, habitat availability and usage, and impacts on vulnerable waterbird species, all of which would be used by managers in coordinated conservation of this unique network of terminal-lake habitats. The terminal lakes of the Great Basin are part of the Basin and Range physiographic province that extends from the Colorado Plateau on the east to the Sierra Nevada on the west, and from the Snake River Plain on the north to the Garlock fault and the Mojave block on the south. The Great Basin is larger than 650,000 square kilometers and encompasses most of the State of Nevada but also extends to western Utah, eastern California, southeastern Idaho, southwestern Wyoming, and southeastern Oregon. The climate is arid to semiarid with a hydrologic regime that is snowmelt dominated, providing as much as 75 percent of total annual runoff for the region. Terminal lakes of the Great Basin occupy the lowest areas of closed (endorheic) drainage basins, such that lake levels and water quality respond rapidly to surface-water inflow. Terminal lakes provide local and regional economic value to the States in the Great Basin, including mineral extraction, aquaculture, public works, and recreational uses. As an example, assessments of Great Salt Lake’s ecological health and economic impact find hemispheric importance for the former and regional importance for the latter. Great Salt Lake creates about 7,000 jobs and $2 billion of economic output per year, most of which would be lost with further declines in lake level. The objectives of this Science Strategy are threefold: (1) to identify how changing water availability affects the quality, diversity, and abundance of habitats supporting continental waterbird populations; (2) to highlight the scientific monitoring and assessment needs of Great Basin terminal lakes; and (3) to support coordinated management and conservation actions to benefit those ecosystems, migratory birds, and other wildlife. There are long-term hydrological, ecological, and societal challenges associated with terminal lakes ecosystems in the Great Basin. This Science Strategy benefits partners by providing a conceptual model, nested at different spatial extents, that identifies key scientific information needs to inform coordinated implementation of management and conservation plans within and among hydrologic basins to address these complex challenges.

California, Idaho, Nevada, Oregon, Utah, Wyoming↗

Linking urbanization to the Biological Condition Gradient (BCG) for stream ecosystems in the Northeastern United States using a Bayesian network approach

Urban development alters important physical, chemical, and biological processes that define urban stream ecosystems. An approach was developed for quantifying the effects of these processes on aquatic biota, and then linking those effects to endpoints that can be used for environmental management. These complex, interacting systems are challenging to model from a scientific standpoint. A desirable model clearly shows the system, simulates the interactions, and ultimately predicts results of management actions. Traditional regression techniques that calculate empirical relations between pairs of environmental factors do not capture the interconnected web of multiple stressors, but urban development effects are not yet understood at the detailed scales required to make mechanistic modeling approaches feasible. Therefore, in contrast to a fully deterministic or fully statistical modeling approach, a Bayesian network model provides a hybrid approach that can be used to represent known general associations between variables while acknowledging uncertainty in predicted outcomes. It does so by quantifying an expert-elicited network of probabilistic relations between variables. Advantages of this modeling approach include (1) flexibility in accommodating many model specifications and information types; (2) efficiency in storing and manipulating complex information, and to parameterize; and (3) transparency in describing the relations using nodes and arrows and in describing uncertainties with discrete probability distributions for each variable. In realization of the aforementioned advantages, a Bayesian network model was constructed to characterize the effect of urban development on aquatic macroinvertebrate stream communities through three simultaneous, interacting ecological pathways affecting stream hydrology, habitat, and water quality across watersheds in the Northeastern United States. This model incorporates both empirical data and expert knowledge to calculate the probabilities of attaining desired aquatic ecosystem conditions under different urban stress levels, environmental conditions, and management options. Ecosystem conditions are characterized in terms of standardized Biological Condition Gradient (BCG) management endpoints. This approach to evaluating urban development-induced perturbations in watersheds integrates statistical and mechanistic perspectives, different information sources, and several ecological processes into a comprehensive description of the system that can be used to support decision making. The completed model can be used to infer which management actions would lead to the highest likelihood of desired BCG tier achievement. For example, if best management practices (BMP) were implemented in a highly urbanized watershed to reduce flashiness to medium levels and specific conductance to low levels, the stream would have a 70-percent chance of achieving BCG Tier 3 or better, relative to a 24-percent achievement likelihood for unmanaged high urban land cover. Results are reported probabilistically to account for modeling uncertainty that is inherent in sources such as natural variability and model simplification error.

Scientific Investigations Report↗