In computational simulations of an idealized subtropical eastern boundary upwelling current system, similar to the California Current, a submesoscale transition occurs in the eddy variability as the horizontal grid scale is reduced to O(1) km. This first paper (in a series of three) describes the transition in terms of the emergent flow structure and the associated time-averaged eddy fluxes. In addition to the mesoscale eddies that arise from a primary instability of the alongshore, wind-driven currents, significant energy is transferred into submesoscale fronts and vortices in the upper ocean. The submesoscale arises through surface frontogenesis growing off upwelled cold filaments that are pulled offshore and strained in between the mesoscale eddy centers. In turn, some submesoscale fronts become unstable and develop submesoscale meanders and fragment into roll-up vortices. Associated with this phenomenon are a large vertical vorticity and Rossby number, a large vertical velocity, relatively flat horizontal spectra (contrary to the prevailing view of mesoscale dynamics), a large vertical buoyancy flux acting to restratify the upper ocean, a submesoscale energy conversion from potential to kinetic, a significant spatial and temporal intermittency in the upper ocean, and material exchanges between the surface boundary layer and pycnocline. Comparison with available observations indicates that submesoscale fronts and instabilities occur widely in the upper ocean, with characteristics similar to the simulations.

This is the second of three papers investigating the regime transition that occurs in numerical simulations for an idealized, equilibrium, subtropical, eastern boundary, upwelling current system similar to the California Current. The emergent upper-ocean submesoscale fronts are analyzed from phenomenological and dynamical perspectives, using a combination of composite averaging and separation of distinctive subregions of the flow. The initiating dynamical process for the transition is near-surface frontogenesis. The frontal behavior is similar to both observed meteorological surface fronts and solutions of the approximate dynamical model called surface dynamics (i.e., uniform interior potential vorticity q and diagnostic force balance) in the intensification of surface density gradients and secondary circulations in response to a mesoscale strain field. However, there are significant behavioral differences compared to the surfacedynamics model. Wind stress acts on fronts through nonlinear Ekman transport and creation and destruction of potential vorticity. The strain-induced frontogenesis is disrupted by vigorous submesoscale frontal instabilities that in turn lead to secondary frontogenesis events, submesoscale vortices, and excitation of even smaller-scale flows. Intermittent, submesoscale breakdown of geostrophic and gradient-wind force balance occurs during the intense frontogenesis and frontal-instability events.

The authors examine the turbulent properties of a baroclinically unstable oceanic flow using primitive equation (PE) simulations with high resolution (in both horizontal and vertical directions). Resulting dynamics in the surface layers involve large Rossby numbers and significant vortical asymmetries. Furthermore, the ageostrophic divergent motions associated with small-scale surface frontogenesis are shown to significantly alter the nonlinear transfers of kinetic energy and consequently the time evolution of the surface dynamics. Such impact of the ageostrophic motions explains the emergence of the significant cyclone-anticyclone asymmetry and of a strong restratification in the upper layers, which are not allowed by the quasigeostrophic (QG) or surface quasigeostrophic (SQG) theory. However, despite this strong ageostrophic character, some of the main surface properties are surprisingly still close to the surface quasigeostrophic equilibrium. They include a noticeable shallow (Ϸk Ϫ2) velocity spectrum as well as a conspicuous local spectral relationship between surface kinetic energy, sea surface height, and density variance over a large range of scales (from 400 to 4 km). Furthermore, surface velocities can be remarkably diagnosed from only the surface density using SQG relations. This suggests that the validity of some specific SQG relations extends to dynamical regimes with large Rossby numbers. The interior dynamics, on the other hand, strongly differ from the surface dynamics, involving a small Rossby number, a steep (Ϸk Ϫ4 ) velocity spectrum, and a somewhat steeper density spectrum. The compensation of the surface restratification by a destratification at depth confirms a connection between the surface and the interior induced by the small-scale divergent motions.

This is the last of a suite of three papers about the transition that occurs in numerical simulations for an idealized equilibrium, subtropical, eastern-boundary upwelling current system similar to the California Current. The transition is mainly explained by the emergence of ubiquitous submesoscale density fronts and ageostrophic circulations about them in the weakly stratified surface boundary layer. Here the highresolution simulations are further analyzed from the perspective of the kinetic energy (KE) spectrum shape and spectral energy fluxes in the mesoscale-to-submesoscale range in the upper ocean. For wavenumbers greater than the mesoscale energy peak, there is a submesoscale power-law regime in the spectrum with an exponent close to Ϫ2. In the KE balance an important conversion from potential to kinetic energy takes place at all wavenumbers in both mesoscale and submesoscale ranges; this conversion is the energetic counterpart of the vertical restratification flux and frontogenesis discussed in the earlier papers. A significant forward cascade of KE occurs in the submesoscale range en route to dissipation at even smaller scales. This is contrary to the inverse energy cascade of geostrophic turbulence and it is, in fact, fundamentally associated with the horizontally divergent (i.e., ageostrophic) velocity component. The submesoscale dynamical processes of frontogenesis, frontal instability, and breakdown of diagnostic force balance are all essential elements of the energy cycle of potential energy conversion and forward KE cascade.

[1] Eddy detection and tracking algorithms are applied to both satellite altimetry and a high-resolution (dx = 5 km) climatological model solution of the U.S. West Coast to study the properties of surface and undercurrent eddies in the California Current System. Eddy properties show remarkable similarity in space and time, and even somewhat in polarity. Summer and fall are the most active seasons for undercurrent eddy generation, while there is less seasonal variation at surface. Most of the eddies have radii in the range of 25-100 km, sea level anomaly amplitudes of 1-4 cm, and vorticity normalized by f amplitudes of 0.025-0.2. Many of the eddies formed near the coast travel considerable distance westward with speeds about 2 km/day, consistent with the b effect. Anticyclones and cyclones show equatorward and poleward displacements, respectively. Long-lived surface eddies show a cyclonic dominance. The subsurface California Undercurrent generates more long-lived anticyclones than cyclones through instabilities and topographic/coastline effects. In contrast, surface eddies and subsurface cyclones have much more widely distributed birth sites. The majority of the identified eddies have lifetimes less than a season. Eddies extend to 800-1500 m depth and have distinctive vertical structures for cyclones and anticyclones. Eddies show high nonlinearity (rotation speed higher than propagation speed) and hence can be efficient in transporting materials offshore.

International audienceThe Peru-Chile current System (PCS) is a region of persistent biases in global climate models. It has strong coastal upwelling, alongshore boundary currents, and mesoscale eddies. These oceanic phenomena provide essential heat transport to maintain a cool oceanic surface underneath the prevalent atmospheric stratus cloud deck, through a combination of mean circulation and eddy flux. We demonstrate these behaviors in a regional, quasi-equilibrium oceanic model that adequately resolves the mesoscale eddies with climatological forcing. The key result is that the atmospheric heating is large (>50 W m^-2) over a substantial strip >500 km wide off the coast of Peru, and the balancing lateral oceanic flux is much larger than provided by the offshore Ekman flux alone. The atmospheric heating is weaker and the coastally influenced strip is narrower off Chile, but again the Ekman flux is not sufficient for heat balance. The eddy contribution to the oceanic flux is substantial. Analysis of eddy properties shows strong surface temperature fronts and associated large vorticity, especially off Peru. Cyclonic eddies moderately dominate the surface layer, and anticyclonic eddies, originating from the nearshore poleward Peru-Chile Undercurrent (PCUC), dominate the subsurface, especially off Chile. The sensitivity of the PCS heat balance to equatorial intra-seasonal oscillations is found to be small. We demonstrate that forcing the regional model with a representative, coarse-resolution global reanalysis wind product has dramatic and deleterious consequences for the oceanic circulation and climate heat balance, the eddy heat flux in particular

[1] The spatial structure of nearshore wind, as measured by aircraft and analyzed in high-resolution atmospheric models (e.g., COAMPS), has strong influence on the patterns of upwelling circulation, surface temperature, and biogeochemical processes in coastal regions. Using a regional oceanic model for both the Southern California Bight and the Central California Coast, we demonstrate the nature of this upwelling sensitivity and infer that present wind analyses do not adequately determine the most important wind properties, viz., the strength of the nearshore curl and the speed drop-off near the coast.

International audienceThe oceanic general circulation is forced at large scales and is unstable to mesoscale eddies. Large-scale currents and eddy flows are approximately in geostrophic balance. Geostrophic dynamics is characterized by an inverse energy cascade except for dissipation near the boundaries. In this paper, we confront the dilemma of how the general circulation may achieve dynamical equilibrium in the presence of continuous large-scale forcing and the absence of boundary dissipation. We do this with a forced horizontal flow with spatially uniform rotation, vertical stratification and vertical shear in a horizontally periodic domain, i.e. a version of Eady's flow carried to turbulent equilibrium. A direct route to interior dissipation is presented that is essentially non-geostrophic in its dynamics, with significant submesoscale frontogenesis, frontal instability and breakdown, and forward kinetic energy cascade to dissipation. To support this conclusion, a series of simulations is made with both quasigeostrophic and Boussinesq models. The quasigeostrophic model is shown as increasingly inefficient in achieving equilibration through viscous dissipation at increasingly higher numerical resolution (hence Reynolds number), whereas the non-geostrophic Boussinesq model equilibrates with only weak dependence on resolution and Rossby number