Projected Changes in Climate Extremes Using CMIP6 Simulations Over SREX Regions

This paper presents projected changes in extreme temperature and precipitation events by using Coupled Model Intercomparison Project phase 6 (CMIP6) data for mid-century (2036–2065) and end-century (2070–2099) periods with respect to the reference period (1985–2014). Four indices namely, Annual maximum of maximum temperature (TXx), Extreme heat wave days frequency (HWFI), Annual maximum consecutive 5-day precipitation (RX5day), and Consecutive Dry Days (CDD) were investigated under four socioeconomic scenarios (SSP1-2.6; SSP2-4.5; SSP3-7.0; SSP5-8.5) over the entire globe and its 26 Special Report on Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (SREX) regions. The projections show an increase in intensity and frequency of hot temperature and precipitation extremes over land. The intensity of the hottest days (as measured by TXx) is projected to increase more in extratropical regions than in the tropics, while the frequency of extremely hot days (as measured by HWFI) is projected to increase more in the tropics. Drought frequency (as measured by CDD) is projected to increase more over Brazil, the Mediterranean, South Africa, and Australia. Meanwhile, the Asian monsoon regions (i.e., South Asia, East Asia, and Southeast Asia) become more prone to extreme flash flooding events later in the twenty-first century as shown by the higher RX5day index projections. The projected changes in extremes reveal large spatial variability within each SREX region. The spatial variability of the studied extreme events increases with increasing greenhouse gas concentration (GHG) and is higher at the end of the twenty-first century. The projected change in the extremes and the pattern of their spatial variability is minimum under the low-emission scenario SSP1-2.6. Our results indicate that an increased concentration of GHG leads to substantial increases in the extremes and their intensities. Hence, limiting CO2 emissions could substantially limit the risks associated with increases in extreme events in the twenty-first century.


Introduction
Over the last 50 years , precipitation related extremes have resulted in many casualties, and caused enormous socioeconomic damage. The death toll from storms was over 577,000; with a further 58,700 from floods; 650,000 from droughts, and 55,700 from extreme temperatures (WMO 2021). Moreover, the last-two decades were the warmest in the entire 140-year long National Oceanic and Atmospheric Administration (NOAA) climate record (NOAA 2021). In the last 20 years, the globe has experienced many extreme events including heatwaves (e.g., Europe 2003; Russia 2010), floods (e.g., USA 2011; India 2015), fires (e.g., Australia 2019, 2020), flash flooding in the Arabian deserts (e.g., Saudi Arabia 2009, 2011, and coastal flooding (e.g., Venice 2019, Mumbai 2021, Karachi 2020), which caused catastrophic situations in these regions (e.g., Almazroui et al. 2016Almazroui et al. ,2018Philips 2020;Zhang et al. 2020;Herring et al. 2019). Are these events a prelude to what is in store in a warmer world, where more intense and more frequent extremes occur as part of the changing climate? Perhaps this is already evident in 2021, which has seen climate extremes break long-standing records by large margins from east to west and north to south. Antarctica recently calved-off the largest iceberg in the world (Climate and Environment 2021). The deadly recent flooding in central Europe, especially in Germany and Belgium where nearly 200 people lost their lives, surprised the region, and is considered a centennial scale (100-year) extreme event (Samenow 2021;Morris et al. 2021;Cornwall 2021). Parts of the Maritime continent, including Indonesia, faced deadly flooding (df Afonson 2021), as did western India (India monsoon 2021). Recent flooding in eastern and central China (China floods 2021) were, according to Gan and Yeung (2021), stand out events in a 1000 years. Scorching heat waves in Canada and the United States of America (US and Canada heatwave 2021) followed by unusual wildfires, along with heat and fire in the southern Europe and recent fires in Turkey and Greece (Turkey 2021), are likely precursors of what is to be expected in a warmer future world. While attribution studies on those recent extreme events have yet to be completed, previous generations of climate models have long projected strong increases in such extremes. As new cycles of Intergovernmental Panel on Climate Change (IPCC) models become available, it is important to assess their potential to simulate future extreme weather events across the globe. This is what is done in the current study.
The IPCC's Special Report for Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (SREX) indicated a certainty about the increase in the frequency and magnitude of climate extremes IPCC 2012). They report that warm daily temperature extremes, more heavy rainfall episodes, and a decrease in cold extremes during the twenty-first century over several regions of the globe. The occurrence probability of hot extremes and their relationship with rainfall deficits at the global scale was investigated by Mueller and Seneviratne (2012). These extreme events often result in destruction of infrastructure, loss of human lives and livestock as well as damage to natural ecosystems, coastal and land areas (Singh et al. 2021;Batibeniz et al. 2020;Chen et al. 2020). The reliable projections and predictions of extreme events are of great importance for implementing adaptation and mitigation strategies (e.g., Bourdeau-Goulet and Hassanzadeh 2021; Abdullah et al. 2020;Ehsan et al. 2019;Travis 2014;Nõges et al. 2010;Abid et al. 2018), in order to avoid human casualties as well as infrastructure and economic losses of billions of dollars (Smith and Matthews 2015).
Global Climate Models (GCMs) from the Coupled Model Intercomparison Project (CMIP) have been used as a primary tool for examining the past and future changes in the mean and extreme of climate at global and regional scales. A plethora of studies (e.g., Ge et al. 2021;Klutse et al. 2021;Li et al. 2021a;Saeed et al. 2021a, b;Akinsanola et al. 2020;Almazroui 2020;Kim et al. 2020;Rastogi et al. 2020;Srivastava et al. 2020;Wehner et al. 2020;Betts et al. 2018;Diffenbaugh et al. 2017) concluded that under the continuously rising global temperature, catastrophic extreme climatic events are happening with ever greater intensity and frequency across the globe. However, representing the characteristics of climate extremes in the past as well as in projections using multiple GCMs has remained an interesting topic in climate change research (e.g., Chen et al. 2021;Sillmann et al. 2013). For instance, Sillmann et al. (2013) compared the performance of the CMIP phases 3 and 5 (CMIP3 and CMIP5) multi-model ensemble means for temperature and precipitation extremes over different global regions, where not much difference was noted for temperature extremes, while a significant difference was noted for precipitation indices, with the CMIP5 models tending to simulate intense precipitation extremes closer to observations. Recently,  considered a large number of temperature extremes over land regions where they showed that CMIP6 models capture the spatiotemporal variations of the observed temperature extremes well compared to a CMIP5 multi-model ensemble (MME), because of the improved spatial mean and skill of the CMIP6 models relative to CMIP5 models. Generally, there is ample evidence that GCMs exhibit large biases both in temperature and precipitation simulations as compared to observation (e.g., Trenberth 2011), yet it is encouraging to see a gradual improvement in the CMIP models, and that the latest 1 3 Published in partnership with CECCR at King Abdulaziz University generation of models is performing better than the earlier ones in simulating climate extremes (Li et al. 2021b;Kim et al. 2020).
The present study aims to evaluate the projected changes in multiple temperature and precipitation extreme indices based on an MME from 21 GCMs that participated in Phase 6 of CMIP over 26 sub-regions, including the entire global landmass over the SREX regions. The CMIP6 models utilize the latest scenarios, namely Shared Socioeconomic Pathways (SSP) (Eyring et al. 2016). Compared with their ancestors (e.g., CMIP3 and CMIP5), the models in the CMIP6 suite are of high resolution (both in space and time) and have comprehensive physical and dynamical schemes enabling the climate research community to tackle a variety of scientific questions related to increasing global warming and its impacts. CMIP6 datasets were also employed to investigate the large-scale picture of the most intense temperature and precipitation extremes by Li et al. (2021b), and a similar regional breakdown over the 41 sub-regions defined for IPCC Sixth Assessment Report (AR6) was conducted by Kim et al. (2020). The present study is complementary to theirs. As climate extremes impact on human and natural systems in different parts of the globe, particularly in recent times, one of the most important tasks is to understand how climate extremes might behave in future periods under different emission scenarios, both globally and for the SREX regions, using CMIP6 multi-model datasets.

Data and Methodology
The CMIP6 simulations used are from 21 GCMs (Table 1) for historical  and two future time frames representing mid-century (2036-2065) and end-century (2070-2099) for four SSP scenarios i.e., SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5. Each model was given equal weighting with one ensemble member per model. While some of our results are presented as global maps, we also illustrate the spatial variation and inter-model spread for regions defined in the Special Report on Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (SREX; Seneviratne et al. 2012;Vogel et al. 2020). All the extreme indices are calculated at each GCM's native resolution as presented in Table1, followed by a regridding to 1° × 1° using bi-linear interpolation in order to carry out ensemble-based analyses and statistics. All the indices are calculated for each GCM using daily data.
Among the many available climate indices, we focus our analysis here on just four of them, each representing a measure of extreme weather. To determine hot extremes, we compute the annual hottest day (TXx) by calculating the yearly maximum value of the maximum daily temperature for each grid cell. The heat wave frequency index (HWFI) is calculated by counting the number of days for each event per year when daily maximum temperature remains higher than the 90th percentile of the historical period continuously for at least six consecutive days. For wet extremes, we compute the maximum five-day precipitation (RX5day) as the maximum sum of precipitation falling during five consecutive days for each year. Similarly for dry extremes, consecutive dry days (CDD) are computed as the maximum number of consecutive days with precipitation less than 1 mm per day for each historical as well as future 30-year time frame.
For future evaluation of all four extreme indices, spatial plots are presented based on criteria of robustness and significance. The robustness of the projected signal is based on the agreement among different model projections on the direction of changes, where the change in the projected signal, whether an increase or decrease, is considered to be robust if at least 66% of the models agree on the direction of change (Haensler et al. 2013). Meanwhile, the significance of the change signal is measured by employing a two-tailed Student's t test based on equal and unequal variance between the future and historical data for each grid box (Kapic et al. 2011). Equal and unequal variance is determined for each grid box using an F test (Archdeacon 1994). In addition to the spatial plots, changes are presented using box and whisker diagrams representing the spatial variability for each SREX region (see too Vogel et al. 2020). These diagrams are plotted by taking temporal average across all the GCMs for two future time horizons as compared to the reference period for four SSP scenarios, where the boxes and whiskers represent the 66% likely range (17th to 83rd percentile) and very likely range (5th to 95th percentile), respectively. Finally, the temporal distribution across all GCMs for the three-time frames are further analyzed using the kernel density distribution (KDD) over each SREX sub-region. Note that prior to applying KDD, the normal distribution, binomial distribution, and Poisson distribution were applied and found to be not suitable for the data used. The basic theory of KDD can be found in Jones (1990) while the similar but updated kernel density estimation (KDE) is presented by Kamalov (2020). The KDD plots are based on M × N degrees of freedom, where "M" is the number of Models and "N" is the number of years. Here, regional KDD is generated by first taking the regional average over all land points, and then employing a bandwidth of 0.1 for the data of 30 years across 21 GCMs. To explain it further, we first calculate the spatial average over the region of interest for each CMIP6 model and then concatenated data from all GCMs to make a sample that consisted of 630 values (i.e., 21 × 30 = 630, where 21 is the number of models, and each model has 30 values).

Projected Changes in Temperature Indices
The temperature extremes in the form of heat waves or cold spells impacts human health, the biophysical environment, ecosystems, and energy consumption (Hales et al. 2003;McGregor et al. 2005;Frank et al. 2015). These extremes occur at various time scales ranging from daily to monthly and annual scales. Many previous studies (e.g., Brown et al. 2008;Christensen et al. 2007;Meehl et al. 2007;Trenberth et al. 2007;Alexander et al. 2006) focused on temperature extremes and the mechanism of changes in these extremes. As mentioned in the previous section, here we examine the projected changes in the annual hottest day temperature (TXx) and the extreme heat wave frequency (HWFI).

Projections of Annual Hottest Day Temperature (TXx)
Trends in temperature extremes can sometimes be different for the most extreme temperature cases (e.g., TXx) than for less extreme events (e.g., warm days) (Brown et al. 2008;Alexander et al. 2006). The reason for this is that the moderate extremes e.g., warm days/nights are generally computed for each day with respect to the long-term statistics for that day. Thus, an increase in warm days from annual analyses does not necessarily imply warming for the very warmest days of the year . Note that the CMIP6 multimodel extreme TXx index is closer to the observations than is the CMIP5 index over most of the globe (Chen et al. 2020;Kim et al. 2020). Freychet et al. (2021) noticed that the best estimate of extremely hot changes by the end of the century could be worse than previously estimated, mostly for tropical and subtropical regions as well as South and East Asia. Therefore, we first examined the projected changes in the annual hottest day temperature over the entire globe and 26 SREX regions (see Figure S1). Figure 1 shows the projected changes in the TXx, which increase significantly over all regions under all future scenarios as compared to the reference period. TXx is projected to increase more over the higher latitudes as compared to the tropical regions. Except for the Amazon region, most of the tropical regions located in Africa and southeast Asia display less increase in TXx values compared to the higher latitudes. This result is in line with Vogel et al. (2020) who analyzed hot extremes using CMIP6 multimodel present day climate along with different levels of additional global warming (+ 1 °C, + 1.5 °C, + 2 °C, and + 3 °C). The projected increase in TXx scales linearly with the future increase in the concentration of greenhouse gases. Figure 1 shows substantial increase in TXx over the Mediterranean region and adjacent parts of Europe and Central Europe, Northern Africa, and Arabian Peninsula. The projected patterns of TXx are consistent among all future scenarios, and also with the results of the CMIP5 analysis reported by Seneviratne and Hauser (2020). The projected changes in the TXx are more pronounced at the end-century as compared to the mid-century. These changes are more intense in the case of the high emission SSP5-8.5 scenario. The projected changes in TXx are robust over global land areas as the majority of CMIP6 models agree with the sign of change.
The spatial variation of the projected TXx over global land and its 26 SREX regions is shown in Fig. 2. The projected change in the median values of TXx over all the SREX regions increases linearly with an increase in the greenhouse gases for the mid-and end-century. Fig. 1 Spatial distribution of future changes (in °C) in the temperature of the annual hottest day (maximum temperature), TXx, for scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and for the two future time slices (2036-2065 and 2070-2099) as compared with the reference period . The backslash and forward slash represent the grid boxes showing significant and robust change, respectively, hence hatching represents the grid boxes having both significant and robust change. Significance is defined based on a two-tailed Student's t test, while robustness is defined when 66% of all models project a climate change signal in the same direction ◂ 1 3 Published in partnership with CECCR at King Abdulaziz University  The globally averaged TXx is projected to increase by 1.5-2.5 °C by mid-century, and by 2-5.5 °C by the end of the century under the four emission scenarios. The projected TXx reveals smallest spatial variability under SSP1-2.6 and largest under the high emission scenario SSP5-8.5. The spatial variability increases towards the end of the twenty-first century as compared to the near future period. The projected spatial variability in the TXx displays marked regional differences. With some exceptions, such as the Amazon and Northeastern Brazil, most of the tropical regions show small spatial variability as compared to the mid and high latitudes. Moreover, the TXx projections for the end-century show significant differences between different SSP scenarios. For example, the TXx displays an increase in median value of 2.5 and 7 °C under SSP1-2.6 and SSP5-8.5, respectively, over the CEU region. Interestingly, for all SREX regions, the difference in the projected TXx signals for mid-and end-century is negligible for the SSP1-2.6, which can be attributed to its peaked emission/concentration pathway during the twentyfirst century. For other future scenarios (e.g., SSP2-4.5, SSP3-7.0, and SSP5-8.5), the TXx is projected to increase more by the end-century as compared to the mid-century.
The KDD curves show a shift in the TXx distribution to the right side for the future periods under SSP2-4.5 and SSP5-8.5 scenarios as compared to the reference period (Fig. 3). The TXx distributions over SREX regions CAM, CAS, EAF, WAS, and WSA show a similar pattern for future periods as noticed in the reference period. Some other regions (e.g., AMZ, CEU, CGI, ENA, NAS etc.) display flatter distributions in the future period as compared to the reference period. The CMIP6 models over some regions (e.g., CNA, WNA) display double peaks in the TXx distribution in the reference period as well as in the future period, except for SSP5-8.5 for the end-century, where this behavior (i.e., two peaks) disappears. In some regions, the CMIP6 models display a drastic increase in TXx as compared to the reference period. For example, SAH and WAS are regions where CMIP6 models, towards their tails, simulate TXx values up to 65 °C in the future under SSP5-8.5. The CMIP6 models in the Amazon region show the highest value of TXx up to 45 °C in the reference period; whereas, in the future it may reach 50 °C under high emission scenarios. In many regions, based on the CMIP6 models, the tails of the distributions show values exceeding 50 °C in the future. The above discussion reveals that extremely hot days will become even hotter in the future period over almost all regions in the world, while the largest increase is likely over the Mediterranean and adjacent European regions, North America, the Amazon, and South Africa. Figure 4 shows the spatial distribution of HWFI for mid-and end-century. Contrary to the TXx, the HWFI increases more in the tropical regions as compared to the higher latitudes. Again, there is a linear increase in the heat wave days with increasing greenhouse gas concentrations in the future period. As expected, the lowest increase in the HWFI is noticed under SSP1-2.6, while the highest increase occurs under SSP5-8.5. Moreover, the projections show a larger increase in the heat wave days towards the end of the twenty-first century as compared to the reference period. By the end of the twentyfirst century, a large part of the globe shows a drastic increase in heat waves under the high emission scenario SSP5-8.5. The largest increase in the HWFI is projected over the Amazon, Sahara, tropical and southern parts of Africa, Arabian Peninsula, south and southeast Asia, and northern Australia. Interestingly, the projections show more spatial variability in the southern hemisphere as compared to the northern hemisphere. The projected HWFI shows enhanced variability in tropical regions such as CAM, AMZ NEB, WAF and SEA (Fig. 5). The largest increase in the projected spatial variability of HWFI occurs over the SEA region where some grid cells show an increase of more than 300 days under the high emission scenario as compared to the reference period. Globally, the spatial variability also increases towards the end of the century. The HWFI is projected to increase from 50 to 200 days as compared to the reference period, under the high emission scenario SSP5-8.5 (Fig. 5).

Extreme Heat Wave Days (HWFI)
KDD plots indicate that future periods are skewed right and become flatter under all scenarios as compared to the reference period (Fig. 6). The distribution of HWFI in future periods reveals a clear increase in the number of heat wave days in the twenty-first century as compared to the reference period under all scenarios over all SREX regions. This increase is more pronounced for the end-century as compared to mid-century and the reference period. By the end of the twenty-first century, the right tail of the distribution is flatter under SSP5-8.5 in some regions (e.g., AMZ, CGI, EAF, NEB, SAF, SAU, SSA, WAF), suggesting a highly pronounced increase in HWFI compared to the reference period. In most regions, the CMIP6 models simulate an increase of more than 100% in the heat wave days under SSP5-8.5 for the end-century period. For example, in the AMZ region, the KDD curve shows a maximum of 150 heat wave days in the reference period, but this number increases to 350 days under SSP5-8.5 for the end-century (Fig. 6). In some regions (e.g., EAF, NEB, SAF, SSA, WAF etc.) the future projections reveal an increase of more than 200% in the heatwave days as compared to the reference period.
The results indicate a clear increase in the heat wave days at global and regional scales. However, their frequency and intensity vary regionally.

Consecutive 5 days Precipitation (RX5day)
Consecutive precipitation amplification may increase the risk of heavy flooding, imposing heavy costs to socio-economic sectors. Excessive precipitation can lead to heavy floods that may also degrade water quality, harming human  health and ecosystems (Marengo et al. 2020a, b). In addition to flooding, extreme consecutive precipitation also increases the risk of landslides (Marengo et al. 2020a, b). Above-normal consecutive precipitation raises the water table and saturates the ground; hence slopes can lose their stability, causing a landslide. We therefore examined the projected changes in annual maximum consecutive 5-day precipitation (RX5day) under different SSP scenarios (Fig. 7). Note that RX5day is reasonably well reproduced by CMIP6 as compared to the observations, in line with  Fig. 3, but for RX5day where horizontal axis represents values of RX5day (in mm). The scales along x-and y-axis are not same for all panels the CMIP5 result (Chen et al. 2020). However, in the precipitation extreme indices climatology patterns, a substantial reduction of dry bias was noticed over South Asia, South America, and east Asia in CMIP6 as compared with CMIP5 (Kim et al. 2020).
The RX5day is projected to increase significantly over most of the globe. Some regions, such as northern Africa and adjacent parts of Mediterranean, western Arabian Peninsula, Caribbean, northeastern and southwestern strips of South America, South Africa, and large areas in Australia show either no change or a reduction in RX5day as compared to the reference period. The RX5day is projected to increase more in tropical Africa, South Asia, East Asia, and Southeast Asia as compared to other regions in the world (Figs. 7, 8). The SAS, EAS and SEA climates are strongly affected by the monsoon regime. Future projections of monsoon rainfall from CMIP3 to CMIP6 models show an intensification over the Asian region with a corresponding increase in precipitation extremes, which has many implications for socioeconomic sectors in these regions (Ge et al. 2021;Almazroui et al. 2020a;Grose et al. 2020;Narsey et al. 2020;Scoccimarro and Gualdi 2020;Zhou et al. 2014). Like other indices, the projected changes in RX5day precipitation intensify with future increases in the concentration of greenhouse gases over all SREX regions except Australia where negative change over central parts disappear in the case of high emission scenarios. The projected changes are more significant under higher emissions than under lower emission scenarios. The projected spatial variability of RX5day is shown in Fig. 8. The tropical regions display large spatial variability of RX5day as compared to the higher latitudes. Further, the RX5day exhibits large spatial variability in the southern hemisphere as compared to the northern hemisphere. Moreover, the spatial variability in RX5day increases with increasing greenhouse gas concentrations. Spatial variability is projected to increase more under the high emission scenario SSP5-8.5 by the end of the twentyfirst century. The largest spatial variability in RX5day is projected over the SAS region, where some grid-cells may receive over 80 mm more than the reference period. The large spatial variability in RX5day could be related to enhanced land surface heterogeneity and the occurrence of monsoons. The KDD of RX5day displays nearly the same pattern over most of the SREX regions during present and future periods (Fig. 9). However, the SAS region is an exception where the distributions are skewed right for both future periods and scenarios. This further indicates the SAS region is more prone to extreme floods in future under high emission scenarios.

Consecutive Dry Days (CDD)
Similar to heavy floods, persistent drought-like conditions severely impact various sectors of life by influencing agriculture, vegetation, and groundwater levels. Therefore, we further examined the projected changes in the longest streak of consecutive dry days (CDD) (Fig. 10) in the two future time slices (2035-2065 and 2069-2100) as compared with the reference period  for scenarios SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5. The future simulations project an intensification of CDD, with increases over the Mediterranean and adjacent areas of northern Africa and Europe, northwestern Arabian Peninsula, South America, South Africa, and Australia. Decreases are projected over Sahara, southern Arabian Peninsula, Greenland, and northeastern parts of Russia. The Mediterranean region is the climate change hotspot (Giorgi 2006), where south of the Mediterranean the largest increase is noted (more than 80 days as compared to the reference period's episode), while over the European region it reaches up to 30 days, which indicates that more prolonged episodes of droughts may occur over this region in the future. South Africa is another region prone to droughts (Pascale et al. 2020). The Cape Town region in South Africa recently experienced a severe multiyear drought with strong rainfall anomalies (Sousa et al. 2018;Mahlalela et al. 2019). The CMIP6 multi-model mean climate projections also show strong dryness over South Africa (Almazroui et al. 2020b). The projected increase in the CDD over South Africa indicates prolongation of future drought-like conditions over South Africa, including Cape Town. Overall, the hotspot regions for increased hot and dry events are projected to be over the Amazon and Brazil, the Mediterranean, Southern Africa, and parts of Australia. These are in line with the hotspot regions mentioned in Vogel et al. (2020), who used CMIP6 present-day climate and different levels of additional global warming. They also mentioned as hotspots the Sahel region and Indonesia under + 1.5 °C global warming, along with the Andes and Himalayan for + 2.0 °C global warming.
The CDD projections also display strong spatial variability over global land and 26 SREX regions (Fig. 11). CMIP6 models tend to overestimate the CDD over some regions of Africa (e.g., northeast Africa, central east Africa) and South America (e.g., south Amazonia, north South America, northeast South America). However, they underestimate the CDD over Sahara, Southwest Africa, and the Tibetan-Plateau (Kim et al. 2020). In addition, they stated that CMIP6 based CDD is larger than CMIP5 CDD over some African Fig. 10 Spatial distribution of future changes (in number of days) in the longest streak of consecutive dry days (CDD) in the two future time slices (2036-2065 and 2070-2099) as compared with the reference period  for scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5). The dots represent the grid boxes showing robust change which is defined when 66% of all models project a climate change signal in the same direction ◂ Published in partnership with CECCR at King Abdulaziz University sub-regions. The projected CDD spatial variability increases more by the end-century as compared to mid-century. Moreover, the projected variability is higher under the high emission scenario. For the low emission scenario SSP1-2.6, the projected changes in the median values of CDD are nearly the same for mid-and end-century with respect to the reference period. Furthermore, the spatial variability of CDD does not show any remarkable difference between mid-and end-century periods. The above results indicate that the increase in future extreme events can be minimized by limiting the global greenhouse emissions according to the Paris Agreement.

Summary and Conclusions
In this study, we examined projected change in hot and wet extreme events that occur once in a year using a 21-member CMIP6 multi-model daily dataset. The projected changes are analyzed by using four indices for mid-century (2036-2065) and end-century (2070-2099) with respect to the reference period  under four emission scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5) for the entire global land and 26 SREX regions. The results indicate a significant and robust increase of hot temperature and precipitation extremes over most of the land areas. Our results show that the annual hottest day temperature (annual maximum of maximum daily temperature) increases more in the extratropical regions as compared to the tropical areas, while the frequency of extreme heat wave days is projected to increase more in the tropical regions. Consecutive dry days (CDD) are projected to increase more in the Amazon and Brazil, Mediterranean, South Africa, and Australia (see Fig. 10). In general, the hotspot regions for increased hot and dry events are projected over Amazon and Brazil, Mediterranean, Southern Africa and parts of Australia.
On the other hand, the Asian monsoon regions (e.g., South Asia, East Asia, and Southeast Asia) are more prone to extreme flooding, associated with the projected increase in the annual maximum consecutive 5-day precipitation. The projected changes in extremes display large spatial variability within each SREX region. The spatial variability increases with increasing greenhouse gas concentration, and it is highest at the end-century as compared to the mid-century and reference period. As expected, there is no remarkable change in extremes and their spatial variability under the low emission scenario SSP1-2.6. The results indicate that increased concentration of global greenhouse gases leads to substantial increases in the extreme events and their intensities. This is recently evident in the deadly heatwave in the US and Canada in June 2021, where temperatures rose above 49.6 °C in Lytton, British Columbia, indicating that human-caused climate disruption is making extreme weather even worse than predicted (Schiermeier 2021) and that societies must be prepared for more such events in future. Our results further underline that for some extremes (e.g., RX5day and HWFI), the countries of the global south already having high incidence of poverty, will be at disproportionately higher risk. The lowest emission scenario, SSP1-2.6 in our analysis, is projected to take the world to almost 2 °C above pre-industrial levels by the end-century. Lower increases in extremes for SSP1-2.6 analyzed in our study, highlight the importance of stronger mitigation actions to stabilize the global temperature at 1.5 °C in-line with the Paris Agreement goals. The findings of the present study indicate that a stabilization of CO 2 emissions consistent with the aims of the Paris Agreement would substantially reduce the risk of increase in extreme events.