David Marshall

Abstract:

Particle Image Velocimetry (PIV) measurement data and numerical simulation data for the paper "Energetic Constraints on Baroclinic Eddy Heat Transport with a Beta Effect in the Laboratory".

The experimental data (dataGCM.zip) used the MITgcm (Marshall et al., 1997) with the input adapted from its rotating tank tutorial (model version checkpoint67c, see the manual from Adcroft et al., 2018). The sample is also for the flow regime at 1 rad/s rotation rate and 4 K temperature contrast, including the beta-plane and f-plane cases. The dataset includes temperature and velocity fields, etc. In the cylindrical coordinate, there are 120 zonal grids, 180 radial grids, and 486 vertical grids for the beta-plane case. This data set is limited in the zonal resolution. This limitation should be considered to validate results derived from the simulation data.

The experimental data (dataPIV.zip) comprises velocity fields from the multi-level PIV measurements, see the method of Wordsworth et al., 2008. The data sample is for the flow regime with the experimental control of 1 rad/s rotation rate and 4 K temperature contrast. Both samples of the beta-plane and f-plane cases are included. The software package UVMAT (Sommeria, 2013) had been applied to process PIV images in time series. For further processing, data should be re-scaled in both the spatial and time dimensions for the annulus diameter of 16 cm and the sampling time step of 0.33 seconds. The multi-level information is indexed for files with 1 to 5 at the lower height level 1, 6 to 10 at the lower level 2, 11 to 15 for middle level 3, etc. Indexes in the multiples of 5 are velocities derived from image pairs between distinct height levels and should not be included for further processing, e.g. indexes 5 or 10. Indexes for upper height level 5 in the beta-plane case should not be directly used for further processing because the sloping topography interferes with the sampling horizontal plane. This data set is limited by the 6-second delay in the sampling between the discrete height levels. This limitation should be considered to validate results derived from the experimental data.

References:

Adcroft, A., Campin, J.-M., Dutkiewicz, S., Evangelinos, C., Ferreira, D., Forget, G., . . . Molod, A. (2018). MITgcm user manual. (https://dspace.mit.edu/handle/1721.1/117188, https://mitgcm.readthedocs.io/en/latest/, [Last Accessed: 25 September 2023])

Marshall, J., Adcroft, A., Hill, C., Perelman, L., & Heisey, C. (1997). A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers. Journal of Geophysical Research: Oceans, 102 (C3), 5753-5766. doi: https://doi.org/10.1029/96JC02775

Sommeria, J. (2013). UVMAT. (http://servforge.legi.grenoble-inp.fr/projects/soft-uvmat, [Last Accessed: 29 September 2023])

Wordsworth, R. D., Read, P. L., & Yamazaki, Y. H. (2008). Turbulence, waves, and jets in a differentially heated rotating annulus experiment. Physics of Fluids, 20 (12), 126602.

Abstract:

There is an increasing interest in mesoscale eddy parameterizations that are scale-aware, normally interpreted to mean that a parameterization does not require parameter recalibration as the model resolution changes. Here we examine whether Gent–McWilliams (GM) based version of GEOMETRIC, a mesoscale eddy parameterization that is constrained by a parameterized eddy energy budget, is scale-aware in its energetics. It is generally known that GM-based schemes severely damp out explicit eddies, so the parameterized component would be expected to dominate across resolutions, and we might expect a negative answer to the question of energetic scale-awareness. A consideration of why GM-based schemes damp out explicit eddies leads a suggestion for what we term a splitting procedure: a definition of a “large-scale” field is sought, and the eddy-induced velocity from the GM-scheme is computed from and acts only on the large-scale field, allowing explicit and parameterized components to co-exist. Within the context of an idealized re-entrant channel model of the Southern Ocean, evidence is provided that the GM-based version of GEOMETRIC is scale-aware in the energetics as long as we employ a splitting procedure. The splitting procedure also leads to an improved representation of mean states without detrimental effects on the explicit eddy motions.

Abstract:

Estimates of the Weddell Gyre transport vary widely between climate simulations. Here, we investigate if inter-model variability can originate from differences in the horizontal resolution of the ocean model. We run an idealized model of the Weddell Gyre at eddy-parameterized, eddy-permitting, and eddy-rich resolutions and find that the gyre is strongly sensitive to horizontal resolution. The gyre transport is largest at eddy-permitting resolutions (45 Sv with a noisy bathymetry) and smallest at eddy-parameterized resolutions (12 Sv). The eddy-permitting simulations have the largest horizontal density gradients and the weakest stratification over the gyre basin. The large horizontal density gradients induce a significant thermal wind transport and increase the mean available potential energy for mesoscale eddies. The distribution of eddy kinetic energy indicates that explicit eddies in simulations intensify the bottom circulation of the gyre via non-linear dynamics. If climate models adopt horizontal resolutions that the Weddell Gyre is most sensitive to, then simulations of the Weddell Gyre could become more disparate.

Abstract:

Although the surface eddy kinetic energy (EKE) has been well studied using satellite altimeter and surface drifter observations, our knowledge of EKE in the ocean interior is much more limited due to the sparsity of subsurface current measurements. Here we develop a new approach for estimating EKE over the full depth of the global ocean by combining 20 years of satellite altimeter and Argo float data to infer the vertical profile of eddies. The inferred eddy profiles are surface-intensified at low latitudes and deep-reaching at mid- and high latitudes. They compare favorably to the first empirical orthogonal function obtained from current meter velocities. The global-integrated EKE estimated from the inferred profiles is about 3.1 × 10¹⁸ J, which is close to that estimated from the surface mode (3.0 × 10¹⁸ J) but about 30% smaller than that estimated from the traditional flat bottom modes (4.6 × 10¹⁸ J).