Exploring the limits of ensemble forecasting via solutions of the Liouville equation for realistic geophysical models

The atmosphere is an exemplary case of uncertain system. The state of such systems is described by means of probability density functions which encompass uncertainty information. In this regard, the Liouville equation is the theoretical framework to predict the evolution of the state of uncertain systems. This study analyses the morphological characteristics of the time evolution of probability density functions for some low complexity geophysical systems by solving the Liouville equation in order to obtain tractable solutions which are otherwise unfeasible with currently available computational resources. The current and usual modest approach to overcome these obstacles and estimate the probability density function of the system in realistic weather and climate applications is the use of a discrete and small number of samples of the state of the system, evolved individually in a deterministic, perhaps sometimes stochastic, way. We investigate particular solutions of the shallow water equations and the barotropic model that allow to apply the Liouville formalism to explore its topological characteristics and interpret them in terms of the ensemble prediction system approach. We provide quantitative evidences of the high variability that solutions to Liouville equation may present, challenging currently accepted uses and interpretations of ensemble forecasts.