64
SOLA, 2006, Vol. 2, 064 067, doi:10.2151/sola.2006 017
Changes in Precipitation-based Extremes Indices Due to Global Warming Projected by a Global 20-km-mesh Atmospheric Model Kenji Kamiguchi1, Akio Kitoh1, Takao Uchiyama1, Ryo Mizuta2, and Akira Noda1 1 Meteorological Research Institute, Tsukuba, Japan 2 Advanced Earth Science and Technology Organization, MRI, Tsukuba, Japan
Future changes in extremes indices on precipitation were projected with a 20-km horizontal grid atmospheric general circulation model. At the end of the 21st Century, heavy precipitation was shown to increase enormously in South Asia, the Amazon, and West Africa, while a dry spell was shown to increase in South Africa, south Australia, and the Amazon, suggesting that the risk of water-related disasters will be higher in these regions. In the Asian monsoon region, heavy precipitation increases notably in Bangladesh and in the Yangtze River basin due to the intensified convergence of water vapor flux in summer. In the Amazon, a dry spell greatly increases due to a reduction in the Walker circulation caused by an El Ni no-like change in SST prescribed as boundary condition.
1. Introduction It is feared that the frequency and intensity of extreme precipitation events, such as heavy rain or severe drought, will increase in certain areas in the near future due to changes in the strength or patterns of hydrological cycles as global warming progresses. These extreme events may threaten human life and economic growth. Changes in the frequency and intensity of extreme events are more likely to affect society than changes in the mean climate (e.g., Katz and Brown 1992). Therefore, reliable future projections using climate models are in high demand. However, extreme events are difficult to express with low-resolution models because phenomena such as those for heavy precipitation are local and their values are smoothed in a coarse grid. We conducted a global warming experiment using a 20-km-mesh atmospheric general circulation model (AGCM) and analyzed future changes in precipitationbased “extremes indices” proposed by Frich et al. (2002). This very high-resolution model is expected to be superior to low-resolution ones because the model topography is more realistic, so that topographically constrained precipitation can be well resolved.
conducted by the “time-slice” method. AJ is a present climate simulation using the observed climatological sea surface temperature (SST) as the boundary conditions. The monthly mean climatological SST and sea ice concentration are those obtained by Reynolds and Smith (1994), averaged from November 1981 to December 1993. AK is a global warming simulation forced by the climatological SST plus anomalies. The SST anomalies are the difference between the average from 1979 to 1998 in the 20th Century climate simulation and that from 2080 to 2099 in a scenario simulation using the MRI-CGCM2.3 (Yukimoto et al. 2006). The concentrations of greenhouse gases (CO2, CH4, N2O) and aerosols in the scenario simulation with MRI-CGCM2.3 and in the AJ-run are based on the IPCC Special Report on Emission Scenario (SRES) A1B, which assumes CO2 increases about twice in concentration between the periods. The atmospheric part of the CGCM has a horizontal spectral truncation of T42 corresponding to a horizontal grid spacing of about 270 km and has 30 vertical levels. 2.2 Extremes Indices In the future projection of extreme events using a climate model, the intercomparison of different models is very important to estimate the uncertainty of the results on a common basis. Among several statistical methods to diagnose extreme events on precipitation, we focused on “extremes indices” proposed by Frich et al. (2002), as they are widely used in recent studies and adopted as IPCC standard output data for the upcoming AR4 (IPCC 4th Assessment Report). The name, definition, and unit of the indices are listed in Table 1, and more detailed information is found in the Stardex website (http:// www.cru.uea.ac.uk/cru/projects/stardex/). Pav (annual mean of daily precipitation) is not an extremes index, but it is used as a basic index of precipitation to diagnose climate change. For all five extremes indices, higher values indicate more extremes on precipitation. CDD means the length of dry spell, whereas R5d, R10, SDII, and R95T express the intensity or frequency of heavy precipitation. In Frich’s original definition, R95T uses the 95th percentile of daily precipitation among the wet days ( 1 mm/day) from 1961 to 1990. However, the 10-year data
2. Model and extremes indices Table 1. Extremes indices and a basic index on precipitation.
2.1 Model and experiment The model used in this study is an MRI-JMA (Meteorological Research Institute and Japan Meteorological Agency) AGCM. The triangular truncation wave number is 959, with a linear Gaussian grid (TL959) which corresponds to a horizontal 20-km grid resolution. This is currently the highest resolution of AGCM used in global warming experiments. The number of vertical levels is 60, with a model top at 0.1 hPa. Further details of the model are described in Mizuta et al. (2006). Two 10-year runs, referred to as AJ and AK, were
Index
Unit
Definition
Pav
mm/day
Annual mean of daily precipitation (this is not an extremes index)
CDD
days
Consecutive dry days: The annual maximum number of consecutive dry days with < 1 mm/day
R5d
mm
Maximum 5-day precipitation total: The annual maximum consecutive 5-day precipitation total
R10
days
The number of days in a year with precipitation 10 mm/day
SDII
Corresponding author: Kenji Kamiguchi, Meteorological Research Institute, 1-1 Nagamine, Tsukuba 305-0052, Japan. E-mail:
[email protected]. ©2006, the Meteorological Society of Japan.
R95T
Simple daily intensity index: Total annual precipitamm/day tion divided by the number of days with 1 mm/day %
Fraction of annual total precipitation due to extreme wet days exceeding the 95th percentile of the present-day simulation
Kamiguchi et al., Changes in Precipitation-based Extremes Indices Due to Global Warming 100 OBS (30S-30N) OBS (40S-70S, 40N-70N)
-1
10
frequency
MODEL (30S-30N) MODEL (40S-70S, 40N-70N)
10-2
10-3
10-4
0
10
20
30
daily precipitation [mm/day]
40
Fig. 1. Probability distribution functions (PDFs) of daily precipitation computed from the observation (GPCP-1DD; blue lines) and the model (AJ; red lines). The solid lines and the dashed lines denote the tropical/subtropical (30°S 30°N) and the extratropics (40°S 70°S, 40°N 70°N), respectively.
of AJ is used here as a base period representing the present climate. All indices are derived in each year, and the 10-year averages are compared between the present and future climate simulations. The results over the land grids are presented here.
3. Results 3.1 Model evaluation The model performance on reproducing the present climate was evaluated by comparing the model (AJ) and the observed climatology based on GPCP-1DD (Huffman et al., 2001) for the 7-year data from 1997 to 2003. GPCP-1DD, which is a daily precipitation data with 1 × 1 degree for horizontal resolution, is considered as the best data available to compare global daily precipitation at the present moment. However, the quality should be evaluated further because it is based mainly on satellite observations. The probability distribution functions (PDFs) of daily precipitation both in the tropics/subtropics (30°S 30°N) and in the extratropics (40°S 70°S, 40°N 70°N) are shown in Fig. 1. The model precipitation data is regridded to the same grid as GPCP-1DD. The simulated PDF agrees well with the observation in the extratropics, though it somewhat overestimates the frequency of heavy rain larger than 20 mm/day. Meanwhile, in the tropics/subtropics, the PDF of the model overestimates the frequency of weak rain between 5 and 15 mm/day and underestimates the frequency of heavy rain larger than 15 mm/day. The distribution maps of the six indices defined in Table 1 in the observation and AJ are shown in the left and middle panels of Fig. 2, respectively. All six indices simulated by the model are generally in good agreement with those from the observation. However, the simulated R5d, R10, and SDII are somewhat less than the observed ones in the tropics/subtropics, as reflected in the result of Fig. 1. CDD in the model is slightly larger in number than those from the observation in Antarctica, Brazilian plateau, and Mongolia. R95T is the least well simulated with respect to its spatial patterns among six indices; the model R95T is larger in eastern North America and smaller in the Indian Peninsula than the observed one. 3.2. Future changes The right panels of Fig. 2 show the maps of future change in the indices. The dry spell index (CDD) increases in most of the Southern Hemisphere except in
65
Antarctica and decreases in the high latitudes in the Northern Hemisphere, the Arabian Peninsula, and inland China. On the other hand, heavy precipitation indices (R5d, R10, SDII, and R95T) increase significantly in many areas. Since their geographical distribution patterns are similar, R5d is chosen as a typical heavy precipitation index for the following discussion (however, other indices are also shown in Table 2 for the purpose of reference). a. Future changes in 23 regions of the world Table 2 summarizes the extremes indices for AJ and their change (AK-AJ) in 23 regions. The regional classification follows that of Giorgi and Francisco (2000a,b) (the name and map for each region are shown in Supplement-1). Pav increases in many regions except MED, SAF, CAM, and SSA. It increases more in the eastern part of the continents than in the western part. As for R5d, no regions show a decrease in warmer climate, and a large positive change is found in low latitudes, especially in SAS, AMZ, and WAF. The change in R5d is large in the regions where its value is large in the present-day simulation. CDD increases in the low latitudes and in the Southern Hemisphere except in Antarctica, especially in SAF, SAU, and AMZ. Compared with other indices, the number of areas showing statistically significant changes in CDD is very small due to the large year-to-year variability. If we focus on arid regions in which the present value of CDD exceeds 60, CDD increases in WAF and SAF, while it decreases in severely arid regions such as ANT, SAH, and TIB whose present annual mean daily precipitation is less than 0.5 mm/day. These results imply that, in the future climate, both heavy precipitation and drought tend to be more severe in regions where they are severe in the present. Drought, however, tends to mitigate in extremely dry regions in the present. b. Regional changes Regional changes in the extremes indices are briefly discussed here with consideration of the meteorological mechanisms. Two regions, the Asian monsoon region (SAS, EAS, and SEA) and the Amazon (AMZ), were chosen here because large changes in the indices were found there. Figure 3 shows a future changes in R5d and in the water vapor flux at 925 hPa in the Asian monsoon region in June-August (JJA). A large increase in R5d is found in Bangladesh and China’s Yangtze River basin, where the convergence of the water vapor flux is intensified. Due to moisture build-up by global warming, the water vapor flux increases from the Bay of Bengal to the Bangladesh delta (Ashrit et al. 2005). Moreover, an intensified subtropical anticyclone in the western Pacific and an increased atmospheric moisture result in an enhanced water vapor flux which strengthens Meiyu/Baiu rainband from the Yangtze River basin toward Japan (Kitoh et al. 2005). Among all the regions, the Amazon shows the most drastic changes. As shown in Fig. 4, CDD increases remarkably in the south of 5°S and R5d increases in the north of 5°S. The following two points can be considered as reasons for the increases in CDD. 1: The onset of the monsoon delays (the pentad number of the timing of the monsoon onset shifts from 46.5 to 48.2 at 10°S, 63° W (cross mark in the left panel in Fig. 4); the monsoon onset is defined as the first pentad after pentad 44, in which rainfall is larger than 5 mm/pentad). 2: Even weak rainfall (smaller than 5 mm/pentad), which is sometimes observed in the dry season in the present, almost vanishes in the future. In the present simulation, in many cases when rain takes place in the dry season, equatorial easterly low-level wind from ITCZ hits the Andes and deviates to the south, bringing scattered clouds to the middle of the Amazon from ITCZ. In the future climate, the Walker circulation is weak in connection with the El Ni no-like SST change prescribed as boundary condition. Moreover, the anticyclone over the
66
SOLA, 2006, Vol. 2, 064 067, doi:10.2151/sola.2006 017
Fig. 2. Distribution maps of indices defined in Table 1 in the observation (GPCP-1DD; left), present-run (AJ; middle), and their future change (AK-AJ; right). Pav (a1, a2, and a3), CDD (b1, b2, and b3), R5d (c1, c2, and c3), R10 (d1, d2, and d3), SDII (e1, e2, and e3), and R95T (f1, f2, and f3). In the maps of SDII and R95T in the observation and in the present (e1, e2, f1, and f2), an undefined value is not colored. In the maps of future change (a3, b3, c3, d3, e3, and f3), colored areas indicate significance at the 90% confidence level. The value on the map is a global average of the index on land.
southwest Atlantic is also weak. These two changes contribute to the weakening of the equatorial easterly wind and the suppression of rainfall in the middle of the Amazon. On the other hand, the primary reason for increase in R5d is not peculiar to this region. Each rainfall amount simply increases mainly due to an enhancement in water vapor content in the air.
4. Concluding remarks Future changes in precipitation-based extremes
indices were projected using a 20-km-mesh global model. R5d, a heavy precipitation index, increases in almost all regions, especially in the area where the present value is large. CDD, a dry spell index, also generally increases in the area in which the present value is large. However, CDD decreases in areas with extremely small amount of annual precipitation, such as Antarctica, Sahara, and Tibet. These results imply that although there are some arid areas in which the length of the dry spell decreases, heavy precipitation and dryness generally become more severe in regions in which these phenomena are severe at present.
Kamiguchi et al., Changes in Precipitation-based Extremes Indices Due to Global Warming
67
Table 2. Extremes indices and a basic index in the present (upper) and future changes (lower) in each region. The underlined value is not statistically significant at the 90% confidence level. The name and location of regions are shown in Supplement-1.
change (AK-AJ)
present (AJ)
index NAS ALA GRL NEU MED CAS TIB EAS WNA CAN ENA SAH WAF EAF SAF SAS SEA NAU SAU CAM AMZ SSA ANT Pav
1.5
1.5
1.6
2.2
1.3
1.6
3.6
CDD
32.9 27.9
33
20.2 51.7 70.6 65.4 36.4 37.1 18.0 14.1 168.4 79.8 72.2 106.3 56.8 13.4 95.6
1.6
1.1
0.9
2.7
1.9
3.1
3.4
0.3
3.6
3.5
2.3
3.5
7.2
39
38.7 43.1 40.7 119.7
5.6
3.1
0.6
R5d
51.9 47.3 51.4 59.5 66.3 47.2
R10
10.9
9
SDII
4.6
4.1
37
88.9 69.2 99.7 94.1 29.5 84.6 81.7 84.4 104.3 111 79.5 71.1 93.2 101 103.9 22.2
12.1 18.4 17.2
10
6.5
30.2 16.9 35.4 38.5
2.4
38.5 39.6 27.4 43.7 99.1 15.2 16.9 45.8 67.4 36.3
3.3
4.4
5.3
4.1
7.1
5.7
7.1
9.2
2.9
R95T 18.5 17.5 17.6 18.1 19.8 18.4 15.7 20.5 20.1 22.3 21.7 16.6 16.8 16.8 18.8 17.4 15.2 21.8 21.3 17.6 17.2 20.6
9.6
5.2
6.4
6
8.5
8.3
7.1
7.5
8.2
8.8
8
7.2
7.7
8.2
Pav
0.2
0.3
0.3
0.1
0.0
0.1
0.2
0.4
0.1
0.3
0.3
0.1
0.2
0.3
0.1
0.3
0.3
0.2
0
0.1
0.3
0.1
0.1
CDD
4.9
4.1
4.3
0.9
1.6
5.4
11.4
4.1
0.9
0.0
0.7
16.4
4
0.6
7
2.8
2.6
7.9
6.6
3.5
5
1.2
21.1
R5d
6.6
10.4
6.9
7
2.3
5.8
5.7
11.8
3.7
8.3
8.1
10.3 13.2 12.8 10.1 16.3 11.8 12.2
9.4
6.1
15.2
3.8
3
R10
2.7
3.5
3.1
3.1
0.9
1.8
1.7
4.9
1.4
4.2
3.8
1
6.4
7.8
1.1
3.5
6.2
1.9
0.4
0.5
5.7
1.6
1.1
SDII
0.3
0.4
0.4
0.5
0.4
0.4
0.4
0.8
0.3
0.7
0.6
0.6
0.7
0.8
0.5
0.7
0.6
0.7
0.5
0.3
0.8
0.4
0.3
R95T
3.5
5.4
4.8
5
3.2
4.3
5
5.6
3
4.1
3.9
6.6
5.8
6.9
4.6
5.1
4.9
5
4.8
4.3
6
2.2
4.1
NAS: North Asia, ALA: Alaska, GRL: Greenland, NEU: Northern Europe, MED: Mediterranean, CAS: Central Asia, TIB: Tibet, EAS: East Asia, WNA: Western North America, CNA: Central North America, ENA: Eastern North America, SAH: Sahara, WAF: Western Africa, EAF: Eastern Africa, SAF: Southern Africa, SAS: South Asia, SEA: South East Asia, NAU: Northern Australia, SAU: Southern Australia, CAM: Central America, AMZ: Amazon, SSA: Southern South America, ANT: Antarctica
Models” which is a sub-program of the Research Revolution 2002 (RR2002) funded by the Ministry of Education, Culture, Sports, Science, and Technology (MEXT). The authors would like to thank Prof. H. Kanzawa (editor) and two reviewers for their helpful comments.
EAS
SAS
Supplements SEA
Fig. 3. Future change (AK-AJ) in R5d and in the water vapor flux (kg/m/sec) at 925 hPa in JJA around Asia. The areas showing a significant change in R5d at the 90% confidence level are colored.
CDD
R5d
Fig. 4. Future change (AK-AJ) in CDD and in R5d in Amazon. Areas showing significant changes at the 90% confidence level are colored. See the text about the cross mark.
In this study, annual extremes indices were calculated, and their 10-year averages were used for analysis. However, in order to consider the effect of extreme events on society, it is necessary to consider not only the mean change in extreme indices as discussed in this paper but also changes in their year-to-year variability. A long-term integration more than 10 years with annually changing SSTs will be required to discuss the year-to-year variability. Such a topic will be discussed in a separate paper based on another set of the time-slice experiments performed on the Earth Simulator.
Acknowledgments This study was conducted under the “Development of Super High Resolution Global and Regional Climate
1. A map of 23 regions used for analysis is given in the page 1 of Supplement-1. Extremes indices and a basic index of the observation (GPCP-1DD) and the model (AJ, AK) are given in the page 2 of Supplement-1.
Ashrit, A., A. Kitoh, and S. Yukimoto, 2005: Transient response of ENSO-monsoon teleconnection in MRI-CGCM2.2 climate change simulations. J. Meteor. Soc. Japan, 83, 273 291. Frich, P., L. V. Alexander, P. Della-Marta, B. Gleason, M. Haylock, A. M. G. Klein Tank, and T. Peterson, 2002: Observed coherent changes in climate extremes during the second half of the twentieth century. Clim. Res., 19, 193 212. Giorgi, F., and R. Francisco, 2000a: Uncertainties in regional climate change predictions: a regional analysis of ensemble simulation with the HADCM2 coupled AOGCM. Clim. Dyn., 16, 169 182. Giorgi, F., and R. Francisco, 2000b: Evaluation uncertainties in the prediction of regional climate change. Geophys. Res. Lett., 27, 1295 1298. Huffman, G. J., R. F. Adler, M. M. Morrissey, D. T. Bolvin, S. Curtis, R. Joyce, B. McGavock, and J. Susskind, 2001: Global precipitation at one-degree daily resolution from multi-satellite observations. J. Hydrometeor., 2(1), 36 50. Katz, R. W., and B. G. Brown, 1992: Extreme events in a changing climate: Variability is more important than averages. Clim. Change, 21, 289 302. Kitoh, A., M. Hosaka, Y. Adachi, and K. Kamiguchi, 2005: Future projections of precipitation characteristics in East Asia simulated by the MRI CGCM2. Adv. Atmos. Sci., 22, 467 478. Mizuta, R., K. Oouchi, H. Yoshimura, A. Noda, K. Katayama, S. Yukimoto, M. Hosaka, S. Kusunoki, H. Kawai, and M. Nakagawa, 2006: 20-km-mesh global climate simulations using JMA-GSM model. mean climate states . J. Meteor. Soc. Japan, 84, 165 185. Reynolds, R. W., and T. M. Smith, 1994: Improved global sea surface temperature analyses using optimum interpolation. J. Climate, 7, 929 948. Yukimoto, S., A. Noda, T. Uchiyama, S. Kusunoki, and A. Kitoh, 2006: Climate changes of the twentieth through twenty-first centuries simulated by the MRI-CGCM2.3. Pap. Meteor. Geophys., 56, 9 24. Manuscript received 26 September 2005, accepted 9 April 2006 SOLA: http://www.jstage.jst.go.jp/browse/sola/