Ocean Mixed Layer Research Articles

Large eddy simulation of the combined effect of heat fluxes and wave forcing of summer monsoon on a diurnal ocean mixed layer in the north Arabian Sea

Large eddy simulation of the combined effect of heat fluxes and wave forcing of summer monsoon on a diurnal ocean mixed layer in the north Arabian Sea

The Relative Importance of Ocean Advection and Surface Heat Fluxes During the Madden–Julian Oscillation in a Coupled Ocean–Atmosphere Model

AbstractIntraseasonal variability of ocean surface mixed layer temperature (MLT) in the tropics can be linked to the Madden–Julian Oscillation (MJO), the main source of intraseasonal atmospheric variability in the tropics. Here, we conduct a boreal winter mixed layer heat budget using 10‐day forecast composites of a coupled ocean–atmosphere Numerical Weather Prediction model of the UK Met Office to reveal that ocean advection plays a major role in modulating intraseasonal anomalies of MLT over the tropical Indian Ocean and the Maritime Continent. Large scale ( 10) intraseasonal anomalies of MLT ( 0.1 C) are driven by net surface heat flux. Ocean advection dominates at smaller horizontal scales ( 1), contributing up to 0.5 C to the intraseasonal MLT anomaly. Prior to the development of the enhanced MJO convection in the western Indian Ocean (phases 8 and 1), ocean advection reinforces the net heat flux warming in this region. However, ocean advection opposes the net heat flux warming in the Maritime Continent prior to the development of suppressed MJO convection in this region. When the MJO convection develops over the central Indian Ocean (phase 3), ocean advection (net surface heat flux) drives the intraseasonal MLT anomalies in the western Indian Ocean (central Indian Ocean and the Maritime Continent). Our results demonstrate the importance of ocean dynamics during the initiation and growth of the MJO, and raise implications for models that do not resolve these processes.

The Feedback of Ocean Cold Wakes on Three Successive Northward-Moving Tropical Cyclones in 2022

Abstract In 2022, three tropical cyclones (TCs) successively moved northward in the Northwest Pacific. Both the preceding and subsequent TCs can cause obvious sea surface temperature cooling (SSTC), which may last from a few days to a month. In this study, an uncoupled atmosphere-only model and a coupled ocean-atmosphere model were employed to investigate the effects of the SSTC caused by the preceding TCs on the track and intensity of the subsequent TCs. The results indicated that in the atmosphere-only simulations, the subsequent TC had a northeast shift in the track and a reduced intensity. The ocean-atmosphere coupled simulations showed similar TC track adjustments, but different TC intensity change results. The analysis suggested that the deviation of the TC track is the combined effect of the changes in the steering flow, beta drift, and diabatic heating. For the change of TC intensity, not only the SSTC but also the thickness of the oceanic mixed layer needs to be considered. The related physical mechanisms are analyzed in this paper.

Impact of internal wave drag on Arctic sea ice

Abstract A parameterization of the impact of internal waves on momentum transfer at the sea-ice–ocean interface based on previous work by McPhee has been implemented in a sea-ice model for the first time. The ice–ocean drag from internal waves is relevant for shallow mixed layer depth and the presence of a density jump at the pycnocline and is also a function of the strength of the stratification beneath the ocean mixed layer and geometry of the ice interface. We present results from a coupled sea-ice–ocean model where the parameterization of internal wave drag has been implemented. We conducted simulations spanning the years from 2000 to 2017. We find a deceleration of ice drift by 5–8% in both winter and summer, but with significant spatial and temporal variation reaching seasonal average values of ~10%. The spatial variation of ice transport leads to local impacts on deformed ice of magnitude ~0.05 m (2–5%), and reductions in ocean-to-ice heat fluxes of ~1 W m−2, and a decrease in bottom melt of ~0.02–0.04 cm d−1. There is an increase of up to 15% in thickness and ice concentration in the Canadian Arctic and a 10% overall impact on the total sea-ice volume.

A dynamically adaptive Langmuir turbulence parameterization scheme for variable wind wave conditions: Model application

Langmuir circulations and turbulence (LT) are crucial in the upper ocean mixed layer, significantly affecting the air-sea exchange of momentum, heat, and mass. The development of an appropriate LT parameterization scheme is vital for ocean modeling. This study employed the Large-eddy Simulation (LES) and the Physics-informed Neural Network (PINN) to optimize the KC04 Langmuir turbulence scheme by dynamically adjusting E6 as a key parameter determined by winds and waves. The LES simulations under different wind wave states indicated the PINN-inferred values for E6. Modelling results from GOTM in OCSPapa station demonstrated that the optimized scheme outperformed the original KC04 scheme in simulating the vertical eddy diffusivity and temperature, with an ∼6.24% annual reduction in the root mean square error (RMSE) for the temperature and an ∼8.23% reduction in the RMSE during autumn. Furthermore, the optimized scheme resulted in a thicker mixed layer, reaching 4.9 m. This enhanced LT parameterization scheme exhibited the improved robustness for variable spatiotemporal resolutions, significantly improving the modeling accuracy.

Responses of the Southern Ocean mixed layer depth to the eastern and central Pacific El Niño events during austral winter

Responses of the Southern Ocean mixed layer depth to the eastern and central Pacific El Niño events during austral winter

A Framework for Assessing Ocean Mixed Layer Depth Evolution

AbstractThe ocean surface mixed layer plays a crucial role as an entry or exit point for heat, salt, momentum, and nutrients from the surface to the deep ocean. In this study, we introduce a framework to assess the evolution of the mixed layer depth (MLD) for realistic forcings and preconditioning conditions. Our approach involves a physically‐based parameter space defined by three dimensionless numbers: λs representing the relative contribution of the buoyancy flux and the wind stress at the air‐sea interface, Rh the Richardson number which characterizes the stability of the water column relative to the wind shear, and f/Nh which characterizes the importance of the Earth's rotation (ratio of the Coriolis frequency f and the pycnocline stratification Nh). Four MLD evolution regimes (“restratification,” “stable,” “deepening,” and “strong deepening”) are defined based on the values of the normalized temporal evolution of the MLD. We evaluate the 3D parameter space in the context of 1D simulations and we find that considering only the two dimensions (λs, Rh) is the best choice of 2D projection of this 3D parameter space. We then demonstrate the utility of this two‐dimensional λs − Rh parameter space to compare 3D realistic ocean simulations: we discuss the impact of the horizontal resolution (1°, 1/12°, or 1/60°) and the Gent‐McWilliams parameterization on MLD evolution regimes. Finally, a proof of concept of using observational data as a truth indicates how the parameter space could be used for model calibration.

Prokaryotic morphological features and maintenance activities governed by seasonal productivity conditions.

Prokaryotic maintenance respiration and associated metabolic activities constitute a considerable proportion of the total respiration of carbon to CO2 in the ocean's mixed layer. However, seasonal influences on prokaryotic maintenance activities in terms of morphological and metabolic adaptations at low (winter) and high productivity (summer) are still unclear. To address this, we examined the natural prokaryotic communities at the mesocosm scale to analyse the differences in their morphological features and gene expression at low and high maintenance respiration, experimentally manipulated with the specific growth rate. Here, we showed that morphological features including membrane blebbing, membrane vesicles and cell‒cell connections occurred under high productivity. Metabolic adaptations associated with maintenance activities were observed under low productivity. Several Kyoto Encyclopedia of Genes and Genomes categories related to signal transduction, energy metabolism, and translational machinery supported maintenance activities under simulated winter conditions. Differential abundances of genes related to transporters, osmoregulation, nitrogen metabolism, ribosome biogenesis, and cold stress were observed. Our results demonstrate how specific growth rate in different seasons can influence resource allocation at the levels of morphological features and metabolic adaptations. This motivates further study of morphological features and their ecological role during high productivity, while investigations of metabolic adaptations during low productivity can advance our knowledge about maintenance activities.

Evaluation of Modeled Diurnal Warming Estimates for Application to Producing Sea Surface Temperature Analyses

AbstractAccurate knowledge of the amount of diurnal warming present in sea surface temperature (SST) observations at different times and effective depths is important for multiple applications including the production of blended SST analyses. This work explores the ability of a modified Kantha‐Clayson‐type one‐dimensional mixed layer ocean model with wave effects to accurately simulate the observed diurnal warming amplitude over a global grid when forced with coarse resolution numerical weather prediction (NWP) outputs. The sensitivity of the modeled diurnal amplitudes to multiple adjustable parameters and model configurations is evaluated to determine whether a preferred configuration can be identified that yields reliable predictions. The accuracy of the predictions is determined through comparison against estimates from operational SST retrievals from geostationary satellites. The results show that a single configuration can yield predictions that reproduce the observed range of diurnal warming amplitudes across a range of seasons and locations and an accurate occurrence frequency of the largest amplitude events. Simulated amplitudes fall along the one‐to‐one line with observations but with significant scatter due to factors including positioning of the NWP fluxes. The identified configuration is also shown to favorably reproduce diurnal warming observations from multiple research cruises. Overall uncertainty of the simulated diurnal warming amplitude across the different tests ranges between 0.4 K for all observations to ∼1 K for the largest warming events. While the focus is not on model comparisons, the results show improved performance relative to other models. Use of the model appears warranted but the associated uncertainty must be considered.

Conceptual Models for Exploring Sea-Surface Temperature Variability Vis-à Long-Range Weather Forecasting

This paper analyzes the ability of three conceptual stochastic models (one-box, two-box, and diffusion models) to reproduce essential features of sea surface temperature variability on intra-annual time scales. The variability of sea surface temperature, which is particularly influenced by feedback mechanisms in ocean surface–atmosphere coupling processes, is characterized by power spectral density, commonly used to analyze the response of dynamical systems to random forcing. The models are aimed at studying local effects of ocean–atmosphere interactions. Comparing observed and theoretical power spectra shows that in dynamically inactive ocean regions (e.g., north-eastern part of the Pacific Ocean), sea surface temperature variability can be described by linear stochastic models such as one-box and two-box models. In regions of the world ocean (e.g., north-western Pacific Ocean, subtropics of the North Atlantic, the Southern Ocean), in which the observed sea surface temperature spectra on the intra-annual time scales do not obey the ν−2 law (where ν is a regular frequency), the formation mechanisms of sea surface anomalies are mainly determined by ocean circulation rather than by local ocean–atmosphere interactions. The diffusion model can be used for simulating sea surface temperature anomalies in such areas of the global ocean. The models examined are not able to reproduce the variability of sea surface temperature over the entire frequency range for two primary reasons; first, because the object of study, the ocean surface mixed layer, changes during the year, and second, due to the difference in the physics of processes involved at different time scales.

Correction: Oceanic bottom mixed layer in the Clarion-Clipperton Zone: potential influence on deep-seabed mining plume dispersal

Correction: Oceanic bottom mixed layer in the Clarion-Clipperton Zone: potential influence on deep-seabed mining plume dispersal

Summer Westerly Wind Intensification Weakens Southern Ocean Seasonal Cycle Under Global Warming

AbstractSince the 1950s, observations and climate models show an amplification of sea surface temperature (SST) seasonal cycle in response to global warming over most of the global oceans except for the Southern Ocean (SO), however the cause remains poorly understood. In this study, we analyzed observations, ocean reanalysis, and a set of historical and abruptly quadrupled CO2 simulations from the Coupled Model Intercomparison Project Phase 6 archive and found that the weakened SST seasonal cycle over the SO could be mainly attributed to the intensification of summertime westerly winds. Under the historical warming, the intensification of summertime westerly winds over the SO effectively deepens ocean mixed layer and damps surface warming, but this effect is considerably weaker in winter, thus weakening the SST seasonal cycle. This wind‐driven mechanism is further supported by our targeted coupled model experiments with the wind intensification effects being removed.

Impact of ocean mixed layer depth on tropical cyclone characteristics: a numerical investigation

Impact of ocean mixed layer depth on tropical cyclone characteristics: a numerical investigation

Suitability of the CICE sea ice model for seasonal prediction and positive impact of CryoSat-2 ice thickness initialization

Abstract. The Los Alamos Community Ice CodE (CICE) sea ice model is being tested in standalone mode to identify biases that limit its suitability for seasonal prediction, where it is driven by atmospheric forcings from the NCEP Climate Forecast System Reanalysis (CFSR) and a built-in mixed-layer ocean model in CICE. The initial conditions for the sea ice and mixed-layer ocean are also from CFSR in the control experiments. The simulated sea ice extent agrees well with observations during the warm season at all lead times up to 12 months, in both the Arctic and the Antarctic. This suggests that CICE is able to provide useful sea ice edge information for seasonal prediction. However, the model's initial conditions have ice that is too thick in the Beaufort Sea, resulting in excessive ice extent in the Arctic at 6-month lead forecasts and errors in ice volume at all lead times when compared to available observations. To address this limitation, additional CS2_IC experiments were conducted, where the Arctic ice thickness was initialized using CryoSat-2 satellite observations while keeping all other initial fields the same as in the control experiments. This reduced the positive bias in the ice thickness in the initial conditions, leading to improvements in both the simulated ice edge and the ice thickness at the seasonal timescale. This indicates that CICE has the potential to improve its seasonal forecast skill and provide more accurate predictions of sea ice extent and thickness when initialized with a more realistic sea ice thickness. This study highlights that the suitability of CICE for seasonal prediction depends on various factors, including initial conditions such as sea ice thickness, in addition to sea ice coverage, as well as oceanic and atmospheric conditions.

Evaluating Existing Ocean Glider Sampling Strategies for Submesoscale Dynamics

Abstract Mixing in the upper ocean is important for biological production and the transfer of heat and carbon between the atmosphere and deep ocean, properties commonly targeted by observational campaigns using ocean gliders. We assess the reliability of ocean gliders to obtain a robust statistical representation of submesoscale variability in the ocean mixed layer of the Weddell Sea. A 1/48° regional simulation of the Southern Ocean is sampled with virtual “bow-tie” glider deployments, which are then compared against the reference model output. Sampling biases of lateral buoyancy gradients associated with the arbitrary alignment between glider paths and fronts are formally quantified, and the magnitude of the biases is comparable to observational estimates, with a mean error of 52%. The sampling bias leaves errors in the retrieved distribution of buoyancy gradients largely insensitive to deployment length and the deployment of additional gliders. Notable sensitivity to these choices emerges when the biases are removed by sampling perpendicular to fronts at all times. Detecting seasonal change in the magnitude of buoyancy gradients is sensitive to the glider-orientation sampling bias but the change in variance is not. We evaluate the impact of reducing the number of dives and climbs in an observational campaign and find that small reductions in the number of dive–climb pairs have a limited effect on the results. Lastly, examining the sensitivity of the sampling bias to path orientation indicates that the bias is not dependent on the direction of travel in our deep ocean study site. Significance Statement Recent observational campaigns have focused on using autonomous vehicles to better understand processes responsible for mixing in the surface region of the ocean. There exists uncertainty around how effective these missions are at returning reliable and representative information. This study seeks to quantify the performance of existing strategies in observing mixing processes, and we confirm that strategies are biased to underestimate indicators of mixing. Furthermore, compensating for the bias by increasing the number of resources or changing the manner in which resources are used has limited reward. Our findings are important for decision-making during the planning phase of an observational campaign and display that further innovations are required to account for the sampling bias.

Polar Lows and Their Effects on Sea Ice and the Upper Ocean in the Iceland, Greenland, and Labrador Seas

AbstractUsing two case studies, we analyze the effects of explicitly resolving polar lows in a global climate model (ICON‐Sapphire) with a high resolution of 2.5 km on the upper ocean and sea ice. We aim to understand the mechanism of how polar lows form in a global coupled model and how they interact with the upper ocean and sea ice. When polar lows form at the sea ice edge, they induce marine cold air outbreaks that lead to large heat loss from the ocean. This heat loss contributes to dense water formation in the Iceland and Greenland Seas, which replenishes the climatically important Denmark Strait Overflow Water (DSOW). The high wind speeds of polar lows open leads and polynyas in the sea ice cover, such as the Sirius Water Polynya in northeastern Greenland. Heat loss in polynyas is compensated for by the formation of new ice, and the rejected brine densifies the water on the Greenland shelf. In the Labrador Sea, polar lows intensify cold air outbreaks from the sea ice and rapidly deepen the ocean mixed layer. Resolving polar lows and kinematic features in the sea ice improves the realism of climate models, in particular the surface heat loss and the dense water formation in (sub)polar oceans.

Influence of Radiation Stress on Upper-Layer Ocean Temperature under Geostrophic Condition

Wave-induced radiation stress (RS), as a primary driver of ocean currents influenced by waves, plays an important role in the response of upper ocean temperatures under typhoons. Previous studies have mainly focused on wave-generated currents and coastal currents in nearshore areas. This paper incorporates the geostrophic effect into the wave-induced radiation stress of wave-current interaction, and the effect of waves on the changes in upper ocean temperature (including sea surface temperature (SST) and mixed layer temperature) under typhoon Nanmadol (2022) is studied. The FVCOM-SWAVE model is used to conduct a preliminary numerical study in the western Pacific Ocean. The RS with the geostrophic effect increased the horizontal and vertical components, leading to an enhancement in turbulent mixing and a decrease in SST by up to 1.0 °C to 1.4 °C, which is closer to the SST obtained by OISST remote sensing fusion observation data. In the strong divergence domain, the direction of the vortex flow exhibits a more pronounced turn to the right, accompanied by an increase in water velocity. The vertical temperature profile of the ocean shows that the water below is perturbed by the RS component of the geostrophic effect, and the depth of the mixed layer increases by about 2 m, which is closer to the depth of the mixed layer observed by the Argo floats, indirectly enhancing the vertical mass transport of the ocean. In general, this shows that RS, which takes into account geostrophic effects, enhances the effect of waves on the water below, indirectly leading to lower temperatures in the upper ocean, and the simulated results align more closely with the observed data, offering valuable insights for enhancing marine numerical forecasting accuracy.

The Inter‐Model Uncertainty of Projected Precipitation Change in Northern China: The Modulating Role of North Atlantic Sea Surface Temperature

AbstractPrecipitation changes in northern China are projected to increase in the Coupled Model Inter‐comparison Project Phase 6 (CMIP6) multi‐model ensemble. However, these projections are accompanied by notable inter‐model uncertainty, and the sources of this uncertainty remain largely unexplored. By analyzing 30 CMIP6 models, this research explores the source of inter‐model uncertainty in projected precipitation change and reveals the fundamental mechanism driving uncertainty spread. Following the empirical orthogonal function of inter‐model projected precipitation change, the leading mode displays a seesaw spatial pattern between northwest and north China. This phenomenon predominantly stems from the inter‐model divergence of projected sea surface temperature (SST) warming in the North Atlantic. Further scrutinizing the ocean mixed layer heat budget, we discover that the combined effect of surface sensible heat flux, net surface shortwave radiation flux, and ocean heat transport convergence influences heat flux and SST change of North Atlantic. The multi‐model projections indicate that localized increases in solar radiation and heat convergence warm sea surface, raising SST and initiating convective motion. This convective motion subsequently transforms the 200 hPa teleconnection wave train, leading to an anti‐phase pattern over northern China. This wave pattern modulates total cloud cover percentage, influences surface upward latent heat flux, and adjusts the top of atmosphere outgoing longwave radiation, collectively resulting in the seesaw pattern. Our study underscores the pivotal role of inter‐model disparities in North Atlantic SST warming projection, which is a primary driver of precipitation uncertainty in northern China. These insights offer an essential foundation for refining and diminishing inter‐model uncertainty.

Evaluation of multi-season convection-permitting atmosphere – mixed-layer ocean simulations of the Maritime Continent

Abstract. A multi-season convection-permitting regional climate simulation of the Maritime Continent (MC) using the Met Office Unified Model (MetUM) with 2.2 km grid spacing is presented and evaluated. The simulations pioneer the use of atmosphere–ocean coupling with the multi-column K profile parametrisation (KPP) mixed-layer ocean model in atmospheric convection-permitting climate simulations. Comparisons are made against a convection-parametrised simulation in which it is nested and which in turn derives boundary conditions from the ERA5 reanalysis. This paper describes the configuration, performance of the mean state and variability in the two simulations compared against observational datasets. The models have both minor sea surface temperature (SST) and wet precipitation biases. The diurnal cycle, representation of equatorial waves, and relationship between SST and precipitation are all improved in the convection-permitting model compared to the convection-parametrised model. The Madden–Julian oscillation (MJO) is present in both models with a faster-than-observed propagation speed. However, it is unclear whether fidelity of the MJO simulation is inherent to the model or whether it predominantly arises from the forcing at the boundaries.

Trend and interannual variability of summer marine heatwaves in the tropical Indian ocean: Patterns, mixed layer heat budget, and seasonal prediction

Marine heatwaves (MHWs) are extreme sea surface temperature (SST) events in all ocean basins, with far-reaching impacts on marine ecosystems and socio-economy. The leading patterns, trend, and interannual variability of summer MHWs in the tropical Indian Ocean (TIO) are investigated in this study. The first empirical orthogonal function (EOF) mode of frequency of MHWs exhibits a monopole pattern over the entire basin. This mode is highly associated with the concurrent Indian Ocean Basin warming, indicating a remarkable trend over the past four decades. The linear trend in MHW properties largely relates to increased summer Indian Ocean mean SST. The second EOF mode exhibits a zonal dipole with the MHW numbers increasing in the west and decreasing in the east. On the interannual timescale, the first two EOF modes are remotely affected by antecedent and concurrent El Niño events, respectively. The ocean mixed layer budget is utilized for examining the formation of different summer MHW patterns. During the preceding spring, the surface heat flux is important for the development of MHWs, while the ocean advections play a secondary role in the South Indian Ocean for the MHW monopole. Once the SST anomaly rises in summer, the ocean advections play a dominant role in maintaining the SST. Last, we assess the prediction skill of summer TIO MHWs by performing a bilinear seasonal statistical prediction model. Our results suggest the frequency of summer MHWs in the TIO could be predicted one season in advance. This study has great implications for understanding and predicting ocean extreme events in the TIO.