No abstract

Numerical experiments were performed using a three-dimensional large-eddy simulation model of the ocean surface mixed layer that includes the Craik-Leibovich vortex force [Craik 1977; Leibovich 1977] to parameterize the interaction of surface waves with mean currents. Results from the experiments show that the vortex force generates Langmuir circulations that can dominate vertical mixing. The simulated vertical velocity fields show linear, small-scale, coherent structures near the surface that extend downwind across the model domain. In the interior of the mixed layer, scales of motion increase to eddy sizes that are roughly equivalent to the mixed-layer depth. Cases with the vortex force have stronger circulations near the surface in contrast to cases with only heat flux and wind stress, particularly when the heat flux is positive. Calculations of the velocity variance and turbulence dissipation rates for cases with and without the vortex force, surface cooling, and wind stress indicate that wave-current interactions are a dominant mixing process in the upper mixed layer. Heat flux calculations show that the entrainment rate at the mixed-layer base can be up to two times greater when the vortex force is included. In a case with reduced wind stress, turbulence dissipation rates remained high near the surface because of the vortex force interaction with preexisting inertial currents. In deep mixed layers (-250 m) the simulations show that Langmuir circulations can vertically transport water 145 rn during conditions of surface heating. Observations of turbulence dissipation rates and the vertical temperature structure support the model results.

Submarine melting has been implicated as a driver of glacier retreat and sea level rise, but to date melting has been difficult to observe and quantify. As a result, melt rates have been estimated from parameterizations that are largely unconstrained by observations, particularly at the near‐vertical termini of tidewater glaciers. With standard coefficients, these melt parameterizations predict that ambient melting (the melt away from subglacial discharge outlets) is negligible compared to discharge‐driven melting for typical tidewater glaciers. Here, we present new data from LeConte Glacier, Alaska, that challenges this paradigm. Using autonomous kayaks, we observe ambient meltwater intrusions that are ubiquitous within 400 m of the terminus, and we provide the first characterization of their properties, structure, and distribution. Our results suggest that ambient melt rates are substantially higher (×100) than standard theory predicts and that ambient melting is a significant part of the total submarine melt flux. We explore modifications to the prevalent melt parameterization to provide a path forward for improved modeling of ocean‐glacier interactions.

Observations from a suite of platforms deployed in the coastal ocean are being combined with numerical models and simulations to investigate the processes that couple the atmosphere and ocean.

Cold pools dominate the surface temperature variability observed over the central Indian Ocean (0°, 80°E) for 2 months of research cruise observations in the Dynamics of the Madden–Julian Oscillation (DYNAMO) experiment in October–December 2011. Cold pool fronts are identified by a rapid drop of temperature. Air in cold pools is slightly drier than the boundary layer (BL). Consistent with previous studies, cold pools attain wet-bulb potential temperatures representative of saturated downdrafts originating from the lower midtroposphere. Wind and surface fluxes increase, and rain is most likely within the ~20-min cold pool front. Greatest integrated water vapor and liquid follow the front. Temperature and velocity fluctuations shorter than 6 min achieve 90% of the surface latent and sensible heat flux in cold pools. The temperature of the cold pools recovers in about 20 min, chiefly by mixing at the top of the shallow cold wake layer, rather than by surface flux. Analysis of conserved variables shows mean BL air is composed of 51% air entrained from the BL top (800 m), 22% saturated downdrafts, and 27% air at equilibrium with the ocean surface. The number of cold pools, and their contribution to the BL heat and moisture, nearly doubles in the convectively active phase compared to the suppressed phase of the Madden–Julian oscillation.

A simple quasi-equilibrium analytical model is used to explore hypotheses related to observed spatial correlations between sea surface temperatures and wind stress on horizontal scales of 50-500 km. It is argued that a plausible contributor to the observed correlations is the approximate linear relationship between the surface wind stress and stress boundary layer depth under conditions in which the stress boundary layer has come into approximate equilibrium with steady free-atmospheric forcing. Warmer sea surface temperature is associated with deeper boundary layers and stronger wind stress, while colder temperature is associated with shallower boundary layers and weaker wind stress. Two interpretations of a previous hypothesis involving the downward mixing of horizontal momentum are discussed, and it is argued that neither is appropriate for the warm-to-cold transition or quasi-equilibrium conditions, while one may be appropriate for the cold-to-warm transition. Solutions of a turbulent large-eddy simulation numerical model illustrate some of the processes represented in the analytical model. A dimensionless ratio ␥ A is introduced to measure the relative influence of lateral momentum advection and local surface stress on the boundary layer wind profile. It is argued that when ␥ A Ͻ 1, and under conditions in which the thermodynamically induced lateral pressure gradients are small, the boundary layer depth effect will dominate lateral advection and control the surface stress.

A numerical modelling study is presented focusing on the effects of mesoscale sea-surface temperature (SST) variability on surface fluxes and the marine atmospheric boundary-layer structure. A basic scenario is examined having two regions of SST anomaly with alternating warm/cold or cold/warm water regions. Conditions upstream from the anomaly region have SST values equal to the ambient atmosphere temperature, creating an upstream neutrally stratified boundary layer. Downstream from the anomaly region the SST is also set to the ambient atmosphere value. When the warm anomaly is upstream from the cold anomaly, the downstream boundary layer exhibits a more complex structure because of convective forcing and mixed layer deepening upstream from the cold anomaly. An internal boundary layer forms over the cold anomaly in this case, generating two distinct layers over the downstream region. When the cold anomaly is upstream from the warm anomaly, mixing over the warm anomaly quickly destroys the shallow cold layer, yielding a more uniform downstream boundary-layer vertical structure compared with the warm-tocold case. Analysis of the momentum budget indicates that turbulent momentum flux divergence dominates the velocity field tendency, with pressure forcing accounting for only about 20% of the changes in momentum. Parameterization of surface fluxes and boundary-layer structure at these scales would be very difficult because of their dependence on subgrid-scale SST spatial order. Simulations of similar flow over smaller scale fronts (<5 km) suggest that small-scale SST variability might be parameterized in mesoscale models by relating the effective heat flux to the strength of the SST variance.

The K-profile parameterization of upperocean mixing is tested and extended using observations and large eddy simulations of upper-ocean response to a westerly windburst. A nonlocal momentum flux term is added, and the amplitude of the nonlocal scalar flux is recalibrated. Parameterizations of Stokes drift effects are added following recent work by McWilliams and Sullivan (2001). These changes allow the parameterization to produce both realistic gradients of momentum and scalars in the nocturnal boundary layer and enhanced mixing during stable conditions. The revised parameterization is expected to produce improved representations of lateral advection and sea-surface temperature in large-scale models.