Spatial asynchrony and cross-scale climate interactions in populations of a coldwater stream fish

Climate change affects populations over broad geographic ranges due to spatially autocorrelated abiotic conditions known as the Moran effect. However, populations do not always respond to broad-scale environmental changes synchronously across a landscape. We combined multiple datasets for a retrospective analysis of time-series count data (5-28 annual samples per segment) at 144 stream segments dispersed over nearly 1,000 linear kilometers of range to characterize the population structure and scale of spatial synchrony across the southern native range of a coldwater stream fish (brook trout, Salvelinus fontinalis), which is sensitive to stream temperature and flow variations. Spatial synchrony differed by life stage and geographic region: it was stronger in the juvenile life stage than in the adult life stage and in the northern sub-region than in the southern sub-region. Spatial synchrony of trout populations extended to 100-200 km but was much weaker than that of climate variables such as temperature, precipitation, and stream flow. Early life stage abundance changed over time due to annual variation in summer temperature and winter and spring stream flow conditions. Climate effects on abundance differed between sub-regions and among local populations within sub-regions, indicating multiple cross-scale interactions where climate interacted with local habitat to generate only a modest pattern of population synchrony over space. Overall, our analysis showed higher degrees of response heterogeneity of local populations to climate variation and consequently population asynchrony than previously shown based on analysis of individual, geographically restricted datasets. This response heterogeneity indicates that certain local segments characterized by population asynchrony and resistance to climate variation could represent unique populations of this iconic native coldwater fish that warrant targeted conservation. Advancing the conservation of this species can include actions that identify such priority populations and incorporate them into landscape-level conservation planning. Our approach is applicable to other widespread aquatic species sensitive to climate change.

Life history strategies of stream fishes linked to predictors of hydrologic stability

Life history theory provides a framework to understand environmental change based on species strategies for survival and reproduction under stable, cyclical, or stochastic environmental conditions. We evaluated environmental predictors of fish life history strategies in 20 streams intersecting a national park within the Potomac River basin in eastern North America. We sampled stream sites during 2018–2019 and collected 3801 individuals representing 51 species within 10 taxonomic families. We quantified life history strategies for species from their coordinates in an ordination space defined by trade‐offs in spawning season duration, fecundity, and parental care characteristic of opportunistic, periodic, and equilibrium strategies. Our analysis revealed important environmental predictors: Abundance of opportunistic strategists increased with low‐permeability soils that produce flashy runoff dynamics and decreased with karst terrain (carbonate bedrock) where groundwater inputs stabilize stream flow and temperature. Conversely, abundance of equilibrium strategists increased in karst terrain indicating a response to more stable environmental conditions. Our study indicated that fish community responses to groundwater and runoff processes may be explained by species traits for survival and reproduction. Our findings also suggest the utility of life history theory for understanding ecological responses to destabilized environmental conditions under global climate change.

Individual behaviour and resource use of thermally stressed brook trout Salvelinus fontinalis portend the conservation potential of thermal refugia

Individual aggression and thermal refuge use were monitored in brook trout Salvelinus fontinalis in a controlled laboratory to determine how fish size and personality influence time spent in forage and thermal habitat patches during periods of thermal stress. On average, larger and more exploratory fish initiated more aggressive interactions and across all fish there was decreased aggression at warmer temperatures. Individual personality did not explain changes in aggression or habitat use with increased temperature; however, larger individuals initiated comparatively fewer aggressive interactions at warmer temperatures. Occupancy of forage patches generally declined as ambient stream temperatures approached critical maximum and fish increased thermal refuge use, with a steeper decline in forage patch occupancy observed in larger fish. These findings suggest that larger individuals may be more vulnerable to stream temperature rise. Importantly, even at thermally stressful temperatures, all fish periodically left the thermal refuge to forage. This indicates that the success of refugia at increasing population survival during periods of stream temperature rise may depend on the location of thermal refugia relative to forage locations within the larger habitat mosaic. These results provide insights into the potential for thermal refugia to improve population survival and can be used to inform predictions of population vulnerability to climate change.

Inferring watershed hydraulics and cold-water habitat persistence using multi-year air and stream temperature signals

Streams strongly influenced by groundwater discharge may serve as "climate refugia" for sensitive species in regions of increasingly marginal thermal conditions. The main goal of this study is to develop paired air and stream water annual temperature signal analysis techniques to elucidate the relative groundwater contribution to stream water and the effective groundwater flowpath depth. Groundwater discharge to streams attenuates surface water temperature signals, and this attenuation can be diagnostic of groundwater gaining systems. Additionally, discharge from shallow groundwater flowpaths can theoretically transfer lagged annual temperature signals from aquifer to stream water. Here we explore this concept using multi-year temperature records from 120 stream sites located across 18 mountain watersheds of Shenandoah National Park, VA, USA and a coastal watershed in Massachusetts, USA. Both areas constitute important cold-water habitat for native brook trout (Salvelinus fontinalis). Observed annual temperature signals indicate a dominance of shallow groundwater discharge to streams in the National Park, in contrast to the coastal watershed that has strong, apparently deeper, groundwater influence. The average phase lag from air to stream signals in Shenandoah National Park is 11 d; however, extended lags of approximately 1 month were observed in a subset of streams. In contrast, the coastal stream has pronounced attenuation of annual temperature signals without notable phase lag. To better understand these observed differences in signal characteristics, analytical and numerical models are used to quantify mixing of the annual temperature signals of surface and groundwater. Simulations using a total heat budget numerical model indicate groundwater-induced annual temperature signal phase lags are likely to show greater downstream propagation than the related signal amplitude attenuation. The measurement of multi-seasonal paired air and water temperatures offers great promise toward understanding catchment processes and informing current cold-water habitat management at ecologically-relevant scales.

Dermocystidium sp. Infection in Blue Ridge Sculpin Captured in Maryland

Raised pale cysts were observed on Blue Ridge Sculpin Cottus caeruleomentum during stream fish community surveys in Catoctin Mountain Park, Maryland. When examined histologically, preserved sculpin exhibited multiple cysts containing spherical endospores with a refractile central body characteristic of Dermocystidium spp. Cysts were not observed on the gills or internally. The portion of the watershed in which affected sculpin were observed contained lower than expected numbers of sculpin, raising concerns about the population effects of this infection. A nearby stream lacked sculpin even though they are common in this region, further suggesting the possibility of regional effects. This is the first report of a Dermocystidium infecting any fish species in the eastern United States. Received October 16, 2015; accepted February 14, 2016.

Accounting for groundwater in stream fish thermal habitat responses to climate change

Forecasting climate change effects on aquatic fauna and their habitat requires an understanding of how water temperature responds to changing air temperature (i.e., thermal sensitivity). Previous efforts to forecast climate effects on brook trout (Salvelinus fontinalis) habitat have generally assumed uniform air-water temperature relationships over large areas that cannot account for groundwater inputs and other processes that operate at finer spatial scales. We developed regression models that accounted for groundwater influences on thermal sensitivity from measured air-water temperature relationships within forested watersheds in eastern North America (Shenandoah National Park, Virginia, USA, 78 sites in nine watersheds). We used these reach-scale models to forecast climate change effects on stream temperature and brook trout thermal habitat, and compared our results to previous forecasts based upon large-scale models. Observed stream temperatures were generally less sensitive to air temperature than previously assumed, and we attribute this to the moderating effect of shallow groundwater inputs. Predicted groundwater temperatures from air-water regression models corresponded well to observed groundwater temperatures elsewhere in the study area. Predictions of brook trout future habitat loss derived from our fine-grained models. were far less pessimistic than those from prior models developed at coarser spatial resolutions. However, our models also revealed spatial variation in thermal sensitivity within and among catchments resulting in a patchy distribution of thermally suitable habitat. Habitat fragmentation due to thermal barriers therefore may have an increasingly important role for trout population viability in headwater streams. Our results demonstrate that simple adjustments to air-water temperature regression models can provide a powerful and cost-effective approach for predicting future stream temperatures while accounting for effects of groundwater.

Seasonal weather patterns drive population vital rates and persistence in a stream fish

Climate change affects seasonal weather patterns, but little is known about the relative importance of seasonal weather patterns on animal population vital rates. Even when such information exists, data are typically only available from intensive fieldwork (e.g., mark-recapture studies) at a limited spatial extent. Here, we investigated effects of seasonal air temperature and precipitation (fall, winter, and spring) on survival and recruitment of brook trout (Salvelinus fontinalis) at a broad spatial scale using a novel stage-structured population model. The data were a 15-year record of brook trout abundance from 72 sites distributed across a 170-km-long mountain range in Shenandoah National Park, Virginia, USA. Population vital rates responded differently to weather and site-specific conditions. Specifically, young-of-year survival was most strongly affected by spring temperature, adult survival by elevation and per-capita recruitment by winter precipitation. Low fall precipitation and high winter precipitation, the latter of which is predicted to increase under climate change for the study region, had the strongest negative effects on trout populations. Simulations show that trout abundance could be greatly reduced under constant high winter precipitation, consistent with the expected effects of gravel-scouring flows on eggs and newly hatched individuals. However, high-elevation sites would be less vulnerable to local extinction because they supported higher adult survival. Furthermore, the majority of brook trout populations are projected to persist if high winter precipitation occurs only intermittently (≤3 of 5 years) due to density-dependent recruitment. Variable drivers of vital rates should be commonly found in animal populations characterized by ontogenetic changes in habitat, and such stage-structured effects may increase population persistence to changing climate by not affecting all life stages simultaneously. Yet, our results also demonstrate that weather patterns during seemingly less consequential seasons (e.g., winter precipitation) can have major impacts on animal population dynamics.

Threshold-dependent sample sizes for selenium assessment with stream fish tissue

Natural resource managers are developing assessments of selenium (Se) contamination in freshwater ecosystems based on fish tissue concentrations. We evaluated the effects of sample size (i.e., number of fish per site) on the probability of correctly detecting mean whole-body Se values above a range of potential management thresholds. We modeled Se concentrations as gamma distributions with shape and scale parameters fitting an empirical mean-to-variance relationship in data from southwestern West Virginia, USA (63 collections, 382 individuals). We used parametric bootstrapping techniques to calculate statistical power as the probability of detecting true mean concentrations up to 3 mg Se/kg above management thresholds ranging from 4 to 8 mg Se/kg. Sample sizes required to achieve 80% power varied as a function of management thresholds and Type I error tolerance (α). Higher thresholds required more samples than lower thresholds because populations were more heterogeneous at higher mean Se levels. For instance, to assess a management threshold of 4 mg Se/kg, a sample of eight fish could detect an increase of approximately 1 mg Se/kg with 80% power (given α=0.05), but this sample size would be unable to detect such an increase from a management threshold of 8 mg Se/kg with more than a coin-flip probability. Increasing α decreased sample size requirements to detect above-threshold mean Se concentrations with 80% power. For instance, at an α-level of 0.05, an 8-fish sample could detect an increase of approximately 2 units above a threshold of 8 mg Se/kg with 80% power, but when α was relaxed to 0.2, this sample size was more sensitive to increasing mean Se concentrations, allowing detection of an increase of approximately 1.2 units with equivalent power. Combining individuals into 2- and 4-fish composite samples for laboratory analysis did not decrease power because the reduced number of laboratory samples was compensated for by increased precision of composites for estimating mean conditions. However, low sample sizes (<5 fish) did not achieve 80% power to detect near-threshold values (i.e., <1 mg Se/kg) under any scenario we evaluated. This analysis can assist the sampling design and interpretation of Se assessments from fish tissue by accounting for natural variation in stream fish populations.

Experimental and environmental factors affect spurious detection of ecological thresholds

Threshold detection methods are increasingly popular for assessing nonlinear responses to environmental change, but their statistical performance remains poorly understood. We simulated linear change in stream benthic macroinvertebrate communities and evaluated the performance of commonly used threshold detection methods based on model fitting (piecewise quantile regression [PQR]), data partitioning (nonparametric change point analysis [NCPA]), and a hybrid approach (significant zero crossings [SiZer]). We demonstrated that false detection of ecological thresholds (type I errors) and inferences on threshold locations are influenced by sample size, rate of linear change, and frequency of observations across the environmental gradient (i.e., sample-environment distribution, SED). However, the relative importance of these factors varied among statistical methods and between inference types. False detection rates were influenced primarily by user-selected parameters for PQR (tau) and SiZer (bandwidth) and secondarily by sample size (for PQR) and SED (for SiZer). In contrast, the location of reported thresholds was influenced primarily by SED. Bootstrapped confidence intervals for NCPA threshold locations revealed strong correspondence to SED. We conclude that the choice of statistical methods for threshold detection should be matched to experimental and environmental constraints to minimize false detection rates and avoid spurious inferences regarding threshold location.

Histologic, immunologic and endocrine biomarkers indicate contaminant effects in fishes of the Ashtabula River

The use of fish as sentinels of aquatic ecosystem health is a biologically relevant approach to environmental monitoring and assessment. We examined the health of the Ashtabula River using histologic, immunologic, and endocrine biomarkers in brown bullhead (BB; Ameiurus nebulosus) and largemouth bass (Micropterus salmoides) and compared fish collected from a reference site (Conneaut Creek). Seasonal analysis was necessary to distinguish differences in fish between the two rivers. Overall BB from the Ashtabula River had a lower condition factor and significantly more macrophage aggregates than those from the reference site. Reduced bactericidal and cytotoxic-cell activity was observed in anterior kidney leukocytes from both BB and largemouth bass from the Ashtabula River. Lower plasma thyroxine and triiodo-L-thyronine in both species in the Ashtabula River indicated disruption of the thyroid axis. Differences in physiological biomarker responses were supported by body burden chemical concentrations when data were analyzed on a seasonal basis. The use of two fish species added a level of rigor that demonstrated biological effects were not exclusive to a single species. The results provide strong evidence that contaminants have affected fish in the Ashtabula River, a Great Lakes Area of Concern, and provide a baseline by which to evaluate remediation activities.

Ecological integrity of streams related to human cancer mortality rates

Assessments of ecological integrity have become commonplace for biological conservation, but their role for public health analysis remains largely unexplored. We tested the prediction that the ecological integrity of streams would provide an indicator of human cancer mortality rates in West Virginia, USA. We characterized ecological integrity using an index of benthic macroinvertebrate community structure (West Virginia Stream Condition Index, SCI) and quantified human cancer mortality rates using county-level data from the Centers for Disease Control and Prevention. Regression and spatial analyses revealed significant associations between ecological integrity and public health. SCI was negatively related to age-adjusted total cancer mortality per 100,000 people. Respiratory, digestive, urinary, and breast cancer rates increased with ecological disintegrity, but genital and oral cancer rates did not. Smoking, poverty, and urbanization were significantly related to total cancer mortality, but did not explain the observed relationships between ecological integrity and cancer. Coal mining was significantly associated with ecological disintegrity and higher cancer mortality. Spatial analyses also revealed cancer clusters that corresponded to areas of high coal mining intensity. Our results demonstrated significant relationships between ecological integrity and human cancer mortality in West Virginia, and suggested important effects of coal mining on ecological communities and public health. Assessments of ecological integrity therefore may contribute not only to monitoring goals for aquatic life, but also may provide valuable insights for human health and safety.

River-stream connectivity affects fish bioassessment performance

Stream fish bioassessment methods assume that fish assemblages observed in sample sites reflect responses to local stressors, but fish assemblages are influenced by local factors as well as regional dispersal to and from connected streams. We hypothesized that fish movement to and from refugia and source populations in connected rivers (i.e., riverine dispersal) would weaken or decouple relations between fish community metrics and local environmental conditions. We compared fish-environment relations between streams that flow into large rivers (mainstem tributaries) and streams that lack riverine confluences (headwater tributaries) at multiple spatial grains using data from the USEPA's Environmental Monitoring and Assessment Program in the mid-Atlantic highlands, USA (n = 157 sites). Headwater and mainstem tributaries were not different in local environmental conditions, but showed important differences in fish metric responses to environmental quality gradients. Stream sites flowing into mainstem channels within 10 fluvial km showed consistently weaker relations to local environmental conditions than stream sites that lacked such mainstem connections. Moreover, these patterns diminished at longer distances from riverine confluences, consistent with the hypothesis of riverine dispersal. Our results suggest that (1) the precision of fish bioassessment metrics may be improved by calibrating scoring criteria based on the spatial position of sites within stream networks and (2) the spatial grain of fish bioassessment studies may be manipulated to suit objectives by including or excluding fishes exhibiting riverine dispersal.