3.1 Historical Simulations of CAOs
From 1979 – 2014, CAOs occurred most frequently across North America and Eurasia. Each of the five climate models were able to reproduce the same general spatial distribution of CAOs as observed with the ERA5, however, each model had a warm bias in the North Atlantic, with the CESM2, WACCM, and MPI having the largest bias (Figure 1). This bias can likely be attributed to how each climate model handles the Atlantic Meridional Overturning Circulation (AMOC; Gent, 2018) or air-sea interactions from fluctuations in Arctic sea ice (Kolstad & Bracegirdle, 2008). Climate model simulations have been shown to underestimate the weakening of the AMOC (Hu et al., 2013; Meehl et al., 2020), which favors more CAOs in the North Atlantic. This may account for the simulation of too few CAOs early in the historical period (R4; Figure A1). The CanESM5 has a cold bias across the western United States and like the MPI, a warm bias across the oceans which is largest in the Southern Pacific. Conversely, the MRI has a large cold bias in the Northern Hemisphere (NH), particularly across the Arctic.
Spatial and temporal similarity were calculated to determine which climate model most accurately simulated the spatial distribution and annual frequency of CAOs for each region (Table 2). While each model was able to simulate the general spatial distribution of CAOs, some regions were better modeled than others (SS; Table 2). Moreover, there were large discrepancies between the time series of mean regional annual CAO days simulated by the climate models and the mean regional annual CAO days from the ERA5 (mean absolute error; MAE; Table2). The CESM2 and MPI had a large warm bias across multiple regions and the largest total error in the SS of annual CAO days. The WACCM (full name: CESM2-WACCM), an extension of CESM2 that models the entire atmosphere (Liu et al., 2010), had less overall bias in SS than the CESM2, followed by theCanESM5 and the MRI. The MRI had the lowest MAE with North America (R1 and R2) while the MPI and WACCM had the lowest MAE in Eurasia. In nearly every region, the model ensemble reduces the errors in spatial similarity (SS) and temporal similarity (MAE).
Table 1: Climate model spatial (SS) and temporal (MAE) historical simulation accuracy (1979-2014). Spatial similarity (SS) - difference between regional mean annual CAO days and the observed annual mean CAO days from the ERA5. The mean absolute error (MAE) is calculated for the annual number of CAO days per region in the historical climate model simulations and the observed (ERA5). Red/blue SS shows where the mean annual CAO days is less than/more than the ERA5. A red/yellow MAE shows where the MAE is large/small. Color intensity of the MAE is relative to the region. Total error is the sum of the absolute values of each column.
|
CESM2
|
WACCM
|
MPI
|
MRI
|
CanESM5
|
Ensemble
|
Region
|
SS
|
MAE
|
SS
|
MAE
|
SS
|
MAE
|
SS
|
MAE
|
SS
|
MAE
|
SS
|
MAE
|
R1
|
-1.0
|
3.7
|
0.2
|
4.1
|
-0.5
|
3.9
|
-0.2
|
3.2
|
0.2
|
3.9
|
-0.3
|
2.8
|
R2
|
-0.2
|
5.2
|
-0.2
|
6.0
|
0.0
|
6.3
|
0.2
|
4.6
|
0.0
|
4.7
|
0.0
|
4.2
|
R3
|
0.3
|
5.4
|
0.1
|
5.5
|
-1.2
|
5.0
|
0.2
|
5.3
|
-0.4
|
5.7
|
-0.2
|
4.3
|
R4
|
-1.7
|
4.3
|
-1.6
|
5.0
|
-1.9
|
5.0
|
-1.0
|
5.6
|
-1.0
|
4.1
|
-1.4
|
4.2
|
R5
|
-0.6
|
5.5
|
0.4
|
4.9
|
0.1
|
5.1
|
0.9
|
5.6
|
0.8
|
6.7
|
0.3
|
4.4
|
R6
|
0.2
|
5.4
|
0.0
|
4.7
|
-0.3
|
4.5
|
0.5
|
5.8
|
-0.1
|
5.5
|
0.1
|
3.9
|
R7
|
-0.8
|
4.6
|
-0.4
|
4.8
|
-0.9
|
4.0
|
-1.0
|
3.3
|
-0.5
|
3.2
|
-0.7
|
3.3
|
R8
|
0.0
|
3.6
|
0.6
|
3.6
|
0.6
|
3.4
|
-0.3
|
2.9
|
0.8
|
3.9
|
0.3
|
2.9
|
R9
|
-1.0
|
2.6
|
-0.2
|
2.6
|
0.3
|
2.5
|
0.6
|
3.3
|
-0.3
|
2.2
|
-0.1
|
2.2
|
R10
|
-0.1
|
1.5
|
0.2
|
2.0
|
-0.6
|
1.6
|
0.1
|
1.5
|
0.3
|
1.7
|
0.0
|
1.3
|
Total Error
|
5.7
|
42.0
|
4.0
|
43.3
|
6.5
|
41.2
|
4.9
|
41.2
|
4.4
|
41.8
|
3.5
|
33.6
|
While there were large discrepancies between the observed and simulated annual number of CAO days (MAE), the spatial distribution of the simulated trends matched the observed trends relatively well (Figure 2). Each model shows the largest decreases in annual CAO days across Northern Hemispheric landmasses. The MPI had the smallest historical trends because the simulation produced too few CAOs early in the historical period and too many late in the period for most places. On the other hand, the MRI has a large negative trend because it produced too many CAOs in the Arctic and western Eurasia early in the historical period. Similar to the observed trends from the ERA5, both the MPI and CESM2 had a neutral to positive trend in CAO days in Eurasia. However, the MPI more accurately replicated the location of this positive trend than the CESM2. Like the ERA5, very few simulated trends in the SH were statistically significant, though the MRI and CanESM5 most accurately simulated the positive trend in CAOs across parts of the Southern Ocean.
3.2 Future Projections of CAOs
Similar to (Vavrus et al., 2006), CAOs are expected to continue decreasing across most of the globe over the next few decades. Compared to the historical period, the ensemble of each SSP shows the mean annual number of CAO days between 2015 and 2054 will decrease between 50% and 100% in most locations (Figure 3). The largest decrease in annual CAO days is in North America and Europe where CAOs have historically occurred most frequently. The CESM2, WACCM, and MRI show a large increase in CAOs across the North Atlantic, consistent with previous studies that have shown a continued weakening of the AMOC in climate model projections (Figure A2; Meehl et al., 2020; Zhang et al., 2019). The MPI and MRI also maintain a relatively large number of mean annual CAO days across North America in all three SSPs. While there are generally fewer annual CAO days with SSP245 and SSP585 than in SSP126, SSP245 and SSP585 do not necessarily result in a larger systematic decrease in CAOs. In the MPI model, more CAOs occur in the Southern Atlantic with SSP245 than SSP126. In the CESM2 model, more CAOs occur in the North Atlantic (R4) from SSP245 than SSP126. SSP585 in the MPI, WACCM, and CESM2 also favor more CAOs in Eurasia (R5 and R6) than in SSP245. Moreover, the WACCM SSP245 simulation shows more CAOs in South America (R9) under than the SSP126 simulation.
Climate models simulate the spatial distribution and trends of CAOs well but are unable to accurately model interannual variability. Though a perfect match is not expected, the large discrepancies between historically simulated and observed annual CAO days indicate the models may be simulating the correct trends for the wrong reasons (Luca et al., 2020). These inaccurate representations of historical climate variability in the models can exacerbate errors in future projections of CAOs (Maraun, 2016). As shown with the historical simulations, an ensemble can be used to reduce the magnitude of individual model error, thus an ensemble is also used for each SSP to better estimate changes in CAOs in each region between 2015 and 2054 (Figure 4).
When compared with the observed annual number of CAO days for each region, the ensemble matches the annual variations and trends well (Figure 4). Only R4 and R5 have particularly poor historical simulations. Climate models have been shown to underestimate variability in R4 (W. M. Kim et al., 2018), which may explain why historical simulations simulated too few CAOs early in the historical period. The complex interaction between and amplified Arctic and surface temperatures in Siberia, which is poorly represented in climate models, may account for much of the discrepancy between annual CAO days simulated in R5 (Cohen et al., 2018; Labe et al., 2020).
Future simulations show a consistent decrease in annual CAO days for most regions with several exceptions. All three SSPs simulate a large increase in annual CAO days between 2030 and 2050 in R4. Though historical simulations for R4 where poor, sea ice melt and a weakening of the AMOC supports the notion that the North Atlantic may be a region of large variability in coming decades. In R1 and R2, future simulations show a slight increase in annual CAO days through 2025 and remaining persistent through 2035 before decreasing to approximately zero annual CAO days by 2054. In R3 (Alaska), historical simulations overestimate the annual number of CAO days early in the historical period and underestimate the annual number of CAO days late in the period which results in an overly negative trend. This suggests the models may be misrepresenting variability in the North Pacific, thus the steady decline in annual CAO days in R3, at least in the near-term, may be off-base. Like R3, historical simulations also underestimated the number of CAO days in R6 (Europe) between 2005 and 2015. Since winter extremes in Europe are heavily reliant on North Atlantic circulation (D. M. Smith et al., 2020), a misrepresentation of variability in the North Atlantic may have caused the discrepancies in observed and simulated CAO days in R6. In South America, annual CAO days remain consistent through 2035 in all SSPs before declining to approximately zero annual CAO days in all but SSP126. Across southern Africa, the already infrequent CAO days are shown to steadily decline to approximately zero annual CAO days by 2035.