We develop a parameterization for representing the effects of submesoscale symmetric instability (SI) in the ocean interior. SI may contribute to water mass modification and mesoscale energy dissipation in flow systems throughout the World Ocean. Dense gravity currents forced by surface buoyancy loss over shallow shelves are a particularly compelling test case, as they are characterized by density fronts and shears susceptible to a wide range of submesoscale instabilities. We present idealized experiments of Arctic shelf overflows employing the GFDL‐MOM6 in
Biological productivity in the Southern Ocean is limited by iron availability. Previous studies of iron supply have focused on mixed‐layer entrainment and diapycnal fluxes. However, the Southern Ocean is a region highly energetic mesoscale and submesoscale turbulence. Here we investigate the role of eddies in supplying iron to the euphotic zone, using a flat‐bottom zonally re‐entrant model, configured to represent the Antarctic Circumpolar Current region, that is coupled to a biogeochemical model with a realistic seasonal cycle. Eddies are admitted or suppressed by changing the model's horizontal resolution. We utilize cross spectral analysis and the generalized Omega equation to temporally and spatially decompose the vertical transport attributable to mesoscale and submesoscale motions. Our results suggest that the mesoscale vertical fluxes provide a first‐order pathway for transporting iron across the mixing‐layer base, where diapycnal mixing is weak, and must be included in modeling the open‐Southern Ocean iron budget.
more » « less- NSF-PAR ID:
- 10374570
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 11
- Issue:
- 12
- ISSN:
- 1942-2466
- Page Range / eLocation ID:
- p. 3934-3958
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract z * and isopycnal coordinates. At the highest resolutions, the dense flow undergoes geostrophic adjustment and forms bottom‐ and surface‐intensified jets. The density front along the topography combined with geostrophic shear initiates SI, leading to onset of secondary shear instability, dissipation of geostrophic energy, and turbulent mixing. We explore the impact of vertical coordinate, resolution, and parameterization of shear‐driven mixing on the representation of water mass transformation. We find that in isopycnal and low‐resolutionz * simulations, limited vertical resolution leads to inadequate representation of diapycnal mixing. This motivates our development of a parameterization for SI‐driven turbulence. The parameterization is based on identifying unstable regions through a balanced Richardson number criterion and slumping isopycnals toward a balanced state. The potential energy extracted from the large‐scale flow is assumed to correspond to the kinetic energy of SI which is dissipated through shear mixing. Parameterizing submesoscale instabilities by combining isopycnal slumping with diapycnal mixing becomes crucial as ocean models move toward resolving mesoscale eddies and fronts but not the submesoscale phenomena they host. -
more » « less
Abstract The primary productivity of the Southern Ocean ecosystem is limited by iron availability. Away from benthic and aeolian sources, iron reaches phytoplankton primarily when iron-rich subsurface waters enter the euphotic zone. Here, eddy-resolving physical/biogeochemical simulations of a seasonally-forced, open-Southern-Ocean ecosystem reveal that mesoscale and submesoscale isopycnal stirring effects a cross-mixed-layer-base transport of iron that sustains primary productivity. The eddy-driven iron supply and consequently productivity increase with model resolution. We show the eddy flux can be represented by specific well-tuned eddy parametrizations. Since eddy mixing rates are sensitive to wind forcing and large-scale hydrographic changes, these findings suggest a new mechanism for modulating the Southern Ocean biological pump on climate timescales.
-
null (Ed.)Submesoscale circulations influence momentum, buoyancy and transport of biological tracers and pollutants within the upper turbulent layer. How much and how far into the water column this influence extends remain open questions in most of the global ocean. This work evaluates the behavior of neutrally buoyant particles advected in simulations of the northern Gulf of Mexico by analyzing the trajectories of Lagrangian particles released multiple times at the ocean surface and below the mixed layer. The relative role of meso- and submesoscale dynamics is quantified by comparing results in submesoscale permitting and mesoscale resolving simulations. Submesoscale circulations are responsible for greater vertical transport across fixed depth ranges and also across the mixed layer, both into it and away from it, in all seasons. The significance of the submesoscale-induced transport, however, is far greater in winter. In this season, a kernel density estimation and a detailed vertical mixing analysis are performed. It is found that in the large mesoscale Loop Current eddy, upwelling into the mixed layer is the major contributor to the vertical fluxes, despite its clockwise circulation. This is opposite to the behavior simulated in the mesoscale resolving case. In the “submesoscale soup,” away from the large mesoscale structures such as the Loop Current and its detached eddies, upwelling into the mixed layer is distributed more uniformly than downwelling motions from the surface across the base of the mixed layer. Maps of vertical diffusivity indicate that there is an order of magnitude difference among simulations. In the submesoscale permitting case values are distributed around 10 –3 m 2 s –1 in the upper water column in winter, in agreement with recent indirect estimates off the Chilean coast. Diffusivities are greater in the eastern portion of the Gulf, where the submesoscale circulations are more intense due to sustained density gradients supplied by the warmer and saltier Loop Current.more » « less
-
more » « less
Abstract Global warming may modify submesoscale activity in the ocean through changes in the mixed layer depth (MLD) and lateral buoyancy gradients. As a case study we consider a region in the NE Atlantic under present and future climate conditions, using a time‐slice method and global and nested regional ocean models. The high resolution regional model reproduces the strong seasonal cycle in submesoscale activity observed under present‐day conditions. Focusing on the well‐resolved winter months, in the future, with a reduction in the MLD, there is a substantial reduction in submesoscale activity, an associated decrease in kinetic energy (KE) at the mesoscale, and the vertical buoyancy flux induced by submesoscale activity is reduced by a factor of 2. When submesoscale activity is suppressed, by increasing the parameterized lateral mixing in the model, the climate change induces a larger reduction in winter MLDs while there is less of a change in KE at the mesoscale. A scaling for the vertical buoyancy flux proposed by (Fox‐Kemper et al., 2008; doi:10.1175/2007JPO3792.1) based on the properties of mixed layer instability (MLI), is found to capture much of the seasonal and future changes to the flux in terms of regional averages as well as the spatial structure, although it over predicts the reduction in the flux in the winter months. The vertical buoyancy flux when the mixed layer is relatively shallow is significantly greater than that given by the scaling based on MLI, suggesting during these times other processes (besides MLI) may dominate submesoscale buoyancy fluxes.
-
more » « less
Abstract The mixing of tracers by mesoscale eddies, parameterized in many ocean general circulation models (OGCMs) as a diffusive‐advective process, contributes significantly to the distribution of tracers in the ocean. In the ocean interior, diffusive contribution occurs mostly along the direction parallel to local neutral density surfaces. However, near the surface of the ocean, small‐scale turbulence and the presence of the boundary itself break this constraint and the mesoscale transport occurs mostly along a plane parallel to the ocean surface (horizontal). Although this process is easily represented in OGCMs with geopotential vertical coordinates, the representation is more challenging in OGCMs that use a general vertical coordinate, where surfaces can be tilted with respect to the horizontal. We propose a method for representing the diffusive horizontal mesoscale fluxes within the surface boundary layer of general vertical coordinate OGCMs. The method relies on regridding/remapping techniques to represent tracers in a geopotential grid. Horizontal fluxes are calculated on this grid and then remapped back to the native grid, where fluxes are applied. The algorithm is implemented in an ocean model and tested in idealized and realistic settings. Horizontal diffusion can account for up to 10% of the total northward heat transport in the Southern Ocean and Western boundary current regions of the Northern Hemisphere. It also reduces the vertical stratification of the upper ocean, which results in an overall deepening of the surface boundary layer depth. Finally, enabling horizontal diffusion leads to meaningful reductions in the near‐surface global bias of potential temperature and salinity.
