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

The nonlinear dynamics of time-dependent subcritical baroclinic currents

Author Posting. © American Meteorological Society, 2007. 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 37 (2007): 1001-1021, doi:10.1175/jpo3034.1.The nonlinear dynamics of baroclinically unstable waves in a time-dependent zonal shear flow is considered in the framework of the two-layer Phillips model on the beta plane. In most cases considered in this study the amplitude of the shear is well below the critical value of the steady shear version of the model. Nevertheless, the time-dependent problem in which the shear oscillates periodically is unstable, and the unstable waves grow to substantial amplitudes, in some cases with strongly nonlinear and turbulent characteristics. For very small values of the shear amplitude in the presence of dissipation an analytical, asymptotic theory predicts a self-sustained wave whose amplitude undergoes a nonlinear oscillation whose period is amplitude dependent. There is a sensitive amplitude dependence of the wave on the frequency of the oscillating shear when the shear amplitude is small. This behavior is also found in a truncated model of the dynamics, and that model is used to examine larger shear amplitudes. When there is a mean value of the shear in addition to the oscillating component, but such that the total shear is still subcritical, the resulting nonlinear states exhibit a rectified horizontal buoyancy flux with a nonzero time average as a result of the instability of the oscillating shear. For higher, still subcritical, values of the shear, a symmetry breaking is detected in which a second cross-stream mode is generated through an instability of the unstable wave although this second mode would by itself be stable on the basic time-dependent current. For shear values that are substantially subcritical but of order of the critical shear, calculations with a full quasigeostrophic numerical model reveal a turbulent flow generated by the instability. If the beta effect is disregarded, the inviscid, linear problem is formally stable. However, calculations show that a small degree of nonlinearity is enough to destabilize the flow, leading to large amplitude vacillations and turbulence. When the most unstable wave is not the longest wave in the system, a cascade up scale to longer waves is observed. Indeed, this classically subcritical flow shows most of the qualitative character of a strongly supercritical flow. This result supports previous suggestions of the important role of background time dependence in maintaining the atmospheric and oceanic synoptic eddy field.GRF was supported by NSF Grant OCE-0137023, and JP was supported by NSF Grant OCE- 9901654

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

Time-dependent response to cooling in a beta-plane basin

Author Posting. © American Meteorological Society, 2006. 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 36 (2006): 2185-2198, doi:10.1175/JPO2967.1.The time-dependent response of an ocean basin to the imposition of cooling (or heating) is examined in the context of a quasigeostrophic, two-layer model on the beta plane. The focus is on the structure and magnitude of the vertical motion and its response to both a switch-on forcing and a periodic forcing. The model employed is a time-dependent version of an earlier model used to discuss the intensification of sinking in the region of the western boundary current. The height of the interface of the two-layer model serves as an analog of temperature, and the vertical velocity at the interface consists of a cross-isopycnal velocity modeled in terms of a relaxation to a prescribed interface height, an adiabatic representation of eddy thickness fluxes parameterized as lateral diffusion of thickness, and the local vertical motion of the interface itself. The presence of time dependence adds additional dynamical features to the problem, in particular the emergence of low-frequency, weakly damped Rossby basin modes. If the buoyancy forcing is zonally uniform the basin responds to a switch-on of the forcing by coming into steady-state equilibrium after the passage of a single baroclinic Rossby wave. If the forcing is nonuniform in the zonal direction, a sequence of Rossby basin modes is excited and their decay is required before the basin achieves a steady state. For reasonable parameter values the boundary layers, in which both horizontal and vertical circulations are closed, are quasi-steady and respond to the instantaneous state of the interior. As in the steady problem the flow is sensitive to small nonquasigeostrophic mass fluxes across the perimeter of the basin. These fluxes generally excite basin modes as well. The basin modes will also be weakly excited if the beta-plane approximation is relaxed. The response to periodic forcing is also examined, and the sensitivity of the response to the structure of the forcing is similar to the switch-on problem.This research was supported in part by NSF Grant OCE-9901654