Baroclinic flow around planetary islands in a double gyre : a mechanism for cross-gyre flow

Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 1075-1086, doi:10.1175/2009JPO4375.1.A quasigeostrophic, two-layer model is used to study the baroclinic circulation around a thin, meridionally elongated island. The flow is driven by either buoyancy forcing or wind stress, each of whose structure would produce an antisymmetric double-gyre flow. The ocean bottom is flat. When the island partially straddles the intergyre boundary, fluid from one gyre is forced to flow into the other. The amount of the intergyre flow depends on the island constant, that is, the value of the geostrophic streamfunction on the island in each layer. That constant is calculated in a manner similar to earlier studies and is determined by the average, along the meridional length of the island, of the interior Sverdrup solution just to the east of the island. Explicit solutions are given for both buoyancy and wind-driven flows. The presence of an island of nonzero width requires the determination of the baroclinic streamfunction on the basin’s eastern boundary. The value of the boundary term is proportional to the island’s area. This adds a generally small additional baroclinic intergyre flow. In all cases, the intergyre flow produced by the island is not related to topographic steering of the flow but rather the pressure anomaly on the island as manifested by the barotropic and baroclinic island constants. The vertical structure of the flow around the island is a function of the parameterization of the vertical mixing in the problem and, in particular, the degree to which long baroclinic Rossby waves can traverse the basin before becoming thermally damped.This research was supported in part by NSF Grant OCE 0451086

Baroclinic instability over topography : unstable at any wave number

Author Posting. © Sears Foundation for Marine Research, 2016. This article is posted here by permission of Sears Foundation for Marine Research for personal use, not for redistribution. The definitive version was published in Journal of Marine Research 74 (2016): 1-19, doi:10.1357/002224016818377595.The instability of an inviscid, baroclinic vertically sheared current of uniform potential vorticity, flowing along a uniform topographic slope, becomes linearly unstable at all wave numbers if the flow is in the direction of propagation of topographic waves. The parameter region of instability in the plane of scaled topographic slope versus wave number then extends to arbitrarily large wave numbers at large slopes. The weakly nonlinear treatment of the problem reveals the existence of a nonlinear enhancement of the instability close to one of the two boundaries of this parametrically narrow unstable region. Because the domain of instability becomes exponentially narrow for large wave numbers, it is unclear how applicable the results of the asymptotic, weakly nonlinear theory are given that it must be limited to a region of small supercriticality. This question is pursued in that parameter domain through the use of a truncated model in which the approximations of weakly nonlinear theory are avoided. This more complex model demonstrates that the linearly most unstable wave in the narrow wedge in parameter space is nonlinearly stable and that the region of nonlinear destabilization is limited to a tiny region near one of the critical curves rendering both the linear and nonlinear growth essentially negligible

The response of a weakly stratified layer to buoyancy forcing

Author Posting. © American Meteorological Society, 2009. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 39 (2009): 1060-1068, doi:10.1175/2008JPO3996.1.The response of a weakly stratified layer of fluid to a surface cooling distribution is investigated with linear theory in an attempt to clarify recent numerical results concerning the sinking of cooled water in polar ocean boundary currents. A channel of fluid is forced at the surface by a cooling distribution that varies in the down-channel as well as the cross-channel directions. The resulting geostrophic flow in the central region of the channel impinges on its boundaries, and regions of strong downwelling are observed. For the parameters of the problem investigated, the downwelling occurs in a classical Stewartson layer but the forcing of the layer leads to an unusual relation with the interior flow, which is forced to satisfy the thermal condition on the boundary while the geostrophic normal flow in the interior is brought to rest in the boundary layer. As a consequence of the layer’s dynamics, the resulting long-channel flow exhibits a nonmonotonic approach to the interior flow, and the strongest vertical velocities are limited to the boundary layer whose scale is so small that numerical models resolve the region only with great difficulty. The analytical model presented here is able to reproduce key features of the previous nonlinear numerical calculations.This research was supported in part by NSF Grant OCE 0451086

The nonlinear downstream development of baroclinic instability

Author Posting. © Sears Foundation for Marine Research, 2011. This article is posted here by permission of Sears Foundation for Marine Research for personal use, not for redistribution. The definitive version was published in Journal of Marine Research 69 (2011): 705-722, doi:10.1357/002224011799849363.The downstream development in both space and time of baroclinic instability is studied in a nonlinear channel model on the f-plane. The model allows the development of the instability to be expressed on space and time scales that are long compared to the growth rates and wavelengths of the most unstable wave. The unstable system is forced by time-varying boundary conditions at the origin of the channel and so serves as a conceptual model for the development of fluctuations in currents like the Gulf Stream and Kuroshio downstream of their separation points from their respective western boundaries. The theory is developed for both substantially dissipative systems as well as weakly dissipative systems for which the viscous decay time is of the order of the advective time in the former case and the growth time in the latter case. In the first case a first order equation in time leads to a hyperbolic system for which exact solutions are found in the case of monochromatic forcing. For a finite bandwidth the governing equations are nonlinear and parabolic and could be put in the form of the Real Ginzburg Landau equation first developed by Newell and Whitehead (1969) and Segel (1969) although we show the equation is not pertinent to the downstream development problem. When the dissipation is small a third order system of partial differential equations is obtained. For steady states the system supports chaotic behavior along the characteristics. This produces for the-time dependent problem new features, principally a strong focusing of amplitude in the regions behind the advancing front and the appearance of what might be called “chaotic shocks.“This research was supported in part by NSF Grant OCE 0925061

A note on the western intensification of the oceanic circulation

The purpose of this note is to provide a simple physical explanation for the westward intensification of the oceanic circulation found in the several dynamically different existing theoretical models (e.g. Stommel 1948, Carrier and Robinson 1962). It has the advantage of showing why the western oceanic boundary is singled out as the boundary-layer region that closes the interior Sverdrup solution

An inertial model of the interaction of Ekman layers and planetary islands

Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 1398–1406, doi:10.1175/JPO-D-13-028.1.An adiabatic, inertial, and quasigeostrophic model is used to discuss the interaction of surface Ekman transport with an island. The theory extends the recent work of Spall and Pedlosky to include an analytical and nonlinear model for the interaction. The presence of an island that interrupts a uniform Ekman layer transport raises interesting questions about the resulting circulation. The consequential upwelling around the island can lead to a local intake of fluid from the geostrophic region beneath the Ekman layer or to a more complex flow around the island in which the fluid entering the Ekman layer on one portion of the island's perimeter is replaced by a flow along the island's boundary from a downwelling region located elsewhere on the island. This becomes especially pertinent when the flow is quasigeostrophic and adiabatic. The oncoming geostrophic flow that balances the offshore Ekman flux is largely diverted around the island, and the Ekman flux is fed by a transfer of fluid from the western to the eastern side of the island. As opposed to the linear, dissipative model described earlier, this transfer takes place even in the absence of a topographic skirt around the island. The principal effect of topography in the inertial model is to introduce an asymmetry between the circulation on the northern and southern sides of the island. The quasigeostrophic model allows a simple solution to the model problem with topography and yet the resulting three-dimensional circulation is surprisingly complex with streamlines connecting each side of the island.This research was supported in part by NSF Grant OCE Grant 0925061

A note on interior pathways in the meridional overturning circulation

Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 643-646, doi:10.1175/JPO-D-17-0240.1.A simple oceanic model is presented for source–sink flow on the β plane to discuss the pathways from source to sink when transport boundary layers have large enough Reynolds numbers to be inertial in their dynamics. A representation of the flow as a Fofonoff gyre, suggested by prior work on inertial boundary layers and eddy-driven circulations in two-dimensional turbulent flows, indicates that even when the source and sink are aligned along the same western boundary the flow must intrude deep into the interior before exiting at the sink. The existence of interior pathways for the flow is thus an intrinsic property of an inertial circulation and is not dependent on particular geographical basin geometry.2018-09-1

Symmetric instability of cross-stream varying currents

Author Posting. © Sears Foundation for Marine Research, 2014. This article is posted here by permission of Sears Foundation for Marine Research for personal use, not for redistribution. The definitive version was published in Journal of Marine Research 72 (2014): 31-45, doi:10.1357/002224014812655206.The symmetric instability of a simple shear flow in which the velocity is a linear function of the vertical coordinate but which varies slowly in the cross-stream direction is studied using an asymptotic analytical method. Explicit analytical solutions are found for the evolution of the envelope of the developing linear instability. Although the problem with no lateral variation yields cell-like instabilities growing in place, the lateral variation of the shear produces time dependence and cross-stream propagation of the envelope and accompanying cells. A similar solution is derived for the case of laterally uniform shear in a current whose depth slowly varies exponentially in the cross-stream direction producing similar time dependence to the otherwise stationary cell pattern

The effect of beta on the downstream development of unstable, chaotic baroclinic waves

Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 49(9), (2019): 2337-2343, doi:10.1175/JPO-D-19-0097.1.The weakly unstable, two-layer model of baroclinic instability is studied in a configuration in which the flow is perturbed at the inflow section of a channel by a slow and periodic perturbation. In a parameter regime where the governing equation would be the Lorenz equations for chaos if the development occurs only in time, the solution behavior becomes considerably more complex as a function of time and downstream coordinate. In the absence of the beta effect it has earlier been shown that the chaotic behavior along characteristics renders the solution nearly discontinuous in the slow downstream coordinate of the asymptotic model. The additional presence of the beta effect, although expunging the chaos for large enough values of the beta parameter, also provides an additional mechanism for abrupt spatial change.2020-02-2

A necessary condition for the existence of an inertial boundary layer in a baroclinic ocean

A criterion for the existence of an inertial boundary current in a baroclinic ocean is derived. It is shown that the effect of topographical variations can play a decisive role in determining whether an inertial boundary current can exist, even if the velocity field is highly baroclinic