Dmitry Sidorenko

Designing variable ocean model resolution based on the observed ocean variability

If unstructured meshes are refined to locally represent eddy dynamics in ocean circulation models, a practical question arises on how to vary the resolution and where to deploy the refinement. We propose to use the observed sea surface height variability as the refinement criterion. We explore the utility of this method (i) in a suite of idealized experiments simulating a wind-driven double gyre flow in a stratified circular basin and (ii) in simulations of global ocean circulation performed with FESOM. Two practical approaches of mesh refinement are compared. In the first approach the uniform refinement is confined within the areas where the observed variability exceeds a given threshold. In the second one the refinement varies linearly following the observed variability. The resolution is fixed in time. For the double gyre case it is shown that the variability obtained in a high-resolution reference run can be well captured on variable-resolution meshes if they are refined where the variability is high and additionally upstream the jet separation point. The second approach of mesh refinement proves to be more beneficial in terms of improvement downstream the midlatitude jet. Similarly, in global ocean simulations the mesh refinement based on the observed variability helps the model to simulate high variability at correct locations. The refinement also leads to a reduced bias in the upper-ocean temperature

Long-term ocean simulations in FESOM: evaluation and application in studying the impact of Greenland Ice Sheet melting

The Finite Element Sea-ice Ocean Model (FESOM) is formulated on unstructured meshes and offers geometrical flexibility which is difficult to achieve on traditional structured grids. In this work, the performance of FESOM in the North Atlantic and Arctic Ocean on large time scales is evaluated in a hindcast experiment. A water-hosing experiment is also conducted to study the model sensitivity to increased freshwater input from Greenland Ice Sheet (GrIS) melting in a 0.1-Sv discharge rate scenario. The variability of the Atlantic Meridional Overturning Circulation (AMOC) in the hindcast experiment can be explained by the variability of the thermohaline forcing over deep convection sites. The model also reproduces realistic freshwater content variability and sea ice extent in the Arctic Ocean. The anomalous freshwater in the water-hosing experiment leads to significant changes in the ocean circulation and local dynamical sea level (DSL). The most pronounced DSL rise is in the northwest North Atlantic as shown in previous studies, and also in the Arctic Ocean. The released GrIS freshwater mainly remains in the North Atlantic, Arctic Ocean and the west South Atlantic after 120 model years. The pattern of ocean freshening is similar to that of the GrIS water distribution, but changes in ocean circulation also contribute to the ocean salinity change. The changes in Arctic and sub-Arctic sea level modify exchanges between the Arctic Ocean and subpolar seas, and hence the role of the Arctic Ocean in the global climate. Not only the strength of the AMOC, but also the strength of its decadal variability is notably reduced by the anomalous freshwater input. A comparison of FESOM with results from previous studies shows that FESOM can simulate past ocean state and the impact of increased GrIS melting well.

Can sea surface height be used to estimate oceanic transport variability

The relation between the sea surface height and the meridional transport across a zonal section at 26.5°N in the North Atlantic is studied by using an eddy resolving ocean state estimate simulated with the Massachusetts Institute of Technology general circulation model. It is shown that the correlation between the zonal sea surface height difference and transport can be substantially increased if the steric height contribution from the seasonal thermocline is removed. The latter explains a substantial part of sea surface height variability, but its effect on transport is weak. It is also found that the zonal steric height difference correlates well with the transport after the contribution of the seasonal thermocline has been removed. There is a similar agreement (with correlation coefficient of 0.63 for the full signal and 0.89 for the mean seasonal cycle) between the meridional transport and steric height based on observations from the Rapid Climate Change (RAPID) project.

Towards multi-resolution global climate modeling with ECHAM6-FESOM. Part II: climate variability

This study forms part II of two papers describing ECHAM6-FESOM, a newly established global climate model with a unique multi-resolution sea ice-ocean component. While part I deals with the model description and the mean climate state, here we examine the internal climate variability of the model under constant present-day (1990) conditions. We (1) assess the internal variations in the model in terms of objective variability performance indices, (2) analyze variations in global mean surface temperature and put them in context to variations in the observed record, with particular emphasis on the recent warming slowdown, (3) analyze and validate the most common atmospheric and oceanic variability patterns, (4) diagnose the potential predictability of various climate indices, and (5) put the multi-resolution approach to the test by comparing two setups that differ only in oceanic resolution in the equatorial belt, where one ocean mesh keeps the coarse ~1° resolution applied in the adjacent open-ocean regions and the other mesh is gradually refined to ~0.25°. Objective variability performance indices show that, in the considered setups, ECHAM6-FESOM performs overall favourably compared to five well-established climate models. Internal variations of the global mean surface temperature in the model are consistent with observed fluctuations and suggest that the recent warming slowdown can be explained as a once-in-one-hundred-years event caused by internal climate variability; periods of strong cooling in the model (‘hiatus’ analogs) are mainly associated with ENSO-related variability and to a lesser degree also to PDO shifts, with the AMO playing a minor role. Common atmospheric and oceanic variability patterns are simulated largely consistent with their real counterparts. Typical deficits also found in other models at similar resolutions remain, in particular too weak non-seasonal variability of SSTs over large parts of the ocean and episodic periods of almost absent deep-water formation in the Labrador Sea, resulting in overestimated North Atlantic SST variability. Concerning the influence of locally (isotropically) increased resolution, the ENSO pattern and index statistics improve significantly with higher resolution around the equator, illustrating the potential of the novel unstructured-mesh method for global climate modeling.

Ocean Modeling on a Mesh With Resolution Following the Local Rossby Radius

We discuss the performance of the Finite Element Ocean Model (FESOM) on locally eddy-resolving global unstructured meshes. In particular, the utility of the mesh design approach whereby mesh horizontal resolution is varied as half the Rossby radius in most of the model domain is explored. Model simulations on such a mesh (FESOM-XR) are compared with FESOM simulations on a smaller-size mesh, where refinement depends only on the pattern of observed variability (FESOM-HR). We also compare FESOM results to a simulation of the ocean model of the Max Planck Institute for Meteorology (MPIOM) on a tripolar regular grid with refinement toward the poles, which uses a number of degrees of freedom similar to FESOM-XR. The mesh design strategy, which relies on the Rossby radius and/or the observed variability pattern, tends to coarsen the resolution in tropical and partly subtropical latitudes compared to the regular MPIOM grid. Excessive variations of mesh resolution are found to affect the performance in other nearby areas, presumably through dissipation that increases if resolution is coarsened. The largest improvement shown by FESOM-XR is a reduction of the surface temperature bias in the so-called North-West corner of the North Atlantic Ocean where horizontal resolution was increased dramatically. However, other biases in FESOM-XR remain largely unchanged compared to FESOM-HR. We conclude that resolving the Rossby radius alone (with two points per Rossby radius) is insufficient, and that careful use of a priori information on eddy dynamics is required to exploit the full potential of ocean models on unstructured meshes.

Improved Estimates of the Oceanic Circulation Using the CHAMP Geoid

Ocean general circulation models which are constrained by altimetry data usually assimilate only temporal sea-surface height anomalies. It is known that this is not enough to correct the mean ocean state. Here we present first results of assimilating the full (absolute) dynamical topography into a steady state version of a finite element ocean model (FEOM) for the North Atlantic. This makes it possible to notably reduce the model-data misfit especially in the western part of the basin and in the southern Labrador Sea.

Influence of a Salt Plume Parameterization in a Coupled Climate Model

Sea ice formation is accompanied by the rejection of salt which in nature tends to be mixed vertically by the formation of convective plumes. Here we analyze the influence of a salt plume parameterization (SPP) in an atmosphere-sea ice-ocean model. Two 330 years long simulations have been conducted with the AWI Climate Model. In the reference simulation, the rejected salt in the Arctic Ocean is added to the upper-most ocean layer. This approach is commonly used in climate modelling. In another experiment, employing SPP, the rejected salt is vertically redistributed within the mixed layer based on a power law profile that mimics the penetration of salt plumes. We discuss the effects of this redistribution on the simulated mean state and on atmosphere-ocean linkages associated with the intensity of deep water formation. We find that the salt plume parametrization leads to simultaneous increase of sea ice (volume and concentration) and decrease of sea surface salinity in the Arctic. The SPP considerably alters the interplay between the atmosphere and the ocean in the Nordic Seas. The parameterization modifies the ocean ventilation; however, resulting changes in temperature and salinity largely compensate each other in terms of density so that the overturning circulation is not significantly affected.

The Relative Influence of Atmospheric and Oceanic Model Resolution on the Circulation of the North Atlantic Ocean in a Coupled Climate Model

It is often unclear how to optimally choose horizontal resolution for the oceanic and atmospheric components of coupled climate models, which has implications for their ability to make best use of available computational resources. Here we investigate the effect of using different combinations of horizontal resolutions in atmosphere and ocean on the simulated climate in a global coupled climate model (Alfred Wegener Institute Climate Model [AWI‐CM]). Particular attention is given to the Atlantic Meridional Overturning Circulation (AMOC). Four experiments with different combinations of relatively high and low resolutions in the ocean and atmosphere are conducted. We show that increases in atmospheric and oceanic resolution have clear impacts on the simulated AMOC, which are largely independent. Increased atmospheric resolution leads to a weaker AMOC. It also improves the simulated Gulf Stream separation; however, this is only the case if the ocean is locally eddy resolving and reacts to the improved atmosphere. We argue that our results can be explained by reduced mean winds caused by higher cyclone activity. Increased resolution of the ocean affects the AMOC in several ways, thereby locally increasing or reducing the AMOC. The finer topography (and reduced dissipation) in the vicinity of the Caribbean basin tends to locally increase the AMOC. However, there is a reduction in the AMOC around 45°N, which relates to the reduced mixed layer depth in the Labrador Sea in simulations with refined ocean and changes in the North Atlantic current pathway. Furthermore, the eddy‐induced changes in the Southern Ocean increase the strength of the deep cell.

Sensitivity of deep ocean biases to horizontal resolution in prototype CMIP6 simulations with AWI-CM1.0

CMIP5 models show substantial biases in the deep ocean that are larger than the level of natural variability and the response to enhanced greenhouse gas concentrations. Here we analyse the influence of horizontal resolution in a hierarchy of five multi-resolution simulations with the AWI Climate Model (AWI-CM), which employs a sea ice-ocean model component formulated on unstructured meshes. The ocean grid sizes considered range from a nominal resolution of∼ 1◦ (CMIP5-type) up to locally eddy-resolving. We show that increasing ocean resolution locally to resolve ocean eddies leads to a major reduction 5 in deep ocean biases. A detailed diagnosis of the simulations allows to identify the origins of the biases. We find that two major sources at the surface are responsible for the deep bias in the Atlantic Ocean. Furthermore, the Southern Ocean density structure is equally improved with locally explicitly resolved eddies compared to parameterized eddies. Part of the bias reduction can be traced back towards improved surface biases over outcropping regions, which are in contact with deeper ocean layers along isopycnal surfaces. Our prototype simulations provide guidance for the optimal choice of ocean grids for AWI-CM to be used 10 in the final runs for phase 6 of the ’Coupled Model Intercomparison Project’ (CMIP6) and for the related flagship simulations in the ’High Resolution Model Intercomparison Project’ (HighResMIP). Quite remarkably, retaining resolution only in areas of high eddy activity along with excellent scalability characteristics of the unstructured-mesh sea ice-ocean model enables us to perform the multi-centennial climate simulations needed in a CMIP context at (locally) eddy-resolving resolution with a throughput of 5–6 simulated years per day. 15

Arctic Sea Ice Decline Significantly Contributed to the Unprecedented Liquid Freshwater Accumulation in the Beaufort Gyre of the Arctic Ocean

The Beaufort Gyre (BG) is the largest liquid freshwater reservoir of the Arctic Ocean. The liquid freshwater content (FWC) significantly increased in the BG in the 2000s during an anticyclonic wind regime and remained at a high level despite a transition to a more cyclonic state in the early 2010s. It is not well understood to what extent the rapid sea ice decline during this period has modified the trend and variability of the BG liquid FWC in the past decade. Our numerical simulations show that about 50% of the liquid freshwater accumulated in the BG in the 2000s can be explained by the sea ice decline caused by the Arctic atmospheric warming. Among this part of the FWC increase, 60% can be attributed to surface freshening associated with the reduction of the net sea ice thermodynamic growth rate, and 40% to changes in ocean circulation, which makes freshwater more accessible to the BG for storage. Thus, the rapid increase of the BG FWC in the 2000s was due to the concurrence of the anticyclonic wind regime and the high freshwater availability. We also find that if the Arctic sea ice had not declined, the liquid FWC in the BG would have shown a stronger decreasing tendency at the beginning of the 2010s owing to the cyclonic wind regime. From our results we argue that changes in sea ice conditions should be adequately taken into account when it comes to understanding and predicting variations of BG liquid FWC in a changing climate.