SESTET: A spatially explicit stream temperature model based on equilibrium temperature

Stream-water temperature is a key variable controlling chemical, biological, and ecological processes in freshwater environments. Most models focus on a single river cross-section; however, temperature gradients along stretches and tributaries of a river network are crucial to assess ecohydrological features such as aquatic species suitability, growth and feeding rates, or disease transmission. We propose SESTET, a deterministic, spatially explicit stream temperature model for a whole river network, based on water and energy budgets at a reach scale and requiring only commonly available spatially distributed datasets, such as morphology and air temperature, as input. Heat exchange processes at the air–water interface are modelled via the widely used equilibrium temperature concept, whereas the effects of network structure are accounted for through advective heat fluxes. A case study was conducted on the prealpine Wigger river (Switzerland), where water temperatures have been measured in the period 2014–2018 at 11 spatially distributed locations. The results show the advantages of accounting for water and energy budgets at the reach scale for the entire river network, compared with simpler, lumped formulations. Because our approach fundamentally relies on spatially distributed air temperature fields, adequate spatial interpolation techniques that account for the effects of both elevation and thermal inversion in air temperature are key to a successful application of the model. SESTET allows the assessment of the magnitude of the various components of the heat budget at the reach scale and the derivation of reliable estimates of spatial gradients of mean daily stream temperatures for the whole catchment based on a limited number of conveniently located (viz., spanning the largest possible elevation range) measuring stations. Moreover, accounting for mixing processes and advective fluxes through the river network allows one to trust regionalized values of the parameters controlling the relationship between equilibrium and air temperature, a key feature to generalize the model to data-scarce catchments.     

Using LM‐OSL of quartz to distinguish sediments derived from surface‐soil and channel erosion

This study describes the use of linearly modulated optically stimulated luminescence (LM‐OSL) to distinguish surface‐soil derived sediments from those derived from channel bank erosion. LM‐OSL signals from quartz extracted from 15 surface‐soil and five channel bank samples were analysed and compared to signals from samples collected from two downstream river sites. Discriminant analysis showed that the detrapping probabilities of fast, first slow and second slow components of the LM‐OSL signal can be used to differentiate between the samples collected from the channel bank and surface‐soil sources. We show that for each of these source end members these components are all normally distributed. These distributions are then used to estimate the relative contribution of surface‐soil derived and channel bank derived sediment to the river bed sediments. The results indicate that channel bank derived sediments dominate the sediment sources at both sites, with 90.1 ± 3% and 91.9 ± 1.9% contributions. These results are in agreement with a previous study which used measurements of ¹³⁷Cs and ²¹⁰Pb_ex fallout radionuclides to estimate the relative contribution from these two sources. This result shows that LM‐OSL may be a useful method, at least in the studied catchment, to estimate the relative contribution of surface soil and channel erosion to river sediments. However, further research in different settings is required to test the difference of OSL signals in distinguishing these sediment sources. And if generally acceptable, this technique may provide an alternative to the use of fallout radionuclides for source tracing. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrological regulation of chemical weathering and dissolved inorganic carbon biogeochemical processes in a monsoonal river

To better understand the mechanisms relating to hydrological regulations of chemical weathering processes and dissolved inorganic carbon (DIC) behaviours, high-frequency sampling campaigns and associated analyses were conducted in the Yu River, South China. Hydrological variability modifies the biogeochemical processes of dissolved solutes, so major ions display different behaviours in response to discharge change. Most ions become diluted with increasing discharge because of the shortened reactive time between rock and water under high-flow conditions. Carbonate weathering is the main source of major ions, which shows strong chemostatic behaviour in response to changes in discharge. Ions from silicate weathering exhibit a significant dilution effect relative to the carbonate-sourced ions. Under high temperatures, the increased soil CO₂ influx from the mineralisation of organic material shifts the negative carbon isotope ratios of DIC (δ¹³C_DIC) during the high-flow season. The δ¹³C_DIC values show a higher sensitivity than DIC contents in response to various hydrological conditions. Results from a modified isotope-mixing model (IsoSource) demonstrate that biological carbon is a dominant source of DIC and plays an important role in temporal carbon dynamics. Furthermore, this study provides insights into chemical weathering processes and carbon dynamics, highlighting the significant influence of hydrological variability to aid understanding of the global carbon cycle.     

A method for evaluating the longitudinal functional connectivity of a river–lake–marsh system and its application in China

The longitudinal functional connectivity of a river–lake–marsh system (RLMS) refers to the actual water-mediated transport of material from upstream to downstream areas along a spatial gradient and is fundamental to understand hydrological and biogeochemical cycles. However, due to a lack of consensus on appropriate data and methods, the quantification of connectivity is still a challenge, especially at the catchment scale. We developed a new method to evaluate longitudinal functional connectivity based on fluxes of materials (water, sediment, and chemicals) along a RLMS. The calculation of fluxes is based on the longitudinal pattern of terrain gradient, which influences transport efficiency, and on contributions from hillslopes, which set the initial spatial template of material loading to the RLMS. We evaluate the contributions from hillslopes to RLMS based on a new modified version of the index of sediment connectivity (IC) proposed by Borselli et al. (2008) and revised by Chartin et al. (2017).We applied this method to the Baiyangdian Basin covering an area of 3.4 × 10⁴ km² in China and quantified longitudinal functional connectivity during normal, wet, and dry periods(April, July and December) in year 2016. We found that areas with good structural connectivity exhibited poor functional connectivity during the normal and dry periods. Modelling testing with discharge data from hydrological stations and measured chemicals from Baiyangdian Lake was satisfactory in test periods. We conclude that public data and Digital Elevation Model-derived information can be used to reliably map the longitudinal functional connectivity of RLMSs. The proposed method provides a useful tool for monitoring and restoring the longitudinal functional connectivity of RLMSs and our results indicate that efforts aimed at restoring functional connectivity in RLMSs should take into account landscape patterns that can greatly influence fluxes in the watershed.     

Hydraulics of overland flow influenced by litter incorporation under extreme rainfall

Water erosion on hillslopes is a worldwide environmental problem, which is a rainfall‐induced process, especially extreme rainfall. The great intensity of extreme rainfall strongly enhances the power of overland flow to detach soil and transport sediment. Plant litter is one of the most important constituents of ecosystems that often covers the soil surface and can be incorporated into topsoil. However, little attention has been paid to its effect on flow hydraulics owing to the veiled nature. This study aimed to examine the effects of incorporated litter on the hydraulic properties under extreme rainfall condition. To reach this goal, six litter rates of 0, 0.05, 0.10, 0.20, 0.35, and 0.50 kg m^?2 and four litter types collected from deciduous trees, coniferous trees, shrubs, and herbs were incorporated into topsoil. Then, simulated rainfall experiments were performed on five slope gradients (5°, 10°, 15°, 20°, and 25°) with an extreme rainfall intensity of 80 mm h^?1. The results showed that Froude number and flow velocity of the overland flow decreased, whereas flow resistance increased exponentially with litter incorporation rate. Litter type had an influence on flow hydraulics, which can mainly be attributed to the variations in surface coverage of the exposed litter and the litter morphology. Flow velocity and Darcy–Weisbach coefficient increased markedly with slope gradient. However, the variation of slope gradient did not modify the relationships between flow hydraulics and incorporated litter rate. The random roughness, resulting from heterogeneous erosion due to the uneven protection of surface exposed litter, increased linearly with litter incorporated rate. As rainfall proceeded, flow hydraulics varied with incorporated litter rate and slope gradient complicatedly due to the increases in flow rate and coverage of the exposed litter and the modification of soil surface roughness.     

Impacts of climate variability on water quality with best management practices in sub‐tropical climate of USA

Efficiency of non‐point source pollution control methods may be altered in future climate. This study investigated climate change impacts on sediment and nutrient transport, and efficiency of best management practices (BMPs), in the Upper Pearl River Watershed (UPRW) in Mississippi. The Soil and Water Assessment Tool was applied to the UPRW using observed flow, sediment and nutrient data. Water quality samples were collected at three US geological survey gauging stations. The model was successfully calibrated and validated for daily time steps (Nash Sutcliffe efficiency and coefficient of determination – R² up to 0.7) using manual and automatic (sequential uncertainty fitting version 2) methods from February 2010 to May 2011. Future weather scenarios were simulated using the LARS‐WG model, a stochastic weather generator, with Community Climate System Model, global climate model, which was developed by the National Center for Atmospheric Research in the USA. On the basis of the Special Report on Emissions Scenarios A1B, A2 and B1 of the Intergovernmental Panel on Climate Change, climate change scenarios were simulated for the mid (2046–2065) and late (2080–2099) century. Effectiveness of four BMPs (Riparian buffer, stream fencing, sub‐surface manure applications and vegetative filter strips) on reducing sediment and nutrient were evaluated in current and future climate conditions. Results show that sediment, nitrogen and phosphorus loadings will be increased up to a maximum of 26.3%, 7.3% and 14.3%, respectively, in future climate conditions. Furthermore, the effectiveness of BMPs on sediment removal will be reduced in future climate conditions, and the efficiency of nitrogen removal will be increased, whereas phosphorus removal efficiency will remain unchanged. Copyright © 2013 John Wiley & Sons, Ltd.     

River discharge,land use change,and surface water quality in the Xiangjiang River,China

To compare the impacts of river discharge on the surface water quality of the Xiangjiang River in China, 12 surface water quality parameters recorded at 31 sampling sites from January 1998 to December 2008 along the river and its main tributaries were analyzed. Significantly higher concentrations of total nitrogen, ammoniacal nitrogen, and total phosphorus, and biochemical oxygen demand were observed during low‐flow periods than during high‐flow periods, implying a higher risk to local residents drinking untreated water during low‐flow periods. Pollution indexes, including the inorganic pollution index and integrated pollution index (IPI), were negatively related to impervious surface area (ISA) and cropland area (CLA) when ISA (CLA) was less than 160 (3000) km². However, the relationship was positive when ISA (CLA) was larger than 160 (3000) km², which provided a reasonable explanation for the observed spatial patterns of water quality. Distinct increasing temporal trends for two kinds of pollution indexes were also found. The annual ISA was significantly related to the rapid degradation of water quality from 1998 to 2008, with correlation coefficient (r) values of 0.816 (p = 0.002) and 0.711 (p = 0.014) for the organic pollution index (OPI) and IPI, respectively. However, annual rainfall was negatively correlated with the two indexes with r values of 0.785 (p = 0.002) and 0.448 (p = 0.093) for OPI and IPI, respectively. Our study highlights that decision makers should be more aware of recent increases in the pollution of the Xiangjiang River, especially at downriver sites and during low‐flow periods. Copyright © 2013 John Wiley & Sons, Ltd.     

Temporal variation of suspended sediment transport in the Koga catchment,North Western Ethiopia and environmental implications

Event sediment transport and yield were studied for 45 events in the upstream part of the 260 km² agricultural Koga catchment that drains to an irrigation reservoir. Discharge and turbidity data were collected over a period of more than a year, accompanied by grab sampling. Turbidity was very well correlated with the sediment concentrations from the samples (r = 0.99), which allowed us to estimate the temporal patterns of sediment concentrations within events. The hysteresis patterns between discharge and sediment concentrations were analysed to provide insight into the different sediment sources. Anticlockwise patterns are the dominant hysteresis patterns in the area, suggesting smaller contributions of suspended sediment from the river channels than from the hillslopes and agricultural areas. Complicated types of hysteresis patterns were mostly observed for long events with multiple peaks. For a given discharge, sediment yields in August and September, when the catchment was almost completely covered with vegetation, were much smaller than during the rest of the rainy season. The hysteresis patterns and timing suggest that the sediment availability from the agricultural areas and hillslopes affects sediment yields more strongly than does peak discharge. Two distinct types of sediment rating curves were observed for the season when the agricultural land was covered with vegetation and when it was not, indicating the dominating contribution of land use/cover to sediment yields in the catchment. The rate of suspended sediment transport in the area was estimated as 25.6 t year^?1 ha^?1. Copyright © 2013 John Wiley & Sons, Ltd.     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

Rainstorm-generated sediment yield model based on soil moisture proxies (SMP)

This study develops improved Soil Moisture Proxies (SMP) based suspended sediment yield (SMPSY) models corresponding to three antecedent moisture conditions (AMCs) (i.e., AMC-I-AMC-III) by coupling the improved initial abstraction (I_a-λ) model, the SMA procedure and the SMP concept for modelling the rainfall generated suspended sediment yield. The SMPSY models specifically incorporate a watershed storage index (S) model to accentuate the transformation from storm to storm and to avoid the sudden jumps in sediment yield computation. The workability of the SMPSY models is tested using a large dataset of rainfall and sediment yield (98 storm events) from twelve small watersheds and a comparison has been made with the existing MSY model. The goodness-of-fit (GOF) statistics is evaluated in terms of the Nash Sutcliffe efficiency (NSE), and error indices, i.e., root mean square error (RMSE), normalized root mean square error (nRMSE), standard error (SE), mean absolute error (MAE), and RMSE-observations standard deviation ratio (RSR). The NSE values vary from 74.31% to 96.57% and from 75.21% to 91.78%, respectively for the SPMSY and MSY model. The NSE statistics indicate that the SMPSY model has lower uncertainty in simulating sediment yield as compared to the MSY model. The error indices are lower for the SMPSY model than the MSY model for most of the watersheds. These results show that the SMPSY model has less uncertainty and performs better than the MSY model. A sensitivity analysis of the SMPSY model shows that the parameter β is most sensitive followed by parameter S, α and A. Overall, the results show that the characterization of soil moisture variability in terms of SMPs and incorporation of improved delivery ratio and runoff coefficient relationship improves the simulation of the erosion and sediment yield generation process.