Lagrangian descriptors with uncertainty

Lagrangian descriptors provide a global dynamical picture of the geometric structures for arbitrarily time-dependent flows with broad applications. This paper develops a mathematical framework for computing Lagrangian descriptors when uncertainty appears. The uncertainty originates from estimating the underlying flow field as a natural consequence of data assimilation or statistical forecast. It also appears in the resulting Lagrangian trajectories. The uncertainty in the flow field directly affects the path integration of the crucial nonlinear positive scalar function in computing the Lagrangian descriptor, making it fundamentally different from many other diagnostic methods. Despite being highly nonlinear and non-Gaussian, closed analytic formulae are developed to efficiently compute the expectation of such a scalar function due to the uncertain velocity field by exploiting suitable approximations. A rapid and accurate sampling algorithm is then built to assist the forecast of the probability density function (PDF) of the Lagrangian trajectories. Such a PDF provides the weight to combine the Lagrangian descriptors along different paths. Simple but illustrative examples are designed to show the distinguished behavior of using Lagrangian descriptors in revealing the flow field when uncertainty appears. Uncertainty can either completely erode the flow3 structure or barely affect the underlying geometry of the flow field. The method is also applied for eddy identification, indicating that uncertainty has distinct impacts on detecting eddies at different time scales. Finally, when uncertainty is incorporated into the Lagrangian descriptor for inferring the source target, the likelihood criterion provides a very different conclusion from the deterministic methods.