• Allen, J. S., and P. A. Newberger, 1998: On symmetric instabilities in oceanic bottom boundary layers. J. Phys. Oceanogr., 28, 11311151, https://doi.org/10.1175/1520-0485(1998)028<1131:OSIIOB>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Armi, L., and R. C. Millard, 1976: The bottom boundary layer of the deep ocean. J. Geophys. Res., 81, 49834990, https://doi.org/10.1029/JC081i027p04983.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Benthuysen, J., and L. N. Thomas, 2012: Friction and diapycnal mixing at a slope: Boundary control of potential vorticity. J. Phys. Oceanogr., 42, 15091523, https://doi.org/10.1175/JPO-D-11-0130.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Benthuysen, J., L. N. Thomas, and S. J. Lentz, 2015: Rapid generation of upwelling at a shelf break caused by buoyancy shutdown. J. Phys. Oceanogr., 45, 294312, https://doi.org/10.1175/JPO-D-14-0104.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blumsack, S. L., and P. J. Gierasch, 1972: Mars: The effects of topography on baroclinic instability. J. Atmos. Sci., 29, 10811089, https://doi.org/10.1175/1520-0469(1972)029<1081:MTEOTO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Boccaletti, G., R. Ferrari, and B. Fox-Kemper, 2007: Mixed layer instabilities and restratification. J. Phys. Oceanogr., 37, 22282250, https://doi.org/10.1175/JPO3101.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., 2012: Baroclinic instability of an idealized tidal mixing front. J. Mar. Res., 70, 661688, https://doi.org/10.1357/002224012805262716.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., 2013: Instability of a tidal mixing front in the presence of realistic tides and mixing. J. Mar. Res., 71, 227251, https://doi.org/10.1357/002224013807719473.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., 2016: Continental shelf baroclinic instability. Part I: Relaxation from upwelling or downwelling. J. Phys. Oceanogr., 46, 551568, https://doi.org/10.1175/JPO-D-15-0047.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., and D. A. Cherian, 2013: Instability of an idealized tidal mixing front: Symmetric instabilities and frictional effects. J. Mar. Res., 71, 425450, https://doi.org/10.1357/002224013812587582.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., and H. Seo, 2016: Continental shelf baroclinic instability. Part II: Oscillating wind forcing. J. Phys. Oceanogr., 46, 569582, https://doi.org/10.1175/JPO-D-15-0048.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Burns, K. J., G. M. Vasil, J. S. Oishi, D. Lecoanet, and B. Brown, 2016: Dedalus: Flexible framework for spectrally solving differential equations. Astrophysics Source Code Library, http://ascl.net/phpBB3/viewtopic.php?t=33963.

  • Callies, J., 2018: Restratification of abyssal mixing layers by submesoscale baroclinic eddies. J. Phys. Oceanogr., 48, 19952010, https://doi.org/10.1175/JPO-D-18-0082.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., and R. Ferrari, 2018: Baroclinic instability in the presence of convection. J. Phys. Oceanogr., 48, 4560, https://doi.org/10.1175/JPO-D-17-0028.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., R. Ferrari, J. M. Klymak, and J. Gula, 2015: Seasonality in submesoscale turbulence. Nat. Commun., 6, 6862, https://doi.org/10.1038/ncomms7862.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., G. Flierl, R. Ferrari, and B. Fox-Kemper, 2016: The role of mixed-layer instabilities in submesoscale turbulence. J. Fluid Mech., 788, 541, https://doi.org/10.1017/jfm.2015.700.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Costello, M. J., A. Cheung, and N. De Hauwere, 2010: Surface area and the seabed area, volume, depth, slope, and topographic variation for the world’s seas, oceans, and countries. Environ. Sci. Technol., 44, 88218828, https://doi.org/10.1021/es1012752.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • de Szoeke, R. A., 1983: Baroclinic instability over wavy topography. J. Fluid Mech., 130, 279298, https://doi.org/10.1017/S0022112083001093.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Dewar, W. K., J. C. McWilliams, and M. J. Molemaker, 2015: Centrifugal instability and mixing in the California Undercurrent. J. Phys. Oceanogr., 45, 12241241, https://doi.org/10.1175/JPO-D-13-0269.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Eady, E. T., 1949: Long waves and cyclone waves. Tellus, 1, 3352, https://doi.org/10.3402/tellusa.v1i3.8507.

  • Ferrari, R., A. Mashayek, T. J. McDougall, M. Nikurashin, and J.-M. Campin, 2016: Turning ocean mixing upside down. J. Phys. Oceanogr., 46, 22392261, https://doi.org/10.1175/JPO-D-15-0244.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fox-Kemper, B., R. Ferrari, and R. Hallberg, 2008: Parameterization of mixed layer eddies. Part I: Theory and diagnosis. J. Phys. Oceanogr., 38, 11451165, https://doi.org/10.1175/2007JPO3792.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Garrett, C., 1991: Marginal mixing theories. Atmos.–Ocean, 29, 313339, https://doi.org/10.1080/07055900.1991.9649407.

  • Garrett, C., P. MacCready, and P. Rhines, 1993: Boundary mixing and arrested Ekman layers: Rotating stratified flow near a sloping boundary. Annu. Rev. Fluid Mech., 25, 291323, https://doi.org/10.1146/annurev.fl.25.010193.001451.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gula, J., M. J. Molemaker, and J. C. McWilliams, 2015: Topographic vorticity generation, submesoscale instability and vortex street formation in the Gulf Stream: Topographic vorticity generation. Geophys. Res. Lett., 42, 40544062, https://doi.org/10.1002/2015GL063731.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gula, J., M. J. Molemaker, and J. C. McWilliams, 2016: Topographic generation of submesoscale centrifugal instability and energy dissipation. Nat. Commun., 7, 12811, https://doi.org/10.1038/ncomms12811.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Haine, T. W. N., and J. Marshall, 1998: Gravitational, symmetric, and baroclinic instability of the ocean mixed layer. J. Phys. Oceanogr., 28, 634658, https://doi.org/10.1175/1520-0485(1998)028<0634:GSABIO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hetland, R. D., 2017: Suppression of baroclinic instabilities in buoyancy-driven flow over sloping bathymetry. J. Phys. Oceanogr., 47, 4968, https://doi.org/10.1175/JPO-D-15-0240.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hoskins, B. J., 1974: The role of potential vorticity in symmetric stability and instability. Quart. J. Roy. Meteor. Soc., 100, 480482, https://doi.org/10.1002/qj.49710042520.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Isachsen, P. E., 2011: Baroclinic instability and eddy tracer transport across sloping bottom topography: How well does a modified Eady model do in primitive equation simulations? Ocean Modell., 39, 183199, https://doi.org/10.1016/j.ocemod.2010.09.007.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Johnson, L., C. M. Lee, and E. A. D’Asaro, 2016: Global estimates of lateral springtime restratification. J. Phys. Oceanogr., 46, 15551573, https://doi.org/10.1175/JPO-D-15-0163.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ledwell, J. R., E. T. Montgomery, K. L. Polzin, L. C. St. Laurent, R. W. Schmitt, and J. M. Toole, 2000: Evidence for enhanced mixing over rough topography in the abyssal ocean. Nature, 403, 179182, https://doi.org/10.1038/35003164.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lévy, M., P. Klein, A.-M. Tréguier, D. Iovino, G. Madec, S. Masson, and K. Takahashi, 2010: Modifications of gyre circulation by sub-mesoscale physics. Ocean Modell., 34, 115, https://doi.org/10.1016/j.ocemod.2010.04.001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • MacCready, P., and P. B. Rhines, 1991: Buoyant inhibition of Ekman transport on a slope and its effect on stratified spin-up. J. Fluid Mech., 223, 631661, https://doi.org/10.1017/S0022112091001581.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • MacCready, P., and P. B. Rhines, 1993: Slippery bottom boundary layers on a slope. J. Phys. Oceanogr., 23, 522, https://doi.org/10.1175/1520-0485(1993)023<0005:SBBLOA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mahadevan, A., and A. Tandon, 2006: An analysis of mechanisms for submesoscale vertical motion at ocean fronts. Ocean Modell., 14, 241256, https://doi.org/10.1016/j.ocemod.2006.05.006.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mahadevan, A., E. D’Asaro, C. Lee, and M. J. Perry, 2012: Eddy-driven stratification initiates North Atlantic spring phytoplankton blooms. Science, 337, 5458, https://doi.org/10.1126/science.1218740.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McDougall, T. J., and R. Ferrari, 2017: Abyssal upwelling and downwelling driven by near-boundary mixing. J. Phys. Oceanogr., 47, 261283, https://doi.org/10.1175/JPO-D-16-0082.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., 2016: Submesoscale currents in the ocean. Proc. Roy. Soc. London, A472, 20160117, https://doi.org/10.1098/rspa.2016.0117.

    • Search Google Scholar
    • Export Citation
  • Mechoso, C. R., 1980: Baroclinic instability of flows along sloping boundaries. J. Atmos. Sci., 37, 13931399, https://doi.org/10.1175/1520-0469(1980)037<1393:BIOFAS>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Molemaker, M. J., J. C. McWilliams, and W. K. Dewar, 2015: Submesoscale instability and generation of mesoscale anticyclones near a separation of the California Undercurrent. J. Phys. Oceanogr., 45, 613629, https://doi.org/10.1175/JPO-D-13-0225.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Moum, J. N., A. Perlin, J. M. Klymak, M. D. Levine, T. Boyd, and P. M. Kosro, 2004: Convectively driven mixing in the bottom boundary layer. J. Phys. Oceanogr., 34, 21892202, https://doi.org/10.1175/1520-0485(2004)034<2189:CDMITB>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nakamura, N., 1988: Scale selection of baroclinic instability—Effects of stratification and nongeostrophy. J. Atmos. Sci., 45, 32533268, https://doi.org/10.1175/1520-0469(1988)045<3253:SSOBIO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nikurashin, M., and R. Ferrari, 2013: Overturning circulation driven by breaking internal waves in the deep ocean. Geophys. Res. Lett., 40, 31333137, https://doi.org/10.1002/grl.50542.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Pedlosky, J., 1979: Geophysical Fluid Dynamics. Springer, 710 pp.

    • Crossref
    • Export Citation
  • Pedlosky, J., 2016: Baroclinic instability over topography: Unstable at any wave number. J. Mar. Res., 74, 119, https://doi.org/10.1357/002224016818377595.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Polzin, K. L., 1997: Spatial variability of turbulent mixing in the abyssal ocean. Science, 276, 9396, https://doi.org/10.1126/science.276.5309.93.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rhines, P. B., 1975: Waves and turbulence on a beta-plane. J. Fluid Mech., 69, 417443, https://doi.org/10.1017/S0022112075001504.

  • Ripa, P., 1999: On the validity of layered models of ocean dynamics and thermodynamics with reduced vertical resolution. Dyn. Atmos. Oceans, 29, 140, https://doi.org/10.1016/S0377-0265(98)00056-6.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ripa, P., 2001: Waves and resonance in free-boundary baroclinic instability. J. Fluid Mech., 428, 387408, https://doi.org/10.1017/S0022112000002731.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ruan, X., A. F. Thompson, M. M. Flexas, and J. Sprintall, 2017: Contribution of topographically generated submesoscale turbulence to Southern Ocean overturning. Nat. Geosci., 10, 840845, https://doi.org/10.1038/ngeo3053.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Sasaki, H., P. Klein, B. Qiu, and Y. Sasai, 2014: Impact of oceanic-scale interactions on the seasonal modulation of ocean dynamics by the atmosphere. Nat. Commun., 5, 5636, https://doi.org/10.1038/ncomms6636.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Solodoch, A., A. L. Stewart, and J. C. McWilliams, 2016: Baroclinic instability of axially symmetric flow over sloping bathymetry. J. Fluid Mech., 799, 265296, https://doi.org/10.1017/jfm.2016.376.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Spall, M. A., 2004: Boundary currents and watermass transformation in marginal seas. J. Phys. Oceanogr., 34, 11971213, https://doi.org/10.1175/1520-0485(2004)034<1197:BCAWTI>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • St. Laurent, L. C., J. M. Toole, and R. W. Schmitt, 2001: Buoyancy forcing by turbulence above rough topography in the abyssal Brazil basin. J. Phys. Oceanogr., 31, 34763495, https://doi.org/10.1175/1520-0485(2001)031<3476:BFBTAR>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stahr, F. R., and T. B. Sanford, 1999: Transport and bottom boundary layer observations of the North Atlantic Deep Western Boundary Current at the Blake Outer Ridge. Deep-Sea Res. II, 46, 205243, https://doi.org/10.1016/S0967-0645(98)00101-5.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stone, P. H., 1966: On non-geostrophic baroclinic stability. J. Atmos. Sci., 23, 390400, https://doi.org/10.1175/1520-0469(1966)023<0390:ONGBS>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stone, P. H., 1971: Baroclinic stability under non-hydrostatic conditions. J. Fluid Mech., 45, 659671, https://doi.org/10.1017/S0022112071000260.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Su, Z., J. Wang, P. Klein, A. F. Thompson, and D. Menemenlis, 2018: Ocean submesoscales as a key component of the global heat budget. Nat. Commun., 9, 775, https://doi.org/10.1038/s41467-018-02983-w.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taylor, J. R., 2016: Turbulent mixing, restratification, and phytoplankton growth at a submesoscale eddy: Submesoscales, turbulence, and phytoplankton growth. Geophys. Res. Lett., 43, 57845792, https://doi.org/10.1002/2016GL069106.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taylor, J. R., and R. Ferrari, 2009: On the equilibration of a symmetrically unstable front via a secondary shear instability. J. Fluid Mech., 622, 103113, https://doi.org/10.1017/S0022112008005272.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taylor, J. R., and R. Ferrari, 2010: Buoyancy and wind-driven convection at mixed layer density fronts. J. Phys. Oceanogr., 40, 12221242, https://doi.org/10.1175/2010JPO4365.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., 2005: Destruction of potential vorticity by winds. J. Phys. Oceanogr., 35, 24572466, https://doi.org/10.1175/JPO2830.1.

  • Thomas, L. N., and P. B. Rhines, 2002: Nonlinear stratified spin-up. J. Fluid Mech., 473, 211244, https://doi.org/10.1017/S0022112002002367.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., and J. R. Taylor, 2010: Reduction of the usable wind-work on the general circulation by forced symmetric instability. Geophys. Res. Lett., 37, L18606, https://doi.org/10.1029/2010GL044680.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., A. Tandon, and A. Mahadevan, 2008: Submesoscale processes and dynamics. Ocean Modeling in an Eddying Regime, Geophys. Monogr., Vol. 177, Amer. Geophys. Union, 17–38, https://doi.org/10.1029/177GM04.

    • Crossref
    • Export Citation
  • Thomas, L. N., J. R. Taylor, R. Ferrari, and T. M. Joyce, 2013: Symmetric instability in the Gulf Stream. Deep-Sea Res. II, 91, 96110, https://doi.org/10.1016/j.dsr2.2013.02.025.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Toole, J., R. Curry, T. Joyce, M. McCartney, and B. Peña-Molino, 2011: Transport of the North Atlantic Deep Western Boundary Current about 39°N, 70°W: 2004–2008. Deep-Sea Res. II, 58, 17681780, https://doi.org/10.1016/j.dsr2.2010.10.058.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Umlauf, L., and H. Burchard, 2011: Diapycnal transport and mixing efficiency in stratified boundary layers near sloping topography. J. Phys. Oceanogr., 41, 329345, https://doi.org/10.1175/2010JPO4438.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Umlauf, L., H. Burchard, and K. Bolding, 2005: GOTM—Scientific documentation, version 3.2. Marine Science Rep. 63, Institut für Ostseeforschung Warnemünde, 274 pp., https://www.io-warnemuende.de/tl_files/forschung/meereswissenschaftliche-berichte/mebe63_2005-gotm.pdf.

  • Umlauf, L., W. D. Smyth, and J. N. Moum, 2015: Energetics of bottom Ekman layers during buoyancy arrest. J. Phys. Oceanogr., 45, 30993117, https://doi.org/10.1175/JPO-D-15-0041.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Vanneste, J., 1993: The Kelvin-Helmholtz instability in a non-geostrophic baroclinic unstable flow. Math. Comput. Model., 17, 149154, https://doi.org/10.1016/0895-7177(93)90099-K.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wenegrat, J. O., and M. J. McPhaden, 2016: Wind, waves, and fronts: Frictional effects in a generalized Ekman model. J. Phys. Oceanogr., 46, 371394, https://doi.org/10.1175/JPO-D-15-0162.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wenegrat, J. O., L. N. Thomas, J. Gula, and J. C. McWilliams, 2018: Effects of the submesoscale on the potential vorticity budget of ocean mode waters. J. Phys. Oceanogr., 48, 21412165, https://doi.org/10.1175/JPO-D-17-0219.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • View in gallery

    (a) Schematic of a BBL over a linear topographic slope, with solid contours indicating buoyancy surfaces and the dashed line indicating the top of a weakly stratified BBL. The rotated coordinate system is indicated in the lower left. (b) Simplified domain used in the basic parameter space exploration (section 2c). Throughout, the background velocity is assumed to be a function only of the slope-normal coordinate z. The relationships between the slope-normal derivative and derivatives in the nonrotated frame are also indicated.

  • View in gallery

    Example idealized vertical structure of a stratified interior overlying a weakly stratified BBL.

  • View in gallery

    Growth rates, normalized by f −1, for the idealized profile shown in Fig. 2 for various slope angles (legend). The low-wavenumber, full-depth baroclinic modes are indicated by the dashed box. Growth rates of the low-wavenumber instability are highly modulated by the slope angle, whereas the high wavenumber features are robust to changing slope angle.

  • View in gallery

    Structure of the fastest-growing perturbation (l ≈ 6 × 10−4 m−1) for θ = 5 × 10−3. Velocities are normalized relative to the maximum across-slope velocity and buoyancy relative to the maximum buoyancy perturbation. Kinetic energy tendency terms are defined in appendix A, with buoyancy production given by VBP and shear production by the sum of the lateral and vertical shear production terms (LSP + VSP).

  • View in gallery

    Baroclinic instability growth rates from (18) as a function of slope parameter α and wavenumber normalized by the deformation radius NH/f . Growth rates are normalized by f−1, and it is assumed that Ri = 100. Shown are (left) a horizontal upper boundary and (right) an upper boundary that matches the topographic slope, as shown schematically for α > 0, and α < 0 in each plot.

  • View in gallery

    Growth rates for Ri = (left) 10 and (right) 1, as a function of slope parameter α and wavenumber normalized by the deformation radius NH/f . Growth rates are normalized by f−1, and Λ = 5f.

  • View in gallery

    Growth rates for a fixed topographic slope angle (θ = 5 × 10−3) and vertical shear Λ = 10−3 s−1, with varying N2. Wavenumbers are normalized by ΛH/f, growth rates by f−1.

  • View in gallery

    Solution to (19)(21) with parameters as given in (25) and (26). In the upper-left panel the along-slope velocity V is shown in orange and the across-slope velocity U is shown in blue.

  • View in gallery

    Growth rate of the baroclinic mode for the base state shown in Fig. 8.

  • View in gallery

    Structure of the fastest-growing baroclinic mode for the base state shown in Fig. 8, discussed in section 3. Velocities are normalized relative to the maximum across-slope velocity and buoyancy relative to the maximum buoyancy perturbation. Kinetic energy tendency terms are defined in appendix A, with buoyancy production given by VBP, shear production by the sum of the lateral and vertical shear production terms (LSP + VSP), and dissipation by DKE.

  • View in gallery

    BBL response to a downwelling-favorable interior flow, = 0.1 m s−1, with f = 10−4 s−1, = 3.5 × 10−3 s−2, and θ = 10−2, as in MacCready and Rhines (1993). Viscosity ν follows κ closely and hence is not shown here.

  • View in gallery

    Hovmöller plot of the baroclinic growth rates as a function of wavenumber and time, normalized by the inertial period, for the case of downwelling favorable interior flow (Fig. 11). Growth rates are normalized by f−1.

  • View in gallery

    Structure of the fastest-growing baroclinic mode for the downwelling-favorable Ekman layer (Fig. 11) at t = 10 inertial periods, discussed in section 4. Velocities are normalized relative to the maximum across-slope velocity, buoyancy relative to the maximum buoyancy perturbation. Kinetic energy tendency terms are defined in appendix A, with buoyancy production given by VBP, shear production by the sum of the lateral and vertical shear production terms (LSP + VSP), and dissipation by DKE.

  • View in gallery

    As in Fig. 11, but for upwelling favorable interior flow, = −0.1 m s−1.

  • View in gallery

    As in Fig. 12, but for the upwelling favorable case shown in Fig. 14.

  • View in gallery

    Initial conditions for along-slope velocity (color) and buoyancy (gray contours) for the nonlinear simulation discussed in section 5 and appendix B.

  • View in gallery

    Snapshots of (top) buoyancy and (bottom) slope-normal velocity (w) evaluated at 50 m above the bottom.

  • View in gallery

    (top) Domain-integrated EKE compared to the prediction of linear theory (calculated as discussed in section 5) and (bottom) domain-averaged EKE tendency (bottom) with components as indicated in the legends.

  • View in gallery

    Example regime diagram for baroclinic instability (BI), symmetric instability (SI), centrifugal instability (CI), and mixed instability (SI/CI) in the ocean BBL. The red line indicates q = 0, the green line indicates , and the dotted black lines indicate α = (1, 2)α. The dashed blue line is the theoretical expectation for LSP/VSP = 1, given by (33), and the solid blue line is calculated using numerical solutions of (6)(10) with δRi−1/2 = 0.1, l = 0, and k = 20f/(NH). See section 6 for further explanation.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 538 538 23
PDF Downloads 358 358 21

Submesoscale Baroclinic Instability in the Bottom Boundary Layer

View More View Less
  • 1 Earth System Science, Stanford University, Stanford, California
  • | 2 Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, California
  • | 3 Earth System Science, Stanford University, Stanford, California
© Get Permissions
Full access

Abstract

Weakly stratified layers over sloping topography can support a submesoscale baroclinic instability mode, a bottom boundary layer counterpart to surface mixed layer instabilities. The instability results from the release of available potential energy, which can be generated because of the observed bottom intensification of turbulent mixing in the deep ocean, or the Ekman adjustment of a current on a slope. Linear stability analysis suggests that the growth rates of bottom boundary layer baroclinic instabilities can be comparable to those of the surface mixed layer mode and are relatively insensitive to topographic slope angle, implying the instability is robust and potentially active in many areas of the global oceans. The solutions of two separate one-dimensional theories of the bottom boundary layer are both demonstrated to be linearly unstable to baroclinic instability, and results from an example nonlinear simulation are shown. Implications of these findings for understanding bottom boundary layer dynamics and processes are discussed.

© 2018 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Jacob O. Wenegrat, jwenegrat@stanford.edu

Abstract

Weakly stratified layers over sloping topography can support a submesoscale baroclinic instability mode, a bottom boundary layer counterpart to surface mixed layer instabilities. The instability results from the release of available potential energy, which can be generated because of the observed bottom intensification of turbulent mixing in the deep ocean, or the Ekman adjustment of a current on a slope. Linear stability analysis suggests that the growth rates of bottom boundary layer baroclinic instabilities can be comparable to those of the surface mixed layer mode and are relatively insensitive to topographic slope angle, implying the instability is robust and potentially active in many areas of the global oceans. The solutions of two separate one-dimensional theories of the bottom boundary layer are both demonstrated to be linearly unstable to baroclinic instability, and results from an example nonlinear simulation are shown. Implications of these findings for understanding bottom boundary layer dynamics and processes are discussed.

© 2018 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Jacob O. Wenegrat, jwenegrat@stanford.edu

1. Introduction

An important development in our understanding of the upper ocean, occurring largely over the last decade, has been the recognition that submesoscale (horizontal scales of approximately 0.1–10 km) fronts, eddies, and instabilities are a prominent component of the dynamics of the ocean surface mixed layer (Boccaletti et al. 2007; Callies et al. 2015). These submesoscale processes alter the classic one-dimensional picture of boundary layer dynamics, giving rise to a host of new physical processes that have been shown to affect both large-scale ocean processes (Lévy et al. 2010; Wenegrat et al. 2018) and the turbulence properties of the boundary layer itself (Taylor and Ferrari 2010; Thomas and Taylor 2010; Taylor 2016). However, despite the large body of literature that has developed on submesoscale processes in the surface boundary layer [as reviewed in Thomas et al. (2008) and McWilliams (2016)], these types of processes have received much less attention in the bottom boundary layer (BBL),1 despite BBLs over sloping topography exhibiting key similarities with surface mixed layers at a front, particularly the existence of available potential energy in the form of a horizontal buoyancy gradient (Fig. 1).

Fig. 1.
Fig. 1.

(a) Schematic of a BBL over a linear topographic slope, with solid contours indicating buoyancy surfaces and the dashed line indicating the top of a weakly stratified BBL. The rotated coordinate system is indicated in the lower left. (b) Simplified domain used in the basic parameter space exploration (section 2c). Throughout, the background velocity is assumed to be a function only of the slope-normal coordinate z. The relationships between the slope-normal derivative and derivatives in the nonrotated frame are also indicated.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Recently, however, observations and numerical modeling have begun to suggest that the BBL does indeed support an active submesoscale turbulence field, with dynamical implications potentially as broad as for submesoscale processes in the surface boundary layer. Submesoscale processes in BBLs along topography have been shown to enhance cross-shelf exchanges (Gula et al. 2015), contribute to interior water mass transformation (Ruan et al. 2017), generate long-lived submesoscale coherent vortices (Molemaker et al. 2015), and provide another route for cascading energy from the large-scale flow into unbalanced motions and eventual dissipation (Gula et al. 2016). Particular focus has been given to centrifugal (CI) and symmetric instabilities (SI) in the BBL, both of which are associated with negative potential vorticity (Hoskins 1974; Haine and Marshall 1998), which can occur in BBLs over steep topography (Dewar et al. 2015) or in the presence of interior flows driving downslope Ekman transport in the boundary layer (Allen and Newberger 1998). The role of baroclinic instability in the BBL is currently less clear, despite the fact that surface mixed layer baroclinic instabilities are one of the more thoroughly studied aspects of submesoscale dynamics, with demonstrated effects on boundary layer restratification (Boccaletti et al. 2007; Fox-Kemper et al. 2008), surface potential vorticity fluxes (Wenegrat et al. 2018), boundary layer turbulence (Taylor 2016), biological productivity (Mahadevan et al. 2012), and the mesoscale eddy field (Sasaki et al. 2014). Understanding the conditions under which baroclinic instability can be expected to be active in the BBL is therefore a key step toward understanding submesoscale turbulence in the BBL.

In this article, we extend a variety of earlier works on baroclinic instability over topography (e.g., Blumsack and Gierasch 1972; Mechoso 1980; Pedlosky 2016; Solodoch et al. 2016) to the nongeostrophic regime appropriate for the low Richardson numbers typical of the BBL. Using a linear stability analysis, we explore the parameter dependence of baroclinic instabilities in a weakly stratified lower layer over sloping topography, which can be considered as a bottom counterpart to the surface mixed layer instability (Boccaletti et al. 2007). The submesoscale BBL mode is shown to be more robust to topographic slope than mesoscale baroclinic instability, and we demonstrate that the solutions to two classic one-dimensional theories of the BBL structure are unstable to BBL baroclinic instability. An example idealized nonlinear simulation further suggests that at finite amplitude the BBL baroclinic instability mode can generate vertical buoyancy fluxes, and vertical velocities, similar in magnitude to those associated with surface mixed layer instabilities. It is therefore expected that submesoscale baroclinic instability, along with the symmetric and centrifugal modes (Allen and Newberger 1998; Molemaker et al. 2015; Dewar et al. 2015; Gula et al. 2016), plays a leading-order role in the dynamics of both coastal (Brink 2016; Brink and Seo 2016; Hetland 2017), and deep ocean BBLs, with implications for closing the upwelling limb of the abyssal circulation (Ferrari et al. 2016; Callies 2018).

The article is organized as follows. In section 2 the equations governing linear perturbations in the BBL are developed, relevant quasigeostrophic results are briefly reviewed, and the parameter dependence of the baroclinic mode in the nongeostrophic limit is explored numerically. In section 3 a steady one-dimensional solution for a BBL in the presence of bottom-intensified turbulence, with parameters similar to those observed along the mid-Atlantic ridge, is shown to support growing baroclinic instability in a weakly stratified outer layer on the order of hundreds of meters thick. In section 4 numerical solutions of the time-dependent, one-dimensional, Ekman adjustment problem for flow along a slope are likewise shown to be unstable to submesoscale baroclinic instability in a well-mixed bottom boundary layer that is on the order of tens of meters thick. An example nonlinear simulation is presented in section 5, allowing us to comment briefly on the finite-amplitude behavior, and in section 6 the baroclinic mode is discussed in relation to symmetric and centrifugal instability, both of which can also be present in BBLs over sloping topography. Major findings and broader implications are summarized in section 7.

2. Baroclinic instability in the bottom boundary layer

a. Theory

We consider the stability characteristics of a linearized Boussinesq system, in a coordinate system rotated to align with a linear slope in the direction (Fig. 1), where throughout the paper we will indicate quantities in the conventional, nonrotated, coordinate system using a caret notation, background fields using uppercase variables, and perturbation quantities using lowercase variables. The slope-normal coordinate is defined positive upward. An arbitrary background flow can be expressed in the rotated frame as , where θ is the angle of the topographic slope (Fig. 1). The equations governing the linear evolution of the perturbations are then
e1
e2
e3
e4
e5
where f is the Coriolis parameter, ν is the viscosity, κ is the diffusivity, ∇χ = (χx, χy, χz) is the gradient operator in the rotated frame, and subscripts denote differentiation. Note that, consistent with Stone (1971), we retain nonhydrostatic effects but make the traditional approximation for the Coriolis force, ignoring terms involving the horizontal components of the Coriolis parameter. This is an accurate approximation for the problems considered here, and further allows the problem to remain invariant under rotation of the slope direction. Boundary conditions depend on the particular problem configuration and will be discussed in the relevant following sections.

For simplicity, we will only consider background flows that are invariant in the rotated along- and across-slope direction, which excludes slope-normal background flows (i.e., W = 0). These restrictions are consistent with the conceptual model of a BBL of uniform thickness along a linear slope (sections 3 and 4) and still allow for a vertically and horizontally sheared background flow in the nonrotated coordinate system (where the horizontal shear is a consequence of the sloping topography, as shown schematically in Fig. 1b). We will also retain only the slope-normal component of the turbulent diffusion terms, which for an isotropic turbulent diffusivity can be understood as an assumption of small aspect ratio, H2/L2 ≪ 1, where H and L scale the slope-normal and rotated horizontal dimensions, respectively. This choice does not affect the conclusions of this article; however, we emphasize that the properties of diapycnal and isopycnal turbulent mixing in the BBL are currently poorly constrained from available observations, and it is possible that under some circumstances—for example, along very steep topography—the full three-dimensional diffusion operator may be important.

Assuming perturbations of the form gives the eigenvalue problem,
e6
e7
e8
e9
e10
Given a specified background state, the eigenvalues of (6)(10) can be found numerically. We do this by projecting the governing equations onto a Chebyshev basis in z (n = 256 in all cases) and then finding the eigenvalues of the resulting matrix using the Dedalus equation-solving framework (Burns et al. 2016).
Later in the article, general profiles of background velocity and buoyancy will be considered using (6)(10) directly; however, to better illustrate the problem parameter dependence, it is useful to assume a background flow in the along-slope direction, with uniform vertical shear, , and stratification, , such that V = Λz secθ, and B = N2(z cosθ + x sinθ). Equations (6)(10) can then be nondimensionalized following Stone (1971) using , , , , , , , and , with primes denoting nondimensional quantities. Dropping grave and prime accents for clarity then gives
e11
e12
e13
e14
e15
The relevant nondimensional parameters are the Richardson number of the background flow, Ri = N22; a nonhydrostatic parameter, δ = f/Λ; the Ekman number, E = Aυ/(fH2); and a slope parameter,
e16
which gives the ratio of the topographic slope to minus the isopycnal slope.2 Alternately, the slope parameter can be expressed as α = (SRi)1/2, where S = (N tanθ/f)2 is the slope Burger number. For the analyses of this paper, we retain nonhydrostatic terms (except for the nonlinear simulation of section 5); however, the nonhydrostatic parameter is generally small, even for weakly stratified layers, and hence the inviscid stability characteristics of this system of equations are primarily determined by the Richardson number and α.

To motivate the following sections, we first consider the baroclinic growth rates for an idealized profile of stratification (N2) and shear (Λ), with varying slope angle. Figure 2 shows profiles of background velocity and buoyancy, where we assume there is an inviscid along-slope flow, V(z) = Λz secθ, and a stratified interior overlying a weakly stratified BBL. Note that these profiles are held fixed in the nonrotated frame, hence the isopycnal slope in the rotated frame changes as a function of the slope angle. The inviscid problem only requires boundary conditions on the slope-normal velocity, which are given by w = 0 at z = 0 and w = −u tanθ at z = 1000 m, representing a rigid horizontal upper boundary 1000 m from the bottom. Growth rates of the baroclinic mode (k = 0) are calculated numerically and are shown as a function of increasing slope angle in Fig. 3. The full-depth baroclinic mode can be seen at wavelengths of ~30 km. This mode is strongly modulated by slope angle, with maximum growth rates reduced, and shifted to higher wavenumber, as the slope angle increases, consistent with the expectation from quasigeostrophic (QG) theory (discussed further below). The fastest-growing instability, however, is found at wavelengths near 1500 m, with perturbation structure and energetics shown in Fig. 4. This is a bottom-enhanced submesoscale baroclinic instability, resulting from the interaction of counterpropagating Rossby waves along the lower-boundary and on the interface of changing stratification between the boundary layer and the interior. These ocean bottom modes are thus similar to baroclinic instabilities in the atmospheric boundary layer (Nakamura 1988) and mirror the structure of surface mixed layer instabilities (Boccaletti et al. 2007; Callies et al. 2016). As shown in Fig. 3, the growth rates of the submesoscale mode are relatively insensitive to increasing slope angle, suggesting that submesoscale baroclinic instability may be a robust feature of the bottom boundary layer. To understand this behavior, we will first briefly review relevant results from the QG limit and then explore parameter space in the non-QG limit relevant for submesoscale instabilities.

Fig. 2.
Fig. 2.

Example idealized vertical structure of a stratified interior overlying a weakly stratified BBL.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Fig. 3.
Fig. 3.

Growth rates, normalized by f −1, for the idealized profile shown in Fig. 2 for various slope angles (legend). The low-wavenumber, full-depth baroclinic modes are indicated by the dashed box. Growth rates of the low-wavenumber instability are highly modulated by the slope angle, whereas the high wavenumber features are robust to changing slope angle.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Fig. 4.
Fig. 4.

Structure of the fastest-growing perturbation (l ≈ 6 × 10−4 m−1) for θ = 5 × 10−3. Velocities are normalized relative to the maximum across-slope velocity and buoyancy relative to the maximum buoyancy perturbation. Kinetic energy tendency terms are defined in appendix A, with buoyancy production given by VBP and shear production by the sum of the lateral and vertical shear production terms (LSP + VSP).

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

b. Quasigeostrophic baroclinic instability over sloping topography

An analytical solution for the growth rate of inviscid baroclinic instabilities over sloping topography in the QG limit (assuming small Rossby number, large Richardson number, and small slope angle) was first provided by Blumsack and Gierasch (1972), considering instabilities in the Martian atmosphere. This work was later extended by Mechoso (1980) to include a sloping upper boundary, introducing an additional parameter,
e17
the ratio of the upper boundary slope θub, to the negative of the isopycnal slope. The case where αub = α can therefore be considered as a simple prototype model of baroclinic instability in a bottom boundary layer that is of uniform thickness in the across-slope direction, ignoring deformation of the upper boundary of the BBL (Fig. 1b). The use of a rigid lid approximation is accurate as long as the stratification contrast between the lower-layer and interior is large (Eady 1949; Boccaletti et al. 2007; Callies et al. 2016), in which case including a deformable upper boundary adds an additional free parameter and leads only to small quantitative modifications. Later in the paper, the assumption of a rigid boundary on the lower-layer will be removed. For the case of uniform stratification and vertical shear, the instability growth rate is given by (Mechoso 1980)
e18
where ωi is the growth rate normalized by f−1, and we have introduced an alternate normalization of the wavenumber using the deformation radius NH/f instead of ΛH/f (such that ), as is appropriate for quasigeostrophic dynamics. This alternate normalization will prove useful for comparing the stability properties with varying Ri (below).

These growth rates are shown as a function of α in Fig. 5 for the case of a horizontal upper boundary (left) and for a sloping upper boundary, αub = α (right). In both cases there is a region of instability for any α > −1, with increasing positive values of α associated with reduced growth rates, and an increasingly narrow region of instability at high wavenumbers, as seen for the low-wavenumber modes in Fig. 3. While linear theory predicts finite growth rates at arbitrarily high α, in reality nonlinear (Pedlosky 2016), viscous, and non-QG (section 2c) effects may provide an upper bound on physically relevant wavenumbers. With an upper boundary parallel to the topography, there is no longer a barotropic potential vorticity (PV) gradient, and hence there is no longer a long-wave cutoff for α < 0 (Fig. 5, right). The sloping boundaries also lead to a destabilization of low wavenumbers for α < 0, with maximal growth rates for α = −0.5. This is a consequence of parcel trajectories being forced by the boundary slope to cross isopycnals in a manner that achieves maximal extraction of available potential energy (Mechoso 1980). The increase of growth rates as can be understood as a resonance between baroclinic and barotropic Rossby waves (Ripa 1999, 2001), associated with a uniform downslope flow extracting available potential energy. These results will now be compared to the non-QG limit.

Fig. 5.
Fig. 5.

Baroclinic instability growth rates from (18) as a function of slope parameter α and wavenumber normalized by the deformation radius NH/f . Growth rates are normalized by f−1, and it is assumed that Ri = 100. Shown are (left) a horizontal upper boundary and (right) an upper boundary that matches the topographic slope, as shown schematically for α > 0, and α < 0 in each plot.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

c. Nongeostrophic baroclinic instabilities over sloping topography

To explore parameter space in the non-QG limit, using (11)(15) in the rotated frame, we continue to assume the flow is inviscid, with uniform background vertical stratification, N2, and along-slope flow with uniform slope-normal shear. This allows the background velocity field to be written in the rotated frame as V(z) = Λz secθ, and in the nonrotated frame as , which has both vertical and horizontal vorticity (section 6). We also assume a rigid upper boundary, aligned parallel to the topography (i.e., αub = α, such that the boundary conditions are w = 0 at z = 0, H), a configuration that approximates a boundary layer of uniform thickness across the slope (Fig. 1, lower panel). The focus here is on the baroclinic mode, and hence we assume perturbation variables are uniform in the across-slope direction (k = 0); other instability modes (k ≠ 0) are discussed in section 6. This simplified model can thus be thought of as Stone’s (1966) nongeostrophic baroclinic instability problem adapted for a slope. In the inviscid and hydrostatic limit, the stability characteristics of (11)(15) are determined by Ri and α, which are varied below.

In the nondimensionalization used here, l acts as a Rossby number, quantifying advection, vis-à-vis Doppler shifting of the perturbations by the mean flow, relative to the Coriolis term. Non-QG effects will therefore become important when , which indicates that the flow is evolving quickly relative to the inertial time scale, violating the QG assumption of slow dynamics. This can be seen in Fig. 6 for Ri = 10 and Ri = 1. In each case the low wavenumber growth rates are similar to the expectation from QG theory, whereas wavenumbers with (or in the plots ) have growth rates reduced by non-QG effects. These effects are most relevant for positive slope angles, α > 0, where growth rates in the high-wavenumber tail are strongly reduced. Considering the instability energetics, topographic reduction in growth rates for α > 0 can result from a reduction of the vertical buoyancy production by the across-slope advection of buoyancy perturbations (appendix A). This reduction of the vertical buoyancy production scales as α/Ri and is therefore negligible in the QG limit, but can become significant for small Ri.

Fig. 6.
Fig. 6.

Growth rates for Ri = (left) 10 and (right) 1, as a function of slope parameter α and wavenumber normalized by the deformation radius NH/f . Growth rates are normalized by f−1, and Λ = 5f.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

An alternate, physically intuitive, way to consider the parameter dependence of the instabilities is to hold the slope angle and horizontal buoyancy gradient fixed (Fig. 7). Varying the stratification thus varies the Richardson number and α, approximating interior stratification encountering deep topography, and adjusting within a turbulent boundary layer (sections 3 and 4). At low Richardson numbers the growth rates increase, and higher wavenumbers become unstable, consistent with the decreasing stratification allowing greater vertical penetration and coupling of the boundary waves. Assuming values typical of interior flows (V ~ 0.1 m s−1, f = 10−4 s−1) the most unstable modes for Ri < 10 have a horizontal scale of O(5) km and growth rates on the order of 0.5–3 inertial periods.

Fig. 7.
Fig. 7.

Growth rates for a fixed topographic slope angle (θ = 5 × 10−3) and vertical shear Λ = 10−3 s−1, with varying N2. Wavenumbers are normalized by ΛH/f, growth rates by f−1.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

It is important to note that for isopycnals that intersect the topography at right angles, a reasonable first approximation to the BBL structure, the slope parameter is given by α = tan2θ, which for small slope angle can be approximated as αθ2. It is estimated that approximately 95% of the world’s ocean has (Costello et al. 2010), implying small α values will be most physically relevant in the BBL. Further, for a BBL with Ri = 1, as is anticipated for a bottom boundary layer that has adjusted after turbulent mixing (Haine and Marshall 1998; Fox-Kemper et al. 2008; Taylor and Ferrari 2009), α = S1/2. Observations of weakly stratified BBLs in both coastal and deep oceans suggest that S ≪ 1 (Armi and Millard 1976; Stahr and Sanford 1999; Moum et al. 2004), and hence α ≪ 1 may be common in the BBL. Finally, the change in the slope parameter with increasing slope angle is given by ∂α/∂θ = NRi1/2 sec2θ/f, which is small for BBLs with low stratification and low Richardson numbers. This suggests that the properties of the BBL instability are less sensitive to changes in slope than the full-depth instability, for which stratification and Richardson number are larger. This is consistent with the robustness of the submesoscale mode shown in Fig. 3. Thus, despite the reduction of growth rates at large α evident in Fig. 6, the growth rates of baroclinic instability can remain significant in weakly stratified layers over topography, where both Ri and α are small (sections 3 and 4), and the submesoscale mode is anticipated to be a robust feature of the ocean bottom boundary layer.

Up until this point we have considered inviscid flows with greatly simplified vertical structures for the background buoyancy and shear fields, as an aid to interpreting the problem parameter dependencies. In the following sections, we will consider the stability characteristics of the bottom boundary layer structure predicted by two somewhat more realistic dynamical theories of the ocean bottom boundary layer. The first results from bottom intensified turbulent mixing, the second from Ekman adjustment of the bottom boundary layer to an imposed interior flow.

3. Bottom-intensified turbulent mixing

Available potential energy near the bottom boundary can be generated by turbulent mixing associated with processes such as the breaking of internal waves over sloping topography. Observations show that turbulence over rough topography is strongly bottom enhanced, suggesting there should be a dipole of vertical velocities along topography, with downwelling in an outer layer a few hundred meters thick, where diapycnal buoyancy fluxes are divergent, and upwelling in an inner layer adjacent to the bottom, where diapycnal buoyancy fluxes are convergent (e.g., Polzin 1997; Ledwell et al. 2000; St. Laurent et al. 2001; Ferrari et al. 2016; McDougall and Ferrari 2017). However, 1D boundary layer solutions, with turbulent viscosities based on observations, imply outer-layer stratification much weaker than observed, and hence weak vertical velocity dipoles (Callies 2018). Baroclinic instabilities in the mixing layer would alter the 1D buoyancy budget, with important implications for closing the upwelling branch of the abyssal circulation, and hence here we consider the stability characteristics of 1D boundary layer solutions with turbulent diffusivities based on observations.

To do this we solve the governing equations for the steady background state over sloping topography, with uniform interior stratification, and no across-slope variations of the boundary layer quantities,
e19
e20
e21
where denotes a specified interior stratification far from the boundary, and the total background buoyancy field is given by , such that represents the portion of the buoyancy not associated with the imposed interior stratification. Boundary conditions are
e22
e23
Practically, the upper-boundary conditions are applied at a finite height, sufficiently far from the lower boundary to not influence the dynamics, which here we set to be z = 2500 m. These equations are solved numerically, using Dedalus, with an idealized, bottom-enhanced profile of turbulent mixing,
e24
Parameters are chosen to approximate observations taken during the Brazil Basin Tracer Release Experiment (Ledwell et al. 2000):
e25
e26
These solutions are used to define the background fields in the eigenvalue problem, (6)(10), with perturbation boundary conditions given by
e27
e28
The free-slip upper-boundary conditions on momentum are necessary to satisfy a no-stress boundary condition at the surface (Thomas and Rhines 2002; Wenegrat and McPhaden 2016).

The solution for the background fields is shown in Fig. 8. There is an inner layer of thickness , where isopycnals are approximately normal to the slope (Callies 2018). This layer is approximately 6 m thick for these parameters. Above this is an outer layer, with thickness , approximately 1000 m for these parameters. Stratification in the outer layer is reduced by a factor of from the interior values, where is the interior slope Burger number (Garrett 1991). The across-slope momentum balance in the outer layer is approximately , which, when combined with the force balance, , gives , that is, the along-slope flow is in geostrophic balance due to the projection of the slope-normal gradient on the true horizontal. The lower portion of the outer layer has Ri ~ O(10–100) and hence is stable to symmetric instability, but, as shown in Fig. 9, it can support growing baroclinic instability. The fastest-growing instability for these parameters has a wavelength of ~5 km and an e-folding time scale of ~5 days. For comparison, assuming an inviscid layer with a rigid lid at H = 500 m, ΛH = −0.015 m s−1, and Ri = 25, the simplified analysis of section 2c would predict a fastest-growing wavelength of ~6 km and a growth rate of ~3 days.

Fig. 8.
Fig. 8.

Solution to (19)(21) with parameters as given in (25) and (26). In the upper-left panel the along-slope velocity V is shown in orange and the across-slope velocity U is shown in blue.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Fig. 9.
Fig. 9.

Growth rate of the baroclinic mode for the base state shown in Fig. 8.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

The perturbation structure and kinetic energy tendency of the fastest-growing mode are shown in Fig. 10. Phase lines are inclined into the mean shear (note f < 0), and the instabilities grow by releasing available potential energy in the outer layer, with maximum vertical buoyancy fluxes near z = 250 m. Perturbations quantities are enhanced near the bottom boundary, except for in the thin inner layer where dissipation dominates the energetics. At finite amplitude the eddy vertical buoyancy fluxes will provide a restratifying tendency, altering the simple 1D balance given by (19)(21), and potentially modifying the structure of the vertical velocity dipole in the boundary layer. The net secondary circulation along sloping topography may therefore involve both eddy and Eulerian components, with implications for the dynamics of the abyssal circulation (Callies 2018).

Fig. 10.
Fig. 10.

Structure of the fastest-growing baroclinic mode for the base state shown in Fig. 8, discussed in section 3. Velocities are normalized relative to the maximum across-slope velocity and buoyancy relative to the maximum buoyancy perturbation. Kinetic energy tendency terms are defined in appendix A, with buoyancy production given by VBP, shear production by the sum of the lateral and vertical shear production terms (LSP + VSP), and dissipation by DKE.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

4. Stability of the bottom Ekman layer

The classic analysis of the laminar bottom boundary layer holds that, in the presence of an interior flow along the slope, a bottom Ekman layer must also develop to satisfy the bottom boundary condition on the momentum (Pedlosky 1979). Interior flow thus gives rise to across-slope flow in the Ekman layer that advects the across-slope buoyancy gradient, modifying the depth and stratification of the boundary layer (MacCready and Rhines 1991, 1993). For interior flow in the direction of topographic wave propagation (topography sloping upward to the right of the direction of flow in the Northern Hemisphere) the resulting Ekman flow is downslope, advecting light water under dense, causing convective mixing and deepening of the BBL (Garrett et al. 1993; Moum et al. 2004), an analog to the surface “Ekman buoyancy flux” associated with downfront winds at a surface mixed layer front (Thomas 2005). For interior flow in the opposite direction, the Ekman flow is upslope, increasing boundary layer stability by advecting dense water under light. In each case this adjustment process generates available potential energy in the BBL (Umlauf et al. 2015), and we demonstrate here that this available potential energy can allow for growing baroclinic instability in the BBL.

Following MacCready and Rhines (1993), we consider the time-dependent adjustment of a BBL over a linear slope. The initial condition is a barotropic, geostrophically balanced interior flow and uniform stable vertical stratification . The dimensional along-slope velocity can then be written as , with governing equations,
e29
e30
e31
where is defined in section 3 following (21). Boundary conditions on the background flow are given by (22) and (23), and on the perturbations by (27) and (28), with the upper boundary applied at z = 150 m, which is sufficiently far from the lower boundary to not modify the characteristics of the BBL solution. Solutions to (29)(31) are found using the General Ocean Turbulence Model (GOTM; Umlauf et al. 2005), modified to solve the rotated equations of motion. The barotropic interior flow is applied as an initial condition and is maintained during integration by imposing a steady barotropic pressure gradient force in the across-slope momentum equation. Turbulence closure is provided using a second-moment kε scheme, with a prognostic equation for the turbulence kinetic energy, and background viscosity and diffusivity of ν = κ = 10−5 m2 s−1. The bottom boundary conditions on momentum are implemented using a log-layer formulation. Detailed discussion of the use of this turbulence parameterization in a BBL over topography can be found in Umlauf and Burchard (2011).

a. Downwelling favorable interior flow

The numerical solution for the case of downwelling-favorable interior flow with magnitude typical of a deep western boundary current (Toole et al. 2011), = 0.1 m s−1, is shown in Fig. 11. A well-mixed turbulent bottom boundary layer is rapidly established by convective mixing due to the across-slope buoyancy advection. This leads to a growing BBL, which reaches ~30 m thickness after only approximately five inertial periods. As the BBL thickness increases, the along-slope flow approaches thermal wind balance, and U(z) → 0, a state known as the “arrested” Ekman layer. An estimate of the Ekman shutdown time scale can be formed as (MacCready and Rhines 1993), where Mx is the across-slope transport, which for this simulation gives τEk ~ 10 inertial periods.

Fig. 11.
Fig. 11.

BBL response to a downwelling-favorable interior flow, = 0.1 m s−1, with f = 10−4 s−1, = 3.5 × 10−3 s−2, and θ = 10−2, as in MacCready and Rhines (1993). Viscosity ν follows κ closely and hence is not shown here.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Instability growth rates (Fig. 12) are calculated using snapshots of the numerical solutions for U, V, B, κ, and ν as the background fields in the eigenvalue problem, (6)(10). The baroclinic mode emerges after approximately five inertial periods, with growth rates increasing and wavenumbers decreasing as the boundary layer deepens in time. Growth of the baroclinic mode becomes faster than the Ekman adjustment time scale after only t ≈ 0.5τEk, emphasizing that these instabilities grow rapidly relative to the shutdown process. High-frequency variability is associated with weak inertial oscillations modifying the boundary layer structure. The structure of the fastest-growing mode at t = 10 inertial periods is shown in Fig. 13. Vertical buoyancy production of kinetic energy dominates the instability growth, with dissipation reducing the kinetic energy tendency by a factor of approximately 1/2. Shear production is a small contribution to the vertically integrated kinetic energy tendency and is dominated by the lateral shear production, although there are locally significant contributions from both vertical and lateral shear production terms at the top of the boundary layer.

Fig. 12.
Fig. 12.

Hovmöller plot of the baroclinic growth rates as a function of wavenumber and time, normalized by the inertial period, for the case of downwelling favorable interior flow (Fig. 11). Growth rates are normalized by f−1.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Fig. 13.
Fig. 13.

Structure of the fastest-growing baroclinic mode for the downwelling-favorable Ekman layer (Fig. 11) at t = 10 inertial periods, discussed in section 4. Velocities are normalized relative to the maximum across-slope velocity, buoyancy relative to the maximum buoyancy perturbation. Kinetic energy tendency terms are defined in appendix A, with buoyancy production given by VBP, shear production by the sum of the lateral and vertical shear production terms (LSP + VSP), and dissipation by DKE.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

It is important to emphasize that, for downwelling-favorable interior flow, the BBL can become unstable to symmetric or centrifugal instability, with growth rates that exceed the baroclinic mode (Allen and Newberger 1998). Future work will therefore consider the full nonlinear Ekman adjustment process; however, the findings presented here suggest that fast-growing submesoscale baroclinic instability likely plays an important role in the 3D adjustment of the BBL to interior flows, emerging after symmetric/centrifugal instability restratifies the boundary layer to a state of marginal stability (discussed further in section 6).

b. Upwelling favorable interior flow

When the interior flow is in the opposite direction, = −0.1 m s−1, the Ekman flow is upslope, which tends to restratify the BBL (Fig. 14). Initially in the simulation, shear production of turbulent kinetic energy leads to the development of a thin turbulent BBL, with a maximum depth of about 15 m after three inertial periods. A weak, low-wavenumber baroclinic mode is evident in the stability analysis (Fig. 15), with a peak growth rate of approximately 0.03f at t = 7.5 inertial periods. The growth rate, however, decreases in time as the upslope Ekman flow restratifies the thin BBL. Thus, aside from a slight initial transient destabilization of the BBL, upwelling-favorable interior flow leads to a BBL that is stable to the baroclinic mode. The transient behavior in this simple numerical integration does, however, suggest that more complex, time-dependent interior flows may lead to baroclinically unstable BBLs, even in the case of time-mean flows that are upwelling favorable (Brink 2016; Brink and Seo 2016).

Fig. 14.
Fig. 14.

As in Fig. 11, but for upwelling favorable interior flow, = −0.1 m s−1.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Fig. 15.
Fig. 15.

As in Fig. 12, but for the upwelling favorable case shown in Fig. 14.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

5. Nonlinear simulation

In this section we present the results of an example of idealized nonlinear simulation. Solutions are found using Dedalus, in a computational domain that is doubly periodic in the rotated horizontal coordinate system. This is achieved by imposing a fixed across-slope buoyancy gradient, (where is an assigned interior stratification), and solving for horizontally periodic departures from the background fields, similar to the “frontal zone” configuration often used in studies of the surface boundary layer (e.g., Taylor and Ferrari 2010). Full details of the model configuration are given in appendix B.

The model is initialized with a uniformly stratified interior (NI = 3.4 × 10−3 s−1) and a barotropic along-slope interior flow of 0.1 m s−1. The BBL is initialized in thermal wind balance, with a thickness of 100 m, such that Ri = 1.5 in the BBL (Fig. 16). The slope angle is θ = 10−2, giving a slope parameter in the BBL of α ≈ 0.15. For simplicity we consider the case of weak slope-normal viscosity and diffusivity, ν = κ = 10−5 m2 s−1, and, for numerical stability, a biharmonic diffusivity in the rotated horizontal (νh = 2 × 105 m4 s−1) (appendix B). A more complete exploration of parameter space using nonlinear simulations, including the use of physically realistic turbulence closures, will be the subject of future work.

Fig. 16.
Fig. 16.

Initial conditions for along-slope velocity (color) and buoyancy (gray contours) for the nonlinear simulation discussed in section 5 and appendix B.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

Snapshots of the instability evolution are shown in Fig. 17. By day 10 the baroclinic mode has reached finite amplitude and is clearly evident in both the buoyancy and slope-normal velocity fields. The along-slope wavelength of the dominant instability, λ ≈ 6.4 km, is similar to that predicted by linear stability analysis of the domain-averaged buoyancy and velocity fields at day 1.5, λ ≈ 5.5 km. Importantly, the slope-normal velocities exceed 250 m day−1 (Fig. 17, bottom row), comparable to values associated with baroclinic instability in the surface boundary layer. These large-magnitude vertical velocities suggest that eddy fluxes in the BBL have the potential to affect a wide range of BBL processes, including, for example, nutrient fluxes between the BBL and interior, sediment transport, and coastal hypoxia.

Fig. 17.
Fig. 17.

Snapshots of (top) buoyancy and (bottom) slope-normal velocity (w) evaluated at 50 m above the bottom.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

The energetics of the nonlinear simulation are shown in Fig. 18. In the top panel the domain-integrated eddy kinetic energy (EKE) is shown, defined as , where primes indicate departure from the rotated horizontal average. An estimate of EKE can be calculated from linear theory as , where to = 1.5 days, the hat notation indicates the wavenumber spectrum (calculated from the model output at t = to), and ωi is determined by linear stability analysis of the domain-averaged profiles at t = to. This estimate thus accounts for the distribution of energy across multiple wavenumbers with varying growth rates, and agrees closely with the growth rate of EKE in the nonlinear simulation over the first 10 days. During this time the instability growth is dominated by the buoyancy production of EKE, as indicated in Fig. 18 (bottom), typical of baroclinic instability. Later, after day 10, the available potential energy of the initial condition has been depleted, EKE is saturated, and buoyancy and shear production of EKE are largely compensated by dissipation (Fig. 18, bottom).

Fig. 18.
Fig. 18.

(top) Domain-integrated EKE compared to the prediction of linear theory (calculated as discussed in section 5) and (bottom) domain-averaged EKE tendency (bottom) with components as indicated in the legends.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

In this simulation, domain-averaged vertical buoyancy fluxes in the BBL reach ~ O(10−8) m2 s−3. This value is comparable to baroclinic instability in the surface boundary layer (Boccaletti et al. 2007), and importantly, is several orders of magnitude larger than values typically associated with the global breaking of internal waves over rough topography, where it is often assumed ~ O(10−10) m2 s−3 (Nikurashin and Ferrari 2013; Ferrari et al. 2016). This simulation thus suggests that finite-amplitude baroclinic instability in the BBL has the potential to both modify the dynamics of the BBL, and to affect the interior dynamics of both coastal and deep oceans, through large, restratifying, eddy vertical buoyancy fluxes.

6. Instability regimes

In addition to the baroclinic mode, (6)(10) also admit symmetric and centrifugal instabilities (Haine and Marshall 1998), each of which may be active in the ocean BBL under certain conditions (Allen and Newberger 1998; Brink 2012; Molemaker et al. 2015; Gula et al. 2016). While the focus of this work is primarily on baroclinic instability, it is useful to also briefly summarize where in parameter space each of these modes is anticipated to be dominant.

An illustrative example can be formed by assuming a geostrophically balanced along-slope flow in the direction of topographic wave propagation, with uniform stable stratification, and uniform slope-normal shear, such that (Fig. 1b). Note that the assumed linear background shear profile excludes the Kelvin–Helmholtz mode, otherwise expected at Ri < 0.25 (Vanneste 1993). The fastest-growing instability type will therefore depend on the sign of f times the Ertel PV of the background state (Haine and Marshall 1998), , which for this simplified configuration can be written nondimensionally as
e32
For q′ < 0, the fastest-growing instabilities will be associated with overturning circulations in the slope-normal and across-slope plane, l = 0, which grow by extracting kinetic energy from the background flow (Hoskins 1974). These instabilities can be of the symmetric, centrifugal, or mixed types, with the ratio of the lateral shear production (LSP) to vertical shear production (VSP) providing a useful discriminator (appendix A; Thomas et al. 2013),
e33
In the case that LSP/VSP < 1 the instabilities are of the symmetric type, for LSP/VSP > 1 and α/Ri > 1 (negative absolute vertical vorticity) the instabilities are of the centrifugal type, and for LSP/VSP > 1 and α/Ri < 1 (positive absolute vertical vorticity) the instabilities are a mixed SI/CI instability (Thomas et al. 2013). For q′ > 0 the baroclinic mode k = 0 will be the fastest-growing instability, with perturbation kinetic energy increasing through the extraction of the background state available potential energy (buoyancy production).

An example regime diagram is shown in Fig. 19, as a function of the slope Burger number and the Richardson number. The q′ = 0 line divides the domain (dashed red line), with the entire area to the right of the line linearly unstable to baroclinic instability, although instabilities at small α are likely the most physically relevant (section 2c). For q′ < 0 the domain is further divided by the energetics, with LSP/VSP = 1 shown from theory [(33); dashed blue line in Fig. 19], assuming δ/Ri1/2 = f/N = 0.1. Below this line the symmetric mode grows by extracting mean kinetic energy through the vertical shear production. Above this line the energetics are dominated by the lateral shear production, either in a pure centrifugal instability mode, to the left of the green line indicating zero absolute vertical vorticity (1 − α/Ri = 0), or in a mixed symmetric/centrifugal type mode, for positive absolute vertical vorticity to the right of the green line. Numerical solutions of (6)(10) [with l = 0, k = 20f/(NH)] were also used to calculate (solid blue line), demonstrating that the theory provides an accurate approximation to the energetics.

Fig. 19.
Fig. 19.

Example regime diagram for baroclinic instability (BI), symmetric instability (SI), centrifugal instability (CI), and mixed instability (SI/CI) in the ocean BBL. The red line indicates q = 0, the green line indicates , and the dotted black lines indicate α = (1, 2)α. The dashed blue line is the theoretical expectation for LSP/VSP = 1, given by (33), and the solid blue line is calculated using numerical solutions of (6)(10) with δRi−1/2 = 0.1, l = 0, and k = 20f/(NH). See section 6 for further explanation.

Citation: Journal of Physical Oceanography 48, 11; 10.1175/JPO-D-17-0264.1

For small Ri, the dominant instability types in the BBL are therefore likely to be of the centrifugal, symmetric, or mixed types, with the centrifugal mode becoming dominant for large S, where the topography begins to act as an effectively vertical boundary (Dewar et al. 2015). For and small S, the symmetric mode is likely most relevant (Allen and Newberger 1998). In the case that , the baroclinic mode will dominate. Further, when considering the nonlinear dynamics of the BBL, it is expected that the time-dependent adjustment of the background fields by finite-amplitude instabilities will lead to a traversal of parameter space as the shear and stratification are adjusted by eddy fluxes. Thus, even for initial background conditions that have fastest-growing modes of the overturning type, the baroclinic mode may emerge later as the instabilities restratify the boundary layer. For example, during the Ekman adjustment to a downwelling-favorable interior flow, potential vorticity is removed at the boundary (Benthuysen and Thomas 2012), and in the BBL the potential vorticity can become negative. The fastest-growing modes of the 1D solutions of section 4a (Fig. 11) will therefore likely be of the symmetric or centrifugal type. However, when these instabilities reach finite amplitude they will adjust the boundary layer stratification such that Ri ≈ 1, at which point the baroclinic mode will dominate (Haine and Marshall 1998; Brink and Cherian 2013). Likewise, even in the case where the boundary layer is forced to remain near q = 0, as, for instance, might occur in the case of continuing downslope Ekman transport, growing baroclinic instability can still be present (Callies and Ferrari 2018).

7. Summary

In this article, we considered baroclinic instability over sloping topography, extending earlier QG results to the nongeostrophic regime appropriate for the ocean BBL. Importantly, weakly stratified BBLs can support a submesoscale baroclinic instability, a BBL counterpart to the surface baroclinic mixed layer instability (Boccaletti et al. 2007). In the BBL these submesoscale instabilities are relatively insensitive to topographic slope, suggesting they are likely a robust part of the dynamics of the ocean BBL, along with the symmetric (Allen and Newberger 1998) and centrifugal (Molemaker et al. 2015) modes.

Two 1D theories of the BBL structure over topography were also shown to result in solutions that are susceptible to growing baroclinic instability, with important implications for our understanding of the dynamics of both the BBL and the interior. In the case of a BBL generated by turbulent mixing, a thick outer layer supports a rapidly growing baroclinic mode that releases available potential energy and provides a restratifying tendency that is not accounted for in the 1D formulation. This suggests that deep submesoscale baroclinic instabilities could alter the stratification and across-slope flow along topographic features such as the mid-Atlantic ridge, with important implications for closure of the abyssal circulation (Callies 2018). Likewise, interior flow along a slope generates available potential energy in the BBL through an across-slope Ekman flow, leading to baroclinic growth rates that quickly exceed the rate of Ekman adjustment. These instabilities thus have the potential to affect both the BBL dynamics and the exchange between the boundary layer and interior, in both the coastal and deep ocean.

Extensive work on the surface boundary layer has demonstrated the leading-order role that submesoscale baroclinic instability plays in boundary layer restratification (Fox-Kemper et al. 2008; Johnson et al. 2016; Su et al. 2018), vertical exchange between the boundary layer and interior (Mahadevan and Tandon 2006), and the surface flux of potential vorticity (Wenegrat et al. 2018). The similarities between the BBL and surface mode (section 2), and the large-amplitude eddy fluxes found in the nonlinear simulation of section 5, suggest that many of the same results may hold in the BBL. However, there are also important differences between the surface and bottom boundary layer that require further consideration. For example, as noted in section 2, topography shapes parcel trajectories, modifying the instability energetics and potentially the associated eddy fluxes of buoyancy and tracers (Spall 2004; Isachsen 2011; Brink 2013). Further, variations in topography, not considered here, can impose an external length scale on instabilities (de Szoeke 1983), generate slope-normal flow and secondary circulations (Benthuysen et al. 2015), and modify instability growth rates (de Szoeke 1983; Solodoch et al. 2016), all of which may alter the role of baroclinic instability in BBL dynamics. Relaxing the assumption of across-slope uniformity in boundary layer height or structure may also lead to an arrest of the inverse cascade due to the topographic beta effect (Rhines 1975; Brink 2012). Full nonlinear simulations to explore the finite-amplitude behavior, and parameter dependence, of the submesoscale BBL baroclinic mode will therefore be the subject of a future paper.

Acknowledgments

The authors thank the contributors to the Dedalus Project (www.dedalus-project.org) and GOTM (www.gotm.net), both of which were extremely useful for the analyses in this paper. Authors JOW and LNT were supported during this work by NSF Grant OCE-1459677. Comments from three anonymous reviewers are also gratefully acknowledged.

APPENDIX A

Energetics

In this appendix we derive several results related to the instability energetics. For simplicity we assume hydrostatic dynamics, with background fields that are uniform in the along-slope direction, and diffusion only in the slope-normal direction. The eddy potential energy equation can be formed by taking b/N2 × (5), giving
ea1
where
ea2
is the conversion of mean potential energy to eddy potential energy through the horizontal buoyancy flux,
ea3
is the conversion of eddy potential energy to EKE through the vertical buoyancy production, and
ea4
is the dissipation of eddy potential energy. In this appendix, angle brackets indicate spatial averaging, both vertically and horizontally over a wavelength.
The EKE tendency can be formed by taking u × (1) + υ × (2), and using (3) and (4), giving
ea5
where
ea6
and
ea7
give the conversion of mean kinetic energy to EKE through vertical and lateral shear production, respectively. Dissipation of EKE is given by
ea8
For the baroclinic mode the EKE source is the vertical buoyancy production (VBP), which consists of two terms related to the slope-normal and across-slope advection of buoyancy. The ratio of these terms scale as
ea9
Noting that PEt > 0 requires 〈ub〉 < 0, the across-slope advection of buoyancy can be seen to reduce (increase) VBP for α > 0 (α < 0).
Finally, for the purposes of discriminating between centrifugal- and symmetric-type instabilities (section 6), it is useful to find the ratio of the lateral and vertical shear production terms. Following Thomas et al. (2013), and assuming a geostrophically balanced along-slope background flow, this ratio can be written as
ea10
where
ea11
Using , this can be written as
ea12
and
ea13
For the portion of parameter space considered in section 6 (Fig. 19), the fastest-growing mode has n = 0, and using small angle approximations,
ea14

APPENDIX B

Configuration of the Nonlinear Model

The nonlinear simulation of section 5 is performed using Dedalus (Burns et al. 2016). The domain is 32 km in each rotated horizontal direction and 1000 m in the slope-normal direction and is decomposed using 64 Fourier components in each horizontal direction and 256 Chebyshev modes in the vertical, giving horizontal and vertical resolutions of Δx = Δy = 500 m, and Δz = 0.075–6 m. The rotated equations are solved in a horizontally periodic domain by assuming a fixed mean across-slope gradient in buoyancy, such that , where is an imposed value of the interior stratification. This configuration is thus similar to the frontal-zone setup commonly used for studying surface mixed layer fronts in periodic domains (e.g., Taylor and Ferrari 2010). The initial condition is assumed to be an along-slope flow in hydrostatic and geostrophic balance that varies only in the slope-normal direction V(z). The rotated equations of motion can then be written in the hydrostatic limit, that is, and , as
eb1
eb2
eb3
eb4
eb5
In the above, uppercase quantities indicate the imposed initial conditions, and time-dependent, horizontally periodic quantities are denoted with the tilde notation. Specifically, the total velocity is decomposed as , the total pressure divided by the reference density as , and the buoyancy as . Horizontal mixing is parameterized using a biharmonic diffusion operator for numerical stability, . Boundary conditions are given by
eb6
eb7
All perturbation quantities are initialized from zero, except for , which is initialized with weak random noise. Note that no restriction is placed on the amplitude of the perturbations, and (B1)(B5) are fully nonlinear.

REFERENCES

  • Allen, J. S., and P. A. Newberger, 1998: On symmetric instabilities in oceanic bottom boundary layers. J. Phys. Oceanogr., 28, 11311151, https://doi.org/10.1175/1520-0485(1998)028<1131:OSIIOB>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Armi, L., and R. C. Millard, 1976: The bottom boundary layer of the deep ocean. J. Geophys. Res., 81, 49834990, https://doi.org/10.1029/JC081i027p04983.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Benthuysen, J., and L. N. Thomas, 2012: Friction and diapycnal mixing at a slope: Boundary control of potential vorticity. J. Phys. Oceanogr., 42, 15091523, https://doi.org/10.1175/JPO-D-11-0130.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Benthuysen, J., L. N. Thomas, and S. J. Lentz, 2015: Rapid generation of upwelling at a shelf break caused by buoyancy shutdown. J. Phys. Oceanogr., 45, 294312, https://doi.org/10.1175/JPO-D-14-0104.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blumsack, S. L., and P. J. Gierasch, 1972: Mars: The effects of topography on baroclinic instability. J. Atmos. Sci., 29, 10811089, https://doi.org/10.1175/1520-0469(1972)029<1081:MTEOTO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Boccaletti, G., R. Ferrari, and B. Fox-Kemper, 2007: Mixed layer instabilities and restratification. J. Phys. Oceanogr., 37, 22282250, https://doi.org/10.1175/JPO3101.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., 2012: Baroclinic instability of an idealized tidal mixing front. J. Mar. Res., 70, 661688, https://doi.org/10.1357/002224012805262716.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., 2013: Instability of a tidal mixing front in the presence of realistic tides and mixing. J. Mar. Res., 71, 227251, https://doi.org/10.1357/002224013807719473.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., 2016: Continental shelf baroclinic instability. Part I: Relaxation from upwelling or downwelling. J. Phys. Oceanogr., 46, 551568, https://doi.org/10.1175/JPO-D-15-0047.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., and D. A. Cherian, 2013: Instability of an idealized tidal mixing front: Symmetric instabilities and frictional effects. J. Mar. Res., 71, 425450, https://doi.org/10.1357/002224013812587582.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brink, K. H., and H. Seo, 2016: Continental shelf baroclinic instability. Part II: Oscillating wind forcing. J. Phys. Oceanogr., 46, 569582, https://doi.org/10.1175/JPO-D-15-0048.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Burns, K. J., G. M. Vasil, J. S. Oishi, D. Lecoanet, and B. Brown, 2016: Dedalus: Flexible framework for spectrally solving differential equations. Astrophysics Source Code Library, http://ascl.net/phpBB3/viewtopic.php?t=33963.

  • Callies, J., 2018: Restratification of abyssal mixing layers by submesoscale baroclinic eddies. J. Phys. Oceanogr., 48, 19952010, https://doi.org/10.1175/JPO-D-18-0082.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., and R. Ferrari, 2018: Baroclinic instability in the presence of convection. J. Phys. Oceanogr., 48, 4560, https://doi.org/10.1175/JPO-D-17-0028.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., R. Ferrari, J. M. Klymak, and J. Gula, 2015: Seasonality in submesoscale turbulence. Nat. Commun., 6, 6862, https://doi.org/10.1038/ncomms7862.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., G. Flierl, R. Ferrari, and B. Fox-Kemper, 2016: The role of mixed-layer instabilities in submesoscale turbulence. J. Fluid Mech., 788, 541, https://doi.org/10.1017/jfm.2015.700.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Costello, M. J., A. Cheung, and N. De Hauwere, 2010: Surface area and the seabed area, volume, depth, slope, and topographic variation for the world’s seas, oceans, and countries. Environ. Sci. Technol., 44, 88218828, https://doi.org/10.1021/es1012752.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • de Szoeke, R. A., 1983: Baroclinic instability over wavy topography. J. Fluid Mech., 130, 279298, https://doi.org/10.1017/S0022112083001093.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Dewar, W. K., J. C. McWilliams, and M. J. Molemaker, 2015: Centrifugal instability and mixing in the California Undercurrent. J. Phys. Oceanogr., 45, 12241241, https://doi.org/10.1175/JPO-D-13-0269.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Eady, E. T., 1949: Long waves and cyclone waves. Tellus, 1, 3352, https://doi.org/10.3402/tellusa.v1i3.8507.

  • Ferrari, R., A. Mashayek, T. J. McDougall, M. Nikurashin, and J.-M. Campin, 2016: Turning ocean mixing upside down. J. Phys. Oceanogr., 46, 22392261, https://doi.org/10.1175/JPO-D-15-0244.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fox-Kemper, B., R. Ferrari, and R. Hallberg, 2008: Parameterization of mixed layer eddies. Part I: Theory and diagnosis. J. Phys. Oceanogr., 38, 11451165, https://doi.org/10.1175/2007JPO3792.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Garrett, C., 1991: Marginal mixing theories. Atmos.–Ocean, 29, 313339, https://doi.org/10.1080/07055900.1991.9649407.

  • Garrett, C., P. MacCready, and P. Rhines, 1993: Boundary mixing and arrested Ekman layers: Rotating stratified flow near a sloping boundary. Annu. Rev. Fluid Mech., 25, 291323, https://doi.org/10.1146/annurev.fl.25.010193.001451.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gula, J., M. J. Molemaker, and J. C. McWilliams, 2015: Topographic vorticity generation, submesoscale instability and vortex street formation in the Gulf Stream: Topographic vorticity generation. Geophys. Res. Lett., 42, 40544062, https://doi.org/10.1002/2015GL063731.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gula, J., M. J. Molemaker, and J. C. McWilliams, 2016: Topographic generation of submesoscale centrifugal instability and energy dissipation. Nat. Commun., 7, 12811, https://doi.org/10.1038/ncomms12811.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Haine, T. W. N., and J. Marshall, 1998: Gravitational, symmetric, and baroclinic instability of the ocean mixed layer. J. Phys. Oceanogr., 28, 634658, https://doi.org/10.1175/1520-0485(1998)028<0634:GSABIO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hetland, R. D., 2017: Suppression of baroclinic instabilities in buoyancy-driven flow over sloping bathymetry. J. Phys. Oceanogr., 47, 4968, https://doi.org/10.1175/JPO-D-15-0240.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hoskins, B. J., 1974: The role of potential vorticity in symmetric stability and instability. Quart. J. Roy. Meteor. Soc., 100, 480482, https://doi.org/10.1002/qj.49710042520.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Isachsen, P. E., 2011: Baroclinic instability and eddy tracer transport across sloping bottom topography: How well does a modified Eady model do in primitive equation simulations? Ocean Modell., 39, 183199, https://doi.org/10.1016/j.ocemod.2010.09.007.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Johnson, L., C. M. Lee, and E. A. D’Asaro, 2016: Global estimates of lateral springtime restratification. J. Phys. Oceanogr., 46, 15551573, https://doi.org/10.1175/JPO-D-15-0163.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ledwell, J. R., E. T. Montgomery, K. L. Polzin, L. C. St. Laurent, R. W. Schmitt, and J. M. Toole, 2000: Evidence for enhanced mixing over rough topography in the abyssal ocean. Nature, 403, 179182, https://doi.org/10.1038/35003164.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lévy, M., P. Klein, A.-M. Tréguier, D. Iovino, G. Madec, S. Masson, and K. Takahashi, 2010: Modifications of gyre circulation by sub-mesoscale physics. Ocean Modell., 34, 115, https://doi.org/10.1016/j.ocemod.2010.04.001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • MacCready, P., and P. B. Rhines, 1991: Buoyant inhibition of Ekman transport on a slope and its effect on stratified spin-up. J. Fluid Mech., 223, 631661, https://doi.org/10.1017/S0022112091001581.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • MacCready, P., and P. B. Rhines, 1993: Slippery bottom boundary layers on a slope. J. Phys. Oceanogr., 23, 522, https://doi.org/10.1175/1520-0485(1993)023<0005:SBBLOA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mahadevan, A., and A. Tandon, 2006: An analysis of mechanisms for submesoscale vertical motion at ocean fronts. Ocean Modell., 14, 241256, https://doi.org/10.1016/j.ocemod.2006.05.006.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mahadevan, A., E. D’Asaro, C. Lee, and M. J. Perry, 2012: Eddy-driven stratification initiates North Atlantic spring phytoplankton blooms. Science, 337, 5458, https://doi.org/10.1126/science.1218740.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McDougall, T. J., and R. Ferrari, 2017: Abyssal upwelling and downwelling driven by near-boundary mixing. J. Phys. Oceanogr., 47, 261283, https://doi.org/10.1175/JPO-D-16-0082.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., 2016: Submesoscale currents in the ocean. Proc. Roy. Soc. London, A472, 20160117, https://doi.org/10.1098/rspa.2016.0117.

    • Search Google Scholar
    • Export Citation
  • Mechoso, C. R., 1980: Baroclinic instability of flows along sloping boundaries. J. Atmos. Sci., 37, 13931399, https://doi.org/10.1175/1520-0469(1980)037<1393:BIOFAS>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Molemaker, M. J., J. C. McWilliams, and W. K. Dewar, 2015: Submesoscale instability and generation of mesoscale anticyclones near a separation of the California Undercurrent. J. Phys. Oceanogr., 45, 613629, https://doi.org/10.1175/JPO-D-13-0225.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Moum, J. N., A. Perlin, J. M. Klymak, M. D. Levine, T. Boyd, and P. M. Kosro, 2004: Convectively driven mixing in the bottom boundary layer. J. Phys. Oceanogr., 34, 21892202, https://doi.org/10.1175/1520-0485(2004)034<2189:CDMITB>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nakamura, N., 1988: Scale selection of baroclinic instability—Effects of stratification and nongeostrophy. J. Atmos. Sci., 45, 32533268, https://doi.org/10.1175/1520-0469(1988)045<3253:SSOBIO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nikurashin, M., and R. Ferrari, 2013: Overturning circulation driven by breaking internal waves in the deep ocean. Geophys. Res. Lett., 40, 31333137, https://doi.org/10.1002/grl.50542.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Pedlosky, J., 1979: Geophysical Fluid Dynamics. Springer, 710 pp.

    • Crossref
    • Export Citation
  • Pedlosky, J., 2016: Baroclinic instability over topography: Unstable at any wave number. J. Mar. Res., 74, 119, https://doi.org/10.1357/002224016818377595.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Polzin, K. L., 1997: Spatial variability of turbulent mixing in the abyssal ocean. Science, 276, 9396, https://doi.org/10.1126/science.276.5309.93.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rhines, P. B., 1975: Waves and turbulence on a beta-plane. J. Fluid Mech., 69, 417443, https://doi.org/10.1017/S0022112075001504.

  • Ripa, P., 1999: On the validity of layered models of ocean dynamics and thermodynamics with reduced vertical resolution. Dyn. Atmos. Oceans, 29, 140, https://doi.org/10.1016/S0377-0265(98)00056-6.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ripa, P., 2001: Waves and resonance in free-boundary baroclinic instability. J. Fluid Mech., 428, 387408, https://doi.org/10.1017/S0022112000002731.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ruan, X., A. F. Thompson, M. M. Flexas, and J. Sprintall, 2017: Contribution of topographically generated submesoscale turbulence to Southern Ocean overturning. Nat. Geosci., 10, 840845, https://doi.org/10.1038/ngeo3053.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Sasaki, H., P. Klein, B. Qiu, and Y. Sasai, 2014: Impact of oceanic-scale interactions on the seasonal modulation of ocean dynamics by the atmosphere. Nat. Commun., 5, 5636, https://doi.org/10.1038/ncomms6636.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Solodoch, A., A. L. Stewart, and J. C. McWilliams, 2016: Baroclinic instability of axially symmetric flow over sloping bathymetry. J. Fluid Mech., 799, 265296, https://doi.org/10.1017/jfm.2016.376.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Spall, M. A., 2004: Boundary currents and watermass transformation in marginal seas. J. Phys. Oceanogr., 34, 11971213, https://doi.org/10.1175/1520-0485(2004)034<1197:BCAWTI>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • St. Laurent, L. C., J. M. Toole, and R. W. Schmitt, 2001: Buoyancy forcing by turbulence above rough topography in the abyssal Brazil basin. J. Phys. Oceanogr., 31, 34763495, https://doi.org/10.1175/1520-0485(2001)031<3476:BFBTAR>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stahr, F. R., and T. B. Sanford, 1999: Transport and bottom boundary layer observations of the North Atlantic Deep Western Boundary Current at the Blake Outer Ridge. Deep-Sea Res. II, 46, 205243, https://doi.org/10.1016/S0967-0645(98)00101-5.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stone, P. H., 1966: On non-geostrophic baroclinic stability. J. Atmos. Sci., 23, 390400, https://doi.org/10.1175/1520-0469(1966)023<0390:ONGBS>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stone, P. H., 1971: Baroclinic stability under non-hydrostatic conditions. J. Fluid Mech., 45, 659671, https://doi.org/10.1017/S0022112071000260.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Su, Z., J. Wang, P. Klein, A. F. Thompson, and D. Menemenlis, 2018: Ocean submesoscales as a key component of the global heat budget. Nat. Commun., 9, 775, https://doi.org/10.1038/s41467-018-02983-w.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taylor, J. R., 2016: Turbulent mixing, restratification, and phytoplankton growth at a submesoscale eddy: Submesoscales, turbulence, and phytoplankton growth. Geophys. Res. Lett., 43, 57845792, https://doi.org/10.1002/2016GL069106.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taylor, J. R., and R. Ferrari, 2009: On the equilibration of a symmetrically unstable front via a secondary shear instability. J. Fluid Mech., 622, 103113, https://doi.org/10.1017/S0022112008005272.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taylor, J. R., and R. Ferrari, 2010: Buoyancy and wind-driven convection at mixed layer density fronts. J. Phys. Oceanogr., 40, 12221242, https://doi.org/10.1175/2010JPO4365.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., 2005: Destruction of potential vorticity by winds. J. Phys. Oceanogr., 35, 24572466, https://doi.org/10.1175/JPO2830.1.

  • Thomas, L. N., and P. B. Rhines, 2002: Nonlinear stratified spin-up. J. Fluid Mech., 473, 211244, https://doi.org/10.1017/S0022112002002367.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., and J. R. Taylor, 2010: Reduction of the usable wind-work on the general circulation by forced symmetric instability. Geophys. Res. Lett., 37, L18606, https://doi.org/10.1029/2010GL044680.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., A. Tandon, and A. Mahadevan, 2008: Submesoscale processes and dynamics. Ocean Modeling in an Eddying Regime, Geophys. Monogr., Vol. 177, Amer. Geophys. Union, 17–38, https://doi.org/10.1029/177GM04.

    • Crossref
    • Export Citation
  • Thomas, L. N., J. R. Taylor, R. Ferrari, and T. M. Joyce, 2013: Symmetric instability in the Gulf Stream. Deep-Sea Res. II, 91, 96110, https://doi.org/10.1016/j.dsr2.2013.02.025.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Toole, J., R. Curry, T. Joyce, M. McCartney, and B. Peña-Molino, 2011: Transport of the North Atlantic Deep Western Boundary Current about 39°N, 70°W: 2004–2008. Deep-Sea Res. II, 58, 17681780, https://doi.org/10.1016/j.dsr2.2010.10.058.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Umlauf, L., and H. Burchard, 2011: Diapycnal transport and mixing efficiency in stratified boundary layers near sloping topography. J. Phys. Oceanogr., 41, 329345, https://doi.org/10.1175/2010JPO4438.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Umlauf, L., H. Burchard, and K. Bolding, 2005: GOTM—Scientific documentation, version 3.2. Marine Science Rep. 63, Institut für Ostseeforschung Warnemünde, 274 pp., https://www.io-warnemuende.de/tl_files/forschung/meereswissenschaftliche-berichte/mebe63_2005-gotm.pdf.

  • Umlauf, L., W. D. Smyth, and J. N. Moum, 2015: Energetics of bottom Ekman layers during buoyancy arrest. J. Phys. Oceanogr., 45, 30993117, https://doi.org/10.1175/JPO-D-15-0041.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Vanneste, J., 1993: The Kelvin-Helmholtz instability in a non-geostrophic baroclinic unstable flow. Math. Comput. Model., 17, 149154, https://doi.org/10.1016/0895-7177(93)90099-K.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wenegrat, J. O., and M. J. McPhaden, 2016: Wind, waves, and fronts: Frictional effects in a generalized Ekman model. J. Phys. Oceanogr., 46, 371394, https://doi.org/10.1175/JPO-D-15-0162.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wenegrat, J. O., L. N. Thomas, J. Gula, and J. C. McWilliams, 2018: Effects of the submesoscale on the potential vorticity budget of ocean mode waters. J. Phys. Oceanogr., 48, 21412165, https://doi.org/10.1175/JPO-D-17-0219.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
1

Throughout this paper, we will use the term bottom boundary layer to denote any weakly stratified lower layer, a definition that includes both very well mixed layers adjacent to topography and the weakly stratified overlying layer that can appear in some solutions to the 1D dynamics (e.g., section 3).

2

Note that we choose to follow the notation of Stone (1966), and hence we denote the slope parameter using α and the nonhydrostatic parameter using δ (cf. Hetland 2017). As in Hetland (2017) our definition of the slope parameter is of the opposite sign to that used in Blumsack and Gierasch (1972) and other preceding work, such that α > 0 implies isopycnal and topographic slopes of opposite sign (see, e.g., Fig. 5).

Save