Estimation of daily stream water temperatures with a Bayesian regression approach

Stream water temperature plays a significant role in aquatic ecosystems where it controls many important biological and physical processes. Reliable estimates of water temperature at the daily time step are critical in managing water resources. We developed a parsimonious piecewise Bayesian model for estimating daily stream water temperatures that account for temporal autocorrelation and both linear and nonlinear relationships with air temperature and discharge. The model was tested at 8 climatically different basins of the USA and at 34 sites within the mountainous Boise River Basin (Idaho, USA). The results show that the proposed model is robust with an average root mean square error of 1.25 °C and Nash–Sutcliffe coefficient of 0.92 over a 2‐year period. Our approach can be used to predict historic daily stream water temperatures in any location using observed daily stream temperature and regional air temperature data.     

Do southern Appalachian Mountain summer stream temperatures respond to removal of understory rhododendron thickets?

Dense understory thickets of the native evergreen shrub Rhododendron maximum expanded initially following elimination of American chestnut by the chestnut blight, and later in response to loss of the eastern hemlock due to hemlock woolly adelgid invasion. Rhododendron thickets often blanket streams and their riparian zones, creating cool, low-light microclimates. To determine the effect of such understory thickets on summer stream temperatures, we removed riparian rhododendron understory on 300 m reaches of two southern Appalachian Mountain headwater streams, while leaving two 300 m reference reaches undisturbed. Overhead canopy was left intact in all four streams, but all streams were selected to have a significant component of dead or dying eastern hemlock in the overstory, creating time-varying canopy gaps throughout the reach. We continuously monitored temperatures upstream, within and downstream of treatment and reference reaches. Temperatures were monitored in all four streams in the summer before treatments were imposed (2014), and for two summers following treatment (2015, 2016). Temperatures varied significantly across and within streams prior to treatment and across years for the reference streams. After rhododendron removal, increases in summer stream temperatures were observed at some locations within the treatment reaches, but these increases did not persist downstream and varied by watershed, sensor, and year. Significant increases in daily maxima in treatment reaches ranged from 0.9 to 2.6°C. Overhead canopy provided enough shade to prevent rhododendron removal from increasing summer temperatures to levels deleterious to native cold-water fauna (average summer temperatures remained below 16°C), and local temperature effects were not persistent.     

Prediction of river water temperature: a comparison between a new family of hybrid models and statistical approaches

River water temperature is a key physical variable controlling several chemical, biological and ecological processes. Its reliable prediction is a main issue in many environmental applications, which however is hampered by data scarcity, when using data‐demanding deterministic models, and modelling limitations, when using simpler statistical models. In this work we test a suite of models belonging to air2stream family, which are characterized by a hybrid formulation that combines a physical derivation of the key equation with a stochastic calibration of parameters. The air2stream models rely solely on air temperature and streamflow, and are of similar complexity as standard statistical models. The performances of the different versions of air2stream in predicting river water temperature are compared with those of the most common statistical models typically used in the literature. To this aim, a dataset of 38 Swiss rivers is used, which includes rivers classified into four different categories according to their hydrological characteristics: low‐land natural rivers, lake outlets, snow‐fed rivers and regulated rivers. The results of the analysis provide practical indications regarding the type of model that is most suitable to simulate river water temperature across different time scales (from daily to seasonal) and for different hydrological regimes. A model intercomparison exercise suggests that the family of air2stream hybrid models generally outperforms statistical models, while cross‐validation conducted over a 30‐year period indicates that they can be suitably adopted for long‐term analyses. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing the effects of catchment‐scale urban green infrastructure retrofits on hydrograph characteristics

Run‐off from impervious surfaces has pervasive and serious consequences for urban streams, but the detrimental effects of urban stormwater can be lessened by disconnecting impervious surfaces and redirecting run‐off to decentralized green infrastructure. This study used a before–after‐control‐impact design, in which streets served as subcatchments, to quantify hydrologic effectiveness of street‐scale investments in green infrastructure, such as street‐connected bioretention cells, rain gardens and rain barrels. On the two residential treatment streets, voluntary participation resulted in 32.2% and 13.5% of parcels having green infrastructure installed over a 2‐year period. Storm sewer discharge was measured before and after green infrastructure implementation, and peak discharge, total run‐off volume and hydrograph lags were analysed. On the street with smaller lots and lower participation, green infrastructure installation succeeded in reducing peak discharge by up to 33% and total storm run‐off by up to 40%. On the street with larger lots and higher participation, there was no significant reduction in peak or total stormflows, but on this street, contemporaneous street repairs may have offset improvements. On the street with smaller lots, lag times increased following the first phase of green infrastructure construction, in which streetside bioretention cells were built with underdrains. In the second phase, lag times did not change further, because bioretention cells were built without underdrains and water was removed from the system, rather than just delayed. We conclude that voluntary green infrastructure retrofits that include treatment of street run‐off can be effective for substantially reducing stormwater but that small differences in design and construction can be important for determining the level of the benefit. Copyright © 2015 John Wiley & Sons, Ltd.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

The response in floodplain respiration of an alpine river to experimental inundation under different temperature regimes

The respiratory potential [i.e. electron transport system activity (ETSA)] of soils and sediments from five floodplain habitats (channel, gravel, islands, riparian forest and grassland) of the Urbach River, Switzerland, and actual respiration rate (R) of the same samples exposed to experimental inundation were measured. Measurements were carried out at three incubation temperatures (4°C, 12°C and 20°C), and ETSA/R ratios (i.e. exploitation of the overall metabolic capacity) were investigated to better understand the effects of temperature and inundation on floodplain functional heterogeneity. Furthermore, ETSA/R ratios obtained during experimental inundation were compared with ETSA/R ratios from field measurements to investigate the exploitation in total metabolic potential at different conditions. Lowest ETSA and R were measured in samples from channel and gravel habitats, followed by those from islands. Substantially higher values were measured in soils from riparian forest and grassland. Both ETSA and R increased with increasing temperature in samples from all habitats, while the ETSA/R ratio decreased because of a rapid response in microbial community respiration to higher temperatures. The metabolic capacity exploitation (i.e. ETSA/R) during experimental inundation was lowest in predominantly terrestrial samples (riparian forest and grassland), indicating the weakest response to wetted conditions. Comparison of experimentally inundated and field conditions revealed that in rarely flooded soils, the metabolic capacity was less exploited during inundation than during non‐flooded conditions. The results suggest high sensitivity in floodplain respiration to changes in temperature and hydrological regime. ETSA/R ratios are considered good indicators of changes in metabolic activity of floodplain soils and sediments, and thus useful to estimate the impact of changes in hydrological regime or to evaluate success of floodplain restoration actions. Copyright © 2015 John Wiley & Sons, Ltd.     

Trends and multi‐annual variability of water temperatures in the river Danube,Serbia

The relationship between air (Ta) and water temperature (Tw) is very important because it shows how the temperature of a water body might respond to future changes in surface Ta. Mean monthly Tw records of three gauging stations (Bezdan, Bogojevo i Veliko Gradi?te) were analysed alongside mean monthly discharge (Q) for the same stations. Additionally, Ta series from two meteorological stations (Sombor and Veliko Gradi?te) were correlated with Tw variations over the period 1950–2012. Locally weighted scatter point smoothing (LOWESS) was used to investigate long‐term trends in the raw data, alongside the Mann–Kendall (MK) trend test. Trend significance was established using Yue–Pilon's pre‐whitening approaches to determine trends in climate data. Also, the rescaled adjusted partial sums (RAPS) method was used to detect dates of possible changes in the time series. Statistically significant warming trends were observed for annual and seasonal minimum and maximum Tw at all investigated sites. The strongest warming was observed at Bogojevo gauging station for seasonal maximum Tw, with +0.05 °C per year on average. RAPS established that the trend began in the 1980s. This behaviour is linked to climate patterns in the North and East Atlantic which determine the amount of heat advected onto mainland Europe. Statistically significant correlations were found for all Tw on an annual basis. Overall, the strongest correlations (p < 0.01) between Tw residuals and the North Atlantic Oscillation (NAO) were recorded for the winter period. These findings suggest possible predictability of Tw over seasonal time‐scales. Copyright © 2016 John Wiley & Sons, Ltd.     

The tectono‐metamorphic evolution of the very low‐grade hangingwall constrains two‐stage gneiss dome formation in the Montagne Noire (Southern France)

The Montagne Noire in the southernmost French Massif Central is made of an ENE‐elongated gneiss dome flanked by Palaeozoic sedimentary rocks. The tectonic evolution of the gneiss dome has generated controversy for more than half a century. As a result, a multitude of models have been proposed that invoke various tectonic regimes and exhumation mechanisms. Most of these models are based on data from the gneiss dome itself. Here, new constraints on the dome evolution are provided based on a combination of very low‐grade petrology, K–Ar geochronology, field mapping and structural analysis of the Palaeozoic western Mont Peyroux and Faugères units, which constitute part of the southern hangingwall of the dome. It is shown that southward‐directed Variscan nappe‐thrusting (D₁) and a related medium‐P metamorphism (M₁) are only preserved in the area furthest away from the gneiss dome. The regionally dominant pervasive tectono‐metamorphic event D₂/M₂ largely transposes D₁ structures, comprises a higher metamorphic thermal gradient than M₁ (transition low‐P and medium‐P metamorphic facies series) and affected the rocks between c. 309 and 300 Ma, post‐dating D₁/M₁ by more than 20 Ma. D₂‐related fabrics are refolded by D₃, which in its turn, is followed by dextral‐normal shearing along the basal shear zone of both units at c. 297 Ma. In the western Mont Peyroux and Faugères units, D₂/M₂ is largely synchronous with shearing along the southern dome margin between c. 311 and 303 Ma, facilitating the emplacement of the gneiss dome into the upper crust. D₂/M₂ also overlaps in time with granitic magmatism and migmatization in the Zone Axiale between c. 314 and 306 Ma, and a related low‐P/high‐T metamorphism at c. 308 Ma. The shearing that accompanied the exhumation of the dome therefore was synchronous with a peak in temperature expressed by migmatization and intrusion of melts within the dome, and also with the peak of metamorphism in the hangingwall. Both, the intensity of D₂ fabrics and the M₂ metamorphic grade within the hangingwall, decrease away from the gneiss dome, with grades ranging from the anchizone–epizone boundary to the diagenetic zone. The related zonation of the pre‐D₃ metamorphic field gradients paralleled the dome. These observations indicate that D₂/M₂ is controlled by the exhumation of the Zone Axiale, and suggest a coherent kinematic between the different crustal levels at some time during D₂/M₂. Based on integration of these findings with regional geological constraints, a two‐stage exhumation of the gneiss dome is proposed: during a first stage between c. 316 and 300 Ma dome emplacement into the upper crust was controlled by dextral shear zones arranged in a pull‐apart‐like geometry. The second stage from 300 Ma onwards was characterized by northeast to northward extension, with exhumation accommodated by north‐dipping detachments and hangingwall basin formation along the northeastern dome margin.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

River evaporation,condensation and heat fluxes within a first‐order tributary of Catamaran Brook (New Brunswick,Canada)

Stream temperature plays an important role in many biotic and abiotic processes, as it influences many physical, chemical and biological properties in rivers. As such, a good understanding of the thermal regime of rivers is essential for effective fisheries management and the protection aquatic habitats. Moreover, a thorough understanding of underlying physical processes and river heat fluxes is essential in the development of better and more adaptive water temperature models. Very few studies have measured river evaporation and condensation and subsequently calculated corresponding heat fluxes in small tributary streams, mainly because microclimate data (data collected within the stream environment) are essential and rarely available. As such, the present study will address these issues by measuring river evaporation and condensation in tributary 1 (Trib 1, a small tributary within Catamaran Brook) using floating minipans. The latent heat flux and other important fluxes were calculated. Results showed that evaporation was low within the small Trib 1 of Catamaran Brook, less than 0.07 mm day^?1. Results showed that condensation played an important role in the latent heat flux. In fact, condensation was present during 34 of 92 days (37%) during the summer, which occurred when air temperature was greater than water temperature by 4–6 °C. Heat fluxes within this small stream showed that solar radiation dominated the heat gains and long‐wave radiation dominated the heat losses. © 2015 Her Majesty the Queen in Right of Canada. Hydrological Processes. © 2015 John Wiley & Sons, Ltd.