Abstract

The aim of this study was to evaluate temperature and sea-level pressure (SLP) fields and to analyse a related anomalous flow over midlatitudes simulated by the CNRM-CM5.1 global climate model (GCM). Simulated flow over midlatitudes of the Northern Hemisphere was assessed through flow indices, classified into 11 circulation types. Reference data were taken from the NOAA-CIRES 20th Century Reanalysis, version 2c. CNRM-CM5.1 exhibited analogous temperature biases to those reported for the mean of the CMIP5 GCMs’ ensemble. The most prominent features were an erroneous temperature dipole pattern in the Atlantic Ocean and a warm bias over regions of deep water upwelling (locally exceeding 5°C). The latter feature was associated with negative SLP biases in those regions. Too low pressure was found over midlatitudes of the Northern Hemisphere, and CNRM-CM5.1 simulated too frequent zonal flow in these latitudes. The usage of three ensemble members with different initial conditions did not improve model’s outputs because the bias is found to be considerably larger compared to the ensemble members’ spread. The study showed that temperature and SLP biases are connected in certain regions, suggesting that improvement of GCMs and development of bias correction methods should be carried out with a complex insight.

1. Introduction

Projections of a possible future climate and advanced understanding of the climate system largely depend on global climate models (GCMs). Although considerable effort has been made on their improvement over past decades, current generation of GCMs still suffers from errors that lower the credibility of their simulations and bring considerable uncertainties in future climate projections [1, 2]. Wang et al. [3] reported that CMIP5 GCMs tend to underestimate an intensity of Atlantic meridional overturning circulation (AMOC) [4] and linked this deficiency to cold biases in the Northern Hemisphere and warm biases over the Southern Ocean. Majahan [5] suggested that the decreased AMOC is related to an increased sea ice extent in the Arctic and may have contributed to a too slow decline of the Arctic sea ice since 1979 simulated by CMIP5 GCMs [6]. In addition, CMIP5 GCMs are not able to reproduce the changes of a sea ice around the Antarctic properly [7].

Possible connections between these large-scale oceanic factors and climate in distant regions have been widely analysed. Many studies assessed possible links between the observed rapid Arctic sea ice loss and modified atmospheric circulation in the Northern midlatitudes [810], with consequent impacts on wintertime temperatures [11, 12]. European temperatures seem to be also related to sea surface temperature (SST) anomalies in the North Atlantic Ocean [13], and even a more distant connection was found by Zhou and Wu [14] who concluded that heat waves over Eurasia are influenced also by El Niño-Southern Oscillation.

Differences in SST among individual GCMs can explain the majority of the intermodel variability in intensity changes of Walker and Hadley circulation cells [15]. This uncertainty propagates to projections of atmospheric circulation and other meteorological variables, such as precipitation and temperature. Moreover, Shepherd [16] concluded that the effect of increased global temperature on atmospheric circulation is not clear, suggesting that climate change in certain regions may possibly be dominated rather by circulation changes than background warming.

Besides these models’ drawbacks, climate projections are affected by uncertainties originating from the choice of the emission scenario and from the internal variability of climate [17]. Emission scenarios represent possible ways of human society’s development that alter the radiative forcing on climate [18], mainly through a modification of atmospheric chemistry and land-cover changes. In midlatitudes, the uncertainty related to the internal climate variability is primarily associated with atmospheric dynamics, dominated by teleconnection modes [19], such as the North Atlantic Oscillation [20].

The aim of this study was to evaluate global patterns of sea-level pressure (SLP) in relation to a simulated temperature bias. In order to better understand possible dynamic causes of biases in temperature over Northern midlatitudes, we also analysed differences between observed and simulated flow indices over three regions there. Historical simulations of the CNRM-CM5.1 GCM [2] were used, and spatial characteristics of SLP and temperature fields are investigated in three ensemble members (initial condition ensemble), which allows analysing the natural variability simulated by the model. Simulated flow and temperature characteristics are evaluated against the NOAA-CIRES 20th Century Reanalysis, version 2c, which is available since 1851.

2. Data and Methods

2.1. CNRM-CM5.1 Global Climate Model

CNRM-CM5.1 GCM simulations performed by the Czech Hydrometeorological Institute were used. The simulations are available for the extended 1740–2005 period at 6-hour temporal resolution. We used daily means of SLP (averaged 4 daily values) and monthly means of air temperature 2 meters above the surface. The CNRM-CM5.1 GCM was developed by CNRM-GAME (Centre National de Recherches Météorologiques-Groupe d’études de l’Atmosphère Météorologique) and Cerfacs (Centre Européen de Recherche et de Formation Avancée). The GCM contains the ARPEGE-Climat atmosphere component, the NEMO oceanic model, and the GELATO sea-ice component. Three ensemble members of the ARPEGE-Climat atmosphere model with different initial conditions are used. Their horizontal spatial resolution is roughly 150 km (T127) with global coverage, and the model has 31 vertical levels. Although the model performs considerably better compared to its predecessor, substantial errors in seasonal precipitation and cloud radiative forcing are still present, including the double intertropical convergence zone, the critical underestimation of low clouds on the eastern side of the tropical ocean basins, or the lack of cloudiness over the Northern Hemisphere continents [2]. In addition, Huszar et al. [21] revealed too large Arctic sea ice extent in this GCM. More detailed information about CNRM-CM5.1 is available in the study by Voldoire et al. [2].

2.2. NOAA-CIRES 20th Century Reanalysis

Model simulations were evaluated against the NOAA-CIRES 20th Century Reanalysis V2c. It provides a global analysis of the most likely state of the atmosphere since 1851 at 6-hourly temporal resolution, and it is available in the 192 × 94 Gaussian grid [22]. Identical variables (temperature 2 meters above the surface and SLP) were taken in order to evaluate the CNRM-CM5.1 simulations. Seasonal and annual averages of SLP and temperature from CNRM-CM5.1 were interpolated to the aforementioned 192 × 94 Gaussian grid using the ordinary kriging method, allowing a direct comparison with the reanalysis.

2.3. Assessment of Model Biases and Circulation Types

We analysed both annual and seasonal global patterns of SLP and temperature bias for the 50-year-long 1956–2005 period. Moreover, globally averaged temperature biases for various time periods (since 1851) were also calculated. Besides the biases, differences among individual ensemble members were also visualised.

SLP biases were associated with anomalous flows in the simulated data, which were quantified through flow indices (flow strength, direction, and vorticity) [23, 24]. The flow indices were calculated identically for both CNRM-CM5.1 and NOAA-CIRES 20th Century Reanalysis V2c using SLP in 16 points, evenly distributed over analysed regions (North America, Europe, and East Asia; Figure 1).

The flow strength (STR) is a vector sum of zonal () and meridional () flow components (1), and represents the latitude of the centre (50°N). SLP values (hPa) in individual points from Figure 1 are indicated by p1–p16 in the following equations:

The flow direction (DIR) is calculated using an function (2). It is similar to the inverse tangent, except that the signs of both arguments are used to determine the quadrant of the result:

Finally, the total vorticity () is the sum of its zonal and meridional components of vorticity (3):

Based on these indices, 11 circulation types were calculated. If the STR (flow strength) and the absolute value of VORT (total vorticity) are lower than 3, a pressure pattern is unclassified (U). When the absolute value of VORT is at least four times larger compared to STR, a pressure pattern is classified as cyclonic (C, if the vorticity is positive) or anticyclonic (A, if the vorticity is negative). When the aforementioned conditions are not met, the pressure pattern is classified based on DIR (the flow direction) into one of the 8 directional types (N, NE, E, SE, S, SW, W, and NW).

3. Global Patterns of Temperature and Sea-Level Pressure Biases

The negative bias of annual global temperature in CNRM-CM5.1 (approximately 2.8°C) was found during all examined time periods (Table 1), and its magnitude was almost identical among the individual ensemble members.

Spatial patterns of annual temperature in CNRM-CM5.1 for the 1956−2005 period are shown in Figure 2(a). Compared to the NOAA-CIRES 20th Century Reanalysis, the CNRM-CM5.1 ensemble mean captured main features reasonably well, but biases were too pronounced in certain regions (Figure 2(b)). Over land, a negative bias prevails, especially over polar regions (both the Arctic and Antarctic) and high-altitude areas. The magnitude of this bias was roughly 5°C, but it reached 10°C in several grid points over the Tibetan Plateau and the Antarctic. Larger terrestrial areas with a positive bias (1−3°C) were present only over eastern parts of North America. Oceans generally exhibit a cold bias in the Northern Hemisphere and a warm bias in the south. The largest positive bias (around 3°C, locally exceeding 5°C) was found over the eastern Pacific Ocean, the southeastern Atlantic Ocean, and the Southern Ocean and in the Baffin and Hudson Bays. A bias with a similar magnitude but opposite sign was located over northern Pacific Ocean, northern Atlantic Ocean, and Arctic Ocean. The described overall temperature biases were similar for all ensemble members, especially in the Tropics, where the difference between the “warmest” and “coldest” ensemble members was only 0.1–0.3°C (Figure 2(c)). The agreement among the ensemble members was somewhat weaker in the midlatitudes, but differences larger than 1°C were still quite rare. The largest variability was found over polar regions, mainly in the Arctic where the difference between the “warmest” and “coldest” ensemble members exceeded 3°C in the Barents Sea.

Turning to SLP, CNRM-CM5.1 was able to simulate locations of atmospheric action centres (Figure 3(a)) quite well; however, some notable biases were found. Over land, negative SLP bias was located in midlatitudes of the Northern Hemisphere (analysed in more detail in Section 4) and in the Antarctic (Figure 3(b)). The most distinctive negative SLP biases over oceans were found in the eastern Pacific Ocean and the southeastern Atlantic Ocean, and the biases were related to overestimated temperature in these regions. By contrast, positive SLP bias was found over western parts of the United States, Mexico, the Tibetan Plateau, the Arctic Ocean, and the Southern Ocean. Differences between ensemble members with highest/lowest SLP (in individual grid points; Figure 3(c)) are analogous to those obtained from temperature output (Figure 2(c)). The best agreement between individual ensemble members was found in the Tropics (especially over the Indian Ocean), while the largest differences were over polar regions and in Azores and North Pacific anticyclones.

Focusing on individual seasons, the smallest global temperature bias was found during the boreal summer (−2.4°C in the 1980–2005 period; Table 2). In the remaining seasons, the biases were somewhat larger (roughly −2.9°C) and quite similar. Analogously to the annual bias, seasonal biases were stable during different time periods and comparable among individual ensemble members.

The global pattern of temperature bias varied across individual seasons (Figure 4). The substantial negative bias in Greenland and Arctic, which was especially pronounced in boreal spring and winter (Figures 4(a) and 4(d)), almost disappeared during boreal summer (Figure 4(b)). Analogously to the Arctic, the strong negative bias over the Antarctic was reduced in austral summer (Figure 4(d)). The aforementioned positive bias over eastern North America was most pronounced during the boreal autumn, when its values exceeded 5°C in some grid points. The oceanic temperature bias dipole was least distinctive during boreal autumn. On the contrary, the eastern Pacific Ocean and the southeastern Atlantic Ocean regions with warm biases were relatively stable during all seasons.

No considerable differences in the magnitude of the SLP bias were found among individual seasons, except for the Antarctic Ocean and the Southern Ocean (Figure 5). The underestimated SLP in Northern midlatitudes and over the upwelling regions in the eastern Pacific Ocean and the southeastern Atlantic Ocean was also fairly stable during the year. The positive SLP bias over the western parts of the United States, Mexico, and the Tibetan Plateau was a year-round feature; however, it was least distinctive during the boreal summer.

4. Northern Hemispheric Wintertime Circulation Types over Midlatitudes

In this section, we evaluated the capability of CNRM-CM5.1 to simulate seasonal frequencies of flow types over three midlatitude regions (North America, Europe, and East Asia). Over North America during winter and spring, the most frequent northerly advection (occurring on about one-fifth of days) was underestimated in all three model integrations, in which the frequency of the northern circulation type varied between 15% and 17% (Figure 6(a)). The occurrence of the types with advection from northeast was also suppressed in the simulations. By contrast, advection from west and northwest and the frequency of the cyclonic type were overestimated, indicating a too frequent zonal flow in simulated data during these seasons. These synoptic-scale anomalies can be clearly linked to biases in the mean SLP gradient over North America, that is, to the positive SLP bias over western parts of the United States and Mexico and the negative bias over Alaska and Canada (Figures 5(a) and 5(d)). Furthermore, these errors in atmospheric dynamics can be linked to positive winter temperature biases (Figure 4(d)), since the types that are associated with the advection of relatively cold Arctic air or the development of cold continental stagnant air occur with lower frequency in the model outputs. The underestimation of northerly types was found also during summer and autumn. By contrast, CNRM-CM5.1 simulated too frequent easterly, south-easterly, and southerly advection in these seasons (Figure 6(a)). The largest temperature anomaly during autumn was not clearly related to biases in large-scale flow, and this error would be possibly associated with land-atmosphere interactions.

Too frequent westerly advection was also found over Europe (except during the summer season; Figure 6(b)). Its observed frequency during winter (20%) was considerably overestimated in CNRM-CM5.1, ranging from 25% to 27%. This overestimation is less pronounced in the transient seasons (spring and autumn), and the summertime frequency of the westerly flow type is underestimated in CNRM-CM5.1. The model also simulated too frequent south-westerly advection, while the frequency of easterly and south-easterly circulation types was underestimated (except in summer). Interestingly, these errors did not result in a notable wintertime positive temperature bias over Europe, probably due to underestimated temperature in the North Atlantic Ocean in the model simulations.

Over East Asia, CNRM-CM5.1 simulated the prevailing wintertime monsoon-related northerly advection quite well (Figure 6(c)); however, the frequencies of the northwestern and western circulation types were markedly overestimated at the expense of the northeastern circulation type. The latter is linked to the bias in the zonal SLP gradient (Figure 5(d)). By contrast, CNRM-CM5.1 failed to simulate the prevailing southerly advection associated with monsoon during summer. The frequency of the southerly type was almost halved in CNRM-CM5.1 compared to the reanalysis.

In all, too frequent zonal flow was found across the Northern midlatitudes. This bias is clearly associated with overestimated/underestimated SLP in two belt-shaped patterns, one between about 20 and 30°N and the other around about 55°N (Figure 5(d)). Furthermore, the underestimated SLP in midlatitudes is linked to the overestimated (underestimated) frequency of strongly cyclonic (anticyclonic) circulation types in all examined regions (Figure 6).

5. Discussion

CNRM-CM5.1 exhibited similar temperature biases as have been shown for the CMIP5 multimodel mean, with respect to their spatial patterns and magnitude [25]. Pronounced warm biases were found over upwelling regions in the southeastern Atlantic Ocean and Pacific Ocean. This bias was already reported for the CMIP3 GCMs, and it was linked to erroneous southward displacement of the intertropical convergence zone, subsequent modification of westerly winds, and deepened thermocline, which reduce upwelling of cold water. In addition, the bias is related also to improper simulation of low-altitude clouds in these regions [26]. An analogous issue was reported for CMIP5 simulations [27], and Wang et al. [3] revealed other interhemispheric links between these warm biases and distant large-scale factors, such as weaker monsoons in the Northern Hemisphere. Therefore, a reduction of these biases cannot be based on correction of regional processes only but has to be carried out with a complex insight.

Another notable deficiency of CMIP5 GCMs (including CNRM-CM5.1) is the oceanic bias dipole, probably originating from a weakened northward heat transport due to a reduced AMOC [3], resulting in the warmer Southern Ocean and the colder North Atlantic Ocean compared to the reference data. The cold bias over the North Atlantic may possibly be related to the persisting negative phase of Atlantic multidecadal oscillation [28] in the model outputs. CMIP3 GCMs, however, have difficulties to simulate this feature, mainly due to a low signal-to-noise ratio and a long period of oscillation [29], and it is possible that CNRM-CM5.1 suffers from the same deficiency as well. The negative temperature bias in northern Atlantic and Europe may also possibly be linked to improper simulation of Gulf Stream, which carries warm water erroneously towards the Canadian Archipelago with consequent impacts on the simulated North American climate. Huszar et al. [21] reported an overestimated sea ice extent in the CNRM-CM5.1, which may explain the strong cold bias over the Arctic Ocean in spring and winter found in our study.

Besides polar regions, considerable negative temperature biases over land were found over mountainous areas, which is in accordance with IPCC [25] and Su et al. [30] who evaluated the performance of CMIP5 GCMs over the Tibetan Plateau and suggested that the negative bias may be linked to snow-cover feedback processes in the models. Besides these physical processes, however, the bias may be affected by a simple representation of orography or a lack of observed data in these remote areas that affects the performance of the reanalysis. In addition, it should be noted that the magnitude of biases may be influenced by a selection of reference data (i.e., the NOAA-CIRES 20th Century Reanalysis V2c), inasmuch as Stryhal and Huth [31] showed that the choice of different reanalyses can have a profound effect on GCM validation over Europe in winter.

Temperature differences between individual ensemble members are relatively small on the global scale, which is in accordance with Kharin et al. [32], who showed that a spread between individual GCMs was larger compared to the difference between ensemble members of one particular GCM. In our study, the difference between the ensemble members in simulating temperature was smallest over equatorial areas and gradually increased poleward. The largest discrepancies were found over polar oceans, in accordance with Kharin et al. [32], and may be related to different boundaries of sea ice extent in individual ensemble members that would considerably alter a surface climate.

Small differences among the individual ensemble members were found also when assessing differences in sea-level pressure fields (except polar regions) and frequencies of circulation types over North America, Europe, and East Asia. The common deficiency in these three regions was an overestimated frequency of zonal flow, which is probably associated with the known issue of GCMs not being able to correctly simulate the occurrence of atmospheric blocking [33, 34].

6. Conclusions

We evaluated temperature and SLP biases and frequency of circulation types in three historical ensemble members of CNRM-CM5.1 GCM against the NOAA-CIRES 20th Century Reanalysis. The main conclusions are summarised as follows:(i)The CNRM-CM5.1 model had a relatively good skill in reproducing global temperature patterns; however, it exhibited drawbacks similar to those reported for other CMIP5 GCMs. An oceanic bias dipole pattern is present in CNRM-CM5.1, resulting in a too warm Southern Ocean and an excessively cold North Atlantic Ocean and Arctic Ocean.(ii)Another deficiency shared with the majority of CMIP5 models is a positive temperature bias in upwelling regions along western coasts of Africa and Southern America (locally exceeding 5°C). This temperature bias was associated with too low sea-level pressure in these regions.(iii)The most apparent biases were present in all three ensemble members, regardless of the member’s initial conditions, and, therefore, are rather a result of model errors than simulated natural variability. This holds not only for mean seasonal and annual temperature and SLP patterns but even for the considerably more chaotic and variable synoptic-scale circulation over three northern midlatitude regions.(iv)CNRM-CM5.1 overestimated the frequency of zonal flow over North America, Europe, and East Asia. The considerable positive wintertime temperature bias in North America was associated with suppressed northerly advection of cold air from the Arctic in CNRM-CM5.1 and was found in all ensemble members.

Data Availability

The outputs of the CNRM-CM5.1 global climate model can be obtained from the ESFG data node (https://esg-dn1.nsc.liu.se/projects/esgf-liu/). NOAA-CIRES 20th Century Reanalysis V2c is available through the ESRL website (https://www.esrl.noaa.gov/psd/data/gridded/data.20thC_ReanV2c.html).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was supported by the Ministry of Education, Youth and Sports of the Czech Republic within the National Sustainability Program I (NPU I), Grant no. LO1415. The authors acknowledge the CNRM-CM5.1 global climate model developed by CNRM-GAME (Centre National de Recherches Météorologiques-Groupe d’études de l’Atmosphère Météorologique) and Cerfacs (Centre Européen de Recherche et de Formation Avancée). The authors also acknowledge the Twentieth Century Reanalysis Project version 2c dataset that is provided by the U.S. Department of Energy, Office of Science Biological and Environmental Research (BER), and by the National Oceanic and Atmospheric Administration Climate Program Office. The authors thank Jan Stryhal, Institute of Atmospheric Physics, Prague, for his valuable comments on the draft of the paper.