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.

In this study, uncoupled and coupled ocean-atmosphere simulations are carried out over the 5 California Upwelling System to assess the dynamic ocean-atmosphere interactions, viz.,the 6 ocean surface current feedback to the atmosphere. We show the current feedback, by mod-7 ulating the energy transfer from the atmosphere to the ocean, controls the oceanic Eddy 8 Kinetic Energy (EKE), and for the first time, we demonstrate the current feedback has an 9 opposite effect on the surface stress and on the wind itself. The current feedback acts as 10 an oceanic eddy killer, reducing by half the Surface EKE, and by 27% the depth-integrated