Search

The mixing of water in the ocean is a key influence on our climate and ecosystems, but its importance is often under-recognized. Mixing in the seas around Antarctica—a key process that redistributes heat, carbon, nutrients, plankton, and all other things in the sea, with profound consequences—also affects the stability of the continent's glaciers and ice sheets, with consequences for sea level rise globally. A recent discovery showed that when the fronts of glaciers disintegrate, they "calve" huge chunks of ice that can cause underwater tsunamis in the ocean, which can spread for many miles and cause strong bursts of mixing when they break. This surprising finding is changing the way we think about mixing close to Antarctica, what causes it, and how it matters. [ABSTRACT FROM AUTHOR]

Shear instability is a dominant mechanism for mixing in the stratified oceans and coastal seas. For the first time, we present fine‐scale, direct measurements of shear instabilities in the bottom front generated by the Kuroshio intrusion into the Changjiang (Yangtze) river plume. Shear instabilities were identified using a shipboard echo‐sounder and the resulting turbulent mixing was quantified using a turbulence microstructure profiler. The shear instabilities generate vigorous turbulent mixing with dissipation rate and vertical diffusivity up to O (10−4 m2 s−3) and O (10−1 m2 s−1), respectively, comparable to values associated with shear instabilities observed in river plumes and western boundary currents but several orders of magnitude larger than typical values in the open ocean. The enhanced turbulence may contribute significantly to mixing between the Kuroshio water and coastal water and thereby alter the coastal biogeochemistry cycles. Plain Language Summary: Strong velocity shear across a density interface can produce instability, causing the interface to roll up to form a train of vortices ("billows") that subsequently break down into turbulence. This process is important in the vertical mixing of the oceanic interior and coastal seas. Using acoustic images, we find that this flow instability is induced by the interaction between the Kuroshio current and the Yangtze river plume. We also quantify the resulting turbulent mixing using high‐resolution measurements. The turbulence efficiently mixes the Kuroshio and river plume‐influence water and can change biogeochemistry cycles in the coastal seas. Key Points: Interaction between Kuroshio and river plume‐influenced coastal waters generates shear instabilityShear instability elevates turbulence level by 2–3 orders of magnitude compared to the fluid above and belowThe enhanced turbulence effectively mixes oceanic and coastal waters [ABSTRACT FROM AUTHOR]

We evaluate the vertical turbulent-kinetic-energy (TKE) mixing scheme of the NEMO-SI3 ocean–sea-ice model in sea-ice-covered regions of the Arctic Ocean. Specifically, we assess the parameters involved in TKE mixed-layer-penetration (MLP) parameterization. This ad hoc parameterization aims to capture processes that impact the ocean surface boundary layer, such as near-inertial oscillations, ocean swells, and waves, which are often not well represented in the default TKE scheme. We evaluate this parameterization for the first time in three regions of the Arctic Ocean: the Makarov, Eurasian, and Canada basins. We demonstrate the strong effect of the scaling parameter that accounts for the presence of sea ice. Our results confirm that TKE MLP must be scaled down below sea ice to avoid unrealistically deep mixed layers. The other parameters evaluated are the percentage of energy penetrating below the mixed layer and the length scale of its decay with depth. All these parameters affect mixed-layer depth and its seasonal cycle, surface temperature, and salinity, as well as underlying stratification. Shallow mixed layers are associated with stronger stratification and fresh surface anomalies, and deeper mixed layers correspond to weaker stratification and salty surface anomalies. Notably, we observe significant impacts on sea-ice thickness across the Arctic Ocean in two scenarios: when the scaling parameter due to sea ice is absent and when the TKE mixed-layer-penetration process vanishes. In the former case, we observe an increase of several meters in mixed-layer depth, along with a reduction in sea-ice thickness ranging from 30 to 40 cm, reflecting the impact of stronger mixing. Conversely, in the latter case, we notice that a shallower mixed layer is accompanied by a moderate increase in sea-ice thickness, ranging from 10 to 20 cm, as expected from weaker mixing. Additionally, interannual variability suggests that experiments incorporating a scaling parameter based on sea-ice concentration display an increased mixed-layer depth during periods of reduced sea ice, which is consistent with observed trends. These findings underscore the influence of enhanced ocean mixing, through specific parameterizations, on the physical properties of the upper ocean and sea ice. [ABSTRACT FROM AUTHOR]

The relationships between metallogenic capacity and geochemical features of adakitic rocks along the South Tan-Lu Fault (STLF) remain unclear. In this study, the ore-barren adakitic rocks (Tuncang, Guandian and Wawuliu) exhibit higher K2O/Na2O ratios and lower Sr/Y ratios than the ore-bearing adakitic rocks (Chuzhou and Shangyaopu). These differences strongly suggest that the ore-barren adakitic rocks originated from the thickened lower continental crust (LCC), whereas the ore-bearing adakitic rocks were derived from oceanic slabs. Notably, the Tuncang granite exhibits higher Y/Yb and (Ho/Yb)N ratios than the Guandian granodiorite and Wawuliu intrusion. Accordingly, we suggest that the Tuncang granite likely originated from a delaminated eclogitic LCC, whereas the Guandian and Wawuliu intrusions were derived from a thickened basaltic LCC sources. The occurrence of diorite and gabbro mafic microgranular enclaves (MMEs) within the Tuncang granite strongly suggests a magma-mixing process. Considering their MgO contents and εNd(t) and (87Sr/86Sr)i values, we suggest that the gabbro MMEs were likely derived from an enriched mantle source previously metasomatized by subduction-related components and that the diorite MMEs were subsequently formed by magma mixing. Due to the slightly younger ages of the ore-bearing adakitic rocks, we propose a model in which the ore-barren adakitic rocks formed through LCC delamination at 130 Ma and the ore-bearing adakitic rocks formed through oceanic slab remelting at 125 Ma. Consequently, the exploration of Cu–Au mineralization along the STLF should focus on younger oceanic slab-derived adakitic rocks. [ABSTRACT FROM AUTHOR]

Deficiencies in upper ocean vertical mixing parameterizations contribute to tropical upper ocean biases in global coupled general circulation models, affecting their simulated ocean heat uptake and ENSO variability. To better understand these deficiencies, we develop a suite of ocean model experiments including both idealized single column models and realistic global simulations. The vertical mixing parameterizations are first evaluated using large eddy simulations as a baseline to assess uncertainties and evaluate their implied turbulent mixing. Global models are then developed following NOAA/GFDL's 0.25° nominal horizontal grid spacing OM4 (uncoupled) configuration of the MOM6 ocean model, with various modifications that target biases in the original model. We test several enhancements to the existing mixing schemes and evaluate them against observational constraints from Tropical Atmosphere Ocean moorings and Argo floats. In particular, we find that we can improve the diurnal variability of mixing in OM4 via modifications to its surface boundary layer mixing scheme, and can improve the net mixing in the upper thermocline by reducing the background vertical viscosity, allowing for more realistic, less diffuse currents. The improved OM4 model better represents the mixing, leading to improved diurnal deep‐cycle variability, a more realistic time‐mean tropical thermocline structure, and a better Pacific Equatorial Undercurrent. Plain Language Summary: Computational models of the oceanic and atmospheric circulation are critical tools for understanding and projecting changes in the Earth's climate. These models have errors that can arise from many sources, including model formulation or the choices in applying the model. One of the more well known sources of error is the representation of turbulent mixing. In this work we consider specially designed small‐scale models that simulate turbulent mixing, and use their results to improve the representation of turbulence and its induced mixing in large‐scale models. In particular, we investigate how the intensity of mixing varies throughout the day, considering the progression from deep mixing during cooler nighttime surface conditions to shallower mixing in the presence of strong solar heating during the day. We find some modifications to the mixing scheme in the ocean climate model that can improve the model solutions when compared to the real ocean. Key Points: Large eddy simulation results are used to evaluate the diurnal cycle of equatorial turbulent mixing in the OM4 ocean modelReducing vertical viscosity in an ocean model increases shear near the Equatorial Undercurrent (EUC) and can result in increased vertical mixingVertical grid spacing of a few meters helps to resolve shear mixing events within and below the EUC in ocean models [ABSTRACT FROM AUTHOR]

The 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. Plain Language Summary: Vertical mixing of water near the ocean surface occurs when cold air temperatures create dense cold water at the surface that tends to sink in the ocean or when a strong wind induces turbulence at the ocean surface. These processes mix heat and salt and create a layer at the top of the ocean that has a uniform temperature and salinity and that is called the "mixed layer." This mixed layer plays a fundamental role in the Earth climate system, and the representation of its evolution in ocean models hence needs to be assessed. For this purpose, we propose to map the mixed layer evolution in a three‐dimensional space where the first axis is related to the wind and the surface heat flux, the second axis to the stability of the water column, and the third axis to the Earth's rotation. We show that this tool performs statistically well and we present how to use it in the context of realistic ocean models. Key Points: A parameter space is proposed to assess the evolution of the mixed layer depth for realistic forcings and preconditioning conditionsAn evaluation of a collection of 1D simulations shows a statistically good performance of the parameter spaceTwo applications demonstrate the utility of the parameter space for assessing and comparing realistic 3D simulations [ABSTRACT FROM AUTHOR]

Motivated by the importance of mixing arising from dissipating internal waves (IWs), vertical profiles of internal‐wave dissipation from a high‐resolution regional ocean model are compared with finestructure estimates made from observations. A horizontal viscosity scheme restricted to only act on horizontally rotational modes (such as eddies) is introduced and tested. At lower resolutions with horizontal grid spacings of 2 km, the modeled IW dissipation from numerical model agrees reasonably well with observations in some cases when the restricted form of horizontal viscosity is used. This suggests the possibility that if restricted forms of horizontal viscosity are adopted by global models with similar resolutions, they could be used to diagnose and map IW dissipation distributions. At higher resolutions with horizontal grid spacings of ∼250 m, the dissipation from vertical shear and horizontal viscosity act much more strongly resulting in dissipation overestimates; however, the vertical‐shear dissipation itself is found to agree well with observations. Plain Language Summary: Oceanic mixing impacts circulation, stratification (layering by density), and the uptake and transport of heat and nutrients. Over most of the ocean, mixing is caused by the breaking (turnover) of internal waves lying on the interfaces of density layers. Most ocean models do not contain a resolved internal wavefield, and therefore must parameterize internal wave (IW) mixing based upon external information. Recently developed high‐resolution ocean models with credible representations of internal waves may make it possible to map and understand global IW mixing without use of external information. Here we compare vertical profiles (profiles in depth) of IW dissipation in a regional model, which can be used to understand sensitivities to numerical schemes and grid spacings. With grid spacings that are attainable in global models, modeled dissipation profiles lie somewhat close to observed profiles, as long as certain choices are made within the numerical schemes. One numerical dissipation scheme is designed to realistically remove energy from eddy fields, which are the non‐wavelike motions in the ocean, and we have adapted this scheme to act less strongly on internal waves. Using this modified scheme, we find that high‐resolution global models may already be able to map IW dissipation. Key Points: Vertical profiles of internal wave (IW) dissipation in high‐resolution regional ocean simulations are compared with observed profilesProfiles in runs with a restricted form of horizontal viscosity and resolutions attainable in global models are close to observationsResults suggest high‐resolution global models can be used to map IW dissipation after numerical sensitivities are tested in regional models [ABSTRACT FROM AUTHOR]

The background potential energy (BPE) is the only reservoir that double diffusive instabilities can tap their energy from when developing from an unforced motionless state with no available potential energy (APE). Recently, Middleton and Taylor linked the extraction of BPE into APE to the sign of the diapycnal component of the buoyancy flux, but their criterion can predict only diffusive convection instability, not salt finger instability. Here, we show that the problem can be corrected if the sign of the APE dissipation rate is used instead, making it emerge as the most fundamental criterion for double diffusive instabilities. A theory for the APE dissipation rate for a two-component fluid relative to its single-component counterpart is developed as a function of three parameters: the diffusivity ratio, the density ratio, and a spiciness parameter. The theory correctly predicts the occurrence of both salt finger and diffusive convection instabilities in the laminar unforced regime, while more generally predicting that the APE dissipation rate for a two-component fluid can be enhanced, suppressed, or even have the opposite sign compared to that for a single-component fluid, with important implications for the study of ocean mixing. Because negative APE dissipation can also occur in stably stratified single-component and doubly stable two-component stratified fluids, we speculate that only the thermodynamic theory of exergy can explain its physics; however, this necessitates accepting that APE dissipation is a conversion between APE and the internal energy component of BPE, in contrast to prevailing assumptions. [ABSTRACT FROM AUTHOR]

Oceanic biogeochemistry plays a pivotal role in regulating Earth's climate system by governing the cycling of key elements such as carbon, oxygen, and nutrients. Various metocean processes including wind, tides, currents, waves, and eddies significantly influence the dynamics of this system. In particular, ocean surface waves contribute to this intricate interplay by facilitating the exchange of heat, gas, and momentum between the atmosphere and the ocean. Although wave‐coupled effects are substantial, studies on their impacts on oceanic biogeochemistry, particularly on phytoplankton abundance are missing in present‐day research. Additionally, wave‐coupled effects cannot be disregarded in regions like the Southern Ocean (SO), where wind and waves activities are prominent. Addressing this gap, we incorporated a parameterization of surface wave mixing into a global ocean biogeochemical model to investigate its effects on upper ocean and biogeochemical parameters. Our results show that surface wave mixing has significant impacts on sea surface temperature (SST), mixed layer depth (MLD), and nutrient distribution—key factors that influence phytoplankton growth. Additionally, we observed significant improvements in model biases against the observations. During austral summer, additional mixing from surface waves can significantly lower SST by 0.5°C, deepen MLD by 13 m, and enhance Chlorophyll‐a (Chl‐a) concentration, an index of phytoplankton population, by 8% in the SO. This observed increase in Chl‐a concentration is mainly driven by enhanced dissolved iron levels resulting from wave‐induced mixing. Our findings underscore the significance of incorporating surface wave mixing in ocean biogeochemistry studies, an aspect that is currently overlooked. Plain Language Summary: Ocean biogeochemistry encompasses the complex interactions of the ocean's biological, geological, and chemical processes and affects the Earth's climate system. This system is highly impacted by various physical forcings such as wind, currents, waves, tides, and eddies. Among these forcings, surface waves play a crucial role by exchanging heat, gas, and momentum between the atmosphere and ocean. Understanding the roles of surface waves on upper ocean mixing and oceanic biogeochemistry is paramount, especially in regions like the Southern Ocean (SO), where wind and waves are prominent. However, present‐day studies lack research on wave‐coupled effects on ocean biogeochemistry, particularly in the SO. Addressing this research gap, we investigated the influence of surface wave mixing on oceanic biogeochemistry by adding a parameterization of surface wave mixing in a global ocean biogeochemical model. From the investigation, we found significant effects of wave coupling on upper ocean and biogeochemical parameters, such as cooling of sea surface temperature, deepening mixed layer depth, enhancing nutrient levels, and consequently, increasing phytoplankton distribution during austral summer in the SO. This study contributes to a deeper understanding of the complex interplay between ocean surface waves and oceanic biogeochemistry, especially in the SO. Key Points: Present‐day studies lack research on the wave‐coupled effects on ocean biogeochemistry, particularly in the Southern OceanA surface wave mixing parameterization is used to investigate the wave‐coupled effects on the upper ocean and biogeochemical parametersWave coupling influences ocean temperature, mixed layer depth, and nutrient levels which are essential for phytoplankton growth [ABSTRACT FROM AUTHOR]

A holistic review is given of the Southern Ocean dynamic system, in the context of the crucial role it plays in the global climate and the profound changes it is experiencing. The review focuses on connections between different components of the Southern Ocean dynamic system, drawing together contemporary perspectives from different research communities, with the objective of closing loops in our understanding of the complex network of feedbacks in the overall system. The review is targeted at researchers in Southern Ocean physical science with the ambition of broadening their knowledge beyond their specific field, and aims at facilitating better‐informed interdisciplinary collaborations. For the purposes of this review, the Southern Ocean dynamic system is divided into four main components: large‐scale circulation; cryosphere; turbulence; and gravity waves. Overviews are given of the key dynamical phenomena for each component, before describing the linkages between the components. The reviews are complemented by an overview of observed Southern Ocean trends and future climate projections. Priority research areas are identified to close remaining loops in our understanding of the Southern Ocean system. Plain Language Summary: The United Nations has identified 2021–2030 as the Decade of Ocean Science, with a goal to improve predictions of ocean and climate change. Improved understanding of the Southern Ocean is crucial to this effort, as it is the central hub of the global ocean. The Southern Ocean is the formation site for much of the dense water that fills the deep ocean, sequesters the majority of anthropogenic heat and carbon, and controls the flux of heat to Antarctica. The large‐scale circulation of the Southern Ocean is strongly influenced by interactions with sea ice and ice shelves, and is mediated by smaller scale processes, including eddies, waves, and mixing. The complex interplay between these dynamic processes remains poorly understood, limiting our ability to understand, model and predict changes to the Southern Ocean, global climate and sea level. This article provides a holistic review of Southern Ocean processes, connecting the smallest scales of ocean mixing to the global circulation and climate. It seeks to develop a common language and knowledge‐base across the Southern Ocean physical science community to facilitate knowledge‐sharing and collaboration, with the aim of closing loops on our understanding of one of the world's most dynamic regions. Key Points: Contemporary perspectives are reviewed on the different components of the Southern Ocean dynamic system from distinct research communitiesKey connections between different components of Southern Ocean dynamics are highlightedCross‐cutting priorities for future Southern Ocean physical science are identified [ABSTRACT FROM AUTHOR]

This study examines the impact of wave-induced processes (WIPs) in modulating thermosteric sea-level changes, highlighting the need to include these processes in future sea-level rise assessments and climate projections. The impact of wave-induced processes on thermosteric sea-level changes is investigated using coupled ocean-wave simulations. These simulations include the effects of Stokes-Coriolis forcing, sea-state dependent surface stress and energy fluxes, and wave-induced mixing. The experiments use a high-resolution configuration of the Geesthacht COAstal Model SysTem (GCOAST), covering the Northeast Atlantic, the North Sea and the Baltic Sea. The GCOAST system uses the Nucleus for European Modelling of the Ocean (NEMO) ocean model to account for wave-ocean interactions and ocean circulation. It is fully coupled with the WAM spectral wind wave model. The aim is to accurately quantify the sea state contribution to thermosteric sea level variability and trends over a 26-year period (1992–2017). The ability of wave-ocean coupled simulations to reveal the contribution of sea state to sea level variability and surge is demonstrated. It is clear that wave-induced processes (WIPs) play a significant role in sea surface dynamics, ocean mixing (mixed layer thickness) and modulation of air-sea fluxes (e.g. heat flux) in both winter (10–20%) and summer (10%), which in turn affect thermosteric sea level variability. The North Atlantic (in summer) and the Norwegian Trench (in winter) show significant contributions (40%) to the thermosteric sea-level variability due to wave-induced processes. The influence of WIPs on thermosteric sea level trends in the North Atlantic is up to the order of 1 mm yr-1 in both winter and summer, in the open ocean and at the shelf break. Smaller contributions are observed over the shelf areas of the North Sea. This study underscores the crucial role of WIPs in modulating sea-level changes and highlights the importance of including these processes in future sea-level rise assessments and climate projections. [ABSTRACT FROM AUTHOR]

Wave-induced mixing can enhance vertical mixing in the upper ocean, facilitating the exchange of heat and momentum between the surface and deeper layers, thereby influencing ocean circulation and climate patterns. Building on previous research, this study proposes a wave-induced mixing parameterization scheme (referred to as EXP3) specifically designed for typhoon periods. This scheme was integrated into the fully coupled ocean–wave–atmosphere model COAWST and applied to analyze Typhoon In-Fa (2021) as a case study. The simulation results were validated against publicly available data, demonstrating a good overall match with observed phenomena. Subsequently, a comparative analysis was conducted between the EXP3 scheme, the previous scheme (EXP2) and the original model scheme (EXP1). Validation against Argo and Drifter buoy data revealed that both EXP2 and EXP3, which include wave-induced mixing effects, resulted in a decrease in the simulated mixed layer depth (MLD) and mixed layer temperature (MLT), with EXP3 showing closer alignment with the observed data. Compared to the other two experiments, EXP3 enhanced vertical motion in the ocean due to intensified wave-induced mixing, leading to increased upper-layer water divergence and upwelling, a decrease in sea surface temperature and accelerated rightward deflection of surface currents. This phenomenon not only altered the temperature structure of the ocean surface layer but also significantly impacted the regional ocean dynamics. [ABSTRACT FROM AUTHOR]

The dissolved nitrate is one of the major essential nutrients for primary production in the tropical ocean and it is brought to the surface though mixing. The depth of nitracline determines how much of nitrate enters to the upper ocean through mixing. The depth of nitracline is traditionally estimated using nitrate concentrations measured at standard depths that introduces significant error due to interpolation of data. Based nitrate profiles measured at 5 m interval using nitrate sensors onboard Argo, the exact depth of nitracline was derived in the Bay of Bengal that displayed a significant linear relationship with depth of 26°C isotherm (D26). Based on climatological D26, the temporal and spatial variations in the depth of nitracline was estimated for the entire Bay of Bengal. The depth of nitracline varied between 5 and 80 m with large spatial and temporal variability in the Bay of Bengal and it is 5–20 m deeper than simulations of numerical models. The relationship between the depth of nitracline and photic zone integrated primary production indicates that 7.5 ± 3 mgC m−2 d−1 of primary production increases due to shallowing of 1 m of depth of nitracline. Therefore, models seem to be over estimating the photic zone integrated primary production by 5%–25% in the Bay of Bengal. The numerical models may improve the simulation of primary production and carbon cycling by accounting the accurate estimation of depth of nitracline in the model initialization. Plain Language Summary: In the tropical ocean, nitrate mainly controls the primary production and identifying the depth of nitracline, where nitrate concentration was 1 μM from that of surface, is important in simulating primary production by numerical models. All the numerical models are using nutrietns profiles collected at the standard depth and estimated depth of nitracline. Since the data collected at standard depth, and interpolation between depth may result in large error in the computation of nitracline depth. We used high depth resolution nitrate data collected by sensors and estimated nitracline depth and the same compared with modeled data. It was observed that models are estimating shallow nitracline than actual as a result higher biological response is simulated than reality. This exercise would be useful in improving model simulation in future. Key Points: The depth of nitracline is strongly controlled by depth of 20°C isothermDeeper nitracline depth was noticed during spring and shallowers during summerModels estimate shallower nitracline depth by 5–20 m than measured [ABSTRACT FROM AUTHOR]

Near-inertial waves (NIWs), a special form of internal waves with a frequency close to the local Coriolis frequency, are ubiquitous in the ocean. NIWs play a crucial role in ocean mixing, influencing energy transport, climate change, and biogeochemistry. This manuscript briefly reviews the generation and propagation of NIWS in the oceans. NIWs are primarily generated at the surface by wind forcing or through the water column by nonlinear wave-wave interaction. Especially at critical latitudes where the tidal frequency is equal to twice the local inertial frequency, NIWs can be generated by a specific subclass of triadic resonance, parametric subharmonic instability (PSI). There are also other mechanisms, including lee wave and spontaneous generation. NIWs can propagate horizontally for hundreds of kilometers from their generating region and radiate energy far away from their origin. NIWs also penetrate deep into the ocean, affecting nutrient and oxygen redistribution through altering mixing. NIW propagation is influenced by factors such as mesoscale eddies, background flow, and topography. This review also discussed some recent observational evidence of interactions between NIWs from different origins, suggesting a complicated nonlinear interaction and energy cascading. Despite the long research history, there are still many areas of NIWs that are not well defined. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Internal tides are pivotal dynamic processes enhancing the mixing of oceanic waters and facilitating energy transfer across various scales within the ocean. In recent years, the proliferation of satellite altimetry observations has enabled global predictions of the elevation and phase of internal tides. This study, leveraging the advanced global internal tide prediction model known as the Multivariate Inversion of Ocean Surface Topography-Internal Tide Model (MIOST-IT), employs a two-dimensional Fourier bandpass filtering approach to decompose the internal tides in the Luzon Strait, thereby addressing the east–west directional blind zones inherent in along-track satellite altimetry-based modal decomposition. To further elucidate the propagation trajectories of individual tidal modes in different directions, we introduce the directional Fourier filter method to characterize the spatial distribution features of each modal internal tide in the vicinity of the Luzon Strait. This work significantly enhances the accuracy and reliability of extracting parameters for distinct modal internal tides, furnishing a scientific basis for subsequent studies on internal tide dynamics and model refinement. [ABSTRACT FROM AUTHOR]

The first operational version of the coupled Hurricane Analysis and Forecast System (HAFSv1) launched in 2023 consists of the HYbrid Coordinate Ocean Model (HYCOM) and finite-volume cubed-sphere (FV3) dynamic atmosphere model. This system is a product of efforts involving improvements and updates over a 4-year period (2019-2022) through extensive collaborations between the Environmental Modeling Center at the US National Centers for Environmental Prediction (NCEP) and NOAA Atlantic Oceanography and Meteorology Laboratory. To provide two sets of numerical guidance, the initial operational capability of HAFSv1 was configured to two systems -HFSA and HFSB. In this study, we present in-depth analysis of the forecast skills of the upper ocean that was co-evolved by the HFSA and HFSB. We chose hurricane Laura (2020) as an example to demonstrate the interactions between the storm and oceanic mesoscale features. Comparisons performed with the available in situ observations from gliders as well as Argos and National Data Buoy Center moorings show that the HYCOM simulations have better agreement for weak winds than high winds (greater than Category 2). The skill metrics indicate that the model sea-surface temperature (SST) and mixed layer depth (MLD) have a relatively low correlation. The SST, MLD, mixed layer temperature (MLT), and ocean heat content (OHC) are negatively biased. For high winds, SST and MLT are more negative, while MLD is closer to the observations with improvements of about 8%-19%. The OHC discrepancy is proportional to predicted wind intensity. Contrarily, the mixed layer salinity (MLS) uncertainties are smaller and positive for higher winds, probably owing to the higher MLD. The less-negative bias of MLD for high winds implies that the wind-force mixing is less effective owing to the higher MLD and high buoyancy stability (approx. 1.5-1.7 times) than the observations. The heat budget analysis suggests that the maximum heat loss by hurricane Laura was 0(< 3°C per day). The main contributor here is advection, followed by entrainment, which act against or with each other depending on the storm quadrant. We also found relatively large unaccountable heat residuals for the in-storm period, and the residuals notably led the heat tendency, meaning that further improvements of the subscale simulations are warranted. In summary, HYCOM simulations showed no systematic differences forced by either HFSA or HFSB. [ABSTRACT FROM AUTHOR]

The atmosphere‐ocean is a highly coupled system with significant diurnal and hourly variations. However, current coupled models usually lack sub‐diurnal scale processes at the air‐sea interface due to the finite vertical resolution for ocean discretization. Previous modeling studies showed that sub‐diurnal scale air‐sea interaction processes are important for ocean mixing. Here, by designing an integrated sub‐diurnal parameterization (ISDP) scheme which combines different temperature profiling functions, we stress sub‐diurnal air‐sea interactions to better represent the local ocean mixing. This scheme has been implemented into two coupled models which contributed to the Climate Model Intercomparison Project (CMIP), referenced by the Intergovernmental Panel on Climate Change—Community Earth System Model and Coupled Model version 2. The results show that the ISDP scheme improves model simulations with better climatology and more realistic spectra, especially in the tropics and North Pacific Ocean. With the scheme, the tropical cold tongue bias is significantly relaxed by reducing the overestimation of ocean upper mixing, and the cold bias of North Pacific Ocean is reduced due to the improvement on currents and net heat fluxes. Our scheme may help better the simulation and prediction skills of coupled models when their horizontal resolution becomes fine but vertical resolution remains relatively coarse as it describes high‐frequency air‐sea interactions more realistically. Plain Language Summary: The atmosphere and ocean interact with each other. In these interactions, changes may occur over different time periods. Some changes take years or months, while others happen in days or hours. The atmosphere changes quickly throughout the day, such as air temperature, wind, and rain. These quick changes in the atmosphere can affect the ocean rapidly. Similarly, there are quick changes in the ocean. These quick changes also affect the atmosphere in return. Understanding these quick changes is important. However, it is difficult to perfectly capture the quick changes in the climate models because they do not always represent the ocean behavior accurately. In this study, we developed a new scheme (named the sub‐diurnal scale parameterization, ISDP) to better represent these quick changes. We added this new scheme into two widely‐used climate models. These models are important tools for studying climate change. Our new scheme represents better the ocean mixing and interaction based on the local weather conditions. When we use ISDP, the climate models get better at simulating climate. They're more accurate in showing ocean temperatures in the tropics. These models make the tropical oceans seem colder than they really are, but our scheme, to some extent, fixes this problem. The results also show the more accurate ocean temperatures in the North Pacific Ocean. Our new scheme has a great potential to make climate models more accurate. Key Points: An integrated sub‐diurnal scale parameterization (ISDP) scheme has a more appropriate representation for local ocean mixingThe ISDP scheme has been implemented into two models which contributed to the Climate Model Intercomparison Project, referenced by IPCC: Community Earth System Model and Coupled Model version 2The model tropical cold tongue bias with ISDP is relaxed by reducing the overestimation of ocean upper mixing [ABSTRACT FROM AUTHOR]

Recent studies have challenged the validity of the Redfield ratio. It is proposed that physical and biogeochemical processes govern the geographical variations in carbon : nitrogen : phosphorus (C : N : P) ratios. However, this proposal remains to be examined through concurrent observations of C : N : P ratios with physical and biogeochemical processes in various marine reservoirs. Here, we sampled the Arabian Sea for its C, N, and P content in organic and inorganic pools during the winter monsoon. We analyzed the role of convective mixing, eddies, and N2 fixation to explain the variation in observed elemental ratios. Convective mixing injected the cold water and enhanced the supply of N and P nutrients in the top layer (surface to 50–75 m) of the northern Arabian Sea. This led to a decrease in the N : P and C : P ratios in the particulate organic matter in the northern region, but C : N : P increased equatorward, averaging 245 : 32 : 1 in the top layer of the Arabian Sea. The variation in the elemental ratios in the top layer is best explained by the changes in water temperature. N2 fixation contributed negligibly to the N : P ratio of the export flux. The substantial decrease in N : P ratios of nutrients in the subsurface waters is most likely caused by the denitrifying conditions in the Arabian Sea. As the processes of convective mixing and eddies are are prevalent oceanic processes, our observations underpin that the interplay of these processes leads to changes in the elemental ratios globally. [ABSTRACT FROM AUTHOR]

The interactions between the internal tide and the mesoscale circulation are studied from the internal tide energy budget perspective. To that end, the modal energy budget of the internal tide is diagnosed using a high resolution numerical simulation covering the North Atlantic. Compared to the topographic contribution, the advection of the internal tide by the low‐frequency flow component and the horizontal and vertical shear are found to be significant at global scale, while the buoyancy contribution is important locally. The advection of the internal tide by the low‐frequency currents is responsible for a net energy transfer from the large scale to smaller scale internal tide, without exchanges with the low‐frequency flow. On the opposite, the shear of the mesoscale circulation and the buoyancy field are responsible for exchanges between the internal tide and the low‐frequency flow. The importance of the shear increases in the northernmost part of the domain, and a partial compensation between the buoyancy and the shear contributions is found in some areas of the North Atlantic, such as in the Gulf Stream region. In addition, the temporal variability of these energy transfers is investigated. In contrast to topographic scattering, for which the spring‐neap cycle is the dominant frequency, the energy transfer terms driven by low‐frequency motions in areas of strong mesoscale activity are also modulated by variations of the low‐frequency current spatial distribution. Plain Language Summary: Internal tides are waves generated when the tidal motions, induced by the Moon and Sun gravitational attraction on the ocean, interact with topographical features. These waves, which are ubiquitous in the ocean, then propagate inside of the ocean interior and interact with other types of currents. This study investigates the wave‐currents interactions and their impact on the energy exchanges in the ocean using a high‐resolution realistic numerical simulation. We find that, although the scattering of internal tides by the topography is dominant in their energy budget, the interaction with the currents are significant at the basin scale. They lead to transfers of energy from large scale internal waves to smaller ones, and between currents and internal tides in the Labrador sea and Faroe islands vicinity. This transfer of energy is likely to impact the ocean mixing driven by the small scale internal tide, and thereby to affect the global ocean circulation. The temporal variability of these current‐internal tides interactions shows the same period as the variability of the semi‐diurnal astronomical tide, with a modulation by the temporal variability of the currents. Key Points: Although topographic scattering is dominant, low‐frequency flow has a strong impact on the internal tide energy cycleThe dominant effect of low‐frequency flow on internal tide is a transfer of energy toward smaller scaleLow‐frequency flow induced interactions show an imprint of the spring neap cycle modulated by variations in the currents [ABSTRACT FROM AUTHOR]

As one kind of submesoscale instabilities, symmetric instability (SI) with spatiotemporal scales of O (100) m–O (1) km and O (1) hour exerts significant effects on vertical material transports and forward energy cascade in the ocean. The potential vorticity (PV) is an important conservative parameter controlling quasi‐geostrophic flows, whose budget can be modulated by SI. However, due to the small spatial scale of SI which is hardly resolved by most current observations and regional models, how SI affects the PV budget and how big the effect is remain unclear. In this work, the effect of SI on the PV budget in the surface mixed layer (SML) of the Kuroshio Extension region is quantitatively analyzed based on high‐resolution simulations by applying an existing SI parameterization scheme. Compared with the case without SI effects, negative PV is found to be eliminated in the SML in the SI‐parameterized case. The negative‐PV likelihood in the SI‐parameterized case is decreased by up to 12% due to SI. Analysis of the PV budget indicates that SI contributes to the PV budget mainly by modulating the friction term. The friction term tends to generate negative PV but its magnitude is decreased by 35% due to SI. Apart from the frictional term, both advection and non‐adiabatic terms are also found to be modulated by SI. This work sheds light on the contribution of SI in the PV budget in the ocean mixed layer and suggests a significant role of SI in quasi‐geostrophic PV dynamics. Plain Language Summary: Potential vorticity (PV) is a conservative but dynamically active tracer and an important indicator to characterize instability in the ocean. Traditionally, PV is applied to quasi‐geophysical dynamics, such as large‐scale circulation (O (1,000) km) and mesoscale currents (O (100) km). As the link between meso‐ and micro‐scales, submesoscale symmetric instability with spatial scales of O (100) m–O (1) km contribute to the PV budget, potentially modulating the quasi‐geostrophic dynamics. This work quantitatively evaluates the effect of SI on the PV budget by applying a parameterization scheme. Key Points: The effect of symmetric instability (SI) on potential vorticity is evaluated by applying a parameterization schemeElimination of negative potential vorticities is mainly caused by the friction effect with magnitude decreased by 35%The advection and non‐adiabatic effect are also found to be modulated by SI [ABSTRACT FROM AUTHOR]

Tropical Instability Waves (TIWs) have been shown to modulate upper ocean mixing. However, previous studies on the modulation of TIW related mixing are based on small numbers of TIWs and have not considered temporal variability, which can lead to discrepancies in the findings. In this study, using a 12‐year simulation carried out with a comprehensive, global, high‐resolution ocean model, we present evidence of seasonally modulated shear instabilities at TIW fronts in the Atlantic Ocean that reach down to the thermocline depth, potentially inducing mixing below the mixed layer depth. We find that, regardless of whether TIWs are present earlier in the year, frontal instabilities and potential mixing primarily occur in boreal summer, coinciding with a vertical shear maximum between the mean zonal currents. We argue that in the Atlantic Ocean, vertical shear at TIW fronts does usually not suffice to cause frontal instabilities below the mixed layer depth. Instead, the background shear needs to be sufficiently large in addition to TIW shear, to overcome the stability. The background shear in turn varies seasonally and is strongly driven by the variability of the northern branch of the South Equatorial Current (nSEC). As such, the variability of the nSEC strongly contributes to the generation and modulation of instabilities at TIW fronts that reach below the mixed layer depth and have the potential to induce mixing. Our results highlight the importance of seasonal variability when studying TIW impacts and their effect on mixing. Plain Language Summary: Tropical Instability Waves (TIWs) create alternating warm and cold patterns at the ocean's surface, leading to sharp temperature fronts. Past studies have shown that mixing can occur along these fronts, but how this mixing changes over time hasn't been explored. Our research focuses on how instabilities below the surface at TIW fronts vary over time. Studying such instabilities helps to better understand the processes that can impact the regional climate. We analyzed over 10 years of data from a detailed ocean model and found that instabilities at TIW fronts in the Atlantic Ocean can extend down to 100 m and follow a clear seasonal cycle. These instabilities are most common in summer when the difference in zonal current speeds is greatest. We discovered that shear at TIW fronts alone isn't enough to cause these deep instabilities; the overall background shear also needs to be strong. This background shear changes with the seasons, mainly influenced by the northern branch of the South Equatorial Current (nSEC). Thus, the nSEC plays a crucial role in modulating instabilities at TIW fronts and potential ocean mixing. Key Points: Tropical instability waves exhibit shear instabilities at the front that reach the thermocline depth and follow a distinct seasonal cycleThe seasonality of frontal shear instabilities results from a superposition of shear from the instability wave and background current shearVariability in the northern South Equatorial Current dominates the modulation of shear instabilities at Tropical Instability Wave fronts [ABSTRACT FROM AUTHOR]

Introduction: The impact of upper-ocean temperature on tropical cyclone (TC) activity is an open issue. Compared to the attention devoted to the effect of seasurface temperature (SST) on TC activities, much less is known about the effect of ocean mixed layer depth (OMLD) on TC activities, which is determined by the ocean temperature below the surface. Methods: In this study, a series of idealized numerical experiments were conducted to investigate the possible responses of TC activities to OMLD. Results: It was found that while OMLD exerts a minor influence on TC track, it evidently affects TC intensity, size, and destructiveness before reaching a certain OMLD threshold (approximately 15 m). Once the OMLD exceeds the threshold, changes in TC intensity, size, and destructiveness become marginal with further increase in OMLD. The threshold of OMLD is largely determined by TC intensity, which in turn is dictated by surface wind speed. Discussion: Specifically, before reaching the threshold of OMLD, the surface wind, namely TC-related surface wind, may bring the cold water from below the OMLD, and effectively decreases the upper ocean temperature (including the SST). As OMLD increases, the effect of surface wind on SST cooling gradually decreases, leading to an increase of SST below the TC. Subsequently, the SST increase leads to more surface enthalpy flux (SEF) input into the TC by increasing air-sea temperature and moisture differences. By altering TC's thermodynamic and dynamic structures, the increase of SEF eventually results in the increase of TC intensity and size, and thus its destructiveness. [ABSTRACT FROM AUTHOR]

Turbulent mixing in the ocean exerts an important control on the rate and structure of the overturning circulation. However, the balance of processes underpinning this mixing is subject to significant uncertainties, limiting our understanding of the overturning's deep upwelling limb. Here, we investigate the hitherto primarily neglected role of tens of thousands of seamounts in sustaining deep-ocean upwelling. Dynamical theory indicates that seamounts may stir and mix deep waters by generating lee waves and topographic wake vortices. At low latitudes, stirring and mixing are predicted to be enhanced by a layered vortex regime in the wakes. Using three realistic regional simulations spanning equatorial to middle latitudes, we show that layered wake vortices and elevated mixing are widespread around seamounts. We identify scalings that relate mixing rate within seamount wakes to topographic and hydrographic parameters. We then apply such scalings to a global seamount dataset and an ocean climatology to show that seamount-generated mixing makes an important contribution to the upwelling of deep waters. Our work thus brings seamounts to the fore of the deep-ocean mixing problem and urges observational, theoretical, and modeling efforts toward incorporating the seamounts' mixing effects in conceptual and numerical ocean circulation models. [ABSTRACT FROM AUTHOR]

Changes in sea ice conditions directly impact the way the wind transfers energy to the Arctic Ocean. The thinning and increasing mobility of sea ice is expected to change the size and speed of ridges on the underside of ice floes, called ice keels, which cause turbulence and impact upper-ocean stratification. However, the effects of changing ice keel characteristics on below-ice mixing are difficult to determine from sparse observations and have not been directly investigated in numerical or laboratory experiments. Here, for the first time, we examine how the size and speed of an ice keel affect the mixing of various upper-ocean stratifications using 16 two-dimensional numerical simulations of a keel moving through a two-layer flow. We find that the irreversible ocean mixing and the characteristic depth over which mixing occurs each vary significantly across a realistic parameter space of keel sizes, keel speeds, and ocean stratifications. Furthermore, we find that mixing does not increase monotonically with ice keel depth and speed but instead depends on the emergence and propagation of vortices and turbulence. These results suggest that changes to ice keel speed and depth may have a significant impact on below-ice mixing across the Arctic Ocean and highlight the need for more realistic numerical simulations and observational estimates of ice keel characteristics. [ABSTRACT FROM AUTHOR]

Based on Argo observations and a coupled atmosphere–ocean–wave model, the upper ocean responses to the tropical cyclone (TC) Mekunu (2018) were investigated, and the role of a pre-existing cold eddy in modulating the temperature response to TC Mekunu was quantified by employing numerical experiments. With TC Mekunu's passage, the mixed layer depth (MLD) on both sides of its track significantly deepened. Moreover, two cold patches (<26 °C) occurred, where the maximum cooling of the mixed layer temperature (MLT) reached 6.62 °C and 6.44 °C. Both the MLD and MLT changes exhibited a notable rightward bias. However, the changes in the mixed layer salinity (MLS) were more complex. At the early stage, the MLS on both sides of the track increased by approximately 0.5 psu. When TC Mekunu made landfall, the MLS change around the track was asymmetric. Significantly, a cold eddy pre-existed where the second cold patch emerged, and this eddy was intensified after TC Mekunu's passage, with an average sea surface height reduction of approximately 2.7 cm. By employing the stand-alone ocean model, the numerical experimental results demonstrated that the pre-existing cold eddy enhanced TC-induced MLT cooling by an average of approximately 0.41 °C due to steeper temperature stratification at the base of mixed layer. Moreover, heat budget analysis indicated that the pre-existing cold eddy also enhanced subsurface temperature cooling mainly through zonal advection. [ABSTRACT FROM AUTHOR]

Diapycnal mixing in the ocean interior is largely fueled by internal tides. Mixing schemes that represent the breaking of internal tides are now routinely included in ocean and earth system models applied to the modern and future. However, this is more rarely the case in climate simulations of deep‐time intervals of the Earth, for which estimates of the energy dissipated by the tides are not always available. Here, we present and analyze two IPSL‐CM5A2 earth system model simulations of the Early Eocene made under the framework of DeepMIP. One simulation includes mixing by locally dissipating internal tides, while the other does not. We show how the inclusion of tidal mixing alters the shape of the deep ocean circulation, and thereby of large‐scale biogeochemical patterns, in particular oxygen distributions. In our simulations, the absence of tidal mixing leads to a relatively stagnant and poorly ventilated deep ocean in the North Atlantic, which promotes the development of a basin‐scale pool of oxygen‐deficient waters, at the limit of complete anoxia. The absence of large‐scale anoxic records in the deep ocean after the Cretaceous anoxic events suggests that such an ocean state most likely did not occur at any time across the Paleogene. This highlights how crucial it is for climate models applied to the deep‐time to integrate the spatial variability of tidally driven mixing as well as the potential of using biogeochemical models to exclude aberrant dynamical model states. Key Points: Inclusion of realistic near‐field tidal mixing substantially modifies global deep ocean circulation in the Early EoceneThese tidally driven changes yield significantly different biogeochemical properties of water masses, in particular in the AtlanticThe simulation that includes tidal mixing compares more favorably to inferences from the O2 proxy record [ABSTRACT FROM AUTHOR]

This study represents the first use of the artificial radionuclides 129I and 236U, released into the ocean mainly from Nuclear Reprocessing Plants, as a dual tracer in the vicinity of Iceland with novel estimation of ocean circulatory pathways and mixing in the region. Iceland lies at the gateway to the Arctic where warm, saline Atlantic waters interact with waters of Arctic origin in ways that have critical consequences for the strength and stability of the Atlantic Meridional Overturning Circulation. Many of these interactions are not yet fully understood, such as how Atlantic water circulates around the Arctic Ocean and Nordic Seas and the composition and fate of the major overflows of the Greenland‐Scotland Ridge. Using new and previous measurements of 129I and 236U in seawater, we present a new method of appraising water mass provenance and mixing in the form of the 129I–236U dual mixing plot. With this method, we estimate that at least half the Atlantic‐origin water entering the Arctic Ocean circulates around the Canada Basin before exiting at Fram Strait and that this outflow is increased by about 40% by mixing with Return Atlantic Water "short‐circuiting" the Arctic Ocean at Fram Strait. We present tracer‐based evidence that water carried by the East Greenland Current has an unbroken pathway to the Faroe‐Shetland Channel and that Iceland‐Scotland Overflow Water (ISOW) entrains 60% Labrador Sea Water during transit past southeast Iceland. We present an unambiguous way to differentiate ISOW from DSOW after they partially merge in the Irminger Sea. Plain Language Summary: This study is the first to use a pair of man‐made nuclear products to study ocean currents near Iceland. Understanding these currents is important because Iceland sits where warm Atlantic waters meet colder Arctic waters, affecting a key ocean circulation system and the global climate. However, many details about how these waters move and mix are still unclear. Using new and previous measurements of these tracers in seawater, we demonstrate a new method to estimate the origins of currents and how they mix. We estimate that the outflow of Atlantic‐origin water from the Arctic Ocean at Fram Strait is increased by about 40% by Atlantic water bypassing the Arctic Ocean altogether at this location. Some of this Atlantic‐origin water then flows all the way from Iceland to Shetland in an unbroken pathway and the water that spills over the Iceland‐Scotland ridge increases by 60% by mixing with water from the Labrador Sea southeast of Iceland. Finally, we present a new way to distinguish this Iceland‐Scotland Overflow Water from similar water that overflows the ridge at Denmark Strait. Key Points: Recirculating Atlantic Water contributes about 40% to Atlantic layer of the East Greenland CurrentAtlantic‐origin waters can be traced from the East Greenland Current to the Faroe‐Shetland ChannelIceland‐Scotland Overflow Water and Denmark Strait Overflow Water are traceable separately after joining Deep Western Boundary Current [ABSTRACT FROM AUTHOR]

North Atlantic subtropical mode water (NASTMW) serves as a major conduit for dissolved carbon to penetrate into the ocean interior by its wintertime outcropping events. Prior research on NASTMW has concentrated on its physical formation and destruction, as well as Lagrangian pathways and timescales of water into and out of NASTMW. In this study, we examine how dissolved inorganic carbon (DIC) concentrations are modified along Lagrangian pathways of NASTMW on subannual timescales. We introduce Lagrangian parcels into a physical‐biogeochemical model and release these parcels annually over two decades. For different pathways into, out of, and within NASTMW, we calculate changes in DIC concentrations along the path (ΔDIC), distinguishing contributions from vertical mixing and biogeochemical processes. The strongest ΔDIC is during subduction of water parcels (+101 μmol L−1 in 1 year), followed by transport out of NASTMW due to increases in density in water parcels (+10 μmol L−1). While the mean ΔDIC for parcels that persist within NASTMW in 1 year is relatively small at +6 μmol L−1, this masks underlying dynamics: individual parcels undergo interspersed DIC depletion and enrichment, spanning several timescales and magnitudes. Most DIC enrichment and depletion regimes span timescales of weeks, related to phytoplankton blooms. However, mixing and biogeochemical processes often oppose one another at short timescales, so the largest net DIC changes occur at timescales of more than 30 days. Our new Lagrangian approach complements bulk Eulerian approaches, which average out this underlying complexity, and is relevant to other biogeochemical studies, for example, on marine carbon dioxide removal. Plain Language Summary: Mode waters are relatively thick water masses with homogeneous properties, such as temperature and salinity. The North Atlantic subtropical mode water (NASTMW), found in the Sargasso Sea, is one such water mass. Lying underneath the ocean surface, it comes into contact with the atmosphere during winter, when the surface layer is vigorously mixed due to strong winds, causing the mixed layer to connect with NASTMW. This way, NASTMW can buffer atmospheric temperature and carbon anomalies during the summer, when there is no surface connection. It is also a conduit for carbon to penetrate beneath the ocean's upper mixed layer, with the potential to sequester it. We study NASTMW from the viewpoint of a water parcel that moves with the currents and see how carbon concentrations in the water parcels change along different NASTMW pathways. For each pathway, the carbon concentration changes due to an interplay of vertical mixing and biogeochemical processes, for example, related to plankton growth and decay. These processes can unfold over different timescales and may counteract or enhance themselves or one another. The largest change in carbon concentration is found when a parcel moves from the upper ocean mixed layer into NASTMW, mostly due to vertical mixing. Key Points: Carbon transformations along pathways of North Atlantic Subtropical Mode Water are split into mixing and biogeochemical contributionsAlong paths into, within, and out of this mode water, mixing and biogeochemistry alter carbon in water parcels over a range of timescalesEnrichment is highest during mixed layer subduction, which few parcels undergo annually; persistence in mode water is the dominant pathway [ABSTRACT FROM AUTHOR]

Small‐scale turbulent mixing supplies potential energy for the upwelling of deep waters in the abyssal ocean, a key component of the global overturning circulation. This process is particularly significant in critical regions such as the Northwestern Pacific where the upwelling structure of deep waters remains poorly understood due to limited knowledge of deep ocean mixing. Here, we investigate the full‐depth spatiotemporal variability of turbulent mixing in the deep Northwestern Pacific based on hydrographic data collected over repeated surveys. Nineteen‐year‐average diapycnal diffusivity of 1.42 × 10−4 m2 s−1 is reveled in the deep Philippine Sea, indicating significantly stronger mixing compared to the stratified ocean interior. Spatially, turbulent mixing strengthens toward the bottom and intensifies westward from the open Pacific to the Philippine Sea due to rough topography. At certain mixing hotspots, enhanced mixing can penetrate up to 2,500 m above the bottom, suggesting a substantial potential for upwelling. Below 2,000 m, turbulent mixing exhibits pronounced seasonal variation that deep mixing is more intense in summer (winter) than in winter (summer) in the West Caroline Basin (the Parece Vela Basin). This spatially varying seasonality may be attributed to the inhomogeneous internal tidal energy dissipation in the Northwestern Pacific. Our study will serve to clarify the modulation of turbulent mixing to deep‐water mass transformation and circulation in the Northwestern Pacific. Plain Language Summary: The global overturning circulation plays an important role in transporting and redistributing climate‐sensitive matters such as oxygen and carbon. This circulation is enclosed by the upwelling of deep waters facilitated by turbulence in the deep ocean, which mixes dense waters with the lighter waters above. In the Northwestern Pacific, a significant portion of deep waters from the Southern Ocean continuously intrudes into the Philippine Sea, rendering it a critical region for global circulation. However, the upwelling process of deep water remains largely unknown due to the limited understanding of turbulent mixing in the deep Northwestern Pacific. To fill this gap, we characterize the spatiotemporal distribution of turbulent mixing in the Northwestern Pacific based on long‐term observations during 2004–2022. We reveal that turbulent mixing here is highly enhanced by rough topography, 10 to 100 times stronger than that in the open ocean, particularly near some key terrains such like Kyushu‐Palau Ridge. Besides, a region‐dependent seasonal variation of deep mixing is also observed. Our study contributes to a further understanding of the abyssal mixing and circulation in the Northwestern Pacific. Key Points: A full‐depth distribution of turbulent mixing in the Northwestern Pacific is obtained based on 19‐year repeated hydrographic observationsRough topography in the Philippine Sea enhances the deep mixing, with diffusivities even up to O (10−2) m2s−1 near the bottomLatitude dependent seasonality of deep ocean turbulent mixing is revealed, which could be related to spatial variations of internal tides [ABSTRACT FROM AUTHOR]

Analysis of dissolved oxygen (O2) in the Arctic's surface ocean provides insights into gas transfer between the atmosphere‐ice‐ocean system, water mass dynamics, and biogeochemical processes. In the Arctic Ocean's Canada Basin mixed layer, higher O2 concentrations are generally observed under sea ice compared to open water regions. Annual cycles of O2 and O2 saturation, increasing from summer through spring and then sharply declining to late summer, are tightly linked to sea ice cover. The primary fluxes that influence seasonal variability of O2 are modeled and compared to Ice‐Tethered Profiler O2 observations to understand the relative role of each flux in the annual cycle. Findings suggest that sea ice melt/growth dominates seasonal variations in mixed layer O2, with minor contributions from vertical entrainment and atmospheric exchange. While the influence of biological activity on O2 variability cannot be directly assessed, indirect evidence suggests relatively minor contributions, although with significant uncertainty. Past studies show that O2 molecules are expelled from sea ice during brine rejection; sea ice cover can then inhibit air‐sea gas exchange resulting in winter mixed layers that are super‐saturated. Decreasing mixed layer O2 concentrations and saturation levels are observed during winter months between 2007 and 2019 in the Canada Basin. Only a minor portion of the decreasing trend in wintertime O2 can be attributed to decreased solubility. This suggests the O2 decline may be linked to more efficient air‐sea exchange associated with increased open water areas in the winter sea ice pack that are not necessarily detectable via satellite observations. Plain Language Summary: Dissolved oxygen (O2) is a valuable ocean property that allows us to better understand the exchange of gases between the different ocean layers, sea ice, and atmosphere, and the physical and biological processes that control its variability. Understanding how and why O2 concentrations in the Arctic Ocean mixed layer vary spatially and seasonally is crucial for interpreting its evolution over timescales of years to decades that are influenced by global warming. We use physical and thermodynamical relationships to model the main factors that influence O2 concentrations in the mixed layer of the Arctic Ocean's Canada Basin, which we compare to observations made by Ice‐Tethered Profilers. Model results indicate that seasonal variations in O2 concentrations are dominated by the effects of sea ice growth and melt. Other processes that modulate mixed layer O2, including air‐sea exchange and ocean mixing, have a lesser influence. Between 2007 and 2019, mixed‐layer O2 has decreased in winter months, which we attribute to more openings in the sea‐ice pack during wintertime in the Canada Basin. Key Points: Spatial and seasonal distributions of O2 concentrations in the Canada Basin mixed layer are linked to the seasonal evolution of sea iceModeled fluxes suggest brine rejection and meltwater dilution during sea ice melt/formation dominate seasonal variability of mixed layer O2Decreases in mixed‐layer O2 during winter (over 2007–2019) suggest outgassing, likely driven by changes in the wintertime sea ice pack [ABSTRACT FROM AUTHOR]

The Arctic Ocean and adjacent seas are undergoing increased freshwater influx due to enhanced glacial and sea ice melt, precipitation, and runoff. Accurate delineation of these freshwater sources is vital as they critically modulate ocean composition and circulation with widespread and varied impacts. Despite this, the delineation of freshwater sources using physical oceanographic measurements (e.g., temperature, salinity) alone is challenging and there is a requirement to improve the partitioning of ocean water masses and their mixing relationships. Here, we complement traditional oceanographic measurements with continuous surface seawater isotopic analysis (δ18O and deuterium excess) across a transect extending from coastal Alaska to Baffin Bay and the Labrador Sea conducted from the US Coast Guard Cutter Healy in Autumn 2021. We find that the diverse isotopic signatures of Arctic freshwater sources, coupled with the high freshwater proportion in these marine systems, facilitates detailed fingerprinting and partitioning. We observe the highest freshwater composition in the Beaufort Sea and Amundsen Gulf regions, with heightened freshwater content in eastern Baffin Bay adjacent to West Greenland. We apply isotopic analysis to delineate freshwater sources, revealing that in the Western Arctic freshwater inputs are dominated by meteoric water inputs—specifically the Mackenzie River—with a smaller sea ice meltwater component and in Baffin Bay the primary sources are local precipitation and glacial meltwater discharge. We demonstrate that such freshwater partitioning cannot be achieved using temperature‐salinity relationships alone, and highlight the potential of seawater isotopic tracers to assess the roles and importance of these evolving freshwater sources. Plain Language Summary: Freshwater inputs to the Arctic seas, including glacial and sea ice meltwater, precipitation, and river runoff, are increasing as the Arctic warms. The impacts of these changing freshwater influxes are varied depending on the type of freshwater source, and thus it is important to delineate and trace these different freshwater sources, which represents a significant challenge using only traditional physical oceanographic measurements (e.g., temperature, salinity). In this study, we utilize a new approach to identify and trace freshwater sources using continuous seawater isotopic measurements during a cruise extending from coastal Alaska, through the Canadian Archipelago, and across Baffin Bay and the Labrador Sea. We show that these isotopic measurements, which have been commonly used in other media (e.g., precipitation, water vapor, ice cores), hold important and distinct information about the source and mixing of different freshwater sources. We use these measurements to identify the freshwater sources (e.g., Mackenzie vs. Yukon River) contributing to ocean surface waters across the Arctic region. Key Points: Seawater isotopic measurements (δ18O, δ2H, deuterium excess) show heightened freshwater content in the Beaufort Sea and Baffin BayIsotopic observations enable freshwater source delineation not feasible from traditional physical oceanographic methodsFreshwater source delineation includes the Mackenzie and Yukon Rivers around coastal Alaska and glacial meltwater in Baffin Bay [ABSTRACT FROM AUTHOR]

Anthropogenic perturbations from fossil fuel burning, nuclear bomb testing, and chlorofluorocarbon (CFC) use have created useful transient tracers of ocean circulation. The atmospheric 14C/C ratio (∆14C) peaked in the early 1960s and has decreased now to pre‐industrial levels, while atmospheric CFC‐11 and CFC‐12 concentrations peaked in the early 1990s and early 2000s, respectively, and have now decreased by 10%–20%. We present the first analysis of a decade of new observations (2007 to 2018–2019) and give a comprehensive overview of the changes in ocean ∆14C and CFC concentration since the WOCE surveys in the 1990s. Surface ocean ∆14C decreased at a nearly constant rate from the 1990–2010s (20‰/decade). In most of the surface ocean ∆14C is higher than in atmospheric CO2 while in the interior ocean, only a few places are found to have increases in ∆14C, indicating that globally, oceanic bomb 14C uptake has stopped and reversed. Decreases in surface ocean CFC‐11 started between the 1990 and 2000s, and CFC‐12 between the 2000–2010s. Strong coherence in model biases of decadal changes in all tracers in the Southern Ocean suggest ventilation of Antarctic Intermediate Water was enhanced from the 1990 to the 2000s, whereas ventilation of Subantarctic Mode Water was enhanced from the 2000 to the 2010s. The decrease in surface tracers globally between the 2000 and 2010s is consistently stronger in observations than in models, indicating a reduction in vertical transport and mixing due to stratification. Plain Language Summary: The ocean contains many dissolved gases that can be measured by sampling ocean water from ships. Some of these dissolved gases are changing over time because their atmospheric concentrations are changing. By measuring these changes in the ocean, we can learn about ocean transport and mixing, which is an important factor regulating climate change. We describe ocean measurements from the past three decades of radiocarbon in dissolved inorganic carbon, that was affected by nuclear bomb testing, and of chlorofluorocarbon gases, CFCs, that were banned in the 1990s due to their destructive effects on the ozone layer. The measurements show how their oceanic distributions have changed, as a result of their atmospheric histories and the patterns of ocean uptake and transport. By comparing the measurements with ocean models, we can identify biases in the models that affect how well future climate change can be predicted. In the most recent decade, the presence of all three tracers is shown to be decreasing in the upper ocean, a reversal of the preceding decades. Also in the most recent decade, the decrease in surface tracers is stronger in the observations than the models, suggesting a change to upper ocean mixing and stratification. Key Points: Recent trends in upper ocean Δ14C, pCFC‐11, and pCFC‐12 are negative, reflecting their decreasing atmospheric trendsIncreases in ∆14C are only observed in a few places over 2000–2010s, showing oceanic bomb 14C uptake has stopped and reversedModel‐data Δ14C and chlorofluorocarbon differences could be consistent with decadal ventilation changes and large‐scale stratification or model errors [ABSTRACT FROM AUTHOR]

Using 9 years of full‐depth profiles from 55 Deep Argo floats in the Southwest Pacific Basin collected between 2014 and 2023, we find consistent warm anomalies compared to a long‐term climatology below 2,000 m ranging between 11 ± 2 to 34 ± 2 m°C, most pronounced between 3,500 and 5,000 m. Over this period, a cooling trend is found between 2,000 and 4,000 m and a significant warming trend below 4,000 m with a maximum rate of 4.1 ± 0.31 m°C yr−1 near 5,000 m, with a possible acceleration over the second half of the period. The integrated Steric Sea Level expansion below 2,000 m was 7.9 ± 1 mm compared to the climatology with a trend of 1.3 ± 1.6 mm dec−1 over the Deep Argo era, contributing significantly to the local sea level budget. We assess the ability to close a full Sea Level Budget, further demonstrating the value of a full‐depth Argo array. Plain Language Summary: Cold, dense waters formed near polar regions in both hemispheres, sink to great depths and fill‐up the majority of the world's deep ocean. Compilation of sparse observations of temperature from global ship‐based surveys at roughly 10‐year intervals worldwide have shown that sequestration of excess atmospheric heat into the deep ocean has caused these waters to warm steadily since the 1990's into the Present. Not only does this warming have implications for changes in large scale ocean circulation, but is also associated with warming‐induced sea level rise. Using a new data set collected between 2014 and 2023 from 55 freely drifting robotic floats (Deep Argo) which gather crucial bimonthly temperature and salinity data between the surface ocean and the ocean floor, we find the greatest warming trend at a depth of 5,000 m of 4 ± 0.3 m°C yr−1 and an associated sea level rise rate below 2,000 m of 1.3 ± 1.6 mm dec−1. Deep Argo data being collected in ocean basins worldwide are crucial in providing high resolution data of the warming deep ocean and its implications on global sea level, ocean mixing and large‐scale ocean circulation. Key Points: Nine years of Deep Argo data in the S.W. Pacific reveals continued warming in the abyss while the mid‐depths cooledWaters below 4,000 m show an accelerated warming trend with a maximum overall warming rate of 4.1 ± 0.31 m°C yr−1 at 5,000 mDeep ocean steric expansion contributed 1.3 ± 1.6 mm dec−1 to total the local sea level [ABSTRACT FROM AUTHOR]

Amplification of sea surface temperature (SST) seasonality in response to global warming is a robust feature in climate model projections but season‐dependent regional disparities in this amplification and the associated mechanisms are not well addressed. Here, by analyzing large ensemble simulations using Community Earth System Model version 2, we investigate detailed spatiotemporal characteristics of the amplification of SST seasonality focusing on the North Pacific and North Atlantic, where robust changes are projected to emerge around 2050 under SSP3‐7.0 scenario. Our results indicate that atmosphere‐ocean coupled processes shape regional changes in SST seasonality differently between warm (MAMJJAS) and cold seasons (ONDJF). During the warm season, the projected warming tendency is mainly due to increased net surface heat flux and weakening of vertical mixing. On the other hand, in the cold season, the projected cooling tendency is driven by strengthened vertical mixing over the North Pacific associated with the northward shift of storm tracks but weakened horizontal advection and mixing due to changes in ocean currents over the North Atlantic. Plain Language Summary: In addition to the increase in annual mean sea surface temperature (SST), the amplitude of SST seasonal cycle is projected to intensify in response to greenhouse warming over many parts of the global ocean due to more warming during local summer than winter. By analyzing large ensemble simulations, we show how different atmospheric and ocean processes influence spatiotemporal differences in SST seasonality changes in the North Pacific and North Atlantic. Our findings indicate that during local summer, the projected warming tendency intensifies with the increased heat input into the ocean and weakening of ocean mixing. Meanwhile, during local winter, the projected cooling tendency intensifies with the strengthening of vertical mixing associated with the enhanced storm track activities in the North Pacific, but the weakening of ocean horizontal advection and mixing attributable to the weakened ocean circulation over the North Atlantic. They together contribute to the amplification of SST seasonality, which could have implications for marine ecosystems, including the plankton phenology. Key Points: Using CESM2 large ensemble simulations, this study explores the robustness of future sea surface temperature (SST) seasonality changes and the associated atmospheric and surface ocean processesIn the warm season, the future SST warming is mainly driven by the reduced upper ocean mixing and the increased net surface heat influx into the oceanIn the cold season, the future SST cooling is attributable to strengthened vertical mixing over the North Pacific but weakened horizontal advection and mixing in the North Atlantic [ABSTRACT FROM AUTHOR]

Mesoscale eddies are known to enhance diapycnal mixing in the ocean, yet direct observation of this effect remains a significant challenge, especially in the robust Antarctic Circumpolar Current (ACC). To quantify the diapycnal mixing induced by mesoscale eddies in the standing meander of the ACC, satellite altimeter and Argo profile data were combined to composite eddies, where the 1.6 m dynamic height contour was used for the first time instead of the climatological Northern Sub-Antarctic Front (SAFN) to define the northern boundary of the ACC to eliminate the influence of frontal shift. The 3D structures of the composite anticyclonic/cyclonic eddy (CAE/CCE) were obtained. Both the CAE and CCE were similar in shape to Taylor columns, from sea surface to the neutral surface of 28.085 kgm − 3 (1689 ± 66 dbar) for the CAE, and from sea surface to 28.01 kgm − 3 (1491 ± 202 dbar) for the CCE. On the same neutral surface, the diffusivity (κ) inside the CCE was one to two orders of magnitude higher than that inside the CAE. Vertically, the maximum influence depth of the CCE on κ reached 1200 dbar, while for the CAE, it reached 800 dbar, where κ exceeded O (10 − 4)   m 2 s − 1 , and κ gradually decreased from these depths downwards. [ABSTRACT FROM AUTHOR]

The non‐local model of mixing based on internal wave breaking, IDEMIX, is implemented as an enhancement of a turbulent kinetic energy closure model in three non‐eddy resolving general circulation ocean models that differ in the discretization and choice of computational grids. In IDEMIX internal wave energy is generated by an energy flux resulting from near‐inertial waves induced by wind forcing at the surface, and at the bottom, by an energy flux that parameterizes the transfer of energy between baroclinic and barotropic tides. In all model simulations with IDEMIX, the mixing work is increased compared to the reference solutions without IDEMIX, reaching values in better agreement with finestructure observations. Furthermore, the horizontal structure of the mixing work is more realistic as a consequence of the heterogeneous forcing functions. All models with IDEMIX simulate deeper thermocline depths related to stronger shallow overturning cells in the Indo‐Pacific. In the North Atlantic, deeper mixed layers in simulations with IDEMIX are associated with an increased Atlantic overturning circulation and an increase of northward heat transports toward more realistic values. The response of the deep Indo‐Pacific overturning circulation and the weak bottom cell of the Atlantic to the inclusion of IDEMIX is incoherent between the models, suggesting that additional unidentified processes and numerical mixing may confound the analysis. Applying different tidal forcing functions leads to simulation differences that are small compared to differences between the different models or between simulations with IDEMIX and without IDEMIX. Plain Language Summary: Waves in the ocean interior play a fundamental role for ocean dynamics since they can carry energy over long distances and, once they break, lead to turbulent mixing. This turbulent mixing can cause dense water masses to rise from the deep ocean with a direct impact on large‐scale currents. The wave dynamics occur on spatial scales that cannot be resolved in global ocean or climate models. To account for these processes, we apply the new parameterization IDEMIX that describes internal wave generation, propagation, and mixing. Using three different ocean models with and without IDEMIX ensures that we can identify model‐specific effects of the parameterization and discriminate them from those independent of the model. We find that the simulated mixing patterns agree better with observations once IDEMIX is applied. Large‐scale currents and the vertical temperature distribution are substantially affected by the internal wave parameterization. Whether this leads to an improved agreement with observed currents and water mass properties depends on the specific model and on numerical effects. In most cases, simulations with IDEMIX are not very sensitive to details of how the internal wave model is driven by tidal energy input. Key Points: the IDEMIX closure for a consistent representation of internal wave‐induced mixing is evaluated in three state‐of‐the‐art ocean modelsonly in simulations with IDEMIX can the models reproduce the magnitude and spatial variability of the observed mixing workmost changes with IDEMIX can be attributed to stronger mixing, but some effects are confounded by other processes and numerical mixing [ABSTRACT FROM AUTHOR]

Radium‐226(226Ra) and barium (Ba) exhibit similar chemical behaviors and distributions in the marine environment, serving as valuable tracers of water masses, ocean mixing, and productivity. Despite their similar distributions, these elements originate from distinct sources and undergo disparate biogeochemical cycles, which might complicate the use of these tracers. In this study, we investigate these processes by analyzing a full‐depth ocean section of 226Ra activities (T1/2 = 1,600 years) and barium concentrations obtained from samples collected along the US GEOTRACES GP15 Pacific Meridional Transect during September–November 2018, spanning from Alaska to Tahiti. We find that surface waters possess low levels of 226Ra and Ba due to export of sinking particulates, surpassing inputs from the continental margins. In contrast, deep waters have higher 226Ra activities and Ba concentrations due to inputs from particle regeneration and sedimentary sources, with 226Ra inputs primarily resulting from the decay of 230Th in sediments. Further, dissolved 226Ra and Ba exhibit a strong correlation along the GP15 section. To elucidate the drivers of the correlation, we used a water mass analysis, enabling us to quantify the influence of water mass mixing relative to non‐conservative processes. While a significant fraction of each element's distribution can be explained by conservative mixing, a considerable fraction cannot. The balance is driven using non‐conservative processes, such as sedimentary, rivers, or hydrothermal inputs, uptake and export by particles, and particle remineralization. Our study demonstrates the utility of 226Ra and Ba as valuable biogeochemical tracers for understanding ocean processes, while shedding light on conservative and myriad non‐conservative processes that shape their respective distributions. Key Points: Water mass analysis shows significant shallow deficits and deep excesses of dissolved 226Ra and Ba in the Pacific Ocean along the GP15Shallow deficits seem driven by BaSO4 precipitation, deep excesses reflect benthic inputs that are more significant for 226Ra than for BaBudgets suggest that regional barium fluxes nearly balance, and the Pacific is likely a significant source of radium‐226 to the global ocean [ABSTRACT FROM AUTHOR]

Turbulent mixing in the ocean, lakes and reservoirs facilitates the transport of momentum, heat, nutrients, and other passive tracers. Turbulent fluxes are proportional to the rate of turbulent kinetic energy dissipation per unit mass, ε. A common method for ε measurements is using microstructure profilers with shear probes. Such measurements are now widespread, and a non-expert practitioner will benefit from best practice guidelines and benchmark datasets. As a part of the Scientific Committee on Oceanographic Research (SCOR) working group on "Analysing ocean turbulence observations to quantify mixing" (ATOMIX), we compiled a collection of five benchmark data of ε from measurements of turbulence shear using shear probes. The datasets are processed using the ATOMIX recommendations for best practices documented separately. Here, we describe and validate the datasets. The benchmark collection is from different types of instruments and covers a wide range of environmental conditions. These datasets serve to guide the users to test their ε estimation methods and quality-assurance metrics, and to standardize their data for archiving. [ABSTRACT FROM AUTHOR]

Oceanic mesoscale eddy mixing plays a crucial role in Earth's climate system by redistributing heat, salt, and carbon. For many ocean and climate models, mesoscale eddies still need to be parameterized. This is often done via an eddy diffusivity K , which sets the strength of turbulent downgradient tracer fluxes. A well-known effect is the modulation of K in the presence of background potential vorticity (PV) gradients, which suppresses cross-PV gradient mixing. Topographic slopes can induce such suppression through topographic PV gradients. However, this effect has received little attention, and topographic effects are often not included in parameterizations for K. In this study, we show that it is possible to describe the effect of topography on K analytically in a barotropic framework, using a simple stochastic representation of eddy–eddy interactions. We obtain an analytical expression for the depth-averaged K as a function of the bottom slope, which we validate against diagnosed eddy diffusivities from a numerical model. The obtained analytical expression can be generalized to any constant barotropic PV gradient. Moreover, the expression is consistent with empirical parameterizations for eddy diffusivity over topography from previous studies and provides a physical rationalization for these parameterizations. The new expression helps to understand how eddy diffusivities vary across the ocean, and thus how mesoscale eddies impact ocean mixing processes. Significance Statement: Large oceanic "whirls," called eddies, can mix and transport ocean properties such as heat, salt, carbon, and nutrients. Mixing plays an important role for oceanic ecosystems and the climate system. In numerical simulations of Earth's climate, eddy mixing is typically represented using a simplified expression. However, an effect that is often not included is that eddy mixing is weaker over a sloping seafloor. In most areas of the ocean the bottom slope is steep enough for this effect to be significant. In this study we derive an expression for eddy mixing that accounts for oceanic bottom slopes. The present effort provides a physical basis for eddy mixing over oceanic bottom slopes, justifying their use in climate models. [ABSTRACT FROM AUTHOR]

Strongly coupled data assimilation allows observations of one Earth system component (e.g., the ocean) to directly update another component (e.g., the atmosphere). The majority of the information transfer in strongly coupled atmosphere‐ocean systems is passed through vertical correlations between atmospheric boundary layer and ocean mixed layer fields. In this work we use correlations from a global, coupled model to study vertical observation‐space localization techniques for strongly coupled data assimilation. We generate target correlations using a bootstrapping approach from a single 24 hr forecast from a realistic global, weakly coupled atmosphere‐ocean cycling system with an 80‐member ensemble, which is the ensemble size currently used by the NOAA operational global data assimilation system. We compare data assimilation methods with different localization schemes using single‐update, offline experiments. We develop a new strategy for optimal observation space localization, called Empirical Optimal R‐localization (EORL), to give an upper bound on the improvement we can expect with any localization scheme. We then evaluate Gaspari‐Cohn localization, which is a commonly used parametric localization function and review its performance with respect to the optimal localization scheme. We investigate how the performance of these localization strategies changes with increasing ensemble sizes. Our results show that strongly coupled data assimilation has the potential to be an improvement over weakly coupled data assimilation when large ensembles are used. We also show that the Gaspari‐Cohn localization function does not appear to be a particularly good choice for cross‐fluid vertical localization. Plain Language Summary: Accurate Earth system forecasts rely on accurate estimates of the current state of the system. This initial state is estimated through a process called data assimilation, which combines the previous forecast with current observations. This process relies on accurate estimates of uncertainty in both the model and the observations. Inaccurate estimates of model and observation uncertainty can result in a degradation of future forecasts. A technique called localization is widely used to minimize the impact of distant observations due to unreliable uncertainty estimates over long distances. In this work we study localization in the context of a global, coupled atmosphere‐ocean model. In particular, we look at strongly coupled data assimilation, where observations of the ocean are used to update the atmospheric model state and vice versa. This is in contrast to weakly coupled data assimilation, where observations of the ocean are only used to update the ocean model state. We investigate different types of localization and find that a commonly used localization function does not perform particularly well in strongly coupled model initialization. Key Points: We investigate vertical localization for strongly coupled atmosphere‐ocean data assimilation in a realistic global modelStrong coupling can improve data assimilation effectiveness over weak coupling when large ensembles are usedWe present a method for optimal observation space localization, called EORL, and demonstrate its performance in offline experiments [ABSTRACT FROM AUTHOR]

Tropical sea surface temperature (SST) biases can cause atmospheric biases on global scales, hence SST needs to be represented well in climate models. A major source of uncertainties is the representation of turbulent mixing in the oceanic boundary layer, or mixed layer (ML). In the present study we focus on near‐inertial wave (NIW) induced mixing. The performance of two mixing schemes, Turbulent Kinetic Energy and K‐profile parameterization (KPP), is assessed at two sites (11.5°N, 23°W and 15°N, 38°W) in the tropical Atlantic. At 11.5°N, turbulence observations (eddy diffusivities, shear and stratification) are available for comparison. We find that the schemes differ in their representation of NIWs, but both under‐represent the observed enhanced diffusivities below the observed ML. However, we find that the models do mix below the ML at 15°N when a storm passes nearby. The near‐inertial oscillations remain below the ML for the following 10 days. Near‐inertial kinetic energy (NIKE) biases in the models are not directly correlated with the wind speed, the MLD biases, or the stratification at the ML base. Instead, NIKE biases are sensitive to the vertical mixing scheme parameterization. NIKE biases are lowest when the KPP scheme is used. Plain Language Summary: The surface temperature of the ocean is highly dependent on the depth of the mixed layer (ML), the uppermost layer in the water column, where density, temperature and salinity are approximately constant. In climate models, the vertical mixing processes cannot be resolved, and instead they are computed with the use of vertical mixing schemes. We assess how well two of such schemes can represent the mixing induced by a specific type of ocean waves, near‐inertial waves (NIWs). We compare recent observations of turbulent mixing induced by NIWs in the tropical Atlantic with numerical simulations that resolve storms. Our results show that the models are able to reproduce the observed NIWs, but underestimate their mixing and amplitude. Our analysis also shows that NIWs are a driver of mixing below the uppermost ocean layer in the models. The strength of the near‐inertial currents is sensitive to the vertical mixing parameterization. Key Points: Observations of inertial oscillations are used to evaluate the performance of two vertical mixing schemes in two high‐resolution modelsBoth the K‐profile parameterization and the Turbulent Kinetic Energy closure underestimate the NIW‐induced mixingNear‐inertial kinetic energy biases are sensitive to the vertical mixing parameterization [ABSTRACT FROM AUTHOR]

The Gregg‐Henyey‐Polzin (GHP) finescale parameterization is widely employed to infer turbulent mixing when microstructure measurements are unavailable. However, the strain‐only GHP scaling is the only option when shear information is lacking, in which case the shear‐to‐strain ratio (Rω) is commonly treated as a constant. Since Rω has been reported to vary spatiotemporally, using a fixed value of Rω might result in a significant bias in inferring turbulent mixing. In this study, we present a global map of Rω, based on the microstructure database contributed by the Climate Process Team on internal wave‐driven ocean mixing, which covers the Indian, Pacific, and Atlantic oceans across both high and low latitudes. Then, we propose a parameterization of Rω by considering buoyancy frequency, Coriolis frequency, and topographic features. Compared to the GHP scaling with constant Rω, the turbulent dissipation rate values inferred from the GHP scaling with parameterized Rω are more accurate. In this manner, this study provides a reference for choosing optimal Rω when using the strain‐only GHP scaling to explore turbulent mixing in the open ocean. Plain Language Summary: Although diapycnal mixing is the smallest‐scale motion in the ocean, understanding its global distribution is crucial because it can modulate large‐scale circulation and hence influence climate change. However, field observations of ocean mixing are costly and therefore limited in a global sense. Given that internal wave breaking is the dominant energy source for ocean mixing, estimating internal wave‐driven mixing using easily available temperature and salinity data becomes an effective way to understand its spatiotemporal variability. Here, we modify an essential parameter in such an estimation method by changing the constant to a variable that depends on stratification, latitude, and topography. This modification significantly improves the accuracy of the method for estimating internal wave‐driven mixing, making it more consistent with observations compared to using a constant. Key Points: Shear‐to‐strain ratio varies significantly with spaceA parameterization of shear‐to‐strain ratio is proposed based on buoyancy frequency, Coriolis frequency, and topographic featuresBy using a parameterized shear‐to‐strain ratio, instead of a fixed value, we show the finescale parameterization performs better [ABSTRACT FROM AUTHOR]

We show that steric sea‐level varies with a period of 18.6 years along the western European coast. We hypothesize that this variation originates from the modulation of semidiurnal tides by the lunar nodal cycle and associated changes in ocean mixing. Accounting for the steric sea level changes in the upper 400 m of the ocean solves the discrepancy between the nodal cycle in mean sea level observed by tide gauges and the theoretical equilibrium nodal tide. Namely, by combining the equilibrium tide with the nodal modulation of steric sea level, we close the gap with the observations. This result supports earlier findings that the observed phase and amplitude of the 18.6‐year cycle do not always correspond to the equilibrium nodal tide. Plain Language Summary: The orbital position of the moon and the gravity pull it exerts on the earth varies with a period of 18.6 years. This cycle is called the lunar nodal cycle and it results in small variations of yearly averaged sea level (∼1–2 cm). Understanding this variability is important because it allows, for example, to quickly detect an acceleration in local sea‐level rise due to global warming. Here we show that the lunar nodal cycle also has an influence on the temperature and salinity in the surface 400m of the ocean. As a result, the ocean density changes and amplifies sea level variations along the western European coast. We make the hypothesis that since the lunar nodal cycle also influences the amplitude of the semidiurnal tides, and since those tides are known to be responsible for a large part of ocean mixing, a change in ocean mixing could be the cause of the ocean density variability that we observe. Key Points: Steric sea level changes are influenced by the 18.6‐year lunar nodal cycle along the western European coastThis influence could result from the modulation of semidiurnal tides by the lunar nodal cycle and the associated change in ocean mixingThis finding is a step toward resolving the long‐standing discrepancy between the theoretical long‐period nodal tide and observed signal [ABSTRACT FROM AUTHOR]

The main sources of redox gradients supporting high-productivity life in the Europan and other icy ocean world oceans were proposed to be photolytically derived oxidants, such as reactive oxygen species (ROS) from the icy shell, and reductants (Fe(II), S(-II), CH4, H2) from bottom waters reacting with a (ultra)mafic seafloor. Important roadblocks to maintaining life, however, are that the degree of ocean mixing to combine redox species is unknown, and ROS damage biomolecules. Here, we envisage a unique solution using an acid mine drainage (AMD)-filled pit lakes analog system for the Europan ocean, which previous models predicted to be acidic. We hypothesize that surface-generated ROS oxidize dissolved Fe(II) resulting in Fe(III) (hydr)oxide precipitates, that settle to the seafloor as "iron snow." The iron snow provides a respiratory substrate for anaerobic microorganisms ("breathing iron"), and limits harmful ROS exposure since they are now neutralized at the ice-water interface. Based on this scenario, we calculated Gibbs energies and maximal biomass productivities of various anaerobic metabolisms for a range of pH, temperatures, and H2 fluxes. Productivity by iron reducers was greater for most environmental conditions considered, whereas sulfate reducers and methanogens were more favored at high pH. Participation of Fe in the metabolic redox processes is largely neglected in most models of Europan biogeochemistry. Our model overcomes important conceptual roadblocks to life in icy ocean worlds and broadens the potential metabolic diversity, thus increasing total primary productivity, the diversity and volume of habitable environmental niches and, ultimately, the probability of biosignature detection. [ABSTRACT FROM AUTHOR]

Temporal and spatial variations in the ocean surface mixed layer are important for the climate and ecological systems. During 1980–2019, the Southern Indian Ocean (SIO) mixed layer depth (MLD) displays a basin‐wide shoaling trend that is absent in the other basins within 40°S–40°N. The SIO MLD shoaling is mostly prominent in austral winter with deep climatology MLD, substantially weakening the MLD seasonality. Moreover, the SIO MLD changes are primarily caused by a southward shift of the subtropical anticyclonic winds and hence ocean gyre, associated with a strengthening of the Southern Annular Mode, in recent decades for both winter and summer. However, the poleward‐shifted subtropical ocean circulation preferentially shoals the SIO MLD in winter when the meridional MLD gradient is sharp but not in summer when the gradient is flat. This highlights the distinct subtropical MLD response to meridional mitigation in winds due to different background oceanic conditions across seasons. Plain Language Summary: The ocean surface mixed layer (ML) is a well‐mixed layer with uniform physical and chemical properties and is key for the ocean in exchanging materials and energy with the atmosphere. The present study shows that during 1980–2019, the Southern Indian Ocean (SIO) ML depth (MLD) displays a basin‐wide decreasing trend, which is absent in the other basins within 40°S–40°N. The SIO MLD shoaling primarily appears in austral winter when the climatology ML is deep but is insignificant in summer, substantially weakening the MLD seasonality. The SIO MLD changes are primarily explained by the ocean dynamical adjustment driven by the surface zonal wind changes. Specifically, the strengthened Southern Annular Mode in recent decades drives southward shifts of the subtropical anticyclonic winds and ocean gyre year‐round. However, the poleward‐shifted ocean gyre preferentially decreases the SIO winter MLD as the meridional MLD gradient is sharp and thus efficiently reduces the deep ML water converging from the Southern Ocean into the SIO. In contrast, the SIO MLD displays negligible change in summer when its meridional gradient is flat. The results highlight that despite under nearly identical southward‐shifted subtropical winds, the winter and summer MLD responses are distinct due to different background oceanic conditions. Key Points: The seasonality of the Southern Indian Ocean surface mixed layer (ML) depth prominently weakens during 1980–2019The weakened seasonality mainly results from a pronounced winter ML shoalingThe southward shift of the subtropical ocean gyre driven by the strengthened Southern Annular Mode dominates the ML shoaling [ABSTRACT FROM AUTHOR]

Through a combination of laboratory experiments and theoretical models, we investigate the interaction of a mean upwelling through a closed basin with a vertical buoyancy flux. The fluid is mixed by a horizontally oscillating rake, which either traverses the whole basin or which oscillates just near one vertical boundary. We first review the steady state and demonstrate that, in both mixing regimes, the vertical density profile across the basin is controlled by the steady-state balance between the upward advective and diffusive fluxes of salinity as described by the classical model introduced by Munk (Deep-Sea Res. , vol. 13, issue 4, 1966, pp. 707–730). However, with boundary mixing, we show that both the upwelling and the buoyancy transport are localised to the mixing zone near the boundary, and the interior fluid is stagnant. We then develop a model to describe the transient evolution of the system if there is either a discrete increase or gradual decrease to the buoyancy flux. In the boundary mixing case, the change in the buoyancy flux at the lower boundary leads to a change in the buoyancy of the fluid in the boundary mixing region, and this induces a transient, buoyancy-driven flow in the boundary region in addition to the steady upwelling. In turn, an equal and opposite vertical flow develops in the interior, and this leads to a change in the density stratification of the interior fluid as the system adjusts to a new equilibrium. However, in our experiments, there is no vertical mixing in the interior and interior fluid may upwell or downwell dependent on the change to the buoyancy forcing. We discuss the implications of our results for the transport and mixing in the deep ocean, and the associated interpretation of field experiments. [ABSTRACT FROM AUTHOR]

The impacts of Stokes drift and sea-state-dependent Langmuir turbulence (LT) on the three-dimensional ocean response to a tropical cyclone in the Bohai Sea are studied through two-way coupled wave-current simulations. The Stokes drift is calculated from the simulated wave spectrum of the wave model, Simulating Waves Nearshore (SWAN), and then input to the Princeton Ocean Model with the generalized coordinate system (POMgcs) to represent the Langmuir effect. The Langmuir circulation is included in the vertical mixing of the ocean model by adding the Stokes drift to the shear of the vertical mean current and by including LT enhancements to the Mellor-Yamada 2.5 turbulent closure submodel. Simulations are assessed through the case study of Typhoon Masta in 2005 with a set of diagnostic experiments that incorporated different terms of Stokes production (SP) respectively. It is shown that with the consideration of SP, a deeper mixed layer, an enhanced vertical mixing coefficient KMS, and a more accurate representation of the vertical temperature distribution could be derived. Moreover, the effect of LT in elevating the turbulence mixing is stronger than that of Coriolis Stokes force (CSF) and Craik-Leibovich vortex force (CLVF). LT has a greater influence on the vertical mixing during typhoon than that in normal weather. [ABSTRACT FROM AUTHOR]

The response of the oceanic surface layer to atmospheric shallow convection is explored using realistic atmospheric large eddy simulations coupled with an oceanic 1D model with high vertical resolution. The effects of cloud shading, rain, and enhanced heat loss due to gust fronts on the edge of cold pools and their interactions are investigated in a case study of the Cooperative Indian Ocean Experiment on Intraseasonal Variability/Dynamics of the Madden–Julian Oscillation experiment in the tropical Indian Ocean, during a suppressed phase of the Madden–Julian Oscillation. Conditions of low surface wind and strong solar heating result in diurnal warming of the oceanic surface of 2 °C over a depth of 1 m. Analysis of specific periods covering the diurnal cycle shows the contrasting effects of cloud shading, rain, and turbulent heat fluxes under the cold pools on the sea temperature at the surface and below. On the one hand, decreasing the solar radiation (cloud shading) results in slight cooling extended horizontally and penetrating down to 1–2 m depth, depending on the time of the day. On the other hand, turbulent heat fluxes enhanced up to 300 W·$$ \cdotp $$m −2$$ {}^{-2} $$ by gusts and freshwater lenses due to rain act together and more locally. They isolate and strongly cool a thin inner layer at the surface, which eventually destabilizes the surface layer and propagates the cooling downward. The exact relative part and efficiency of these processes depend on the time evolution of the thermal stratification and vertical turbulent mixing in the oceanic upper layer. Surface cooling of up to −$$ - $$0.5 °C may occur in a few tens of minutes and last for several hours, significantly mitigating the effects of diurnal warming over large extents. [ABSTRACT FROM AUTHOR]