This study investigates the variability of water mass transformation (WMT) within the Weddell Gyre (WG). The WG serves as a pivotal site for the meridional overturning circulation (MOC) and ocean ventilation because it is the primary origin of the largest volume of water mass in the global ocean, Antarctic Bottom Water (AABW). Recent mooring data suggest substantial seasonal and interannual variability of AABW properties exiting the WG, and studies have linked the variability to the large-scale climate forcings affecting wind stress in the WG region. However, the specific thermodynamic mechanisms that link variability in surface forcings to variability in water mass transformations and AABW export remain unclear. This study explores WMT variability via WMT volume budgets derived from Walin’s classic WMT framework, using three state-of-the-art, data-assimilating ocean reanalyses: Estimating the Circulation and Climate of the Ocean state estimate (ECCOv4), Southern Ocean State Estimate (SOSE) and Simple Ocean Data Assimilation (SODA). From the model outputs, we diagnose a closed form of the water mass budget for AABW that explicitly accounts for transport across the WG boundary, surface forcing, interior mixing, and numerical mixing. We examine the annual mean climatology of the WMT budget terms, the seasonal climatology, and finally the interannual variability. In ECCO and SOSE, we see strong interannual variability in AABW volume budget. In SOSE, we find an accelerating loss of AABW, driven largely by interior mixing and changes in surface salt fluxes. ECCO shows a similar trend during a 3-yr time period beyond what is covered in SOSE, but also reveals such trends to be part of interannual variability over a much longer time period. Overall, ECCO provides the most useful timeseries for understanding the processes and mechanisms that drive WMT and export variability. SODA, in contrast, displays unphysically large variability in AABW volume, which we attribute to its data assimilation scheme. We examine correlations between the WMT budgets and large-scale climate indices, including ENSO and SAM; no strong relationships emerge, suggesting that these reanalysis products may not reproduce the AABW export pathways and mechanisms hypothesized from observations.
more »
« less
Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models
Abstract Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under‐utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite‐Volume Generalized Vertical Coordinate (FV‐GVC) ocean models, including the identification of a non‐negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute‐force, comparing time‐averaged (“offline”) diagnostics on various vertical grids to timestep‐averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray‐enabled Python packages for evaluating WMT budgets in FV‐GVC models (culminating in the newxwmbpackage), which is intended to be model‐agnostic and available for community use and development.
more »
« less
- Award ID(s):
- 2332379
- PAR ID:
- 10626708
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 17
- Issue:
- 3
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Indo-Pacific Ocean appears exponentially stratified between 1- and 3-km depth with a decay scale on the order of 1 km. In his celebrated paper “Abyssal recipes,” W. Munk proposed a theoretical explanation of these observations by suggesting a pointwise buoyancy balance between the upwelling of cold water and the downward diffusion of heat. Assuming a constant upwelling velocity w and turbulent diffusivity κ , the model yields an exponential stratification whose decay scale is consistent with observations if κ ∼ 10 −4 m 2 s −1 . Over time, much effort has been made to reconcile Munk’s ideas with evidence of vertical variability in κ , but comparably little emphasis has been placed on the even stronger evidence that w decays toward the surface. In particular, the basin-averaged w nearly vanishes at 1-km depth in the Indo-Pacific. In light of this evidence, we consider a variable-coefficient, basin-averaged analog of Munk’s budget, which we verify against a hierarchy of numerical models ranging from an idealized basin-and-channel configuration to a coarse global ocean simulation. Study of the budget reveals that the decay of basin-averaged w requires a concurrent decay in basin-averaged κ to produce an exponential-like stratification. As such, the frequently cited value of 10 −4 m 2 s −1 is representative only of the bottom of the middepths, whereas κ must be much smaller above. The decay of mixing in the vertical is as important to the stratification as its magnitude . Significance Statement Using a combination of theory and numerical simulations, it is argued that the observed magnitude and shape of the global ocean stratification and overturning circulation appear to demand that turbulent mixing increases quasi-exponentially toward the ocean bottom. Climate models must therefore prescribe such a vertical profile of turbulent mixing in order to properly represent the heat and carbon uptake accomplished by the global overturning circulation on centennial and longer time scales.more » « less
-
Abstract The abyssal Southwest Pacific Basin has warmed significantly between 1992-2017, consistent with warming along the bottom limb of the meridional overturning circulation seen throughout the global oceans. Here we present a framework for assessing the abyssal heat budget that includes the time-dependent unsteady effects of decadal warming and direct and indirect estimates of diapycnal mixing from microscale temperature measurements and finescale parameterizations. The unsteady terms estimated from the decadalwarming rate are shown to be within a factor of 3 of the steady state terms in the abyssal heat budget for the coldest portion of the water column and therefore, cannot be ignored. We show that a reduction in the lateral heat flux for the coldest temperature classes compensated by an increase in warmer waters advected into the basin has important implications for the heat balance and diffusive heat fluxes in the basin. Finally, vertical diffusive heat fluxes are estimated in different ways: using the newly available CTD-mounted microscale temperature measurements, a finescale strain parameterization, and a vertical kinetic energy parameterization from data along the P06 transect along 32.5°S. The unsteady-state abyssal heat budget for the basin shows closure within error estimates, demonstrating that (i) unsteady terms have become consequential for the heat balance in the isotherms closest to the ocean bottom and (ii) direct and indirect estimates from full depth GO-SHIP hydrographic transects averaged over similarly large spatial and temporal scales can capture the basin-averaged abyssal mixing needed to close the deep overturning circulation.more » « less
-
Abstract In classic models of the tidally averaged gravitationally driven estuarine circulation, denser salty oceanic water moves up the estuary near the bottom, while less dense riverine water flows toward the ocean near the surface. Traditionally, it is assumed that the associated pressure gradient forces and salt advection are balanced by vertical mixing. This study, however, demonstrates that lateral (across the estuary width) transport processes are essential for maintaining the estuarine circulation. This is because for realistic estuarine bathymetry, the depth-integrated salt transport up the estuary is enhanced in the deeper estuary channel. A closed salt budget then requires the lateral transport of this excess salt in the deeper channel toward the estuarine flanks. To understand how such lateral transport affects the estuarine salt and momentum balances, we devise an idealized model with explicit lateral transport focusing on tidally averaged lateral mixing effects. Solutions for the along-estuary velocity and salinity are nondimensionalized to depend only on one single nondimensional parameter, referred to as the Fischer number, which describes the relative importance of lateral to vertical tidal mixing. For relatively strong lateral tidal mixing (greater Fischer number), salinity and velocity variations are predominantly vertical. For relatively weak lateral tidal mixing (smaller Fischer number), salinity and velocity variations are predominantly lateral. Overall, lateral transport greatly affects the estuarine circulation and controls the estuarine salinity intrusion length, which is demonstrated to scale inversely with the Fischer number.more » « less
-
Due to their limited resolution, numerical ocean models need to be interpreted as representing filtered or averaged equations. How to interpret models in terms of formally averaged equations, however, is not always clear, particularly in the case of hybrid or generalized vertical coordinate models, which limits our ability to interpret the model results and to develop parameterizations for the unresolved eddy contributions. We here derive the averaged hydrostatic Boussinesq equations in generalized vertical coordinates for an arbitrary thickness‐weighted average. We then consider various special cases and discuss the extent to which the averaged equations are consistent with existing ocean model formulations. As previously discussed, the momentum equations in existing depth‐coordinate models are best interpreted as representing Eulerian averages (i.e., averages taken at fixed depth), while the tracer equations can be interpreted as either Eulerian or thickness‐weighted isopycnal averages. Instead we find that no averaging is fully consistent with existing formulations of the parameterizations in semi‐Lagrangian discretizations of generalized vertical coordinate ocean models such as MOM6. A coordinate‐following average would require “coordinate‐aware” parameterizations that can account for the changing nature of the eddy terms as the coordinate changes. Alternatively, the model variables can be interpreted as representing either Eulerian or (thickness‐weighted) isopycnal averages, independent of the model coordinate that is being used for the numerical discretization. Existing parameterizations in generalized vertical coordinate models, however, are not always consistent with either of these interpretations, which, respectively, would require a three‐dimensional divergence‐free eddy tracer advection or a form‐stress parameterization in the momentum equations.more » « less
An official website of the United States government

