Evaluation of oceanic and atmospheric trajectory schemes in the
TRACMASS trajectory model v6.0
Abstract. Two different trajectory schemes for oceanic and atmospheric general circulation models are compared in two different experiments. The theories of the two trajectory schemes are presented showing the differential equations they solve and why they are mass conserving. One scheme assumes that the velocity fields are stationary for a limited period of time and solves the trajectory path from a differential equation only as a function of space, i.e. "stepwise stationary". The second scheme uses a continuous linear interpolation of the fields in time and solves the trajectory path from a differential equation as a function of both space and time, i.e. "time-dependent". A special case of the "stepwise-stationary" scheme, when velocities are assumed constant between GCM outputs, is also considered, named "fixed GCM time step". The trajectory schemes are tested "off-line", i.e. using the already integrated and stored velocity fields from a GCM. The first comparison of the schemes uses trajectories calculated using the velocity fields from an eddy-resolving ocean general circulation model in the Agulhas region. The second comparison uses trajectories calculated using the wind fields from an atmospheric reanalysis. The study shows that using the "time-dependent" scheme over the "stepwise-stationary" scheme greatly improves accuracy with only a small increase in computational time. It is also found that with decreasing time steps the "stepwise-stationary" scheme becomes more accurate but at increased computational cost. The "time-dependent" scheme is therefore preferred over the "stepwise-stationary" scheme. However, when averaging over large ensembles of trajectories the two schemes are comparable, as intrinsic variability dominates over numerical errors. The "fixed GCM time step" is found to be less accurate than the "stepwise-stationary" scheme, even when considering averages over large ensembles.