Selecting Between One‐Dimensional and Two‐Dimensional Hydrodynamic Models for Ecohydraulic Analysis

Aquatic habitat assessment and river restoration design require geospatially explicit maps of hydraulic conditions. Diverse mechanistic ecohydraulic models compute spatially explicit depth and velocity results to evaluate habitat suitability spatially as a function of these abiotic conditions. This study compared depth and velocity results from two‐dimensional (2D) and one‐dimensional (1D) hydraulic models with algorithms that laterally discretize 1D velocity and interpolate depth and velocity spatially based on the Laplacian heat mapping approach. These ‘conveyance distributed’ methods constitute ‘best 1D modelling practice’ and were compared with 2D results for the first time. The 1D and 2D models were applied to three morphologically distinct reaches (leveed, meandering, and anastomosing) for three flows (base, bankfull, and flood flows) of the partially regulated, gravel/cobble lower Yuba River in north–central California. The test metrics were the coefficient of determination (R²) and the median absolute residual ( ). These metrics quantified the incremental uncertainty 1D approximation incurs, results which make explicit cost–benefit processes of model selection possible. Finally, velocity residual maps were analysed to identify regions and processes where residuals were high, indicating divergence from the 1D assumptions. Paired data (1D–2D) fell between 0.94 ≥ R² ≥ 1.00 (R²_mean = 0.98 and R²_median = 0.99) for depth and median absolute residuals were all 3.8 ≤ ≤ 7.2% (i.e. 50% of residuals are approximately within ±1.7 to 3.6%). Higher flows and lower gradient reaches had lower residuals and higher R². Velocity diverged more, particularly for base flow in anastomosing reaches (0.42 < R² < 0.58). One‐dimensional, conveyance distributed, assumptions performed better for other channel types, where 0.69 < R² < 0.81 (R²_mean = 0.75 and R²_median = 0.77), with median absolute residuals between 9.6% > > 22.4% (i.e. ~ ± 4.6 to ±11.2%), where _mean = 14.2% and _median = 13% (~ ±7.1 and 6.5%). The conveyance distributed 1D velocity model performed best, where the orthogonal flow assumptions obtained and where side channels did not transition from backwater to conveying area between flows. Copyright © 2015 John Wiley & Sons, Ltd.     

Role of Local Flow Conditions in River Biofilm Colonization and Early Growth

Direct numerical simulations of a turbulent boundary layer flow over a bed of hemispheres of height h are performed using an immersed boundary method for comparison with river biofilm growth experiments performed in a hydraulic flume. Flow statistics above the substrates are shown to be in agreement with measurements performed by laser Doppler velocimetry and particle image velocimetry in the experiments. Numerical simulations give access to flow components inside the roughness sublayer, and biofilm colonization patterns found in the experiments are shown to be associated with low shear stress regions on the hemisphere surface. Two bed configurations, namely staggered and aligned configurations, lead to different colonization patterns because of differences in the local flow topology. Dependence with the Reynolds number of the biofilm distribution and accrual 7 days after inoculum is shown to be associated to local flow topology change and shear stress intensity. In particular, the shear stress τ on the surface of the hemispheres is found to scale as , where Re_t = u^*h/ν, with u^* as the log law friction velocity and ν as the fluid kinematic viscosity. This scaling is due to the development of boundary layers along the hemisphere surface. Associated with a critical shear stress for colonization and early growth, it explains the increasing delay in biomass accrual for increasing flow velocities in the experiments. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling movements of lake charr eggs on spawning grounds

Because little is known about the effects of currents on lake charr (Salvelinus namaycush) deposition and fate of eggs, a method to determine their role in spawning success was developed. It is suggested that an ideal strategy for lake charr egg survival is to select hydrodynamic conditions that allow for eggs to roll on lakebed and into void spaces in the substratum, whilst avoiding to be lifted from it. On that basis, egg movement on spawning grounds were modeled using an approach that combines fluid mechanics and biology. Wave modeling combined with drag and lift forces and egg parameters were used to propose both critical entrainment and resuspension velocities as well as a modular range of depths that should be ideal for lake charr spawning given a certain fetch. Modeling suggested that critical velocities of 4.1 and 10.3 cm/s would cause egg entrainment and resuspension respectively. The model was partially validated using in situ sampling on three spawning sites in Lake Témiscouata (Québec, Canada). Then, further model applications were conducted on twelve characterized Great Lakes spawning sites using data found in the literature. Relation of egg survival and the difference between resuspension depth and spawning depth at those sites was found to fit $y=66.78e^{0.13x}$ with a correlation coefficient $r^2=0.89$. The findings provide the first explanation of the impact of currents on lake charr spawning success. They provide insight to help researchers and authorities in their restorative and protection efforts for lake charr populations.     

Large wood addition for aquatic habitat rehabilitation in an incised,sand‐bed stream,Little Topashaw Creek,Mississippi

Large wood (LW) is a key component of stream habitats, and degraded streams often contain little wood relative to less‐impacted ones. Habitat rehabilitation and erosion control techniques that emphasize addition of natural wood in the form of individual elements or structures are increasingly popular. However, the efficacy of wood addition, especially in physically unstable, warmwater systems is not well established. The effects of habitat rehabilitation of Little Topashaw Creek, a sinuous, sand‐bed stream draining 37 km² in northwest Mississippi are described herein. The rehabilitation project consisted of placing 72 LW structures along eroding concave banks of a 2‐km reach and planting 4000 willow cuttings in sandbars opposite or adjacent to the LW structures. Response was measured by monitoring flow, channel geometry, physical aquatic habitat and fish populations in treated and untreated reaches for 2 years before and 4 years after rehabilitation. Initially, LW structures reduced high flow velocities at concave bank toes. Progressive failure of the LW structures and renewed erosion began during the second year after rehabilitation, with only 64% of the structures and about 10% of the willow plantings surviving for 3 years. Accordingly, long‐term changes in physical habitat attributable to rehabilitation were limited to an increase in LW density. Fish biomass increased in the treated reach, and species richness approximately doubled in all reaches after rehabilitation, suggesting the occurrence of some sort of stressful event prior to our study. Fish community composition shifted toward one typical of a lightly degraded reference site, but similar shifts occurred in the untreated reaches downstream, which had relatively high levels of naturally occurring LW. Large wood is a key component of sand‐bed stream ecosystems, but LW addition for rehabilitation should be limited to sites with more stable beds and conditions that foster rapid woody plant colonization of sediment deposits. Published in 2006 by John Wiley & Sons, Ltd.     

The role of riparian vegetation and woody debris in the development of macroinvertebrate assemblages in streams

Riparian vegetation development and macroinvertebrate assemblages were studied in 16 streams formed between 35 and 230 years ago, following glacial recession in Glacier Bay National Park, southeast Alaska. Riparian vegetation established most rapidly in streams where flow variation in downstream reaches was buffered by a lake. Riparian vegetation development was positively correlated with lower bank stability, but was independent of stream age. Roots and branches of riparian vegetation trailing into streams (trailing riparian habitat—TRH) were shown to be an important habitat for a number of macroinvertebrate taxa. In young and unstable streams, TRH was colonized mainly by Plecoptera whereas in more stable lake‐influenced streams Simuliidae dominated. Significant coarse woody debris (CWD) accumulations were not observed until after approximately 130 years of stream development had occurred when certain channel features, such as gravel bars, were stabilized by dead wood. Where dead wood was present, opportunistic wood taxa were abundant, even in the younger streams. However, a xylophagous species, Polypedilum fallax, was not recorded until streams were over 100 years old. Two‐way indicator species analysis (TWINSPAN) using presence/absence of macroinvertebrate taxa on TRH, initially divided streams into lake and non‐lake systems, but subsequent divisions were consistent with differences in stream age. TWINSPAN of macroinvertebrate assemblages on dead wood again highlighted differences in stream age. Canonical correspondence analysis indicated that bed stability and stream age were the most important environmental variables influencing macroinvertebrate distribution on TRH. Trailing riparian habitat was most abundant in moderately unstable streams where it facilitates invertebrate colonization. CWD contributes markedly to channel stabilization, provides habitat for invertebrate xylophages, and confers additional habitat complexity. Maximum levels of CWD are predicted to occur in non‐lake streams after approximately 300 years, but at least a further 100 years will be required in stable streams below lakes where dead wood entrainment is not enhanced by flooding, channel migration and bank undercutting. A conceptual model summarizing the role of TRH and CWD on stream development in Glacier Bay is presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Lake Urmia restoration success story: A natural trend or a planned remedy?

Lake Urmia is the second-largest hypersaline lake in the world. There has been a drastic water level drop of 7.2 m from 1995 to 2016. Beginning in October 2013, the Lake Urmia Restoration Plan (LURP) launched a 10-year program. An increase in water level and a relative improvement in Lake Urmia condition has been observed since 2017. It is an undecided and controversial issue whether the recent positive trend of Lake Urmia has been due to the LURP activities or it is a natural contribution of climate factors variations. To shed some light on this issue, we examine three other lakes, adjacent to the Lake Urmia basin, with similar rainfall variability to investigate their status during the same period. Van (Turkey), Mosul, and Tharthar (both in Iraq), are evaluated as well as Lake Urmia. Three decades of remotely sensed data including precipitation ($P$), water level ($\mathit{WL}$), and lake extent ($A$) were considered for the mentioned lakes. A significant correlation was observed between standardized $\mathit{WL}-P$, and $A-P$ over the long-term period, especially for the recent three years ($R$² = 0.63–0.87). We show that the cumulative precipitation in the antecedent months played a major role in the improvement of these lakes' situation but with different time lags (up to 6 months for Van and Mosul lakes and up to 36 months for Lake Urmia and Tharthar lake). These findings could inform the planners of LURP to adopt strategies for achieving a sustainable state of Lake Urmia based on a more realistic outlook.     

Nitrate removal in the hyporheic zone of a salmon river in Alaska

Pacific boreal streams and riparian zones are believed to receive significant N loads that are derived from the ocean in the form of decaying sockeye salmon (Oncorhynchus nerka). Using a small stream in south‐central Alaska we examined whether the associated riparian forest could take up the pulse of marine‐derived nitrogen (MDN) entering the hyporheic zone from spawning and dying sockeye salmon. We evaluate the relative importance of riparian uptake and denitrification in nitrate‐N removal in hyporheic sediment. We found that maximum biological removal of nitrate peaked within 1 h of water entering the hyporheic zone, decreasing exponentially with subsurface flow duration. Plant and microbial uptake reached 14 µg NO‐N L^?1 min^?1 and denitrification reached 4 µg NO‐N L^?1 min^?1 during the initial 2 h of transit time. Our results reinforce the hypothesis that MDN from Pacific salmon can be transferred to riparian zone via hyporheic flow. Most nitrate‐N removal along hyporheic flow paths is by plant and microbial uptake (the respective contributions could not be determined). Denitrifying bacteria are present and active in the hyporheic zones of this well‐oxygenated Alaskan stream but their contribution to the nitrate‐N removal is small compared to plant and microbial uptake in such nitrate‐N poor environment. Copyright © 2008 John Wiley & Sons, Ltd.     

Large wood debris recruitment on differing riparian landforms along a Gulf Coastal Plain (USA) stream: a comparison of large floods and average flows

In southeastern Coastal Plain streams, wood debris can be very abundant and is recruited from extensive forested floodplains. Despite importance of wood debris, there have been few opportunities to examine recruitment and redistribution of wood in an undisturbed setting, particularly in the southeastern Coastal Plain. Following extensive flooding in 1994, measurements of individual downed trees (species, dbh, orientation, distance from base‐flow channel and condition) were made across replicated riparian landforms in a Gulf Coastal Plain 5th‐order stream. Annually, the fate of these trees was determined and newly recruited trees were noted. More than 300 downed trees have been recorded. Recruitment varied across landforms with more constrained reaches having greatest mortality. Total tree mortality varied substantially across years. Generally, tree recruitment was greatest in years with substantial floods (1994 and 1998). For each riparian landform type, tree mortality was correlated with the maximum daily flow during the period preceding annual debris surveys. This relationship was particularly strong for sand ridges (r² = 0.942) and low terraces (r² = 0.915), but was significant for floodplains (r² = 0.413). Greatest rates of debris recruitment per maximum daily flow were observed for sand ridges followed by low terraces. Flood characteristics also influenced debris recruitment. The 1994 flood was caused by a tropical storm and resulted in a rapid rise in streamflow. Much of the debris recruited during this flood was from toppled trees and was oriented parallel to the stream channel. In contrast, the 1998 flood was preceded by a wetter than average winter with more gradually rising flows and there was no relationship between riparian landform and debris characteristics. These results indicate that wood recruitment dynamics in Coastal Plain streams are complex. Wood recruitment rates are controlled by cyclical variations in climate interacting with riparian geomorphology. Infrequent high flows appear critical in the maintenance of the instream debris pool. Copyright © 2007 John Wiley & Sons, Ltd.     

Experimental study on discharge coefficient of a gear-shaped weir

This study focused on hydraulic characteristics around a gear-shaped weir in a straight channel. Systematic experiments were carried out for weirs with two different gear heights and eight groups of geometrical parameters. The impacts of various geometrical parameters of gear-shaped weirs on the discharge capacity were investigated. The following conclusions are drawn from the experimental study: (1) The discharge coefficient ($m_{\mathrm{c}}$) was influenced by the size of the gear: at a constant discharge, the weir with larger values of a/b (a is the width of the gear, and b is the width between the two neighboring gears) and a/c (c is the height of the gear) had a smaller value of $m_{\mathrm{c}}$. The discharge capacity of the gear-shaped weir was influenced by the water depth in the weir. (2) For type C1 with a gear height of 0.01 m, when the discharge was less than 60 m³/h and $H_{\mathrm{1}}/P$ < 1.0 ($H_{\mathrm{1}}$ is the water depth at the low weir crest, and P is the weir height), $m_{\mathrm{c}}$ significantly increased with the discharge and $H_{\mathrm{1}}/P$; with further increases of the discharge and $H_{\mathrm{1}}/P$, $m_{\mathrm{c}}$ showed insignificant decreases and fluctuated within small ranges. For type C2 with a gear height of 0.02 m, when the discharge was less than 60 m³/h and $H_{\mathrm{1}}/P$ < 1.0, $m_{\mathrm{c}}$ significantly increased with the discharge and $H_{\mathrm{1}}/P$; when the discharge was larger than 60 m³/h and $H_{\mathrm{1}}/P$ > 1.0, $m_{\mathrm{c}}$ slowly decreased with the increases of the discharge and $H_{\mathrm{1}}/P$ for $a/b$ ≤ 1.0 and $a/c$ ≤ 1.0, and slowly increased with the discharge and $H_{\mathrm{1}}/P$ for $a/b$ > 1.0 and $a/c$ > 1.0. (3) A formula of $m_{\mathrm{c}}$ for gear-shaped weirs was established based on the principle of weir flow, with consideration of the water depth in the weir, the weir height and width, and the height of the gear.     

Phosphorus removal by adsorbent based on poly-aluminum chloride sludge

Phosphorus adsorption tests were carried out using poly-aluminum chloride sludge (PACS), which was collected from a water treatment plant in Nanjing. The amount of phosphorus adsorbed by PACS increased quickly within the first hour and reached equilibrium after about 48 h. The adsorption behavior of PACS for phosphorus is consistent with the Langmuir adsorption isotherm equation (R² > 0.99) and parallel first-order kinetic equation (R² > 0.98). With the increase of the PACS concentration, the adsorption capacity of PACS for phosphorus decreased, and the removal rate increased. The results of batch tests showed that the adsorption capacities of PACS for phosphorus ranged from 1.64 to 1.13 mg/g when the pH value varied from 4 to 10. However, the adsorption capacity of PACS was not evidently influenced by temperature. In comparison with the ion exchange resin, the adsorption capacity of PACS was barely inhibited by competitive ions, such as ${\text{SO}}_4^{2-}$, ${\text{NO}}_3^{-}$, and Cl⁻. The PACS surface after adsorption became smooth, and the vibration peaks of Al–O and Al–OH shifted. Both HCl and NaOH have a strong desorption effect on PACS after adsorption saturation, and with higher concentrations of HCl and NaOH, the desorption effect was stronger. Results of column adsorption experiments showed that with lower phosphorus and hydraulic loads, the adsorption column took longer to reach saturation. This indicated that PACS could be used as an efficient material for removal of phosphorus from water. This study provides a new treatment method with PACS.     

Large wood aids spawning Chinook salmon (Oncorhynchus tshawytscha) in marginal habitat on a regulated river in California

To determine whether large wood (LW, ≥1‐m length, ≥10‐cm diameter) plays a role in Chinook salmon (Oncorhynchus tshawytscha) redd (i.e. egg nest) placements in a regulated, Mediterranean‐climate, medium‐sized river (where channel width is less than the upper quartile of length of potential instream wood pieces), characteristics of 527 large wood pieces, locations of 650 redds, and mesohabitat delineations (riffle, run, glide, pool) were collected during a spawning season along a 7.7 km reach directly below Camanche Dam on the Mokelumne River, CA. LW was regularly distributed across the study reach an average 70 LW pieces km^‐1. Some LW clustering was evident at islands and meander bends. Spawners built 85% of redds within one average channel width (31 m) of LW. Spawners utilized LW within a 10 m radius 36% of the time in the upper 3 km rehabilitated reach, and 44% of the time in the lower 4.7 km marginal habitat reach. A greater percentage of LW was utilized in riffles in the upper 3 km reach where 90% of redds were built, while a larger percentage of spawners used LW in riffles in the lower 4.7 km reach. LW‐redd interactions occurred at greater rates than by random chance alone in the lower 4.7 km reach, which implies that LW aids spawning in marginal habitats. River managers and salmonid spawning habitat rehabilitation (SHR) projects should take LW additions into consideration as an important component of river rehabilitation. Copyright © 2010 John Wiley & Sons, Ltd.     

Three-dimensional modeling of sediment resuspension in a large shallow lake

We simulated bottom resuspension events in Lake Erie, using a coupled three-dimensional hydrodynamic and water quality model. Key parameters in the model, including critical bottom shear stress (${\tau }_{\mathit{cr}}$) and resuspension rate (α) were calibrated and validated by comparing the model output to observations. These included total suspended solid (TSS) concentrations in the bottom boundary layer (RMSE = 0.74 $\text{mg}{\text{L}}^{\text{-1}}$) and water column (RMSE = 0.81 $\text{mg}{\text{L}}^{\text{-1}}$), and to time series of acoustic backscatter signal (R² > 0.8) and turbidity (R² ≈ 0.4) from long-term moorings near the lakebed in 2008–09 and 2013. Signals from phytoplankton, in spring and summer, caused discrepancies between modeled TSS and the observed turbidity data. Although common practice, we show that literature-based or field-observed critical shear stress should not be directly applied in large-scale Reynolds-averaged sediment model as this will likely underestimate resuspension. In agreement with the literature, the model reproduced more frequent and intensive surface-wave driven resuspension in the shallow regions (< ~20 m), particularly in the western basin, compared to the deeper central and eastern basins, where bottom stresses induced by mean currents (${\tau }_c$) were comparable with those due to surface waves (${\tau }_w$). However, on the north-shore of the eastern basin, ${\tau }_c$ often predominated over ${\tau }_w$. We simulated thermocline motion, including up- and down-welling events and swashing of the internal Poincaré wave, to contribute to ${\tau }_c$ in the central basin and form nepheloid layers.     

Improvement of LCM model and determination of model parameters at watershed scale for flood events in Hongde Basin of China

Considering the fact that the original two-parameter LCM model can only be used to investigate rainfall losses during the runoff period because the initial abstraction is not included, the LCM model was redefined as a three-parameter model, including the initial abstraction coefficient λ, the initial abstraction I_a, and the rainfall loss coefficient R. The improved LCM model is superior to the original two-parameter model, which only includes r and R, where r is the initial rainfall loss index and can be calculated with λ using the Soil Conservation Service curve number(SCS-CN) method, with r=(1+λ). The trial method was used to determine the parameter values of the improved LCM model at the watershed scale for 15 flood events in the Hongde Basin in China. The results show that larger r values are associated with smaller R values, and the parameter R ranges widely from 0.5 to 2.0. In order to improve the practicability of the LCM model, r = 0.833 with λ= 0.2 is reasonable for simplifying calculation. When the LCM model is applied to arid and semi-arid regions, rainfall without yielding runoff should be deducted from the total rainfall for more accurate estimation of rainfall-runoff.