Numerical simulation of the impact of sediment supply and streamflow variations on channel grain sizes and Chinook salmon habitat in mountain drainage networks

Climatically driven changes in streamflow and hillslope sediment supply could potentially alter stream surface grain size distribution patterns and thereby impact habitat for a number of threatened and endangered in‐stream fish species. Relatively little is known about hydrograph (shape, peak flow) influence or the relative importance of chronic and episodic hillslope inputs on channel conditions. To better understand these external drivers, we calculated sediment routing through a gravel‐bedded river network using a one‐dimensional (1D) bedload transport model. We calculated changes in grain sizes and estimated Chinook salmon habitat suitability caused by a dry year and an extreme flood hydrograph, and chronic (diffusive, overland flow) or pulse (landslide, debris flow) hillslope sediment supplies. To obtain accurate channel conditions, a relatively high reference Shields stress, representative of steep mountain streams, was needed. An extreme event flood without any hillslope sediment inputs caused widespread bed coarsening and a decrease in aquatic habitat. Chronic sediment input combined with this hydrograph eliminated any changes in grain size and habitat, although when combined with a dry year flow, caused systematic bed fining. The influence of a given hydrograph therefore highly depends on the hillslope sediment supply. Regardless of the flow hydrograph or sediment pulse timing, grain size distribution or location, pulse sediment inputs did not cause widespread grain size changes despite being 100 times the total chronic input volume. Widespread and continuous hillslope sediment inputs may influence channel grain sizes and aquatic habitat more than a single discrete sediment pulse. Depending on the magnitudes of flow hydrograph and sediment supply alterations, climate change may induce no differences in grain sizes or very dramatic changes with significant consequences for long‐term sustainability. Copyright © 2013 John Wiley & Sons, Ltd.     

The impact of exceptional events on erosion,bedload transport and channel stability in a step‐pool channel

Sediment transport in the Erlenbach, a small stream with step‐pool morphology in the canton of Schwyz, Switzerland, has been monitored for more than 20 years. During this time three exceptional events (events with high sediment yield and long return times that have a large effect on channel morphology) have impacted the stream and partly or completely rearranged the existing step‐pool morphology. In the aftermath of the events, sediment transport rates at a given discharge and total sediment yield remained elevated for about a year or longer. For the last event, dated on the 20 June 2007, observations of boulder mobility and step destruction were used to interpret channel stability. Boulders with median diameters of up to 135 cm and estimated weights of more than 2·5 tons have moved during the 2007 event. Using hydraulic observations and shear stress calculations boulders up to 65 cm in diameter were predicted to have been fully mobile in peak conditions, even if form resistance and increased critical stresses needed for the initiation of motion in steep streams were taken into account. For two of the events, estimated peak shear stresses at the bed exceeded 1000 Pa, calculated both from observations of the flow hydraulics and from boulder mobility. This suggests that highly energetic flows occur relatively frequently in small, steep streams and that large boulders can be transported by fluvial processes in such streams. The observations have potential significance for hazard risk mitigation, stream engineering and restoration. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluating the performance of topobathymetric LiDAR to support multi-dimensional flow modelling in a gravel-bed mountain stream

Stream biophysical processes are commonly studied using multi-dimensional numerical modelling that quantifies flow hydraulics from which parameters such as habitat suitability, stream carrying capacity, and bed mobility are derived. These analyses would benefit from accurate high-resolution stream bathymetries spanning tens of kilometres of channel, especially in small streams or where navigation is difficult. Traditional ground-based survey methods are limited by survey time, dense vegetation and stream access, and are usually only feasible for short reaches. Conversely, airborne topobathymetric LiDAR surveys may overcome these limitations, although limited research is available on how errors in LiDAR-derived digital elevation models (DEMs) might propagate through flow models. This study investigated the performance of LiDAR-derived topobathymetry in support of multi-dimensional flow modelling and ecohydraulics calculations in two gravel-bedded reaches (approximately 200 m long), one morphologically complex and one morphologically simple, and at the segment scale (32 km-long stream segment) along a 15 m-wide river in central Idaho, USA. We compared metre and sub-metre-resolution DEMs generated from RTK-GPS ground and Experimental Advanced Airborne Research LiDAR-B (EAARL-B) surveys and water depths, velocities, shear stresses, habitat suitability, and bed mobility modelled with two-dimensional (2D) hydraulic models supported by LiDAR and ground-surveyed DEMs. Residual statistics, bias (B), and standard deviation (SD) of the residuals between depth and velocity predicted from the model supported by LiDAR and ground-survey topobathymetries were up to −0.04 (B) and 0.09 m (SD) for depth and −0.09 (B) and 0.20 m s⁻¹ (SD) for velocity. The accuracy (B = 0.05 m), precision (SD = 0.09 m), and point density (1 point m⁻²) of the LiDAR topobathymetric survey (regardless of reach complexity) were sufficient to support 2D hydrodynamic modelling and derivative stream habitat and process analyses, because these statistics were comparable to those of model calibration with B = 0 m and SD = 0.04 m for water surface elevation and B = 0.05 m s⁻¹ and SD = 0.22 m s⁻¹ for velocity in our investigation. © 2020 John Wiley & Sons, Ltd.     

Drivers of peatland water table dynamics in the central Andes,Bolivia and Peru

Cushion plant dominated peatlands are key ecosystems in tropical alpine regions of the Andes in South America. The cushion plants have formed peat bodies over thousands of years that fill many valley bottoms, and the forage produced by the plants is critical for native and nonnative domesticated mammals. The sources and flow paths of water supporting these peatlands remain largely unknown. Some studies have suggested that glacier meltwater streams support some peatlands, and that the ongoing loss of glaciers and their meltwaters could lead to the loss or diminishment of peatlands. We analysed the hydrologic regime of 10 peatlands in four mountain regions of Bolivia and Peru using groundwater monitoring. Groundwater levels in peatlands were relatively stable and within 20 cm of the ground surface during the rainy season, and many sites had water tables 40–90 cm below the ground surface in the dry season. Topographic and groundwater elevations in the peatlands demonstrated that the water source of all 10 peatlands was hillslope groundwater flowing from lateral moraines, talus, colluvium, or bedrock aquifers into the peatlands. There was little to no input from streams, whether derived from glacier melt or other sources, and glacier melt could not have recharged the hillslope aquifers supporting peatlands. We measured the stable water isotopes in water samples taken during different seasons, distributed throughout the catchments, and the values are consistent with this interpretation. Our findings indicate that peatlands in the study region are recharged by hillslope groundwater discharge rather than stream water and may not be as vulnerable to glacial decline as other studies have indicated. However, both glaciers and peatlands are susceptible to changing thermal and precipitation regimes that could affect the persistence of peatlands.     

Application of 210Pbex inventories to measure net hillslope erosion at burned sites

Fallout radionuclides, including lead‐210 excess (²¹⁰Pb_ex), have been broadly and successfully used to quantify net hillslope sediment transport in agricultural, pastoral and forested landscapes but have only recently been applied in burned terrain. Quantifying post‐fire erosion is important because fires can amplify hillslope erosion, impacting terrestrial and aquatic habitat and water quality. However, we lack a basic understanding of the fate of ²¹⁰Pb_ex in fires. To address this knowledge gap, we collected over 400 soil samples from unburned, moderately and severely burned forested sites in central Idaho. We measured soil ²¹⁰Pb_ex content at stable reference and eroding sites and in mineral and organic soil components. At all sites, organic matter had the highest concentration of ²¹⁰Pb_ex, representing 30% to 73% of the total activity. At the severely and moderately burned sites, ²¹⁰Pb_ex reference inventories were lower by 58% and 41%, with about 40% less organic mass, relative to the unburned site. These results indicate that most ²¹⁰Pb_ex in our semi‐arid, forested sites was bound to organic matter, and that a substantial portion of this lead was lost due to forest fires. These losses likely occurred through volatilization and wind transport of smoke and ash. In the moderately burned site, ²¹⁰Pb_ex losses were more spatially variable, potentially due to spatially uneven fire intensity and effects. Despite equal percent losses of ²¹⁰Pb_ex, lower inventories at the burned sites produced lower calculated net erosion rates relative to the unburned site. Thus, given methodological uncertainties, ²¹⁰Pb_ex losses due to fire, and the subsequent sensitivity of calculated net erosion rates to these lower ²¹⁰Pb_ex inventories, we suggest this method should not be used in burned terrain to calculate absolute net erosion and deposition rates. However, within a given burned site, ²¹⁰Pb_ex inventories still provide useful information describing relative soil losses and storage across the landscape. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison of alternative representations of hydraulic-conductivity anisotropy in folded fractured-sedimentary rock: modeling groundwater flow in the Shenandoah Valley (USA)

A numerical representation that explicitly represents the generalized three-dimensional anisotropy of folded fractured-sedimentary rocks in a groundwater model best reproduces the salient features of the flow system in the Shenandoah Valley, USA. This conclusion results from a comparison of four alternative representations of anisotropy in which the hydraulic-conductivity tensor represents the bedrock structure as (model A) anisotropic with variable strikes and dips, (model B) horizontally anisotropic with a uniform strike, (model C) horizontally anisotropic with variable strikes, and (model D) isotropic. Simulations using the US Geological Survey groundwater flow and transport model SUTRA are based on a representation of hydraulic conductivity that conforms to bedding planes in a three-dimensional structural model of the valley that duplicates the pattern of folded sedimentary rocks. In the most general representation, (model A), the directions of maximum and medium hydraulic conductivity conform to the strike and dip of bedding, respectively, while the minimum hydraulic-conductivity direction is perpendicular to bedding. Model A produced a physically realistic flow system that reflects the underlying bedrock structure, with a flow field that is significantly different from those produced by the other three models.     

Origin of halite brine in the Onondaga Trough near Syracuse,New York State,USA: modeling geochemistry and variable-density flow

Halite brine (saturation ranging from 45 to 80%) lies within glacial sediments that fill the Onondaga Trough, a bedrock valley deepened by Pleistocene glaciation near Syracuse, New York State, USA. The most concentrated brine occupies the northern end of the trough, about 10 km downgradient of the northern limit of halite beds in the Silurian Salina Group, the assumed source of salt. The chemical composition of the brine and its radiocarbon age suggest that the brine originally formed about 16,700  years ago through dissolution of halite by glacial melt water and later mixed with saline bedrock water. Two hypotheses regarding the formation of the brine pool were tested through variable-density flow simulations using SEAWAT. Simulation results supported the first hypothesis that the brine pool was derived from a source in the glacial sediments and then migrated to its current position, where it has persisted for over 16,000  years. A second hypothesis that the brine pool formed through steady accumulation of brine from upward flow of a source in the underlying bedrock was not supported by simulation results, because the simulated age distribution was much younger than the age estimated from geochemical modeling.     

Carbon to nitrogen (C:N) stoichiometry of the spring–summer phytoplankton bloom in the North Water Polynya (NOW)

The carbon to nitrogen (C:N) stoichiometry of phytoplankton production varied significantly during the spring–summer bloom in the North Water Polynya (NOW), from April through July 1998. The molar ratio of particulate organic carbon (POC) to nitrogen (PON) production by phytoplankton (ΔPOC:ΔPON) increased from 5.8 during April through early June to 8.9 in late June and July. The molar dissolved inorganic carbon (DIC) to nitrate+nitrite (NO₃) drawdown ratio (ΔDIC: ΔNO₃) increased from 6.7 in April and May, to 11.9 in June (no estimate for July because of ice melting). The discrepancy between ΔPOC:ΔPON and ΔDIC:ΔNO₃ was likely due to dissolved organic carbon (DOC) production. Increased ΔPOC:ΔPON of phytoplankton and surface water ΔDIC:ΔNO₃ throughout the phytoplankton blooms resulted from changes in physical properties of the upper water column, such as reduced thickness of the surface mixed layer that exposed phytoplankton to increased photosynthetically available radiation (PAR), accompanied by NO₃ depletion. This is expected to have significant effects on the cycling of carbon (C) and nitrogen (N) in pelagic ecosystems, as the increased C:N ratio of organic matter decreases its quality as substrate for grazers and microbial communities. Based on ΔPOC:ΔPON, the ratio of POC to chlorophyll a (Chl) production (ΔPOC:ΔChl) and the relationship between Chl yields and NO₃ depletion, we estimate that 71±17% and 46±20% of the depleted NO₃ went to PON production in the euphotic zone over the polynya from April to early June, and late June to July, respectively. The remaining NO₃ was likely channelled to dissolved organic nitrogen (DON) and heterotrophic bacteria, which were not returned to the dissolved inorganic nitrogen (DIN) pool through recycling during the course of the study. Hence, the autotrophic production of organic N and its recycling by the microbial food web were not coupled temporally.