## 1. Introduction

The conceptual paradigm for planetary boundary layers is horizontal homogeneity, and often stationarity as well. Besides embodying simplicity, these assumptions express a space- and time-scale separation between the boundary layer turbulence and the interior flow and the boundary characteristics. Yet these assumptions are often false in nature. In particular, the boundary itself need not be homogeneous, which is the focus of this paper.

The most important new dynamical element with inhomogeneity is three-dimensional secondary circulations (SCs), also called standing eddies, that can extend throughout the boundary layer from the surface to the inversion layer adjacent to the vertical interior, and thereby also modulate the more familiar transient turbulent eddies. This behavior is relevant for both atmospheric and oceanic boundary layers. We are particularly concerned with inhomogeneities on intermediate scales, comparable to or somewhat larger than the boundary layer thickness. Random, smaller-scale surface irregularities can be considered as analogous to surface roughness in the momentum stress formula.

Previously we have examined the effect of sea surface temperature (SST) inhomogeneity in large-eddy simulations (LES) of a weakly convective atmospheric boundary layer (ABL) for configurations of a SST front and filament and for different geostrophic wind directions (Sullivan et al. 2020, 2021, hereafter S20 and S21, respectively). This is a thermodynamic ocean–atmosphere coupling effect through the surface heat flux. Here we examine the complementary effect of an inhomogeneous oceanic current acting through the surface stress, and the current has the idealized shape of an isolated circular eddy.

There is an extensive literature for the ABL above inhomogeneous surfaces, especially land-forms (e.g., Garratt 1992; Raasch and Harbusch 2001; Owinoh et al. 2005; Patton et al. 2005). Over the ocean SST gradients modulate the gravitational stability of the ABL and induce local changes in the surface wind, in particular at the oceanic mesoscale (∼100 km horizontally) (Chelton et al. 2011); this in turn forces inhomogeneous Ekman pumping in the ocean and impacts plankton distributions Gaube et al. (2013, 2015). Surface current gradients modulate the surface stress (Seo et al. 2016; Renault et al. 2016). Its clearest effect is to extract kinetic energy from the mesoscale currents (i.e., “eddy killing”), but its influence has also been shown for larger-scale currents (Renault et al. 2019). The feedback on the ABL from these surface couplings is less well understood, but it seems clear that they can be significant in low and high winds (Edson et al. 2007; Chen et al. 2007; Minobe et al. 2008; Frenger et al. 2013); this is an active area of research (Robinson et al. 2019; Ayet et al. 2021).

Much of the previous coupled air–sea modeling has relied on single-column boundary layer parameterizations, rather than turbulence-resolving LES, as here. Furthermore, most of the documented interactions are at the oceanic mesoscale (∼10–100 km). However, the oceanic surface is full of submesoscale currents and SST fronts and filaments (∼0.1–10 km) McWilliams (2016, 2020), and their coupling effects have only just begun to be explored (Renault et al. 2018); it is likely that significant impacts will continue to this smaller scale.

A framework for interpreting the consequences of stress inhomogeneity is uniform-density, rotating, laminar boundary layers. For swirling flows or a rotating-boundary, radial outflows or inflows occur near the surface that are closed by secondary circulations (e.g., Schlichting 1968; Rotunno 2013, 2014). In our LES the eddy current differs from solid-body rotation, the boundary layer is turbulent, and buoyancy effects are active. Nevertheless, in this article we also describe a linearized Ekman-layer flow model coupled to an oceanic eddy that helps guide the interpretation of the computational results.

The outline of the paper is as follows: An analytic model for Ekman-layer winds over an oceanic eddy is developed in section 2. The LES equations and setup of the LES experiments for varying stability are described in section 3. The coupling of the ABL with an oceanic eddy and the surface boundary conditions are described in section 4. Results with winds alone (no surface heating), free convection (no mean winds), and winds plus convection are discussed in three separate sections (sections 5, 6, and 7, respectively). Section 8 is a summary of the findings.

## 2. Ekman-layer winds over an oceanic eddy

*z ϵ*[0, ∞], where

**u**

*= (*

^{h}*u*,

*υ*) is the horizontal velocity,

*w*is the vertical velocity with a surface boundary condition of

*w*(0) = 0, and

*f*and

*ν*are the Coriolis frequency and vertical eddy viscosity, respectively, both assumed constant. For our purposes we choose a surface stress boundary condition:

*ρ*is the atmospheric density and

_{a}

*τ**is the stress of the underlying ocean on the atmosphere above. The interior boundary condition is ∂*

_{s}

_{z}**u**

*→ 0 as*

^{h}*z*→ ∞. The problem is well posed as a function of

*z*for any

*τ**(*

_{s}*x*,

*y*,

*t*).

*x*and

*y*components of the vector as its real and imaginary parts. In particular,

*λ*is negative, the first term decreases with

*z*, and because the imaginary part is negative, its velocity direction rotates clockwise while ascending (a.k.a. the Ekman spiral). The second term is a constant interior velocity, presumed in geostrophic balance with an interior pressure gradient. The surface wind is rotated clockwise relative to the surface stress by

*π*/4. The complex notation facilitates recognizing the various vector directions. Additionally the horizontal transport is

*π*/2 relative to the surface stress. From (1) it is related to the vertical velocity at the top of the boundary layer (aka Ekman pumping) by

*C*:

_{D}*x*,

*y*,

*t*)-averaged value for boundary

**u**

*at the top of the layer and rotated slightly counterclockwise (i.e., by about*

_{g}*ϕ*= 10° or 0.17 radians), see section 5. Thus, in fitting an Ekman solution to this LES we choose a stress of

_{a}*z*is the boundary layer height, which in the LES solution is stratification limited by the interior inversion layer. As an example (Table 1), for

_{i}*f*= 10

^{−4}s

^{−1},

**u**

*= 5 m s*

_{g}^{−1}to the east, and

*z*= 310 m, fitted Ekman parameters are

_{i}*ν*= 5 m

^{2}s

^{−1}so that the Ekman vertical decay scale

*z*. Of course, this is only a gross fitting, because the Ekman model neglects the more realistic effects of a variable

_{i}*ν*(

*z*), buoyancy and advection. The orientations of these various quantities are shown in Fig. 1.

LES parameters for free convection, mixed shear–convective, neutral, and Ekman simulations where Free, Mixed, Neutral, and Ekman mark the simulation family, respectively. *θ*_{BC} = 0 fixed SST and *θ*_{BC} = 1 fixed

*ϕ*, consistent with (8) above. The current is assumed to be a circular eddy with amplitude

_{a}*U*(positive for a cyclone), radial shape

_{o}*F*(

*r*), and radial size

*r*; i.e., in cylindrical coordinates:

_{m}*ε*=

*U*/

_{o}*U*is a small parameter used to make this first-order approximation to (9). At leading order, this formula is the same as horizontally homogeneous (8), and at first order we see the inhomogeneous effects of the oceanic eddy ∝

_{g}*U*. The first current term ∝ −exp(

_{o}*iϕ*) has its stress in the same direction as the average stress (Fig. 1), and the second term ∝ −exp[

_{a}*i*(

*ϕ*−

*π*2)] = exp[

*i*(

*ϕ*+

*π*2)] has its stress directed along the eddy current, hence its transport

**V**outward from the eddy center for a cyclone, with downwelling in the eddy core, and vice versa for an anticyclone. Because of the imprecision in evaluating

**u**

*in the bulk formula, we will not attempt to make an accurate estimate of the magnitude of the eddy-induced atmospheric circulation but rather focus on its spatial patterns; in terms of the parameters in Table 1 for the LES cases with*

_{a}*u*= 5 m s

_{g}^{−1},

*ε*= 0.466 with a maximum eddy current speed of 1 m s

^{−1}.

*ε*. Its spatial pattern is in Fig. 2 for a cyclonic eddy. Apart from the small counterclockwise rotation by

*ϕ*, it is a dipole with enhanced stress to the north and reduced stress to the south. This is easily understood as a higher stress where the eddy current opposes the wind and a lower stress where it is in the wind direction.

_{a}*w*(

_{e}*x*,

*y*) in (6). Using

The interior downwelling region extends further outward in the east–west direction, and there are local maxima in the exterior upwelling in the north–south direction. For an anticyclonic eddy, the Ekman model predicts a sign reversal for the *w _{e}* and

Finally, although the Ekman model has no buoyancy effects, we can anticipate that the bulk formula for surface heat flux, with *r* = *r _{m}* for both cyclones and anticyclonic eddies because there is a positive stress magnitude anomaly for both.

## 3. LES details

### a. LES equations

**u**≡

*u*= (

_{i}*u*,

*υ*,

*w*) denote the Cartesian velocity components,

*θ*is virtual potential temperature and

*p*is the pressure variable normalized by air density

*ρ*. The three Cartesian coordinates are

_{a}**x**≡

*x*= (

_{i}*x*,

*y*,

*z*) also referred to as (streamwise, spanwise, vertical) directions, respectively. The set of LES equations that describe rotating stratified turbulent flow in a ABL under the incompressible Boussinesq approximation are

The above equation set includes transport equations: for momentum *ρ _{a}*

**u**, (18a), and for virtual potential temperature,

*θ*(18b);

*θ*is referred to as “temperature” in the narrative. The divergence free (incompressible) condition (18c) determines the elliptic pressure variable

*p*. Equation set (18) also includes: geostrophic winds

**u**

*= (*

_{g}*u*,

_{g}*υ*), rotation vector

_{g}**f**= (0, 0,

*f*) with Coriolis parameter

*f*, unit vector

*β*=

*g*/

*θ*where

_{r}*g*is gravity and

*θ*is the reference potential temperature. The subgrid-scale (SGS) momentum and temperature fluxes (

_{r}**T**,

**B**) are estimated using turbulent eddy viscosity prescriptions. Other details are found in S20 and S21.

### b. Mesh resolution and data sampling

Four families of LES experiments with varying geostrophic winds, domain size, and surface boundary conditions are carried out: 1) free convection, 2) mixed shear and convection, 3) shear driven neutral, and 4) Ekman see Table 1. The simulations span the stability range −*z*/*L* = [0, ∞], where *z* is the vertical coordinate and *L* is the Monin–Obukhov stability length. The LES setup including the number of grid points used for each type of simulation are listed in Table 1. The horizontal mesh resolution for all simulations is Δ(*x*, *y*) = 6.5 m, except for the Ekman simulation where Δ(*x*, *y*) = 13 m. The vertical mesh is a smoothly stretched grid built with a near unity stretching factor *K* = Δ*z _{k}*

_{+1}/Δ

*z*= 1.000 99. Here

_{k}*k*is the vertical index, and the first vertical

*w*grid point is located at Δ

*z*

_{1}= 3 m, (Fig. 5). The vertical domain is

*L*= 1.4 km and is more than 2 times the final boundary layer depth

_{z}*z*. It is convenient to position the origin of the horizontal coordinates (

_{i}*x*,

*y*) at the center of the LES domain: thus

*x*= [−

*L*/2,

_{x}*L*/2] and

_{x}*y*= [−

*L*/2,

_{y}*L*/2] where (

_{y}*L*,

_{x}*L*) denote the lengths of the horizontal domain.

_{y}All the simulations are started using the same recipe. First, a small domain LES (5, 5, 1.4) km with mesh points (*N _{x}*,

*N*,

_{y}*N*) = (768, 768, 384) is run starting with random initial conditions, fixed temperature flux

_{z}*x*and

*y*to span the larger domains in Table 1. At the same time the number of mesh points is increased, i.e., the mesh resolution is not changed. Finally the large domain fine mesh LES is started but now including the eddy currents. The simulations are run for approximately 30 000 additional steps. Data are analyzed over the last 2.8 h of the simulation when the fields are quasi-steady. Variables in Table 1 include the wind speed at the 10 m level

*u*

_{10}, friction velocity

*z*(Sullivan et al. 1998), temperature flux

_{i}*r*. To expose the eddy induced motions we use low-pass filtering, specifically a Gaussian filter applied in spectral space with filter scale

_{m}*δ*(e.g., Pope 2000, p. 563).

## 4. LES and current coupling

### a. Oceanic eddy

*F*(

*r*) and azimuthal speed given by

*r*,

*ϕ*), the Cartesian components used in the LES are (

*u*,

_{o}*υ*) =

_{o}*U*[−sin(

_{s}*ϕ*), cos(

*ϕ*)]. The parameter

*U*controls the eddy speed and rotational direction and

_{o}*r*is a measure of the eddy width. The eddy speed

_{m}*U*increases linearly with small

_{s}*r*, i.e., with solid body rotation, and is exponentially damped in the far field. The eddy speed reaches a maximum

*U*, normalized by the maximum current

_{s}*U*

_{s}_{,}

*= 1 m s*

_{m}^{−1}. To illustrate the symmetries of the circular eddy, Fig. 4 also shows the variation of the Cartesian velocity components

**u**

*= (*

_{o}*u*,

_{o}*υ*) slicing through the eddy center;

_{o}*u*for varying

_{o}*y*/

*r*at constant

_{m}*x*= 0, and

*υ*

_{0}for varying

*x*/

*r*at constant

_{m}*y*= 0. Note (

*u*< 0,

_{o}*u*> 0) in the upper and lower halves of the vortex while (

_{o}*υ*> 0,

_{o}*υ*< 0) in the right and left halves of the vortex. The sinusoidal oscillation of the (

_{o}*u*,

_{o}*υ*) velocity components in the vortex has an important impact on the results. In all simulations the eddy center is located at (

_{o}*x*,

*y*) = (0, 0). For the simulations in Table 1 the parameter

*ε*=

*U*/

_{o}*u*= [∞, 0.932, 0.466[, and

_{g}*ε*(>0, <0) for (cyclonic, anticyclonic) rotation.

### b. Surface boundary conditions

Figure 5 illustrates the staggered vertical arrangement of flow variables at two levels above the water surface used to implement the surface boundary conditions in the LES. Variables (*u*, *υ*, *p*, *θ*) are located at the (*z*_{1/2}, *z*_{3/2})-levels while variables (*w*, *e*) are located at the (*z*_{1}, *z*_{2})-levels; the (*z*_{1/2}, *z*_{3/2})-levels are centered midway between the (*z*_{1}, *z*_{2})-levels. We adopt rough-wall Monin–Obukhov (M–O) similarity formulas for the calculation of momentum and temperature fluxes at the water surface (e.g., Businger 1972). Variables (*u*, *υ*, *θ*)_{1/2} stored at the *z*_{1/2}-level and variables (*u*, *υ*, *θ*)* _{o}* stored at the water surface are both used in the M–O formulas. The imposed roughness is

*z*= 2 × 10

_{o}^{−4}m characteristic of a low or moderate wind marine boundary layer (e.g., Large and Pond 1981). The M–O formulas for friction velocity

*z*

_{1/2}/

*L*where the M–O length

*κ*= 0.4. Given (

*u*,

*υ*,

*θ*,

*z*)

_{1/2}and

*z*, a robust bisection iteration scheme solves for (

_{o}*u*,

_{o}*υ*) and imposed SST

_{o}*θ*.

_{o}

*τ*_{s},

*q*)/

_{s}*ρ*are used, for clarity density is omitted in Fig. 5 and in Eqs. (22) and (23). The velocity and temperature differences between the

_{a}*z*

_{1/2}-level and the water surface are

**u**, and the average

*C*,

_{d}*C*) are found in the bisection iteration scheme, additional details about the scheme can be found in the appendix of Sullivan et al. (2014). We emphasize that Eqs. (20a) and (20b) are coupled and nonlinear rules for surface fluxes. Although we primarily impose SST as the surface boundary condition, we also consider a boundary condition with imposed surface temperature flux. The boundary condition for

_{h}*w*is the usual no-flow condition

*w*= 0 at

*z*= 0. The M–O formulas are applied at every

*x*–

*y*grid point in the LES.

### c. Currents and turbulent motions in the boundary layer interior

*z*

_{1/2},

*z*

_{1})-levels in Fig. 5:

*u*,

*υ*,

*θ*) fields. The pressure

*p*

_{1/2}, and thus the horizontal and vertical pressure gradients, further depend on the surface fluxes as shown below.

*p*(

**x**) at each substep that satisfies the divergence free condition. To expose the role of the surface boundary conditions consider the Poisson pressure equation and its source term near the water surface in Fig. 5. At the (

*z*

_{1/2},

*z*

_{3/2})-levels the pressure satisfies the semidiscrete equation:

*w*(

*z*= 0) = 0, for further details see Sullivan et al. (2014). Inspection of (23) shows the source term in the pressure stencil depends on the horizontal gradients of the surface stress

*z*

_{1}grid level where

*θ*

_{1}= (

*θ*

_{3/2}+

*θ*

_{1/2})/2; because of (22d) temperature flux

*q*is implicit in

_{s}*θ*

_{1}and hence

*p*

_{1/2}. The eddy current

**u**

*and SST*

_{o}*θ*thus impact the pressure because of surface fluxes (20). Pressure gradients then transfer the current impact to the dynamical fields (

_{o}*u*,

*υ*,

*w*) in the interior of the boundary layer.

What are the impacts of specifying surface SST versus surface flux in (20b)? If we specify SST then *q _{s}* inherits all the variability in the surface winds and temperature, and the underlying eddy current, i.e.,

*θ*

_{1/2}according to (22d), and hence the pressure. Specifying constant surface temperature flux, however, short circuits variability in the temperature field. There are fluctuations in the exchange coefficients (

*C*,

_{d}*C*) since we use M–O similarity theory to diagnose them, but

_{h}*C*and Δ

_{h}**u**variability is balanced by changes in

*θ*in (20b), and

_{o}*θ*does not appear in (22d). Boundary conditions constant

_{o}*θ*and constant

_{o}*q*thus produce different levels of variability in

_{s}*θ*

_{1/2}, the pressure field, and hence the boundary layer. With imposed SST, both wind and temperature fluctuations (Δ

**u**, Δ

*θ*) result in heterogeneous

*q*. Broadly, incorporating oceanic currents in the LES introduces both mechanical and thermal couplings within the overlying atmospheric boundary layer.

_{s}## 5. Results: Neutral winds

We first examine the impact of winds and currents on boundary layer dynamics under neutral surface conditions, i.e., where −*z*/*L* → 0. In contrast to the simulations in sections 6 and 7, these LES experiments impose a temperature flux boundary condition *q _{s}*(

*x*,

*y*) = 0. This boundary condition isolates the coupling between surface momentum and temperature fluxes in the presence of eddy currents, and also allows a comparison with the linearized Ekman model discussed previously. In simulation k03 the geostrophic wind

*u*= 5 m s

_{g}^{−1}, the eddy radius

*r*= 4 km, and

_{m}*ϵ*= 0.466. The LES domain size, number of mesh points, and bulk parameters are provided in Table 1. Note the smaller values of friction velocity

*z*in the absence of convective forcing. Thus the ratio

_{i}*z*/

_{i}*r*is smaller than in the mixed shear-convective cases in section 7.

_{m}A robust feature of all simulations is the development of organized variability in the surface momentum flux, a flux dipole shown in Fig. 1 which is also predicted by the Ekman model. To compare with (14), we process the LES output and compute the average stress anomaly (21); to further smooth the results, a low-pass filter is applied with filter *δ* = 2 km. The theory-simulation comparison in Fig. 1 shows good qualitative and quantitative agreement. The peak magnitude of the stress anomaly occurs near *x* = 0 and *y*/*r _{m}* ∼ ±0.707, i.e., where

*U*is maximum. Simulation k01 with

_{s}*u*= 2.5 m s

_{g}^{−1}and

*ε*= 0.932 with smaller values of

*z*generates a similar flux dipole as in Fig. 1 with (positive, negative) perturbations over the (north, south) quadrants of the eddy. Anticyclonic eddies

_{i}*ε*< 0 in the LES and also in the Ekman theory reverse the sign of the anomaly. At

*z*= 7.5 m, the near surface horizontal velocity anomalies

**u**

*exhibits a cyclonic swirling motion about the location (*

_{a}*x*,

*y*)/

*r*= (0.6, 0.2) at

_{m}*z*= 7.5 m and at (

*x*,

*y*)/

*r*= (1.4, 0.5) at

_{m}*z*= 59 m, see supplementary material. With increasing height, advection continues to shift the horizontal and vertical anomaly patterns further to the northeast, and at the same time the magnitude of the anomaly decays. For

*z*> 59 m the coherent signature of the eddy disappears in the background turbulence and the horizontal motions morph into a streamwise oriented pattern.

A demanding test compares the vertical velocity *w* from the Ekman model and the LES. One qualitative difference is that the unstratified Ekman model has *w* → *w _{e}* as

*z*increases, and it extends to an arbitrary height. In contrast in the stratified LES

*w*does not penetrate beyond the inversion layer, partly because of the

*w*= 0 boundary condition at the top of the LES domain; therefore, we compare its

*w*field at an intermediate level within the ABL. Figure 3 depicts the pumping velocity

*w*from (15) with the vertical velocity anomaly computed from simulation k03 at

_{e}*z*= 60.5 m. The fields in Fig. 3 are each normalized by

*z*= 60.5 m (

*z*/

*z*= 0.19), the LES predicts the center of the downdraft is offset to the northeast and is surrounded by weaker updrafts in qualitative agreement with the Ekman model. However, the magnitude of the updraft and downdrafts is considerably weaker in LES. Tracking the vertical path of the downdraft–updraft pattern in the LES shows it is offset from the eddy center at all

_{i}*z*with maximum positive and negative magnitudes near

*z*/

*z*= 0.19. Vertical cutting planes

_{i}*x*–

*z*at

*y*/

*r*= 0.7 and

_{m}*y*–

*z*at

*x*/

*r*= 0.8 show the vertical organization of

_{m}*w*induced by the eddy is confined to the near surface. For case k01 with

*ε*= 0.933 the magnitude of the vertical velocity anomaly increases by more than a factor of 2 compared to case k03 with

*ε*= 0.466; and the pattern of downdrafts and updrafts is more coherent over the depth of the ABL. Changing the sign of eddy rotation in simulation k04 simply reversed the updraft and downdraft patterns compared to k03 but otherwise did not alter the magnitude of the vertical velocity.

The results from simulations k03 and k04 motivated further investigation. An Ekman simulation m01 with *u _{g}* = 5 m s

^{−1}and

*r*= 4 km was carried out to remove stable stratification impacts on the eddy induced motions: Recall k01, k03 and k04 are stably stratified at all

_{m}*z*but especially at

*z*levels near the ABL inversion. The LES Ekman simulation m01 develops vigorous turbulence over the depth of the ABL

*z*∼ 700 m but we readily identify the same eddy induced flux dipole and vertical velocity anomaly as in Figs. 2 and 3. Inspection of the LES Ekman solution also shows a complex transition from the downdraft–updraft pattern near the surface to an elongated streamwise pattern for

_{i}*z*> 100 m. The amplitude of the horizontal motions at

*z*> 100 m are large, nearly a factor of 5 greater than the eddy induced vertical velocity. Horizontal advection that varies with

*z*appears to be the main source of disagreement between the linearized Ekman model and the LES, not stratification influences.

## 6. Results: Free convection

Idealized free convection, no mean wind (e.g., Schmidt and Schumann 1989), is the asymptotic stability regime −*z*/*L* → ∞ frequently used to examine the impact of land surface heterogeneity on the atmospheric boundary layer (e.g., Raasch and Harbusch 2001; Owinoh et al. 2005; Patton et al. 2005). Flow visualization and analysis of our LES runs also shows unique features in the free convection regime when coupled to an oceanic eddy with heterogeneous currents.

As discussed in sections 4b and 4c even with fixed SST, the surface temperature flux depends on the eddy currents and overlying turbulent winds because of coupling in the surface flux formulas (20). A cyclonic eddy with radius *r _{m}* = 400 m interacting with convective turbulence generates rich variability in the surface flux as shown in the left panel of Fig. 7. Close inspection shows a high concentration of positive temperature flux at (

*x*,

*y*) = (0, 0), the instantaneous flux at the eddy center frequently exceeds 60% of the average flux 〈

*q*〉. In regions outside of the eddy the flux is often less than 50% of 〈

_{s}*q*〉. Coupling the currents with the overlying winds also results in spatially intermittent momentum flux. The right panel of Fig. 7 shows the distribution of radial and azimuthal momentum fluxes (

_{s}*τ*,

_{r}*τ*) near the eddy center. Notice (

_{ϕ}*τ*,

_{r}*τ*), computed from

_{ϕ}Eddy currents coupled to surface fluxes in free convection leave an imprint on the full atmospheric boundary layer particularly in the pressure and temperature fields and to a lesser degree the vertical velocity field, see Figs. 8 and 9. The (*p*, *θ*) fields show horizontal organization at a scale larger the boundary layer depth *z _{i}* and eddy radius

*r*. Based on the decay of the eddy current in (19),

_{m}*U*/

_{s}*U*

_{s}_{,}

*≈ 0.06 at*

_{m}*r*/

*r*= 2.1 and thus we expect to observe an impact on

_{m}*θ*and

*p*at scales with an upper bound

*r*∼ 1 km. At the water surface, the coherent negative pressure at the eddy center is nearly 1 km wide and is surrounded by an outer ring of less organized positive pressure. The coupling between large horizontal scales in

*p*and

*θ*at the surface is contained in (23). We find that the LES requires a horizontal domain of 10 km in order to resolve the pressure field organization; coherent pressure first begins to develop midway through the simulation. An

*x*–

*z*plane cutting through the eddy center, see the bottom panel of Fig. 8, shows that the pressure organization extends over the full depth of the boundary layer, a sign change (positive pressure) occurs in the middle and upper boundary layer.

At the surface the *θ* field exhibits two scales, a wide slightly warm core overlain with a web of warmer narrow intersecting filaments; the former is induced by the oceanic eddy while the latter is a characteristic convection feature found over a homogeneous surface (e.g., Schmidt and Schumann 1989; Sullivan and Patton 2011). Near the boundary layer top a strong and wide central updraft, induced by a secondary circulation, is surrounded by clusters of smaller plumes arranged in a ring of radius ∼ 2 km. The no flow boundary condition on vertical velocity prevents *w* from developing the larger horizontal scale observed in the *θ* field near the surface. Depending on the horizontal position of an *x*–*z* or *y*–*z* cutting plane we also observe intense splats in *u* or *υ* at the ABL inversion. Similar results are obtained for cyclonic and anticyclonic eddies, i.e., changing the sign of the eddy rotation does not change the sign of the ( *p*, *θ*) fields. Apparently oceanic eddies can act as an agent to organize structures over the depth of the boundary layer in free convection.

It is interesting to contrast our convective boundary layer results with analytic models of neutral rotating boundary layers (e.g., Schlichting 1968; Rotunno 2013, 2014). To facilitate the comparison the LES surface fluxes *u*, *υ*) are first converted into radial and azimuthal components (*τ _{r}*,

*τ*,

_{ϕ}*υ*,

_{r}*υ*) at each (

_{ϕ}*x*,

*y*) location. These variables and

*w*are then interpolated to a fine polar (

*r*,

*ϕ*) mesh, and integrated around a ring

*ϕ*= [0, 2

*π*] at fixed values of

*r*. Finally, these results are averaged in time; square brackets [ ] in the narrative below indicate a time average in (

*r*,

*ϕ*) polar coordinates.

Results for the surface fluxes [*τ _{r}*,

*τ*,

_{ϕ}*q*] are shown in Fig. 10 for cyclonic and anticyclonic eddies. Notice the temperature and momentum fluxes peak near

_{s}*r*= 0.7

*r*where the eddy current is maximum. Fluxes [

_{m}*τ*,

_{ϕ}*q*] decay with increasing

_{s}*r*>

*r*. Noteworthy is that the flux

_{m}*τ*tends to

_{r}*r*, thus the radial flow near the surface is inward toward the low pressure region in Fig. 8. The spatial variation of the fluxes is similar for cyclonic and anticyclonic eddies, but with a sign change in

*τ*, the sign of

_{ϕ}*τ*matches the eddy rotation. The magnitude of the temperature flux at

_{ϕ}*r*/

*r*= 1 depends on the eddy rotation; [

_{m}*q*] is noticeably larger for the anticyclonic case compared to its cyclonic counterpart. Recall the linear Ekman model [Eq. (15), section 2] predicts a central updraft for anticyclonic rotation. In the LES there is cooperative behavior between surface convection and anticyclonic eddy rotation that enhances the surface temperature flux.

_{s}Average velocity fields [*υ _{r}*,

*υ*,

_{ϕ}*w*](

*r*,

*z*) are shown in Fig. 11 at four different vertical levels

*z*= (7.5, 184.0, 409.0, 646.0) m for simulations with cyclonic and anticyclonic eddies. Inspection of the results shows that the oceanic eddy has generated a large scale coherent SC under free convection. The radial flow

*υ*< 0 is inward toward the eddy center over the bulk of the boundary layer

_{r}*r*/

*r*= [0, 12] and

_{m}*z*/

*z*< 0.5. Vertical velocity [

_{i}*w*] > 0 is concentrated in an updraft near the eddy center, the strength of the updraft increases with vertical distance from the surface as expected. In the middle of the boundary layer, [

*w*] > (

*υ*,

_{r}*υ*) as a result of the potent SC. In the inversion layer the updraft

_{ϕ}*w*is quenched and the resulting splat induces outward radial motion

*υ*> 0. The sign of the azimuthal velocity

_{r}*υ*in the lower boundary layer depends on the direction of the eddy rotation. The surface stress

_{ϕ}*τ*acts to

_{ϕ}*drag*the surface winds in the same direction as the eddy rotation, i.e., similar to a plane Couette flow. In the middle and upper boundary layer

*υ*is weak but greater than zero. In the anticyclonic case

_{ϕ}*υ*switches sign near

_{ϕ}*z*/2, and we speculate the sign switch is a consequence of positive Coriolis rotation. The sign of the central updraft is independent of the sign of the rotation, i.e., [

_{i}*w*] > 0 independent of the rotation, markedly different than the neutral regime in section 5. However, in the interior of the ABL, vertical velocity [

*w*] for 0 <

*r*/

*r*< 1.5 is noticeably enhanced with anticyclonic eddy rotation. As mentioned previously, this is a signature of cooperative behavior between surface convection and the central updraft induced by the eddy rotation.

_{m}Thus a circular eddy of modest strength generates a SC in the free-convective atmospheric boundary layer that draws in fluid at the bottom and expels fluid vertically upward into the interior of the boundary layer. Near the top of the boundary layer the rising plume splashes against the stably stratified inversion and the splat generates a radial outflow. The outflow is deflected downward and the descending shell away from the eddy center feeds the radially inward flow at the surface completing the circulation. Thus the surface eddy and convection act in concert to create a secondary circulation. In a neutral cyclonic rotating boundary layer discussed in section 2 and by Schlichting (1968), Rotunno (2013, 2014) the SC features a downdraft at the center with outward radial flow near the surface. Convective boundary layer turbulence coupled to an oceanic eddy produces a SC that is opposite in sign and different compared to a neutral rotating boundary layer. The results depend on the sign of the eddy rotation and the surface boundary condition for temperature.

In a two-step process, eddy currents first generate spatial heterogeneity in the surface temperature flux with fixed SST. Then the resulting heterogeneity in *q _{s}*(

*x*,

*y*) acts to drive organized depth-filling secondary circulations similar to atmospheric flows with heterogeneous SST fronts as shown in S20 and S21. A surface boundary condition that holds

*q*spatially constant eliminates surface heterogeneity and SC even in the presence of currents. To test these ideas simulation f05 imposes a flux boundary condition with no currents, but uses a point-by-point time average flux

_{s}*p*,

*θ*,

*w*) fields from f05 match those obtained from f01 (not shown). In other words, spatial heterogeneity in

*q*(

_{s}*x*,

*y*) generated by coupling eddy currents in the surface flux formulas is responsible for the observed organization in free convection. Simulation f08 with a 10-times-larger radius

*r*= 4 km generates a similar impact on the ABL as

_{m}*r*= 0.4 km; see supplementary material.

_{m}## 7. Results: Winds and convection

To examine the combined impact of winds, convection, and currents on boundary layer dynamics in an intermediate regime of −*z*/*L* a suite of LES experiments are carried out varying the geostrophic wind *u _{g}* = (2.5, 5) m s

^{−1}and eddy radius

*r*= (0.4, 1.2, 2, 4) km. In these simulations the SST is adjusted to generate an average temperature flux

_{m}*q*similar to that in the free convection experiments. The LES details and bulk parameters are provided in Table 1. We focus on simulations with

_{s}*r*= 4 km, but touch on the impact of smaller

_{m}*r*. For discussion purposes it is convenient to divide the circular eddy into four quadrants (left, upper, right, lower) referred to as (west, north, east, south) quadrants, respectively.

_{m}### a. Surface flux anomalies

In the intermediate regime of −*z*/*L* an oceanic eddy generates flux anomalies that differ from the coherent structures in a homogeneous flow, and are also different than the anomalies in the free convection and neutral regimes discussed previously. Typical results for (*τ _{a}*,

*q*) from cases i04 (no eddy) and i03 (with an eddy) are compared in Fig. 12. In both simulations the winds are light

_{a}*u*= 2.5 m s

_{g}^{−1}and the average fluxes

*x*axis. In i03 the eddy rotation is cyclonic,

*r*= 2 km, and the ratio of the boundary layer depth to eddy radius

_{m}*z*/

_{i}*r*= 0.35. To emphasize the impact of the oceanic eddy on the magnitude and scale content of the anomalies, the average fluxes computed from (21), shown in Fig. 12, are further low-pass filtered using a Gaussian filter with scale

_{m}*δ*= 500 m.

In the horizontally homogeneous case i04 (no eddy), the momentum and temperature flux perturbations are streaky anisotropic structures mainly aligned with the streamwise wind at all scales. ABL turbulence generates flux fluctuations mainly at scales less than 500 m; (*τ _{a}*,

*q*) are less than 10% of

_{a}*r*.

_{m}The momentum flux organization is explained by the wind and current alignment. Given *u*_{1/2} > 0, in the north quadrant *u _{o}* < 0 and then

*u*> 0 and then

_{o}*υ*

_{1/2}> 0, Δ

*υ*is positive and large in the west quadrant and small or negative in the east quadrant.

The new feature in the mixed shear-convection regime is the formation of a temperature flux dipole. Even with spatially uniform *θ _{o}*(

*x*,

*y*) a temperature flux dipole develops because of wind–current coupling

*θ*< 0,

*q*in the (north, south) quadrant is (larger, smaller), respectively. In the south quadrant wind–current coupling nearly switches the temperature flux to stable forcing, i.e., locally the M–O stability parameter −

_{s}*z*/

*L*is very unstable in the north quadrant while in the south quadrant −

*z*/

*L*tends to neutral. The eddy radius

*r*plays an important role in the dipole pattern as the latter is enhanced with larger

_{m}*r*, i.e., as

_{m}*r*/

_{m}*z*increases the impact of the heterogeneous currents also increases. The sign and spatial dipole patterns in the bottom panels of Fig. 12 reverse with an anticyclonic eddy (not shown). Thus oceanic eddy currents can modify both momentum and temperature fluxes in an intermediate regime of −

_{i}*z*/

*L*. The surface fluxes are (largest, smallest) where winds and currents are (opposed, aligned). For the low wind case i03, the normalized peak (stress, temperature flux) anomalies are approximately (0.7, 0.3) or about (0.0074 N m

^{−2}, 4.8 W m

^{−2}), respectively.

All simulations with an eddy generate flux patterns qualitatively similar to those shown in the bottom panels of Fig. 12. To compare the main features of the flux dipoles a spanwise cut through the anomalies at *x* = 0 is collected from each simulation and presented in Fig. 13; the cut at *x* = 0 passes near the maximum and minimum in the dipoles. In each case, the maximum (positive) magnitude occurs near *y* = 0.7 *r _{m}*, i.e., where the eddy current

*u*/

_{o}*U*< 0 is a negative maximum (Fig. 4). The magnitude of

_{o}*q*is consistently smaller by a factor of 3 compared to

_{a}*τ*. For a fixed eddy radius, the magnitude of

_{a}*τ*increases with decreasing wind speeds. For example comparing cases (i07, i05) with

_{a}*r*= 2 km the magnitude of

_{m}*τ*increases from 0.4 to 0.8 as

_{a}*u*decreases from 5 to 2.5 m s

_{g}^{−1}. This approximate linear dependence is predicted by the linearized Ekman theory, recall the eddy impact varies with

*ε*=

*U*/

_{o}*u*. Convection does not markedly alter the observed patterns, an 83% increase in surface convection in case h04 leads to (normalized) anomalies nearly identical to i07 (not shown). For a fixed wind speed, notice the scale content of (

_{g}*τ*,

_{a}*q*) widens with increasing

_{a}*r*; the LES are mainly in the regime

_{m}*z*/

_{i}*r*≤ 1. In other words the scale content of the surface flux which the ABL “feels” widens as

_{m}*r*increases.

_{m}It is important to notice that the spanwise gradients *r _{m}* at fixed wind speed, and thus we anticipate that the impact of current heterogeneity also increases as described in S20 and S21. With eastward

*u*, the surface wind is mainly aligned with the

_{g}*x*direction but has a small ageostrophic component in the

*y*direction. This combination of surface winds and spanwise gradients in the anomalies is in the “alongfront” regime discussed in S21. With alongfront winds, pressure gradients are equally or more important than mean horizontal advection and the induced secondary circulation is in a plane perpendicular to the alongfront wind direction. However, the sign of the secondary flow in the present simulations depends on the sign of the eddy rotation, see Fig. 13.

### b. Boundary layer interior

With mixed shear and convection an oceanic eddy generates organized motions in the turbulent fields near the surface but also persistent patterns in the interior of the ABL often extending up to the inversion. The eddy imprint in the shear-convection regime is not as dramatic as in free convection but is clearly visible in the temperature field. Eddy induced patterns strengthen with increasing *r _{m}* and decreasing wind

*u*and are noticeably different than observed in free convection. Visualization in horizontal planes near the surface and ABL inversion shows two scales in

_{g}*θ*, an extended streamwise streaky pattern with spanwise width ≈

*r*. And the streamwise pattern is punctuated with smaller scale updrafts and updrafts associated with boundary layer convection. In the case of a cyclonic eddy, the updrafts and downdrafts are stronger and weaker over the north and south sides of the eddy, see Fig. 12. Very near the surface the (

_{m}*u*,

*υ*,

*θ*,

*p*) patterns mimic the closed dipoles observed in the surface fluxes, e.g., Fig. 12. With increasing

*z*the eddy induced motions weaken against the background turbulence as the spatial patterns broaden in

*x*and morph into elongated ellipses centered over the surface anomalies. The major axis of the ellipse is aligned with the geostrophic wind, and tends to span the entire

*x*extent of the LES domain.

*y*–

*z*plane, i.e., in the plane perpendicular to the geostrophic wind. An anomaly for a field

*h*is computed from

*t*and

*x*over the limited interval [−0.8, 1.5]

*r*: 〈

_{m}*h*〉

*is the far-field upstream average. To extract the induced motions a Gaussian filter with scale*

_{f}*δ*= 1 km is finally applied in the

*y*direction at each

*z*. Results for simulations with

*u*= 2.5 m s

_{g}^{−1}, i05 cyclonic eddy

*r*= 4 km, and i06 anticyclonic eddy

_{m}*r*= 4 km are shown in the four-panel plots in Figs. 14 and 15, respectively. Each figure shows four fields (

_{m}*υ*,

*w*,

*θ*,

*p*)

*. For easy comparison the figures use the same color bar ranges, and the fields are dimensionless using proper combinations of (*

_{a}*x*–

*y*plane. The impact of an oceanic eddy is clearly visible in Figs. 14 and 15. An eddy with horizontal scale

*r*>

_{m}*z*sets the horizontal scale of persistent anomalies (

_{i}*υ*,

*θ*,

*p*)

*that extend over the depth of the ABL; the eddy also appears to provide a mechanism to excite stronger vertical motions in mid-ABL. Notice the anomalies are responsive to the sign of the eddy rotation, for example, near the surface*

_{a}*p*and ∂

_{a}*reverse sign for cyclonic and anticyclonic eddies. Closer inspection of the results shows that the pressure gradient −∂*

_{y}p_{a}*drives a (positive, negative) spanwise anomaly*

_{y}p_{a}*υ*near the (surface, inversion), respectively. The pressure field responds to the baroclinic changes in temperature. Notice the underlying eddy motion produces large-scale horizontal variability in

_{a}*θ*near the inversion which can influence entrainment. As mentioned previously, the orientation of the geostrophic wind relative to the eddy-induced temperature heterogeneity in the present flows is analogous to the “alongfront wind” regime examined in S21. With alongfront winds pressure gradients are of importance in generating organized secondary circulations. The eddy induced motions decrease with

_{a}*ε*=

*U*/

_{o}*u*, e.g., the motions are weak in simulation i07 with

_{g}*u*= 5 m s

_{g}^{−1}which is anticipated given the variation of surface fluxes in Fig. 13.

Two-dimensional streamtraces (average particle paths) in a *y*–*z* plane at *x*/*r _{m}* = 4.375 constructed from anomalies (

*υ*,

*w*)

*illustrate the formation of a persistent SC in the interior of the ABL downstream of the eddy center; see supplementary material. The center of the SC is positioned off axis above the north side of the eddy for both cyclonic and anticyclonic eddies, and the spanwise extent of the eddy, approximately 4 km, matches the eddy radius*

_{a}*r*. Weaker circulations (not shown) are also present left and right of the dominant SC and are believed to be residuals from limited averaging.

_{m}### c. Eddy “wake”

The vertical coherence of the eddy-induced motions found in section 7b is surprising and motivated additional investigation focusing on their horizontal coherence. Because of the horizontal periodicity in the LES, simulations in Table 1 are essentially flow over a periodic train of eddies. To examine horizontal coherence in the ABL above an isolated (nonperiodic) oceanic eddy we utilize the fringe method introduced in S20 and S21. The fringe technique utilizes small and large LES domains. The turbulent inflow to the large LES domain is nudged to match the homogeneous fields in the small LES domain effectively eliminating periodicity effects, see S20 for details. The forcing conditions are geostrophic wind *u _{g}* = 2.5 m s

^{−1}with SST set to match simulation i03. In the large LES domain a cyclonic eddy with

*r*= 2 km is located 3

_{m}*r*downstream (to the right) of the inflow. The horizontal dimensions of the large domain are (20, 10) km utilizing the spatial resolution in Table 1, i.e., a mesh with (3072, 1536) points.

_{m}A key result is the development of an eddy “temperature wake” in the boundary layer interior, see Fig. 16. The field is noisy as only point-by-point time averaging is used to construct the anomaly. Near the ABL top *z* = 398 m the temperature wake is first visible aft of the eddy center. At *z* = 28.6 m the magnitude of the anomaly pattern above the eddy is large and exhibits a dipole structure similar to the surface flux pattern in Fig. 12. However, the *θ* field at *z* = 28.6 m also exhibits a coherent trailing wake downstream of the eddy. The downstream vertical tilt of the wake between the surface and ABL inversion is a consequence of *z*-varying horizontal advection. Flow visualization shows the temperature wake persists for at least 5 radii downstream of the eddy. A vertical slice in an *y*–*z* plane at *x*/*r _{m}* = 6.875

*r*, shows clear rotation and a persistent secondary circulation generated by the vector field (

_{m}*υ*,

*w*)

*, see supplementary material. Notice the center of the rotation is asymmetric about*

_{a}*y*= 0 matching the asymmetric heterogeneity in the surface temperature flux,

*q*is stronger on the north side of the eddy. The streamtraces are found from anomaly fields averaged over 2.5 km in

_{s}*x*. The streamwise anomaly

*y*–

*z*plane and the eddy wake is three-dimensional. The structure and coherence of eddy induced wakes in the ABL are a topic for future investigation.

### d. Wind work

*F*,

_{a}K_{a}*F*). Thus, let

_{o}K_{o}*denotes wind work in the far field. By construction these anomalies vanish away from the eddy. The opposite sign in*

_{f}*F*assures that positive values in both wind-work expressions have the sense of energy transfer into, respectively, the atmosphere or the ocean. In (25a), the surface wind

_{o}K_{o}**u**

*is evaluated at the first model level*

_{a}*z*

_{1/2}(Fig. 5). The above correlations are next rewritten in terms of a temporal mean and fluctuation, i.e.,

**u**

*is constant in time; therefore, the turbulence term in*

_{o}*F*is zero.

_{o}K_{o}Contour plots of the eddy contribution to (*F _{a}K_{a}*,

*F*) from simulation i03 are shown in Fig. 17, where the fields are normalized by

_{o}K_{o}*x*–

*y*location are obtained from 22 surface volumes spanning approximately 6000 s. Results are smoothed with a spatial low-pass Gaussian filter with cutoff

*δ*= 500 m. For reference the normalized average far-field atmospheric wind work

*F*on the north side of the cyclonic eddy clearly show the enhanced “eddy-killing” mechanism described by Renault et al. (2016) where the ocean on average is losing its eddy kinetic energy back into the atmosphere. The computed far-field wind work in the ocean for our circular eddy is near zero as expected with no far-field currents. The eddy impact on wind work in the atmosphere is interesting. Inside the eddy, wind work is large with peak values approaching ±50% of the far-field atmospheric wind work, which translates to 0.012 W m

_{o}K_{o}^{−2}, a small value compared to the temperature flux anomaly of 4.8 W m

^{−2}. Note that the wind work is (from, to) the atmosphere on the (north, south) sides of the eddy. In contrast to the ocean, the (positive, negative) values of wind work in the atmosphere are approximately in balance, and the net wind work is nearly zero. The peak magnitude in |

*F*| ≈ 14 is not very different from the peak magnitude in |

_{a}K_{a}*F*| ≈ 19, and the patterns in both are dipolar.

_{o}K_{o}## 8. Summary and discussion

High Reynolds-number, large-eddy simulation (LES) is used to examine stratified marine atmospheric boundary layers (ABL) forced by a combination of geostrophic winds and convection above an oceanic eddy of varying radius *r _{m}* = [0.4–4] km; the maximum azimuthal eddy speed is 1 m s

^{−1}. The simulations span the unstable regime −1/

*L*= [0, ∞] where

*L*is the Monin–Obukhov (M–O) stability parameter. Eddy currents are coupled to the ABL in the M–O surface flux formulas.

A linearized Ekman model is developed to help guide and interpret the LES. In contrast to the classical Ekman model, the present formulation couples ABL winds to a nonlinear rough-wall stress boundary condition that depends on the air–sea relative velocity Δ**u**. As a result, the Ekman model shows how surface stress develops a dipole anomaly under the action of neutral winds. The dipole is understood as a consequence of surface winds aligned or opposing the surface currents. For a cyclonic eddy the stress anomaly is (positive, negative) in the (north, south) quadrants of the eddy. In good agreement with the Ekman model, the LES also develops a dipole in the surface stress. However, the Ekman model for a cyclonic eddy predicts a strong downdraft magnitude for vertical velocity in the ABL interior which is not found in the simulations. In the absence of an accompanying convective surface flux, the LES predicts weaker downdrafts and an asymmetric *w* pattern because of horizontal wind advection; the updraft–downdraft pattern for *w* is reversed for an anticyclonic eddy with neutral winds in both the Ekman model and LES, see sections 2 and 5.

In the free-convection regime (no mean wind), an oceanic eddy with *r _{m}* = [0.4, 4] km is particularly impactful on the ABL. A wide vigorous updraft is found above the eddy center and persists over the depth of the ABL. The radial wind is inward at the surface and outward at the ABL top. The azimuthal wind aligns with the eddy rotation either cyclonic or anticyclonic, and there is a positive

**u**

*= 0) predicts a central updraft (downdraft) for a cyclonic (anticyclonic) eddy; clearly the surface heat flux feedback in free convection overwhelms this predicted neutral circulation. There is clear cooperative behavior between surface convection and the central updraft induced by anticyclonic eddy rotation; in the region near the eddy core 0 <*

_{g}*r*/

*r*< 1.5 both the temperature flux and the vertical velocity in the ABL interior are enhanced compared to simulations with a cyclonic eddy, see section 6.

_{m}The new feature in the mixed wind and convection regime is the development of a dipole in surface temperature flux that coexists with the surface stress anomaly. Even with fixed (constant) SST a temperature flux anomaly appears because of current coupling *θ _{a}*, pressure

*p*, and spanwise velocity

_{a}*υ*anomaly fields are sensitive metrics of an oceanic eddy impact. For large

_{a}*r*the temperature field is persistently heterogeneous at scales larger than the ABL depth

_{m}*z*, i.e.,

_{i}*z*/

_{i}*r*< 1. With geostrophic winds oriented in the

_{m}*x*direction the surface temperature flux develops persistent gradients ∂

*, which in turn drive pressure gradients ∂*

_{y}θ_{a}*near the surface and inversion. ABL winds over an oceanic eddy are in the “alongfront” wind regime Sullivan et al. (2021). Well organized patterns in the spanwise anomaly*

_{y}p_{a}*υ*are found near the surface and at the top of the boundary layer in the low-wind regime with

_{a}*u*= 2.5 m s

_{g}^{−1}combined with convection at

*r*> 2 km. The spanwise horizontal scale of

_{m}*υ*matches the variation in the anomaly (

_{a}*θ*,

_{a}*p*) fields. Simulations past an isolated eddy develop a downstream “wake” and secondary circulations that fill the boundary layer vertically and persists for more than 5 radii downstream of the eddy. The anomalous wind work associated with the current feedback acts to deplete the oceanic kinetic energy (“eddy killing”) but only spatially rearranges the atmospheric kinetic energy balance, see section 7.

_{a}In the submesoscale ocean regime with horizontal scales 0.1 < ℓ < 10 km, an ocean eddy impacts the ABL in low-wind convective conditions: Edson et al. (2007) documents typical low-wind conditions. The strength and scale of the induced secondary circulations and modification to the surface fluxes grows with increasing eddy radius and speed. The LES and the linearized Ekman model show the importance of coupling currents in both surface momentum and temperature flux rules. Also, the LES results provide a target for focused surface layer and boundary layer observations with varying atmospheric stability above a heterogeneous ocean. The impact of mesoscale eddies with *r _{m}* > 20 km using LES and coupled atmosphere–ocean models is a topic for future work.

## Acknowledgments.

PPS and JCM were supported by the Office of Naval Research through the Physical Oceanography Program Awards N00014-17-1-2334 and N00014-18-1-2599, and by the National Oceanic and Atmospheric Administration Award NA19OAR4310378. PPS acknowledges support from the National Science Foundation and the Geophysical Turbulence Program at the National Center for Atmospheric Research (NCAR). This research benefited greatly from computer resources provided by the Department of Defense High Performance Computing Modernization Program and NCAR’s Computational and Information Systems Laboratory (doi:10.5065/D6RX99HX) sponsored by the National Science Foundation.

## Data availability statement.

Simulation data are available on request.

## REFERENCES

Ayet, A., N. Rascle, B. Chapron, F. Couvreux, and L. Terray, 2021: Uncovering air-sea interaction in oceanic submesoscale frontal regions using high-resolution satellite observations.

,*U.S. Clivar Var.***19**, 10–17.Businger, J. A., 1972: Turbulent transfer in the atmospheric surface layer.

, D. A. Haugen, Ed., Boston, MA, Amer. Meteor. Soc., 67–100.*Workshop on Micrometeorology*Chelton, D. B., M. G. Schlax, and R. M. Samelson, 2011: Global observations of nonlinear mesoscale eddies.

,*Prog. Oceanogr.***91**, 167–216, https://doi.org/10.1016/j.pocean.2011.01.002.Chen, S. S., J. F. Price, W. Zhao, M. A. Donelan, and E. J. Walsh, 2007: The CBLAST-Hurricane program and the next-generation fully coupled atmosphere–wave–ocean models for hurricane research and prediction.

,*Bull. Amer. Meteor. Soc.***88**, 311–318, https://doi.org/10.1175/BAMS-88-3-311.Conzemius, R. J., and E. Fedorovich, 2006: Dynamics of sheared convective boundary layer entrainment. Part I: Methodological background and large-eddy simulations.

,*J. Atmos. Sci.***63**, 1151–1178, https://doi.org/10.1175/JAS3691.1.Deardorff, J. W., 1970: A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers.

,*J. Fluid Mech.***41**, 453–480, https://doi.org/10.1017/S0022112070000691.Edson, J., and Coauthors, 2007: The Coupled Boundary Layers and Air-Sea Transfer Experiment in low winds.

,*Bull. Amer. Meteor. Soc.***88**, 341–356, https://doi.org/10.1175/BAMS-88-3-341.Frenger, I., N. Gruber, R. Knutti, and M. Münnich, 2013: Imprint of Southern Ocean eddies on winds, clouds and rainfall.

,*Nat. Geosci.***6**, 608–612, https://doi.org/10.1038/ngeo1863.Garratt, J. R., 1992:

Cambridge University Press, 316 pp.*The Atmospheric Boundary Layer.*Gaube, P., D. B. Chelton, P. G. Strutton, and M. J. Behrenfield, 2013: Satellite observations of chlorophyll, phytoplankton biomass, and Ekman pumping in nonlinear mesoscale eddies.

,*J. Geophys. Res. Oceans***118**, 6349–6370, https://doi.org/10.1002/2013JC009027.Gaube, P., D. B. Chelton, R. M. Samelson, M. G. Schlax, and L. W. O’Neil, 2015: Satellite observations of mesoscale eddy-induced Ekman pumping.

,*J. Phys. Oceanogr.***45**, 104–132, https://doi.org/10.1175/JPO-D-14-0032.1.Large, W. G., and S. Pond, 1981: Open ocean flux measurements in moderate to strong winds.

,*J. Phys. Oceanogr.***11**, 324–336, https://doi.org/10.1175/1520-0485(1981)011<0324:OOMFMI>2.0.CO;2.Marchesiello, P., J. C. McWilliams, and A. Shchepetkin, 2003: Equilibrium structure and dynamics of the California Current System.

,*J. Phys. Oceanogr.***33**, 753–783, https://doi.org/10.1175/1520-0485(2003)33<753:ESADOT>2.0.CO;2.McWilliams, J. C., 2016: Submesoscale currents in the ocean.

,*Proc. Roy. Soc. London***472**, 1–32, https://doi.org/10.1098/rspa.2016.0117.McWilliams, J. C., 2020: Oceanic frontogenesis.

,*Annu. Rev. Mar. Sci.***13**, 227–253, https://doi.org/10.1146/annurev-marine-032320-120725.Minobe, S., A. Kuwano-Yoshida, N. Komori, S. Xie, and R. J. Small, 2008: Influence of the Gulf Stream on the troposphere.

,*Nature***452**, 206–209, https://doi.org/10.1038/nature06690.Moeng, C.-H., and P. P. Sullivan, 1994: A comparison of shear and buoyancy driven planetary-boundary-layer flows.

,*J. Atmos. Sci.***51**, 999–1022, https://doi.org/10.1175/1520-0469(1994)051<0999:ACOSAB>2.0.CO;2.Owinoh, A., J. Hunt, A. Orr, P. Clark, R. Klein, H. Fernando, and F. Nieuwstadt, 2005: Effects of changing surface heat flux on atmospheric boundary-layer flow over flat terrain.

,*Bound.-Layer Meteor.***116**, 331–361, https://doi.org/10.1007/s10546-004-2819-z.Patton, E. G., P. P. Sullivan, and C.-H. Moeng, 2005: The influence of idealized heterogeneity on wet and dry planetary boundary layers coupled to the land surface.

,*J. Atmos. Sci.***62**, 2078–2097, https://doi.org/10.1175/JAS3465.1.Pope, S. B., 2000:

Cambridge University Press, 771 pp.*Turbulent Flows.*Raasch, S., and G. Harbusch, 2001: An analysis of secondary circulations and their effects caused by small-scale surface in homogeneities using large-eddy simulation.

,*Bound.-Layer Meteor.***101**, 31–59, https://doi.org/10.1023/A:1019297504109.Renault, L., M. J. Molemaker, J. C. McWilliams, A. F. Shchepetkin, F. Lemarié, D. Chelton, S. Illig, and A. Hall, 2016: Modulation of wind work by oceanic current interaction with the atmosphere.

,*J. Phys. Oceanogr.***46**, 1685–1704, https://doi.org/10.1175/JPO-D-15-0232.1.Renault, L., J. McWilliams, and J. Gula, 2018: Dampening of submesoscale currents by air-sea stress coupling in the Californian Upwelling System.

,*Sci. Rep.***8**, 13388, https://doi.org/10.1038/s41598-018-31602-3.Renault, L., P. Marchesiello, S. Masson, and J. C. McWilliams, 2019: Remarkable control of western boundary currents by eddy killing, a mechanical air-sea coupling process.

,*Geophys. Res. Lett.***46**, 2743–2751, https://doi.org/10.1029/2018GL081211.Robinson, W. P., P. Chang, E. Chassignet, and S. Speich, 2019: Ocean mesoscale eddy interactions with the atmosphere. U.S. CLIVAR and CLIVAR Workshop Tech. Rep. 2019-2, 22 pp., https://doi.org/10.5065/ebjm-5q77.

Rotunno, R., 2013: The fluid dynamics of tornados.

,*Annu. Rev. Fluid Mech.***45**, 59–84, https://doi.org/10.1146/annurev-fluid-011212-140639.Rotunno, R., 2014: Secondary circulations in rotating-flow boundary layers.

,*Aust. Meteor. Oceanogr. J.***64**, 27–35, https://doi.org/10.22499/2.6401.004.Schlichting, H., 1968:

6th ed. McGraw-Hill, 748 pp.*Boundary-Layer Theory.*Schmidt, H., and U. Schumann, 1989: Coherent structure of the convective boundary layer derived from large-eddy simulations.

,*J. Fluid Mech.***200**, 511–562, https://doi.org/10.1017/S0022112089000753.Seo, H., A. J. Miller, and J. R. Norris, 2016: Eddy–wind interaction in the California current system: Dynamics and impacts.

,*J. Phys. Oceanogr.***46**, 439–459, https://doi.org/10.1175/JPO-D-15-0086.1.Sullivan, P. P., and E. G. Patton, 2011: The effect of mesh resolution on convective boundary-layer statistics and structures generated by large-eddy simulation.

,*J. Atmos. Sci.***68**, 2395–2415, https://doi.org/10.1175/JAS-D-10-05010.1.Sullivan, P. P., C.-H. Moeng, B. Stevens, D. H. Lenschow, and S. D. Mayor, 1998: Structure of the entrainment zone capping the convective atmospheric boundary layer.

,*J. Atmos. Sci.***55**, 3042–3064, https://doi.org/10.1175/1520-0469(1998)055<3042:SOTEZC>2.0.CO;2.Sullivan, P. P., J. C. McWilliams, and E. G. Patton, 2014: Large eddy simulation of marine boundary layers above a spectrum of moving waves.

,*J. Atmos. Sci.***71**, 4001–4027, https://doi.org/10.1175/JAS-D-14-0095.1.Sullivan, P. P., J. C. McWilliams, J. C. Weil, E. G. Patton, and H. J. S. Fernando, 2020: Marine atmospheric boundary layers above heterogeneous SST: Across-front winds.

,*J. Atmos. Sci.***77**, 4251–4275, https://doi.org/10.1175/JAS-D-20-0062.1.Sullivan, P. P., J. C. McWilliams, J. C. Weil, E. G. Patton, and H. J. S. Fernando, 2021: Marine atmospheric boundary layers above heterogeneous SST: Alongfront winds.

,*J. Atmos. Sci.***78**, 3297–3315, https://doi.org/10.1175/JAS-D-21-0072.1.