Improving the stability of temporal statistics in transition path theory with sparse data

Ulam's method is a popular discretization scheme for stochastic operators that involves the construction of a transition probability matrix controlling a Markov chain on a set of cells covering some domain. We consider an application to satellite-tracked undrogued surface-ocean drifting buoy trajectories obtained from the National Oceanic and Atmospheric Administration Global Drifter Program dataset. Motivated by the motion of Sargassum in the tropical Atlantic, we apply Transition Path Theory (TPT) to drifters originating off the west coast of Africa to the Gulf of Mexico. We find that the most common case of a regular covering by equal longitude-latitude side cells can lead to a large instability in the computed transition times as a function of the number of cells used. We propose a different covering based on a clustering of the trajectory data that is stable against the number of cells in the covering. We also propose a generalization of the standard transition time statistic of TPT that can be used to construct a partition of the domain of interest into weakly dynamically connected regions.