Passive Tracer Research Articles

Experimental Visualization and Modeling of the Transport Behaviors of Monofilament Microplastic Fibers Through an Idealized Porous Media

AbstractMicroplastic fibers (MPF) are the largest fraction of microplastics in the environment by mass. The endpoints of these contaminants' movement is generally known at large‐scale (i.e., their origins and where they end up), but the mechanics of how they get to those sinks remains poorly understood. The objective of this work was to improve understanding of the mechanisms driving MPF migration through terrestrial systems by directly imaging their motion through idealized representations of porous media. Monofilament line with 0.3 mm diameter was passed through a bench‐scale, pseudo‐2d flow cell to capture trajectories of MPFs of three different lengths and trajectories of passive micro‐bead tracers were also captured. Video processing and automated image analysis converted the video of the experiments into a database of trajectories, allowing comparison of the experimental data to various numerical models. Simple advection‐dispersion models were adequate for modeling the passive tracer but did not provide a good description of MPF transport. A physics‐based, distributed model was able to generate realistic trajectories through the domain, but the speeds of the fibers in the initial simulation were too fast, despite working well for the passive tracer. Adding a delay (waiting time) process resulted in good description of the trajectories and travel times. The specifics of the delay process could not be deduced from these experiments, but its overall impact on transport provides mechanistic insights. These direct observation of the trajectories and speeds of MPFs moving through porous media show that MPFs likely have strong interactions with their surroundings.

Vertical Carbon Export During a Phytoplankton Bloom in the Chukchi Sea: Physical Setting and Frontal Subduction

AbstractIn order to quantify pelagic‐benthic coupling on high‐latitude shelves, it is imperative to identify the different physical mechanisms by which phytoplankton are exported to the sediments. In June–July 2023, a field program documented the evolution of an under‐ice phytoplankton bloom on the northeast Chukchi shelf. Here, we use in situ data from the cruise, a simple numerical model, historical water column data, and ocean reanalysis fields to characterize the physical setting and describe the dynamically driven vertical export of chlorophyll associated with the bloom. A water mass front separating cold, high‐nutrient winter water in the north and warmer summer waters to the south—roughly coincident with the ice edge—supported a baroclinic jet which is part of the Central Channel flow branch that veers eastward toward Barrow Canyon. A plume of high chlorophyll fluorescence extending from the near‐surface bloom in the winter water downwards along the front was measured throughout the cruise. Using a passive tracer to represent phytoplankton in the model, it was demonstrated that the plume is the result of subduction due to baroclinic instability of the frontal jet. This process, in concert with the gravitational sinking, pumps the chlorophyll downwards an order of magnitude faster than gravitational sinking alone. Particle tracking using the ocean reanalysis fields reveals that a substantial portion of the chlorophyll away from the front is advected off of the northeast Chukchi shelf before reaching the bottom. This highlights the importance of the frontal subduction process for delivering carbon to the sea floor.

Onset of global instability in a premixed annular V-flame

We investigate self-excited axisymmetric oscillations of a lean premixed methane–air V-flame in a laminar annular jet. The flame is anchored near the rim of the centrebody, forming an inverted cone, while the strongest vorticity is concentrated along the outer shear layer of the annular jet. Consequently, the reaction and vorticity dynamics are largely separated, except where they coalesce near the flame tip. The global eigenmodes corresponding to the linearised reacting flow equations around the steady base state are computed in an axisymmetric setting. We identify an arc branch of eigenmodes exhibiting strong oscillations at the flame tip. The associated eigenvalues are robust with respect to domain truncation and numerical discretisation, and they become destabilised as the Reynolds number increases. The frequency of the leading eigenmode is found to correspond to the Lagrangian disturbance advection time from the nozzle outlet to the flame tip. The essential role of this convective mechanism is also supported by resolvent analysis, which finds that the same flame-tip disturbance structure and frequency are optimally amplified when the flame is subjected to external white noise forcing. Strong non-modal effects in the form of pseudo-resonance are not found. Nonlinear time-resolved simulation further reveals notable hysteresis phenomena in the subcritical regime prior to instability. Hence, even when the flame is linearly stable, perturbations of sufficient amplitude can trigger limit-cycle oscillations and higher-dimensional dynamics sustained by nonlinear feedback. A Monte Carlo simulation of passive tracers in the unsteady flame suggests a nonlinear non-local instability mechanism. Notably, linear analysis of the subcritical time-averaged limit-cycle state yields eigenvalues that do not match the nonlinear periodic oscillation frequencies. This mismatch is attributed to the fundamentally nonlinear dynamics of the subcritical V-flame instability, where the dichromatic, non-local interaction between the heat release rate along the flame surface and the vortex dynamics in the jet shear layer cannot be approximated as a simple distortion of the mean flow.

Recommended coupling to global meteorological fields for long-term tracer simulations with WRF-GHG

Abstract. Atmospheric transport models are often used to simulate the distribution of greenhouse gases (GHGs). This can be in the context of forward modeling of tracer transport using surface–atmosphere fluxes or flux estimation through inverse modeling, whereby atmospheric tracer measurements are used in combination with simulated transport. In both of these contexts, transport errors can bias the results and should therefore be minimized. Here, we analyze transport uncertainties in the commonly used Weather Research and Forecasting (WRF) model coupled with the greenhouse gas module (WRF-GHG), enabling passive tracer transport simulation of CO2 and CH4. As a mesoscale numerical weather prediction model, WRF's transport is constrained by global meteorological fields via initialization and at the lateral boundaries of the domain of interest. These global fields were generated by assimilating various meteorological data to increase the accuracy of modeled fields. However, in limited-domain models like WRF, the winds in the center of the domain can deviate considerably from these driving fields. As the accuracy of the wind speed and direction is critical to the prediction of tracer transport, maintaining a close link to the observations across the simulation domain is desired. On the other hand, a link that is too close to the global meteorological fields can degrade performance at smaller spatial scales that are better represented by the mesoscale model. In this work, we evaluated the performance of strategies for keeping WRF's meteorology compatible with meteorological observations. To avoid the complexity of assimilating meteorological observations directly, two main strategies of coupling WRF-GHG with ERA5 meteorological reanalysis data were tested over a 2-month-long simulation over the European domain: (a) restarting the model daily with fresh initial conditions (ICs) from ERA5 and (b) nudging the atmospheric winds, temperatures, and moisture to those of ERA5 continuously throughout the simulation period, using WRF's built-in four-dimensional data assimilation (FDDA) in grid-nudging mode. Meteorological variables and simulated mole fractions of CO2 and CH4 were compared against observations to assess the performance of the different strategies. We also compared planetary boundary layer height (PBLH) with radiosonde-derived estimates. Either nudging or daily restarts similarly improved the meteorology and GHG transport in our simulations, with a small advantage of using both methods in combination. However, notable differences in soil moisture were found that accumulated over the course of the simulation when not using frequent restarts. The soil moisture drift had an impact on the simulated PBLH, presumably via changing the Bowen ratio. This is partially mitigated through nudging without requiring daily restarts, although not entirely alleviated. Soil moisture drift did not have a noticeable impact on GHG performance in our case, likely because it was dominated by other errors. However, since the PBLH is critical for accurately simulating GHG transport, we recommend transport model setups that tie soil moisture to observations. Our method of frequently re-initializing simulations with meteorological reanalysis fields proved suitable for this purpose.

Impact of ageostrophic dynamics on the predictability of Lagrangian trajectories in surface-ocean turbulence

Turbulent flows at the surface of the ocean deviate from geostrophic equilibrium on scales smaller than about 10 km. These scales are associated with important vertical transport of active and passive tracers, and should play a prominent role in the heat transport at climate scales and for plankton dynamics. Measuring velocity fields on such small scales is notoriously difficult but new, high-resolution satellite altimetry is starting to reveal them. However, the satellite-derived velocities essentially represent the geostrophic flow component, and the impact of unresolved ageostrophic motions on particle dispersion needs to be understood to properly characterize transport properties. Here, we investigate ocean fine-scale turbulence using a model that represents some of the processes due to ageostrophic dynamics. We take a Lagrangian approach and focus on the predictability of the particle dynamics, comparing trajectories advected by either the full flow or by its geostrophic component only. Our results indicate that, over long times, relative dispersion is marginally affected after filtering out the ageostrophic component. Nevertheless, advection by the geostrophic-only flow leads to an overestimation of the typical pair-separation rate, and to a bias on trajectories (in terms of displacement from the actual ones), whose importance grows with the Rossby number. We further explore the intensity of the transient particle clustering induced by ageostrophic motions and find that it can be significant, even for small flow compressibility. Indeed, we show that clustering is here due to the interplay between compressibility and persistent flow structures that trap particles, enhancing their aggregation. Published by the American Physical Society 2024

The formation and ventilation of an oxygen minimum zone in a simple model for latitudinally alternating zonal jets

Abstract. An advection–diffusion model coupled to a simple dynamical ocean model is used to illustrate the formation and ventilation of an oxygen minimum zone. The advection–diffusion model carries a passive tracer with a source at the western boundary and a Newtonian damping term to mimic oxygen consumption. The dynamical model is a non-linear one-and-a-half-layer reduced-gravity model. The latter is forced by an annually oscillating mass flux confined to the near-equatorial band that, in turn, leads to the generation of mesoscale eddies and latitudinally alternating zonal jets at higher latitudes. The model uses North Atlantic geometry and develops a tracer minimum zone remarkably similar in location to the observed oxygen minimum zone in the eastern tropical North Atlantic (ETNA). This is despite the absence of wind forcing for a subtropical gyre and the shadow zone predicted by the ventilated thermocline theory. Although the model is forced only at the annual period, the model nevertheless exhibits decadal and multidecadal variability in its spun-up state. The associated trends are comparable to observed trends in oxygen within the ETNA oxygen minimum zone. Notable exceptions are the multi-decadal decrease in oxygen in the lower-oxygen minimum zone and the sharp decrease in oxygen in the upper-oxygen minimum zone between 2006 and 2013.

Nonlinear changes in urban heat island intensity, urban breeze intensity, and urban air pollutant concentration with roof albedo

This study systematically examines how the urban heat island (UHI) and urban breeze circulation (UBC) respond to an increase in roof albedo (αr) and its influence on urban air pollutant dispersion. For this, idealized ensemble simulations are performed using the Weather Research and Forecasting (WRF) model. The increase in αr from 0.20 to 0.65 decreases the UHI intensity, UBC intensity, and urban planetary boundary layer (PBL) height in the daytime (from 1200 to 1700 LST) by 47%, 36%, and 6%, respectively. As both UBC intensity and urban PBL height decrease, the daytime urban near-surface passive tracer concentration increases by 115%. The daytime UHI intensity, UBC intensity, and urban tracer concentration nonlinearly change with αr: For 0.10 ≤ αr < 0.80, the rates of changes in the UHI intensity, UBC intensity, and urban tracer concentration with αr overall increase as αr increases. For αr ≥ 0.80, the daytime roof surface temperature is notably lower than the daytime urban near-surface air temperature, the UHI intensity, UBC intensity, and urban tracer concentration very slightly changing with αr. This study provides insights into the associations between changes in roof surface temperature and roof surface energy fluxes with αr and those in UHI intensity.

Toward production zonal flow allocation using novel completion tracer dosing line and diffuser: Prototype laboratory flow loop experiments and computational fluid dynamics studies

Complex completion schemes with inflow control valves (ICVs) in multi-lateral wells are now the industry standard for modern fields. Understanding the production contributions from different production zones/laterals is of great importance in optimizing the zonal flow rates to prolong the field lifetimes, reduce water cuts and the associated costs, and achieve production targets. Furthermore, the improved operational efficiency can reduce the producers’ energy usages and carbon footprints for a more sustainable operation. For the past decade, inflow tracer technology based on embedding tracers in degradable polymer matrices, which are then placed in downhole completion zones at initial installation, has been used in understanding production zonal flow contributions in field demonstrations as well as in commercial projects. However, the current technology based on degradable polymers has the following limitations: (1) there are limited lifetimes for the degradable polymer matrices that will lose their functionality past their expiration dates, and (2) the tracers are released from the polymer matrices at a constant rate; so in order to create a transient for tracer flowback measurements the wells need to be shut-in for a period of time which will interfere with normal production and possibly, with the long-term productivity of the well. Here, we introduce an alternative approach for deploying inflow tracers that has no such limitations mentioned above. Specifically, we propose the installation of passive chemical tracer dosing lines connected to diffuser chambers at the annular regions downhole during completion. Pulses of tracers can be injected into the different production zones and monitored for their flowback characteristics on-demand. The passive tracer dosing lines and diffusers only require minimum upfront investments and have no limitations of lifetime compared to the degradable polymer matrices. Furthermore, the injection of tracer pulses into the annulus near production zones creates appropriate transient responses for their flowback measurements, which do not require well shut-in and will not affect normal production.Laboratory flow loop experiments and computational fluid dynamics (CFD) studies were performed to demonstrate the proposed methodology and to understand the underlying physics. Half-scale and quarter-scale flow loops made of transparent polyvinyl chloride (PVC) pipes that resemble annular regions near the ICVs were constructed with tracer dosing lines and diffusers attached within. Colored dyes were injected into the flow loops to visualize their flush out and flowback behaviors. In CFD studies, a two-step approach was used with first solving for the turbulent fluid flow profiles, and then solving for the tracer convection diffusion behaviors in the turbulent fields. The CFD model parameters were carefully matched with flow loop experiments to predict the possible real field responses. With appropriate diffuser design, we can achieve desirable zonal flow-rate-dependent tracer flowback responses with reasonable sampling times and injected tracer quantities.

Diagnosing tracer transport in convective penetration of a stably stratified layer

We use large-eddy simulations to study the penetration of a buoyant plume carrying a passive tracer into a stably stratified layer with constant buoyancy frequency. Using a buoyancy-tracer volume distribution, we develop a method for objectively partitioning plume fluid in buoyancy-tracer space into three regions, each of which corresponds to a coherent region in physical space. Specifically, we identify a source region where undiluted plume fluid enters the stratified layer, a transport region where much of the transition from undiluted to mixed fluid occurs in the plume cap and an accumulation region corresponding to a radially spreading intrusion. This method enables quantification of different measures of turbulence and mixing within each of the three regions, including potential energy and turbulent kinetic energy dissipation rates, an activity parameter and the instantaneous mixing efficiency. We find that the most intense buoyancy gradients lie in a thin layer at the cap of the penetrating plume. This provides the primary stage of mixing between plume and environment and exhibits a mixing efficiency around 50 %. Newly generated mixtures of environmental and plume fluid join the intrusion and experience relatively weak turbulence and buoyancy gradients. As the intrusion spreads radially, environmental fluid surrounding the intrusion is mixed into the intrusion with moderate mixing efficiency. This dominates the volume of environmental fluid entrained into the region containing plume fluid. However, the ‘strongest’ entrainment, as measured by the specific entrainment rate, is largest in the plume cap, where the most buoyant environmental fluid is entrained.

Enhanced Diffusion and Non-Gaussian Displacements of Colloids in Quasi-2D Suspensions of Motile Bacteria

In the real world, active agents interact with surrounding passive objects, thus introducing additional degrees of complexity. The relative contributions of far-field hydrodynamic and near-field contact interactions to the anomalous diffusion of passive particles in suspensions of active swimmers remain a subject of ongoing debate. We constructed a quasi-two-dimensional microswimmer–colloid mixed system by taking advantage of Serratia marcescens’ tendency to become trapped at the air–water interface to investigate the origins of the enhanced diffusion and non-Gaussianity of the displacement distributions of passive colloidal tracers. Our findings reveal that the diffusion behavior of colloidal particles exhibits a strong dependence on bacterial density. At moderate densities, the collective dynamics of bacteria dominate the diffusion of tracer particles. In dilute bacterial suspensions, although there are multiple dynamic types present, near-field contact interactions such as collisions play a major role in the enhancement of colloidal transport and the emergence of non-Gaussian displacement distributions characterized by heavy exponential tails in short times. Despite the distinct types of microorganisms and their diverse self-propulsion mechanisms, a generality in the diffusion behavior of passive colloids and their underlying dynamics is observed.

Uncertainties in exposure predictions arising from point measurements of carbon dioxide in classroom environments.

Predictions of airborne infection risk can be made based on the fraction of rebreathed air inferred from point measurements of carbon dioxide (CO[Formula: see text]). We investigate the extent to which environmental factors, particularly spatial variations due to the ventilation provision, affect the uncertainty in these predictions. Spatial variations are expected to be especially problematic in naturally ventilated spaces, which include the majority of classrooms in the UK. An idealized classroom, broadly representative of the physics of (buoyancy-driven) displacement ventilation, is examined using computational fluid dynamics, with different ventilation configurations. Passive tracers are used to model both the CO[Formula: see text] generated by all 32 occupants and the breath of a single infectious individual (located in nine different regions). The distribution of infected breath is shown to depend strongly on the distance from the release location but is also affected by the pattern of the ventilating flow, including the presence of stagnating regions. However, far-field exposure predictions based on single point measurements of CO[Formula: see text] within the breathing zone are shown to rarely differ from the actual exposure to infected breath by more than a factor of two-we argue this uncertainty is small compared with other uncertainties inherent in modelling airborne infection risk.

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

Direct Entrainment as a Measure of Dilution

Abstract The turbulent transport of mass, energy, moisture, and momentum between the clouds and the surrounding environment plays a central role in determining the vertical structure of the troposphere. This article investigates the connection between net entrainment, defined here as the net transport of air into individual clouds, and net dilution, defined as the tendencies of passive tracers such as static energy or total water mixing ratio. Entrainment and detrainment rates for 2.6 × 106 individual cloud samples are obtained from a large-eddy simulation of shallow convective boundary layer atmosphere that explicitly calculates the turbulent fluxes across the cloud boundaries. The equations describing the tendencies of cloud mass and tracer concentrations are derived as a function of directly calculated entrainment and detrainment rates of the individual clouds, and used to calculate net entrainment and dilution rates. Directly calculated net entrainment and dilution rates agree well with cloud mass and tracer tendencies and give a dilution time scale of 13 min. In contrast, the traditional bulk-plume approximation overestimates the effect of entrainment and detrainment on the dilution of cloud field properties due to the differential tracer transport through the moist shell surrounding the cloud. Using direct measures of entrainment and detrainment for individual clouds separates different processes that influence the turbulent mass transport between the clouds and their environment.

Horizontal and vertical dispersion in a wind-driven oceanic gyre model

This study addresses the horizontal and vertical dispersion of passive tracers in idealized wind-driven subtropical gyres. Synthetic particles within a closed basin are numerically advected to analyze their dispersion under different theoretical velocity fields. Horizontal dispersion simulations incorporate the classic wind-driven Stommel circulation along with (i) surface Ekman drift associated with the Stommel wind field and (ii) inertial effects due to particle size and buoyancy. Results reveal that the Ekman drift inhibits particle dispersion across the entire domain leading to tracer concentration in a quasi-stable distribution skewed toward the western side of the basin. Similar behavior is observed with inertial particles. The equilibrium state is quantified for different diffusivity values, particle sizes, and buoyancies. For vertical dispersion, simulations incorporate the three-dimensional Ekman velocity, which includes a negative vertical component, while ignoring inertial effects. Initially, surface particles accumulate around the gyre center while slowly sinking, but they disperse across the basin once they surpass the Ekman layer and are free from surface effects. Tracers sink more on the western side of the basin, regardless of horizontal diffusivity. On average, ignoring inertial effects, particles sink less with higher diffusivity and more with lower diffusivity, suggesting a potential for high horizontal distribution of sunken tracers in the ocean.

The effect of Southern Ocean topography on the global MOC and abyssal water mass distribution

Abstract We investigate the role of Southern Ocean topography and wind stress in the deep and abyssal ocean overturning and water mass composition using a suite of idealized global ocean circulation models. Specifically, we address how the presence of a meridional ridge in the vicinity of Drake Passage and the formation of an associated Southern Ocean gyre influences the water mass composition of the abyssal cell. Our experiments are carried out using a numerical representation of the global ocean circulation in an idealized two-basin geometry under varying wind-stress and Drake Passage ridge height. In the presence of a low Drake Passage ridge the overall strength of the meridional overturning circulation is primarily influenced by wind stress, with a topographically-induced weakening of the mid-depth cell and concurrent strengthening of the abyssal cell occurring only after ridge height passes 2500m. Passive tracer experiments show that a strengthening mid-depth cell leads to increased abyssal ventilation by North Atlantic water masses, as more North Atlantic Deep Water (NADW) enters the Southern Ocean and then spreads into the Indo-Pacific. We repeat our tracer experiments without restoring in the high-latitude Southern Ocean in order to identify the origin of water masses that circulate through the Southern Ocean before sinking into the abyss as Antarctic Bottom Water. Our results from these “exchange” tracer experiments show that an increasing ridge height in Drake Passage and the concurrent gyre spin-up lead to substantially decreased NADW-origin waters in the abyssal ocean, as more surface waters from north of the ACC are transferred into the Antarctic Bottom Water formation region.

Trough‐Scale Slope Countercurrent Over the East China Sea Continental Slope Driven by Upwelling Divergence

AbstractObservations have revealed the existence of persistent slope countercurrents (SCCs) that flow southwestward beneath the Kuroshio Current at several locations over the East China Sea (ECS) continental slope. It was not clear whether these flows are localized circulation features or segments of a trough‐scale circulation system in the Okinawa Trough (OT). We demonstrate that there indeed exists a potentially continuous trough‐scale SCC along the ECS slope that is associated with an OT‐wide cyclonic circulation using high‐resolution model simulations and physical interpretations. The detailed features of the deep OT circulation are illustrated by the trajectories of the Lagrangian drifters and the time‐varying distributions of passive tracers. The SCC in the ECS is characterized by its weak yet persistent nature, typically located in narrow sloping regions at the isopycnal layer of 26.6–27.3 kg m−3. It exhibits a characteristic speed of approximately O‐(1) cm s−1. Analyses and experiments suggest that the divergence of upwelling in the SCC layer (26.6–27.3 σθ surface) gives rise to lateral potential vorticity transport, ultimately driving the deep cyclonic circulation. Furthermore, the SCC also displays a substantial connection with the onshore intrusion of the Kuroshio Current, particularly to the northeast of Taiwan Island. The SCC may potentially play a crucial role in the transport of heat and nutrients, as well as in regulating sediment distributions within the deep OT. This mechanism offers fresh insights into explaining the presence of undercurrents in semi‐enclosed marginal seas.

Large eddy simulation analysis of a model reactive tracer through spatial filtering

Large eddy simulations (LES) provide a methodology for both analyzing and simulating multi-scale flows when the smallest scales of motion cannot be resolved. Within environmental flows there exist numerous biogeochemical processes involving tracers undergoing reactions. In this study, we perform an a posteriori LES analysis on a direct numerical simulation of an idealized model reactive tracer subjected to three-dimensional turbulence induced by a Rayleigh–Taylor instability. The governing equations, including an advection–diffusion–reaction equation for the reactive tracer, are filtered, and the resulting sub-filter-scale terms are expressed in terms of interactions between scales. The procedure is demonstrated for a generalized degree N polynomial reaction function. Various spectral filters are applied to the data and compared. The preferential choice is to use the widest filter possible with a smoothed cutoff. The sub-filter-scale reaction term that results from filtering the reaction function is considered for each of the filter choices. When using a particularly harsh filter, local balances are found for the resolved scale and cross-scale components of the sub-filter-scale reaction term. The same result is shown for the vertical sub-filter-scale flux for both a reactive and a passive tracer. The components of the sub-filter-scale reaction and vertical flux terms involving interactions at the sub-filter-scale do not show any evidence of local balances and are distributed around the fine turbulent structures in the flow. This suggests that parameterizations for the sub-filter-scale terms would benefit from considering event specific dynamics.

Seasonal Tracer Subduction in the Subpolar North Atlantic Driven by Submesoscale Fronts

AbstractSubmesoscale flows (0.1–10 km) are often associated with large vertical velocities, which can have a significant impact on the transport of surface tracers, such as carbon. However, global models do not adequately account for these small‐scale effects, which still require a proper parameterization. In this study, we introduced a passive tracer into the surface mixed layer (ML) of a northern Atlantic Ocean simulation based on the primitive‐equation model CROCO, with a horizontal resolution of Δx = 800 m, aiming to investigate the seasonal submesoscale effects on vertical transport. Using surface vorticity and strain rate criteria, we identified regions with submesoscale fronts and quantified the associated subduction, that is the export of tracer below the ML depth. The results suggest that the tracer vertical distribution and the contribution of frontal subduction can be estimated from surface strain and vorticity. Notably, we observed significant seasonal variations. In winter, the submesoscale fronts contribute up to 40% of the total vertical advective transport of tracer below the ML, while representing only 5% of the domain. Conversely, in summer, fronts account for less than 1% of the domain and do not contribute significantly to the transport below the ML. The findings of this study contribute to a better understanding of the seasonal water subduction due to fronts in the region.

Fusion of Lagrangian drifter data and numerical model outputs for improved assessment of turbulent dispersion

Abstract. Transport and dispersion processes in the ocean are crucial, as they determine the lifetime and fate of biological and chemical quantities drifting with ocean currents. Due to the complexity of the coastal ocean environment, numerical circulation models have difficulties to accurately simulate highly turbulent flows and dispersion processes, especially in highly energetic tidal basins such as the eastern English Channel. A method of improving the results of coastal circulation modeling and tracer dispersion in the Dover Strait is proposed. Surface current velocities derived from Lagrangian drifter measurements in November 2020 and May 2021 were optimally interpolated in time and space to constrain a high-resolution coastal circulation MARS model, with careful attention given to selecting ensemble members composing the model covariance matrix. The space–time velocity covariances derived from model simulations were utilized by the optimal interpolation algorithm to determine the most likely evolution of the velocity field under constraints provided by Lagrangian observations and their error statistics. The accuracy of the velocity field reconstruction was evaluated at each time step. The results of the fusion of model outputs with surface drifter velocity measurements show a significant improvement (by ∼ 50 %) of the model capability to simulate the drift of passive tracers in the Dover Strait. Optimized velocity fields were used to quantify the absolute dispersion in the study area. The implications of these results are important, as they can be used to improve existing decision-making support tool or design new tools for monitoring the transport and dispersion in a coastal ocean environment.

Flow structure around a microswimmer at fluid–fluid interface

Living organisms such as bacteria and algae often form biofilms at air–liquid and/or liquid–liquid interfaces. Therefore, it is important to understand the hydrodynamic interaction between the fluid–fluid interface and microorganisms. In the present study, a squirmer model is employed to study the flow field structures around a symmetrically trapped spherical microswimmer at an interface separating two fluids with viscosity contrast. In these simulations, Reynolds ( Re ) and Capillary ( Ca ) numbers are very small, and hence the contribution of inertia and interface deformations are neglected. The squirmer model is validated against the analytical solution obtained by considering a source dipole. It is observed that the flow structure and vorticity distribution around the microswimmers are strongly influenced by the squirmer parameter (β) and viscosity contrast (λ). Furthermore, the interplay between force-dipole and source-dipole along with viscosity contrast leads to a range of flow structures such as symmetric four-lobe and asymmetric quadrupolar flow fields. This flow structure asymmetry around the swimmer is quantified with different steady state orientations (φ). The effect of flow profile on the passive tracer particle transport is explored in terms of trajectories. Specifically, larger values of φ result in more profound curly trajectories.