  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 coherent 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.
