This paper aims to stimulate discussion based on the experiences derived from the QUICS project (Quantifying Uncertainty in Integrated Catchment Studies). First it briefly discusses the current state of knowledge on uncertainties in sub-models of integrated catchment models and the existing frameworks for analysing uncertainty. Furthermore, it compares the relative approaches of both building and calibrating fully integrated models or linking separate sub-models. It also discusses the implications of model linkage on overall uncertainty and how to define an acceptable level of model complexity. This discussion includes, whether we should shift our attention from uncertainties due to linkage, when using linked models, to uncertainties in model structure by necessary simplification or by using more parameters. This discussion attempts to address the question as to whether there is an increase in uncertainty by linking these models or if a compensation effect could take place and that overall uncertainty in key water quality parameters actually decreases. Finally, challenges in the application of uncertainty analysis in integrated catchment water quality modelling, as encountered in this project, are discussed and recommendations for future research areas are highlighted.

Muddy coasts provide ecological habitats, supply food and form a natural coastal defence. Relative sea level rise, changing wave energy and human interventions will increase the pressure on muddy coastal zones. For sustainable coastal management it is key to obtain information on the geomorphology of and historical changes along muddy areas. So far, little is known about the distribution and behaviour of muddy coasts at a global scale. In this study we present a global scale assessment of the occurrence of muddy coasts and rates of coastline change therein. We combine publicly available satellite imagery and coastal geospatial datasets, to train an automated classification method to identify muddy coasts. We find that 14% of the world’s ice-free coastline is muddy, of which 60% is located in the tropics. Furthermore, the majority of the world’s muddy coasts are eroding at rates exceeding 1 m/yr over the last three decades.

Current riverine plastic monitoring best practices mainly consider surface observations, thus neglecting the underlying distribution of plastics in the water column. Bias on plastic budgets estimations hinders advances on modelling and prediction of plastics fate. Here, we experimentally disclose the structure of plastics transport in surface water flows by investigating how thousands of samples of plastics commonly found in fluvial environments travel in turbulent river flows. We show for the first time that surface tension plays a key role in the transport of plastics since its effects can be of the same magnitude as buoyancy and turbulence, therefore holding a part of the dispersed buoyant plastics captive by the water surface. We investigate two types of transport; surfaced plastics (surface tension-turbulence-buoyancy dominated), in contact with the free surface, and suspended plastics (turbulence-buoyancy dominated). We prove that this duality in transport modes is a major source of error in the estimation of plastic budgets, which can be underestimated by 90 % following current, well-established monitoring protocols if sampling is conducted solely in the water surface. Based on our empirical findings, we optimize physics-driven monitoring strategies for plastic fluxes in rivers, thereby achieving over a ten-fold reduction of the bias and uncertainty of riverine plastic pollution estimates.

Receiving water quality simulation in highly urbanised areas requires the integration of several processes occurring at different space-time scales. These integrated catchment models deliver results with a significant uncertainty level associated. Still, uncertainty analysis is seldom applied in practice and the relative contribution of the individual model elements is poorly understood. Often the available methods are applied to relatively small systems or individual sub-systems, due to limitations in organisational and computational resources. Consequently this work presents an uncertainty propagation and decomposition scheme of an integrated water quality modelling study for the evaluation of dissolved oxygen dynamics in a large-scale urbanised river catchment in the Netherlands. Forward propagation of the measured and elicited uncertainty input-parametric distributions was proposed and contrasted with monitoring data series. Prior ranges for river water quality-quantity parameters lead to high uncertainty in dissolved oxygen predictions, thus the need for formal calibration to adapt to the local dynamics is highlighted. After inferring the river process parameters with system measurements of flow and dissolved oxygen, combined sewer overflow pollution loads became the dominant uncertainty source along with rainfall variability. As a result, insights gained in this paper can help in planning and directing further monitoring and modelling efforts in the system. When comparing these modelling results to existing national guidelines it is shown that the commonly used concentration-duration-frequency tables should not be the only metric used to select mitigation alternatives and may need to be adapted in order to cope with uncertainties.

In urban hydrological models, rainfall is the main input and one of the main sources of uncertainty. To reach sufficient spatial coverage and resolution, the integration of several rainfall data sources, including rain gauges and weather radars, is often necessary. The uncertainty associated with rain gauge measurements is dependent on rainfall intensity and on the characteristics of the devices. Common spatial interpolation methods do not account for rain gauge uncertainty variability. Kriging for Uncertain Data (KUD) allows the handling of the uncertainty of each rain gauge independently, modelling space- and time-variant errors. The applications of KUD to rain gauge interpolation and radar-gauge rainfall merging are studied and compared. First, the methodology is studied with synthetic experiments, to evaluate its performance varying rain gauge density, accuracy and rainfall field characteristics. Subsequently, the method is applied to a case study in the Dommel catchment, the Netherlands, where high-quality automatic gauges are complemented by lower-quality tipping-bucket gauges and radar composites. The case study and the synthetic experiments show that considering measurement uncertainty in rain gauge interpolation usually improves rainfall estimations, given a sufficient rain gauge density. Considering measurement uncertainty in radar-gauge merging consistently improved the estimates in the tested cases, thanks to the additional spatial information of radar rainfall data but should still be used cautiously for convective events and low-density rain gauge networks.

A methodology is presented which can be used in the evaluation of parametric uncertainty in urban flooding simulation. Due to the fact that such simulations are time consuming, the following methodology is proposed: (a) simplification of the description of the physical process; (b) derivation of a training data set; (c) development of a data-driven surrogate model; (d) use of a forward uncertainty propagation scheme. The simplification comprises the following steps: (a) unit hydrograph derivation using a 2D hydrodynamic model; (b) calculation of the losses in order to determine the effective rainfall depth; (c) flood event simulation using the principle of the proportionality and superposition. The above methodology was implemented in an urban catchment located in the city of Athens, Greece. The model used for the first step of the simplification was FLOW-R2D, whereas the well-known SWMM software (US Environmental Protection Agency, Washington, DC, USA) was used for the second step of the simplification. For the training data set derivation, an ensemble of 100 Unit Hydrographs was derived with the FLOW-R2D model. The parameters which were modified in order to produce this ensemble were the Manning coefficients in the two friction zones (residential and urban open space areas). The surrogate model used to replicate the unit hydrograph derivation, using the Manning coefficients as an input, was based on the Polynomial Chaos Expansion technique. It was found that, although the uncertainties in the derived results have to be taken into account, the proposed methodology can be a fast and efficient way to cope with dynamic flood simulation in an urban catchment.