Thomas Rackow

Effects of Stochastic Ice Strength Perturbation on Arctic Finite Element Sea Ice Modeling

AbstractThe ice strength parameter P* is a key parameter in dynamic/thermodynamic sea ice models that cannot be measured directly. Stochastically perturbing P* in the Finite Element Sea Ice–Ocean Model (FESOM) of the Alfred Wegener Institute aims at investigating the effect of uncertainty pertaining to this parameterization. Three different approaches using symmetric perturbations have been applied: 1) reassignment of uncorrelated noise fields to perturb P* at every grid point, 2) a Markov chain time correlation, and 3) a Markov chain time correlation with some spatial correlation between nodes. Despite symmetric perturbations, results show an increase of Arctic sea ice volume and a decrease of Arctic sea ice area for all three approaches. In particular, the introduction of spatial correlation leads to a substantial increase in sea ice volume and mean thickness. The strongest response can be seen for multiyear ice north of the Greenland coast. An ensemble of eight perturbed simulations generates a spread ...

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.

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

Modeling ENSO with ECHAM6-FESOM: Influence of the ocean resolution

A new climate model supporting multi-resolution meshes in the ocean component has been established at the Alfred Wegener Institute (AWI) in Bremerhaven. The atmospheric component is ECHAM6 with T63L47 setting, while the ocean is simulated by the AWI multi-resolution model FESOM, supporting triangular unstructured meshes. Two multi-century simulations with ECHAM6-FESOM, REF and TRO, document the beneficial role of an increased tropical ocean resolution for ENSO simulations. REF features a tropical ocean resolution of about 1°, TRO employs more than 0.25° in a narrow equatorial band, with resolution gradually decreasing to 1° as in REF. Outside the tropical belt (15°N to 15°S), both meshes are identical. REF and TRO simulate a mean climate comparable to some of the best CMIP5 models. In TRO, however, both the cold tongue SST bias and the western Pacific SST standard deviation bias appear to improve along with the Nino-3 index statistics. Also, advanced ENSO diagnostics including the Nino-3.4 seasonal variance, the annual cycle representation, and its interaction with ENSO tend to improve. The robustness of these improvements is analyzed and their physical explanations are explored.