Observed and projected functional reorganization of riverine fish assemblages from global change.

Climate and land-use/land-cover change ("global change") are restructuring biodiversity, globally. Broadly, environmental conditions are expected to become warmer, potentially drier (particularly in arid regions), and more anthropogenically developed in the future, with spatiotemporally complex effects on ecological communities. We used functional traits to inform Chesapeake Bay Watershed fish responses to future climate and land-use scenarios (2030, 2060, and 2090). We modeled the future habitat suitability of focal species representative of key trait axes (substrate, flow, temperature, reproduction, and trophic) and used functional and phylogenetic metrics to assess variable assemblage responses across physiographic regions and habitat sizes (headwaters through large rivers). Our focal species analysis projected future habitat suitability gains for carnivorous species with preferences for warm water, pool habitats, and fine or vegetated substrates. At the assemblage level, models projected decreasing habitat suitability for cold-water, rheophilic, and lithophilic individuals but increasing suitability for carnivores in the future across all regions. Projected responses of functional and phylogenetic diversity and redundancy differed among regions. Lowland regions were projected to become less functionally and phylogenetically diverse and more redundant while upland regions (and smaller habitat sizes) were projected to become more diverse and less redundant. Next, we assessed how these model-projected assemblage changes 2005-2030 related to observed time-series trends (1999-2016). Halfway through the initial projecting period (2005-2030), we found observed trends broadly followed modeled patterns of increasing proportions of carnivorous and lithophilic individuals in lowland regions but showed opposing patterns for functional and phylogenetic metrics. Leveraging observed and predicted analyses simultaneously helps elucidate the instances and causes of discrepancies between model predictions and ongoing observed changes. Collectively, results highlight the complexity of global change impacts across broad landscapes that likely relate to differences in assemblages' intrinsic sensitivities and external exposure to stressors.

Explainable machine learning improves interpretability in the predictive modeling of biological stream conditions in the Chesapeake Bay Watershed, USA.

Anthropogenic alterations have resulted in widespread degradation of stream conditions. To aid in stream restoration and management, baseline estimates of conditions and improved explanation of factors driving their degradation are needed. We used random forests to model biological conditions using a benthic macroinvertebrate index of biotic integrity for small, non-tidal streams (upstream area ≤200 km2) in the Chesapeake Bay watershed (CBW) of the mid-Atlantic coast of North America. We utilized several global and local model interpretation tools to improve average and site-specific model inferences, respectively. The model was used to predict condition for 95,867 individual catchments for eight periods (2001, 2004, 2006, 2008, 2011, 2013, 2016, 2019). Predicted conditions were classified as Poor, FairGood, or Uncertain to align with management needs and individual reach lengths and catchment areas were summed by condition class for the CBW for each period. Global permutation and local Shapley importance values indicated percent of forest, development, and agriculture in upstream catchments had strong impacts on predictions. Development and agriculture negatively influenced stream condition for model average (partial dependence [PD] and accumulated local effect [ALE] plots) and local (individual condition expectation and Shapley value plots) levels. Friedman's H-statistic indicated large overall interactions for these three land covers, and bivariate global plots (PD and ALE) supported interactions among agriculture and development. Total stream length and catchment area predicted in FairGood conditions decreased then increased over the 19-years (length/area: 66.6/65.4% in 2001, 66.3/65.2% in 2011, and 66.6/65.4% in 2019). Examination of individual catchment predictions between 2001 and 2019 showed those predicted to have the largest decreases in condition had large increases in development; whereas catchments predicted to exhibit the largest increases in condition showed moderate increases in forest cover. Use of global and local interpretative methods together with watershed-wide and individual catchment predictions support conservation practitioners that need to identify widespread and localized patterns, especially acknowledging that management actions typically take place at individual-reach scales.

Time marches on, but do the causal pathways driving instream habitat and biology remain consistent?

Stream ecosystems are complex networks of interacting terrestrial and aquatic drivers. To untangle these ecological networks, efforts evaluating the direct and indirect effects of landscape, climate, and instream predictors on biological condition through time are needed. We used structural equation modeling and leveraged a stream survey program to identify and compare important predictors driving condition of benthic macroinvertebrate and fish assemblages. We used data resampled 14 years apart at 252 locations across Maryland, USA. Sample locations covered a wide range of conditions that varied spatiotemporally. Overall, the relationship directions were consistent between sample periods, but their relative strength varied temporally. For benthic macroinvertebrates, we found that the total effect of natural landscape (e.g., elevation, longitude, latitude, geology) and land use (i.e., forest, development, agriculture) predictors was 1.4 and 1.5 times greater in the late 2010s compared to the 2000s. Moreover, the total effect of water quality (e.g., total nitrogen and conductivity) and habitat (e.g., embeddedness, riffle quality) was 1.2 and 4.8 times lower in the 2010s, respectively. For fish assemblage condition, the total effect of land use-land cover predictors was 2.3 times greater in the 2010s compared to the 2000s, while the total effect of local habitat was 1.4 times lower in the 2010s, respectively. As expected, we found biological assemblages in catchments with more agriculture and urban development were generally comprised of tolerant, generalist species, while assemblages in catchments with greater forest cover had more-specialized, less-tolerant species (e.g., Ephemeroptera, Plecoptera, and Trichoptera taxa, clingers, benthic and lithophilic spawning fishes). Changes in the relative importance of landscape and land-use predictors suggest other correlated, yet unmeasured, proximal factors became more important over time. By untangling these ecological networks, stakeholders can gain a better understanding of the spatiotemporal relationships driving biological condition to implement management practices aimed at improving stream condition.

Disentangling the potential effects of land-use and climate change on stream conditions.

Land-use and climate change are significantly affecting stream ecosystems, yet understanding of their long-term impacts is hindered by the few studies that have simultaneously investigated their interaction and high variability among future projections. We modeled possible effects of a suite of 2030, 2060, and 2090 land-use and climate scenarios on the condition of 70,772 small streams in the Chesapeake Bay watershed, United States. The Chesapeake Basin-wide Index of Biotic Integrity, a benthic macroinvertebrate multimetric index, was used to represent stream condition. Land-use scenarios included four Special Report on Emissions Scenarios (A1B, A2, B1, and B2) representing a range of potential landscape futures. Future climate scenarios included quartiles of future climate changes from downscaled Coupled Model Intercomparison Project - Phase 5 (CMIP5) and a watershed-wide uniform scenario (Lynch2016). We employed random forests analysis to model individual and combined effects of land-use and climate change on stream conditions. Individual scenarios suggest that by 2090, watershed-wide conditions may exhibit anywhere from large degradations (e.g., scenarios A1B, A2, and the CMIP5 25th percentile) to small degradations (e.g., scenarios B1, B2, and Lynch2016). Combined land-use and climate change scenarios highlighted their interaction and predicted, by 2090, watershed-wide degradation in 16.2% (A2 CMIP5 25th percentile) to 1.0% (B2 Lynch2016) of stream kilometers. A goal for the Chesapeake Bay watershed is to restore 10% of stream kilometers over a 2008 baseline; our results suggest meeting and sustaining this goal until 2090 may require improvement in 11.0%-26.2% of stream kilometers, dependent on land-use and climate scenario. These results highlight inherent variability among scenarios and the resultant uncertainty of predicted conditions, which reinforces the need to incorporate multiple scenarios of both land-use (e.g., development, agriculture, etc.) and climate change in future studies to encapsulate the range of potential future conditions.