Discovery and analysis of submesoscale variability O(0.3–30) km on the continental shelf is made possible by a high-resolution (Δx = 75 m) Regional Oceanic Modeling System (ROMS) simulation of the Southern California Bight (SCB). This variability is manifest in ubiquitous yet ephemeral coherent structures: fronts, filaments, and vortices. Similar to their open-ocean counterparts, fronts and filaments on the shelf are identified by their strong vertical velocity, surface convergence, cyclonic vorticity, and horizontal density gradient. Life cycles of these features typically last 3–5 days, with the formation dominated by a horizontal advective tendency that increases density and velocity gradients (i.e., frontogenesis). The shape of the coastline and depth of the water column both influence the abundance and spatial orientation of shallow-water fronts and filaments. Closer to shore, fronts and filaments often align themselves parallel to isobaths, and headlands often act as sites of intense vorticity generation through bottom stress. A quasi-steady, approximate momentum balance among rotation, pressure gradient, and vertical mixing—known as turbulent thermal wind (TTW)—often is valid in the strong secondary circulations local to fronts and filaments. However, front and filament circulations subject to strong diurnal variation in surface heating and vertical mixing are inconsistent with steady-state TTW balance. The secondary circulations can induce ephemeral material trapping and substantial vertical heat fluxes on the shelf.

Realistic simulation of nearshore (from the shoreline to approximately 10‐km offshore) Lagrangian material transport is required for physical, biological, and ecological investigations of the coastal ocean. Recently, high‐resolution simulations of the coastal ocean have revealed a shelf populated with small‐scale, rapidly evolving currents that arise at resolutions ⪅100 m. However, many historical and recent investigations of coastal connectivity utilize circulation models with ≈1‐km resolution. Here we show a resolution sensitivity to simulated Lagrangian transport and coastal connectivity with a hierarchy of Regional Oceanic Modeling System simulations of the Santa Barbara Channel at Δx= 1, 0.3, 0.1, and 0.036 km. At higher resolution ( normalΔx⪅ 100 m), rapid alongshore and vertical transport occurs in regions less than 1 km from the shoreline due to submesoscale shelf currents that open up new transport pathways on the shelf: submesoscale fronts and filaments, topographic wakes, and narrow alongshore jets. Shallow‐water fronts and filaments induce early time downwelling and subsequent dispersal at depth of surface material; this is not captured at coarser resolution (Δx= 1 km). Differences in three‐dimensional and two‐dimensional transport are explored in a higher‐resolution simulation: In general, three‐dimensional trajectories are more dispersive than two‐dimensional, due to a separation in their respective trajectories.

The local circulation of submesoscale fronts and filaments can be partly understood through a horizontal momentum balance of Coriolis, a horizontal pressure gradient, and vertical diffusivity in a turbulent boundary layer, known as the turbulent thermal wind balance (TTW or T2W). T2W often reproduces the instantaneous relative vorticity and divergence of submesoscale circulations in open-ocean and shelf settings. However, a diurnal cycle in submesoscale vorticity and divergence is characterized by a non-T2W phasing: a maximum in divergence magnitude lags the maximum in vertical diffusivity (with vorticity lagging divergence). Here, an idealized model is used to solve the transient turbulent thermal wind (T3W) equations and to investigate the diurnal evolution of front and filament circulation in a 2D plane. Relative to a steady-state circulation, transient evolution can cause both instantaneous strengthening and a weaker diurnal average of the secondary circulation. The primary mechanisms controlling the diurnal variability exist in a 1D Ekman layer that imprints onto the 2D circulation. In midlatitudes, acceleration during separate phases of the diurnal cycle (from night to day and from day to night) is dominated by distinct inertial oscillation and vertically diffusive mechanisms, respectively. However, the manifestation of these dual accelerations is sensitive to latitude, boundary layer depth, and the strength of the forcing. A simple 1D model predicts the diurnal phasing of submesoscale divergence and vorticity in realistic primitive equation simulations of the southwestern Pacific and coastal California.

The diurnal cycling of submesoscale circulations in vorticity, divergence, and strain is investigated using drifter data collected as part of the Lagrangian Submesoscale Experiment (LASER) experiment, which took place in the northern Gulf of Mexico during winter 2016, and ROMS simulations at different resolutions and degree of realism. The first observational evidence of a submesoscale diurnal cycle is presented. The cycling is detected in the LASER data during periods of weak winds, whereas the signal is obscured during strong wind events. Results from ROMS in the most realistic setup and in sensitivity runs with idealized wind patterns demonstrate that wind bursts disrupt the submesoscale diurnal cycle, independently of the time of day at which they happen. The observed and simulated submesoscale diurnal cycle supports the existence of a shift of approximately 1–3 h between the occurrence of divergence and vorticity maxima, broadly in agreement with theoretical predictions. The amplitude of the modeled signal, on the other hand, always underestimates the observed one, suggesting that even a horizontal resolution of 500 m is insufficient to capture the strength of the observed variability in submesoscale circulations. The paper also presents an evaluation of how well the diurnal cycle can be detected as function of the number of Lagrangian particles. If more than 2000 particle triplets are considered, the diurnal cycle is well captured, but for a number of triplets comparable to that of the LASER analysis, the reconstructed diurnal cycling displays high levels of noise both in the model and in the observations.

Offshore aquaculture has the potential to expand the macroalgal industry. However, moving into deeper waters requires suspended structures that will present novel farm-environment interactions. Here, we present a computational modeling framework, the Macroalgal Cultivation Modeling System (MACMODS), to explore within-farm modifications to light, seawater flow, and nutrient fields across time and space scales relevant to macroalgae. A regional ocean model informs the site-specific setting, the Santa Barbara Channel in the Southern California Bight. A fine-scale hydrodynamic model predicts modified flows and turbulent mixing within the farm. A spatially resolved macroalgal growth model, parameterized for giant kelp, Macrocystis pyrifera, predicts kelp biomass. Key findings from model integration are that regional ocean conditions set overall farm performance, while fine-scale within-farm circulation and nutrient delivery are important to resolve variation in within-farm macroalgal performance. Therefore, we conclude that models resolving within-farm dynamics can provide benefit to farmers with insight on how farm design and regional ocean conditions interact to influence overall yield. Here, the presence of repeating longlines aligned with the mean current generate flow diversions around the farm as well as attached Langmuir circulations and increased turbulence intensity. These flow-induced phenomena lead to less biomass in the interior portion of the farm relative to the edges. We also find that there is an effluent “footprint” that extends as much as 20 km beyond the farm. In this regard, MACMODS can be used to not only evaluate farm design and cultivation practices that maximize yield but also explore interactions between the farm and ecosystem in order to minimize impacts.