Spatial and Temporal Variability in Stream Thermal Regime Drivers for Three River Networks During the Summer Growing Season.

Many cold-water dependent aquatic organisms are experiencing habitat and population declines from increasing water temperatures. Identifying mechanisms which drive local and regional stream thermal regimes facilitates restoration at ecologically relevant scales. Stream temperatures vary spatially and temporally both within and among river basins. We developed a modeling process to identify statistical relationships between drivers of stream temperature and covariates representing landscape, climate, and management-related processes. The modeling process was tested in 3 study areas of the Pacific Northwest USA during the growing season (May [start], August [warmest], September [end]). Across all months and study systems, covariates with the highest relative importance represented the physical landscape (elevation [1st], catchment area [3rd], main channel slope [5th]) and climate covariates (mean monthly air temperature [2nd] and discharge [4th]). Two management covariates (ground water use [6th] and riparian shade [7th]) also had high relative importance. Across the growing season (for all basins) local reach slope had high relative importance in May, but transitioned to a regional main channel slope covariate in August and September. This modeling process identified regionally similar and locally unique relationships among drivers of stream temperature. High relative importance of management-related covariates suggested potential restoration actions for each system.

Riparian vegetation shade restoration and loss effects on recent and future stream temperatures.

River temperatures are expected to increase this century harming species requiring cold-water habitat unless restoration activities protect or improve habitat availability. Local shading by riparian vegetation can cool water temperatures, but uncertainty exists over the scaling of this local effect to larger spatial extents. We evaluate this issue using a regional spatial stream network temperature model with covariates representing shade effects to predict mean August stream temperatures across 78,195 km of tributaries flowing into the Columbia River in the northwestern US. We evaluate nine scenarios predicting stream temperatures for three riparian shade conditions (current, restored, and no riparian vegetation) within three different climate periods (2000s, 2040s, and 2080s). Results suggest riparian shade restoration (2000s climate) could decrease mean August stream temperatures by 0.62°C across the study network. Under the same restored shade conditions, temperature predictions for tributaries at their confluence with the Columbia River range from 0.02-2.08°C cooler than under current shade conditions. The climate warming effect predicted for the 2040s and 2080s, however, is greater than the cooling effect from restoring riparian shade. Streams less than 10m bankfull width cooled more frequently with riparian shade restoration. In Oregon, the proportion of fish habitat for salmon and trout rearing and migration that meet temperature numeric water quality criteria could be increased by 20% under restored shade conditions although net habitat declines may still occur in the future. We conclude riparian vegetation restoration could partially mitigate future warming and help maintain cold-water habitats that function as thermal refuges if implemented strategically.

Integrating thermal infrared stream temperature imagery and spatial stream network models to understand natural spatial thermal variability in streams.

Under a warmer future climate, thermal refuges could facilitate the persistence of species relying on cold-water habitat. Often these refuges are small and easily missed or smoothed out by averaging in models. Thermal infrared (TIR) imagery can provide empirical water surface temperatures that capture these features at a high spatial resolution (<1 m) and over tens of kilometers. Our study examined how TIR data could be used along with spatial stream network (SSN) models to characterize thermal regimes spatially in the Middle Fork John Day (MFJD) River mainstem (Oregon, USA). We characterized thermal variation in seven TIR longitudinal temperature profiles along the MFJD mainstem and compared them with SSN model predictions of stream temperature (for the same time periods as the TIR profiles). TIR profiles identified reaches of the MFJD mainstem with consistently cooler temperatures across years that were not consistently captured by the SSN prediction models. SSN predictions along the mainstem identified ~80% of the 1-km reach scale temperature warming or cooling trends observed in the TIR profiles. We assessed whether landscape features (e.g., tributary junctions, valley confinement, geomorphic reach classifications) could explain the fine-scale thermal heterogeneity in the TIR profiles (after accounting for the reach-scale temperature variability predicted by the SSN model) by fitting SSN models using the TIR profile observation points. Only the distance to the nearest upstream tributary was identified as a statistically significant landscape feature for explaining some of the thermal variability in the TIR profile data. When combined, TIR data and SSN models provide a data-rich evaluation of stream temperature captured in TIR imagery and a spatially extensive prediction of the network thermal diversity from the outlet to the headwaters.

Seagrass habitat metabolism increases short-term extremes and long-term offset of CO2 under future ocean acidification.

The role of rising atmospheric CO2 in modulating estuarine carbonate system dynamics remains poorly characterized, likely due to myriad processes driving the complex chemistry in these habitats. We reconstructed the full carbonate system of an estuarine seagrass habitat for a summer period of 2.5 months utilizing a combination of time-series observations and mechanistic modeling, and quantified the roles of aerobic metabolism, mixing, and gas exchange in the observed dynamics. The anthropogenic CO2 burden in the habitat was estimated for the years 1765-2100 to quantify changes in observed high-frequency carbonate chemistry dynamics. The addition of anthropogenic CO2 alters the thermodynamic buffer factors (e.g., the Revelle factor) of the carbonate system, decreasing the seagrass habitat's ability to buffer natural carbonate system fluctuations. As a result, the most harmful carbonate system indices for many estuarine organisms [minimum pHT, minimum Ωarag, and maximum pCO2(s.w.)] change up to 1.8×, 2.3×, and 1.5× more rapidly than the medians for each parameter, respectively. In this system, the relative benefits of the seagrass habitat in locally mitigating ocean acidification increase with the higher atmospheric CO2 levels predicted toward 2100. Presently, however, these mitigating effects are mixed due to intense diel cycling of CO2 driven by aerobic metabolism. This study provides estimates of how high-frequency pHT, Ωarag, and pCO2(s.w.) dynamics are altered by rising atmospheric CO2 in an estuarine habitat, and highlights nonlinear responses of coastal carbonate parameters to ocean acidification relevant for water quality management.

Toxicity of atmospheric aerosols on marine phytoplankton.

Atmospheric aerosol deposition is an important source of nutrients and trace metals to the open ocean that can enhance ocean productivity and carbon sequestration and thus influence atmospheric carbon dioxide concentrations and climate. Using aerosol samples from different back trajectories in incubation experiments with natural communities, we demonstrate that the response of phytoplankton growth to aerosol additions depends on specific components in aerosols and differs across phytoplankton species. Aerosol additions enhanced growth by releasing nitrogen and phosphorus, but not all aerosols stimulated growth. Toxic effects were observed with some aerosols, where the toxicity affected picoeukaryotes and Synechococcus but not Prochlorococcus. We suggest that the toxicity could be due to high copper concentrations in these aerosols and support this by laboratory copper toxicity tests preformed with Synechococcus cultures. However, it is possible that other elements present in the aerosols or unknown synergistic effects between these elements could have also contributed to the toxic effect. Anthropogenic emissions are increasing atmospheric copper deposition sharply, and based on coupled atmosphere-ocean calculations, we show that this deposition can potentially alter patterns of marine primary production and community structure in high aerosol, low chlorophyll areas, particularly in the Bay of Bengal and downwind of South and East Asia.