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.

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.

A terrain-based paired-site sampling design to assess biodiversity losses from eastern hemlock decline.

Biodiversity surveys are often hampered by the inability to control extraneous sources of variability introduced into comparisons of populations across a heterogenous landscape. If not specifically accounted for a priori, this noise can weaken comparisons between sites, and can make it difficult to draw inferences about specific ecological processes. We developed a terrain-based, paired-site sampling design to analyze differences in aquatic biodiversity between streams draining eastern hemlock (Tsuga canadensis) forests, and those draining mixed hardwood forests in Delaware Water Gap National Recreation Area (USA). The goal of this design was to minimize variance due to terrain influences on stream communities, while representing the range of hemlock dominated stream environments present in the park. We used geographic information systems (GIS) and cluster analysis to define and partition hemlock dominated streams into terrain types based on topographic variables and stream order. We computed similarity of forest stands within terrain types and used this information to pair hemlock-dominated streams with hardwood counterparts prior to sampling. We evaluated the effectiveness of the design through power analysis and found that power to detect differences in aquatic invertebrate taxa richness was highest when sites were paired and terrain type was included as a factor in the analysis. Precision of the estimated difference in mean richness was nearly doubled using the terrain-based, paired site design in comparison to other evaluated designs. Use of this method allowed us to sample stream communities representative of park-wide forest conditions while effectively controlling for landscape variability.