Global Ocean General Circulation Model Research Articles

Simulated inventory and distribution of 137Cs released from multiple sources in the global ocean

Radioactive cesium (137Cs) is distributed in the world's oceans as a result of global fallout from atmospheric nuclear weapons tests, releases from fuel reprocessing plants, and inputs from nuclear power plant accident. In order to detect future radionuclide contamination, it is necessary to establish a baseline global distribution of radionuclides such as 137Cs and to understand the ocean transport processes that lead to that distribution. In order to aid in the interpretation of the observed database, we have conducted a suite of simulations of the distribution of 137Cs using a global ocean general circulation model (OGCM). Simulated 137Cs radioactivity concentrations agree well with observations, and the results were used to estimate the changes in inventories for each ocean basin. 137Cs activity concentration from atmospheric nuclear weapons tests are expected to be detectable in the world ocean until at least 2030.

Seamount-induced mixing revealed through idealized experiments and its parameterization in an Oceanic General Circulation Model

Though hidden under hundreds if not thousands of meters of seawater, seamounts and the vertical mixing they induce are crucial components of mixing in the global ocean. To gain a deeper understanding of the influence of different seamount topographies on vertical mixing, and distributional characteristics and causes of seamount-induced vertical mixing in different grid scales, this study builds a three-dimensional idealized seamount in a regional ocean model. Results illustrate that, given a constant westward inflow, seamount-induced vertical mixing is mainly distributed on the north and south sides of the summit, the foot of the north and east sides, and downstream from the seamount. Mixing on the north and south sides of the seamount summit, and at its foot, is caused by density fronts induced by friction between ambient currents and the seamount. By contrast, mixing downstream of the seamount is due to convective instability, lee waves, and seamount wake eddies. Vertical mixing enhancement at subgrid scales mainly occurs on the southern and eastern sides of the seamount summit, at its foot, and at the seabed. Here, vertical mixing is primarily caused by friction in the bottom boundary layer. Sensitivity experiments reveal that mixing intensity is affected by seamount topographic roughness. As the seamount height and slope increases, vertical mixing at the seamount summit is enhanced, and mixing on its slope gradually shifts to the seamount base. Due to limitations of global ocean model resolutions, seamount-induced vertical mixing is often ignored. Consequently, a global background mixing coefficient (BMC) relationship is built from an Argo dataset-estimated BMC combined with seafloor topography roughness, and then applied to a global Oceanic General Circulation Model (OGCM). Results show that simulations of deep-sea water temperature and salinity improved after modifying the BMC in seamount-rich areas, with a simultaneous reduction of root-mean-square error by 17% and 12%, respectively.

Local Drivers of Extreme Upper Ocean Marine Heatwaves Assessed Using a Global Ocean Circulation Model

The growing threat of Marine heatwaves (MHWs) to ecosystems demands that we better understand their physical drivers. This information can be used to improve the performance of ocean models in predicting major events so more appropriate management decisions can be made. Air-sea heat fluxes have been found to be one of the dominant drivers of MHWs but their impact are expected to decrease for MHWs extending deeper into the water column. In this study, we examine the most extreme MHWs occurring within an upper ocean layer and quantify the relative contributions of oceanic and atmospheric processes to their onset and decay phases. The base of the upper ocean layer is defined as the local winter mixed layer depth so that summer events occurring within a shallower mixed layer are also included. We perform a local upper ocean heat budget analysis at each grid point of a global ocean general circulation model. Results show that in 78% of MHWs, horizontal heat convergence is the main driver of MHW onset. In contrast, heat fluxes dominate the formation of MHWs in 11% of cases, through decreased latent heat cooling and/or increased solar radiation. These air-sea heat flux driven events occur mostly in the tropical regions where the upper ocean layer is shallow. In terms of MHW decay, heat advection is dominant in only 31% of MHWs, while heat flux dominance increases to 23%. For the majority of remaining events, advection and air-sea heat flux anomalies acted together to dissipate the excessive heat. This shift toward a comparable contribution of advection and air-sea heat flux is a common feature of extreme MHW decay globally. The anomalous air-sea heat flux cooling is mostly due to an increased latent heat loss feedback response to upper ocean temperature anomalies. Extreme upper ocean MHWs coincided with SST MHWs consistently, but with lower intensity in extra-tropical regions, where the upper ocean layer is deeper. This suggests that the upper ocean heat accumulation may pre-condition the SST MHWs in these regions. Our analysis provides valuable insights into the local physical processes controlling the onset and decay of extreme MHWs.

Incorporating tides and internal gravity waves within global ocean general circulation models: A review

Incorporating tides and internal gravity waves within global ocean general circulation models: A review

The Spectral Diagram as a new tool for model assessment in the frequency domain: Application to a global ocean general circulation model with tides

Here we introduce a new tool, the Spectral Diagram (SD), for the comparison of time series in the frequency domain. The SD provides a novel way to display the coherence function, power, amplitude, phase, and skill score of discrete frequencies of two time series. Each SD summarises these quantities in a single plot for multiple targeted frequencies. The versatility of SDs is demonstrated through a series of sea-level comparisons between observations from tide gauges and the model results from a global eddy-permitting ocean general circulation model (MOM5) with explicit tidal forcing. Phase information for the eight principal lunisolar constituents (M2, S2, N2, K2, K1, O1, P1, Q1) is added to the default configuration of MOM5. Inaccurate estimation of phase information is an important source of barotropic errors in ocean modelling, thereby compromising the skill scores in regions where amplitudes are close to the tidal gauge datasets. The greatest contribution of SD analysis is the indication that some diurnal estimates can be improved by adjusting the phase lag in the model as severe underestimation of semidiurnal amplitudes is the main reason for lower skill scores despite higher coherence. Although the SD has been designed for tidal analysis, it is a powerful tool for detecting co-oscillating patterns in multi-scale analyses, and this approach might provide guidance in devising skill scores for inter-comparing model results.

Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain

Abstract. The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.

A Barotropic Solver for High-Resolution Ocean General Circulation Models

High-resolution global ocean general circulation models (OGCMs) play a key role in accurate ocean forecasting. However, the models of the operational forecasting systems are still not in high resolution due to the subsequent high demand for large computation, as well as the low parallel efficiency barrier. Good scalability is an important index of parallel efficiency and is still a challenge for OGCMs. We found that the communication cost in a barotropic solver, namely, the preconditioned conjugate gradient (PCG) method, is the key bottleneck for scalability due to the high frequency of the global reductions. In this work, we developed a new algorithm—a communication-avoiding Krylov subspace method with a PCG (CA-PCG)—to improve scalability and then applied it to the Nucleus for European Modelling of the Ocean (NEMO) as an example. For PCG, inner product operations with global communication were needed in every iteration, while for CA-PCG, inner product operations were only needed every eight iterations. Therefore, the global communication cost decreased from more than 94.5% of the total execution time with PCG to less than 63.4% with CA-PCG. As a result, the execution time of the barotropic modes decreased from more than 17,000 s with PCG to less than 6000 s with CA-PCG, and the total execution time decreased from more than 18,000 s with PCG to less than 6200 s with CA-PCG. Besides, the ratio of the speedup can also be increased from 3.7 to 4.6. In summary, the high process count scalability when using CA-PCG was effectively improved from that using the PCG method, providing a highly effective solution for accurate ocean simulation.

Depth Structure of Ningaloo Niño/Niña Events and Associated Drivers

AbstractMarine heatwaves along the coast of Western Australia, referred to as Ningaloo Niño, have had dramatic impacts on the ecosystem in the recent decade. A number of local and remote forcing mechanisms have been put forward; however, little is known about the depth structure of such temperature extremes. Utilizing an eddy-active global ocean general circulation model, Ningaloo Niño and the corresponding cold Ningaloo Niña events are investigated between 1958 and 2016, with a focus on their depth structure. The relative roles of buoyancy and wind forcing are inferred from sensitivity experiments. Composites reveal a strong symmetry between cold and warm events in their vertical structure and associated large-scale spatial patterns. Temperature anomalies are largest at the surface, where buoyancy forcing is dominant, and extend down to 300-m depth (or deeper), with wind forcing being the main driver. Large-scale subsurface anomalies arise from a vertical modulation of the thermocline, extending from the western Pacific into the tropical eastern Indian Ocean. The strongest Ningaloo Niños in 2000 and 2011 are unprecedented compound events, where long-lasting high temperatures are accompanied by extreme freshening, which emerges in association with La Niñas, that is more common and persistent during the negative phase of the interdecadal Pacific oscillation. It is shown that Ningaloo Niños during La Niña phases have a distinctively deeper reach and are associated with a strengthening of the Leeuwin Current, while events during El Niño are limited to the surface layer temperatures, likely driven by local atmosphere–ocean feedbacks, without a clear imprint on salinity and velocity.

Global tide simulations with ICON-O: testing the model performance on highly irregular meshes

The global tide is simulated with the global ocean general circulation model ICON-O using a newly developed tidal module, which computes the full tidal potential. The simulated coastal M2 amplitudes, derived by a discrete Fourier transformation of the output sea level time series, are compared with the according values derived from satellite altimetry (TPXO-8 atlas). The experiments are repeated with four uniform and sixteen irregular triangular grids. The results show that the quality of the coastal tide simulation depends primarily on the coastal resolution and that the ocean interior can be resolved up to twenty times lower without causing considerable reductions in quality. The mesh transition zones between areas of different resolutions are formed by cell bisection and subsequent local spring optimisation tolerating a triangular cell’s maximum angle up to 84°. Numerical problems with these high-grade non-equiangular cells were not encountered. The results emphasise the numerical feasibility and potential efficiency of highly irregular computational meshes used by ICON-O.

The IOD Impacts on the Indian Ocean Carbon Cycle

AbstractThis paper examines the impacts of Indian Ocean Dipole (IOD) Mode on the upper ocean carbon cycle and its variability in the Indian Ocean with available biogeochemical observations and 60 years (1960–2019) of model outputs from a global ocean biogeochemical general circulation model. The upper ocean carbon cycle variability of the Indian Ocean is faithfully reproduced by the model when compared with available observations. The IOD leads to a substantial sea‐to‐air CO2 flux variability in the southeastern tropical Indian Ocean over a broad region (70–105°E, 0–20°S), with more focus near the coast of Java‐Sumatra due to the prevailing upwelling dynamics and associated westward propagating anomalies. The sea‐to‐air CO2 fluxes, surface ocean partial pressure of CO2 (pCO2), the concentration of dissolved inorganic carbon (DIC), and ocean alkalinity (ALK) range as much as ±1.0 mole m−2 year−1, ±20 μatm, ±35 μmole kg−1, and ±22 μmole kg−1 within 80–105°E, 0–10°S due to IOD. The DIC and ALK are significant drivers of pCO2 variability associated with IOD. The roles of temperature (T) and biology are found negligible. A relatively warm T and extremely high freshwater forcing make the southeastern tropical Indian Ocean carbon cycle variability submissive to DIC and ALK evolutions in contrast to the tropical eastern Pacific where changes in DIC and T dominate the pCO2 interannual variability. For the first time, this study provides a most comprehensive and extended analysis for the region while highlighting significant differences in carbon cycle variability of the eastern tropical Indian Ocean compared to that of the other parts of the global oceans.

Is interactive air sea coupling relevant for simulating the future climate of Europe?

The majority of regional climate change assessments for the Euro-CORDEX region is based on high resolution atmosphere models. These models use prescribed lower boundary conditions, such as sea surface temperatures (SST) from global ocean General Circulation Models (GCMs), that do not respond to changes simulated by the regional atmosphere model, thus lacking an important feedback to the atmosphere. However, research during the past decade indicated that the use of coupled atmosphere–ocean models can lead to significantly altered model solutions compared to standalone atmosphere models for the present day climate imposing some uncertainty on the widely used uncoupled future scenarios. We here present the first multi-model and multi scenario (RCP2.6, RCP4.5, RCP8.5) ensemble of future climate change scenarios downscaled with a coupled atmosphere—ocean model in which sea surface temperature and sea ice fields are explicitly simulated by a coupled state-of-the-art high resolution ocean model and communicated to the atmosphere at 3-hourly time steps. Our ensemble generally confirms results of previous uncoupled ensembles over land areas implying that the coupling effect is restricted mainly to the coupled area and the adjacent coastal zone. By contrast, over the North Sea and Baltic Sea small scale processes point to important coupling effects that mediate the response to climate change and that can not be simulated by uncoupled models. Our results therefore impose general uncertainty on the usage of regional climate change data from uncoupled ensembles over marine areas such as for purposes of offshore wind or mussel farming, the planing of marine protected areas, and marine recreation along the coastal zone. It further sets in question the usage of uncoupled scenario data (such as Euro-CORDEX) to force high resolution ocean models. Comparing coupled and uncoupled hindcast simulations reveals that the coupling effect over land is most pronounced during the warm season when prescribed and modelled sea surface temperatures (SST) differ strongest. In addition, a generally weaker wind regime in summer damps the heat dispersion in the atmosphere so that air temperature anomalies can extent further over land compared to winter. Future projections are discussed under consideration of land-sea warming characteristics for selected climate indices as well as mean seasonal climate change. At the end of the century a clear land-sea pattern is seen in all scenarios with stronger warming over land than over open sea areas. On average land areas warm at a rate 1.5 times faster than areas over the open ocean. Over the coupled area, i.e. the North Sea and Baltic Sea tropical nights are impacted strongest and the Baltic Sea turns out to be a hot spot in future climate. This has been unrecognized in previous studies using high resolution atmosphere models with prescribed SSTs from global models which do not represent small scale ocean processes in the Baltic Sea adequately.

Resolution dependency of sinking Lagrangian particles in ocean general circulation models.

Any type of non-buoyant material in the ocean is transported horizontally by currents during its sinking journey. This lateral transport can be far from negligible for small sinking velocities. To estimate its magnitude and direction, the material is often modelled as a set of Lagrangian particles advected by current velocities that are obtained from Ocean General Circulation Models (OGCMs). State-of-the-art OGCMs are strongly eddying, similar to the real ocean, providing results with a spatial resolution on the order of 10 km on a daily frequency. While the importance of eddies in OGCMs is well-appreciated in the physical oceanographic community, other marine research communities may not. Further, many long term climate modelling simulations (e.g. in paleoclimate) rely on lower spatial resolution models that do not capture mesoscale features. To demonstrate how much the absence of mesoscale features in low-resolution models influences the Lagrangian particle transport, we simulate the transport of sinking Lagrangian particles using low- and high-resolution global OGCMs, and assess the lateral transport differences resulting from the difference in spatial and temporal model resolution. We find major differences between the transport in the non-eddying OGCM and in the eddying OGCM. Addition of stochastic noise to the particle trajectories in the non-eddying OGCM parameterises the effect of eddies well in some cases (e.g. in the North Pacific gyre). The effect of a coarser temporal resolution (once every 5 days versus monthly) is smaller compared to a coarser spatial resolution (0.1° versus 1° horizontally). We recommend to use sinking Lagrangian particles, representing e.g. marine snow, microplankton or sinking plastic, only with velocity fields from eddying Eulerian OGCMs, requiring high-resolution models in e.g. paleoceanographic studies. To increase the accessibility of our particle trace simulations, we launch planktondrift.science.uu.nl, an online tool to reconstruct the surface origin of sedimentary particles in a specific location.

Simulation of Deep Cycle Turbulence by a Global Ocean General Circulation Model

AbstractDeep cycle turbulence (DCT) is a diurnally oscillating turbulence that penetrates into a stratified shear layer below the surface mixed layer, which is often observed in the eastern Pacific and Atlantic above the Equatorial Undercurrent (EUC). Here we present the simulation of DCT by a global ocean general circulation model (OGCM) for the first time. As the k‐ε vertical mixing scheme is used in the OGCM, the simulation of observed DCT structure based on in situ microstructure measurements can be explicitly demonstrated. The simulated DCT is found in all equatorial ocean basins, and its characteristics agree very well with observations. Zonal and meridional variations of DCT in the entire equatorial Pacific and Atlantic are described through constructing the composite diurnal cycle. In the central Pacific where the maximum shear associated with EUC is deep, the separation of DCT from the surface mixed layer is much more prominent than other areas.

A non-breaking-wave-generated turbulence mixing scheme for a global ocean general circulation model

A novel vertical mixing scheme to describe the influence of the non-breaking surface waves in ocean general circulation models is proposed based on the second-order turbulence closure model and a mathematical relationship fitted between the in situ observations of the turbulence kinetic energy dissipation rate and the calculated velocity shear module of non-breaking surface waves. The non-breaking-wave-generated turbulence mixing coefficients can be calculated empirically in terms of the wave spectral parameters including ΦSW, K, and ω, where ΦSW is the wave number spectrum of the significant wave height, K is the wave number, ω is the surface wave frequency. The effect of the new mixing scheme on global ocean circulation simulation was tested using the MArine Science and NUmerical Modeling ocean model. The results indicate significant improvement in the upper-ocean temperature structure and the mixed layer depth, i.e., the non-breaking-wave-generated turbulence mixing plays an important role in the upper ocean at high latitudes in summer for both hemispheres compared with traditional Mellor-Yamada 2.5 scheme.

Energy Flux Observations in an Internal Tide Beam in the Eastern North Atlantic

AbstractLow‐mode internal waves propagate over large distances and provide energy for turbulent mixing when they break far from their generation sites. A realistic representation of the oceanic energy cycle in ocean and climate models requires a consistent implementation of their generation, propagation, and dissipation. Here we combine the long‐term mean energy flux from satellite altimetry with results from a 1/10° global ocean general circulation model that resolves the low modes of internal waves and in situ observations of stratification and horizontal currents to study energy flux and dissipation along a 1000 km internal tide beam in the eastern North Atlantic. Internal wave fluxes were estimated from twelve 36‐ to 48‐hr stations in along‐ and across‐beam direction to resolve both the inertial period and tidal cycle. The observed internal tide energy fluxes range from 5.9 kW m−1 near the generation sites to 0.5 kW m−1 at distant stations. Estimates of energy dissipation come from both finestructure and upper ocean microstructure profiles and range, vertically integrated, from 0.5 to 3.3 mW m−2 along the beam. Overall, the in situ observations confirm the internal tide pattern derived from satellite altimetry, but the in situ energy fluxes are more variable and decrease less monotonically along the beam. Internal tides in the model propagate over shorter distances compared to results from altimetry and in situ measurements, but more spatial details close the main generation sites are resolved.

Evaluation of an Online Grid‐Coarsening Algorithm in a Global Eddy‐Admitting Ocean Biogeochemical Model

AbstractIn order to explore the effects of mesoscale eddies on marine biogeochemistry over climate timescales, global ocean biogeochemical general circulation models (OBGCMs) need at least to be run at a horizontal resolution of a 0.25°, the minimal resolution admitting eddies. However, their use is currently limited because of a prohibitive computational cost and storage requirements. To overcome this problem, an online coarsening algorithm is evaluated in the oceanic component (NEMO‐GELATO‐PISCES) of CNRM‐ESM2‐1. This algorithm allows to compute biogeochemical processes at a coarse resolution (0.75°) while inheriting most of the dynamical characteristics of the eddy‐admitting OBGCM (0.25°). Through the coarse‐graining process, the effective resolution of the ocean dynamics seen by the biogeochemical model is higher than that which would be obtained from an OBGCM run at 0.75°. In this context, we assess how much the increase from low (1°) to coarse‐grained horizontal resolution impacts the ocean dynamics and the marine biogeochemistry over long‐term climate simulations. The online coarsening reduces the computational cost by 60% with respect to that of the eddy‐admitting OBGCM. In addition, it improves the representation of chlorophyll, nutrients, oxygen, and sea‐air carbon fluxes over more than half of the open ocean area compared to the 1° OBGCM. Most importantly, the coarse‐grained OBGCM captures the physical‐biogeochemical coupling between sea‐air carbon fluxes and sea surface height and between oxygen minimum zone boundaries and eddies, as produced by the eddy‐admitting OBGCM. Such a cost‐efficient coarsening algorithm offers a good trade‐off to conduct process‐based studies over centennial timescales at higher resolution.

Improving the Upper‐Ocean Temperature in an Ocean Climate Model (FESOM 1.4): Shortwave Penetration Versus Mixing Induced by Nonbreaking Surface Waves

AbstractAs the first mature global ocean general circulation model based on unstructured‐mesh methods, the multiresolution Finite Element Sea ice‐Ocean Model (FESOM) has shown great capability in reconstructing the ocean and sea ice in both standalone and coupled simulations at a relatively low computational cost. Parameterizations of some important processes, including the vertical mixing induced by surface waves, however, are still missing, contributing to temperature biases in the upper ocean. In this work we incorporate the vertical mixing induced by nonbreaking surface waves derived from a wave model into FESOM and compare its effect with that of shortwave penetration, another key process to vertically redistribute the heat in the upper ocean. Numerical experiments with and without the shortwave penetration scheme and the nonbreaking surface wave mixing reveal that both processes ameliorate the simulation of upper‐ocean temperature in middle and low latitudes mainly on the summer hemisphere. The role of nonbreaking surface waves is more pronounced in decreasing the mean cold biases at 50 m (by 1.0 °C, in comparison to 0.5 °C achieved by applying shortwave penetration). We conclude that the incorporation of mixing induced by nonbreaking surface waves into FESOM is practically very helpful and suggest that it needs to be considered in other ocean climate models as well.

Recent Sea Ice Decline Did Not Significantly Increase the Total Liquid Freshwater Content of the Arctic Ocean

Abstract The freshwater stored in the Arctic Ocean is an important component of the global climate system. Currently the Arctic liquid freshwater content (FWC) has reached a record high since the beginning of the last century. In this study we use numerical simulations to investigate the impact of sea ice decline on the Arctic liquid FWC and its spatial distribution. The global unstructured-mesh ocean general circulation model Finite Element Sea Ice–Ocean Model (FESOM) with 4.5-km horizontal resolution in the Arctic region is applied. The simulations show that sea ice decline increases the FWC by freshening the ocean through sea ice meltwater and modifies upper ocean circulation at the same time. The two effects together significantly increase the freshwater stored in the Amerasian basin and reduce its amount in the Eurasian basin. The salinification of the upper Eurasian basin is mainly caused by the reduction in the proportion of Pacific Water and the increase in that of Atlantic Water (AW). Consequently, the sea ice decline did not significantly contribute to the observed rapid increase in the Arctic total liquid FWC. However, the changes in the Arctic freshwater spatial distribution indicate that the influence of sea ice decline on the ocean environment is remarkable. Sea ice decline increases the amount of Barents Sea branch AW in the upper Arctic Ocean, thus reducing its supply to the deeper Arctic layers. This study suggests that all the dynamical processes sensitive to sea ice decline should be taken into account when understanding and predicting Arctic changes.

Seasonal Cycle and Annual Reversal of the Somali Current in an Eddy‐Resolving Global Ocean Model

AbstractThe annual cycle of the Somali Current is studied using a forced strongly eddy‐active global ocean general circulation model. The annual cycle and timing of the reversal of the Somali Current differ depending on the location along the coast of northeastern Africa. A momentum budget analysis calculated from the model output shows that both annual Rossby waves and the life cycle of the Great Whirl impact the northern part of the Somali Current (roughly 5°N to 10°N). In particular, the nonlinear term is of leading order importance when the Great Whirl is present. This is in contrast to the leading‐order geostrophic balance between 2°N and 5°N. South of 2°N, equatorial dynamics prevail, and the nonlinear term is again important, and its reversal is strongly affected by the northward flowing East African Coastal Current. Alongshore wind forcing is influential at all latitudes, but only the transition in flow direction during boreal spring between 2°N and 5°N directly follows that of the wind.

Transport and storage of anthropogenic C in the North Atlantic Subpolar Ocean

Abstract. The North Atlantic Ocean is a major sink region for atmospheric CO2 and contributes to the storage of anthropogenic carbon (Cant). While there is general agreement that the intensity of the meridional overturning circulation (MOC) modulates uptake, transport and storage of Cant in the North Atlantic Subpolar Ocean, processes controlling their recent variability and evolution over the 21st century remain uncertain. This study investigates the relationship between transport, air–sea flux and storage rate of Cant in the North Atlantic Subpolar Ocean over the past 53 years. Its relies on the combined analysis of a multiannual in situ data set and outputs from a global biogeochemical ocean general circulation model (NEMO–PISCES) at 1∕2∘ spatial resolution forced by an atmospheric reanalysis. Despite an underestimation of Cant transport and an overestimation of anthropogenic air–sea CO2 flux in the model, the interannual variability of the regional Cant storage rate and its driving processes were well simulated by the model. Analysis of the multi-decadal simulation revealed that the MOC intensity variability was the major driver of the Cant transport variability at 25 and 36∘ N, but not at OVIDE. At the subpolar OVIDE section, the interannual variability of Cant transport was controlled by the accumulation of Cant in the MOC upper limb. At multi-decadal timescales, long-term changes in the North Atlantic storage rate of Cant were driven by the increase in air–sea fluxes of anthropogenic CO2. North Atlantic Central Water played a key role for storing Cant in the upper layer of the subtropical region and for supplying Cant to Intermediate Water and North Atlantic Deep Water. The transfer of Cant from surface to deep waters occurred mainly north of the OVIDE section. Most of the Cant transferred to the deep ocean was stored in the subpolar region, while the remainder was exported to the subtropical gyre within the lower MOC.