## 1. Introduction

A large portion of upper-ocean density fronts reside in the mixed layer, directly exposed to atmospheric forcing. Many fronts in the ocean, such as the subpolar fronts of the Gulf Stream and Kuroshio systems and fronts composing the Antarctic Circumpolar Current, are subject to powerful wind stress *τ** _{w}* blowing mostly in the direction of the frontal currents, yielding large, positive kinetic energy input

*τ**·*

_{w}**(Oort et al. 1994). Often the winds are accompanied by strong surface buoyancy fluxes. In the case of subpolar fronts in the Northern Hemisphere, these winds originate from the continent and in the winter bring cold dry air in contact with warm water of southern origin, extracting huge amounts of heat (up to 400 W m**

*u*^{−2}) from the surface ocean (da Silva et al. 1994). The combination of the atmospheric forcing and out-cropping of the pycnocline at these fronts permits the formation and subduction of intermediate and mode waters, water masses that retain the heat, salinity, potential vorticity, and chemical properties that were set in the mixed layer. Upwelling and interaction with the atmosphere of these water masses at some point during their journey through the ocean gyres may play an important role in the decadal variability of the ocean–atmosphere climate system (e.g., Latif and Barnett 1994). The first step in understanding this important process is to determine the physical mechanisms responsible for subduction and the intensification of ocean fronts in the presence of strong atmospheric forcing.

Many theoretical and modeling studies on the dynamics of ocean fronts (e.g., MacVean and Woods 1980; Bleck et al. 1988; Wang 1993; Spall 1995, 1997) are founded in the inviscid, adiabatic frontogenesis theory of Hoskins and Bretherton (1972), and therefore neglect mixing by surface fluxes. The frontal model of Hoskins and Bretherton (1972) describes the manner in which a confluent geostrophic flow (which for an oceanic application might represent the collision of western boundary currents at gyre boundaries or, on the mesoscale, confluence by eddy circulations) intensifies an initially weak baroclinic zone via its horizontal deformation field. This process involves the generation of an ageostrophic secondary circulation (ASC) whose convergent flow augments the confluence and leads to the formation of an infinitely strong front in a finite time. Not only is the ASC responsible for frontogenesis, but its downwelling branch determines the subduction rate. Therefore, understanding the mechanics of the ASC is crucial to frontal dynamics.

Ageostrophic secondary circulations arise at fronts to keep the alongfront geostrophic flow in geostrophic balance over subinertial time scales, as is required in the semigeostrophic approximation (Hoskins 1982). Advection of density and momentum by confluent flow tends to push the frontal jet out of a thermal-wind balance and hence induces an ASC whose spatial structure is governed by the *ω* equation (Hoskins et al. 1978). Like confluent flow, redistribution of momentum or buoyancy by small-scale turbulent mixing disrupts the geostrophic balance and, therefore, likewise drive a geostrophy-restoring ASC (Eliassen 1951).

As mentioned previously, winds over many fronts in the ocean are mostly down-front (i.e., they blow along the frontal jet). Down-front winds destabilize the water column as Ekman flow advects dense water over light water for this wind orientation. Therefore, at these fronts the winds as well as the destabilizing surface buoyancy flux will lead to gravitational instability. Mixing by gravitational instability redistributes buoyancy and can drive an ASC. If this ASC increases the horizontal density contrast across the front, the potential for the down-front winds to destabilize the water column is strengthened. This gives rise to the possibility of the following frontogenetic scenario: down-front winds drive convection, mixing buoyancy, and disrupting the geostrophic balance; the subsequent geostrophy-restoring ASC strengthens the front, further enhancing the wind-driven gravitational instability, buoyancy mixing, ASC, and frontal intensification. It will be shown in this paper using an analytic theory and high-resolution nonhydrostatic numerical simulations that this frontogenetic mechanism does indeed occur and is an efficient means for frontal intensification and subduction.

The outline of the paper is as follows. First, the semigeostrophic equations for a baroclinic zone forced by wind stress and surface buoyancy flux will be formulated. Next, the method and solution for the analytic theory of the wind/buoyancy-driven ASC and its frontogenetic effects are presented in section 3. Following this, the nonhydrostatic numerical simulations will be detailed in section 4. In section 5 the implications of wind-driven frontogenesis for the subpolar front of the Japan/East Sea forced by cold-air outbreaks is touched upon using observations. The paper is concluded in section 6.

## 2. Semigeostrophic dynamics of a wind-forced baroclinic zone

*u*,

_{g}*υ*). As often assumed in frontal studies, the alongfront length scale is supposed to be much larger than the cross-front length scale so that the geostrophic flow is nearly two dimensional, dominated by the alongfront flow (which is arbitrarily chosen to be oriented in the zonal direction); that is,

_{g}*u*≫

_{g}*υ*. The properties of the baroclinic zone are concisely described by three parameters: the square of the buoyancy frequency

_{g}*N*

^{2}≡ ∂

*b*/∂

*z*[where

*b*= −(

*g*/

*ρ*)

_{o}*ρ*is the buoyancy,

*g*is the gravitational acceleration,

*ρ*is a reference density, and

_{o}*ρ*is the density of the fluid]; a quantity proportional to the vertical component of the absolute vorticity of the geostrophic flow

*f*

^{ 2}≡

*f*(

*f*− ∂

*u*/∂

_{g}*y*), where

*f*is the Coriolis parameter; and a quantity describing the baroclinicity of the zone expressed through the thermal-wind balance

*υ*,

*w*) is the ageostrophic secondary circulation in the (

*y*−

*z*) plane, and

*X*and

*Q*are generic momentum and buoyancy source terms, whose specific form will be described below. A single equation for the streamfunction

*ψ*(where

*υ*≡ ∂

*ψ*/∂

*z*and

*w*≡ −∂

*ψ*/∂

*y*) describing the ASC can be obtained by using (1) to eliminate the rate of change and geostrophic advection terms of (2) and (3), yielding the time-independent form of the Eliassen–Sawyer equation (Eliassen 1951; Sawyer 1956),

*ω*equation (Hoskins 1982).

*ω**is the absolute vorticity vector in the (*

_{a}*y*–

*z*) plane (Hoskins 1982). The sign of the PV is determined by the relative slopes of surfaces of constant density to those of constant absolute momentum

*ψ*decay away from regions of forcing. Solutions for the ASC take on a cellular nature when isopycnals are steeper than

*τ*is the alongfront, zonal wind stress applied to the baroclinic zone) is modified by the surface vertical vorticity of the geostrophic flow

^{x}_{w}*ζ*= –∂

_{s}*u*/∂

_{g}*y*|

_{z}_{=0}. This result, which was first derived by Stern (1965) and Niiler (1969), is a consequence of nonlinear advection of momentum by the meridional Ekman flow

*υ*. In deriving (9) it has been assumed that the momentum source/sink

_{e}*X*represents a wind-induced frictional force confined to the Ekman layer of thickness

*δ*. It has also been assumed that the Ekman layer is much thinner than the vertical scale

_{e}*D*of the geostrophic flow; that is,

*h*are considered. With these vertical scales in mind, approximation (10) becomes

*δ*≪

_{e}*h*. Although it is conventionally assumed that the mixed layer and Ekman layer depths are equivalent, this need not be true for fronts with deep mixed layers, as is the case for fronts associated with the Antarctic Circumpolar Current, where the mixed layer depth can reach 500 m (Rintoul and England 2002). For mixed layers less than 100 m deep, direct observations of wind-derived ageostrophic transport indicate that wind-driven momentum can penetrate through the base of the mixed layer (Chereskin and Roemmich 1991; Wijffels et al. 1994; Lee and Eriksen 1996). Therefore, the theoretical analysis and use of approximation (10) may be more applicable to fronts with deep mixed layers. Having stated this caveat however, it should be pointed out that no approximation for the relative sizes of the mixed layer and Ekman layer depths is required for the numerical simulations of section 4 so that numerical experiments with Ekman depths of the same order as the mixed layer depth can and will be performed to see how critical approximation (10) is to the frontogenetic mechanism detailed in the theoretical analysis.

*ψ*,

*F*is the buoyancy flux associated with vertical mixing in the Ekman layer. In (12), the lateral buoyancy gradient has been evaluated at the surface owing to its negligible vertical variation across the Ekman layer in the limit where

^{B}_{e}*δ*≪

_{e}*D*. Integrating (12) down from the surface through the Ekman layer, one obtains an expression for the flux of buoyancy exiting the base of the Ekman layer

*F*

^{B}_{atm}is the buoyancy loss/gain to the atmosphere. The second term in (13), a result which has been independently derived in the context of thermocline ventilation (Marshall and Nurser 1992), wind effects on convection (Straneo et al. 2002), and nonlinear stratified spinup (Thomas and Rhines 2002), is a wind-driven buoyancy flux (WDBF) representing the flux of buoyancy exiting (entering) the fluid beneath the Ekman layer that is required to balance the tendency of lateral advection of buoyancy by the flow in the Ekman layer to reduce (increase) the buoyancy when the wind stress is down-front (up-front).

*O*(

*N*

^{−1}) and is characterized by overturning cells of small horizontal width of order the mixed layer depth. Therefore, the effects of gravitational instability on the larger-scale, subinertial ASC must be parameterized. Convection is parameterized by a buoyancy sink/source equal to the convergence of the convective buoyancy flux

*F*=

^{B}_{c}*w*and

_{c}*b*are the vertical velocity and buoyancy perturbations of the convection, and the bar represents an ensemble average), that is,

_{c}*K*profile parameterization” (KPP) of Large et al. (1994) incorporates this nonlocal behavior into their mixing scheme for unstable conditions, parameterizing the turbulent buoyancy flux in the convective limit by the following form:

*F*is the surface buoyancy flux and

^{B}_{s}*G*is a function describing the vertical structure of the flux in the mixed layer of thickness

*h*. In the limit that

*δ*≪

_{e}*h*, the surface buoyancy flux can be approximated by the effective buoyancy flux of (13); that is,

*F*→

^{B}_{s}*F*

^{B}_{eff}.

The novelty of the present study is the incorporation of nonlinear Ekman effects [see (9), (11), (12), and (13)] into the dynamics of the ASC. In doing this, processes that drive an ASC, that is, lateral gradients in vertical mixing of buoyancy and Ekman pumping, are no longer solely determined by the atmospheric forcing, but instead are functions of the lateral density contrast and vertical vorticity of the evolving front. In this way, Ekman pumping and the buoyancy source/sink term of (4) are analogous to the geostrophic forcing, and like the geostrophic forcing, it will be shown that they too drive a frontogenetic ASC.

## 3. Analytic theory

### a. Simple model of a baroclinic zone

*(*u

*z*) in a thermal wind balance with a buoyancy field

*(*b

*y*,

*z*), that is,

S

^{2}≡

*fd*u /

*dz*= −∂

*/∂*b

*y*(where

S

^{2}is a constant for simplicity; i.e.,

*has uniform vertical shear), embedded in a mixed layer of thickness*u

*h*, bounded below by a pycnocline (see Fig. 1). For the subinertial time scales appropriate to semigeostrophic dynamics, the vertical stratification in the mixed layer is nonzero because of restratification by geostrophic adjustment (Tandon and Garrett 1994). The degree to which small-scale turbulent processes mix zonal momentum determines the mixed layer vertical stratification

*N*

^{2}

_{ml}, yielding values ranging from

*N*

^{2}

_{ml}=

S

^{4}/

*f*

^{ 2}for complete mixing of momentum, to

*N*

^{2}

_{ml}= 0 if the geostrophic flow is unaffected. With this range of stratification, the potential vorticity (PV) in the mixed layer

*q*

_{ml}≡

*N*

^{2}

_{ml}

*f*

^{ 2}−

S

^{4}is less than or equal to zero, that is, −

S

^{4}≤

*q*

_{ml}≤ 0, owing to the geostrophic shear and lateral density gradient. Within the mixed layer

N

^{2}= ∂

*/∂*b

*z*is constant, while beneath the mixed layer, it is specified that

N

^{2}grows linearly with depth, increasing by an amount

S

^{4}/

*f*

^{ 2}over a depth

*δ*; that is,

*H*is the thickness of the mixed layer at

*t*= 0. This is meant to be a simple representation of a pycnocline whose stiffness is reflected in the smallness of

*δ*. For simplicity, the limit of a rigid pycnocline, that is,

*δ*/

*H*→ 0 will be considered. In this case,

*h*=

*H*for all times, and the vertical velocity must vanish at the base of the mixed layer. Solutions for flexible pycnoclines can be found in Thomas (2003) and differ only in the details from solutions for the ridgid pycnocline limit discussed below.

Added to this laterally homogeneous flow, which will be referred to as the “basic state,” is a meridionally varying geostrophic disturbance with zonal velocity *ũ*(*y*, *z*, *t*) and buoyancy *b̃*(*y*, *z*, *t*) meant to represent an incipient frontal jet. To make the analysis amenable to analytical solutions, the disturbance is considered weak relative to the basic state, i.e., *ũ* ≪ * u* and

*b̃*≪

*. The objective of the analysis is to determine if the frontal jet grows when exposed to a spatially uniform down-front wind stress*b

*τ*=

^{x}_{w}*τ*> 0. So as to isolate the effects of atmospheric forcing, no geostrophic confluence is imposed; that is,

_{o}*υ*= 0.

_{g}*t*= 0 the atmospheric forcing is initiated, driving Ekman pumping/suction (11) and an effective surface buoyancy flux (13), which for weak disturbances take the following linearized forms:

*M*=

^{o}_{e}*τ*/

_{o}*ρ*is the magnitude of the classical Ekman transport,

_{o}f*ζ̃*= −∂

_{s}*ũ*/∂

*y*|

_{z=0}, and

*S̃*

^{2}

_{s}=

*f*∂

*ũ*/∂

*z*|

_{z=0}= −∂

*b̃*/∂

*y*|

_{z=0}. The effective surface buoyancy flux (18) is mostly spatially uniform but is modulated by the “dynamic cooling” (19), that is, the part of the WDBF that evolves with the disturbance’s surface lateral buoyancy gradient and vertical vorticity, the later functionality arising from the dependence of the WDBF on the nonlinear Ekman transport (9), which decreases in magnitude in regions of positive vertical vorticity.

*F*

^{B}_{eff}is positive—that is, destabilizing—convection is induced. The convective buoyancy flux is parameterized using (15), yet a choice must be made on its vertical structure function

*G*. In the mixed layer,

*F*typically decreases linearly with depth, turning negative if the convection is penetrative as it entrains fluid from the pycnocline (Large et al. 1994). For simplicity, the convection is approximated to be nonpenetrative—that is, its buoyancy flux is zero at the base of the mixed layer:

^{B}_{c}### b. Governing equations for weak disturbances

*Q*=

*Q*—and that momentum sources/sinks are confined to the Ekman layer (hence their effects are captured in the Ekman dynamics), yielding

_{c}*w*is equal to the Ekman pumping/suction at the surface [see (17)] and

*w*=

*ψ*= 0 for

*z*≤ −

*H*(a consequence of the rigid pycnoline).

### c. Solution for the ageostrophic secondary circulation

*ζ̃*and the dynamic cooling

_{s}*F*

_{dyn}. To obtain a solution for the ASC, a meridional structure for

*ζ̃*and

_{s}*F*

_{dyn}must be chosen. The meridional structure of the disturbance is chosen to be sinusoidal with a wavelength

*L:*

*l*= 2

*π*/

*L*.

*ψ*and Ekman pumping

_{c}*ψ*

_{ep}, that is,

*ψ*= 0. The Ekman pumping driven ASC is a boundary driven homogeneous solution to (23) having the form

*γ*=

*f*

^{ 2}/

S

^{2}. The

*ansatz*(26) transforms the homogeneous form of (23) into an ordinary differential equation for the vertical structure function

*q*

_{ml}> 0, (27) is elliptic, yielding solutions which decay away from the surface. In this case, the depth over which the ASC penetrates is the semigeostrophic Prandtl depth

*δ*

_{SG}=

*Lf*/(2

*π*

*N*

^{2}

_{ml}−

S

^{4}/

*f*

^{ 2}

*q*

_{ml}< 0, (27) is hyperbolic and solutions are cellular so that the ASC does not necessarily decay away from the surface. Mathematically summarizing the Ekman pumping part of the solution:

*m*=

*πL*/(2

_{o}*HL*), where

*ψ*solves (23) subject to an upper boundary condition of

_{c}*w*= 0 and is as follows:

*satisfies the following relationship*

_{p}*ψ*is an ASC that balances convective mixing through vertical advection of buoyancy.

_{c}### d. Evolution equations and stability analysis

Although the spatial structure of the secondary circulation is known, its temporal evolution is as yet unknown. The solution for *ψ* depends only parametrically on time through the surface vorticity and dynamic cooling, therefore it is the evolution equations of *ζ̃ _{s}* and

*F*

_{dyn}that determine the time-dependence of

*ψ*and the disturbance.

*ζ̃*and

_{s}*S̃*

^{2}

_{s}, the latter of which is governed by the following equation:

*M*S

^{o}_{e}^{2}/

*f*× (32) from

*M*× (33) yields an equation for the evolution of

^{o}_{e}*F*

_{dyn}

*F*

_{dyn}.

*ζ̂*and

_{s}*F̂*:

*ζ̂*,

_{s}*F̂*) = (

*ζ̂*,

^{o}_{s}*F̂*)

^{o}*e*

^{σt}into (35) generates an eigenvalue problem for the growth rate

*σ*of the disturbance. In this way, a stability analysis can be performed that determines whether or not a geostrophic disturbance of a particular wavelength introduced to the baroclinic zone grows or decays.

The dependence on *L* of the growth rate of the fastest growing eigenvector for a basic state with the following parameters: ^{2} = 8.8 × 10^{−7} s^{−2}, *f* = 9.3 × 10^{−5} s^{−1}, *H* = 50 m, *q*_{ml} = ±0.1^{4}, and *τ _{o}* = 0.1 N m

^{−2}was calculated and is plotted in Fig. 2. The growth rate is complex,

*σ*=

*σ*+

_{r}*iσ*, and varies greatly with the wavelength of the disturbance and the PV of the mixed layer. For

_{i}*q*

_{ml}< 0, disturbances with wavelengths

*L*≃

*L*/(2

_{o}*n*+ 1), with

*n*being an integer greater than or equal to zero, are the fastest growing and are quasi-stationary, that is, |

*σ*| ≫ |

_{r}*σ*|, while propagating disturbances travel with or against the Ekman transport and grow more slowly. Rapidly propagating disturbances have wavelengths

_{i}*L*≃

*L*/(2

_{o}*n*) and small

*σ*. Unlike the

_{r}*q*

_{ml}< 0 case, all disturbances in positive PV mixed layers propagate in the direction of the Ekman transport (southward), that is,

*σ*> 0, and the fastest growing disturbances are those with the shortest wavelengths—that is, there is no scale selection.

_{i}*q*

_{ml}/

S

^{4}is fixed, the shape of the growth rate spectrum is unaffected by changes in

S

^{2},

*H*, and

*τ*; however, its scale is. Increasing the wind stress

_{o}*τ*thinning the mixed layer thickness

_{o,}*H*, or reducing the background lateral buoyancy gradient

S

^{2}all amplify the growth rate by an amount consistent with an advective time scale scaling (that is,

*HL*/

*M*) for the evolution of the disturbance: that is,

^{o}_{e}### e. Implications for frontogenesis

The objective of this study was to determine if destabilization of the water column by down-front winds could drive frontogenetic ASCs. It has been shown in the previous section that geostrophic frontal jets introduced into a baroclinic zone forced by down-front winds can propagate and grow. In this section, the physical mechanisms responsible for the propagation and growth of the frontal jets and the implication of these results for frontogenesis are summarized.

Propagation of the frontal jets is caused by vortex stretching associated with the Ekman pumping driven ASC *ψ*_{ep}. Vortex stretching forced by Ekman pumping (17) is in quadrature with the surface vertical vorticity and hence results in the propagation of *ζ̃ _{s}*. As demonstrated in Fig. 3 the vertical structure of

*ψ*

_{ep}, which is set by the PV of the mixed layer and the wavelength of the disturbance, determines the propagation speed and direction. For mixed layers with positive PV,

*ψ*

_{ep}decays in the vertical over a Prandtl depth, inducing vortex stretching that lags

*ζ̃*by 90° in

_{s}*y*and forces southward propagation, in the direction of the Ekman transport. Reducing

*L*for

*q*

_{ml}> 0 decreases the Prandtl depth, makes vortex stretching stronger, and hence results in faster propagation (i.e., larger

*σ*, see Fig. 2). Secondary circulations in mixed layers with negative PV form cellular structures and can decay, remain unchanged, or increase with depth depending on

_{i}*L.*For

*q*

_{ml}< 0 and

*L*=

*L*

_{0}[more generally,

*L*=

*L*

_{0}/(2

*n*+ 1)], surface flow is parallel to

*q*

_{ml}< 0 and

*L*=

*L*

_{0}/1.75,

*ψ*

_{ep}increases with depth making vortex stretching lead

*ζ̃*, driving northward propagation (i.e.,

_{s}*σ*< 0), counter to the Ekman transport.

_{i}The components of the analytical solution that are key to the growth of the frontal jets are shown in Fig. 4 using the solution for a disturbance with *L* = *L _{o}* and initial surface vorticity of magnitude 0.0025

*f*introduced into a mixed layer with negative PV. Seven inertial periods after the introduction of the disturbance, frontal interfaces have developed in the buoyancy field as a result of the convergent flow field of the ASC. Convective mixing is most intense beneath the fronts and drives an ASC

*ψ*characterized by downwelling on the dense side of the front, upwelling along the frontal interface, and northward surface flow coincident with the eastward flow of the frontal jets. The vertical vorticity of the frontal jets modifies the Ekman transport and results in Ekman pumping on the dense side of the front and suction/upwelling along the frontal interface. The spatial structure of the zonal velocity reflects the vertical circulation, with high momentum fluid subducted on the dense side of the front and low momentum fluid upwelled along the frontal interface, resulting in a reduction of the vertical vorticity at the front.

_{c}A schematic diagram illustrating the role of each key component of the analytical solution in the steps that lead to intensification of ocean fronts by down-front winds is shown in Fig. 5. These steps to frontogenesis are as follows. First, Ekman flux induces a wind-driven buoyancy flux and convective mixing of buoyancy that is concentrated at the front where the lateral buoyancy gradient is largest. Second, localized mixing drives an ASC with northward surface flow at the frontal outcrop, northward flow that accelerates down-wind frontal jets via the Coriolis force. Third, spinup of the frontal jets and associated surface vertical vorticity strengthens the Ekman pumping/suction, forcing water to downwell to the north of the fronts and upwell along the frontal interface. Last, this differential vertical motion tilts the meridional component of the vorticity downward, reducing the vertical vorticity, and consequently enhancing the Ekman transport at the front. The differential vertical motion also reduces the lateral density gradient; however, for mixed layers with negative potential vorticity, the enhancement of the Ekman transport outweighs the slight reduction in the lateral density gradient so that the wind-driven buoyancy flux and convective mixing experiences a net increase at the front. As a result, the mixing driven ASC is strengthened, the spinup of the frontal jets is accelerated, the Ekman pumping/suction is magnified, and the frontogenetic convergent flow of the total ASC is intensified. Repetition of these steps leads to the exponential growth of the disturbance and frontogenesis.

## 4. Nonhydrostatic numerical simulations

To test and extend the results of the theory outlined above, two-dimensional, nonhydrostatic, high-resolution numerical simulations run with the same basic state configuration of section 3a and forced by a spatially uniform down-front wind stress and atmospheric cooling were performed. A series of five experiments designed to cover a range of values of *L _{o}*, through variation of

S

^{2},

*H*,

*f*, and

*N*

^{2}

_{ml}were performed (see Table 1). All experiments were forced by the same wind stress

*τ*= 0.1 N m

_{o}^{−2}and atmospheric buoyancy flux

*F*= 6.3 × 10

_{o}^{−8}m

^{2}s

^{−3}. The exact configuration of the numerical model is detailed in the appendix.

Listed in Table 1 is the thickness of the Ekman layer *δ _{e}* for each of the experiments. Notice that for none of the experiments is

*δ*very much smaller than the initial depth of the mixed layer and hence the experiments do not strictly satisfy approximation (10) used in the theory. Therefore, a comparison of the numerical and analytical solutions can be used to determine if (10) is critical to the frontogenesis mechanism outlined in the theory.

_{e}A Hovmöller plot of −*ψ* evaluated at the base of the Ekman layer for experiment E is shown in Fig. 6 and illustrates the typical *y* − *t* structure of the near-surface ASC of the numerical simulations. Small-scale convective overturning cells of width the depth of the mixed layer *H* fill the domain after a half an inertial period. By two inertial periods, out of this chaos emerge three large-scale secondary circulations characterized by intense downwelling centered at fronts that do not travel with the Ekman flow, but remain stationary.

Meridional sections of the density and streamfunction for experiment E reveal that frontal interfaces form, convection is concentrated beneath and to the north of the fronts, and Ekman transport (southward flow in the Ekman layer, i.e., upper 10 m) is channeled down the frontal interface (Fig. 7a). The zonal velocity of the numerical solution, like the analytical solution (cf. Figs. 4b and 7b) reflects the vertical circulation, with high momentum subducted down the dense side of the front and low momentum upwelled at the base of the frontal interface. Down-wind surface jets develop just to the south of the front in a region where the Ekman flow stagnates or reverses.

The lateral structure of the surface expression of the fronts is illustrated in Fig. 8. So as to highlight the robust features of the fronts of experiment E, a composite front has been constructed and is used in the figure. The method of constructing this composite front is as follows: all fronts in a particular meridional section are aligned relative to their maximum lateral density gradient and the particular variable of interest is then averaged both spatially over the number of fronts in the section and temporally over the last two inertial periods of the experiment. The intensity of this composite front is evidenced by the precipitous drop in buoyancy and zonal velocity (which has been averaged in the vertical over the Ekman layer) crossing the front from south to north (Figs. 8a and 8b). The tremendous lateral shear in the zonal velocity gives rise to positive vertical vorticity at the front that greatly exceeds *f.* With such large vorticity, nonlinear Ekman theory [see (9)] predicts that the front should effectively form a barrier to the Ekman transport. In Fig. 8c, the Ekman transport derived from (9) is compared with −*ψ* evaluated at the base of the Ekman layer. The Ekman transport estimated using the vertical vorticity calculated from a smoothed version of the zonal velocity agrees well with −*ψ* north of the front. On account of the singularity of (9) at *ζ _{s}* = −

*f*, the Ekman transport calculated from the zonal velocity without smoothing is scattered, yet drops to zero at the center of the front, similar to −

*ψ*. This sharp shutdown of the Ekman transport generates powerful Ekman pumping and leads to the downturn of

*ψ*at the fronts shown in Fig. 7a.

South of the front the Ekman transport is lower than −*ψ* by 0.5–1.0 m^{2} s^{−1}. To account for this discrepancy, there must be a northward transport added to the Ekman flow. This presumed northward flow is coincident with the down-wind surface jets. The Coriolis force associated with the northward flow is in the right sense to accelerate the frontal jet. Hence, in analogy with the analytical theory of section 3e, this flow may be attributable to the convective mixing driven ASC *ψ _{c}* that accomplishes this task in the theory. If this is the case, the downwelling of the northward flow at the front (where it must vanish as the discrepancy between Ekman transport and −

*ψ*is small) should be a location of enhanced convective mixing. To test this, the convective buoyancy flux

*F*=

^{B}_{c}To calculate the convective buoyancy flux, *w* and *b* were high-pass filtered in *y* to separate *w _{c}* and

*b*from the larger-scale flows associated with the ASC. The product of

_{c}*w*and

_{c}*b*was then averaged using the compositing method outlined above to obtain

_{c}*F*. As illustrated in Fig. 9a the wind-driven buoyancy flux and

^{B}_{c}*F*evaluated at

^{B}_{c}*z*= −7.8 m, near the base of the Ekman layer, scale well with each other within ∼700 m of the front. Both buoyancy fluxes are concentrated sharply at the front where they far exceed the strength of the atmospheric buoyancy flux

*F*. An example of overturning cells associated with gravitational instability is shown in Fig. 9b. Convection is characterized by a strong overturning cell locked to the front that vigorously mixes down dense water of northern origin advected toward the front by the Ekman flow. A vertical profile of

_{o}*F*associated with this convection cell locked to the front (see Fig. 9c) matches the wind-driven buoyancy flux at the base of the Ekman layer and decays linearly with depth through the mixed layer, as was assumed in the analytical theory [see (15)]. However, over the depth range

^{B}_{c}*z*< −25 m and

*z*> −35 m

*F*is negative, indicating that the convection is penetrative and entrains fluid from the pycnocline. A train of weaker convection cells extends ∼700 m north of the front, yielding a relatively homogeneous meridional distribution of

^{B}_{c}*F*. Beyond ∼700 m north of the front, even though convection cells are not evident in the plot of the high-passed streamfunction shown in Fig. 9b, convection is occurring as

^{B}_{c}*F*≠ 0. In this region,

^{B}_{c}*F*decays nearly exponentially in

^{B}_{c}*y*, resulting in vertical mixing of buoyancy by convection too weak to homogenize the density inversions that are being generated by the Ekman flow. A local linear stability analysis of the flow in this region (which has been omitted for brevity) reveals that gravitational instabilities both grow and propagate to the south, with the Ekman flow. The combination of growth and southward propagation of convection leads to the exponential increase of

*F*with decreasing

^{B}_{c}*y*. The distance it takes the instability to reach finite amplitude depends on the strength of noise in the buoyancy and flow fields which sets the initial amplitude of the disturbance: the weaker the noise the longer the distance. When convection reaches finite amplitude,

*F*scales with the wind-driven buoyancy flux and density inversions are rapidly mixed away. Therefore, the finite distance over which density inversions persist in the experiment, evident in Figs. 7 and 9b, is a consequence of the weak noise level in the numerical model.

^{B}_{c}A key quantity predicted by the theory is the wavelength of the fastest growing disturbance. As shown in Fig. 6, there is a clear wavenumber selection for the ASC in experiment E. This was true for all the experiments listed in Table 1. To assess whether the distance of separation between fronts and wavelength of the ASC scales with the predicted length scale *L _{o}*, spectra in

*y*of the streamfunction at the base of the Ekman layer were calculated, and are shown in Fig. 10. For all the experiments, the spacing of the fronts is proportional to

*L*, with a multiplier of 2/3. It was found that for all the experiments, the mixed layer depth shoaled from its initial value. An example of this can be seen in the depth of penetration of the convective buoyancy flux (see Fig. 9c), which is 30 m versus the initial mixed layer depth of 50 m. Since

_{o}*L*is proportional to the mixed layer depth, it is conceivable that the reduced wavelengths of the ASCs in the numerical experiments are associated with the thinner mixed layers. The shoaling of the mixed layer was due to the large vertical diffusivity of buoyancy, which resulted in an upward diffusion of the pycnocline.

_{o}## 5. Implications for the subpolar front of the Japan/East Sea during cold-air outbreaks

^{−2}) and powerful, northwesterly wind stress (up to 0.5 N m

^{−2}) with a significant down-front component, making the front a prime testing ground for the theory. During January and February 2000, as part of the Japan Sea Initiative of the Office of Naval Research, a series of high-resolution SeaSoar (towed profiler instrumented to measure temperature, salinity, depth, and chlorophyll fluorescence) and ship-board acoustic Doppler current profiler (ADCP) surveys sampled the oceanic response to three cold-air outbreaks. During the outbreaks, the two-dimensional PV of the flow

*u*, was negative at the front owing to the strong vertical shear and horizontal density gradient (see Fig. 11a). Beneath the front pycnostads were found with tracers of northern, surface waters: low PV, freshwater, and high chlorophyll fluorescence (not shown), a sign of active subduction. The alongfront velocity (Fig. 11b) was characterized by a high Rossby number (|

*ζ*|/

_{s}*f*∼ 0.7), downwind jet at the surface, and an upwind flow at the base of the frontal interface, a flow pattern reminiscent of the upwelled low-momentum fluid of the analytical and numerical solutions (see Figs. 4b and 7b). The Ekman transport (9) estimated using the surface vertical vorticity and the zonal component of the wind stress derived from shipboard meteorological measurements differs dramatically from its linear counterpart, with strong divergence at the center of the front midway between the cyclonic and anticyclonic vorticity of the jet (Fig. 11d). Utilizing the Ekman transport to calculate the wind-driven buoyancy flux of (13), shows that the destabilizing buoyancy flux is concentrated at the front, reaching a peak value of 2.5 × 10

^{−6}m

^{2}s

^{−3}(Fig. 11c). In units of heat flux, this wind-driven buoyancy flux would be equivalent to a heat loss of 4900 W m

^{−2}, nearly an

*order of magnitude larger*than the

*atmospheric heat loss.*The negative PV, evidence for subduction, subsurface upwind flow, Ekman suction up the frontal interface, and strong, concentrated wind-driven buoyancy flux all suggest that wind-driven secondary circulation and frontal intensification are active at the front during cold-air outbreaks.

A comparison between the observed width of the subpolar front of the Japan/East Sea to *L _{o}*, the length scale of the fastest growing frontogenetic ASC, can be used to qualitatively evaluate the applicability of the theory to the observations. Substituting an average value for the negative PV in the core of the front based on the observations

*q*

_{ml}∼ −2 × 10

^{−13}s

^{−4}, a mixed layer depth

*H*∼ 50 m, and the Coriolis parameter

*f*= 9.35 × 10

^{−5}s

^{−1}into (29) yields

*L*≈ 10 km, a value comparable to the observed width of the front.

_{o}Because the Japan/East Sea subpolar front exhibits energetic meanders and eddies, two-dimensional theories cannot fully characterize its dynamics. Nonetheless, the theoretical model developed in this study can provide an upper bound on the strength of vertical circulation at the front. For wind-driven frontogenesis, the frontal vertical velocity scales with the Ekman pumping and downwelling balancing convective mixing. As shown in Fig. 11d, the Ekman transport changes by about 5 m^{2} s^{−1} in 5 km so that the Ekman pumping/suction would be 0.1 cm s^{−1}. Using the scaling for *ψ _{c}* (31) to estimate the strength of the downwelling balancing convective mixing, that is,

*w*≈

*F*

^{B}_{eff}/(

*N*

^{2}

*H*), with

*F*

^{B}_{eff}≈ 2 × 10

^{−6}m

^{2}s

^{−3},

*N*

^{2}∼ 2 × 10

^{−5}s

^{−2}(based on an average value for the vertical stratification at the core of the front), and

*H*∼ 50 m, yields a value of 0.2 cm s

^{−1}. Such large vertical velocities of 0.1–0.2 cm s

^{−1}(86–170 m day

^{−1}) will not necessarily cause the mixed layer to deepen by ∼100 m in a day because vertical motions are not purely perpendicular to isopycnals on account of the slanted nature of the frontal interface. If the ASC at the subpolar front were similar to the analytical and numerical solutions plotted in Figs. 4 and 7, isopycnals would be distorted by the ASC in the following manner. Upwelling centered at the front by the Ekman divergence lifts the frontal interface, while beneath the Ekman layer convergence intensifies the lateral buoyancy gradient of the front. On the dense side of the front, water downwells on a slanted path, tucking surface fluid from the north under the frontal interface. This “lift-and-tuck” flow configuration of the ASC is conducive to subduction and might explain the existence of the boluses of low PV freshwater found underneath the front.

## 6. Conclusions

A spatially uniform wind blowing over a baroclinic zone in the direction of the surface currents leads to the formation of strong frontogenetic ageostrophic secondary circulations and multiple fronts within the zone.

Using an analytic theory, it has been shown that for mixed layers with negative potential vorticity, a particular length scale *L _{o}* [see (29)] is selected by the most frontogenetic, stationary ASCs and determines the separation distance between fronts. Frontogenesis is ultimately a consequence of wind-driven gravitational instability and nonlinear Ekman pumping. Ekman flux of down-front winds advects dense water over light triggering convection centralized to the front. Mixing of buoyancy by this convection drives an ASC that accelerates the frontal jet. The vorticity contrast of the jet induces Ekman pumping/suction that enhances the Ekman flux at the front, hence strengthening the destabilizing density advection, subsequent convective mixing, and jet-accelerating ASC. Repetition of this process leads to frontal intensification within several inertial periods, with stronger winds producing faster frontogenesis. This frontogenesis mechanism does not require that the wind stress have negative curl (or any curl for that matter) as is the case with classic Ekman frontogenesis (e.g., Roden 1980; de Szoeke 1980; Cushman-Roisin 1981). Nor does it require a horizontal deformation field and thus is distinct from frontogenesis mechanisms based on the frontal model of Hoskins and Bretherton (1972).

The vertical circulation associated with the ASCs is characterized by subduction on the dense side of the front and upwelling along the frontal interface. Vertical velocity is scaled by both nonlinear Ekman pumping and downwelling balancing wind-driven convective mixing. The magnitude of the latter is set by a wind-driven buoyancy flux proportional to the product of the Ekman transport with the lateral buoyancy gradient at the front, a quantity that can be much larger than the atmospheric buoyancy loss owing to the large frontal density contrast.

Nonhydrostatic numerical simulations of a baroclinic zone forced by down-front winds capture the process of frontal intensification by the formation of frontogenetic ASCs in a manner consistent with the mechanism presented in the analytic theory. The experiments verify several key predictions of the theory: 1) fronts form that do not propagate with the Ekman flow, but remain stationary; 2) the stationary fronts are separated by a distance that scales with *L _{o}*; 3) convection is concentrated in the proximity of the front and leads to a convective buoyancy flux that scales with the wind-driven buoyancy flux; and 4) frontal jets develop with strong positive vertical vorticity that effectively form a barrier to the Ekman transport, causing the Ekman transport to be channeled down the dense side of the front. This agreement between theory and numerical experiments, experiments with Ekman and mixed layer depths of comparable size, indicates that assumption (10) used in the theory that the Ekman layer is much thinner than the mixed layer is not crucial to frontal intensification by down-front winds.

High-resolution cross-front sections of density, horizontal velocity, and PV made at the subpolar front of the Japan/East Sea during strong atmospheric forcing by northwesterly winds with a significant down-front component associated with cold-air outbreaks suggest that wind-driven frontogenesis is active at the front. Evidence to support this hypothesis was: the occurrence of negative PV at the front; the discovery of waters beneath the frontal interface of northern origin apparently having been recently subducted; the estimation of a wind-driven buoyancy flux focused at the front peaking at a value equivalent to a heat loss of 4900 W m^{−2}, an order of magnitude larger than the atmospheric heat flux; and the calculation of significant cross-front variation of the Ekman transport owing to the vertical vorticity contrast of the frontal jet. The width of the front was found to scale with *L _{o}*. Upper bound estimates of the vertical velocity based on the theory (86–170 m day

^{−1}) were quite large in magnitude, indicating that neglecting wind-forcing and using the quasigeostrophic

*ω*equation to infer vertical velocity at wind-forced fronts could lead to significant errors. However, assuming that the front is purely two-dimensional by neglecting the effects of geostrophic forcing resulting from meanders and eddies could also lead to errors in the estimation of the vertical velocity. The Eliassen–Sawyer equation (4) accounts for both geostrophic and wind forcing and thus provides an appropriate theoretical model for describing vertical frontal circulations at wind-driven fronts with low to negative PV. In a future study, the complete set of data taken at the subpolar front of the Japan/East Sea, including both cross and alongfront measurements, will be used to obtain solutions of the Eliassen–Sawyer equation. These solutions will be utilized to both estimate the ASC at the front and to determine the relative importance of geostrophic versus wind forcing in the driving of the ASC.

The dominance of the destabilizing wind-driven buoyancy flux over the atmospheric buoyancy flux in the analytical and numerical solutions as well as in the estimate from the subpolar front of the Japan/East Sea suggests that Ekman advection of buoyancy could play a major role in the formation of mode waters at fronts forced by down-front winds. This is especially true for fronts in the Southern Ocean where there are no cold-air outbreaks on account of the lack of continents. Evidence pointing to the importance of Ekman dynamics in the formation of mode waters in the Southern Ocean is given by Rintoul and England (2002) who show that the magnitude of temperature and salinity variation in the Subantarctic Mode Water cannot be explained by air–sea fluxes of heat and freshwater, but can be accounted for by advection of cold, freshwater across the subantarctic front by northward Ekman transport. As shown in the example of the subpolar front of the Japan/East Sea forced by cold-air outbreaks, it is not necessary that the atmospheric buoyancy flux be weak for the wind-driven buoyancy flux to dominate. Like the subpolar front of the Japan/East Sea, the Kuroshio is forced by cold dry air during the winter. The formation of mode waters associated with the Kuroshio has been attributed to wintertime convection driven by heat loss. However, Yasuda and Hanawa (1997) conclude that decadal variations in the temperature of the North Pacific central mode water are due to both decadal changes in the air–sea heat flux and Ekman advection of heat, with the latter contributing twice as much to the decadal variability. Both of these examples suggest that wind-driven frontogenesis may be active at these large-scale fronts of global importance. The connection with mode water formation further implies that knowledge of the processes involved in wind-driven fronts is important for understanding the ocean–atmosphere climate system. In this respect, care should be taken in designing global circulation models so that in frontal regions the wind-driven buoyancy flux, ensuing gravitational instability, and frontogenetic ASCs are well parameterized.

## Acknowledgments

We thank Peter Rhines, LuAnne Thompson, Parker MacCready, and Jonathan Lilly, who improved this manuscript with helpful comments and fruitful discussions. We also thank the members of the Woods Hole SeaSoar Group—Frank Bahr, Jerry Dean, Paul Fucile, Allan Gordon, and Craig Marquette, whose efforts made the Japan/East Sea observations possible. Captain Christopher Curl and the crew of the R/V *Revelle* provided professional, enthusiastic assistance through difficult wintertime operating conditions. They were a pleasure to work with, and we are grateful for both their exceptional skill and good nature. L. T. was supported under NSF grants OCE-00-95971 and OCE-03-51191. The Japan/East Sea program was supported by the Office of Naval Research under Grant N00014-98-1-0370 (CML).

## REFERENCES

Bleck, R., R. Onken, and J. D. Woods, 1988: A two-dimensional model of mesoscale frontogenesis in the ocean.

,*Quart. J. Roy. Meteor. Soc.***114****,**347–371.Chereskin, T. K., and D. Roemmich, 1991: A comparison of measured and wind-derived Ekman transport at 11°N in the Atlantic Ocean.

,*J. Phys. Oceanogr.***21****,**869–878.Cushman-Roisin, B., 1981: Effects of horizontal advection on upper ocean mixing: A case of frontogenesis.

,*J. Phys. Oceanogr.***11****,**1345–1356.da Silva, A. M., C. C. Young, and S. Levitus, 1994:

*Algorithms and Procedures*. Vol. 1,*Atlas of Surface Marine Data 1994,*NOAA Atlas NESDIS 6, 83 pp.de Szoeke, R. A., 1980: On the effects of horizontal variability of wind stress on the dynamics of the ocean mixed layer.

,*J. Phys. Oceanogr.***10****,**1439–1454.Eliassen, A., 1951: Slow thermally or frictionally controlled meridional circulation in a circular vortex.

,*Astrophys. Norvegia***5****,**19–60.Hoskins, B. J., 1982: The mathematical theory of frontogenesis.

,*Annu. Rev. Fluid Mech.***14****,**131–151.Hoskins, B. J., and F. P. Bretherton, 1972: Atmospheric frontogenesis models: Mathematical formulation and solution.

,*J. Atmos. Sci.***29****,**11–37.Hoskins, B. J., and I. Draghici, 1977: The forcing of ageostrophic motion according to the semi-geostrophic equations and in an isentropic coordinate model.

,*J. Atmos. Sci.***34****,**1859–1867.Hoskins, B. J., I. Draghici, and H. C. Davies, 1978: A new look at the

*ω*-equation.,*Quart. J. Roy. Meteor. Soc.***104****,**31–38.Large, W. G., J. C. McWilliams, and S. C. Doney, 1994: Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parameterization.

,*Rev. Geophys.***32****,**363–403.Latif, M., and T. P. Barnett, 1994: Causes of decadal climate variability over the North Pacific and North America.

,*Science***266****,**634–637.Lee, C. M., and C. C. Eriksen, 1996: The subinertial momentum balance of the North Atlantic subtropical convergence zone.

,*J. Phys. Oceanogr.***26****,**1690–1704.MacVean, M. K., and J. D. Woods, 1980: Redistribution of scalars during upper ocean frontogenesis: A numerical model.

,*Quart. J. Roy. Meteor. Soc.***106****,**293–311.Marshall, J. C., and A. J. G. Nurser, 1992: Fluid dynamics of oceanic thermocline ventilation.

,*J. Phys. Oceanogr.***22****,**583–595.Niiler, P., 1969: On the Ekman divergence in an oceanic jet.

,*J. Geophys. Res.***74****,**7048–7052.Oort, A. H., L. A. Anderson, and J. P. Peixoto, 1994: Estimates of the energy cycle of the oceans.

,*J. Geophys. Res.***99****,**7665–7688.Rintoul, S. R., and M. H. England, 2002: Ekman transport dominates local air–sea fluxes in driving variability of Subantarctic Mode Water.

,*J. Phys. Oceanogr.***32****,**1308–1321.Roden, G. I., 1980: On the subtropical frontal zone north of Hawaii during winter.

,*J. Phys. Oceanogr.***10****,**342–362.Sawyer, J. S., 1956: The vertical circulation at meterological fronts and its relation to frontogenesis.

,*Proc. Roy. Soc. London***A234****,**346–362.Spall, M. A., 1995: Frontogenesis, subduction, and cross-front exchange at upper ocean fronts.

,*J. Geophys. Res.***100****,**2543–2557.Spall, M. A., 1997: Baroclinic jets in confluent flow.

,*J. Phys. Oceanogr.***27****,**1054–1071.Stern, M. E., 1965: Interaction of a uniform wind stress with a geostrophic vortex.

,*Deep-Sea Res.***12****,**355–367.Straneo, F., M. Kawase, and R. S. Pickart, 2002: Effects of wind on convection in strongly and weakly baroclinic flows with application to the Labrador Sea.

,*J. Phys. Oceanogr.***32****,**2603–2618.Tandon, A., and C. Garrett, 1994: Mixed layer restratification due to a horizontal density gradient.

,*J. Phys. Oceanogr.***24****,**1419–1424.Thomas, L. N., 2003: Nonlinear dynamics of wind-driven oceanic circulation and fronts. Ph.D. thesis, University of Washington, 134 pp.

Thomas, L. N., and P. B. Rhines, 2002: Nonlinear stratified spin-up.

,*J. Fluid Mech.***473****,**211–244.Wang, D. P., 1993: Model of frontogenesis: Subduction and upwelling.

,*J. Mar. Res.***51****,**497–513.Wijffels, S., E. Firing, and H. Bryden, 1994: Direct observation of the Ekman balance at 10°N in the Pacific.

,*J. Phys. Oceanogr.***24****,**1666–1679.Yasuda, T., and K. Hanawa, 1997: Decadal changes in the mode waters in the midlatitude North Pacific.

,*J. Phys. Oceanogr.***27****,**858–870.

## APPENDIX

### Numerical Model Equations and Configuration

*ψ*is a streamfunction; that is,

*υ*≡ ∂

*ψ*/∂

*z*and

*w*≡ −∂

*ψ*/∂

*y*,

*ξ*≡ −∇

^{2}

*ψ*is the zonal component of the vorticity,

*J*(

*A*,

*B*) ≡

*A*−

_{y}B_{z}*A*is the Jacobian,

_{z}B_{y}*κ*= 0.005 m

_{υ}^{2}s

^{−1}(

*κ*= 0.1 m

_{H}^{2}s

^{−1}) is the vertical (horizontal) diffusivity of momentum and buoyancy, and the total zonal velocity and buoyancy fields have been decomposed as

*u*= (

S

^{2}/

*f*)

*z*+

*u*′(

*y*,

*z*,

*t*) and

*b*= −

S

^{2}

*y*+

*b*′(

*y*,

*z*,

*t*), where

S

^{2}is a constant. The vertical stratification of the basic state (16) is incorporated into the initial condition for

*b*′

*u*′,

*ξ,*and

*ψ*—are set to zero at

*t*= 0.

*u*′,

*ξ,*and

*ψ*at the corners of the grid. Vertical and horizontal grid spacings of the model (Δ

*z*= 2 m and Δ

*y*= 10 m) are designed to be fine enough so that the Ekman layer (which unlike the idealized Ekman layer used in the theory is of finite thickness) and narrow convective plumes are well resolved. An Adams–Bashforth time-stepping scheme is used with a 6-s time step. Lateral boundary conditions at

*y*= 0 and

*y*= 12 km are periodic. On the top and bottom boundaries of the domain a no normal flow condition is used:

*ψ*= 0 at

*z*= 0, −

*D*, where the total depth of the fluid

*D*= 250 m is much deeper than the Prandtl depth calculated from the anticipated length scale

*L*of the ASCs. For the buoyancy and zonal velocity, flux and stress boundary conditions are specified

_{o}*τ*is the wind stress and

_{o}*F*is the atmospheric buoyancy flux. At

_{o}*t*= 0, the wind stress and buoyancy flux is turned on impulsively and the buoyancy in the top grid is seeded with a weak perturbation of random noise to trigger gravitational instability.

Experimental parameters for numerical simulations.