Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “submesoscale”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Isopycnal Submesoscale Stirring Crucially Sustaining Subsurface Chlorophyll Maximum in Ocean Cyclonic Eddies

Abstract Mesoscale and submesoscale processes have crucial impacts on ocean biogeochemistry, importantly enhancing the primary production in nutrient‐deficient ocean regions. Yet, the intricate biophysical interplay still holds mysteries. Using targeted high‐resolution in situ observations in the South China Sea, we reveal that isopycnal submesoscale stirring serves as the primary driver of vertical nutrient transport to sustain the dome‐shaped subsurface chlorophyll maximum (SCM) within a long‐lived cyclonic mesoscale eddy. Density surface doming at the eddy core increased light exposure for phytoplankton production, while along‐isopycnal submesoscale stirring disrupted the mesoscale coherence and drove significant vertical exchange of tracers. These physical processes play a crucial role in maintaining the elevated phytoplankton biomass in the eddy core. Our findings shed light on the universal mechanism of how mesoscale and submesoscale coupling enhances primary production in ocean cyclonic eddies, highlighting the pivotal role of submesoscale stirring in structuring marine ecosystems.

Cao, Haijin↗

Response of the Atmospheric Boundary Layer to Submesoscale Sea Surface Temperature Fronts

Submesoscale sea surface temperature fronts are ubiquitous throughout much of the global ocean; however, the response of the marine atmospheric boundary layer (MABL) to the ocean submesoscale is not well understood. Here large-eddy simulation is used to explore the time-dependent response of the MABL to idealized submesoscale sea surface temperature fronts, with an emphasis on how the dynamics of the MABL determine the strength and position of gradients in wind speed and air temperature. Results suggest that horizontal mixing only becomes important in response to frontogenesis by horizontally convergent ageostrophic flows, contrary to the common assumption that the MABL response will be strongly dependent on horizontal turbulent mixing. The fronts that develop in the MABL are also associated with large vertical relative vorticity, suggesting the possibility that submesoscale fronts may induce inertial instability in the MABL. Finally, these results provide guidance for high-resolution ocean and atmosphere modeling and for interpreting observations.

58 GEOSCIENCES↗

Frontal Maintenance in Submesoscale Flows

Classic deformation theory includes parameters—divergence, total strain, and vorticity—that are invariant to changes in the coordinate system. However, these parameters are sometimes ambiguous with respect to characterizing how fronts are formed and maintained because the presence of a front imposes a reference coordinate system. To help remedy this ambiguity, we propose a framework in frontal coordinates based on along- and cross-front velocity gradients to better characterize frontal maintenance, which can also be used to define divergence and normal strain in frontal coordinates. The framework with these four parameters (along-, cross-front velocity gradients, divergence, and normal strain in frontal coordinate) defines eight characteristic flow types at a front, providing a complete characterization of the flow that strengthens or weakens a front. This framework highlights the importance of the “strain efficiency” concept, which unambiguously defines the contribution of total strain to frontogenesis. Two examples, one based on a realistic simulation of submesoscales in the northern Gulf of Mexico and the other based on an idealized model with similar flow characteristics, are provided to demonstrate how this framework can be used to enhance our understanding of frontal dynamics in submesoscale flows.

frontogenesis↗

The Scale and Activity of Symmetric Instability Estimated from a Global Submesoscale-Permitting Ocean Model

Abstract Symmetric instability (SI) extracts kinetic energy from fronts in the surface mixed layer (SML), potentially affecting the SML structure and dynamics. Here, a global submesoscale-permitting ocean model named MITgcm LLC4320 simulation is used to examine the Stone linear prediction of the maximum SI scale to estimate grid spacings needed to begin resolving SI. Furthermore, potential effects of SI on the usable wind work are estimated roughly: this estimate of SI “activity” is useful for assessing if these modes should be resolved or parameterized. The maximum SI scale varies by latitude with median values from 568 to 23 m. Strong seasonality is observed in the SI scale and activity. The median scale in winter is 188 m globally, 2.5 times of that of summer (75 m). SI is more active in winter: 15% of the time compared with 6% in summer. The strongest SI activity is found in the western Pacific, western Atlantic, and Southern Oceans. The required grid spacings for a global model to begin resolving SI eddies in the SML are 24 m (50% of regions resolved) and 7.9 m (90%) in winter, decreasing to 9.4 m (50%) and 3.6 m (90%) in summer. It is also estimated that SI may reduce usable wind work by an upper bound of 0.83 mW m −2 globally, or 5% of the global magnitude. The sensitivity of these estimates to empirical thresholds is provided in the text.

Dong, Jihai↗

The DOE E3SM version 2.1: overview and assessment of the impacts of parameterized ocean submesoscales

The U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 2.1 builds on E3SMv2 with several changes, with the most notable being the addition of the Fox-Kemper et al. (2011) mixed-layer eddy parameterization. This parameterization captures the effect of finite-amplitude, mixed-layer eddies as an overturning streamfunction and has the primary function of restratification. Herein, we outline the changes to the mean climate state of E3SM that were introduced by the addition of this parameterization. Overall, the presence of the submesoscale parameterization improves the fidelity of the v2.1 simulation by reducing the ocean surface biases in the North Atlantic present in v2, as illustrated by changes in the climatological sea surface temperature and salinity and the Arctic sea-ice extent. Other impacts include a slight shoaling of the mixed-layer depths in the North Atlantic and a small improvement in the Atlantic Meridional Overturning Circulation (AMOC). We note that the expected shoaling due to the parameterization is regionally dependent in our coupled configuration. In addition, we investigate why the parameterization and its impacts on mixed-layer depth have little impact on the simulated AMOC: despite increased dense-water formation in the Norwegian Sea, only a small fraction of the water formed makes its way south into the North Atlantic basin. Version 2.1 also exhibits small improvements in the atmospheric climatology, with smaller biases in many notable quantities and modes of variability.

54 ENVIRONMENTAL SCIENCES↗

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗

Amplified Mesoscale and Submesoscale Variability and Increased Concentration of Precipitation under Global Warming over Western North America

Abstract Cold-season precipitation statistics in simulations from the storm-resolving WRF Model at 6-km and 1-h resolution over western North America are analyzed. Pseudo–global warming future simulations for the 2041–80 period, constrained by GCMs under the RCP8.5 scenario, are compared to the 1981–2020 historical simulation. The analysis focuses on the dynamical properties of precipitation time series at subdaily scales and on the morphology of storms. The statistical distribution of precipitation intensities in each pixel of the simulation domain is characterized through nonparametric statistical indicators: frequency of wet hours, mean wet-hour precipitation intensity, and Gini coefficient as a measure of the temporal concentration of the precipitation volume. Additionally, the temporal and spatial Fourier power spectra of precipitation time series and precipitation fields are analyzed. The half-power period (HPP) and half-power wavelength (HPW) are defined as spectral measures of the characteristic scales of precipitation’s temporal and spatial patterns. The results show statistically significant increases in the mean wet-hour precipitation intensity and in the Gini coefficient in 99% of the pixels, indicating that the seasonal precipitation volume becomes more concentrated within a smaller number of hours with higher precipitation intensity. The statistics of change in the frequency of wet hours are more contrasted across the simulation domain. The changes are also reflected in the power spectra, which show the spatial and temporal variability increasing proportionally more with finer spatial and temporal scales and the HPW and HPP decreasing. These projected changes are expected to have consequences, not only in terms of hydrologic impacts but also in terms of the predictability of precipitation patterns. Significance Statement The precipitation characteristics of winter storms over the western United States and southwestern Canada are analyzed in future climate simulations for the 2041–80 period. As compared to present-day climate, the most intense parts of the storms are projected to produce a higher rainfall volume, with increased concentration over smaller areas and shorter time intervals. The propensity of rainfall intensity to vary rapidly over time will be enhanced in the future according to the simulations. These model predictions imply an increased risk of rapid flooding in small basins. They also suggest that predicting several hours ahead the time and location at which a storm will produce maximum rainfall may become more challenging in the future.

Climate change↗

Surface waves affect frontogenesis

Here, this paper provides a detailed analysis of momentum, angular momentum, vorticity, and energy budgets of a submesoscale front undergoing frontogenesis driven by an upper–ocean, submesoscale eddy field in a Large Eddy Simulation (LES). The LES solves the wave–averaged, or Craik–Leibovich, equations in order to account for the Stokes forces that result from interactions between nonbreaking surface waves and currents, and resolves both submesoscale eddies and boundary layer turbulence down to 4.9 m × 4.9 m × 1.25 m grid scales. It is found that submesoscale frontogenesis differs from traditional frontogenesis theory due to four effects: Stokes forces, momentum and kinetic energy transfer from submesoscale eddies to frontal secondary circulations, resolved turbulent stresses, and unbalanced torque. In the energy, momentum, angular momentum, and vorticity budgets for the frontal overturning circulation, the Stokes shear force is a leading–order contributor, typically either the second or third largest source of frontal overturning. These effects violate hydrostatic and thermal wind balances during submesoscale frontogenesis. The effect of the Stokes shear force becomes stronger with increasing alignment of the front and Stokes shear and with a nondimensional scaling. The Stokes shear force and momentum transfer from submesoscale eddies significantly energize the frontal secondary circulation along with the buoyancy.

58 GEOSCIENCES↗

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗

Frontogenesis and frontal arrest of a dense filament in the oceanic surface boundary layer

The evolution of upper ocean currents involves a set of complex, poorly understood interactions between submesoscale turbulence (e.g. density fronts and filaments and coherent vortices) and smaller-scale boundary-layer turbulence. Here we simulate the lifecycle of a cold (dense) filament undergoing frontogenesis in the presence of turbulence generated by surface stress and/or buoyancy loss. This phenomenon is examined in large-eddy simulations with resolved turbulent motions in large horizontal domains using ${\sim}10^{10}$ grid points. Steady winds are oriented in directions perpendicular or parallel to the filament axis. Due to turbulent vertical momentum mixing, cold filaments generate a potent two-celled secondary circulation in the cross-filament plane that is frontogenetic, sharpens the cross-filament buoyancy and horizontal velocity gradients and blocks Ekman buoyancy flux across the cold filament core towards the warm filament edge. Within less than a day, the frontogenesis is arrested at a small width, ${\approx}100~\text{m}$ , primarily by an enhancement of the turbulence through a small submesoscale, horizontal shear instability of the sharpened filament, followed by a subsequent slow decay of the filament by further turbulent mixing. The boundary-layer turbulence is inhomogeneous and non-stationary in relation to the evolving submesoscale currents and density stratification. The occurrence of frontogenesis and arrest are qualitatively similar with varying stress direction or with convective cooling, but the detailed evolution and flow structure differ among the cases. Thus submesoscale filament frontogenesis caused by boundary-layer turbulence, frontal arrest by frontal instability and frontal decay by forward energy cascade, and turbulent mixing are generic processes in the upper ocean.

Mechanics↗

Langmuir turbulence and filament frontogenesis in the oceanic surface boundary layer

Submesoscale currents, small-scale turbulence and surface gravity waves co-exist in the upper ocean and interact in complex ways. To expose the couplings, the frontogenetic life cycle of an idealized cold dense submesoscale filament interacting with upper ocean Langmuir turbulence is investigated in large-eddy simulations (LESs) based on the incompressible wave-averaged equations. The simulations utilize large domains and fine meshes with $6.4\times 10^{9}$ grid points. Case studies are made with surface winds or surface cooling with waves oriented in across-filament (perpendicular) or down-filament (parallel) directions relative to the two-dimensional filament axis. The currents $u$ , $v$ and $w$ are aligned with the across-filament, down-filament and vertical directions, respectively. Frontogenesis is induced by across-filament Lagrangian secondary circulations in the boundary layer, and it is shown to be strongly impacted by surface waves, in particular the propagation direction relative to the filament axis. In a horizontally heterogeneous boundary layer, surface waves induce both mean and fluctuating Stokes-drift vortex forces that modify a linear, hydrostatic turbulent thermal wind (TTW) approximation for momentum. Down-filament winds and waves are found to be especially impactful, they significantly reduce the peak level of frontogenesis by fragmenting the filament into primary and secondary down-welling sites in a broad frontal zone over a width ${\sim}500~\text{m}$ . At peak frontogenesis, opposing down-filament jets $\langle v\rangle$ overlie each other resulting in a vigorous vertical shear layer $\unicode[STIX]{x2202}_{z}\langle v\rangle$ with large vertical momentum flux $\langle v^{\prime }w^{\prime }\rangle$ . Filament arrest is induced by a lateral shear instability that generates horizontal momentum flux $\langle u^{\prime }v^{\prime }\rangle$ at low wavenumbers. The turbulent vertical velocity patterns, indicative of coherent Langmuir cells, change markedly across the horizontal domain with both across-filament and down-filament winds under the action of submesoscale currents.

Mechanics↗

Intermittent Frontogenesis in the Alboran Sea

Abstract We present a phenomenological description and dynamical analysis of the Alboran fronts using a realistic simulation at submesoscale resolution. The study is focused on east Alboran fronts emerging within relatively strong flows that separate from the Spanish coast into the basin interior. Despite modest lateral shifting associated with the position of the Alboran anticyclonic gyres and variations in intensity, these fronts present a similar structure and dynamical configuration as the climatological Almeria–Oran front. The statistical analysis of our solution shows that strained-induced frontogenesis is a recurrent submesoscale mechanism associated with these fronts, and the process is assessed in terms of the advective Lagrangian frontogenetic tendencies associated with buoyancy and velocity horizontal gradients. Intermittency in their strength and patterns is indicative of high variability in the occurrence of active frontogenesis in association with the secondary (overturning) circulation across the frontal gradient. As a result, we find many episodes with strong surface fronts that do not have much associated downwelling. Frontogenesis and the associated secondary circulation are further explored during two particular frontal events, both showing strong downwelling of (1) cm s −1 extending down into the pycnocline. A frontogenetic contribution of turbulent vertical momentum mixing to the secondary circulation is identified in the easternmost region during the cold season, when the dynamics are strongly influenced by the intrusion of the salty Northern Current. The background vertical velocity fields observed during the analyzed events indicate other currents in the submesoscale range, including tidal and topographic internal waves.

Capó, Esther↗

Seasonal variation of the surface cross-shelf exchange in the northern South China Sea: a Lagrangian perspective

Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.

Hao, Dongliang↗

Quantification of Physical and Numerical Mixing in a Coastal Ocean Model Using Salinity Variance Budgets

Numerical mixing, the spurious mixing primarily generated by the discretization of advection, is often significant in estuarine and coastal models due to sharp, energetic fronts. In this study, we compare on- and offline estimates of numerical mixing in a submesoscale-resolving realistic simulation of the ocean state over the Texas-Louisiana continental shelf. While offline estimates of numerical mixing differ from online estimates, offline methods may be the only analysis available. This study offers insight into the differences between the on- and offline methods. We use two methods to estimate numerical mixing offline, based on salinity squared $s^2$ and volume-mean salinity variance $s^{\prime^2}$. Numerical mixing estimated from the $s^{\prime^2}$ budget is generally within 60\% of the magnitude for the online method but captures the temporal variability well. However, the $s^2$ budget compares poorly due to larger truncation errors associated with the tendency and advection terms, which can be reduced by increasing the model output frequency. We also investigate the effects of horizontal resolution on numerical mixing using a two-way nested grid. The volume-integrated numerical mixing constitutes 57\% of the bulk physical mixing -- the mixing prescribed by the turbulence closure scheme -- in the coarse model and may exceed the physical mixing by half an order of magnitude. We find that numerical mixing is reduced by 35\% on average in the nested model, less than expected based on scaling of the numerical mixing for an upwind advection scheme, likely due to new dynamical processes that emerge in the nested simulation.

54 ENVIRONMENTAL SCIENCES↗

Atmospheric Turbulent Intermittency Over the Arctic Sea‐Ice Surface During the MOSAiC Expedition

Turbulent motions in the Arctic stable boundary layer are characterized by intermittency, but they are rarely investigated due to limited observations, in particular over the sea‐ice surface. In the present study, we explore the characteristics of turbulent intermittency over the Arctic sea‐ice surface using data collected during the Multidisciplinary drifting Observation for the Study of Arctic Climate expedition from October 2019 to September 2020. We first develop a new algorithm, which performs well in identifying the spectral gap over the Arctic sea‐ice surface. Then the characteristics of intermittency are investigated. It is found that the strength of intermittency increases under the conditions of light surface wind speed, small surface wind speed gradient, and strong surface air temperature gradient. The momentum flux, sensible heat flux, and latent heat flux calculated by raw eddy‐covariance fluctuations are overestimated by 3%, 10%, and 24%, respectively, because submesoscale motions are included. Furthermore, the characteristics of the atmospheric boundary layer structure under various intermittency conditions reveal that strong low‐level jets are favorable to surface turbulent motions that result in weak intermittency, while strong temperature inversions above the surface layer suppress surface turbulent motions and lead to strong intermittency.

54 ENVIRONMENTAL SCIENCES↗

Pervasive iron limitation at subsurface chlorophyll maxima of the California Current

Subsurface chlorophyll maximum layers (SCMLs) are nearly ubiquitous in stratified water columns and exist at horizontal scales ranging from the submesoscale to the extent of oligotrophic gyres. These layers of heightened chlorophyll and/or phytoplankton concentrations are generally thought to be a consequence of a balance between light energy from above and a limiting nutrient flux from below, typically nitrate (NO 3 ). Here we present multiple lines of evidence demonstrating that iron (Fe) limits or with light colimits phytoplankton communities in SCMLs along a primary productivity gradient from coastal to oligotrophic offshore waters in the southern California Current ecosystem. SCML phytoplankton responded markedly to added Fe or Fe/light in experimental incubations and transcripts of diatom and picoeukaryote Fe stress genes were strikingly abundant in SCML metatranscriptomes. Using a biogeochemical proxy with data from a 40-y time series, we find that diatoms growing in California Current SCMLs are persistently Fe deficient during the spring and summer growing season. We also find that the spatial extent of Fe deficiency within California Current SCMLs has significantly increased over the last 25 y in line with a regional climate index. Finally, we show that diatom Fe deficiency may be common in the subsurface of major upwelling zones worldwide. Our results have important implications for our understanding of the biogeochemical consequences of marine SCML formation and maintenance.

54 ENVIRONMENTAL SCIENCES↗

Improving Process Level Understanding of Boundary Layer Winds over the Northeast U.S. Shelf: The Third Wind Forecast Improvement Project (WFIP3)

The third Wind Forecast Improvement Project (WFIP3), a U.S. Department of Energy and National Oceanic and Atmospheric Administration sponsored investigation, sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the marine atmospheric boundary layer. WFIP3 focused on mesoscale and submesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high-use coastal zone, using a 3D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the marine atmospheric boundary layer over the ocean was done from an air-sea interaction flux tower and extended deployments of a large autonomous barge platform. Numerous critical forecasting phenomena were observed that are being evaluated within regional coupled and uncoupled modeling systems, including the National Oceanic and Atmospheric Administration's foundational High-Resolution Rapid Refresh forecast model.

Kirincich, Anthony↗