Using hysteresis analysis of high-resolution water quality monitoring data, including uncertainty, to infer controls on nutrient and sediment transfer in catchments

A large proportion of nutrients and sediment is mobilised in catchments during storm events. Therefore understanding a catchment's hydrological behaviour during storms and how this acts to mobilise and transport nutrients and sediment to nearby watercourses is extremely important for effective catchment management. The expansion of available in-situ sensors is allowing a wider range of water quality parameters to be monitored and at higher temporal resolution, meaning that the investigation of hydrochemical behaviours during storms is increasingly feasible. Studying the relationship between discharge and water quality parameters in storm events can provide a valuable research tool to infer the likely source areas and flow pathways contributing to nutrient and sediment transport. Therefore, this paper uses 2years of high temporal resolution (15/30min) discharge and water quality (nitrate-N, total phosphorus (TP) and turbidity) data to examine hysteretic behaviour during storm events in two contrasting catchments, in the Hampshire Avon catchment, UK. This paper provides one of the first examples of a study which comprehensively examines storm behaviours for up to 76 storm events and three water quality parameters. It also examines the observational uncertainties using a non-parametric approach. A range of metrics was used, such as loop direction, loop area and a hysteresis index (HI) to characterise and quantify the storm behaviour. With two years of high resolution information it was possible to see how transport mechanisms varied between parameters and through time. This study has also clearly shown the different transport regimes operating between a groundwater dominated chalk catchment versus a surface-water dominated clay catchment. This information, set within an uncertainty framework, means that confidence can be derived that the patterns and relationships thus identified are statistically robust. These insights can thus be used to provide information regarding transport processes and biogeochemical processing within river catchments

Sensitivity of a hydraulic model to changes in channel erosion during extreme flooding

Recent research into flood modelling has primarily concentrated on the simulation of inundation flow without considering the influences of channel morphology. River channels are often represented by a simplified geometry that is implicitly assumed to remain unchanged during flood simulations. However, field evidence demonstrates that significant morphological changes can occur during floods to mobilise the boundary sediments. Despite this, the effect of channel morphology on model results has been largely unexplored. To address this issue, the impact of channel cross-section geometry and channel long-profile variability on flood dynamics is examined using an ensemble of a 1D-2D hydraulic model (LISFLOOD-FP) of the 1:2102 year recurrence interval floods in Cockermouth, UK, within an uncertainty framework. A series of hypothetical scenarios of channel morphology were constructed based on a simple velocity based model of critical entrainment. A Monte-Carlo simulation framework was used to quantify the effects of channel morphology together with variations in the channel and floodplain roughness coefficients, grain size characteristics, and critical shear stress on measures of flood inundation. The results showed that the bed elevation modifications generated by the simplistic equations reflected a good approximation of the observed patterns of spatial erosion despite its overestimation of erosion depths. The effect of uncertainty on channel long-profile variability only affected the local flood dynamics and did not significantly affect the friction sensitivity and flood inundation mapping. The results imply that hydraulic models generally do not need to account for within event morphodynamic changes of the type and magnitude modelled, as these have a negligible impact that is smaller than other uncertainties, e.g. boundary conditions. Instead morphodynamic change needs to happen over a series of events to become large enough to change the hydrodynamics of floods in supply limited gravel-bed rivers like the one used in this research

A high-resolution global flood hazard model

Floods are a natural hazard that affect communities worldwide, but to date the vast majority of flood hazard research and mapping has been undertaken by wealthy developed nations. As populations and economies have grown across the developing world, so too has demand from governments, businesses, and NGOs for modeled flood hazard data in these data-scarce regions. We identify six key challenges faced when developing a flood hazard model that can be applied globally and present a framework methodology that leverages recent cross-disciplinary advances to tackle each challenge. The model produces return period flood hazard maps at ∼90 m resolution for the whole terrestrial land surface between 56°S and 60°N, and results are validated against high-resolution government flood hazard data sets from the UK and Canada. The global model is shown to capture between two thirds and three quarters of the area determined to be at risk in the benchmark data without generating excessive false positive predictions. When aggregated to ∼1 km, mean absolute error in flooded fraction falls to ∼5%. The full complexity global model contains an automatically parameterized subgrid channel network, and comparison to both a simplified 2-D only variant and an independently developed pan-European model shows the explicit inclusion of channels to be a critical contributor to improved model performance. While careful processing of existing global terrain data sets enables reasonable model performance in urban areas, adoption of forthcoming next-generation global terrain data sets will offer the best prospect for a step-change improvement in model performance.</p

Improving estuary models by reducing uncertainties associated with river flows

© 2018 The Authors To mitigate against future changes to estuaries such as water quality, catchment and estuary models can be coupled to simulate the transport of harmful pathogenic viruses, pollutants and nutrients from their terrestrial sources, through the estuary and to the coast. To predict future changes to estuaries, daily mean river flow projections are typically used. We show that this approach cannot resolve higher frequency discharge events that have large impacts to estuarine dilution, contamination and recovery for two contrasting estuaries. We therefore characterise sub-daily scale flow variability and propagate this through an estuary model to provide robust estimates of impacts for the future. River flow data (35-year records at 15-min sampling) were used to characterise variabilities in storm hydrograph shapes and simulate the estuarine response. In particular, we modelled a fast-responding catchment-estuary system (Conwy, UK), where the natural variability in hydrograph shapes generated large variability in estuarine circulation that was not captured when using daily-averaged river forcing. In the extreme, the freshwater plume from a ‘flash’ flood (lasting < 12 h) was underestimated by up to 100% – and the response to nutrient loading was underestimated further still. A model of a slower-responding system (Humber, UK), where hydrographs typically last 2–4 days, showed less variability in estuarine circulation and good approximation with daily-averaged flow forcing. Our result has implications for entire system impact modelling; when we determine future changes to estuaries, some systems will need higher resolution future river flow estimates

Ensemble evaluation of hydrological model hypotheses

It is demonstrated for the first time how model parameter, structural and data uncertainties can be accounted for explicitly and simultaneously within the Generalized Likelihood Uncertainty Estimation (GLUE) methodology. As an example application, 72 variants of a single soil moisture accounting store are tested as simplified hypotheses of runoff generation at six experimental grassland field-scale lysimeters through model rejection and a novel diagnostic scheme. The fields, designed as replicates, exhibit different hydrological behaviors which yield different model performances. For fields with low initial discharge levels at the beginning of events, the conceptual stores considered reach their limit of applicability. Conversely, one of the fields yielding more discharge than the others, but having larger data gaps, allows for greater flexibility in the choice of model structures. As a model learning exercise, the study points to a “leaking” of the fields not evident from previous field experiments. It is discussed how understanding observational uncertainties and incorporating these into model diagnostics can help appreciate the scale of model structural error

Modular Assessment of Rainfall-Runoff Models Toolbox (MARRMoT) v1.2:an open-source, extendable framework providing implementations of 46 conceptual hydrologic models as continuous state-space formulations

This paper presents the Modular Assessment of Rainfall-Runoff Models Toolbox (MARRMoT): A modular open-source toolbox containing documentation and model code based on 46 existing conceptual hydrologic models. The toolbox is developed in MATLAB and works with Octave. MARRMoT models are based solely on traceable published material and model documentation, not on already-existing computer code. Models are implemented following several good practices of model development: The definition of model equations (the mathematical model) is kept separate from the numerical methods used to solve these equations (the numerical model) to generate clean code that is easy to adjust and debug; the implicit Euler time-stepping scheme is provided as the default option to numerically approximate each model's ordinary differential equations in a more robust way than (common) explicit schemes would; threshold equations are smoothed to avoid discontinuities in the model's objective function space; and the model equations are solved simultaneously, avoiding the physically unrealistic sequential solving of fluxes. Generalized parameter ranges are provided to assist with model inter-comparison studies. In addition to this paper and its Supplement, a user manual is provided together with several workflow scripts that show basic example applications of the toolbox. The toolbox and user manual are available from span classCombining double low line"uri"https://github.com/wknoben/MARRMoT/span (last access: 30 May 2019; a hrefCombining double low line"https://doi.org/10.5281/zenodo.3235664"https://doi.org/10.5281/zenodo.3235664). Our main scientific objective in developing this toolbox is to facilitate the inter-comparison of conceptual hydrological model structures which are in widespread use in order to ultimately reduce the uncertainty in model structure selection

The role of bedrock topography on subsurface storm flow

[1] We conducted a detailed study of subsurface flow and water table response coupled with digital terrain analysis (DTA) of surface and subsurface features at the hillslope scale in Panola Mountain Research Watershed (PMRW), Georgia. Subsurface storm flow contributions of macropore and matrix flow in different sections along an artificial trench face were highly variable in terms of timing, peak flow, recession characteristics, and total flow volume. The trench flow characteristics showed linkages with the spatial tensiometer response defining water table development upslope. DTA of the ground surface did not capture the observed spatial patterns of trench flow or tensiometric response. However, bedrock surface topographic indices significantly improved the estimation of spatial variation of flow at the trench. Point-scale tensiometric data were also more highly correlated with the bedrock surface-based indices. These relationships were further assessed for temporal changes throughout a rainstorm. Linkages between the bedrock indices and the trench flow and spatial water table responses improved during the wetter periods of the rainstorm, when the hillslope became more hydrologically connected. Our results clearly demonstrate that in developing a conceptual framework for understanding the mechanisms of runoff generation, local bedrock topography may be highly significant at the hillslope scale in some catchments where the bedrock surface acts as a relatively impermeable boundary