An Asynchronous Solver for Differential Equations Arising from River Basin Models

Source: ICHE Conference Archive - https://mdi-de.baw.de/icheArchiv

Doing Hydrology Backwards—Analytic Solution Connecting Streamflow Oscillations at the Basin Outlet to Average Evaporation on a Hillslope

The concept of doing hydrology backwards, introduced in the literature in the last decade, relies on the possibility to invert the equations relating streamflow fluctuations at the catchment outlet to estimated hydrological forcings throughout the basin. In this work, we use a recently developed set of equations connecting streamflow oscillations at the catchment outlet to baseflow oscillations at the hillslope scale. The hillslope-scale oscillations are then used to infer the pattern of evaporation needed for streamflow oscillations to occur. The inversion is illustrated using two conceptual models of movement of water in the subsurface with different levels of complexity, but still simple enough to demonstrate our approach. Our work is limited to environments where diel oscillations in streamflow are a strong signal in streamflow data. We demonstrate our methodology by applying it to data collected in the Dry Creek Experimental Watershed in Idaho and show that the hydrology backwards principles yield results that are well within the order of magnitude of daily evapotranspiration fluctuations. Our analytic results are generic and they encourage the development of experimental campaigns to validate integrated hydrological models and test implicit parameterization assumptions

Improving Hillslope Link Model Performance from Non-Linear Representation of Natural and Artificially Drained Subsurface Flows

This study evaluates the potential for a newly proposed non-linear subsurface flux equation to improve the performance of the hydrological Hillslope Link Model (HLM). The equation contains parameters that are functionally related to the hillslope steepness and the presence of tile drainage. As a result, the equation provides better representation of hydrograph recession curves, hydrograph timing, and total runoff volume. The authors explore the new parameterization’s potential by comparing a set of diagnostic and prognostic setups in HLM. In the diagnostic approach, they configure 12 different scenarios with spatially uniform parameters over the state of Iowa. In the prognostic case, they use information from topographical maps and known locations of tile drainage to distribute parameter values. To assess performance improvements, they compare simulation results to streamflow observations during a 17-year period (2002–2018) at 140 U.S. Geological Survey (USGS) gauging stations. The operational setup of the HLM model used at the Iowa Flood Center (IFC) serves as a benchmark to quantify the overall improvement of the model. In particular, the new equation provides better representation of recession curves and the total streamflow volumes. However, when comparing the diagnostic and prognostic setups, the authors found discrepancies in the spatial distribution of hillslope scale parameters. The results suggest that more work is required when using maps of physical attributes to parameterize hydrological models. The findings also demonstrate that the diagnostic approach is a useful strategy to evaluate models and assess changes in their formulations

Improving Hillslope Link Model Performance from Non-Linear Representation of Natural and Artificially Drained Subsurface Flows

This study evaluates the potential for a newly proposed non-linear subsurface flux equation to improve the performance of the hydrological Hillslope Link Model (HLM). The equation contains parameters that are functionally related to the hillslope steepness and the presence of tile drainage. As a result, the equation provides better representation of hydrograph recession curves, hydrograph timing, and total runoff volume. The authors explore the new parameterization’s potential by comparing a set of diagnostic and prognostic setups in HLM. In the diagnostic approach, they configure 12 different scenarios with spatially uniform parameters over the state of Iowa. In the prognostic case, they use information from topographical maps and known locations of tile drainage to distribute parameter values. To assess performance improvements, they compare simulation results to streamflow observations during a 17-year period (2002–2018) at 140 U.S. Geological Survey (USGS) gauging stations. The operational setup of the HLM model used at the Iowa Flood Center (IFC) serves as a benchmark to quantify the overall improvement of the model. In particular, the new equation provides better representation of recession curves and the total streamflow volumes. However, when comparing the diagnostic and prognostic setups, the authors found discrepancies in the spatial distribution of hillslope scale parameters. The results suggest that more work is required when using maps of physical attributes to parameterize hydrological models. The findings also demonstrate that the diagnostic approach is a useful strategy to evaluate models and assess changes in their formulations

On the propagation of diel signals in river networks using analytic solutions of flow equations

Several authors have reported diel oscillations in streamflow records and have hypothesized that these oscillations are linked to evapotranspiration cycles in the watershed. The timing of oscillations in rivers, however, lags behind those of temperature and evapotranspiration in hillslopes. Two hypotheses have been put forth to explain the magnitude and timing of diel streamflow oscillations during low-flow conditions. The first suggests that delays between the peaks and troughs of streamflow and daily evapotranspiration are due to processes occurring in the soil as water moves toward the channels in the river network. The second posits that they are due to the propagation of the signal through the channels as water makes its way to the outlet of the basin. In this paper, we design and implement a theoretical model to test these hypotheses. We impose a baseflow signal entering the river network and use a linear transport equation to represent flow along the network. We develop analytic streamflow solutions for the case of uniform velocities in space over all river links. We then use our analytic solution to simulate streamflows along a self-similar river network for different flow velocities. Our results show that the amplitude and time delay of the streamflow solution are heavily influenced by transport in the river network. Moreover, our equations show that the geomorphology and topology of the river network play important roles in determining how amplitude and signal delay are reflected in streamflow signals. Finally, we have tested our theoretical formulation in the Dry Creek Experimental Watershed, where oscillations are clearly observed in streamflow records. We find that our solution produces streamflow values and fluctuations that are similar to those observed in the summer of 2011