Pathways and time scales of ocean heat uptake and redistribution in a global ocean-ice model

Author(s):  
Alice Marzocchi ◽  
George Nurser ◽  
Louis Clement ◽  
Elaine McDonagh

<p>Changes in regional ocean heat content are not only sensitive to anthropogenic and natural influences, but also substantially impacted by the redistribution of heat, which is in turn driven by changes in ocean circulation and air-sea fluxes. Using a set of numerical simulations with an ocean-sea-ice model of the NEMO framework, we assess where the ocean takes up heat from the atmosphere and how ocean currents transport and redistribute that heat. Here, the strength and patterns of the net uptake of heat by the ocean are treated like a passive tracer, by including simulated sea water vintage dyes, which are released annually between 1958 and 2017. An additional tracer released in year 1800 is also used to investigate longer-term variability. All dye tracers are released from 29 surface patches, representing different water mass production sites, allowing us to identify when and where water masses were last ventilated. The tracers’ distribution and fluxes are shown to capture years of strong and weak convection at deep and mode water formation sites in both hemispheres, when compared to the available observations. Using this approach, which can be applied to any passive tracer in the ocean, we can: (1) assess the relative role of each of the water mass production sites, (2) evaluate the regional and depth distribution of the tracers, and (3) determine their variability on interannual, multidecadal and centennial time scales.</p>

2020 ◽  
Author(s):  
Katherine Hutchinson ◽  
Julie Deshayes ◽  
Jean-Baptiste Sallee ◽  
Julian Dowdeswell ◽  
Casimir de Lavergne ◽  
...  

<p>The physical oceanographic environment, water mass mixing and transformation in the area adjacent to Larsen C Ice Shelf (LCIS) are investigated using hydrographic data collected during the Weddell Sea Expedition 2019. The results shed light on the ocean conditions adjacent to a thinning LCIS, on a continental shelf that is a source region for the globally important water mass, Weddell Sea Deep Water (WSDW). Modified Weddell Deep Water (MWDW), a comparatively warmer water mass of circumpolar origin, is identified on the continental shelf and is observed to mix with local shelf waters, such as Ice Shelf Water (ISW), which is a precursor of WSDW. Oxygen measurements enable the use of a linear mixing model to quantify contributions from source waters revealing high levels of mixing in the area, with much spatial and temporal variability. Heat content anomalies indicate an introduction of heat, presumed to be associated with MWDW, into the area via Jason Trough. Furthermore, candidate parent sources for ISW are identified in the region, indicating the potential for the circulation of continental shelf waters into the ice shelf cavity. This highlights the possibility that offshore climate signals are conveyed under LCIS. ISW is observed within Jason Trough, likely exiting the sub-ice shelf cavity en route to the Slope Current. This onshore-offshore flux of water masses links the region of the Weddell Sea adjacent to northern LCIS to global ocean circulation and Bottom Water characteristics via its contribution to ISW and hence WSDW properties. </p><p>What remains to be clarified is whether MWDW found in Jason Trough has a direct impact on basal melting and thus thinning of LCIS. More observations are required to investigate this, in particular direct observations of ocean circulation in Jason Trough and underneath LCIS. Modelling experiments could also shed light on this, and so preliminary results based on NEMO global simulations explicitly representing the circulation in under-ice shelf seas, will be presented. </p>


2021 ◽  
Author(s):  
Maurice Huguenin ◽  
Ryan Holmes ◽  
Matthew England

<p>Uptake and storage of heat by the ocean plays a critical role in modulating the Earth's climate system. In the last 50 years, the ocean has absorbed over 90% of the additional energy accumulating in the Earth system due to radiative imbalance. However, our knowledge about ocean heat uptake (OHU), transport and storage is strongly constrained by the sparse observational record with large uncertainties. In this study, we conduct a suite of historical 1972–2017 hindcast simulations using a global ocean-sea ice model that are specifically designed to account for a cold start climate and model drift. The hindcast simulations are initialised from an equilibrated control simulation that uses repeat decade forcing over the period 1962-1971. This repeat decade forcing approach is a compromise between an early unobserved period (where our confidence in the forcing is low) and later periods (which would result in a shorter experiment period and a smaller fraction of the total OHU). The simulations are aimed at giving a good estimate of the trajectory of OHU in the tropics, the extratropics and individual ocean basins in recent decades. Many modelling studies that look at recent OHU rates so far use a simpler approach for the forcing. For example, they use repeating cycles of 1950-2010 Coordinated Ocean Reference Experiment (CORE) forcing that is consistent with the Ocean Model Intercomparison Project 2 (OMIP-2). However, this approach cannot account for model drift. The new simulations here highlight the dominant role of the extratropics, and in particular the Southern Ocean in OHU. In contrast, little heat is absorbed in the tropics and simulations forced with only tropical trends in atmospheric forcing show only weak global ocean heat content trends. Almost 50% of the heat taken up from the atmosphere in the Southern Ocean is transported into the Atlantic Ocean. Two-thirds of this Southern Ocean-sourced heat is then subsequently lost to the atmosphere in the North Atlantic but nevertheless this basin gains heat overall. Our results help to estimate the large-scale cycling of anthropogenic heat within the ocean today and have implications for heat content trends under a changing climate.</p>


2014 ◽  
Vol 27 (5) ◽  
pp. 1945-1957 ◽  
Author(s):  
John M. Lyman ◽  
Gregory C. Johnson

Abstract Ocean heat content anomalies are analyzed from 1950 to 2011 in five distinct depth layers (0–100, 100–300, 300–700, 700–900, and 900–1800 m). These layers correspond to historic increases in common maximum sampling depths of ocean temperature measurements with time, as different instruments—mechanical bathythermograph (MBT), shallow expendable bathythermograph (XBT), deep XBT, early sometimes shallower Argo profiling floats, and recent Argo floats capable of worldwide sampling to 2000 m—have come into widespread use. This vertical separation of maps allows computation of annual ocean heat content anomalies and their sampling uncertainties back to 1950 while taking account of in situ sampling advances and changing sampling patterns. The 0–100-m layer is measured over 50% of the globe annually starting in 1956, the 100–300-m layer starting in 1967, the 300–700-m layer starting in 1983, and the deepest two layers considered here starting in 2003 and 2004, during the implementation of Argo. Furthermore, global ocean heat uptake estimates since 1950 depend strongly on assumptions made concerning changes in undersampled or unsampled ocean regions. If unsampled areas are assumed to have zero anomalies and are included in the global integrals, the choice of climatological reference from which anomalies are estimated can strongly influence the global integral values and their trend: the sparser the sampling and the bigger the mean difference between climatological and actual values, the larger the influence.


2021 ◽  
Author(s):  
Ryan Holmes ◽  
Jan Zika ◽  
Stephen Griffies ◽  
Andrew Hogg ◽  
Andrew Kiss ◽  
...  

<p>Numerical mixing, the physically spurious diffusion of tracers due to the numerical discretization of advection, is known to contribute to biases in ocean circulation models. However, quantifying numerical mixing is non-trivial, with most studies utilizing specifically targeted experiments in idealized settings. Here, we present a precise method based on water-mass transformation for quantifying numerical mixing, including its spatial structure, that can be applied to any conserved variable in global general circulation ocean models. The method is applied to a suite of global MOM5 ocean-sea ice model simulations with differing grid spacings and sub-grid scale parameterizations. In all configurations numerical mixing drives across-isotherm heat transport of comparable magnitude to that associated with explicitly-parameterized mixing. Numerical mixing is prominent at warm temperatures in the tropical thermocline, where it is sensitive to the vertical diffusivity and resolution. At colder temperatures, numerical mixing is sensitive to the presence of explicit neutral diffusion, suggesting that much of the numerical mixing in these regions acts as a proxy for neutral diffusion when it is explicitly absent. Comparison of equivalent (with respect to vertical resolution and explicit mixing parameters) 1/4-degree and 1/10-degree horizontal resolution configurations shows only a modest enhancement in numerical mixing at the eddy-permitting 1/4-degree resolution. Our results provide a detailed view of numerical mixing in ocean models and pave the way for future improvements in numerical methods.</p>


2020 ◽  
Vol 13 (11) ◽  
pp. 5465-5483
Author(s):  
Clément Bricaud ◽  
Julien Le Sommer ◽  
Gurvan Madec ◽  
Christophe Calone ◽  
Julie Deshayes ◽  
...  

Abstract. Ocean biogeochemical models are key tools for both scientific and operational applications. Nevertheless the cost of these models is often expensive because of the large number of biogeochemical tracers. This has motivated the development of multi-grid approaches where ocean dynamics and tracer transport are computed on grids of different spatial resolution. However, existing multi-grid approaches to tracer transport in ocean modelling do not allow the computation of ocean dynamics and tracer transport simultaneously. This paper describes a new multi-grid approach developed for accelerating the computation of passive tracer transport in the Nucleus for European Modelling of the Ocean (NEMO) ocean circulation model. In practice, passive tracer transport is computed at runtime on a grid with coarser spatial resolution than the hydrodynamics, which reduces the CPU cost of computing the evolution of tracers. We describe the multi-grid algorithm, its practical implementation in the NEMO ocean model, and discuss its performance on the basis of a series of sensitivity experiments with global ocean model configurations. Our experiments confirm that the spatial resolution of hydrodynamical fields can be coarsened by a factor of 3 in both horizontal directions without significantly affecting the resolved passive tracer fields. Overall, the proposed algorithm yields a reduction by a factor of 7 of the overhead associated with running a full biogeochemical model like PISCES (with 24 passive tracers). Propositions for further reducing this cost without affecting the resolved solution are discussed.


2021 ◽  
Author(s):  
Marti Florence ◽  
Ablain Michaël ◽  
Fraudeau Robin ◽  
Jugier Rémi ◽  
Meyssignac Benoît ◽  
...  

<p>The Earth Energy Imbalance (EEI) is a key indicator to understand climate change. However, measuring this indicator is challenging since it is a globally integrated variable whose variations are small, of the order of several tenth of W.m<sup>-2</sup>, compared to the amount of energy entering and leaving the climate system of ~340 W.m<sup>-2</sup>. Recent studies suggest that the EEI response to anthropogenic GHG and aerosols emissions is 0.5-1 W.m<sup>-2</sup>. It implies that an accuracy of <0.3 W.m<sup>-2</sup> at decadal time scales is necessary to evaluate the long term mean EEI associated with anthropogenic forcing. Ideally an accuracy of <0.1 W.m<sup>-2</sup> at decadal time scales is desirable if we want to monitor future changes in EEI.</p><p>In the frame of the MOHeaCAN project supported by ESA, the EEI indicator is deduced from the global change in Ocean Heat Content (OHC) which is a very good proxy of the EEI since the ocean stores 93% of the excess of heat  gained by the Earth in response to EEI. The OHC is estimated from space altimetry and gravimetry missions (GRACE). This “Altimetry-Gravimetry'' approach is promising because it provides consistent spatial and temporal sampling of the ocean, it samples nearly the entire global ocean, except for polar regions, and it provides estimates of the OHC over the ocean’s entire depth. Consequently, it complements the OHC estimation from the ARGO network. </p><p>The MOHeaCAN product contains monthly time series (between August 2002 and June 2017) of several variables, the main ones being the regional OHC (3°x3° spatial resolution grids), the global OHC and the EEI indicator. Uncertainties are provided for variables at global scale, by propagating errors from sea level measurements (altimetry) and ocean mass content (gravimetry). In order to calculate OHC at regional and global scales, a new estimate of the expansion efficiency of heat at global and regional scales have been performed based on the global ARGO network. </p><p>A scientific validation of the MOHeaCAN product has also been carried out performing thorough comparisons against independent estimates based on ARGO data and on the Clouds and the Earth’s Radiant energy System (CERES) measurements at the top of the atmosphere. The mean EEI derived from MOHeaCAN product is 0.84 W.m<sup>-2</sup> over the whole period within an uncertainty of ±0.12 W.m<sup>-2</sup> (68% confidence level - 0.20 W.m<sup>-2</sup> at the 90% CL). This figure is in agreement (within error bars at the 90% CL) with other EEI indicators based on ARGO data (e.g. OHC-OMI from CMEMS) although the best estimate is slightly higher. Differences from annual to inter-annual scales have also been observed with ARGO and CERES data. Investigations have been conducted to improve our understanding of the benefits and limitations of each data set to measure EEI at different time scales.</p><p><strong>The MOHeaCAN product from “altimetry-gravimetry” is now available</strong> and can be downloaded at https://doi.org/10.24400/527896/a01-2020.003. Feedback from interested users on this product are welcome.</p>


2013 ◽  
Vol 43 (10) ◽  
pp. 2211-2229 ◽  
Author(s):  
Yu Zhang ◽  
Geoffrey K. Vallis

Abstract Ocean heat uptake is explored with non-eddying (2°), eddy-permitting (0.25°), and eddy-resolving (0.125°) ocean circulation models in a domain representing the Atlantic basin connected to a southern circumpolar channel with a flat bottom. The model is forced with a wind stress and a restoring condition for surface buoyancy that is linearly dependent on temperature, both being constant in time in the control climate. When the restore temperature is instantly enhanced regionally, two distinct processes are found relevant for the ensuing heat uptake: heat uptake into the ventilated thermocline forced by Ekman pumping and heat absorption in the deep ocean through meridional overturning circulation (MOC). Temperature increases in the thermocline occur on the decadal time scale whereas, over most of the abyss, it is the millennial time scale that is relevant, and the strength of MOC in the channel matters for the intensity of heat uptake. Under global, uniform warming, the rate of increase of total heat content increases with both diapycnal diffusivity and strengthening Southern Ocean westerlies. In models with different resolutions, ocean responses to uniform warming share similar patterns with important differences. The transfer by mesoscale eddies is insufficiently resolved in the eddy-permitting model, resulting in steep isopycnals in the channel and weak lower MOC, and this in turn leads to weaker heat uptake in the abyssal ocean. Also, the reduction of the Northern Hemisphere meridional heat flux that occurs in a warmer world because of a weakening MOC increases with resolution. Consequently, the cooling tendency near the polar edge of the subtropical gyre is most significant in the eddy-resolving model.


2015 ◽  
Vol 8 (8) ◽  
pp. 6611-6668 ◽  
Author(s):  
B. Tranchant ◽  
G. Reffray ◽  
E. Greiner ◽  
D. Nugroho ◽  
A. Koch-Larrouy ◽  
...  

Abstract. INDO12, a 1/12° regional version of the NEMO physical ocean model covering the whole Indonesian EEZ has been developed and is now running every week in the framework of the INDESO project (Infrastructure Development of Space Oceanography) implemented by the Indonesian Ministry of Marine Affairs and Fisheries. The initial hydrographic conditions as well as open boundary conditions are derived from the operational global ocean forecasting system at 1/4° operated by Mercator Ocean. Atmospheric forcing fields (3 hourly ECMWF analyses) are used to force the regional model. INDO12 is also forced by tidal currents and elevations, and by the inverse barometer effect. The turbulent mixing induced by internal tides is taken into account through a specific parameterization. In this study we evaluate the model skill through comparisons with various datasets including outputs of the parent model, climatologies, in situ temperature and salinity measurements, and satellite data. The simulated and altimeter-derived Eddy Kinetic Energy fields display similar patterns and confirm that tides are a dominant forcing in the area. The volume transport of the Indonesian ThroughFlow is in good agreement with the INSTANT current meter estimates while the transport through Luzon Strait is, on average, westward but probably too weak. Significant water mass transformation occurs along the main routes of the Indonesian Throughflow (ITF) and compares well with observations. Vertical mixing is able to erode the South and North Pacific subtropical waters salinity maximum as seen in TS diagrams. Compared to satellite data, surface salinity and temperature fields display marked biases in the South China Sea. Altogether, INDO12 proves to be able to provide a very realistic simulation of the ocean circulation and water mass transformation through the Indonesian Archipelago. A few weaknesses are also detected. Work is on-going to reduce or eliminate these problems in the second INDO12 version.


2021 ◽  
Vol 21 (2) ◽  
pp. 1357-1373
Author(s):  
Shigeyuki Ishidoya ◽  
Satoshi Sugawara ◽  
Yasunori Tohjima ◽  
Daisuke Goto ◽  
Kentaro Ishijima ◽  
...  

Abstract. Systematic measurements of the atmospheric Ar∕N2 ratio have been made at ground-based stations in Japan and Antarctica since 2012. Clear seasonal cycles of the Ar∕N2 ratio with summertime maxima were found at middle- to high-latitude stations, with seasonal amplitudes increasing with increasing latitude. Eight years of the observed Ar∕N2 ratio at Tsukuba (TKB) and Hateruma (HAT), Japan, showed interannual variations in phase with the observed variations in the global ocean heat content (OHC). We calculated secularly increasing trends of 0.75 ± 0.30 and 0.89 ± 0.60 per meg per year from the Ar∕N2 ratio observed at TKB and HAT, respectively, although these trend values are influenced by large interannual variations. In order to examine the possibility of the secular trend in the surface Ar∕N2 ratio being modified significantly by the gravitational separation in the stratosphere, two-dimensional model simulations were carried out by arbitrarily modifying the mass stream function in the model to simulate either a weakening or an enhancement of the Brewer–Dobson circulation (BDC). The secular trend of the Ar∕N2 ratio at TKB, corrected for gravitational separation under the assumption of weakening (enhancement) of BDC simulated by the 2-D model, was 0.60 ± 0.30 (0.88 ± 0.30) per meg per year. By using a conversion factor of 3.5 × 10−23 per meg per joule by assuming a one-box ocean with a temperature of 3.5 ∘C, average OHC increase rates of 17.1 ± 8.6 ZJ yr−1 and 25.1 ± 8.6 ZJ yr−1 for the period 2012–2019 were estimated from the corrected secular trends of the Ar∕N2 ratio for the weakened- and enhanced-BDC conditions, respectively. Both OHC increase rates from the uncorrected- and weakened-BDC secular trends of the Ar∕N2 ratio are consistent with 12.2 ± 1.2 ZJ yr−1 reported by ocean temperature measurements, while that from the enhanced-BDC is outside of the range of the uncertainties. Although the effect of the actual atmospheric circulation on the Ar∕N2 ratio is still unclear and longer-term observations are needed to reduce uncertainty of the secular trend of the surface Ar∕N2 ratio, the analytical results obtained in the present study imply that the surface Ar∕N2 ratio is an important tracer for detecting spatiotemporally integrated changes in OHC and BDC.


2020 ◽  
Author(s):  
Till Kuhlbrodt ◽  
Aurore Voldoire ◽  
Matthew Palmer ◽  
Rachel Killick ◽  
Colin Jones

<p>Ocean heat content is arguably one of the most relevant metrics for tracking global climate change and in particular the current global heating. Because of its enormous heat capacity, the global ocean stores about 93 percent of the excess heat in the Earth System. Time series of global ocean heat content (OHC) closely track Earth’s energy imbalance as observed as the net radiative balance at the top of the atmosphere. For these reasons simulated OHC time series are a cornerstone for assessing the scientific performance of Earth System models (ESM) and global climate models. Here we present a detailed analysis of the OHC change in simulations of the historical climate (20<sup>th</sup> century up to 2014) performed with four of the current, state-of-the art generation of ESMs and climate models. These four models are UKESM1, HadGEM3-GC3.1-LL, CNRM-ESM2-1 and CNRM-CM6-1. All four share the same ocean component, NEMO3.6 in the shaconemo eORCA1 configuration, and they all take part in CMIP6, the current Phase 6 of the Coupled Model Intercomparison Project. Analysing a small number of models gives us the opportunity to analyse OHC change for the global ocean as well as for individual ocean basins. In addition to the ensemble means, we focus on some individual ensemble members for a more detailed process understanding. For the global ocean, the two CNRM models reproduce the observed OHC change since the 1960s closely, especially in the top 700 m of the ocean. The two UK models (UKESM1 and HadGEM3-GC3.1-LL) do not simulate the observed global ocean warming in the 1970s and 1980s, and they warm too fast after 1991. We analyse how this varied performance across the models relates to the simulated radiative forcing of the atmosphere. All four models show a smaller ocean heat uptake since 1971, and a larger transient climate response (TCR), than the CMIP5 ensemble mean. Close analysis of a few individual ensemble members indicates a dominant role of heat uptake and deep-water formation processes in the Southern Ocean for variability and change in global OHC. Evaluating OHC change in individual ocean basins reveals that the lack of warming in the UK models stems from the Pacific and Indian basins, while in the Atlantic the OHC change 1971-2014 is close to the observed value. Resolving the ocean warming in depth and time shows that regional ocean heat uptake in the North Atlantic plays a substantial role in compensating small warming rates elsewhere. An opposite picture emerges from the CNRM models. Here the simulated OHC change is close to observations in the Pacific and Indian basins, while tending to be too small in the Atlantic, indicating a markedly different role for the Atlantic meridional overturning circulation (AMOC) and cross-equatorial heat transport in these models.</p>


Sign in / Sign up

Export Citation Format

Share Document