Next Article in Journal
Microcosm Experiment to Assess the Capacity of a Poplar Clone to Grow in a PCB-Contaminated Soil
Previous Article in Journal
Comparison of Meteorological- and Agriculture-Related Drought Indicators across Ethiopia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Future Changes in Precipitation and Drought Characteristics over Bangladesh under CMIP5 Climatological Projections

1
Department of Agricultural Engineering, Institute of Agriculture and Life Science, Gyeongsang National University, 52828, 501, Jinju-daero, 13557 Jinju, Gyezongnam, Korea
2
Bangladesh Rice Research Institute, Gazipur 1701, Bangladesh
3
Research Department, APEC Climate Center, Busan 48058, Korea
*
Author to whom correspondence should be addressed.
Water 2019, 11(11), 2219; https://doi.org/10.3390/w11112219
Submission received: 18 September 2019 / Revised: 14 October 2019 / Accepted: 22 October 2019 / Published: 24 October 2019
(This article belongs to the Section Hydrology)

Abstract

:
The impacts of climate change on precipitation and drought characteristics over Bangladesh were examined by using the daily precipitation outputs from 29 bias-corrected general circulation models (GCMs) under the representative concentration pathway (RCP) 4.5 and 8.5 scenarios. A precipitation-based drought estimator, namely, the Effective Drought Index (EDI), was applied to quantify the characteristics of drought events in terms of the severity and duration. The changes in drought characteristics were assessed for the beginning (2010–2039), middle (2040–2069), and end of this century (2070–2099) relative to the 1976–2005 baseline. The GCMs were limited in regard to forecasting the occurrence of future extreme droughts. Overall, the findings showed that the annual precipitation will increase in the 21st century over Bangladesh; the increasing rate was comparatively higher under the RCP8.5 scenario. The highest increase in rainfall is expected to happen over the drought-prone northern region. The general trends of drought frequency, duration, and intensity are likely to decrease in the 21st century over Bangladesh under both RCP scenarios, except for the maximum drought intensity during the beginning of the century, which is projected to increase over the country. The extreme and medium-term drought events did not show any significant changes in the future under both scenarios except for the medium-term droughts, which decreased by 55% compared to the base period during the 2070s under RCP8.5. However, extreme drought days will likely increase in most of the cropping seasons for the different future periods under both scenarios. The spatial distribution of changes in drought characteristics indicates that the drought-vulnerable areas are expected to shift from the northwestern region to the central and the southern region in the future under both scenarios due to the effects of climate change.

1. Introduction

Nowadays, climate change is regarded as a major global issue, and it poses significant challenges to human existence and socio-economic development, particularly in Bangladesh. The dominant features of climate change in Bangladesh detected in the late 19th century include the significant increases in temperature and monsoon and post-monsoon precipitation due to global warming. However, a decreasing trend of precipitation was observed in the winter season [1]. The overall increasing rate of warming in Bangladesh is higher than the present increasing rate of global warming, and this trend is expected to continue over the 21st century [2,3,4]. Notably, Bangladesh has been suffering from frequent natural disasters. Among the most expensive natural disasters, drought is a chronic natural disaster that can have severe and long-lasting impacts on water resources, agriculture, ecosystems, and human societies [5]. The impact of drought events is often worst in developing countries because their economies are driven mainly by agricultural products that are adversely affected by meteorological droughts [6]. Drought may generally be defined as a scarcity of water in a region over a prolonged period of time triggered by a lack of rainfall. Presently, great effort needs to be expended on researching future changes in rainfall patterns, which are the major causes of drought in Bangladesh and can lead to adverse changes in economic and social development. The effects of gradual climate changes and extreme weather events may negatively impact overall socio-economic development in many regions, and therefore, the scientific community and policymakers need more information about the probability of future occurrences of such events [7,8].
The Coupled Model Intercomparison Project (CMIPs) has made available general circulation model (GCM) outputs for the Programme for Climate Model Diagnosis and Intercomparision (PCMDI), and these products are available for research. The climate model outputs from Phase 3 of the CMIP (CMIP3) were broadly used in the Intergovernmental Panel on Climate Change (IPCC) 4th Assessment Report (AR4) [9]. Recently, climate models from Phase 5 of the CMIP (CMIP5) along with greenhouse gas concentration scenarios termed as representative concentration pathways (RCPs) were adopted by the IPCC for its 5th Assessment Report (AR5) [10]. The models from CMIP5 joined with the RCP scenarios have delivered more precise representations of climate outputs than the CMIP3 model results because corrections were made in regard to some key assumptions of climate that were overlooked previously by the model developers [11]. Sperber et al. [12] demonstrated that the CMIP5 models are more competent for capturing numerous features of the Asian monsoon climate compare to the CMIP3 models. Therefore, it is of supreme importance to assess future changes in climatological drought by using the precipitation data from the new sets of CIMP5 GCM projections.
Presently, GCMs are the principal tools for predicting and projecting future climate changes. However, most such global climate models are typically run at coarse resolutions, e.g., more than hundreds of kilometers. Therefore, the GCM outputs are inherently unable to represent regional or local climate features and dynamics at the necessary spatial resolutions for detailed analyses [13]. To overcome this problem, downscale techniques have been developed to obtain local climate change information at the desired scale from coarser-resolution GCM outputs [14,15]. Dynamic downscaling data have been recognized to be more representative of fine-scale physical processes than statistical downscaling data; however, the former technique requires more expensive computing resources than the latter technique [16,17,18]. A large number of simulations with multiple GCM configurations and emissions scenarios can be computed efficiently by using a statistical approach, and such an approach is well accepted in the scientific community; these types of data are widely used in downscaling climate projections. Thus, projections from multiple GCMs from CMIP5 have been used in this study.
In parallel, climate change impact assessments at the regional scale rarely use the raw GCM outputs because climate model data suffer from systematic biases due to the uncertainty in the parameterization of unsolved processes [19]. Therefore, bias-corrected GCM outputs are essential for regional climate impact studies and vulnerability assessments. A vast number of bias correction procedures are in use, such as the monthly mean correction [20], delta change [21], and quantile mapping [22] techniques. The quantile mapping methods are considered to be the most accurate methods in terms of precipitation [23] among all of the other methods. Therefore, a quantile based bias correction approach was used to adjust for the model biases in this study.
One way of assessing the changes in future drought characteristics is to use climate projections from a couple of GCM simulations under different greenhouse gas emission scenarios [24]. According to Coelho and Goddard [25], future changes in precipitation patterns may aggravate drought risk in highly vulnerable tropical areas. In several areas, climate change is expected to primarily affect precipitation, and thus, it is likely to affect the frequency and severity of metrological droughts. A remarkable number of studies have assessed the future drought risks on global and regional scales by using multiple climate model scenarios [26,27,28,29]. For instance, Chen et al. [30] investigated the future change in the drought pattern over the 21st century in China by using global climate models and regional climate models (RCMs) under the SRES (Special Report on Emissions Scenarios) A1B scenario, and the findings indicated that droughts will become less frequent in most areas of China. Conversely, Wang et al. [31] examined future drought in China by using CMIP5 model outputs under the RCP4.5 and RCP8.5 scenarios and revealed that extreme drought events would increase in the future. This discrepancy may have been due to the limitations of drought indices and the lack of realistic climate information. Importantly, the future drought outlook may depend on the indicators used in drought calculations along with realistic climate data. The selection of appropriate drought indices is therefore vital for understanding the future drought characteristics and for planning a drought mitigation strategy for a region [32]. Hence, it is necessary to predict future drought characteristics on a local scale by using realistic climate information and locally applicable and appropriate drought indices.
In the past few years, several studies have projected the average temperatures and annual rainfall for Bangladesh in the future [33,34,35]. For instance, Rahman et al. [34] estimated the annual rainfall by using version 3 of an RCM (RegCM3), and the findings indicated that a 50% decrease in rainfall will occur by 2060. Nowreen et al. [33] found that the amount of total rainfall is likely to increase in the future based on simulations with a 17-member ensemble-driven by the Hadley Centre Coupled Model (HadCM3). A number of studies have suggested that the overall monsoon rainfall will increase and post-monsoon rainfall will decreases in the future throughout most parts of the country [34,35]. Notably, very few studies have been done concerning future drought projections [35,36]. Islam et al. [36] evaluated the drought hazards at current and future climate change conditions in the western region of Bangladesh by using simulated climate data from the outputs of three global climate models for the period between 2041 and 2070. Hasan et al. [35] estimated the future drought conditions by using 25 km high resolution downscaled and projected climate data generated from the RCM known as PRECIS for a continuous period of 1971–2100. The results revealed that droughts will generally decrease in future years but there will be a comparative higher frequency of droughts in the mid 21st century. However, all of these studies used the climate model output from CMIP3 and were driven by the previous IPCC AR4 scenarios.
Recently, Hasan et al. [37] projected the future climate and associated extremes while considering the new RCP4.5 and RCP8.5 scenarios by using RCM results driven by the GCMs over Bangladesh within the new CMIP5. The study revealed that overall precipitation and temperature trends are likely to increase in the future over this region. However, drought characteristics were not evaluated by using these outputs. Mortuza et al. [38] projected only the future drought frequency from 2020 to 2100 by using bias-corrected four GCM (CMIP5) model outputs under the RCP4.5 and RCP8.5 scenarios and showed that the drought frequency will decrease in the future (2020–2100) compared to the past (1961–2010). They also pointed out that more frequent and severe droughts will occur on the west side of the country. However, a complete assessment of future changes in drought characteristics in terms of the frequency, duration, and intensity of meteorological drought events and seasonal drought days mainly based on CMIP5 multimodel simulated data still has not been performed to understand the future conditions of droughts in Bangladesh, and this was the prime motivation for this study. The downscaling techniques of spatial disaggregation were utilized in this study to provides finer-resolution climate projections.
Over the year, a good number of drought indices (DIs) have been proposed and developed to identify the spatiotemporal patterns of droughts and quantify their intensity. The maximum number of drought indices is region-specific, as indices are limited in terms of their applications to different climatic conditions due to the inherent complexity of drought phenomena. In Bangladesh, most of the studies have measured drought severity by using the Standardized Precipitation Index (SPI; [39]), while very few research studies have used other DIs. The SPI is calculated based on the averaged monthly precipitation for a certain period. Therefore, the time steps involved in the SPI tend to produce several different values for the same period. Additionally, the SPI does not take into account the water resources generated by rainfall that may have already been lost due to outflow as well as the effect of evaporation. Moreover, the SPI tends to assign equal weight to temporally different precipitation events, thereby resulting in inaccuracies in predictions of the drought severity [40]. Byun and Wilhite [41] developed a new series of indices to overcome these limitations. Specifically, they used a new concept of effective precipitation (EP), in which the EP represents the summed value of all daily precipitation with a time-dependent reduction function. The EDI was found to be more responsive to drought conditions compared with other drought indices, and it could capture the real essence of the meteorological drought situation in the study area [42,43,44,45]. In particular, Kamruzzaman et al. [46] empirically demonstrated the superiority of the EDI over the SPI when monitoring both long-term and short-term droughts in Bangladesh. Thus, the EDI was employed in this study to interpret the changes in drought patterns to gain knowledge about the potential future changes in drought over Bangladesh.
With the aim of investigating the future changes in precipitation and meteorological drought characteristics over Bangladesh, the drought characteristics were assessed by using the EDI and the bias-corrected CMIP5 GCM precipitation data under the RCP4.5 and RCP8.5 scenarios from 2010 to 2100 were used in this study. Twenty-nine different GCMs and corresponding Multi-Model Ensemble (MME) means were evaluated to assess the performance of the models in reproducing the observed drought characteristics during the historical period and for deriving future changes in drought characteristics in terms of the frequency, duration, and severity. The changes were evaluated in three future time ranges, namely, the 2010s (2010–2039), 2040s (2040–2069), and 2070s (2070–2099), and data were compared to the baseline period (1976–2005).

2. Materials and Methods

2.1. Study Area and Observation Data

The study area covers Bangladesh, which is situated in latitudes between 20°34′ and 26°38′ north and longitudes between 88°01′ and 92°41′ east in South Asia; this region is bordered by India on three sides (i.e., west, north, and northeast), and it is bordered by Myanmar to the southeast. The southern border is demarcated by the Bay of Bengal with a long coastline. Although elevations up to 105 m above sea level occur in the northern part of the country, most of the elevations are less than 10 m above sea level; heights decrease in the coastal south (Figure 1). A predominance of agricultural lands is evident, and such lands comprise three-fourths of the total geographical area followed by forests, including orchards, as shown in Figure 2. Bangladesh is one of the most vulnerable nations to the increasing effects of global climate change. Presently, this country is regularly affected by natural disasters such as floods, tornadoes, droughts, and tidal bores. It has experienced drought conditions recurrently over the past several years; on average, these events have occurred once in 2.5 years [47] and have mainly affected agricultural lands, with huge losses in food grains [48,49]. While drought is a periodic occurrence in many parts of Bangladesh, the northwestern part of the country is the most susceptible to drought due to the high variability of rainfall [50]. During 2006, the average crop production was decreased by 25–30% due to the effects of drought in the northwestern part of Bangladesh [51]. This area is comparatively dry in relation to other areas of the country, as it receives much lower rainfall [52]. Additionally, this area contains sandy soils that have a low moisture retention capacity and a high infiltration rate [53]. Therefore, drought happens in this region regularly.
The climate of Bangladesh is characterized by moderately warm temperatures, high humidity, and subtropical monsoons with wide seasonal variations in rainfall. The following four metrological seasons are generally recognized: a hot, humid pre-monsoon period (March to May); a humid, warm, and rainy monsoon period (June to September); a post-monsoon period (October to November); and a dry winter (December to February). Moreover, the cropping season in Bangladesh is categorized into Pre-Kharif (March–June), Kharif (July–October), and Rabi (November–February). The mean annual temperature within the country is about 25 °C. The temperature varies from month to month. In general, the mean temperature across the country usually ranges between 11 °C and 29 °C and between 21 °C and 34 °C during the summer months. April is the warmest month in most parts of the country [54]. The annual rainfall in the region varies from 1536 mm to 4124 mm during 1976 and 2005, and the 30-year averaged annual rainfall was 2410 mm. (Table 1).
The Bangladesh Meteorological Department (BMD) operates 35 weather stations throughout the country. However, only 25 stations have continuous rainfall records for more than 30 years from 1976 to 2005 (Figure 2). In this study, the daily time series rainfall data were used to evaluate the GCM data for the same period. When working with observation data, missing data are common. Here, there were some missing data in the dataset, but the amount was <2%. The average values of the same date from neighboring stations were used in place of the missing data.

2.2. EDI Calculation

In this study, EDI was used to characterize droughts in terms of the frequency, duration, intensity, and seasonal drought days. The estimation of drought features in this method is the use of daily precipitation. The EDI calculation process was carried out according to the following equations:
E P i = n = 1 i [ ( m = 1 n P m / n ]
D E P = E P M E
E D I = D E P / S T D E P
In Equation (1), EP is the daily cumulative effective precipitation, Pm is the precipitation of m1 days ago, n is the duration of the preceding period, and i is the duration over which precipitation is summed. For example, if i is 4, then the daily EP is P 1 + P 1 + P 2 2 + P 1 + P 2 + P 3 3 + P 1 + P 2 + P 3 + P 4 4 . In this study, i was set to 365 days.
Secondly, the mean EP (MEP) is computed for each calendar day, i.e., from 1 to 365, and then the average of the analyzed 30-year period. The MEP is computed with the results from Equation (1). For instance, the MEP of the 1st day of January is the mean of 30 values for the 1st day of January composed over 30 years.
The third step is to calculate the deviation of the EP (DEP) from the MEP (Equation (2)). The DEP indicates the deficiency (negative DEP) or surplus (positive DEP) of water resources for a particular day. If the dry period is longer than 365 days, i in Equation (1) increases by the number of dry days. For instance, if the negative DEP value continues for 2 days, i is set to be 367 and Equations (1)–(3) are again calculated.
Finally, the standardized value of DEP is calculated, where ST (DEP) denotes the standard deviation of each day’s DEP. The EDI is calculated then by using Equation (3). Originally, the EDI was proposed for monitoring the drought conditions at a daily time step. Then, it was extended for monthly drought monitoring [45,55,56]. In this study, the daily EDI was used for drought calculations.

2.3. Definitions of Drought Characteristics

The EDI value illuminates the characteristics of drought events researched in this study, as shown in Figure 3. A meteorological drought event was considered to have occurred when the EDI values were less than −1, as shown in Figure 3. The drought start date was considered to be the day in which the EDI first indicated that the value of EDI was −1.0 or below, and the end date was considered to be the day when the EDI regained the value of −1.0 and above. Drought durations are the periods between the start and the end date, as shown in Figure 3. The drought severity is the cumulative deficit below the −1.0 level for the duration of a drought event, as shown in Figure 3, and the maximum intensity of a drought is the value of the minimum EDI at each event.
In this study, the classification level of dryness follows the EDI classification proposed by Kim and Byun [57], as shown in Figure 3, and drought events were also categorized based on the drought duration as shown in Table 2.
The inverse distance weight (IDW) algorithm was used to expand the data to the entire study area spatially. Twenty-five (25) weather station location points of the BMD were considered for spatial mapping of the change in drought characteristics as well as precipitation.
In this study, the following different parameters were used for drought characterization: (a) frequency—number of drought events according to the severity and duration over the period of interest, (b) duration—average duration of all drought events over the period of interest, (c) intensity—extreme intensity (minimum value of EDI) among all of the drought events over the period of interest, and (d) seasonal drought days—number of drought days in each category in different cropping seasons over the period of interest.

2.4. CMIP5 GCM Projections

The CMIP5 projection data, including precipitation for the two RCP scenarios (RCP4.5 and RCP8.5), were acquired from the outputs of 29 CMIP5 GCM models. The GCM models used in this study are described in Table 3.
The GCM outputs vary considerably in terms of the spatial resolution and include systematic biases precluding their immediate application to the assessment of climate change impacts. Therefore, certain processes to make up for these limitations are generally required before use. Daily data were downscaled for the period of 1976–2100 and bias-corrected against the observation data from 25 weather stations by using a simple quantile mapping (SQM) method [58]. The SQM technique performs independent refinements by observation points and meteorological variables through empirical quantile mapping. In this study, the following three-step process was used: (1) extract the GCM grid data corresponding to each target station, (2) estimate the biases of the retrospective simulations, and (3) bias-correct the future projections. A single grid covering the target station was extracted for each GCM, and the biases of retrospective simulation outputs for the selected grid were estimated in comparison to observations. The differences between the observed and simulated cumulative distribution functions (CDFs) for the retrospective period were quantified and then applied to the future simulations for a given percentile (Equation (4)):
x p t = x p t + F obs 1 F p . sim x p t F r . sim 1 F p . sim x p t
where x p t and x p t denote the bias-corrected and raw future projections on day t, and F(θ) and F−1(θ) are a CDF of the daily data θ and its inverse, respectively. The subscripts p.sim, r.sim, and obs indicate the future projection, retrospective simulation, and observed daily data, respectively.
The non-parametric empirical equation temporally measures the amount of daily observation data and raw GCM data. According to Gudmundsson [59], the nonparametric methods have shown the best skill in reducing the systematic bias compared to the parametric style because these use the real distribution of the observed and simulated data, without estimating a probability distribution function.

2.5. Climate Indices for Evaluations of GCMs

The Joint Expert Team on Climate Change Detection and Indices (ETCCDI) has defined a set of 27 core climate indices mainly focusing on extreme events, which are derived from the daily climate data. In this study, six extreme rainfall-related indices were selected, as shown in Table 4. These indices were used to quantitatively evaluate the future extremes and drought-related rainfall patterns of 29 CMIP5 GCM data. Analyses were conducted for 30-year periods starting from the 2010s, 2040s, and 2070s relative to the 1976–2005 baseline period.

3. Results and Discussion

3.1. Evaluations for Retrospective Simulations of GCMs

In this study, the reproducibility of output from 29 CMIP5 GCMs for the historical period from 1976 to 2005 was evaluated through comparisons with observed data concerning ETCCDI precipitation extreme indices, as shown in Table 4. To validate the performance of the bias correction algorithms, ETCCDI precipitation extreme indices were used. The statistical characteristics of each GCM model have been compared with an observed median to determine the performance of each model before (raw) and after bias-correction, and these results are presented in Figure 4. It was found that the bias-corrected GCM results reasonably captured the observed patterns of ETCCDI precipitation extremes indices, while raw GCM outputs included significant differences from the observed indices. The CDD, CWD, and R1mm of the raw GCM outputs were overestimated by 0.64%, 142.79%, and 11.02%, while R10mm, PRCPTOT, and SDII were underestimated by 43.34%, 59.51%, and 53.19%, respectively. Definitions of the indices are given in Table 4. After bias-correction, most of the errors were removed and only 11.15% and 55.57% of the errors in CDD and CWD remained, which were related to the temporal pattern of a precipitation event.

3.2. GCM Skills in Reproducing Drought Characteristics

Drought event reproducibility from GCM simulations was evaluated through EDI values calculated from the raw and bias-corrected GCM projected rainfall, which was compared with observed EDI values calculated from the BMD data from 1976 to 2005, as shown in Table 5. In this study, the average events from 29 GCMs were considered as historical drought events over 30 years. Results were compared for different BMD stations.
Table 5 presents the number of moderate, severe, and extreme drought events that occurred in the past based on observed and modeled results. The results showed that the bias-corrected GCM data could capture the averaged frequencies of moderate and severe droughts over the country quite well, while raw model data underestimated the frequency of moderate droughts by 9% and overestimated that of severe droughts by 5% compared to the observed values. However, bias-correction did not improve the spatial variability of droughts as well as the skills in representing extreme drought events. The number of observed moderate drought events varied within a range of 17–102 over the stations, whereas raw and bias-corrected modeling results showed ranges of 46–55 and 51–58, respectively. Similar limitations in representing the spatial variability of drought events were also found in the case of severe droughts. This was likely due to the issue of the coarse resolution of the GCM itself (Table 3). Specifically, the grid spacing of the GCM configuration may not be detailed enough to present the spatial distributions of the variable (i.e., the frequency) on a station basis.
The frequency of observed extreme droughts varied within a range of 0–7 over the stations, while the raw and bias-corrected modeling results showed ranges of 2–3 and 1–2, respectively, during the historical period between 1976 and 2005. As only a few extreme drought events occurred in the past, the model was not able to provide enough indications of such an occurrence. This was likely due to the differences between the observed and GCM temporal patterns of precipitation events such as the CDD and CWD, of which errors remained even after bias-correction (Figure 3). However, bias-corrected GCM data reasonably provided information on severe and moderate droughts. These findings were similarly discussed in the previous study by Hasan et al. [35].

3.3. Changes in Precipitation

This study evaluated how the characteristics of precipitation, which could potentially induce drought events, will change in the future so that we could better illuminate the driving processes underlying changes in droughts. The spatial distributions of averaged precipitation amounts from the MME mean in the 2010s, 2040s, and 2070s were compared to the baseline period (1976–2005) under the RCP4.5 and RCP8.5 scenarios, as shown in Figure 5. The projected percentage change in precipitation was uneven over the entire country, but overall, precipitation tended to increase in the future over Bangladesh.
For the beginning of the century (the 2010s), Bangladesh is expected to face an increase of up to 18.72% and 26.23% in precipitation under RCP4.5 and RCP8.5, respectively, compared to the baseline period. In particular, the maximum increases were located on the border between the Natore, Panbna, and Kustia districts. An increase of up to 14.63% in precipitation was also projected to occur by the middle of the century (the 2040s) under RCP4.5 compared to the baseline period, but the amount was less than that from the 2010s. The highest increases were found in the Shatkhira and Comilla districts, which are located in the southern and eastern regions, respectively. However, the most significant increase of rainfall, up to 35.67%, was projected during the middle of the century (the 2040s) under RCP8.5, especially in the districts in the northern region, and the increasing trend grew more severe with time. Growth of up to 26.14% and 53.95% in precipitation was also projected to occur under RCP4.5 and RCP8.5, respectively, by the end of the century (the 2070s) over Bangladesh. The highest increase in rainfall is expected to happen over the northern region of Bangladesh, especially the Dinajpur and Rajshahi districts.
Overall, the dominant feature detected was that the widespread precipitation will be increased in the 21st century over Bangladesh, especially under the RCP8.5 scenario. This result agreed with the previous study by Fahad et al. [60]. The increasing trend of precipitation may generally imply a decrease in drought occurrence.
The future changes in climate indices (PRCPTOT, CWD, CDD, R10mm, R1mm, and SDII) calculated by using data from the 29 GCMs for RCP4.5 and RCP8.5 were further investigated, as shown in Figure 6. When comparing the median values of the box and whisker diagrams, most of the climate indices showed changes indicative of increasing severity of precipitation (i.e., more precipitation with a lower number of rainy days). According to Figure 6c, the total annual precipitation (PRCPTOT) tended to increase gradually in the three future periods under both RCP scenarios; PRCPTOT increased at the end of the century (the 2070s) by around 31.98% under the RCP8.5 scenario compared to the baseline period followed by about 17.66% under the RCP4.5 scenario in the same period. Conversely, the number of rainy days (R1mm) decreased in the middle (the 2040s) and end of the century (the 2070s) by around 0.51% and 0.61%, respectively, under RCP8.5 and by 0.49% at the beginning of the century under RCP4.5 compared to the baseline period, as shown in Figure 6d. This finding was related to the gradual increase of the simple daily rainfall intensity (SDII) and heavy rainfall events (PRCP ≥ 10 mm), as shown in Figure 6e,f.
Changes in the frequency of precipitation events were also manifest as changes in the duration of dry spells (consecutive dry days, CDD) and wet spells (consecutive wet days, CWD). The CDD will be increased in the three future periods under both scenarios, as shown in Figure 6a. The highest CDD rise occurred at the end of the century (the 2070s) by around 12.08% under the RCP8.5 scenario compared to the baseline period followed by approximately 9.57% in the middle of the century (2040s) under the same scenario. However, the lowest CDD increase occurred at the beginning of the century (the 2010s) by around 4.83% under the RCP8.5 scenario compared to the base period. The CWD is also expected to show a minimum increase in the 21st century with a range of 1.12 to 2.37% compared to the base period as shown in Figure 6b. The highest CWD increase occurred at the end of the century (2070) under the RCP8.5 scenario, while the lowest rise occurred in the middle of the century (the 2040s) under the RCP4.5 scenario. Previously, Hasan et al. [37] also projected that heavy rainfall is expected to increase in the future.
Generally, the change in extreme precipitation index patterns confirms the increasing trend of precipitation in the future. Moreover, a decrease in the number of drought events and a reduction of long-term drought duration can be expected given the related increases in total precipitation in Bangladesh. However, the CDD will be increased, which may lead to an increase in extreme drought events as well as extreme drought days in the future. However, the model used in this study cannot provide the real essence of changes in extreme drought events (Table 5).

3.4. Projections of Future Changes in Drought Characteristics

Future climatological changes in drought frequency, mean duration, and maximum intensity was investigated by using the GCM simulations in the 2010s, 2040s, and 2070s compared to the baseline period (1976–2005) under the RCP4.5 and RCP8.5 scenarios. To compute the differences between baseline and future drought characteristics, the percent change for the three next periods under each scenario were calculated by counting the number of events for each level of drought severity and duration for the three future periods under both scenarios.

3.4.1. Changes in the Drought Frequency

Figure 7 compares the drought frequency for each level of drought severity (i.e., extreme, severe, and moderate drought events) and future periods under the RCP4.5 and RCP8.5 scenarios. Overall, the projected changes in drought frequency for the three future periods showed different characteristics. However, the respective change in climatological severity under RCP4.5 and RCP8.5 exhibited a constant magnitude response for any future period, thus indicating an insignificant effect of emission scenarios on the projected change in extreme and severe droughts. Besides, extreme droughts are not expected to change much in the 21st century under both scenarios.
According to Figure 7, it was found that the MME mean produced the maximum number of moderate drought events compared to the others. For the beginning of the century (the 2010s), moderate droughts decreased by around 9% and 14% under RCP4.5 and RCP8.5, respectively. However, only severe droughts increased by about 15% under both of the scenarios. Moreover, moderate droughts decreased by about 18% in the middle of the century (the 2040s) compared to the base period under both of the scenarios, while severe droughts are expected to remain the same under both scenarios. Furthermore, the moderate and severe droughts events are expected to decrease by around 29% and 12%, respectively, at the end of the century (the 2070s) compared to the base period under the RCP4.5 scenario, while 38% of moderate and 23% of severe drought events will decrease in the same period under the RCP8.5 scenario.
The spatial distributions of climatological changes in drought frequency for the future periods are shown in Figure 8. The projected percentage change in drought frequency was non-uniform over the country. However, the overall drought frequency will likely decrease with time under climate change. The CMIP5 does not well capture the spatial variability in Bangladesh, as discussed in Table 5. Therefore, the spatial distributions analyzed here may contain large uncertainty.
For the beginning of the century (the 2010s), the drought frequency appeared to decrease in most of the country. In particular, the maximum decrease was detected in the eastern mountainous regions, especially the Chittagong and Bandarban districts, with ranges up to 10–15% compared to the base period under the RCP4.5 scenario. However, it increased up to 3% in the districts of Jessore and Satkhira, but only under RCP4.5, by the 2010s compared to the base period. On the other hand, the drought frequency decreased all over the country under the RCP8.5 scenario. The highest decrease was detected in the districts of Noakhali, Hobigonj, and Moulovibazar during the 2010s. Furthermore, a decrease of up to 23% in drought frequency was also projected to occur by the middle of the century (the 2040s) under RCP4.5 compared to the baseline period, which was greater than that from the 2010s. The highest decrease was detected in the districts of Mymensingh and Chittagong. A decrease in frequency, up to 24%, was projected to occur over Bangladesh during the middle of the century (the 2040s) under RCP8.5, especially in the eastern region. A greater decrease in frequency can be expected over Bangladesh for the 2070s than for the 2010s and 2040s. A decrease of up to 33.33% and 46.77% in drought frequency was also projected under RCP4.5 and RCP8.5, respectively, to occur by the end of the century (the 2070s) over Bangladesh relative to the baseline period of 1976–2005. The highest decrease in frequency was projected to occur over the drought-vulnerable northern and northeastern regions of Bangladesh, especially in the Rajshahi and Sylhet districts.

3.4.2. Changes in the Drought Duration

Figure 9 presents the changes in drought duration for the 21st century based on the historical period (1976–2005) concerning long-term, medium-term, short-term, and very short-term drought events under the RCP4.5 and RCP8.5 scenarios. Overall, the projected changes in drought duration for the three future periods (the 2010s, 2040s, and 2070s) showed different characteristics. However, the respective change in the climatological length of droughts under RCP4.5 and RCP8.5 exhibited a stable magnitude for any future period, thus showing a negligible effect of emission scenarios. The medium-term droughts are not expected to change in the 21st century under both scenarios except by the end of the century (the 2070s) under RCP8.5, where a decrease of 55% was detected compared to the base period. Only the long-term droughts will decrease by around 55% over the three future periods under both scenarios compared to the base period.
According to Figure 9, The short-term and very short-term droughts will be decreased in the 21st century under both scenarios except for at the beginning of the century (the 2010s) under RCP4.5. The highest decreasing rates of very short-term and short-term drought events can be expected at the end of the century and will amount to around 26% and 21%, respectively, compared to the base period under the RCP4.5 scenario, while 36% of the very short-term and 31% of the short-term drought events decreased in the same period under the RCP8.5 scenario.
The mean duration of droughts was computed from the average duration of all drought events in the historical and three future periods with the MME mean of the GCMs. The climatological changes in mean drought duration for the future periods are shown in Figure 10. The difference is a percentage that designates the change in the mean duration of droughts, where a negative value indicates a decreasing trend and vice versa. The average length of droughts can be expected to decline in the 21st century all over the country. However, the decreasing rate varied with the projected period, location, and emission scenario. The results indicate that the decreasing rate will be increased with time and the highest decreasing rate can be expected during the end of the century (the 2070s) under RCP4.5.
In particular, the mean duration of droughts is expected to decrease within a range of 0 to 15% over the entire country, whereas the highest decrease will occur in the districts of Chittagong, Cox’s Bazar, and Jessore at the beginning of the century (the 2010s). The values decreased within a range of 6–24% and 14–36% in the middle of the century (the 2040s) and end of the century (2070s), respectively, under the RCP4.5 scenario. The highest decreases occurred in the hilly district of Rangamati in the 2040s and in the Rangpur and Dinajpur districts in the 2070s. On the other hand, the RCP8.5 scenario results tended to show a higher decreasing rate in the northern region, especially in the Bogra district. In particular, the mean duration of droughts decreased within a range of 0 to −16% all over the country at the beginning of the century (the 2010s), and decreases occurred within a range of 6–24% and 10–28% in the middle of the century (the 2040s) and end of the century (2070s), respectively, under the RCP8.5 scenario.

3.4.3. Changes in the Drought Intensity

The intensity of the EDI was computed from the MME means derived from the GCMs, and resulting data were used to assess the extreme intensity trends (minimum values of the EDI) among all of the drought events in the three future periods compared to the historical period (1976–2005). The difference is a percentage designating the change in maximum intensity of droughts, where a positive value indicates an increasing trend and vice versa. The climatological changes in the maximum drought intensity for future periods are shown in Figure 11. Overall, the projected change in maximum drought intensity for the three future periods showed different characteristics. The maximum drought intensity increased over most of the country at the beginning of the century (the 2010s). However, a greater decrease can be expected for the end of the century (the 2070s) under RCP8.5.
In particular, the maximum increase occurred in northern and southwestern coastal regions, especially the Borguna district, with ranges up to 8–10% under the RCP4.5 scenario, whereas a maximum increase up to 6–8% was also projected under the RCP8.5 scenario, in locations including the southwestern and central regions, during the beginning of the century compared to the base period. However, values decreased up to 1–2% in the districts of Rajshahi and Sylhet only under RCP4.5 in the 2010s, and negligible decreases were detected under the RCP8.5 scenario.
Besides, the maximum drought intensity for the middle of the century showed a tendency to decrease under RCP4.5 but an increasing trend under RCP8.5. Precisely, a decrease of up to 6–8% in the drought intensity is projected to occur by the middle of the century (the 2040s) under RCP4.5 compared to the baseline period. The higher decrease was located in the hilly eastern region and southwestern coastal region, especially in the districts of Rangamati and Chittagong. Conversely, an increase in the intensity of up to 6–8% was projected to occur over Bangladesh during the middle of the century (the 2040s) under RCP8.5, especially in the central region.
A decrease of up to 4–6% and 8–10% in the maximum drought intensity was also projected under RCP4.5 and RCP8.5, respectively, by the end of the century (the 2070s) over Bangladesh relative to the baseline period of 1976–2005. The highest decrease in intensity was detected over the drought-vulnerable northern region and northeastern region of Bangladesh, especially in the Rajshahi and Sylhet districts under the RCP8.5 scenario.

3.5. Projections of Seasonal Changes in the Number of Drought Days

In this study, the changes in moderate, severe, and extreme drought days were investigated during the Pre-Kharif, Kharif, and Rabi seasons for the 2010s, 2040s, and 2070s compared to the baseline period (1976–2005) under the RCP4.5 and RCP8.5 scenarios. The different categories of drought days were counted by using EDI values based on the MME means of 29 GCMs under both scenarios. The drought days estimated at different stations are presented by using box plots. The 25th and 75th quartiles of drought days are demarked by the lower and upper lines of the boxes, whereas the median value drought days are indicated by the middle lines of the boxes at different stations over Bangladesh. The changes in drought days for the Pre-Kharif, Kharif, and Rabi drought seasons are shown in Figure 12. Overall, an increase of median values of extreme drought days was observed with some exceptions, whereas the opposite trends in moderate and severe drought days were detected under both scenarios compared to the base period.
Particularly, the maximum decrease in the moderate and severe drought days was found in the Rabi season with changes of around 50% and 44% under RCP4.5 and 44% and 53% under RCP8.5, respectively, during the end of the century (the 2070s) compared to the base period. However, the severe droughts days increased at the beginning of the century in the Pre-Kharif and Kharif seasons by around 28% and 3%, respectively, under RCP4.5, and 12% in the Pre-Kharif season under RCP8.5. The moderate drought days are expected to increase at the beginning of the century in the Pre-Kharif and Kharif seasons by 14% and 7% under RCP4.5, respectively, and 6% in the middle of the century under RCP8.5, respectively. However, values are not expected to increase in the future under the RCP4.5 scenario. Only 30% of the extreme drought days will be decreased in the Kharif season in the middle of the century under RCP4.5, while 17% and 10% decreases are expected in the Kharif and Rabi seasons, respectively, at the end of the century under the RCP8.5 scenario. The highest increase in extreme drought days was found in the Pre-Kharif season, and the increase was around 205% at the end of the century under RCP4.5.

4. Conclusions

Climate change is regarded as a future threat that is expected to change regional hydrological patterns as well as the characteristics of droughts, which could lead to severe drought-related disasters. Hence, it is important to assess the impacts of climate change on potential future droughts. Such work will allow us to respond preventively by providing policy analysis information and useful guidance to policymakers, who can devise a framework for water resource management as well as agricultural adaptation in the future. In this regard, the principal objectives of this study were to evaluate the changes in drought-related precipitation extremes and drought characteristics over Bangladesh by using 29 GCM projections under the RCP4.5 and RCP8.5 scenarios for the period of 2010–2099. The Effective Drought Index (EDI) was used for the characterization of droughts in terms of the frequency, duration, intensity, and number of drought days. To evaluate the reproducibility of the GCMs, retrospective simulations of models for the historical period from 1976 to 2005 were compared to the quantities based on observed climatology. It was found that the bias-corrected GCM results showed an appropriate consistency with observation data and were successfully able to reproduce the real essence of the drought situation compared to the raw GCM data except for extreme drought events. Future changes of drought characteristics were then investigated for the entire area by using the historical period as the calibration period. The differences between climate change scenarios were not recognizable from the results because a statistical trend analysis was not included in the study. The changes in drought characteristics in the future were examined only through comparisons with the historical period. The increasing or decreasing trends of future precipitation and drought characteristics were visually inspected in this study.
The future changes in precipitation tended to show increases under both scenarios; the increasing rate was higher under RCP8.5 than RCP4.5. The projected change in climatological frequencies of drought based on severity and duration under RCP4.5 and RCP8.5 exhibited a consistent magnitude for any future period, thus indicating a negligible effect of emission scenarios on the projected changes in the extreme and severe droughts and long-term and medium-term droughts. Overall, the drought frequency, mean duration, and maximum intensity will likely decrease with time under climate change in relation to the increase of total precipitation. In particular, the occurrence of severe and moderate droughts will be less frequent in the 21st century, according to this research. The frequency of extreme drought events did not show a significant change in the 21st century under both scenarios. However, the extreme drought days are likely to be increased in most of the cropping season and future periods under both scenarios, which may affect agricultural production in the future, as driven by the increasing pattern of CDD. The spatial pattern of change in drought characteristics indicates that the drought-vulnerable areas will be shifted from the northwestern to central and southern coastal regions in the future due to the effects of climate change. However, GCM may predict the trends well but may not reveal the change in a spatial pattern in detail.
Lastly, to our best of knowledge, this study is the first attempt to characterize the future droughts in Bangladesh widely using the EDI. The present results can help resource managers to optimally allocate scarce water resources and develop long-term strategies for protecting communities against natural hazards related to water scarcity. Furthermore, the outcomes of the study are expected to represent important measures for mitigating the losses in agricultural production for drought-prone areas in Bangladesh.

Author Contributions

M.K. and M.-W.J. are equal contributions as a first author; M.K. and S.H. designed the research, analyzed the data, and wrote the manuscript; M.-W.J. helped in the preparation of the manuscript and subsequent revisions; J.C. provided important intellectual content; S.H. supervised the study and providing critical evaluations of the manuscript.

Funding

This research was funded by Research Program for Agricultural Science & Technology Development, National Institute of Agricultural Sciences, Rural Development Administration, Korea, project number PJ01254903.

Acknowledgments

The authors are grateful to the Bangladesh Meteorological Department (BMD) for providing the rainfall data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Basak, J.K.; Titumir, R.A.M.; Dey, N.C. Climate Change in Bangladesh: A Historical Analysis of Temperature and Rainfall Data. J. Environ. 2013, 2, 41–46. [Google Scholar]
  2. Ahmad, Q.K.; Warrick, R.A.; Ericksen, N.J.; Mirza, M.Q. The Implications of Climate and Sea–Level Change for Bangladesh; Warrick, R.A., Ahmad, Q.K., Eds.; Kluwer Academic Publisher: Dordrecht, The Netherlands, 1996; ISBN 978-94-009-0241-1. [Google Scholar]
  3. Shahid, S. Recent trends in the climate of Bangladesh. Clim. Res. 2010, 42, 185–193. [Google Scholar] [CrossRef] [Green Version]
  4. Rahman, M.R.; Lateh, H. Climate change in Bangladesh: A spatiotemporal analysis and simulation of recent temperature and rainfall data using GIS and time series analysis model. Appl. Clim. 2017, 128, 27–41. [Google Scholar] [CrossRef]
  5. Dai, A. Drought under global warming: A review. Wiley Interdiscip. Rev. Clim. Chang. 2011, 2, 45–65. [Google Scholar] [CrossRef]
  6. Brown, C.; Meeks, R.; Hunu, K.; Yu, W. Hydro climatic risk to economic growth in sub-Saharan Africa. Clim. Chang. 2011, 106, 621–647. [Google Scholar] [CrossRef]
  7. Hartmann, D.L.; Tank, A.M.G.K.; Rusticucci, M. Climatic Change 2013: The Physical Science Basis; IPCC Fifth Assessment Report; IPCC: Geneva, Switzerland, 2013; pp. 31–39. [Google Scholar]
  8. Inter-Governmental Panel on Climate Change (IPCC). Summary for Policymakers. In Climate Change 2013: The Physical Science Basis; Contribution of Working Group I to the Fifth Assessment Report of IPCC; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. [Google Scholar]
  9. Solomon, S. (Ed.) Climate Change—The Physical Science Basis: Working Group I Contribution to the Fourth Assessment Report of the IPCC 4; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2007. [Google Scholar]
  10. Van Vuuren, D.P.; Edmonds, J.; Kainuma, M. The representative concentration pathways: An overview. Clim. Chang. 2011, 109, 5–31. [Google Scholar] [CrossRef]
  11. Taylor, K.E.; Stouffer, R.J.; Meehl, G.A. An overview of CMIP5 and the experiment design. Bull. Am. Meteorol. Soc. 2012, 93, 485–498. [Google Scholar] [CrossRef]
  12. Sperber, K.R.; Annamalai, H.; Kang, I.S.; Kitoh, A.; Moise, A.; Turner, A.; Zhou, T. The Asian summer monsoon: An intercomparison of CMIP5 vs. CMIP3 simulations of the late 20th century. Clim. Dyn. 2013, 41, 2711–2744. [Google Scholar] [CrossRef]
  13. Bhaskaran, B.; Ramachandran, A.; Jones, R.; Moufouma-Okia, W. Regional climate model applications on sub-regional scales over the Indian monsoon region: The role of domain size on downscaling uncertainty. J. Geophys. Res. Atmos. 2012, 117. [Google Scholar] [CrossRef]
  14. Xu, C.Y.; Widen, E.; Halldin, S. Modelling hydrological consequences of climate change-progress and challenges. Adv. Atmos. Sci. 2005, 22, 789–797. [Google Scholar] [CrossRef]
  15. Fowler, H.J.; Kilsby, C.G. Using regional climate model data to simulate historical and future river flows in northwest England. Clim. Chang. 2007, 80, 337–367. [Google Scholar] [CrossRef]
  16. Hong, S.Y.; Kanamitsu, M. Dynamical downscaling: Fundamental issues from an NWP point of view and recommendations. Asia Pac. J. Atmos. Sci. 2014, 50, 83–104. [Google Scholar] [CrossRef]
  17. Kang, S.; Hur, J.; Ahn, J.B. Statistical downscaling methods based on APCC multi-model ensemble for seasonal prediction over South Korea. Int. J. Clim. 2014, 34, 3801–3810. [Google Scholar] [CrossRef]
  18. Lee, J.W.; Hong, S.Y. Potential for added value to downscaled climate extremes over Korea by increased resolution of a regional climate model. Theor. Appl. Climatol. 2014, 117, 667–677. [Google Scholar] [CrossRef]
  19. Kunkel, K.E.; Liang, X.Z.; Zhu, J.; Lin, Y. Can CGCMs simulate the twentieth-century “warming hole” in the central United States? J. Clim. 2006, 19, 4137–4153. [Google Scholar] [CrossRef]
  20. Fowler, H.J.; Blenkinsop, S.; Tebaldi, C. Linking climate change modeling to impacts studies: Recent advances in downscaling techniques for hydrological modeling. Int. J. Clim. 2007, 27, 1547–1578. [Google Scholar] [CrossRef]
  21. Hay, L.E.; Wilby, R.L.; Leavesley, G.H. A comparison of delta change and downscaled GCM scenarios for three mountainous basins in the United States. J. Am. Water Resour. Assoc. 2000, 36, 387–397. [Google Scholar] [CrossRef]
  22. Wood, A.W.; Maurer, E.P.; Kumar, A.; Lettenmaier, D.P. Long-range experimental hydrologic forecasting for the eastern United States. J. Geophys. Res. 2002, 107, 4429. [Google Scholar] [CrossRef]
  23. Heo, J.H.; Ahn, H.; Shin, J.Y.; Kjeldsen, T.R.; Jeong, C. Probability Distributions for a Quantile Mapping Technique for a Bias Correction of Precipitation Data: A Case Study to Precipitation Data Under Climate Change. Water 2019, 11, 1475. [Google Scholar] [CrossRef]
  24. Zhao, T.; Dai, A. Uncertainties in historical changes and future projections of drought. Part II: Model-simulated historical and future drought changes. Clim. Chang. 2016, 144, 535–548. [Google Scholar] [CrossRef]
  25. Cohelo, C.A.S.; Goddard, L. El Niño-Induced Tropical Droughts in Climate Change Projections. J. Clim. 2009, 22, 6456–6476. [Google Scholar] [CrossRef]
  26. Rhee, J.; Cho, J. Future Changes in Drought Characteristics: Regional Analysis for South Korea under CMIP5 Projections. J. Hydrometeorol. 2016, 17, 437–451. [Google Scholar] [CrossRef]
  27. Touma, D.; Ashfaq, M.; Nayak, M.A.; Kao, S.C.; Diffenbaugh, N.S. A multi-model and multi-index evaluation of drought characteristics in the 21st century. J. Hydrol. 2015, 526, 196–207. [Google Scholar] [CrossRef] [Green Version]
  28. Dai, A. Increasing drought under global warming in observations and models. Nat. Clim. Chang. 2013, 3, 52–58. [Google Scholar] [CrossRef]
  29. Orlowsky, B.; Seneviratne, S.I. Elusive drought: Uncertainty in observed trends and short- and long-term CMIP5 projections. Hydrol. Earth Syst. Sci. 2013, 17, 1765–1781. [Google Scholar] [CrossRef]
  30. Chen, H.P.; Sun, J.Q.; Chen, X.L. Future changes of drought and flood events in China under a global warming scenario. Atmos. Ocean. Sci. Lett. 2013, 6, 8–13. [Google Scholar] [CrossRef]
  31. Wang, L.; Chen, W.; Zhou, W. Assessment of future drought in Southwest China based on CMIP5 multimodal projections. Adv. Atmos. Sci. 2014, 31, 1035–1050. [Google Scholar] [CrossRef]
  32. Burke, E.J.; Brown, S.J. Evaluating uncertainties in the projection of future drought. J. Hydrometeorol. 2008, 9, 292–299. [Google Scholar] [CrossRef]
  33. Nowreen, S.; Murshed, S.B.; Islam, A.K.M.S. Changes of rainfall extremes around the haor basin areas of Bangladesh using multi-member ensemble RCM. Appl. Clim. 2015, 119, 363–377. [Google Scholar] [CrossRef]
  34. Rahman, M.M.; Islam, M.N.; Ahmed, A.U.; Georgi, F. Rainfall and temperature scenarios for Bangladesh for the middle of 21st century using RegCM. J. Earth. Syst. Sci. 2012, 121, 287–295. [Google Scholar] [CrossRef] [Green Version]
  35. Hasan, M.A.; Saiful Islam, A.K.M.; Bokhtiar, S.M. Future change of the metrological drought over Bangladesh using high-resolution climate scenarios. In Proceedings of the International Conference on Climate Change Impact and Adaptation (I3CIA-2013), Gazipur, Bangladesh, 15–17 November 2013. [Google Scholar]
  36. Islam, A.R.M.T.; Shen, S.; Hu, Z.; Rahman, M.A. Drought Hazard Evaluation in Boro Paddy Cultivated Areas of Western Bangladesh at Current and Future Climate Change Conditions. Adv. Meteorol. 2017, 2017, 3514381. [Google Scholar] [CrossRef]
  37. Hasan, M.A.; Saiful Islam, A.K.M.; Akanda, A.S. Climate projections, and extremes in dynamically downscaled CMIP5 model output over the Bengal delta: A quartile based bias-correction approach with new gridded data. Clim. Dyn. 2018, 51, 2169–2190. [Google Scholar] [CrossRef]
  38. Mortuza, M.R.; Moges, E.; Demissie, Y.; Hong-Yi, L. Historical and future drought in Bangladesh using copula-based bivariate regional frequency analysis. Appl. Clim. 2019, 135, 855–871. [Google Scholar] [CrossRef]
  39. McKee, T.B.; Doesken, N.J.; Kleist, J. The relationship of drought frequency and duration to time scales. Preprints. In Proceedings of the Eighth Conference on Applied Climatology, Anaheim, CA, USA, 17–22 January 1993; pp. 179–184. [Google Scholar]
  40. Byun, H.R.; Kim, D.W. Comparing the Effective Drought Index and the Standardized Precipitation Index. In: López-Francos, A. (comp.), López-Francos, A. (collab.). Zaragoza: CIHEAM/FAO/ICARDA/GDAR/CEIGRAM/MARM. Econ. Drought Drought Prep. A Clim. Chang. Context 2010, 95, 85–89. [Google Scholar]
  41. Byun, H.R.; Wilhite, D.A. Objective quantification of drought severity and duration. J. Clim. 1999, 12, 2747–2756. [Google Scholar] [CrossRef]
  42. Morid, S.; Smakhtin, V.; Moghaddasi, M. Comparison of Seven Meteorological drought indices for drought monitoring in Iran. Int. J. Climatol. 2006, 26, 971–985. [Google Scholar] [CrossRef]
  43. Pandey, R.P.; Dash, B.B.; Mishra, S.K.; Singh, R. Study of indices for drought characterization in KBK districts in Orissa (India). Hydrol. Process. 2008, 22, 1895–1907. [Google Scholar] [CrossRef]
  44. Dogan, S.; Berktay, A.; Singh, V.P. Comparison of multi-monthly rainfall-based drought severity indices, with application to semi-arid Kenya, closed basin, Turkey. J. Hydrol. 2012, 470, 255–268. [Google Scholar] [CrossRef]
  45. Jain, V.K.; Pandey, R.P.; Jain, M.K.; Byun, H.R. Comparison of drought indices for appraisal of drought characteristics in the Ken River Basin. Weather Clim. Extrem. 2015, 8, 1–11. [Google Scholar] [CrossRef] [Green Version]
  46. Kamruzzaman, M.; Cho, J.; Jang, M.; Hwang, S. Comparative Evaluation of Standardized Precipitation Index (SPI) and Effective Drought Index (EDI) for Meteorological Drought Detection over Bangladesh. J. Korean Soc. Agric. Eng. 2019, 61, 143–157. [Google Scholar] [CrossRef]
  47. Mondal, M.A.H.; Ara, I.; Das, S.C. Meteorological Drought Index Mapping in Bangladesh Using Standardized Precipitation Index during 1981–2010. Adv. Meteorol. 2017, 2017, 4642060. [Google Scholar] [CrossRef]
  48. Murad, H.; Islam, A.S. Drought assessment using remote sensing and GIS in north-west region of Bangladesh. In Proceedings of the 3rd International Conference on Water & Flood Management, Dhaka, Bangladesh, 8–10 January 2011; pp. 797–804. [Google Scholar]
  49. Rafiuddin, M.; Dash, B.K.; Khanam, F. Diagnosis of Drought in Bangladesh using Standardized Precipitation Index. In Proceedings of the International Conference on Environment Science and Engineering; IACSIT Press: Singapore, 2011. [Google Scholar]
  50. Shahid, S.; Behrawan, H. Drought risk assessment in the western part of Bangladesh. Nat. Hazards 2008, 46, 391–413. [Google Scholar] [CrossRef]
  51. Rahman, A.A.; Alam, M.; Alam, S.S.; Uzzaman, M.R.; Rashid, M.; Rabbani, G. Risks, Vulnerability, and Adaptation in Bangladesh; UNDP Human Development Report, Bangladesh Centre for Advanced Studies (BCAS): Dhaka, Bangladesh, 2008. [Google Scholar]
  52. Paul, B.K. Coping mechanisms practiced by drought victims (1994/95) in North Bengal, Bangladesh. Appl. Geogr. 1998, 18, 355–373. [Google Scholar] [CrossRef]
  53. Habiba, U.; Hassan, A.W.R.; Shaw, R. Livelihood adaptation in the drought-prone areas of Bangladesh. In Climate Change Adaptation Actions in Bangladesh; Springer: Tokyo, Japan, 2013; pp. 227–252. ISBN 978-4-431-54248-3. [Google Scholar]
  54. Rashid, H.E. Geography of Bangladesh; University Press Limited: Dhaka, Bangladesh, 1991. [Google Scholar]
  55. Deo, R.C.; Sahin, M. Application of the extreme learning machine algorithm for the prediction of monthly Effective Drought Index in eastern Australia. Atmos. Res. 2015, 153, 512–525. [Google Scholar] [CrossRef] [Green Version]
  56. Morid, S.; Smakhtin, V.; Bagherzadeh, K. Drought forecasting using artificial neural networks and time series of drought indices. Int. J. Climatol. 2007, 27, 2103–2111. [Google Scholar] [CrossRef]
  57. Kim, D.W.; Byun, H.R. Future pattern of Asian drought under global warming scenario. Appl. Clim. 2009, 98, 137–150. [Google Scholar] [CrossRef] [Green Version]
  58. Cho, J.; Cho, W.; Jung, I. RSQM: Statistical Downscaling Toolkit for Climate Change Scenario Using Nonparametric Quantile Mapping. 2018. Available online: cran.r-project.org/web/packages/rSQM/index.html (accessed on 24 May 2018).
  59. Gudmundsson, L.; Bremnes, J.B.; Haugen, J.E.; Engen-Skaugen, T. Technical Note: Downscaling RCM precipitation to the station scale using statistical transformations—A comparison of methods. Hydrol. Earth Syst. Sci. 2012, 16, 3383–3390. [Google Scholar] [CrossRef]
  60. Fahad, M.G.R.; Saiful Islam, A.K.M.; Nazari, R.; Hasan, M.A.; Tarekul Islam, G.M.; Bala, S.K. Regional changes of precipitation and temperature over Bangladesh using bias-corrected multi-model ensemble projections considering high-emission pathways. Int. J. Climatol. 2007, 38, 1634–1648. [Google Scholar] [CrossRef]
Figure 1. The elevation map of Bangladesh (data source: https://SRTM.csi.cgiar.org).
Figure 1. The elevation map of Bangladesh (data source: https://SRTM.csi.cgiar.org).
Water 11 02219 g001
Figure 2. Location of meteorological stations and land use map of the study area.
Figure 2. Location of meteorological stations and land use map of the study area.
Water 11 02219 g002
Figure 3. Schematic concept of the drought characteristics (severity, duration) evaluated with the EDI time series.
Figure 3. Schematic concept of the drought characteristics (severity, duration) evaluated with the EDI time series.
Water 11 02219 g003
Figure 4. Comparison of raw (left) and bias-corrected (right) reproducibility with observed data for all meteorological stations in terms of the (a) CDD, (b) CWD, (c) PRCPTOT, (d) R1mm, (e) SDII, and (f) R10mm.
Figure 4. Comparison of raw (left) and bias-corrected (right) reproducibility with observed data for all meteorological stations in terms of the (a) CDD, (b) CWD, (c) PRCPTOT, (d) R1mm, (e) SDII, and (f) R10mm.
Water 11 02219 g004
Figure 5. Percentage (%) change of precipitation over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and (2) RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Figure 5. Percentage (%) change of precipitation over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and (2) RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Water 11 02219 g005
Figure 6. Future changes of climate indices under RCP4.5 (left) and RCP8.5 (right); (a) CDD, (b) CWD, (c) PRCPTOT, (d) R1mm, (e) SDII, and (f) R10mm.
Figure 6. Future changes of climate indices under RCP4.5 (left) and RCP8.5 (right); (a) CDD, (b) CWD, (c) PRCPTOT, (d) R1mm, (e) SDII, and (f) R10mm.
Water 11 02219 g006
Figure 7. Change in the drought intensity based on a historical period (1976–2005) in regard to extreme, severe, and moderate drought events under the (a) RCP4.5 (left) and (b) RCP8.5 (right) scenarios.
Figure 7. Change in the drought intensity based on a historical period (1976–2005) in regard to extreme, severe, and moderate drought events under the (a) RCP4.5 (left) and (b) RCP8.5 (right) scenarios.
Water 11 02219 g007
Figure 8. Percentage (%) change in the drought frequency over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Figure 8. Percentage (%) change in the drought frequency over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Water 11 02219 g008
Figure 9. Change in the drought duration based on a historical period (1976–2005) in regard to long-term, medium-term, short-term, and very short-term drought events under the (a) RCP4.5 (left) and (b) RCP8.5 (right) scenarios.
Figure 9. Change in the drought duration based on a historical period (1976–2005) in regard to long-term, medium-term, short-term, and very short-term drought events under the (a) RCP4.5 (left) and (b) RCP8.5 (right) scenarios.
Water 11 02219 g009
Figure 10. Percentage (%) change in the mean drought duration over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Figure 10. Percentage (%) change in the mean drought duration over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Water 11 02219 g010
Figure 11. Percentage (%) change in the maximum drought intensity over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Figure 11. Percentage (%) change in the maximum drought intensity over Bangladesh for the 2010s (1st column), 2040s (2nd column), and 2070s (3rd column) under the RCP4.5 (upper row) and RCP8.5 (lower row) scenarios relative to the baseline period (1976–2005).
Water 11 02219 g011
Figure 12. Future changes in seasonal drought days under RCP4.5 (left) and RCP8.5 (right); (a) moderate, (b) severe, and (c) extreme droughts.
Figure 12. Future changes in seasonal drought days under RCP4.5 (left) and RCP8.5 (right); (a) moderate, (b) severe, and (c) extreme droughts.
Water 11 02219 g012
Table 1. Mean monthly and annual rainfall at various stations in Bangladesh.
Table 1. Mean monthly and annual rainfall at various stations in Bangladesh.
Station NameData PeriodLocationMean Monthly Rainfall of Multiple Stations During 1976–2005 (mm)Annual Avg. Precip. (mm)
Lon.Lat.JanFebMarAprMayJunJulAugSepOctNovDec
Rangpur1954 to present89.2325.7381125117274453498356383176892318
Dinajpur1948 to present 88.6825.65101114762223584713733701609102082
Bogra1948 to present89.3724.85813217920932040229330915212111828
Rajshahi1964 to present88.7024.371016256514026132725528512714111536
Ishwardi1961 to present89.0524.1362133861942923262302881171691616
Jessore1948 to present89.1723.181426467717631831627426913829131695
Khulna1948 to present89.5322.78123953841933553043292511313781795
Shatkhira1948 to present89.0822.72133942901592953352962801273391717
Barishal1949 to present90.3722.751025551192124184193622821784682132
Patuakhali1973 to present90.3322.33822431152385355404503421855362536
Khepupara1974 to present90.2321.9892449972595136064903692435582721
Bhola1966 to present90.6522.681030531312654714473882911724372308
Maizdi Court1951 to present91.1022.871027761543355747506313841814473172
Swandip1966 to present91.4322.48922681463496998606214362565093526
Dhaka1953 to present90.3823.777226914631834635929832618329122115
Mymensingh1948 to present90.4324.727223814535939545332632221618102309
Hatiya1966 to present91.1022.434144211623754155748432219336122559
Chandpur1964 to present90.7023.27620611392473413753262591383971956
Comilla1964 to present91.1823.437247014932235941131825015434102108
Feni1973 to present91.4223.03628691803695356524963301814592900
Sylhet1956 to present91.8824.973414936757176983360252922228134124
Srimangal1948 to present91.7324.35308922144544237133629616332152445
Chittagong1949 to present91.8122.354235012829256064548622717960132669
Rangamati1957 to present92.2022.535236213933350457544229415457132601
Cox’es Bazar1948 to present91.9721.455183211730181286966835719895143486
Mean8235313128145950840532217337102410
± ±±±±±±±±±±±±±
STDV382761961481761266636203649
Table 2. Drought classification based on the drought duration.
Table 2. Drought classification based on the drought duration.
Duration (Days)Category
Less than or equal to 30Very short-term
31 to 90Short-term
91 to 180Medium-term
Greater than 180Long-term
Table 3. Summary of 29 climate models from CMIP5 used in this study.
Table 3. Summary of 29 climate models from CMIP5 used in this study.
Model NameModeling CenterResolution (Lon × Lat)
bcc-csm1-1Beijing Climate Center, China Meteorological Administration, China2.81° × 2.79°
bcc-csm1-1-m1.13° × 1.12°
CanESM2Canadian Centre for Climate Modelling and Analysis, Canada2.81° × 2.79°
CCSM4National Center for Atmospheric Research, USA1.25° × 0.94°
CESM1-BGCNational Science Foundation, Department of Energy, National Center for Atmospheric Research, USA1.25° × 0.94°
CESM1-CAM5
CMCC-CMCentro Euro-Mediterraneo sui Cambiamenti Climatici, Italy0.75° × 0.75°
CMCC-CMS1.88° × 1.86°
CNRM-CM5Centre National de Recherches Meteorologiques/Centre Europeen de Recherche et Formation Avancees en Calcul Scientifique, France1.41° × 1.40°
CSIRO-Mk3-6-0Commonwealth Scientific and Industrial Research Organisation in collaboration with the Queensland Climate Change Centre of Excellence, Australia1.88° × 1.86°
FGOALS-g2Institute of Atmospheric Physics, Chinese Academy of Sciences; and CESS, Tsinghua University, China2.81° × 3.05°
FGOALS-s2Institute of Atmospheric Physics, Chinese Academy of Sciences, China2.81° × 1.66°
GFDL-CM3Geophysical Fluid Dynamics Laboratory, USA2.50° × 2.00°
GFDL-ESM2G
GFDL-ESM2M
HadGEM2-AONational Institute of Meteorological Research/Korea Meteorological Administration, South Korea1.88° × 1.25°
HadGEM2-CCMet Office Hadley Centre (additional HadGEM2-ES realizations contributed by Instituto Nacional de Pesquisas Espaciais), UK1.88° × 1.25°
HadGEM2-ES
inmcm4Institute of Numerical Mathematics, Russia2° × 1. 5°
IPSL-CM5A-LRInstitut Pierre-Simon Laplace, France3.75° × 1.89°
IPSL-CM5A-MR2.50° × 1.27°
IPSL-CM5B-LR3.75° × 1.89°
MIROC5Atmosphere and Ocean Research Institute (The University of Tokyo), National Institute for Environmental Studies, and Japan Agency for Marine-Earth Science and Technology, Japan1.41°×1.40°
MIROC-ESMJapan Agency for Marine-Earth Science and Technology, Atmosphere and Ocean Research Institute (The University of Tokyo), and National Institute for Environmental Studies, Japan2.81° × 2.79°
MIROC-ESM-CHEM
MPI-ESM-LRMax Planck Institute for Meteorology (MPI-M), Germany1.88° × 1.86°
MPI-ESM-MR
MRI-CGCM3Meteorological Research Institute, Japan1.13° × 1.12°
NorESM1-MNorwegian Climate Centre, Norway2.50° × 1.89°
Table 4. List of the ETCCDI extreme climate indices.
Table 4. List of the ETCCDI extreme climate indices.
IDIndicator NameDefinitionUnit
PRCPTOTAnnual total wet-day precipitationAnnual total PRCP in wet days (RR ≥ 1 mm)mm
CDDConsecutive dry daysMaximum number of consecutive days with RR < 1 mmdays
CWDConsecutive wet daysMaximum number of consecutive days with RR ≥ 1 mmdays
R10mmNumber of heavy precipitation daysThe annual count of days when PRCP ≥ 10 mmdays
R1mmNumber of days above 1 mmThe annual count of days when PRCP ≥ 1 mmdays
SDIISimple daily intensity indexAnnual total precipitation divided by the number of wet days (defined as PRCP ≥ 1.0 mm) in the yearmm/day
Table 5. Comparison of several drought events calculated from observed precipitation data and raw and bias-corrected GCM data for 1976–2005. Values in parentheses refer to the standard deviation of the results over the GCMs.
Table 5. Comparison of several drought events calculated from observed precipitation data and raw and bias-corrected GCM data for 1976–2005. Values in parentheses refer to the standard deviation of the results over the GCMs.
StationsModerate DroughtSevere DroughtExtreme Drought
ObservedHistorical (Raw)Historical (Bias-Corr.)ObservedHistorical (Raw)Historical (Bias-Corr.)ObservedHistorical (Raw)Historical (Bias-Corr.)
Rangpur4352 (±12.3)58 (±13.5)38 (±4.2)8 (±2.8)22 (±1.7)1 (±1.2)
Dinajpur4850 (±13.2)58 (±13.3)117 (±4.4)7 (±2.7)42 (±1.8)1 (±1.2)
Bogra5450 (±12.0)51 (±13.3)118 (±3.7)9 (±2.9)32 (±1.8)1 (±0.8)
Rajsahi5546 (±11.0)56 (±11.9)108 (±4.2)7 (±3.5)32 (±1.6)1 (±1.4)
Ishwardi6048 (±08.7)55 (±11.2)58 (±4.0)8 (±3.1)22 (±1.6)1 (±1.2)
Jessore8647 (±10.1)56 (±10.2)109 (±4.2)8 (±3.4)22 (±1.6)2 (±1.1)
Khulna6549 (±11.4)52 (±10.8)139 (±4.8)8 (±3.5)22 (±1.5)1 (±1.0)
Satkhira6546 (±10.6)54 (±11.3)139 (±4.2)8 (±3.5)53 (±1.8)2 (±1.4)
Barisal7149 (±11.4)51 (±15.1)209 (±5.0)9 (±3.1)42 (±1.6)2 (±1.2)
Patuakhali2148 (±11.8)57 (±12.1)29(±4.0)9 (±3.4)23 (±1.7)2 (±1.4)
Khepupara6048 (±12.4)54 (±11.0)119 (±3.7)9 (±3.4)73 (±1.7)2 (±1.3)
Bhola4348 (±12.2)54 (±12.8)810 (±5.1)9 (±3.7)22 (±1.7)2 (±1.3)
Maizdicourt6852 (±13.7)54 (±13.4)169 (±5.1)9 (±4.1)12 (±1.4)2 (±1.0)
Sawndip6154 (±16.6)54 (±17.1)48 (±4.1)8 (±3.3)33 (±1.8)2 (±1.4)
Dhaka10247 (±10.9)57 (±10.5)109 (±4.2)8 (±3.2)32 (±1.9)2 (±1.3)
Mymensingh4951 (±14.8)57 (±13.4)128 (±4.4)8 (±3.2)22 (±1.6)2 (±1.0)
Hatiya6151 (±15.7)54 (±09.7)29 (±4.6)8 (±3.6)12 (±1.7)2 (±1.1)
Chandpur1748 (±11.0)54 (±11.5)09 (±5.0)7 (±3.9)12 (±1.5)1 (±0.9)
Comilla6452 (±13.6)54 (±13.1)119 (±4.8)8 (±3.3)42 (±1.5)1 (±1.0)
Feni4555 (±14.3)56 (±15.5)69 (±4.9)9 (±3.9)52 (±1.8)2 (±1.3)
Sylhet8454 (±14.4)51 (±11.8)98 (±4.6)9 (±3.6)32 (±1.4)2 (±1.3)
Srimongal5354 (±13.7)57 (±12.8)58 (±4.1)8 (±3.5)22 (±1.7)2 (±1.3)
Rangamati5453 (±15.9)58 (±17.3)98 (±4.1)7 (±2.7)32 (±1.7)2 (±1.3)
Cox’s bazar1946 (±12.9)56 (±13.3)310 (±3.7)9 (±4.2)12 (±1.5)2 (±1.5)
Chittagong2452 (±15.6)58 (±15.5)09 (±4.6)8 (±4.0)02 (±1.9)2 (±1.2)
Bangladesh54.8849.9455.048.168.608.22.682.291.68
In percent (%) Under-estimate 9%Match Over-estimate 5%Match Under-estimate 15%Under-estimate 37%

Share and Cite

MDPI and ACS Style

Kamruzzaman, M.; Jang, M.-W.; Cho, J.; Hwang, S. Future Changes in Precipitation and Drought Characteristics over Bangladesh under CMIP5 Climatological Projections. Water 2019, 11, 2219. https://doi.org/10.3390/w11112219

AMA Style

Kamruzzaman M, Jang M-W, Cho J, Hwang S. Future Changes in Precipitation and Drought Characteristics over Bangladesh under CMIP5 Climatological Projections. Water. 2019; 11(11):2219. https://doi.org/10.3390/w11112219

Chicago/Turabian Style

Kamruzzaman, Mohammad, Min-Won Jang, Jaepil Cho, and Syewoon Hwang. 2019. "Future Changes in Precipitation and Drought Characteristics over Bangladesh under CMIP5 Climatological Projections" Water 11, no. 11: 2219. https://doi.org/10.3390/w11112219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop