We construct a climatology of Lagrangian coherent structures (LCSs)—the concealed skeleton that shapes transport—with a twelve-year-long data-assimilative simulation of the sea-surface circulation in the Gulf of Mexico (GoM). Computed as time-mean Cauchy–Green strain tensorlines of the climatological velocity, the climatological LCSs (cLCSs) unveil recurrent Lagrangian circulation patterns. The cLCSs strongly constrain the ensemble-mean Lagrangian circulation of the instantaneous model velocity, showing that a climatological velocity can preserve meaningful transport information. The quasi-steady transport patterns revealed by the cLCSs agree well with aspects of the GoM circulation described in several previous observational and numerical studies. For example, the cLCSs identify regions of persistent isolation, and suggest that coastal regions previously identified as high-risk for pollution impact are regions of maximal attraction. We also show that cLCSs are remarkably accurate at identifying transport patterns observed during the Deepwater Horizon and Ixtoc oil spills, and during the Grand LAgrangian Deployment (GLAD) experiment. Thus it is shown that computing cLCSs is an efficient and meaningful way of synthesizing vast amounts of Lagrangian information. The cLCS method confirms previous GoM studies, and contributes to our understanding by revealing the persistent nature of the dynamics and kinematics treated therein.

Persistent Lagrangian transport patterns at the ocean surface are revealed from Lagrangian Coherent Structures (LCSs) computed from daily climatological surface current velocities in the northwestern Gulf of Mexico (NWGoM). The velocities are produced by a submesoscale permitting regional ocean model of the Gulf of Mexico. The significance of the climatological LCSs (cLCSs) is supported with ensemble-mean drifter density evolutions from simulated and historical satellite-tracked drifter trajectories. A persistent attracting barrier between the NWGoM shelf and the deep ocean is effectively identified by a hook-like pattern associated with groups of overall strongly attracting cLCSs that extend along the shelf break.Localized reductions in the attraction rate along these overall strongly attracting cLCSs proximal to cross-shore oriented cLCSs identify a pathway for potential transport across the shelf break.Groups of overall weakly-attracting cLCSs are not seen to strongly constrain material transport.Tracers originating over the shelf tend to be trapped there by the hook-like pattern as they spread cyclonically. Tracers originating beyond the shelf tend to be initially attracted to the hook-like pattern as they spread anti-cyclonically and eventually over the deep ocean. The findings have important implications for the mitigation of contaminant accidents such as oil spills.

The Department of Energy’s (DOE’s) National Energy Technology Laboratory’s (NETL’s) Blowout and Spill Occurrence Model (BLOSOM), and the National Oceanic and Atmospheric Administration’s (NOAA’s) General NOAA Operational Modeling Environment (GNOME) are compared. Increasingly complex simulations are used to assess similarities and differences between the two models’ components. The simulations presented here are forced by ocean currents from a Finite Volume Community Ocean Model (FVCOM) implementation that has excellent skill in representing tidal motion, and with observed wind data that compensates for a coarse vertical ocean model resolution. The comprehensive comparison between GNOME and BLOSOM presented here, should aid modelers in interpreting their results. Beyond many similarities, aspects where both models are distinct are highlighted. Some suggestions for improvement are included, e.g., the inclusion of temporal interpolation of the forcing fields (BLOSOM) or the inclusion of a deflection angle option when parameterizing wind-driven processes (GNOME). Overall, GNOME and BLOSOM perform similarly, and are found to be complementary oil spill models. This paper also sheds light on what drove the historical Point Wells spill, and serves the additional purpose of being a learning resource for those interested in oil spill modeling. The increasingly complex approach used for the comparison is also used, in parallel, to illustrate the approach an oil spill modeler would typically follow when trying to hindcast or forecast an oil spill, including detailed technical information on basic aspects, like choosing a computational time step. We discuss our successful hindcast of the 2003 Point Wells oil spill that, to our knowledge, had remained unexplained. The oil spill models’ solutions are compared to the historical Point Wells’ oil trajectory, in time and space, as determined from overflight information. Our hindcast broadly replicates the correct locations at the correct times, using accurate tide and wind forcing. While the choice of wind coefficient we use is unconventional, a simplified analytic model supported by observations, suggests that it is justified under this study’s circumstances. We highlight some of the key oceanographic findings as they may relate to other oil spills, and to the regional oceanography of the Salish Sea, including recommendations for future studies.

The Brazil Current (BC) is a weak western boundary current flowing along the Southwestern Atlantic Ocean. It is frequently described as a flow with intense mesoscale activity and relatively low volume transport between 5.0 and 10.0 Sv. We use a 13-year eddy-resolving primitive-equation simulation to show that the presence of persistent meanders and eddies leads to characteristic quasi-steady Lagrangian transport patterns, aptly extracted through climatological Lagrangian Coherent Structures (cLCS). The cLCS position the surface expression of the BC core along the 2000-m isobath, in excellent visual agreement with high resolution satellite sea-surface temperature and the model Eulerian mean velocity. The cLCS deformation pattern also responds to zonally persistent cross-shelf SSH transition from positive (high) values near coastline to low (negative) values between 200- and 2000-m and back to positive (high) offshore from the 2000-m isobath. Zonally-paired cyclonic and anticyclonic structures are embedded in this transition, also causing the cLCS to deform into chevrons. An efficient transport barrier is identified close to the 200-m isobath confirmed by limited inshore movement of drogued drifters and accurately indicated by an along slope maxima of climatological strength of attraction. We also show that the persistent cyclonic and anticyclonic structures may induce localized cross-shelf transport. Regions of low climatological strength of attraction coincide with large shelves and with stagnant synthetic trajectories. We also show that cLCS accurately depict trajectories initiated at the location of Chevron’s spill (November 2011) as compared to synthetic and satellite-tracked trajectories, and the outline of the oil from that accident. There is also an agreement between the large-scale oil slicks reaching the Brazilian beaches (from August 2019 to February 2020) and the strength of climatological attraction at the coast. Our work also clarifies the influence of persistent mesoscale structures on the regional circulation. The identification and quantitative description of climatological Lagrangian coherent structures is expected to improve the effectiveness of future emergency response to oil spills, contingency planning, rescue operations, larval and fish connectivity assessment, drifter launch strategies, waste pollutant and marine debris dispersion and destination.