Ocean General Circulation Model Research Articles

Vertical Velocity 3D Estimates from the Linear Vorticity Balance in the North Atlantic Ocean

Abstract Large-scale three-dimensional vertical movements of the world’s oceans play a fundamental role in their functioning, yet they are still overlooked. This study aims to enhance our understanding of the relationship between horizontal and vertical flows as a function of depth to derive some major characteristics of the w field from a relatively well-constrained variable. For this purpose, we analyze an eddy-permitting ocean general circulation model simulation, focusing on the linear vorticity balance (LVB) over the North Atlantic basin. For spatial scales greater than 5°, over more than 50-yr average, the LVB holds with less than 10% error over the whole depth range of most of the subtropical upper tachocline–upper ocean and part of the tropical one, away from western boundary currents and zonal tropical flows. The spatial extent of the balance diminishes in the intermediate ocean, particularly in the western half of the basin. Within the deep tachocline, the balance holds over large portions of the abyssal plains. Based on these results, the β-plane geostrophic vertical velocity is derived from the vertical indefinite integral of the geostrophic LVB from the surface. While its time mean aligns only qualitatively with the model reference w, it accurately captures the interannual variability over most of the basin, even close to eastern boundaries. However, it cannot reproduce the variability in the lower tropics, the western boundary, and the North Atlantic Current system. Our results emphasize the dominance of geostrophy in establishing most of the ocean model vertical flow variability.

GRACE Satellite Observations of Antarctic Bottom Water Transport Variability

AbstractAntarctic Bottom Water (AABW) formation and transport constitute a key component of the global ocean circulation. Direct observations suggest that AABW volumes and transport rates may be decreasing, but these observations are too temporally or spatially sparse to determine the cause. To address this problem, we develop a new method to reconstruct AABW transport variability using data from the GRACE (Gravity Recovery and Climate Experiment) satellite mission. We use an ocean general circulation model to investigate the relationship between ocean bottom pressure and AABW: we calculate both of these quantities in the model, and link them using a regularized linear regression. Our reconstruction from modeled ocean bottom pressure can capture 65%–90% of modeled AABW transport variability, depending on the ocean basin. When realistic observational uncertainty values are added to the modeled ocean bottom pressure, the reconstruction can still capture 30%–80% of AABW transport variability. Using the same regression values, the reconstruction skill is within the same range in a second, independent, general circulation model. We conclude that our reconstruction method is not unique to the model in which it was developed and can be applied to GRACE satellite observations of ocean bottom pressure. These advances allow us to create the first global reconstruction of AABW transport variability over the satellite era. Our reconstruction provides information on the interannual variability of AABW transport, but more accurate observations are needed to discern AABW transport trends.

Effects of Horizontal Resolution on Long‐Range Equatorward Radiation of Near‐Inertial Internal Waves in Ocean General Circulation Models

AbstractWind‐generated near‐inertial internal waves (NIWs) are characterized by dominant long‐range equatorward radiation due to the gradient of the planetary vorticity, known as the β‐refraction effect. In this study, we analyze the effects of horizontal model resolution on the long‐range equatorward radiation of NIWs. In a high‐resolution Community Earth System Model (CESM‐HR) with a 0.1° oceanic resolution, about 25% (15%) of NIW energy flux injected downward the surface boundary layer base poleward of 30°N (30°S) radiates into the lower‐latitude region. This ratio decreases to about 15% (8%) in a low‐resolution CESM (CESM‐LR) with a 1° oceanic resolution. The higher long‐range equatorward radiation efficiency in the CESM‐HR than the CESM‐LR is directly attributed to the faster equatorward group velocity of the NIWs of the first three vertical modes, which reflects the better representation of equatorward propagation and beta‐refraction of smaller scale NIWs in the CESM‐HR. The enhancement of equatorward wavenumber induced by the β‐refraction is inhibited in the CESM‐LR, which underrepresent the long‐range equatorward radiation of NIWs. These results underscore the necessity of high‐resolution ocean models in accurately simulating the spatial variabilities of NIWs and their induced turbulent diapycnal mixing in the global ocean.

Interacting internal waves explain global patterns of interior ocean mixing

Across the stable density stratification of the abyssal ocean, deep dense water is slowly propelled upward by sustained, though irregular, turbulent mixing. The resulting mean upwelling determines large-scale oceanic circulation properties like heat and carbon transport. In the ocean interior, this turbulent mixing is caused mainly by breaking internal waves: generated predominantly by winds and tides, these waves interact nonlinearly, transferring energy downscale, and finally become unstable, break and mix the water column. This paradigm, long parameterized heuristically, still lacks full theoretical explanation. Here, we close this gap using wave-wave interaction theory with input from both localized and global observations. We find near-ubiquitous agreement between first-principle predictions and observed mixing patterns in the global ocean interior. Our findings lay the foundations for a wave-driven mixing parameterization for ocean general circulation models that is entirely physics-based, which is key to reliably represent future climate states that could differ substantially from today’s.

Spatially Heterogeneous Responses of Planktonic Foraminiferal Assemblages Over 700,000 Years of Climate Change

ABSTRACTAimTo determine the degree to which assemblages of planktonic foraminifera track thermal conditions.LocationThe world's oceans.Time PeriodThe last 700,000 years of glacial–interglacial cycles.Major Taxa StudiedPlanktonic foraminifera.MethodsWe investigate assemblage dynamics in planktonic foraminifera in response to temperature changes using a global dataset of Quaternary planktonic foraminifera, together with a coupled Atmosphere–Ocean General Circulation Model (AOGCM) at 8000‐year resolution. We use ‘thermal deviance’ to assess assemblage responses to climate change, defined as the difference between the temperature at a given location and the bio‐indicated temperature (i.e., the abundance‐weighted average of estimated temperature optima for the species present).ResultsAssemblages generally tracked annual mean temperature changes through compositional turnover, but thermal deviances are evident under certain conditions. The coldest‐adapted species persisted in polar regions during warming but were not joined by additional immigrants, resulting in minimal assemblage turnover. The warmest‐adapted species persisted in equatorial regions during cooling with similarly minimal assemblage change. Assemblages at mid‐latitudes mostly tracked temperature cooling and warming.Main ConclusionsPlanktonic foraminiferal assemblages were generally able to track or endure temperature changes: as climate warmed or cooled, bio‐indicated temperature also became warmer or cooler, although to a variable degree. At polar sites under warming and at equatorial sites under cooling, the change in bio‐indicated temperature was less than, or even opposite to, what would be expected from estimated environmental change. Nevertheless, all studied species persisted across the study interval, regardless of thermal deviance—a result that highlights the resilience and inertia of planktonic foraminifera on an assemblage level to the last 700,000 years of climate change.

Mixed Rossby Gravity Waves at Middepths of the Equatorial Indian Ocean

Abstract This study examines mixed Rossby gravity (MRG) waves at middepths (500–2000 m) of the Indian Ocean, using multiyear velocity time series obtained from current-meter moorings at 77°, 83°, and 93°E along the equator over the period 2000–19. These data are analyzed in combination with a high-resolution wind-forced ocean general circulation model. The spectrum of observed meridional velocity showed elevated energy over a wide range of periods from about 10 to 100 days with the spectral peak at a period of about 30 days. The model was able to simulate the characteristics of the observed spectrum. Further diagnostics determined that the detected variability is generally consistent with theoretical MRG waves in a resting ocean. Statistical analysis and a model sensitivity experiment identified distinct variations at three periods, where meridional velocity has sizable energy. The 14-day variability is wind driven and has a long zonal (∼3300 km) and vertical (∼4200 m) wavelength. The 28-day variability is excited by the dynamical instability of the background flow in the equatorial western Indian Ocean near the surface and propagates to the study area. It is characterized by a shorter zonal (∼1100 km) and vertical (∼2800 m) wavelength compared to 14-day variability. The 43-day variability has a zonal wavelength (∼820 km) comparable to the 28-day variability but does not show the tendency of propagation and is likely generated in situ through nonlinear interactions. These results show that various processes contribute to the excitation of MRG waves at middepths of the Indian Ocean. Significance Statement Mixed Rossby gravity (MRG) waves are a prominent mode of intraseasonal variability in tropical oceans. They are of climatic importance owing to their effects on turbulent mixing in the deep sea and how they affect the large-scale ocean circulation. This study examines MRG waves at the middepths of the equatorial Indian Ocean using a unique set of in situ current-meter measurements and a high-resolution ocean general circulation model. We describe characteristics such as horizontal and vertical wavelengths, frequency spectra, dispersion properties, and energy sources. Results show that wind forcing at the sea surface, dynamical instability of the background flow near the surface in the western Indian Ocean, and nonlinear interactions generate MRG waves at distinct periods and wavelengths.

Combined physical and biogeochemical assessment of mesoscale eddy parameterisations in ocean models: Eddy-induced advection at eddy-permitting resolutions

Combined physical and biogeochemical assessment of mesoscale eddy parameterisations in ocean models: Eddy-induced advection at eddy-permitting resolutions

Parameterized Internal Wave Mixing in Three Ocean General Circulation Models

AbstractThe non‐local model of mixing based on internal wave breaking, IDEMIX, is implemented as an enhancement of a turbulent kinetic energy closure model in three non‐eddy resolving general circulation ocean models that differ in the discretization and choice of computational grids. In IDEMIX internal wave energy is generated by an energy flux resulting from near‐inertial waves induced by wind forcing at the surface, and at the bottom, by an energy flux that parameterizes the transfer of energy between baroclinic and barotropic tides. In all model simulations with IDEMIX, the mixing work is increased compared to the reference solutions without IDEMIX, reaching values in better agreement with finestructure observations. Furthermore, the horizontal structure of the mixing work is more realistic as a consequence of the heterogeneous forcing functions. All models with IDEMIX simulate deeper thermocline depths related to stronger shallow overturning cells in the Indo‐Pacific. In the North Atlantic, deeper mixed layers in simulations with IDEMIX are associated with an increased Atlantic overturning circulation and an increase of northward heat transports toward more realistic values. The response of the deep Indo‐Pacific overturning circulation and the weak bottom cell of the Atlantic to the inclusion of IDEMIX is incoherent between the models, suggesting that additional unidentified processes and numerical mixing may confound the analysis. Applying different tidal forcing functions leads to simulation differences that are small compared to differences between the different models or between simulations with IDEMIX and without IDEMIX.

Breaking Internal Waves and Ocean Diapycnal Diffusivity in a High‐Resolution Regional Ocean Model: Evidence of a Wave‐Turbulence Cascade

AbstractIt is generally understood that the origin of ocean diapycnal diffusivity is primarily associated with the stratified turbulence produced by breaking internal (gravity) waves (IW). However, it requires significant effort to verify diffusivity values in ocean general circulation models in any particular geographical region of the ocean due to the scarcity of microstructure measurements. Recent analyses of downscaled IW fields from an internal‐wave‐admitting global ocean simulation into higher‐resolution regional configurations northwest of Hawaii have demonstrated a much‐improved fit of the simulated IW spectra to the in‐situ profiler measurements such as the Garrett‐Munk (GM) spectrum. Here, we employ this dynamically downscaled ocean simulation to directly analyze the nature of the IW‐breaking and the wave‐turbulence cascade in this region. We employ a modified version of the Kappa Profile Parameterization (KPP) to infer what the horizontally averaged vertical profile of diapycnal diffusivity should be, and compare this to the background profile that would be employed in the ocean component of a low‐resolution coupled climate model such as the Community Earth System Model (CESM) of the US National Center for Atmospheric Research (NCAR). In pursuing this goal, we also demonstrate that the wavefield in the high‐resolution regional domain is dominated by a well‐resolved spectrum of low‐mode IWs that are predictable by solving an appropriate eigenvalue problem for stratified flow. We finally suggest a new tentative approach to improve the KPP parameterization.

Why Do Oceanic Nonlinearities Contribute Only Weakly to Extreme El Niño Events?

AbstractExtreme El Niño events have outsized global impacts and control the El Niño Southern Oscillation (ENSO) warm/cold phases asymmetries. Yet, a consensus regarding the relative contributions of atmospheric and oceanic nonlinearities to their genesis remains elusive. Here, we isolate the contribution of oceanic nonlinearities by conducting paired experiments forced with opposite wind stress anomalies in an oceanic general circulation model, which realistically simulates extreme El Niño events and oceanic nonlinearities thought to contribute to ENSO skewness (Tropical Instability Waves (TIWs), Nonlinear Dynamical Heating (NDH)). Our findings indicate a weak contribution of oceanic nonlinearities to extreme El Niño events in the eastern Pacific, owing to compensatory effects between lateral (NDH and TIWs) and vertical processes. These results hold across different vertical mixing schemes and modifications of the upper‐ocean heat budget mixed layer criterion. Our study reinforces previous research underscoring the pivotal role of atmospheric nonlinearities in shaping extreme El Niño events.

An ensemble-based data assimilation system for forecasting variability of the Northwestern Pacific ocean

An adjoint-free four-dimensional variational (a4dVar) data assimilation (DA) is implemented in an operational ocean forecast system based on an eddy-resolving ocean general circulation model for the Northwestern Pacific. Validation of the system against independent observations demonstrates that fitting the model to time-dependent satellite altimetry during a 10-day DA window leads to substantial skill improvements in the succeeding 60-day hindcast. The a4dVar corrects representation of the Kuroshio path variation south of Japan by adjusting the dynamical balance between amplitude/wavelength of the meander and flow advection. A larger ensemble spread tends to reduce the skill in representing the observed sea surface height anomaly, suggesting that it is possible to use the ensemble information for quantifying the forecast error. The ensemble information is also utilized for modification of the background error covariance (BEC), which improves the accuracy of temperature and salinity distributions. The modified BEC yields the skill decline of the Kuroshio path variation during the 60-day hindcast period, and the ensemble sensitivity analysis shows that changes in the dynamical balance caused by the ensemble BEC result in such skill deterioration.

Evaluating westward and eastward propagating mesoscale eddies using a 1/10° global ocean simulation of CAS-LICOM3

Evaluating westward and eastward propagating mesoscale eddies using a 1/10° global ocean simulation of CAS-LICOM3

An overview of the NEMO modelling for the BaySys project

This article is intended as an introduction to discuss the development of a modelling framework to examine simulated climate change and river discharge regulation and their combined impact on marine conditions in the Hudson Bay Complex as a contribution to BaySys, a collaborative project between Manitoba Hydro, Hydro-Quebec, the University of Manitoba, the University of Alberta, Université Laval and Ouranos. In support of this work, a sea ice and oceanographic model was improved and then used to further study the effects of freshwater loading and ice cover on the circulation of Hudson Bay. This modelling perspective is based on the Nucleus for European Modelling of the Ocean (NEMO) ocean general circulation model coupled to version 2 of the Louvain-la-Neuve sea ice model (LIM2). The goal of the modelling was to provide a framework and tool for simulating projected changes in marine state and dynamic variables, while also enabling an integration of observations and numerical analyses. A key aspect of this work was the climate-hydrologic-ocean model integration aspect. The inclusion of a biogeochemical model and explicit tidal forcing to examine the evolution of a Canadian marginal sea with century-long integrations was also a novel aspect of the work. Overall, this work examines the NEMO modelling configuration used in BaySys, how it is set up and the experiments carried out. A broader picture evaluation of the model output is made including the BaySys mooring observations, showing that the modelling framework is suitable to examine the posed questions on the role of climate change and river regulation.

A Dynamics-Weighted Principal Components Analysis of Dominant Atmospheric Drivers of Ocean Variability with an Application to the North Atlantic Subpolar Gyre

Abstract This paper describes a framework for identifying dominant atmospheric drivers of ocean variability. The method combines statistics of atmosphere–ocean fluxes with physics from an ocean general circulation model to derive atmospheric patterns optimized to excite variability in a specified ocean quantity of interest. We first derive the method as a weighted principal components analysis and illustrate its capabilities in a toy problem. Next, we apply our analysis to the problem of interannual upper ocean heat content (HC) variability in the North Atlantic Subpolar Gyre (SPG) using the adjoint of the MITgcm and atmosphere–ocean fluxes from the ECCOv4-r4 state estimate. An unweighted principal components analysis reveals that North Atlantic heat and momentum fluxes in ECCOv4-r4 have a range of spatiotemporal patterns. By contrast, dynamics-weighted principal components analysis collapses the space of these patterns onto a small subset—principally associated with the North Atlantic Oscillation—that dominates interannual SPG HC variance. By perturbing the ECCOv4-r4 state estimate, we illustrate the pathways along which variability propagates from the atmosphere to the ocean in a nonlinear ocean model. This technique is applicable across a range of problems across Earth system components, including in the absence of a model adjoint. Significance Statement While the oceans have absorbed 90% of the excess heat associated with human-forced climate change, the change in the ocean’s heat content is not steady, with peaks and troughs superimposed upon a general increase. These fluctuations come from chaotic changes in the atmosphere and ocean and can be hard to disentangle. We use this case of ocean heat content variability to introduce a new method for determining the patterns of weather and climate in the atmosphere that are most effective at generating fluctuations in the ocean. To do this, we combine the statistics of recent atmospheric activity with output from a state-of-the-art numerical ocean model that reveals physical processes driving changes in ocean quantities including ocean heat content. This approach suggests that the atmospheric patterns that stimulate the most energetic changes in ocean heat content in the northern North Atlantic are not necessarily the most energetic patterns present in the atmosphere. We test our findings by preventing these patterns from affecting the ocean in our numerical model and measure a strong reduction in ocean heat content fluctuations.

Energetics of eddy-mean flow interactions in the western tropical Pacific Ocean

Energetics of eddy-mean flow interactions in the western tropical Pacific Ocean

Dominant contribution of atmospheric nonlinearities to ENSO asymmetry and extreme El Niño events

Extreme El Niño events have outsized impacts and strongly contribute to the El Niño Southern Oscillation (ENSO) warm/cold phase asymmetries. There is currently no consensus on the respective importance of oceanic and atmospheric nonlinearities for those asymmetries. Here, we use atmospheric and oceanic general circulation models that reproduce ENSO asymmetries well to quantify the atmospheric nonlinearities contribution. The linear and nonlinear components of the wind stress response to Sea Surface Temperature (SST) anomalies are isolated using ensemble atmospheric experiments, and used to force oceanic experiments. The wind stress-SST nonlinearity is dominated by the deep atmospheric convective response to SST. This wind-stress nonlinearity contributes to ~ 40% of the peak amplitude of extreme El Niño events and ~ 55% of the prolonged eastern Pacific warming they generate until the following summer. This large contribution arises because nonlinearities consistently drive equatorial westerly anomalies, while the larger linear component is made less efficient by easterly anomalies in the western Pacific during fall and winter. Overall, wind-stress nonlinearities fully account for the eastern Pacific positive ENSO skewness. Our findings underscore the pivotal role of atmospheric nonlinearities in shaping extreme El Niño events and, more generally, ENSO asymmetry.

Improving Global Barotropic Tides With Sub‐Grid Scale Topography

AbstractIn recent years, efforts have been made to include tides in both operational ocean models as well as climate and earth system models. The accuracy of the barotropic tides is often limited by the model topography, which is in turn limited by model horizontal resolution. In this work, we explore the reduction of barotropic tidal errors in an ocean general circulation model (Modular Ocean Model version 6; MOM6) using sub‐grid scale topography representation. We follow the methodology from Adcroft (2013, https://doi.org/10.1016/j.ocemod.2013.03.002), which utilizes statistics from finer resolution topographic data sets to represent sub‐grid scale features with a light computational cost in a structured finite volume formulation. The geometric effect from sub‐grid scale topography can be introduced to the model with only a few parameters at each grid cell. The porous barriers, which are implemented at the walls of the grid cells, are used to modify transport between grid cells. Our results show that the globally averaged tidal error in lower‐resolution simulations is significantly reduced with the use of porous barriers. We argue this method is a potentially useful tool to improve simulations of tides (and other flows) in low‐resolution simulations.

Accounting for residual errors in atmosphere–ocean background models applied in satellite gravimetry

The Atmosphere and Ocean non-tidal De-aliasing Level-1B (AOD1B) product is widely used in precise orbit determination and satellite gravimetry to correct for transient effects of atmosphere–ocean mass variability that would otherwise alias into monthly mean global gravity fields. The most recent release is based on the global ERA5 reanalysis and ECMWF operational data together with simulations from the general ocean circulation model MPIOM consistently forced with fields from the corresponding atmospheric dataset. As background models are inevitably imperfect, residual errors will consequently propagate into the resulting geodetic products. Accounting for uncertainties of the background model data in a statistical sense, however, has been shown before to be a useful approach to mitigate the impact of residual errors leading to temporal aliasing artefacts. In light of the changes made in the new release RL07 of AOD1B, previous uncertainty assessments are deemed too pessimistic and thus need to be revisited. We here present an analysis of the residual errors in AOD1B RL07 based on ensemble statistics derived from different atmospheric reanalyses, including ERA5, MERRA2 and JRA55. For the oceans, we investigate the impact of both the forced and intrinsic variability through differences in MPIOM simulation experiments. The atmospheric and oceanic information is then combined to produce a new time-series of true errors, called AOe07, which is applicable in combination with AOD1B RL07. AOe07 is further complemented by a new spatial error variance–covariance matrix. Results from gravity field recovery simulation experiments for the planned Mass-Change and Geosciences International Constellation (MAGIC) based on GFZ’s EPOS software demonstrate improvements that can be expected from rigorously implementing the newly available stochastic information from AOD1B RL07 into the gravity field estimation process.

Carbon isotopes in the marine biogeochemistry model FESOM2.1-REcoM3

Abstract. In this paper we describe the implementation of the carbon isotopes 13C and 14C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13CDIC and Δ14CDIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14CDIC.

Subsurface eddies east of the Philippines: Geographic characteristics, vertical structures, volume and thermohaline transport

Subsurface eddies east of the Philippines: Geographic characteristics, vertical structures, volume and thermohaline transport