In this work, we analyze Lagrangian material coherence in dynamic open domains. We derive and apply new theory and efficient schemes to extract material sets in dynamic flow fields that remain the most or the least coherent throughout the time interval of interest, with special attention to realistic ocean domains that have multiple time-dependent inlets and/or outlets. First, the partial differential equation (PDE)-based method of composition is extended to efficiently compute flow maps in open domains, evolving a dynamic mask field without compounding numerical errors. This permits the use of existing grid based PDE solvers to compute flow maps on their dynamic non-regular domain. Inherent parallelization capabilities with accuracy as trajectory-based schemes but importantly with also an optimal grid-based resolution make this method very attractive. Second, we derive a novel approach to compute material sets in dynamic fluid flows that undergo minimal stress throughout the considered time interval. The level sets of the proposed metric, called the ‘extended polar distance’, yield material subdomains that remain rigid (i.e. only undergo translation and rotation) throughout the time interval of interest up to a certain tolerance. This metric and the corresponding persistently coherent sets and incoherent sets are computed using the PDE-based flow map computation. We further relate the extended polar distance and the diffusion barrier strength metric and show that the extended polar distance rigorously cumulates the tendency of a material subdomain to be prone to diffusion and the average strain it undergoes. We utilize the new theory and numerical methods to analyze Lagrangian coherence in analytical and realistic scenarios – an analytical unsteady double gyre flow and a realistic simulation in the Southern Pacific Ocean. The former helps us better understand the proposed theory in practice, and highlights the evolution of coherent, persistently coherent, and incoherent sets. In the latter Southern Pacific Ocean application, we find that the surface regions around Palau island are highly incoherent due to the steep topography and complex interactive dynamics. However, we also find a rigid set advected by the larger-scale currents around the Island, retrieving its shape at the end, as well as a persistently rigid set that approximately maintains its shape throughout the time interval, maximally resisting advective stretching and diffusive transport.

This paper focuses on the extractions of Lagrangian Coherent Sets from realistic velocity fields obtained from ocean data and simulations, each of which can be highly resolved and non volume-preserving. We introduce two novel methods for computing two formulations of such sets. First, we propose a new “diffeomorphism-based” criterion to extract “rigid sets”, defined as sets over which the flow map acts approximately as a rigid transformation. Second, we develop a matrix-free methodology that provides a simple and efficient framework to compute “coherent sets” with operator methods. Both new methods and their resulting rigid sets and coherent sets are illustrated and compared using three numerically simulated flow examples, including a high-resolution realistic, submesoscale to large-scale dynamic ocean current field in the Palau Island region of the western Pacific Ocean.

Vertical transport in the ocean plays a critical role in the exchange of freshwater, heat, nutrients, and other biogeochemical tracers. While there are situations where vertical fluxes are important, studying the vertical transport and displacement of material requires analysis over a finite interval of time. One such example is the subduction of fluid from the mixed layer into the pycnocline, which is known to occur near density fronts. Divergence has been used to estimate vertical velocities indicating that surface measurements, where observational data is most widely available, can be used to locate these vertical transport conduits. We evaluate the correlation between surface signatures derived from Eulerian (horizontal divergence, density gradient, and vertical velocity) and Lagrangian (dilation rate and finite time Lyapunov exponent) metrics and vertical displacement conduits. Two submesoscale resolving models of density fronts and a data-assimilative model of the western Mediterranean were analyzed. The Lagrangian surface signatures locate significantly more of the strongest displacement features and the difference in the expected displacements relative to Eulerian ones increases with the length of the time interval considered. Ensemble analysis of forecasts from the Mediterranean model demonstrates that the Lagrangian surface signatures can be used to identify regions of strongest downward vertical displacement even without knowledge of the true ocean state.

The Red Sea Project (RSP) is based on a coastal lagoon with over 90 pristine islands. The project intends to transform the Red Sea coast into a world-class tourist destination. To better understand the regional dynamics and water exchange scenarios in the lagoon, a high-resolution numerical model is implemented. The general and tidal circulation dynamics are then investigated with a particular focus on the response of the lagoon to strong wind jets. Significant variations in winter and summer circulation patterns are identified. The tidal amplitude inside the lagoon is greater than that outside, with strong tidal currents passing over its surrounding coral reef banks. The lagoon rapidly responds to the strong easterly wind jets that occur mainly in winter; it develops a reverse flow at greater depths, and the coastal water elevation is instantly affected. Lagrangian particle simulations are conducted to study the residence time of water in the lagoon. The results suggest that water renewal is slow in winter. Analysis of the Lagrangian coherent structures (LCS) reveals that water renewal is largely linked to the circulation patterns in the lagoon. In winter, the water becomes restricted in the central lagoon with only moderate exchange, whereas in summer, more circulation is observed with a higher degree of interaction between the central lagoon and external water. The results of LCS also highlight the tidal contribution to stirring and mixing while identifying the hotspots of the phenomenon. Our analysis demonstrates an effective approach for studying regional water mixing and connectivity, which could support coastal management in data-limited regions.

Deep-sea polymetallic nodule mining research activity has substantially increased in recent years, but the expected level of environmental impact is still being established. One environmental concern is the discharge of a sediment plume into the midwater column. We performed a dedicated field study using sediment from the Clarion Clipperton Fracture Zone. The plume was monitored and tracked using both established and novel instrumentation, including acoustic and turbulence measurements. Our field studies reveal that modeling can reliably predict the properties of a midwater plume in the vicinity of the discharge and that sediment aggregation effects are not significant. The plume model is used to drive a numerical simulation of a commercial-scale operation in the Clarion Clipperton Fracture Zone. Key takeaways are that the scale of impact of the plume is notably influenced by the values of environmentally acceptable threshold levels, the quantity of discharged sediment, and the turbulent diffusivity in the Clarion Clipperton Fracture Zone.

Every year, hundreds of people die at sea because of vessel and airplane accidents. A key challenge in reducing the number of these fatalities is to make Search and Rescue (SAR) algorithms more efficient. Here, we address this challenge by uncovering hidden TRansient Attracting Profiles (TRAPs) in ocean-surface velocity data. Computable from a single velocity-field snapshot, TRAPs act as short-term attractors for all floating objects. In three different ocean field experiments, we show that TRAPs computed from measured as well as modeled velocities attract deployed drifters and manikins emulating people fallen in the water. TRAPs, which remain hidden to prior flow diagnostics, thus provide critical information for hazard responses, such as SAR and oil spill containment, and hence have the potential to save lives and limit environmental disasters.

To cleanup marine plastics, accurate modeling is needed. We outline and illustrate a new partial-differential-equation methodology for characterizing and modeling plastic transports in time and space (4D), showcasing results for Massachusetts Bay. We couple our primitive equation model for ocean dynamics with our composition based advection for Lagrangian transport. We show that the ocean physics predictions have skill by comparison with synoptic data. We predict the fate of plastics originating from four sources: rivers, beach and nearshore, local Bay, and remote offshore. We analyze the transport patterns and the regions where plastics accumulate, comparing results with and without plastic settling. Simulations agree with existing debris and plastics data. They also show new results: (i) Currents set-up by wind events strongly affect floating plastics. Winds can for example prevent Merrimack outflows reaching the Bay; (ii) There is significant chaotic stirring between nearshore and offshore floating plastics as explained by ridges of Lagrangian Coherent Structures (LCSs); (iii) With 4D plastic motions and settling, plastics from the Merrimack and nearshore regions can settle to the seabed before offshore advection; (iv) Internal waves and tides can bring plastics downward and out of main currents, leading to settling to the deep bottom. (v) Attractive LCSs ridges are frequent in the northern Cape Cod Bay, west of the South Shore, and southern Stellwagen Bank. They lead to plastic accumulation and sinking along thin subduction zones.

To understand the dynamics and health of marine ecosystems such as lagoons and coral reefs as well as to understand the impact of human activities on these systems, it is imperative to predict the residence times of water masses and connectivity between ocean domains. In the present work, we consider the pristine lagoons and coral reefs of the Red Sea as an example of such sensitive ecosystems, with a large number of marine species, many of which are unique to the region. To study the residence times and connectivity patterns, we make use of recent advances in dynamic three-dimensional Lagrangian analyses using partial differential equations. Specifically, we extend and apply our novel efficient flow map composition scheme to predict the time needed for any particular water parcel to leave the domain of interest (i.e., a lagoon) as well as the time for any particular water parcel to enter that domain. These spatiotemporal residence time fields along with four-dimensional Lagrangian metrics such as finite time Lyapunov exponent (FTLE) fields provide a quantitative description of the Lagrangian pathways and connectivity patterns of lagoons in the Red Sea.

We propose a novel numerical methodology to compute the advective transport and diffusion-reaction of tracer quantities. The tracer advection occurs through flow map composition and is super-accurate, yielding numerical solutions almost devoid of compounding numerical errors, while allowing for direct parallelization in the temporal direction. It is computed by implicitly solving the characteristic evolution through a modified transport partial differential equation and domain decomposition in the temporal direction, followed by composition with the known initial condition. This advection scheme allows a rigorous computation of the spatial and temporal error bounds, yields an accuracy comparable to that of Lagrangian methods, and maintains the advantages of Eulerian schemes. We further show that there exists an optimal value of the composition timestep that yields the minimum total numerical error in the computations, and derive the expression for this value. We develop schemes for the addition of tracer diffusion, reaction, and source terms, and for the implementation of boundary conditions. Finally, the methodology is applied in three flow examples, namely an analytical reversible swirl flow, an idealized flow exiting a strait undergoing sudden expansion, and a realistic ocean flow in the Bismarck sea. New benchmark problems for advection-diffusion-reaction schemes are developed and used to compare and contrast results with those of classic schemes. The results highlight the theoretical properties of the methodology as well as its efficiency, super-accuracy with minimal numerical errors, and applicability in realistic simulations.

Three-dimensional (3D) underwater sound field computations have been used for a few decades to understand sound propagation effects above sloped seabeds and in areas with strong 3D temperature and salinity variations. For an approximate simulation of effects in nature, the necessary 3D sound-speed field can be made from snapshots of temperature and salinity from an operational data-driven regional ocean model. However, these models invariably have resolution constraints and physics approximations that exclude features that can have strong effects on acoustics, example features being strong submesoscale fronts and nonhydrostatic nonlinear internal waves (NNIWs). Here, work to predict NNIW fields to improve 3D acoustic forecasts using an NNIW model nested in a tide-inclusive data-assimilating regional model is reported. The work was initiated under the Integrated Ocean Dynamics and Acoustics project. The project investigated ocean dynamical processes that affect important details of sound-propagation, with a focus on those with strong intermittency (high kurtosis) that are challenging to predict deterministically. Strong internal tides and NNIW are two such phenomena, with the former being precursors to NNIW, often feeding energy to them. Successful aspects of the modeling are reported along with weaknesses and unresolved issues identified in the course of the work.

Recent advances in probabilistic forecasting of regional ocean dynamics, and stochastic optimal path planning with massive ensembles motivate principled analysis of their large datasets. Specifically, stochastic time-optimal path planning in strong, dynamic and uncertain ocean flows produces a massive dataset of the stochastic distribution of exact timeoptimal trajectories. To synthesize such big data and draw insights, we apply machine learning and data mining algorithms. Particularly, clustering of the time-optimal trajectories is important to describe their PDFs, identify representative paths, and compute and optimize risk of following these paths. In the present paper, we explore the use of hierarchical clustering algorithms along with a dissimilarity matrix computed from the pairwise discrete Frechet distance between all the optimal trajectories. We apply the algorithms to two datasets of massive ensembles of vehicle trajectories in a stochastic flow past a circular island and stochastic wind driven double gyre flow. These paths are computed by solving our dynamically orthogonal level set equations. Hierarchical clustering is applied to the two datasets, and results are qualitatively and quantitatively analyzed.

With advances in engineering and technology, mining the deep sea for untapped rare metal resources from the bottom of the ocean has recently become economically viable. However, extracting these metal ores from the seabed creates plumes of fine particles that are deposited at various depths within the ocean, and these may be extremely harmful to the marine ecosystems and its components. Thus, for sustainable management, it is of utmost importance to carefully monitor and predict the impact of such harmful activities including plume dispersion on the marine environment. To forecast the plume dispersion in real-time, data-driven ocean modeling has to be coupled with accurate, efficient, and rigorous sediment plume transport computations. The goal of the present paper is to demonstrate the real-time applications of our coupled 3D-and-time data-driven ocean modeling and plume transport forecasting system. Here, the region of focus is the southern California bight, where the PLUMEX 2018 deep sea mining real-time sea experiment was recently conducted (23 Feb – 5 Mar, 2018). Specifically, we demonstrate the improved capabilities of the multiscale MSEAS primitive equation ocean modeling system to capture the complex oceanic phenomenon in the region of interest, the application of the novel method of composition to efficiently and accurately compute the transport of sediment plumes in 3D+1 domains, and the portability of our software and prediction system to different operational regions and its potential in estimating the environmental impacts of deep sea mining activities, ultimately aiding sustainable management and science-based regulations.

Quantifying the uncertainty of Lagrangian motion can be performed by solving a large number of ordinary differential equations with random velocities, or equivalently a stochastic transport partial differential equation (PDE) for the ensemble of flow-maps. The Dynamically Orthogonal (DO) decomposition is applied as an efficient dynamical model order reduction to solve for such stochastic advection and Lagrangian transport. Its interpretation as the method that applies instantaneously the truncated SVD on the matrix discretization of the original stochastic PDE is used to obtain new numerical schemes. Fully linear, explicit central advection schemes stabilized with numerical filters are selected to ensure efficiency, accuracy, stability, and direct consistency between the original deterministic and stochastic DO advections and flow-maps. Various strategies are presented for selecting a time-stepping that accounts for the curvature of the fixed rank manifold and the error related to closely singular coefficient matrices. Efficient schemes are developed to dynamically evolve the rank of the reduced solution and to ensure the orthogonality of the basis matrix while preserving its smooth evolution over time. Finally, the new schemes are applied to quantify the uncertain Lagrangian motions of a 2D double gyre flow with random frequency and of a stochastic flow past a cylinder.

A pressing environmental question facing the ocean is the potential impact of possible deep-sea mining activities. This work presents our initial results in developing an ocean and plume modeling system for the Bismark Sea where deep sea mining operations will probably take place. We employ the MSEAS modeling system to both simulate the ocean and to downscale initial conditions from a global system (HYCOM) and tidal forcing from the global TPXO-8 Atlas. We found that at least 1.5 km resolution was needed to adequately resolve the multiscale flow fields. In St. Georges channel, the interaction between the tides, background currents, and underlying density fields increased the subtidal flows. Comparing to historical transport estimates, we showed that tidal forcing is needed to maintain the correct subtidal transport through that Channel. Comparisons with past simulations and measured currents all showed good agreement between the MSEAS hindcasts. Quantitative comparisons made between our hindcasts and independent synoptic ARGO profiles showed that the hindcasts beat persistence by 33% to 50%. These comparisons demonstrated that the MSEAS current estimates were useful for assessing plume advection. Our Lagrangian transport and coherence analyses indicate that the specific location and time of the releases can have a big impact on their dispersal. Our results suggest that ocean mining plumes can be best mitigated by managing releases in accord with such ocean modeling and Lagrangian transport forecasts. Real-time integrated mining-modeling-sampling is likely to provide the most effective mitigation strategies.

The observation, computation and study of “Lagrangian Coherent Structures” (LCS) in turbulent geophysical flows have been active areas of research in fluid mechanics for the last 30 years. Growing evidence for the existence of LCSs in geophysical flows (e.g., eddies, oscillating jets, chaotic mixing) and other fluid flows (e.g., separation prole at the surface of an airfoil, entrainment and detrainment by a vortex) generates an increasing interest for the extraction and understanding of these structures as well as their properties. In parallel, realistic ocean modeling with dense data assimilation has developed in the past decades and is now able to provide accurate nowcasts and predictions of ocean flow fields to study coherent structures. Robust numerical methods and sufficiently fast hardware are now available to compute real-time forecasts of oceanographic states and render associated coherent structures. It is therefore natural to expect the direct predictions of LCSs based on these advanced models. The impact of uncertainties on the coherent structures is becoming an increasingly important question for practical applications. The transfer of these uncertainties from the ocean state to the LCSs is an unexplored but intriguing scientific problem. These two questions are the motivation and focus of this presentation. Using the classic formalism of continuous-discrete estimation [1], the spatially discretized dynamics of the ocean state vector x and observations are described by (1a) dx =M(x; t) + d yok (1b) = H(xk; tk) + k where M and H are the model and measurement model operator, respectively. The stochastic forcings d and k are Wiener/Brownian motion processes,   N(0;Q(t)), and white Gaussian sequences, k  N(0;Rk), respectively. In other words, Efd(t)d T (t)g := Q(t) dt. The initial conditions are also uncertain and x(t0) is random with a prior PDF, p(x(t0)), i.e. x(t0) = bx0 + n(0) with n(0) random. Of course, vectors and operators in Eqs. (1a-b) are multivariate which impacts the PDFs: e.g. their moments are also multivariate. The estimation problem at time t consists of combining all available information on x(t), the dynamics and data (Eqs. 1a-b), their prior distributions and the initial conditions p(x(t0)). Defining the set of all observations prior to time t by yt

A data and dynamics driven approach to estimate, decompose, organize and analyze the evolving three-dimensional variability of ocean fields is outlined. Variability refers here to the statistics of the differences between ocean states and a reference state. In general, these statistics evolve in time and space. For a first endeavor, the variability subspace defined by the dominant eigendecomposition of a normalized form of the variability covariance is evolved. A multiscale methodology for its initialization and forecast is outlined. It combines data and primitive equation dynamics within a Monte-Carlo approach. The methodology is applied to part of a multidisciplinary experiment that occurred in Massachusetts Bay in late summer and early fall of 1998. For a 4-day time period, the three-dimensional and multivariate properties of the variability standard deviations and dominant eigenvectors are studied. Two variability patterns are discussed in detail. One relates to a displacement of the Gulf of Maine coastal current offshore from Cape Ann, with the creation of adjacent mesoscale recirculation cells. The other relates to a Bay-wide coastal upwelling mode from Barnstable Harbor to Gloucester in response to strong southerly winds. Snapshots and tendencies of physical fields and trajectories of simulated Lagrangian drifters are employed to diagnose and illustrate the use of the dominant variability covariance. The variability subspace is shown to guide the dynamical analysis of the physical fields. For the stratified conditions, it is found that strong wind events can alter the structures of the buoyancy flow and that circulation features are more variable than previously described, on multiple scales. In several locations, the factors estimated to be important include some or all of the atmospheric and surface pressure forcings, and associated Ekman transports and downwelling/upwelling processes, the Coriolis force, the pressure force, inertia and mixing.