Abstract

This article presents the comparison of two algorithms for data assimilation of two dimensional shallow water flows. The first algorithm is based on a linearization of the model equations and a quadratic programming (QP) formulation of the problem. The second algorithm uses Ensemble Kalman Filtering (EnKF) applied to the non-linear two dimensional shallow water equations. The two methods are implemented on a scenario in which boundary conditions and Lagrangian measurements are available. The performance of the methods is evaluated using twin experiments with experimentally measured bathymetry data and boundary conditions from a river located in the Sacramento Delta. The sensitivity of the algorithms to the number of drifters, low or high discharge and time sampling frequency is studied.

Highlights

  • The modeling and monitoring of river hydraulics are increasingly important as they provide drinkable water for populations as well as irrigation for a variety of crops

  • We present two data assimilation methods applied to river flows, namely, a novel algorithm based on Quadratic Programming (QP), and an algorithm using Ensemble Kalman Filtering (EnKF)

  • Two data assimilation methods were applied to a river in the Sacramento Delta and their respective performances were compared

Read more

Summary

Introduction

The modeling and monitoring of river hydraulics are increasingly important as they provide drinkable water for populations as well as irrigation for a variety of crops. The notations for the EnKF algorithm presented in the rest of section 4 are as follows: θn: Vector of state variables, namely the velocity components (u, v), the water height h for each mesh point and the positions of the drifters xD, yD at a time instant tn. Lagrangian drifter positions are embedded in the state vector in discretized form and the new state vector θn at time tn is written as θn = This approach makes it possible to take the Lagrangian nature of observations into account in a way which does not require deriving Eulerian quantities from Lagrangian measurements as was done in the previous section in equations (18) and (19). For the state space model (20)–(21), the EnKF algorithm can be summarized as in [15, 12]: Algorithm

Measurement update
Results
Conclusion
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.