Ocean Mixing Research Articles

Generation Mechanisms of SST Anomalies Associated with the Canonical El Niño Focusing on Vertical Mixing

Abstract The mean vertical advection of anomalous vertical temperature gradient is considered the dominant generation mechanism of positive sea surface temperature (SST) anomalies associated with the canonical El Niño. However, most past studies had a residual term in their heat budget analysis and/or did not quantify the role of vertical mixing even though active vertical turbulent mixing in the upper ocean is observed in the eastern equatorial Pacific. To quantitatively assess the importance of vertical mixing, a mixed layer heat budget analysis is performed using a hindcast simulation forced by daily mean atmospheric reanalysis data. It is found that when the mixed layer depth is defined as the depth at which potential density increases by 0.125 kg m−3 from the sea surface, the development of positive SST anomalies is predominantly governed by reductions in the cooling by vertical mixing, and their magnitude is much larger than those by vertical advection. The anomalous warming by vertical mixing may be partly explained by an anomalous deepening of the thermocline that leads to a decrease in the vertical temperature gradient, giving rise to suppression of the climatological cooling by vertical mixing. Also, an anomalously thick mixed layer reduces sensitivity to cooling by the mean vertical mixing and contributes to the anomalous SST warming. On the other hand, the dominant negative feedbacks are attributed to both anomalous surface heat loss and anomalous deepening of the mixed layer that weakens warming by the mean surface heat flux.

Linking Ocean Mixing and Overturning Circulation

Linking Ocean Mixing and Overturning Circulation

Spatial and Temporal Variability of Turbulent Mixing in the Deep Northwestern Pacific

AbstractSmall‐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.

Arctic Freshwater Sources and Ocean Mixing Relationships Revealed With Seawater Isotopic Tracing

AbstractThe 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.

Controls on Dissolved Barium and Radium‐226 Distributions in the Pacific Ocean Along GEOTRACES GP15

AbstractRadium‐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.

Parameterization of Shear‐To‐Strain Ratio Used in Finescale Parameterization

AbstractThe 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.

Dynamic causes contribute to the increasing trend of red tides in the east China sea during 2020–2022

Red tide is a marine phenomenon caused by the excessive growth of microscopic algae in the ocean. This study aims to analyze the development trends of red tides in the past 20 years and the dynamic external causes that induce red tides based on existing satellite remote sensing and numerical simulation data. And the changes in dominant species of red tides in different seasons are analyzed. The results show significant temperature fluctuations within the week before the red tide occurs, with an average increase of 1.42 °C. In contrast, the change in salinity is relatively small. Meanwhile, ocean fronts are areas in the ocean where different water masses meet and form boundaries. The average strength of ocean fronts increased by 3.7%, indicating enhanced ocean mixing over a short period of time. Under the combined influence of these factors, the probability of a red tide outbreak in the East China Sea increases rapidly. Therefore, this study has important reference value for further research on the causes of red tides and their response to ocean dynamic changes.

Boundary mixing. Part 2. The impact of ventilation

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.

The Effect of the 18.6‐Year Lunar Nodal Cycle on Steric Sea Level Changes

AbstractWe 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.

Bioenergetics of iron snow fueling life on Europa

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.

Seasonal variability and long-term winter shoaling of the upper mixed layer in the southern Baltic Sea

The upper ocean mixed layer plays a crucial role in regulating the exchange between the ocean and the atmosphere. The Mixed Layer Depth (MLD) is a key parameter affecting the air-sea exchanges of momentum and heat and determining the upper ocean temperature. Numerous previous studies have investigated MLD variability in the global ocean or regional seas but no such studies were carried in the Baltic Sea. In this study, we present the first observational assessment of the MLD and its properties in the Southern Baltic Sea including quantification of its seasonal and long-term changes and identification of the multi-year winter shoaling. We calculated monthly maps of MLD in the southern Baltic Sea using a large number of historical CTD profiles collected in 1995–2021 from a combination of different data sets. To test the robustness of the results we compared the MLDs calculated using different threshold methods. Throughout the southern Baltic Sea, across its three basins, a distinct seasonality is evident in the MLD, with values varying from 12 m in July to 60 m in December–March. During winter the water column is well mixed down to the upper halocline depth and the MLD reaches about 45 m in the Bornholm Basin, 50 m in the Slupsk Furrow, and 60 m in the Gdansk Basin. The observed global warming and decadal changes in the salty inflows from the North Sea to the Baltic have had an impact on stratification by increasing water densities in the intermediate and deep layers. Consequently, density gradients have strengthened with depth while the upper ocean mixing has weakened during the winter season. The results reveal a significant winter shoaling of the mixed layer by 4 m per decade, driven by the increased stratification due to rising temperatures and salinity. These changes could have significant impacts on the dynamics and productivity of marine ecosystems.

Synthetic Lagrangian turbulence by generative diffusion models

Lagrangian turbulence lies at the core of numerous applied and fundamental problems related to the physics of dispersion and mixing in engineering, biofluids, the atmosphere, oceans and astrophysics. Despite exceptional theoretical, numerical and experimental efforts conducted over the past 30 years, no existing models are capable of faithfully reproducing statistical and topological properties exhibited by particle trajectories in turbulence. We propose a machine learning approach, based on a state-of-the-art diffusion model, to generate single-particle trajectories in three-dimensional turbulence at high Reynolds numbers, thereby bypassing the need for direct numerical simulations or experiments to obtain reliable Lagrangian data. Our model demonstrates the ability to reproduce most statistical benchmarks across time scales, including the fat-tail distribution for velocity increments, the anomalous power law and the increased intermittency around the dissipative scale. Slight deviations are observed below the dissipative scale, particularly in the acceleration and flatness statistics. Surprisingly, the model exhibits strong generalizability for extreme events, producing events of higher intensity and rarity that still match the realistic statistics. This paves the way for producing synthetic high-quality datasets for pretraining various downstream applications of Lagrangian turbulence.

Anthropocene North Western Pacific Oceanography Recorded as Seasonal‐Resolution Radiocarbon in Coral From Kikai Island, Japan

AbstractRadiocarbon (14C) in corals can be used as a relatively high‐sensitivity indicator of vertical and horizontal advection of water masses, which contributes to the understanding of ocean circulation. In this study, we reconstruct Kuroshio and Ryukyu current transport with a seasonal resolution Δ14C record spanning 1947–2009. This record covers the beginning of the atomic era and was obtained from a coral on Kikai Island in the south of Japan. The Kikai Δ14C curve features a newly discovered Δ14C spike in July 1955, a rapid increase after 1962, and a steady decrease after 1980. The spike in 1955 may directly reflect ocean current transport. The lack of periodicity in the Δ14C record suggests the existence of mesoscale eddies and the complexity of Kuroshio and Ryukyu current transport. In addition, comparing the high‐resolution Δ14C of Kikai and Ishigaki islands, both situated along the path of the Kuroshio, reveals the influence of Pacific Decadal Oscillation and El Niño‐Southern Oscillation on the Kuroshio and Ryukyu currents. This suggests that seasonally resolved Δ14C in corals along an ocean current can produce a long‐term record of ocean mixing that responds to climate variability.

Modulation of the northward penetration of Antarctica intermediate waters into the eastern equatorial Indian Ocean under glacial and interglacial conditions

The Indo-Pacific warm pool is the warmest and most dynamic ocean–atmosphere-climate system on Earth and was subject to significant climate changes during the Pleistocene glacial-interglacial transitions. This has been shown to significantly affected the strength of surface waters that redistribute heat from the tropics to the southern part of the Indian Ocean. Here we investigate the response of the oceanic circulation at intermediate depth (1200 m) of the eastern equatorial Indian Ocean (EEIO) with neodymium (Nd) isotopes in the context of the climatic oscillation of the last 500 ka. The most striking feature of our new dataset is the seesaw Nd record that mimics glacial-interglacial cycles. While the interglacial periods are characterized by a higher contribution of the less radiogenic neodymium (~ − 7εNd) Antarctic Intermediate Water (AAIW), the glacial periods are characterized by more radiogenic water mass of Pacific origin (~ − 5εNd). To explain the increase in the εNd signature toward a more radiogenic signature as the Indo-Pacific connection is reduced under the low sea level of the glacial periods, we show that under global cooling, the AAIW advances northward into the tropics, which is a consequence of the general slowdown of the thermohaline circulation. Therefore, oceanic mixing at intermediate depth in the eastern tropical Indian intermediate water is modulated by the production rate of the AAIW in the Southern Ocean. Our study provides new evidence for the role that changes in the deep oceanic conditions play in amplifying externally forced climate changes that ultimately lead to drier/moister atmospheric conditions and weaker/stronger monsoons during glacial/interglacial periods over eastern tropical Indian Ocean.

Studies on parameterizations of turbulent mixing in the deep ocean

Studies on parameterizations of turbulent mixing in the deep ocean

Long-term trends and extreme events of marine heatwaves in the Eastern China Marginal Seas during summer

Marine heatwaves (MHWs) are a type of widespread, persistent, and extreme marine warming event that can cause serious harm to the global marine ecology and economy. This study provides a systematic analysis of the long-term trends of MHWs in the Eastern China Marginal Seas (ECMS) during summer spanning from 1982 to 2022, and occurrence mechanisms of extreme MHW events. The findings show that in the context of global warming, the frequency of summer MHWs in the ECMS has increased across most regions, with a higher rate along the coast of China. Areas exhibiting a rapid surge in duration predominantly reside in the southern Yellow Sea (SYS) and southern East China Sea (ECS, south of 28°N). In contrast, the long-term trends of mean and maximum intensities exhibit both increases and decreases: Rising trends primarily occur in the Bohai Sea (BS) and Yellow Sea (YS), whereas descending trends are detected in the northern ECS (north of 28°N). Influenced jointly by duration and mean intensity, cumulative intensity (CumInt) exhibits a notable positive growth off the Yangtze River Estuary, in the SYS and southern ECS. By employing the empirical orthogonal function, the spatio-temporal features of the first two modes of CumInt and their correlation with summer mean sea surface temperature (SST) and SST variance are further examined. The first mode of CumInt displays a positive anomalous pattern throughout the ECMS, with notable upward trend in the corresponding time series, and the rising trend is primarily influenced by summer mean SST warming. Moreover, both of the first two modes show notable interannual variability. Extreme MHW events in the SYS in 2016 and 2018 are examined using the mixed layer temperature equation. The results suggest that these extreme MHW events originate primarily from anomalous atmospheric forcing and oceanic vertical mixing. These processes involve an anomalous high-pressure system over the SYS splitting from the western Pacific subtropical high, augmented atmospheric stability, diminished wind speeds, intensified solar radiation, and reduced oceanic mixing, thereby leading to the accumulation of more heat near the sea surface and forming extreme MHW events.

Influence of planetary rotation on metal-silicate mixing and equilibration in a magma ocean

At a late stage of its accretion, the Earth experienced high-energy planetary impacts. Following each collision, the metal core of the impactor sank into molten silicate magma oceans. The efficiency of chemical equilibration between these silicates and the metal core controlled the composition of the Earth interior and left a signature on geochemical and isotopic data. These data constrain the timing, pressure and temperature of Earth formation, but their interpretation strongly depends on the efficiency of metal-silicate mixing and equilibration. We investigate the role of planetary rotation on the dynamics of the sinking metal and on its chemical equilibration using laboratory experiments of particle clouds settling in a rotating fluid. Our clouds initially sink as spherical turbulent thermals, but after a critical depth, rotation becomes important and they transition to a vortical columnar flow aligned with the rotation axis. Applied to Earth formation, our results predict that rotation strongly affects the fall of metal in the magma ocean for impactors smaller than 459 km in radius on a proto-Earth that rotates twice faster than today. On a proto-Earth spinning 5 times faster than today, rotation is important for any impactor smaller than the Earth itself. In contrast with a thermal that grows in all directions, the vortical column grows vertically but keeps a constant horizontal extent. The slower dilution in vortical columns reduces chemical equilibration compared to previous estimates that neglect planetary rotation. We find that rotation significantly affects the degree of equilibration for highly siderophile elements with partition coefficients larger than 103. In this case, for a planet spinning twice faster than today, the degree of equilibration decreases by up to a factor 2 compared to previous estimates that neglect the effect of rotation. Finally, the ultimate fate of iron drops is to be detrained from the vortical column as an iron rain, reconciling the traditional iron rain scenario with the model of turbulent thermal.

Interannual Variability of Subpolar Mode Water in the Subpolar North Atlantic

AbstractSubpolar Mode Water (SPMW) is an important water mass originating in the eastern North Atlantic. Its formation, subject to modification through oceanic interior mixing, can directly influence the volume of water contributing to the Atlantic meridional overturning circulation. Utilizing observation‐based data sets spanning from 1993 to 2018, we estimated the formation rates and volume of SPMW within isopycnal layers and examined its temporal variability. Two complementary approaches were used to estimate the formation rate: a thermodynamic approach focusing on the air‐sea interactions and a kinematic approach involving volume transport from the mixed layer to the ocean's interior, including the entrainment/detrainment of the mixed layer itself. This is the first time that thermodynamic and kinematic approaches are applied to observation‐based data in the North Atlantic. Our results suggest a substantial role of diapycnal mixing in diluting the dense waters formed by air‐sea fluxes toward the range of SPMW densities. The study reveals a complex interplay of processes, with entrainment being the primary driver of subduction/obduction rates, while advection contributes to the overall small‐scale dynamics. Variations in the volume and location of SPMW formation are observed from year to year. Notably, when SPMW forms extensively in lighter isopycnal layers, the volume occupied by denser isopycnals decreases and vice versa. We attributed this compensation effect to a propagation signal, where formation in the lightest isopycnal bins influences the formation in denser isopycnal bins with a delay of a few years, emphasizing the circulation's role in shaping the SPMW distribution.

Eddy–Internal Wave Interactions and Their Contribution to Cross-Scale Energy Fluxes: A Case Study in the California Current

Abstract Oceanic mixing, mostly driven by the breaking of internal waves at small scales in the ocean interior, is of major importance for ocean circulation and the ocean response to future climate scenarios. Understanding how internal waves transfer their energy to smaller scales from their generation to their dissipation is therefore an important step for improving the representation of ocean mixing in climate models. In this study, the processes leading to cross-scale energy fluxes in the internal wave field are quantified using an original decomposition approach in a realistic numerical simulation of the California Current. We quantify the relative contribution of eddy–internal wave interactions and wave–wave interactions to these fluxes and show that eddy–internal wave interactions are more efficient than wave–wave interactions in the formation of the internal wave continuum spectrum. Carrying out twin numerical simulations, where we successively activate or deactivate one of the main internal wave forcing, we also show that eddy–near-inertial internal wave interactions are more efficient in the cross-scale energy transfer than eddy–tidal internal wave interactions. This results in the dissipation being dominated by the near-inertial internal waves over tidal internal waves. A companion study focuses on the role of stimulated cascade on the energy and enstrophy fluxes.

Enriched Regions of 228Ra Along the U.S. GEOTRACES Pacific Meridional Transect (GP15)

AbstractThe half‐life of 228Ra (5.7 years) aligns well with near‐surface and near‐bottom ocean mixing timescales. Because 228Ra is sourced from sediments, regions of enhanced activity represent water that has recently interacted with sediments on the continental margin or seabed. The GP15 meridional transect from Alaska to Tahiti along152°W encountered several regions in the upper ocean where 228Ra was enriched. These enrichments follow surface and subsurface ocean current patterns and pair with earlier measurements of 228Ra and transient radionuclides to reveal the origins of these enriched regions. An enriched region at Alaska margin stations 1–3 was sourced locally but did not extend to the Alaskan trench at station 4. A large shallow region between 47° and 32°N. was sourced from the west by the North Pacific Current; another shallow enriched region between 11° and 5° N was also sourced from the west by the North Equatorial Countercurrent. Subsurface enrichments (100–400 m) between 18 and 47°N were associated with Central Mode Water and North Pacific Intermediate Water. The 228Ra activities in the upper Pacific were six times lower than activities in the Atlantic. In deep waters the primary enrichment was 27°–47°N. Two stations (32° and 37°N) were especially enriched, having near‐bottom inventories several times greater than other stations. With these two exceptions the remaining Pacific stations exhibited averaged inventories lower than those in the Atlantic. There was one region of enriched 223Ra (half‐life = 11 days) above the Puna Ridge near Hawaii.