GIS‐based modelling of soil erosion processes using the modified‐MMF (MMMF) model in a large watershed having vast agro‐climatological differences

The present study demonstrates a spatially distributed application of a field‐scale annual soil loss model, the modified‐MMF (MMMF), to a large watershed using hydrological routing techniques, remote sensing data and geospatial technologies. In this study, the MMMF model is implemented after incorporating the corrections suggested in recent literature along with appropriate modifications of the model to suit the agro‐climatological conditions prevailing in most parts of India. Sensitivity analysis carried out through an Average Linear Sensitivity approach indicates that the model outputs are highly sensitive to soil moisture (MS), bulk density (BD), effective hydraulic depth (EHD), ground cover (GC) and settling velocity for clay (VS_c). During calibration and validation, the performance evaluation statistics are mostly in the range of very good to satisfactory for both runoff and soil loss at the watershed outlet. Even spatial validation of the results of intermediate processes in the water phase and the sediment phase, although qualitative, seems to be reasonable and rational. Furthermore, the soil erosion severity analysis for different land‐uses existing in the watershed indicates that about 90% of the watershed area, especially that occupied by agricultural lands, is vulnerable to the long‐term effects of soil erosion. Copyright © 2018 John Wiley & Sons, Ltd.     

Gully head modelling: A Mediterranean badland case study

Predicting the location of gully heads in various environments is an important step towards predicting gully erosion rates. So far, field data collection and modelling of topographic thresholds for gully head development has mainly focused on gullies that formed in forested areas, rangelands, pastures and cropland. Such information for gullies in badlands however is very scarce. Therefore, this paper aims to extend the database on gully head topographical thresholds through data collection in a badland area and to improve the prediction of gully heads forming at sites with a very low erosion resistance value. For this, we chose a badland site located in central Italy that is characterized by biancana forms and both active and dormant gullies. The definition of the conditions under which present‐day gully heads developed allowed a better modelling of the gully head threshold equation, with modification of a previous model and the exemplification of how to use the updated model. The model shows that the resistance to gully head retreat depends on slope gradient and drainage area at gully heads, land use at the moment of gully development (as numerically expressed using parameters derived from the Runoff Curve Number method), surface rock fragment cover, presence of joints, pipes, and factors/processes affecting detachment rate. This study attempted to better understand environmental conditions that control the development of gully heads in badlands through a combination of field data collection of gully heads, an analysis of land use changes over 10 centuries, focusing on the period 1820–2005, and land use management through repeat photography and a critical examination of historical documents. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Removal of water‐surface reflection effects with a temporal minimum filter for UAV‐based shallow‐water photogrammetry

The recent development of structure‐from‐motion (SfM) and multi‐view stereo (MVS) photogrammetry techniques has enabled semi‐automatic high‐resolution bathymetry using aerial images taken by consumer‐grade digital cameras mounted on unmanned aerial vehicles (UAVs). However, the applicability of these techniques is sometimes limited by sun and sky reflections at the water surface, which render the point‐cloud density and accuracy insufficient. In this research, we present a new imaging technique to suppress the effect of these water‐surface reflections. In this technique, we order a drone to take a short video instead of a still picture at each waypoint. We then apply a temporal minimum filter to the video. This filter extracts the smallest RGB values in all the video frames for each pixel, and composes an image with greatly reduced reflection effects. To assess the performance of this technique, we applied it at three small shallow‐water sites. Specifically, we evaluated the effect of the technique on the point cloud density and the accuracy and precision of the photogrammetry. The results showed that the proposed technique achieved a far denser point cloud than the case in which a randomly chosen frame was used for each waypoint, and also showed better overall accuracy and precision in estimating water‐bottom elevation. The effectiveness of this new technique should depend on the surface wave state and sky radiance distribution, and this dependence, as well as the applicability to large areas, should be investigated in future research. Copyright © 2018 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of hillslope–channel coupling and implications for the longitudinal profile in a dryland basin

The long‐term evolution of channel longitudinal profiles within drainage basins is partly determined by the relative balance of hillslope sediment supply to channels and the evacuation of channel sediment. However, the lack of theoretical understanding of the physical processes of hillslope–channel coupling makes it challenging to determine whether hillslope sediment supply or channel sediment evacuation dominates over different timescales and how this balance affects bed elevation locally along the longitudinal profile. In this paper, we develop a framework for inferring the relative dominance of hillslope sediment supply to the channel versus channel sediment evacuation, over a range of temporal and spatial scales. The framework combines distinct local flow distributions on hillslopes and in the channel with surface grain‐size distributions. We use these to compute local hydraulic stresses at various hillslope‐channel coupling locations within the Walnut Gulch Experimental Watershed (WGEW) in southeast Arizona, USA. These stresses are then assessed as a local net balance of geomorphic work between hillslopes and channel for a range of flow conditions generalizing decadal historical records. Our analysis reveals that, although the magnitude of hydraulic stress in the channel is consistently higher than that on hillslopes, the product of stress magnitude and frequency results in a close balance between hillslope supply and channel evacuation for the most frequent flows. Only at less frequent, high‐magnitude flows do channel hydraulic stresses exceed those on hillslopes, and channel evacuation dominates the net balance. This result suggests that WGEW exists mostly (~50% of the time) in an equilibrium condition of sediment balance between hillslopes and channels, which helps to explain the observed straight longitudinal profile. We illustrate how this balance can be upset by climate changes that differentially affect relative flow regimes on slopes and in channels. Such changes can push the long profile into a convex or concave condition. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Assessment of soil factors controlling ephemeral gully erosion on agricultural fields

The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion. Copyright © 2018 John Wiley & Sons, Ltd.     

Bacteria‐sediment associations in an alpine,freshwater environment and their effects on particle size,density and settling velocity

This study measures the presence of bacteria‐sediment associations (BSAs) in an alpine, glacier‐fed watershed in the Southern Coast Mountains of British Columbia, Canada. The impact of BSAs on the creation of flocculated particles and their settling velocity are quantified using a laser transmissometer. Results from the study indicate that BSAs are present in the watershed and vary over both space and time. The percentage of bacteria associated with sediment particles was found to range from < 1% to 40%. Major sources of planktonic bacteria such as agricultural land and wastewater treatment outflow co‐occur with large decreases in the BSA ratio. Laboratory analysis demonstrates that an increase in the concentration of bacteria was associated with a decrease in the volume concentration of small particles, and a decrease in both estimated density and measured settling velocity for particles in larger size classes; consistent with flocculated particles of increasing complexity arising from combinations of primary particles and/or BSAs. Paleoenvironmental reconstructions using laminated lake sediments in alpine, glacier‐fed systems benefit from a fuller understanding of the geomorphologic processes by which they formed. While bacteria are noted to enhance formation of flocculated particles in laboratory systems, their impact upon geomorphic processes in natural systems have yet to be fully explored. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Mesocosm experiments identifying hotspots of groundwater upwelling in a water column by fibre optic distributed temperature sensing

Lacustrine groundwater discharge (LGD) can substantially impact ecosystem characteristics and functions. Fibre optic distributed temperature sensing (FO‐DTS) has been successfully used to locate groundwater discharge into lakes and rivers at the sediment–water interface, but locating groundwater discharge would be easier if it could be detected from the more accessible water surface. So far, it is not clear if how and under which conditions the LGD signal propagates through the water column to the water surface–atmosphere interface, and what perturbations and signal losses occur along this pathway. In the present study, LGD was simulated in a mesocosm experiment. Under winter conditions, water with temperatures of 14 to 16 °C was discharged at the bottom of a 10 × 2.8‐m mesocosm. Water within this mesocosm ranged from 4.0 to 7.4 °C. Four layers (20, 40, 60, and 80 cm above the sediment) of the 82 cm deep mesocosm were equipped with FO‐DTS for tracing thermal patterns in the mesocosm. Aims are (a) to test whether the positive buoyancy of relatively warm groundwater imported by LGD into shallow water bodies allows detection of LGD at the lake's water surface–atmosphere interface by FO‐DTS, (b) to analyse the propagation of the temperature signal from the sediment‐water interface through the water column, and (c) to learn more about detectability of the signal under different discharge rates and weather conditions. The experiments supported the benchmarking of scale dependencies and robustness of FO‐DTS applications for measuring upwelling into aquatic environments and revealed that weather conditions can have important impacts on the detection of upwelling at water surface–atmosphere interfaces at larger scales.     

Continuous real‐time analysis of the isotopic composition of precipitation during tropical rain events: Insights into tropical convection

To investigate stable isotopic variability of precipitation in Singapore, we continuously analysed the δ‐value of individual rain events from November 2014 to August 2017 using an online system composed of a diffusion sampler coupled to Cavity Ring‐Down Spectrometer. Over this period, the average value (δ¹⁸O_Avg), the lowest value (δ¹⁸O_Low), and the initial value (δ¹⁸O_Init) varied significantly, ranging from ?0.45 to ?15.54‰, ?0.9 to ?17.65‰, and 0 to ?13.13‰, respectively. All 3 values share similar variability, and events with low δ¹⁸O_Low and δ¹⁸O_Avg values have low δ¹⁸O_Init value. Individual events have limited intraevent variability in δ‐value (Δδ) with the majority having a Δδ below 4‰. Correlation of δ¹⁸O_Low and δ¹⁸O_Avg with δ¹⁸O_Init is much higher than that with Δδ, suggesting that convective activities prior to events have more control over δ‐value than on‐site convective activities. The d‐excess of events also varies considerably in response to the seasonal variation in moisture sources. A 2‐month running mean analysis of δ¹⁸O reveals clear seasonal and interannual variability. Seasonal variability is associated with the meridional movement of the Intertropical Convergence Zone and evolution of the Asian monsoon. El Niño–Southern Oscillation is a likely driver of interannual variability. During 2015–2016, the strongest El Niño year in recorded history, the majority of events have a δ¹⁸O value higher than the weighted average δ¹⁸O of daily precipitation. δ¹⁸O shows a positive correlation with outgoing longwave radiation in the western Pacific and the Asian monsoon region, and also with Oceanic Niño Index. During El Niño, the convection centre shifts eastward to the central/eastern Pacific, weakening convective activities in Southeast Asia. Our study shows that precipitation δ‐value contains information about El Niño–Southern Oscillation and the Intertropical Convergence Zone, which has a significant implication for the interpretation of water isotope data and understanding of hydrological processes in tropical regions.     

Spatio‐temporal tracer variability in the glacier melt end‐member — How does it affect hydrograph separation results?

Geochemical and isotopic tracers were often used in mixing models to estimate glacier melt contributions to streamflow, whereas the spatio‐temporal variability in the glacier melt tracer signature and its influence on tracer‐based hydrograph separation results received less attention. We present novel tracer data from a high‐elevation catchment (17 km², glacierized area: 34%) in the Oetztal Alps (Austria) and investigated the spatial, as well as the subdaily to monthly tracer variability of supraglacial meltwater and the temporal tracer variability of winter baseflow to infer groundwater dynamics. The streamflow tracer variability during winter baseflow conditions was small, and the glacier melt tracer variation was higher, especially at the end of the ablation period. We applied a three‐component mixing model with electrical conductivity and oxygen‐18. Hydrograph separation (groundwater, glacier melt, and rain) was performed for 6 single glacier melt‐induced days (i.e., 6 events) during the ablation period 2016 (July to September). Median fractions (±uncertainty) of groundwater, glacier melt, and rain for the events were estimated at 49±2%, 35±11%, and 16±11%, respectively. Minimum and maximum glacier melt fractions at the subdaily scale ranged between 2±5% and 76±11%, respectively. A sensitivity analysis showed that the intraseasonal glacier melt tracer variability had a marked effect on the estimated glacier melt contribution during events with large glacier melt fractions of streamflow. Intra‐daily and spatial variation of the glacier melt tracer signature played a negligible role in applying the mixing model. The results of this study (a) show the necessity to apply a multiple sampling approach in order to characterize the glacier melt end‐member and (b) reveal the importance of groundwater and rainfall–runoff dynamics in catchments with a glacial flow regime.     

A multicatchment analysis of headwater and downstream temperature effects from contemporary forest harvesting

Stream temperature is a key physical water‐quality parameter, controlling many biological, chemical, and physical processes in aquatic ecosystems. Maintenance of cool stream temperatures during summer is critical for high‐quality aquatic habitat. As such, transmission of warm water from small, nonfish‐bearing headwater streams after forest harvesting could cause warming in downstream fish‐bearing stream reaches with negative consequences. In this study, we evaluate (a) the effects of contemporary forest management practices on stream temperature in small, headwater streams, (b) the transmission of thermal signals from headwater reaches after harvesting to downstream fish‐bearing reaches, and (c) the relative role of lithology and forest management practices in influencing differential thermal responses in both the headwater and downstream reaches. We measured summer stream temperatures both preharvest and postharvest at 29 sites—12 upstream sites (4 reference, 8 harvested) and 17 downstream sites (5 reference, 12 harvested)—across 3 paired watershed studies in western Oregon. The 7‐day moving average of daily maximum stream temperature (T_7DAYMAX) was greater during the postharvest period relative to the preharvest period at 7 of the 8 harvested upstream sites. Although the T_7DAYMAX was generally warmer in the downstream direction at most of the stream reaches during both the preharvest and postharvest period, there was no evidence for additional downstream warming related to the harvesting activity. Rather, the T_7DAYMAX cooled rapidly as stream water flowed into forested reaches ~370–1,420 m downstream of harvested areas. Finally, the magnitude of effects of contemporary forest management practices on stream temperature increased with the proportion of catchment underlain by more resistant lithology at both the headwater and downstream sites, reducing the potential for the cooling influence of groundwater.