Ocean Processes Research Articles

Metrics for quantifying the efficiency of atmospheric CO2 reduction by marine carbon dioxide removal (mCDR)

Abstract Marine carbon dioxide removal (mCDR) is gaining interest as a tool to meet global climate goals. Because the response of the ocean–atmosphere system to mCDR takes years to centuries, modeling is required to assess the impact of mCDR on atmospheric CO2 reduction. Here, we use a coupled ocean–atmosphere model to quantify the atmospheric CO2 reduction in response to a CDR perturbation. We define two metrics to characterize the atmospheric CO2 response to both instantaneous ocean alkalinity enhancement (OAE) and direct air capture (DAC): the cumulative additionality (α) measures the reduction in atmospheric CO2 relative to the magnitude of the CDR perturbation, while the relative efficiency (ϵ) quantifies the cumulative additionality of mCDR relative to that of DAC. For DAC, α is 100% immediately following CDR deployment, but declines to roughly 50% by 100 years post-deployment as the ocean degasses CO2 in response to the removal of carbon from the atmosphere. For instantaneous OAE, α is zero initially and reaches a maximum of 40%–90% several years to decades later, depending on regional CO2 equilibration rates and ocean circulation processes. The global mean ϵ approaches 100% after 40 years, showing that instantaneous OAE is nearly as effective as DAC after several decades. However, there are significant geographic variations, with ϵ approaching 100% most rapidly in the low latitudes while ϵ stays well under 100% for decades to centuries near deep and intermediate water formation sites. These metrics provide a quantitative framework for evaluating sequestration timescales and carbon market valuation that can be applied to any mCDR strategy.

A Digital Twin of the Trondheim Fjord for Environmental Monitoring—A Pilot Case

Digital Twins of the Ocean (DTO) are a rapidly emerging topic that has attracted significant interest from scientists in recent years. The initiative, strongly driven by the EU, aims to create a digital replica of the ocean to better understand and manage marine environments. The Iliad project, funded under the EU Green Deal call, is developing a framework to support multiple interoperable DTO using a federated systems-of-systems approach across various fields of applications and ocean areas, called pilots. This paper presents the results of a Water Quality DTO pilot located in the Trondheim fjord in Norway. This paper details the building blocks of DTO, specific to this environmental monitoring pilot. A crucial aspect of any DTO is data, which can be sourced internally, externally, or through a hybrid approach utilizing both. To realistically twin ocean processes, the Water Quality pilot acquires data from both surface and benthic observatories, as well as from mobile sensor platforms for on-demand data collection. Data ingested into an InfluxDB are made available to users via an API or an interface for interacting with the DTO and setting up alerts or events to support ’what-if’ scenarios. Grafana, an interactive visualization application, is used to visualize and interact with not only time-series data but also more complex data such as video streams, maps, and embedded applications. An additional visualization approach leverages game technology based on Unity and Cesium, utilizing their advanced rendering capabilities and physical computations to integrate and dynamically render real-time data from the pilot and diverse sources. This paper includes two case studies that illustrate the use of particle sensors to detect microplastics and monitor algae blooms in the fjord. Numerical models for particle fate and transport, OpenDrift and DREAM, are used to forecast the evolution of these events, simulating the distribution of observed plankton and microplastics during the forecasting period.

Dynamics of subsurface chlorophyll maxima in the northern Indian Ocean

Subsurface chlorophyll maxima (SCM) significantly contributes to oceanic primary productivity, emphasizing the need to study its dynamics and governing mechanisms. We used datasets from various platforms to investigate relationships between the SCM characteristics (SCM depth (ZSCM), SCM magnitude (Chlmax), SCM thickness (TSCM)) and environmental variables modulated by various physical processes in the Northern Indian Ocean (NIO). In the Arabian Sea (western NIO), seasonal processes like convective mixing and upwelling, primarily regulated the SCM characteristics. In the Bay of Bengal (eastern NIO), SCM characteristics were jointly influenced by fresh water influx, barrier layer formation, presence of eddies, and the propagation of Kelvin and Rossby waves. Any changes in these oceanic processes, potentially driven by climate change, could therefore impact oceanic primary production. Additionally, a positive association obtained between Chlmax and downward CO2 flux, while a shallower ZSCM, associated with higher concentrations of DMS, indicated SCM's role in regulating atmospheric gases.

Local and Remote Atmosphere‐Ocean Coupling During Extreme Warming Events Impacting Subsurface Ocean Temperature in an Antarctic Embayment

AbstractCoastal ocean temperatures can respond to different atmospheric and oceanic processes at local spatial scales or through remote teleconnections. This study focused on subsurface ocean temperatures (subT) at 10 m depth in Maxwell Bay, northern Antarctic Peninsula from February 2017 and January 2022. It investigated extreme warming events during austral summers and their interaction with atmospheric and oceanic conditions regionally and locally. The analysis identified active and extreme Marine Heat Waves (MHWs) in March 2017 and January‐February 2020 associated with a significantly negative Southern Annular Mode index observed 3–4 months before the temperature increase. In March 2017, temperatures exceeded the climatological mean by over 1°C. This anomaly was linked to a strengthened Amundsen Sea Low and a blocking anticyclone moving between the Scotia Sea and the South‐West Atlantic Ocean that deflected westerly winds and facilitated the anomalous transport of warmer northern air masses to the AP. In January‐February 2020, the highest recorded subT was observed (2.97°C), although air‐sea heat fluxes did not show a similar pattern. In February 2020, one of the most intense atmospheric heatwaves ever recorded in West Antarctica was observed. This heatwave corresponded with maximum subT and positive sea surface temperature anomalies extending throughout the western region of the Southern Ocean, related to an extremely negative Southern Annular Mode. This study provides valuable insights into the impact of strong MHWs, a phenomenon that has been less documented in Antarctic coastal regions.

Shoaling Internal Tides Propagating From a Shallow Ridge Modulated by the Kuroshio

AbstractThis study presents remotely generated internal tides that propagate into shallow regions and are modulated by background flows using numerical simulations and field observation data. Numerical results indicate that strongly enhanced semi‐diurnal (∼M2) internal tide energy is generated over a shallow ridge, the Izu‐Ogasawara Ridge. The generated internal tides propagates toward the Kuroshio upstream and shoal toward shallow regions near Cape Shiono‐Misaki. The intensified internal tide energy flux toward the cape is explained by two mechanisms: (a) intensified internal tide generation over the ridge due to the interaction between tides and the Kuroshio and (b) wave energy convergence along the Kuroshio due to wave refractions. A 13‐year field data set obtained from the cape was used to investigate shoaling internal tides influenced by the Kuroshio. The observed results reveal a significant positive correlation between the kinetic energy of semi‐diurnal internal tides and the background flows caused by the Kuroshio, which evidently supports the intensified internal tides attributed to the interaction between background flows and tides, as proposed by recent studies. The intensity of shoaling internal tides is also largely influenced by the path of the Kuroshio and seasonal effects. The magnitude of shoaling internal tides is clearly weaken as the Kuroshio meander occurs. Shoaling internal tides modulated by the Kuroshio can provide new insights into energy transport and mixing processes in coastal oceans.

Deep learning for the super resolution of Mediterranean sea surface temperature fields

Abstract. Sea surface temperature (SST) is one of the essential variables of the Earth's climate system. Being at the air–sea interface, SST modulates heat fluxes in and out of the ocean, provides insight into several upper and interior ocean dynamical processes, and is a fundamental indicator of climate variability potentially impacting the health of marine ecosystems. Its accurate estimation and regular monitoring from space is therefore crucial. However, even if satellite infrared/microwave measurements provide much better coverage than what is achievable from in situ platforms, they cannot sense the sea surface under cloudy and rainy conditions. Large gaps are present even in merged multi-sensor satellite products, and different statistical strategies, mostly based on optimal interpolation (OI) algorithms, have thus been proposed to obtain gap-free (L4) images. These techniques, however, filter out the signals below the space–time decorrelation scales considered, significantly smoothing most of the small mesoscale and submesoscale features. Here, deep learning models, originally designed for single-image super resolution (SR), are applied to enhance the effective resolution of SST products and the accuracy of SST gradients. SR schemes include a set of computer vision techniques leveraging convolutional neural networks to retrieve high-resolution data from low-resolution images. A dilated convolutional multi-scale learning network, which includes an adaptive residual strategy and implements a channel attention mechanism, is used to reconstruct features in SST data at 1/100° spatial resolution starting from 1/16° data over the Mediterranean Sea. The application of this technique shows an improvement in the high-resolution reconstruction, capturing small-scale features and providing a root-mean-squared-difference improvement of 0.02 °C with respect to the L3 ground-truth data.

Linkage of El Niño-Southern Oscillation to astronomic forcing

Abstract Deciphering the causes of El Niño-Southern Oscillation (ENSO) represents one of the greatest scientific and societal challenges because ENSO impacts people’s safety, food, water, health, and economy. Traditionally, ENSO has been considered a phenomenon that is mostly influenced by the interactions between oceanic and atmospheric processes, i.e. by the internal variability in Earth’s ocean-atmosphere system. However, dendrochronological records of climate indices, ENSO among them, have statistically significant variance at periodicities related to solar activity (sunspots) and lunar precessions. Other studies indicate a modulation by the lunar nodal cycle of ocean’s mixing and its implications on water and air temperatures, steric sea levels, coastal flooding, rain, river discharge and heat transport. Those findings suggest that astronomic forces may influence ENSO. In fact, this study shows that a fit to well-established periodicities from lunar precessions, solar activity and their interactions explains 91% of the variance of an ENSO index smoothed at 5 years, and 67% of a 3 yr filtered ENSO index. Provided that the future Earth’s ocean-atmospheric system remains roughly within historical bounds under a changing climate in the next 2–3 decades, one can venture a projection into the future of 3- and 5 yr smoothed ENSO variability. Such projection is given here and might help in preparations for adaptation and mitigation measures caused by ENSO-related coastal hazards.

The late Permian–Early Triassic island-Arc granitic rocks in the western Qaidam Basin, Qinghai–Tibetan Plateau: constraints on subduction processes of Paleo-Tethys Ocean

The late Permian–Early Triassic island-Arc granitic rocks in the western Qaidam Basin, Qinghai–Tibetan Plateau: constraints on subduction processes of Paleo-Tethys Ocean

An Assessment of Convergence of Climate Reanalyses from Two Coupled Data Assimilation Systems with Identical High-Efficiency Filtering

Abstract Coupled data assimilation (CDA) uses coupled model dynamics and physics to extract observational information from measured data in multiple Earth system domains to reconstruct historical states of the Earth system, forming a reanalysis of climate variability. Due to imperfect numerical schemes in modeling dynamics and physics, models are usually biased from the real world. Such model bias is a critical obstacle in the reconstruction of historical variability by combining model and observations and, to some degree, causes divergence of CDA results because of individual model behavior in each system. Here, based on a multitime-scale high-efficiency filtering algorithm which includes a deep-ocean bias relaxing scheme, we first develop a high-efficiency online CDA system with the Community Earth System Model (CESM-MSHea-CDA). Then, together with the other previously established CDA system within the Coupled Model, version 2.1 (CM2-MSHea-CDA), that is developed by the Geophysical Fluid Dynamics Laboratory, we conduct climate reanalysis for the past 4 decades (1978–2018). Evaluations show that due to improved representation for multiscale background statistics and effective deep-ocean model bias relaxing, both CDA systems produce convergent estimation of variability for major climate signals such as variability of basin-scale ocean heat content, El Niño–Southern Oscillation (ENSO), and Pacific decadal oscillation (PDO). Particularly, both CDA systems generate similar time mean of global and Atlantic meridional overturning circulations that converge to the geostrophic velocity estimate from climatological temperature and salinity data. The CDA-estimated mass transport at typical measurement sections is mostly consistent with the observations. Significance Statement A coupled climate model simulates the interactions of the atmosphere, ocean, sea ice, and land processes and derives change and variations of the climate system by combining these components. Due to imperfect numerical schemes, a model usually has systematic biases to the real world. Coupled data assimilation (CDA) combines coupled model dynamics and physics with atmospheric and oceanic observations to reconstruct the historical states of Earth’s climate system, forming the analysis of climate change and variability. However, because of the existence of model errors and individual data assimilation algorithm behaviors, currently, the ocean variability of reanalysis produced by different numerical systems is basically divergent. Here, we first apply an identical coupled data assimilation algorithm and effective model error relaxing scheme to two widely used IPCC coupled models to establish two CDA systems. Then, we use both CDA systems to assimilate historical atmospheric and oceanic observations to form two sets of climate reanalyses for the past 4 decades. We found, by some degree, that the results of these multimodel coupled climate reanalyses converged to what historical states ought to be in observational estimation. These convergent climate reanalyses are expected to significantly increase our understanding of CDA and climate assessment.

An atmosphere–ocean coupled model for simulating physical and biogeochemical state of north Indian Ocean: Customisation and validation

An atmosphere-ocean biogeochemistry coupled model is configured, customised and validated to decipher the role of atmospheric and oceanic processes in the north Indian Ocean (NIO, 3–30° N, 40–100° E), Arabian Sea (AS, 4–25° N, 50–75° E) and Bay of Bengal (BoB, 4–25° N, 76–100° E). The validation of model results with measurements shows that Sea Surface Temperature is better simulated by the coupled model than the standalone ocean model, with an average bias of ±0.2 °C in NIO. The simulated Sea Surface Salinity has a smaller bias with respect to reanalysis in AS (0.2 psu), but slightly higher in BoB (0.5 psu). The vertical distribution of temperature and salinity is also better represented in the coupled model. The atmospheric forcing, such as Long Wave Radiation, wind stress, and net surface heat and salt fluxes are better simulated by the coupled model with comparable variability and seasonality to the reanalysis data. The winter and summer Chlorophyll-a (Chl-a) blooms in AS are also well reproduced by the coupled model. The model also well simulates NO3 within top 100 m in AS and BoB. The nitrate budget analysis indicates that vertical advection and entrainment play big roles in governing the overall budget in upwelling regions of NIO. Nitrate uptake by phytoplankton is the dominant biological process, except in the northern AS, where denitrification is dominant. New production contributes to Net Primary Productivity (NPP) in most upwelling regions in AS except in the northern AS, where regenerated production is higher than the new production. In BoB, the regenerated production dominates except in East Coast of India during monsoon season. Therefore, our study provides new insights on the capability of the coupled ocean-atmosphere models in simulating the physical and biogeochemical processes, and air-sea interactions in NIO.

Satellite and in situ measurements of water vapour in the Brazil‐Malvinas Confluence region

AbstractThis decade‐long study from the INTERCONF programme addresses the data gap on ocean dynamics in the Southern Hemisphere. Focusing on the Brazil‐Malvinas Confluence (BMC), we investigated the effect of temperature differences between the warm Brazil Current (BC) and the cold Malvinas Current (MC) on water vapour in the marine atmospheric boundary layer (MABL). Our results show a clear distinction: warmer BMC waters have 32% more water vapour (2 kg·m−2 on average) compared to the MC. This highlights the direct link between ocean temperatures and atmospheric processes. Analysis of radiosonde data alongside satellite measurements showed better agreement over cooler waters with lower water vapour, leading to lower variability. This suggests less atmospheric turbulence and improved data compatibility, especially for satellite retrievals such as AIRS. Comparisons between reanalysis data (CFSR), satellite sounders (AIRS) and radiosondes (RS) showed consistent air temperature profiles, with average errors within the 10% threshold for satellite measurements. While capturing humidity variations remains a challenge, especially at high concentrations (indicated by higher mean squared error values), our study highlights the reliability of satellite data, particularly over the cold BMC region. This research highlights the importance of studying the interactions between ocean fronts and atmospheric phenomena for a complete picture of Southern Hemisphere ocean dynamics. It offers valuable insights for scientists across disciplines, providing a broad perspective on the results and their significance in different contexts.

Ocean processes south of the Drygalski Ice Tongue, western Ross Sea

We describe the first year-long hydrographic mooring timeseries from a location just to the south of the Drygalski Ice Tongue – the ice margin that forms the southern boundary of the Terra Nova Bay Polynya in the western Ross Sea. The region is where any northward flowing component of the Victoria Land Coastal Current encounters the ice tongue and supports an occasional polynya. The hydrographic mooring was deployed nearby Geikie Inlet from February 2017 through to March 2018, and was coupled with several contemporaneous oceanographic moorings to the north of the Drygalski Ice Tongue. This provides data with which to examine the water column dynamics in the context of local circulation and interaction with the ice tongue. The Terra Nova Bay region is subject to strong katabatic winds, however the polynya to the south of the Drygalski Ice Tongue operates at different times through the annual cycle when compared to the Terra Nova Bay Polynya to the north, as the sea ice in the south-side region is far more constrained in its motion yet, temperature and salinity are broadly consistent north and south of the ice tongue. Sub-surface Ice Shelf Water is observed south of the ice tongue. Transients in near-bed temperature and salinity are observed on both sides of the ice tongue, albeit with the northside leading by ∼8–9 days. Notably, the temperature transient precedes that of salinity by around 40 days. This suggests that, at this near-coastal position, the circulation beneath the ice tongue is primarily southward.

Sediment remobilization over subaqueous sand waves: In-situ observation in the northern South China Sea

High-resolution tripod observations were conducted in a subaqueous sand-wave field in the northern South China Sea. The objective was to gain insight into the dynamic mechanism by which sandy sediments are remobilized by oceanic dynamic processes, in particular internal solitary waves and internal tides. Daily-recurring high suspended sediment concentrations were observed, likely due to lateral transport of sediments resuspended by critical reflection of diurnal internal tides at the shelf break. The passage of extreme internal solitary waves results in resuspension of sediments from the seabed and ejection of these sediments out of the boundary layer. The resuspension of sediments from the seabed is controlled by current-induced strong bed shear stress, while the ejection of sediments out of the bottom boundary layer is driven by the compression and subsequent expansion of the mixing layer, together with the developments of the separation bubbles behind the internal solitary waves. These findings provide new insights into understanding the dynamic mechanism of sediment remobilization by internal solitary waves. Moreover, they contribute to our understanding of the migration of subaqueous sand waves.

Role of atmospheric and oceanic processes on interannual summertime (2016–2017) decrease of sea ice in the Antarctic regions of the Southern Ocean

In this article, the interannual variability of sea ice in the Antarctic sea ice regions between 2013–2018 is studied using a global ocean sea ice coupled model and satellite observation. The numerical model reasonably well simulates satellite observed interannual variability of sea ice concentration (SIC) and sea surface temperature (SST) in the Antarctic regions of Southern Ocean during all four austral seasons; summer (December–February), autumn (March–May), winter (June–August), and spring (September–November).A comparison of satellite and model shows that, during last two decades between 2001–2020, summertime of 2016–2017 had the lowest (highest) SIC (SST) across the Antarctic sea ice regions. Mixed layer heat budget analysis has been performed to comprehend how thermodynamic processes affect changes in SIC and SST in the Antarctic sea ice regions. The strong positive net atmospheric heat flux and the negative ocean vertical entrainment during summertime of 2016–2017 resulted in increased SST compared to other years, which lead to decreased SIC during above years. Also, loss of sea ice during summertime of 2016–2017 in the Antarctic sea ice regions are linked with significant decrease of wind stress magnitude and increase of wind stress curl.

A Study of the Vortex Filament Pool Left by a Super Typhoon

AbstractEnhancing typhoon forecasts hinges on a deeper understanding of the upper‐ocean’s response and feedback mechanisms to typhoons. Presently, our knowledge of typhoon‐ocean interactions is largely derived from low‐resolution numerical simulations (often >10 km) and limited observations, which inadequately capture submesoscale processes (SPs) in the ocean. Connecting extreme typhoons to upper‐ocean SPs remains a challenge. This study reveals the formation of a distinctive vortex filament pool (VFP) in high‐resolution (∼1.2 km) numerical experiments. The experiments show that under specific conditions, typhoons can generate this visually striking phenomenon, displaying SP dynamics and kinematics typical of the upper ocean, with Rossby numbers and the nondimensional strain and divergence rates exceeding 2. The VFP formation is mainly driven by strain‐induced frontogenesis linked to the flow generated by Typhoon Nangka after a major turn. Initially, the pool consists of many near‐parallel filaments, but processes such as merging, stretching, and destabilization subsequently occur lead to numerous small vortices with a mean radius of ∼13 km. While the high‐resolution numerical experiments highlight phenomena requiring observational validation, they suggest the presence of natural processes previously undetectable with low‐resolution models and limited observations. This study underscores the need for enhanced observations and numerical models to better understand refined ocean dynamical processes.

Drivers of the mean biases of the tropical atmospheric circulation in a moist static energy framework

In this study, we apply the moist static energy for first baroclinic mode (MSEB) model to examine the drivers of the mean tropical atmospheric circulation biases over oceanic regions. The model diagnoses the vertical motion in an air column of the tropical regions based on net energy heat flux and advection of moisture or heat into the air column in relation to the stability of the air column due to the gradients in moist static energy. Analysis of Coupled Model Intercomparison Project (CMIP) and Atmospheric Model Intercomparison Project (AMIP) simulations helped to identified errors intrinsic to the atmospheric models or errors due to atmosphere–ocean coupling process. Despite some limitations of the MSEB model, our multi-model mean analysis over the entire tropical ocean reveals that the primary drivers of the tropical circulation biases mostly result from intrinsic atmospheric model errors in top of the atmosphere longwave radiation and surface latent heats fluxes, suggesting a link to biases in the hydrological cycle. Oceanic coupling significantly enhanced some of the biases. Biases in the advection of moist static energy also play an important role, while biases in the gross moist stability profiles play only a minor role. Further, we examine the inter-model variations in four main regional large-scale biases (double-ITCZ, Pacific cold tongue, southward shift of ITCZ over the Atlantic, and dipole bias over the Indian Ocean). The analysis suggests that regional bias patterns across general circulation models are primarily driven by coupling errors, except for the bias in the Indian Ocean, which is intrinsic to the atmospheric model but amplified by coupling. Notably, longwave radiation biases at the top of the atmosphere are prevalent among the four bias patterns, as well as biases in moisture advection over the Atlantic. Our results underscore the significant role of net longwave radiation at the top of the atmosphere, an aspect not sufficiently emphasized in previous studies.

Quantifying the Role of Ocean Dynamics in SST Variability across GCMs and Observations

Abstract Midlatitude SSTs forced by mesoscale oceanic processes can affect the large-scale atmosphere, pointing to the ocean's crucial role outside the tropics. Previous studies have shown oceanic mesoscale processes' effect on global and regional climate variability. This study quantifies the local contribution of ocean dynamics to mixed-layer temperature across the globe by directly estimating the ocean heat flux divergence resolved by state-of-the-art ocean reanalysis, eddy-resolving, and eddy-parameterized versions of two US national climate models and indirectly from air-sea flux satellite-based estimates. Our results show that the eddy-resolving climate simulations resolve mixed-layer temperature variances that are larger and closer to those inferred from observations than both their eddy-parameterized counterparts and ECCO over much of the extratropics. The observations and the eddy-resolving models indicate a more significant role of ocean dynamics in the mixed layer temperature variability than the surface fluxes over most extratropics compared to their eddy-parameterized versions. A frequency domain analysis shows that the better-resolved ocean mesoscale and thermal gradients enhance the variance over timescale from two months to thirty years. Results show agreement in the ocean's contribution among satellite-based estimates, ocean reanalysis products, and ocean eddy-resolving simulations. At the same time, differences emerge for ECCO and the eddy-parameterized models, suggesting that surface fluxes account for a larger fraction of the mixed layer temperature variability in most of the extratropics.

Disruptions in thermohaline staircases caused by subsurface mesoscale eddies in the eastern Caribbean Sea

Thermohaline staircases and mesoscale eddies play crucial roles in the transport of heat, salt, and nutrients in the ocean, yet their complex relationship remains unclear due to the limitations in observational resolution and lateral range. Here, utilizing high spatial resolution seismic images with a total length of 2518 km, we show the widespread occurrence of thermohaline staircases in the eastern Caribbean Sea and demonstrate their three-dimensional distribution. These staircases occupy about 70% of the total length of the seismic lines, with the remaining 30% occupied by six subsurface mesoscale eddies. These staircases are disrupted and interrupted due to the enhanced turbulent diffusivity driven by the vertical shear of these eddies. The numerous high-resolution seismic observations presented herein enhance our comprehension of the interplay between staircases and multiscale oceanic dynamical processes.

Kilometer-scale ocean processes behind the variability of the Island Mass Effect in the Maldives

The Island Mass Effect (IME) around the Maldives is responsible for intense blooms with distinct seasonal patterns. These blooms sustain the fishing industry of the archipelagic nation, a vital source of income that occupies about 30% of the population. Through high resolution ocean simulations, we explore the physical processes responsible for the increased productivity and its observed variability, and their sensitivity to changes in land distribution. Year-round the frictional break of the currents in the presence of shallow bathymetry produces a strong vertical shear in the flow that favors vertical mixing, independently of the currents direction. The impact of this mixing is visible at the ocean surface during March and April, when the mixed-layer is shallow and the ocean currents are generally weak. A different mechanism than observed in spring modulates the IME during the monsoon seasons, in both winter and summer, when the zonal currents forced by the strong winds cross the archipelago: the flow accelerates while squeezing between the atolls, giving rise to intense wakes, and strong upwelling originates in the lees in response to the kilometer-scale flow divergence. This upwelling creates an asymmetric cooling signal in the lee of the islands and obfuscates the effects of the enhanced vertical mixing. The land reclamation efforts being planned on some of the islands may undercut the upwelling that drives the blooms during the monsoon seasons.

Two distinct early Paleozoic subduction zones in the eastern Central Asian Orogenic Belt: Evidence of subduction recycling and arc evolution

This paper presents zircon U–Pb dating and Hf isotopic, and whole-rock geochemical and Sr–Nd isotopic data for early Paleozoic igneous rocks as well as zircon U–Pb dating for the sedimentary rocks of the Songnen Massif within the eastern Central Asian Orogenic Belt (CAOB), aiming to reveal arc magma evolution process. The zircon U–Pb dating results of the siltstone of the Xiaojingou Formation yielded 450 Ma, 471 Ma, 502 Ma, 773 Ma, 818 Ma and 949 Ma age peaks, as well as multi-episodic Paleo-proterozoic ages (2486-1612 Ma), implying its deposition time is younger than ∼ 450 Ma. The zircon U–Pb dating results reveal that the Songnen Massif exists three episodes of volcanic rocks of the middle Cambrian (502 Ma), Middle Ordovician (453-449 Ma) and the Silurian (435-420 Ma). The middle Cambrian basalts are low-K tholeiitic to calc-alkaline series, weakly depleted εHf(t) and εNd(t) values. They were derived from the weakly depleted mantle metasomatized by slab fluid in an immature island arc setting. The Middle Ordovician to Silurian volcanic rocks show SiO2 of 49.4–58.1 wt%, plotting into high-K calc-alkaline series, which are also featured by higher Th/Nd and (Hf/Sm)PM ratios, together with moderate depleted εNd(t) values, suggesting their derivation of the depleted mantle wedge metasomatized by sediment-derived melt. In contrast, coeval early Paleozoic volcanic rocks in the eastern Xing’an Massif with the MORB-like εNd(t) values are characterized by an origin of the depleted mantle wedge and the distinct interaction of oceanic crust melt and mantle peridotite. Combined with the variation of calculated crust thickness, we suggested that there was a continental arc in the eastern Songnen Massif, which developed from immaturity to maturity during middle Cambrian to Silurian, while the early Paleozoic island arc in the eastern Xing’an Massif is characterized by ocean ridge subduction process.