We present a unified approach in which we reformulate various hydrological routing approaches as a cascade of non-linear reservoirs. This covers linear and nonlinear reservoir routing as well as Muskingum-type schemes. Whereas original variable-parameter versions of these schemes, e.g. Muskingum-Cunge, are not mass conservative, the reformulated version guarantees strict mass conservation. An iterative Newton-Raphson scheme integrates the implicit schematization of the reservoir equation in time. The first order sensitivities are derived by the implicit function theorem and the adjoint sensitivity equation at the computational costs of the time integration itself.
The novel framework is applied to an academic test case and a routing network in the upper basin of Main River in Germany. Results show the performance of the routing scheme and verify the mass conservation properties of the approach. Furthermore, we give an outlook at the future use of the model within variational DA and MPC scenarios.