Modelling the barotropic sea level in the Mediterranean Sea using data assimilation

This paper analyses the variability of the sea level barotropic components in the Mediterranean Sea and their reproduction using a hydrodynamic model with and without data assimilation. The impact of data assimilation is considered both in reanalysis and short-forecast simulations. We used a two-dimensional finite element model paired with an ensemble Kalman filter, which assimilated hourly sea level data from 50 stations in the Mediterranean basin. The results brought about a significant improvement given by data assimilation in the reanalysis of the astronomical tide, the surge, and the barotropic total sea level, even in coastal areas and far from the assimilated stations (e.g. the southeastern Mediterranean Sea). As with the reanalysis simulations, the forecast simulations, which start from analysis states, improve, especially on the first day (37 % average error reduction) and when seiche oscillations are triggered. Since seiches are free barotropic oscillations that depend only on the initial state, their reproduction improves very effectively with data assimilation. Finally, we estimate the periods and the energy of these oscillations by means of spectral analysis, both in the Adriatic Sea, where they have been extensively studied, and in the Mediterranean Sea, where the present documentation is scarce. While the periods are well reproduced by the model even without data assimilation, their energy shows a good improvement when using it.