Research Outputs
The influence of pool-riffle morphological features on river mixing
Accurate prediction of pollutant concentrations in a river course is of great importance in environmental management. Mathematical dispersion models are often used to predict the spatial distribution of substances to help achieve these objectives. In practice, these models use a dispersion coefficient as a calibration parameter that is calculated through either expensive field tracer experiments or through empirical equations available in the scientific literature. The latter are based on reach-averaged values obtained from laboratory flumes or simple river reaches, which often show great variability when applied to natural streams. These equations cannot directly account for mixing that relates specifically to spatial fluctuations of channel geometry and complex bed morphology. This study isolated the influence of mixing related to bed morphology and presented a means of calculating a predictive longitudinal mixing equation that directly accounted for pool-riffle sequences. As an example, a predictive equation was developed by means of a three-dimensional numerical model based on synthetically generated pool-riffle bathymetries. The predictive equation was validated with numerical experiments and field tracer studies. The resulting equation was shown to more accurately represent mixing across complex morphology than those relations selected from the literature.
Runoff volume and sediment yield from an endorheic watershed generated by rare rainfall events in the Atacama Desert
Surface runoff, channel activation and sediment transport processes in desert environments have been convincingly shown to strongly depend on the duration and intensity of local and convective rainfall events. Among these environments, the Atacama, considered the driest desert on Earth, is situated in a remote and rugged area, where documentation of historical and recent hydro-sedimentary processes is rare. We characterize the hydrological processes in an endorheic watershed of the Atacama's Altiplano Desert, where the occurrence of flash floods was evaluated on event-based signatures of water and sediment in a small ephemeral playa. Twelve pits were dug in the playa, with five identified event sediment couplets, each corresponding to computed flood volumes that gave rise to transport and depositional events occurring between 1978 and 2019. Detailed topography allowed reconstruction of a 3D terrain model, from which we estimated a 11.3 t/yr/km2 local historic sediment yield. The timing of the older identified events did not match local rainfall records, and proved to be uncorrelated with occurrences of El Niño–Southern Oscillation (ENSO). This suggests a high spatial patchiness of rainfall events in terms of depth and intensity, and implies that small playa records are not necessarily always helpful in reconstructing the regional climatic history of the recent past. The sediment concentration and volume of the reconstructed hypercontracted events are not well corrected with the magnitude of the rainfall event, suggesting the important role played by variable sediment availability and connectivity at basin scale. This spatiotemporal variability plays a major role in understanding the present and historic hydro-sedimentary processes in the Atacama's altiplano.
Identifying dynamic equilibrium of an undeveloped alluvial stream by extremal hypotheses
The indeterminate channel problem arises from uncertainty in finding a closure relation for alluvial channels created by self-organizing erosional and depositional processes. Extremal hypotheses have been proposed as one potential approach to closing the system of governing equations for alluvial channels. Many different extremal hypotheses have been presented, but no substantive evidence has been developed to select which hypothesis may be most appropriate for natural alluvial river systems. This paper evaluates the ability of ten extremal hypotheses to identify dynamic equilibrium across a geomorphic gradient in the remote and undeveloped mid-latitude watershed of Rio Murta, Chile. This study (a) introduces extremal hypotheses, (b) describes the field site and geomorphic conditions, and (c) examines which extremal hypotheses are supported by the field data in identifying the evolutionary trend toward dynamic-equilibrium. The extremal hypotheses that identified dynamic equilibrium within the geomorphic gradient in the field are: (1) minimum kinetic energy, (2) minimum specific stream power, (3) maximum friction factor, and (4) maximum total friction factor, which collectively support minimizing kinetic energy of the system.
Testing the effective-discharge paradigm in gravel-bed river restoration
Restoration of an alluvial wet-meadow system was conducted in the late 1990s to reestablish hydraulic interactions between the river, floodplain, and groundwater to support aquatic–riparian ecosystem function. A single-discharge approach sized the bankfull channel dimensions to the effective discharge (Qe) and three degrees of channel widening relative to the Qe design were explored to identify which design attained dynamic equilibrium in the shortest time. The three experimental channel designs were implemented with bankfull widths of 96%, 157%, and 191%, respectively, of the Qe geometry. Response trajectories were documented for channel dimensions, sediment mobility, channel morphology, floodplain connectivity, and riparian vegetation for the three channel designs, and the efficacy of a single-discharge approach for restoration was examined. Analysis of 20 years of monitoring data and hydraulic modeling revealed that each design responded differently to the imposed initial channel conditions and evolved at substantially different rates. The design with bankfull dimensions most closely approximating Qe reached dynamic equilibrium within four years of restoration, whereas the moderately over-widened channel (57% larger) exhibited slower responses toward dynamic equilibrium for some metrics and did not fully attain the Qe design dimensions within the monitoring period. The extremely-over-widened channel (91% larger) mainly induced slow rates of bed deposition that are projected to take nearly 300 years for the bankfull dimensions to narrow to the Qe width. All reaches had low bed mobility (bankfull Shields stress < 0.03) 14 years after restoration, demonstrating the challenge of reducing the drivers of channel widening, while maintaining sufficient competence for bedload transport and a sustained supply of coarse bed material for salmonid habitat. Restoration that sizes the channel to Qe can provide rapid dynamic equilibrium, but is a first-order simplification of 1) channel dynamics and 2) the range of flows needed for restoring physical and biological processes in wet-meadow systems.
Evaluation of extremal hypotheses in an undeveloped alluvial river
Recent work in the undeveloped Rio Murta Basin, located in Chilean Patagonia, identified an evolutionary trend in the fluvial system as it progresses toward and away from dynamic equilibrium. A location-for-time-substitution model employed over the longitudinal extent of a 16 km study site assessed the performance of extremal hypotheses in identifying dynamic equilibrium conditions. Numerous extremal hypotheses were successful in identifying the spatial trend, but no means were available to discern differences between them. Thus, this study uses field measurements within the evolutionary trend to propose a new metric for evaluating extremal hypotheses. A thorough review and synthesis of the extremal approach are additionally presented. The new method compares theoretical predictions against field-measured values to determine which extremal hypothesis is most effective in identifying the condition of dynamic equilibrium in a gravel-bed river. Channel width and depth are identified as the dependent stream variables that uniquely differentiate most extremal hypotheses from one another. The results indicate that extremal hypotheses based on energy metrics of the flow are most successful, with the strongest support for minimum kinetic energy and minimum specific stream power.
Influence of dam breach parameter statistical definition on resulting rupture maximum discharge
However rare, dam breach occurrences are recently reported and associated with significant damage to life and property. The rupture of the structural dam wall generates severe flow rates that exceed spillway capacity consequently generating unprecedented flooding scenarios. The present research aims to assess the influence of the dam breach statistical configuration on the most relevant parameters to predict the rupture maximum discharge (RMD). McBreach© software was used to provide the necessary inputs for the operation of the HEC-RAS dam breach module. McBreach© automates the process of batch mode simulations providing a Monte Carlo approach to characterize the breach parameters stochastically. Thus, a sensitivity analysis was performed to identify the most influential breach parameters, followed by an uncertainty assessment regarding their statistical definition of the resultant RMD. Analysis showed that the overtopping failure mode discharges are most sensitive to the breach formation time (tf) parameter, followed by the final height breach (Inv) and the final width of the breach (B), which combined are responsible for 85% of the rupture’s maximum discharge. Further results indicated highly variable RMD magnitudes (up to 300%) depending on the breach parameter’s statistical definition (i.e., probability density function and associated statistical parameters). The latter significantly impacts the estimated flood risk associated with the breach, the flood zone delimitation, preparation of emergency action plans (EAP) and scaling of future dam projects. Consequently, there is a plausible need for additional investigations to reduce this uncertainty and, therefore, the risk associated with it.
Urban stream syndrome: Quantifying topographic variation along an urban-rural gradient
Topographic variation within fluvial systems is essential for providing a mosaic of physical habitats and supporting the dynamic hydraulic, geochemical, and biological processes that determine both aquatic and riparian ecosystem function. In highly-modified rivers through both urban and rural settings, the physical heterogeneity of alluvial channels has been diminished by anthropogenic activities. As riparian areas are increasingly under pressure from agricultural and urban development, identifying the geomorphic controls on physical heterogeneity through these environments is critical. In this study, we use the bed coefficient of variation (CV) extracted from a high-resolution bathymetric LiDAR survey as a dimensionless metric for topographic variation and physical heterogeneity over 100 km of the Boise River corridor that spans an urban-rural gradient. Our CV results for both the streambed and channel demonstrate that the average topographic variation of reaches in urban areas is 22–25% lower than reaches located in rural areas along the same river. While these results initially support the application of the urban stream syndrome hypothesis, CV values had similar magnitudes in both urban and rural reaches suggesting there is a dominant control on topographic variation that was not directly related to urban land use. Analysis of CV values relative to normalized levee width indicates that the causative driver of morphologic simplification in the channel was lateral constraints from levees. In the Boise River, topographic variation increased linearly with normalized levee widths that ranged between 50% and >300% of the average channel width. Further, topographic variation was maximized in reaches where flow expansion during high discharge inundated between 1 and 2 times the average channel width (approximately 65–70% of the available floodplain). Our simple and objective watershed-scale approach leverages highresolution topography data to identify reaches of high physical heterogeneity for river conservation, as well as help guide environmental flow releases in managed rivers.
Effects on fluvial geomorphology and vegetation cover following hydroelectric power plant operation: A case study in the Maule river (Chile)
The installation of hydroelectric plants has generated multiple environmental impacts on the world’s river systems. In central Chile, the impacts of hydroelectric reservoir operation have been documented in ecological and hydrologic regime terms. This investigation assesses the changes in channel morphology, vegetation distribution, and flows in the middle section of the Maule River during the period following the start-up of a hydroelectric plant. Changes in fluvial morphology (active area) and land cover are quantified using LANDSAT images, contrasted with a vegetation sampling and flow analysis. The results show a 12% decrease in active areas of the river, indicating a loss of geomorphological diversity. Within the active channel, there was a gradual increase in plant-covered surface area, which reached 159% between 1989 and 2018, mainly due to reductions in water (−61%), active bar (−35%), and bare soil surface areas (−29%). The changes were evident ten years after plant operations began and intensified during the period known as the megadrought in central Chile (2008–2018). The flow magnitudes present a decrease for exceedance probabilities (P) below 85% in the period after 1985, with a slight increase recorded for low flows (P > 85%). In the segments with superior stabilization, invasive species such as Acacia dealbata (silver wattle) predominated, which are specialists at taking advantage of disturbances to settle and stabilize active areas, narrowing the possibilities for morphological change.
Evaluation of products 3B42 v7 and 3IMERG for the hydroclimatic regions of Chile
Data provided by products 3B42 V7 (TRMM) and its successor 3IMERG (GPM) are compared with discrete rainfall information throughout the Chilean territory covering four macro hydroclimatic zones. Precipitation data was obtained from weather stations available on a daily basis from the 1930s to present. A total of 143 stations were chosen and rainfall estimates performed for years 2014 through 2018. Applying the same metrics we showed how GPM performances improve as the temporal aggregation increases. Several drawbacks were detected in the coastal areas, which were characterized by lower accuracy performances than internal areas. However, the 3IMERG product could be a strong source of data to study the impact that climatic disturbances have on the hydrologic cycle in the Central and South zones of Chile. Additionally, its offers a fundamental source of data for remote zones or areas where access is complicated to install weather stations.
The shell phenotypic variability of the keyhole limpet Fissurella latimarginata: Insights from an experimental approach using a water flow flume
Hydrodynamics are a major environmental factor on intertidal rocky shores. Morphological responses to this factor are expected to strongly influence spatial distribution of species across environmental gradients. We here analysed the shell phenotypic variability of the limpet Fissurella latimarginata using geometric morphometric analysis. The limpets were obtained from a sheltered intertidal coastal area and a wave-exposed environment. To determine whether the shell shape variation of the intertidal molluscs is linked to their resistance to differential intertidal wave exposure, mesocosm studies were developed in a hydraulic flume to explore the effects of hydrodynamic forces on this limpet species. A unidirectional current was used to test the impacts of step-by-step increased current flow velocities for each limpet. The phenotypic variability observed in the populations of F. latimarginata was associated with habitats characterized by contrasting wave exposure. Limpets from exposed environments showed a flattened, round to laterally wider and posteriorly narrower shell shape, larger foot and higher full limpet height, and were dislodged at higher velocities. A more laterally compressed and peaked shape was found in limpets from sheltered areas and these showed a lower resistance to wave action by dislodging at lower velocities.