Eddy Potential Vorticity Flux Research Articles

Simulation of eddy-driven deep circulation in the East/Japan Sea by using a three-layer model with wind, throughflow and deep water formation forcings

A three-layer model was used to investigate the impact of grid resolution, bottom topography, wind, throughflow (the Tsushima Warm Current) and deep water formation (DWF) forcing on deep circulation in the East/Japan Sea. High grid resolution (1/36°) and real bottom topography are essential to capture the observed features of surface and deep circulation. An increase of grid resolution from 1/12° to 1/36°, wind and DWF forcing dramatically increase the current velocity and eddy kinetic energy (EKE) of the lower-layer. Of the three forcings, the influence of wind forcing is the most powerful on the abyssal circulation, followed by the order of the throughflow and DWF forcing, respectively. DWF forcing does not largely change the horizontal circulation pattern of the lower-layer in the 1/36° resolution model, but it increases the current velocity by 33.9%. Model results indicate that surface EKE is larger in the southern region than in the northern region of the East/Japan Sea, while deep EKE is larger in the northern area. These results are in good agreement with EKE distributions obtained from the surface drifter and the ARGO float trajectories. The EKE maxima of the lower-layer is geographically concentrated on the bottom slope, indicating that the strong circulation in the lower-layer is due to the strong eddy activity on the slope. The flux calculation shows that the primary eddy fluxes of potential vorticity driving deep circulation are layer thickness flux (LAY) and relative vorticity flux (REL).

Comparison of the Effect of Parameterized Eddy Fluxes of Thickness and Potential Vorticity

AbstractParameterization of mesoscale eddies is an important problem of modern ocean dynamics and modeling. The most widely used scheme is the so-called Gent–McWilliams parameterization, which describes the eddy-induced transport of tracers, including temperature, density, and isopycnal thickness (TH). An alternative scheme, proposed by Green and Welander, deals with parameterizing eddy fluxes of potential vorticity (PV). Many recent studies propose using it, for it includes the effect of eddy Reynolds stresses that may influence mean flows. These two schemes are compared in the simplest configuration of two-layer quasigeostrophic channel flow, which enables analytical solutions for zonal-mean fields. It is shown how the parameterizations shape the zonally averaged zonal velocity profiles, with special attention paid to the role of the Reynolds stresses and momentum conservation. The zonally averaged zonal velocity profiles are sensitive to the amplitude and profiles of TH and PV diffusivities. For small enough diffusivities the TH parameterization may lead to solutions resembling those for the PV parameterization if it uses the diffusivity of the latter; that is, it may mimic the impact of the Reynolds stresses on the mean flow.

Rapid Eddy-Induced Modification of Subtropical Mode Water during the Kuroshio Extension System Study

Abstract From 2004 to 2006 an observational array of current- and pressure-recording inverted echo sounders (CPIES) were deployed as part of the Kuroshio Extension (KEx) System Study (KESS). KESS observed a transition from a weakly meandering (“stable”) to strongly meandering (“unstable”) state (Qiu and Chen). As the KEx made this transition, potential vorticity (PV) observed within the southern recirculation gyre (SRG) rapidly increased from January to July 2005. In this study, the authors diagnose eddy PV fluxes (EPVFs) in isentropic coordinates within the subtropical mode water (STMW) layer from the CPIES data to determine the role of mesoscale eddies in this rapid increase of PV. The rapid increase in PV within the SRG coincided with enhanced cross-front EPVFs and eddy PV flux convergence upstream of a mean trough in the KEx path and adjacent to the SRG. The enhanced cross-front EPVFs were the result of the formation of a cold-core ring (CCR) and the interaction of the jet with a preexisting CCR. Eddy diffusivities are diagnosed for the unstable regime with values that range from 100 to 2000 m2 s−1. The high eddy diffusivities during the unstable regime reflect the nature of mesoscale CCR formation and CCR–jet interaction as efficient mechanisms for stirring and mixing high PV waters from the north side of the KEx into the low PV waters of the SRG where STMW resides. This mechanism for cross-frontal exchange can explain observed increases in the STMW PV in the SRG over the 16 months of KESS observations.

On the formation and maintenance of the stratospheric surf zone as inferred from the zonally averaged potential vorticity distribution

The diagnostic relation between eddy potential vorticity flux and Eliassen–Palm flux convergence in the transformed Eulerian mean shallow‐water system is used to infer the flux convergence needed to establish a stratospheric surf zone comprising a single region of perfectly homogenized potential vorticity, or to maintain the surf zone in steady state against the restoring effect of radiative relaxation. In the transient case, and when wave breaking is assumed to mix potential vorticity on time‐scales shorter than the radiative time‐scale, the required flux convergence is an order of magnitude larger than that required to maintain the flow in steady state away from radiative equilibrium.

The Role of Subtropical Irreversible PV Mixing in the Zonal Mean Circulation Response to Global Warming–Like Thermal Forcing

Abstract The atmospheric circulation response to the global warming–like tropical upper tropospheric heating is revisited using a dry atmospheric general circulation model (AGCM) in light of new diagnostics based on the concept of finite-amplitude wave activity (FAWA) on equivalent latitude. For a given tropical heating profile, the linear Wentzel–Kramers–Brillouin (WKB) wave refraction analysis sometimes gives a very different and even opposite prediction of the eddy momentum flux response to that of the actual full model simulation, exposing the limitation of the traditional linear approach in understanding the full dynamics of the atmospheric response under global warming. The implementation of the FAWA diagnostics reveals that in response to the upper tropospheric heating, effective diffusivity—a measure of the mixing efficiency—increases and advances upward and poleward in the subtropics and the resultant enhancement and the poleward encroachment of eddy potential vorticity mixing leads to a poleward displaced potential vorticity (PV) gradient peak in the upper troposphere. The anomalous eddy PV flux, in balance with the PV dissipation, gives rise to a poleward shift in the eddy-driven jet and eddy-driven mean meridional circulation. Sensitivity experiments show that these irreversible dissipation processes in the upper troposphere are robust, regardless of the width of the tropical heating.

Integral Constraints for Momentum and Energy in Zonal Flows with Parameterized Potential Vorticity Fluxes: Governing Parameters

Abstract Integral constraints for momentum and energy impose restrictions on parameterizations of eddy potential vorticity (PV) fluxes. The impact of these constraints is studied for a wind-forced quasigeostrophic two-layer zonal channel model with variable bottom topography. The presence of a small parameter, given by the ratio of Rossby radius to the width of the channel, makes it possible to find an analytical/asymptotic solution for the zonally and time-averaged flow, given diffusive parameterizations for the eddy PV fluxes. This solution, when substituted in the constraints, leads to nontrivial explicit restrictions on diffusivities. The system is characterized by four dimensionless governing parameters with a clear physical interpretation. The bottom form stress, the major term balancing the external force of wind stress, depends on the governing parameters and fundamentally modifies the restrictions compared to the flat bottom case. While the analytical solution bears an illustrative character, it helps to see certain nontrivial connections in the system that will be useful in the analysis of more complicated models of ocean circulation. A numerical solution supports the analytical study and confirms that the presence of topography strongly modifies the eddy fluxes.

The Impact of Baroclinic Eddy Feedback on the Persistence of Jet Variability in the Two-Layer Model

Abstract Although it is well known that the persistence of extratropical jet shifts is enhanced by a positive eddy feedback, the dynamics of this feedback is still debated. Two types of mechanisms have been proposed: barotropic mechanisms rely on changes in upper-level propagation and baroclinic mechanisms rely on the coupling between barotropic and baroclinic flow. Recent studies have suggested that barotropic models can capture key aspects of the observed jet variability but the role of baroclinic dynamics has been less explored. This study investigates the temporal relations between barotropic and baroclinic anomalies and their eddy forcings during the internal variability of the simple two-layer quasigeostrophic model. A large correlation is found between barotropic and baroclinic anomalies and between the meridional and vertical components of the Eliassen–Palm divergence, especially at low frequency. The low-frequency variability is consistent with the baroclinic mechanism: persistent upper-level eddy momentum convergence is associated with (and precedes) persistent anomalies in the poleward eddy heat flux. In contrast, at high frequency, poleward heat flux anomalies are associated with eddy momentum divergence aloft and both eddy forcings have same-sign contributions to the upper-level eddy potential vorticity (PV) flux. In this limit the eddy PV flux is associated with wave activity transience as effective diffusivity is too small to dissipate the wave–mean flow interaction term. The large correlation between barotropic and baroclinic anomalies implies that the low-frequency variability of barotropic flow may be affected by thermal damping when this damping is sufficiently strong. For example, zonal index persistence drops drastically in our model when baroclinicity shifts are prevented by strong thermal restoration.

Up-gradient eddy fluxes of potential vorticity near the subtropical jet

The role of eddy fluxes in the general circulation is often approached by treating eddies as (macro)turbulence. In this approach, the eddies act to diffuse certain quasiconservative quantities, such as potential vorticity (PV), along isentropic surfaces in the free atmosphere. The eddy fluxes are determined primarily by the eddy diffusivities and are necessarily down-gradient of the basic state PV field. Support for the (macro)turbulence approach stems from the fact that the eddy fluxes of PV in the free atmosphere are generally down-gradient in the long-term mean. Here we call attention to a pronounced and significant region of upgradient eddy PV fluxes on the poleward flank of the jet core in both hemispheres. The region of up-gradient (i.e., notionally “antidiffusive”) eddy PV fluxes is most pronounced during the winter and spring seasons and partially contradicts the turbulence approach described above. Analyses of the PV variance (potential enstrophy) budget suggest that the up-gradient PV fluxes represent local wave decay and are maintained by poleward fluxes of PV variance. Finite-amplitude effects thus represent leading order contributions to the PV variance budget, whereas dissipation is only of secondary importance locally. The appearance of up-gradient PV fluxes in the long-term mean is associated with the poleward shift of the jet—and thus the region of wave decay relative to wave growth—following wave-breaking events.

Dynamical and Physical Processes Leading to Tropical Cyclone Intensification under Upper-Level Trough Forcing

Abstract The rapid intensification of Tropical Cyclone (TC) Dora (2007, southwest Indian Ocean) under upper-level trough forcing is investigated. TC–trough interaction is simulated using a limited-area operational numerical weather prediction model. The interaction between the storm and the trough involves a coupled evolution of vertical wind shear and binary vortex interaction in the horizontal and vertical dimensions. The three-dimensional potential vorticity structure associated with the trough undergoes strong deformation as it approaches the storm. Potential vorticity (PV) is advected toward the tropical cyclone core over a thick layer from 200 to 500 hPa while the TC upper-level flow turns cyclonic from the continuous import of angular momentum. It is found that vortex intensification first occurs inside the eyewall and results from PV superposition in the thick aforementioned layer. The main pathway to further storm intensification is associated with secondary eyewall formation triggered by external forcing. Eddy angular momentum convergence and eddy PV fluxes are responsible for spinning up an outer eyewall over the entire troposphere, while spindown is observed within the primary eyewall. The 8-km-resolution model is able to reproduce the main features of the eyewall replacement cycle observed for TC Dora. The outer eyewall intensifies further through mean vertical advection under dynamically forced upward motion. The processes are illustrated and quantified using various diagnostics.

The Vertical Structure of the Eddy Diffusivity and the Equilibration of the Extratropical Atmosphere

Abstract Observations suggest that the time- and zonal-mean state of the extratropical atmosphere adjusts itself such that the so-called “criticality parameter” (which relates the vertical stratification to the horizontal temperature gradient) is close to one. T. Schneider has argued that the criticality parameter is kept near one by a constraint on the zonal momentum budget in primitive equations. The constraint relies on a diffusive closure for the eddy flux of potential vorticity (PV) with an eddy diffusivity that is approximately constant in the vertical. The diffusive closure for the eddy PV flux, however, depends crucially on the definition of averages along isentropes that intersect the surface. It is argued that the definition favored by Schneider results in eddy PV fluxes whose physical interpretation is unclear and that do not satisfy the proposed closure in numerical simulations. An alternative definition, first proposed by T.-Y. Koh and R. A. Plumb, is preferred. A diffusive closure for the eddy PV flux under this definition is supported by analysis of the PV variance budget and can be used to close the near-surface zonal momentum budget in idealized numerical simulations. Following this approach, it is shown that O(1) criticalities are obtained if the eddy diffusivity decays from its surface value to about zero over the depth of the troposphere, which is likely to be the case in Earth’s atmosphere. Large criticality parameters, however, are possible if the eddy diffusivity decays only weakly in the vertical, consistent with results from quasigeostrophic models. This provides theoretical support for recent numerical studies that have found supercritical mean states in primitive equation models.

Influence of Bottom Topography on Integral Constraints in Zonal Flows with Parameterized Potential Vorticity Fluxes

Abstract An integral constraint for eddy fluxes of potential vorticity (PV), corresponding to global momentum conservation, is applied to two-layer zonal quasigeostrophic channel flow. This constraint must be satisfied for any type of parameterization of eddy PV fluxes. Bottom topography strongly influences the integral constraint compared to a flat bottom channel. An analytical solution for the mean flow solution has been found by using asymptotic expansion in a small parameter, which is the ratio of the Rossby radius to the meridional extent of the channel. Applying the integral constraint to this solution, one can find restrictions for eddy PV transfer coefficients that relate the eddy fluxes of PV to the mean flow. These restrictions strongly deviate from restrictions for the channel with flat bottom topography.

A Framework for Parameterizing Eddy Potential Vorticity Fluxes

AbstractA framework for parameterizing eddy potential vorticity fluxes is developed that is consistent with conservation of energy and momentum while retaining the symmetries of the original eddy flux. The framework involves rewriting the residual-mean eddy force, or equivalently the eddy potential vorticity flux, as the divergence of an eddy stress tensor. A norm of this tensor is bounded by the eddy energy, allowing the components of the stress tensor to be rewritten in terms of the eddy energy and nondimensional parameters describing the mean shape and orientation of the eddies. If a prognostic equation is solved for the eddy energy, the remaining unknowns are nondimensional and bounded in magnitude by unity. Moreover, these nondimensional geometric parameters have strong connections with classical stability theory. When applied to the Eady problem, it is shown that the new framework preserves the functional form of the Eady growth rate for linear instability. Moreover, in the limit in which Reynolds stresses are neglected, the framework reduces to a Gent and McWilliams type of eddy closure where the eddy diffusivity can be interpreted as the form proposed by Visbeck et al. Simulations of three-layer wind-driven gyres are used to diagnose the eddy shape and orientations in fully developed geostrophic turbulence. These fields are found to have large-scale structure that appears related to the structure of the mean flow. The eddy energy sets the magnitude of the eddy stress tensor and hence the eddy potential vorticity fluxes. Possible extensions of the framework to ensure potential vorticity is mixed on average are discussed.

Role of the Rossby Waves in the Broadening of an Eastward Jet

Abstract To investigate the effect of the Rossby waves on an eastward jet such as the Kuroshio or Gulf Stream Extensions, a series of numerical experiments is conducted using a primitive equation model. In these experiments, an inflow and an outflow imposed on the western and eastern boundaries drive an unstable narrow jet and a broad interior flow in the western and eastern regions of the model domain, respectively. The barotropic Rossby waves are radiated from the transient region between the two regions. The eddy potential vorticity flux by the waves tends to compensate for the difference in the mean potential vorticity along mean streamlines between both sides of the transient region. Instability of the jet is insufficient for this compensation and weakens the mean potential vorticity gradient too much. Moreover, as the potential vorticity of the outflow is increased, the Rossby waves are intensified in order to compensate for the increase in the difference in the mean potential vorticity. These features strongly suggest that the Rossby waves are substantial in matching a jet with an interior flow. The speed of the waves and properties of eddies in recirculations of the jet are consistent with a two-layer analytic model, which indicates that the Rossby waves are radiated from eddies in recirculations. These eddies as well as the Rossby waves increase in amplitude with the transport of the recirculation near the surface presumably because of mean advection. Therefore, the mean potential vorticity of the interior flow, the intensity of the Rossby waves, and the transport of the recirculation change consistently with one another.

Properties of Steady Geostrophic Turbulence with Isopycnal Outcropping

Abstract High-resolution simulations of β-channel, zonal-jet, baroclinic turbulence with a three-dimensional quasigeostrophic (QG) model including surface potential vorticity (PV) are analyzed with emphasis on the competing role of interior and surface PV (associated with isopycnal outcropping). Two distinct regimes are considered: a Phillips case, where the PV gradient changes sign twice in the interior, and a Charney case, where the PV gradient changes sign in the interior and at the surface. The Phillips case is typical of the simplified turbulence test beds that have been widely used to investigate the effect of ocean eddies on ocean tracer distribution and fluxes. The Charney case shares many similarities with recent high-resolution primitive equation simulations. The main difference between the two regimes is indeed an energization of submesoscale turbulence near the surface. The energy cycle is analyzed in the (k, z) plane, where k is the horizontal wavenumber. In the two regimes, the large-scale buoyancy forcing is the primary source of mechanical energy. It sustains an energy cycle in which baroclinic instability converts more available potential energy (APE) to kinetic energy (KE) than the APE directly injected by the forcing. This is due to a conversion of KE to APE at the scale of arrest. All the KE is dissipated at the bottom at large scales, in the limit of infinite resolution and despite the submesoscales energizing in the Charney case. The eddy PV flux is largest at the scale of arrest in both cases. The eddy diffusivity is very smooth but highly nonuniform. The eddy-induced circulation acts to flatten the mean isopycnals in both cases.

The Position of the Midlatitude Storm Track and Eddy-Driven Westerlies in Aquaplanet AGCMs

Abstract The sensitivity of the midlatitude storm track and eddy-driven wind to the sea surface temperature (SST) boundary forcing is studied over a wide range of perturbations using both simple and comprehensive general circulation models over aquaplanet lower boundary conditions. Under the single-jet circulation regime similar to the conditions of the present climate in the Northern Hemisphere winter or the Southern Hemisphere summer, the eddy-driven jet shifts monotonically poleward with both the global mean and the equator-to-pole gradient of the SST. The eddy-driven jet can have a reverse relationship to the gradient if it is well separated from the subtropical jet and Hadley cell boundary in a double-jet circulation regime. A simple scaling is put forward to interpret the simulated sensitivity of the storm-track/eddy-driven westerly wind position within the single-jet regime in both models. The rationale for the scaling is based on the notion that the wave activity flux can propagate horizontally away from the source region, resulting in a broader distribution of eddy potential vorticity (PV) flux in the upper troposphere than that of the flux in the opposite direction in the lower troposphere. As a consequence, the position of the maximum of the eddy-driven westerlies tends to be controlled by the profile of the relatively sharp-peaked low-level PV flux, which is dominated by the eddy heat flux component of the Eliassen–Palm (EP) flux. Thus, the position of the eddy-driven surface westerlies may be inferred from the vertical EP flux coming out of the lower troposphere. The vertical EP flux can be parameterized by a measure of baroclinicity, whose latitudinal variations show a linear relationship with the meridional displacement of the eddy-driven westerlies and the storm track. This relationship still holds well within the single-jet regime, even when only the variation of static stability is taken into consideration in estimating the baroclinicity (the temperature gradient component of which is fixed). To the extent that the static stability is deterministically constrained by and hence can be predicted from the given SST conditions through a moist scaling for the midlatitude stratification, one may, given SST perturbations, predict which way the storm track and eddy-driven wind should shift with respect to a chosen reference climate state. The resultant anomaly-wise scaling turns out to be valid for both the idealized and comprehensive models, regardless of the details in the model physics. By corollary, it can be argued that the poleward shift of storm track found in the global warming simulations by fully coupled climate models may be attributed, at least partially, to the increase in the subtropical and midlatitude static stability with global warming.

Finite-Amplitude Wave Activity and Mean Flow Adjustments in the Atmospheric General Circulation. Part I: Quasigeostrophic Theory and Analysis

Abstract A diagnostic relationship between finite-amplitude wave activity and the associated adiabatic adjustments to the zonal-mean zonal wind and temperature is developed in the quasigeostrophic (QG) framework and is applied to a 23-yr segment (1979–2001) of the 40-yr ECMWF Re-Analysis (ERA-40) data. Wave activity is defined in terms of an instantaneous areal displacement of QG potential vorticity (PV) from zonal symmetry. Unlike previous forms, the tendency of wave activity equals exactly the negative of the eddy PV flux (Eliassen–Palm flux divergence) in the conservative limit, even at finite amplitude. This allows one to integrate the transformed Eulerian mean (TEM) theory in time and quantify the departure (adiabatic adjustment) of the zonal-mean state from an eddy-free reference state in terms of the observed wave activity. The structure of wave activity identifies synoptic eddies in the extratropics and planetary waves in the high latitudes of winter-to-spring stratosphere. In addition, a thin layer of high wave activity is found at the top of the lowermost stratosphere (∼17 km) in the summer extratropics. The reference state is constructed by “zonalizing” the PV contours conservatively (preserving area) on the isobaric surface and by inverting the resultant PV gradient for the mean flow. The adjustment associated with wave activity depends on the assumed surface boundary condition for the reference state. With a no-slip condition, the observed zonal-mean temperature is on average ∼33 (90) K higher than the reference state in the troposphere (stratosphere) of the Arctic winter, while the zonal-mean zonal wind is ∼30 m s−1 slower in the upper stratosphere. Since the reference state filters out the advective eddy–mean flow interaction, it fluctuates less than the zonal-mean state, potentially improving the signal-to-noise ratio for climate diagnosis.

Transport and Mixing in the Extratropical Tropopause Region in a High-Vertical-Resolution GCM. Part II: Relative Importance of Large-Scale and Small-Scale Dynamics

Abstract The relative roles of atmospheric motions on various scales, from mesoscale to planetary scale, in transport and mixing in the extratropical tropopause region are investigated using a high-vertical-resolution general circulation model (GCM). The GCM with a vertical resolution of about 300 m explicitly represents the propagation and breaking of gravity waves and the induced transport and mixing. A downward control calculation shows that the Eliassen–Palm (E-P) flux of the gravity waves diverges and induces a mean equatorward flow in the extratropical tropopause region, which differs from the mean poleward flow induced by the convergence of large-scale E-P fluxes. The diffusion coefficients estimated from the eddy potential vorticity flux in tropopause-based coordinates reveal that isentropic motions diffuse air between 20 K below and 10 K above the tropopause from late autumn to early spring, while vertical mixing is strongly suppressed at around 10–15 K above the tropopause throughout the year. The isentropic mixing is mainly caused by planetary- and synoptic-scale motions, while small-scale motions with a horizontal scale of less than a few thousand kilometers largely affect the three-dimensional mixing just above the tropopause. Analysis of the gravity wave energy and atmospheric instability implies that the small-scale dynamics associated with the dissipation and saturation of gravity waves is a significant cause of the three-dimensional mixing just above the tropopause. A rapid increase in the static stability in the tropopause inversion layer is considered to play an important role in controlling the gravity wave activity around the tropopause.

Finite-Amplitude Equilibration of Baroclinic Waves on a Jet

Abstract A two-layer quasigeostrophic model is used to study the equilibration of baroclinic waves. In this model, if the background flow is relaxed toward a jetlike profile, a finite-amplitude baroclinic wave solution can be realized in both supercritical and subcritical regions of the model’s parameter space. Analyses of the model equations and numerical model calculations indicate that the finite-amplitude wave equilibration hinges on the breaking of Rossby waves before they reach their critical latitude. This “jetward” wave breaking results in an increase in the upper-layer wave generation and a reduction in the vertical phase tilt. This change in the phase tilt has a substantial impact on the Ekman pumping, as it weakens the damping on the lower-layer wave for some parameter settings and enables the Ekman pumping to serve as a source of wave growth at other settings. Together, these processes can account for the O(1)-amplitude wave equilibration. From a potential vorticity (PV) perspective, the wave breaking reduces the meridional scale of the upper-layer eddy PV flux, which destabilizes the mean flow. This is followed by a strengthening of the lower-layer eddy PV flux, which weakens the lower-layer PV gradient and constrains the growth of the lower-layer eddy PV. The same jetward wave breaking focuses the upper-layer PV flux toward the jet center where the upper-layer PV gradient is greatest. This results in an intensification of the upper-layer eddy PV relative to lower-layer eddy PV. Because of this large ratio, the upper-layer eddy PV plays the primary role in inducing the upper- and lower-layer eddy streamfunction fields, decreasing the vertical phase tilt. As a result, the Ekman pumping on the eddies is weakened, and for some parameter settings the Ekman pumping can even act as a wave source, contributing toward O(1)-amplitude wave equilibration. By reducing the horizontal shear of the zonal wind, the same wave breaking process weakens the barotropic decay, which also contributes to the wave amplification.

The Kuroshio Extension Northern Recirculation Gyre: Profiling Float Measurements and Forcing Mechanism

Abstract Middepth, time-mean circulation in the western North Pacific Ocean (28°–45°N, 140°–165°E) is investigated using drift information from the profiling floats deployed in the Kuroshio Extension System Study (KESS) and the International Argo programs. A well-defined, cyclonic recirculation gyre (RG) is found to exist north of the Kuroshio Extension jet, confined zonally between the Japan Trench (∼145°E) and the Shatsky Rise (∼156°E), and bordered to the north by the subarctic boundary along ∼40°N. This northern RG, which is simulated favorably in the eddy-resolving OGCM for the Earth Simulator (OFES) hindcast run model, has a maximum volume transport at 26.4 Sv across 159°E and its presence persists on the interannual and longer time scales. An examination of the time-mean x-momentum balance from the OFES hindcast run output reveals that horizontal convergence of Reynolds stresses works to accelerate both the eastward-flowing Kuroshio Extension jet and a westward mean flow north of the meandering jet. The fact that the northern RG is eddy driven is further confirmed by examining the turbulent Sverdrup balance, in which convergent eddy potential vorticity fluxes are found to induce the cyclonic RG across the background potential vorticity gradient field. For the strength of the simulated northern RG, the authors find the eddy dissipation effect to be important as well.

Formation of intrathermocline eddies at ocean fronts by wind-driven destruction of potential vorticity

A mechanism for the generation of intrathermocline eddies (ITEs) at wind-forced fronts is examined using a high resolution numerical simulation. Favorable conditions for ITE formation result at fronts forced by “down-front” winds, i.e. winds blowing in the direction of the frontal jet. Down-front winds exert frictional forces that reduce the potential vorticity (PV) within the surface boundary in the frontal outcrop, providing a source for the low-PV water that is the materia prima of ITEs. Meandering of the front drives vertical motions that subduct the low-PV water into the pycnocline, pooling it into the coherent anticyclonic vortex of a submesoscale ITE. As the fluid is subducted along the outcropping frontal isopycnal, the low-PV water, which at the surface is associated with strongly baroclinic flow, re-expresses itself as water with nearly zero absolute vorticity. This generation of strong anticyclonic vorticity results from the tilting of the horizontal vorticity of the frontal jet, not from vortex squashing. During the formation of the ITE, high-PV water from the pycnocline is upwelled alongside the subducting low-PV surface water. The positive correlation between the ITE’s velocity and PV fields results in an upward, along-isopycnal eddy PV flux that scales with the surface frictional PV flux driven by the wind. The relationship between the eddy and wind-induced frictional PV flux is nonlocal in time, as the eddy PV flux persists long after the wind forcing is shut off. The ITE’s PV flux affects the large-scale flow by driving an eddy-induced transport or bolus velocity down the outcropping isopycnal layer with a magnitude that scales with the Ekman velocity.