2.1 Study area description
Bangladesh is one of the world's most drought-prone countries. The most severe drought in Bangladesh's history happened in the western half of the country. It is separated into two parts: The north-western region (240.30ꞌ- 26040ꞌ N, 880.01ꞌ-890.90ꞌ E) and the south-western region (220.52ꞌ-230.90ꞌ N, 880.20ꞌ-90030ꞌE) (Islam et al. 2017). The northwestern stations are Rangpur Bogra, Rajshahi and southwestern stations are Faridpur, Satkhira, Khulna, Barishal. Aman rice production is the highest in this region of the country. With an average rainfall of less than 1500 mm annually and summertime humidity of less than 50%, Bangladesh's western zone is the driest. In the summer, it is the warmest climatic zone. The typical maximum temperature in summer exceeds 35°C, while rainfall is substantially lower. It also gets a lot less rainfall than the rest of the country. The excesses of the zones to the north are softened in the south-western zone. The annual rainfall varies from 1,500 to 1,800 mm (Islam et al. 2019).
2.2 Data source and quality check
Climate data
Seven stations in Bangladesh's western area provided daily meteorological data (Fig. 1). These seven stations have been selected as the study area because of their climatic characteristics. The study area was chosen for a variety of reasons, including rainfall patterns, temperature variability, hydrological characteristics, and historical drought records. This study used daily maximum, minimum, and average temperatures, solar radiation, sunshine hours, rainfall, static pressure, and wind velocity from 1980 to 2017. It should be emphasized that the research is based on data from the period 1980 to 2017. The information was obtained from the Bangladesh Meteorological Department (BMD). The quality control of these datasets was verified using quality control tests such as the Chi-Squire test and the Grub test. BMD staff confirmed the accuracy and validity of the dataset. In the case of missing data on rainfall, temperature, and humidity, those were collected from the nearest station, and a well-developed method was applied to make them accurate.
Crop data
To get information on the Aman rice growing season, daily crop data such as farmed area, crop yield, and pheno-phase were gathered from the Bangladesh Bureau of Statistics (BBS 2019). Phenophase data of Aman rice was collected from agrometeorological stations in Bangladesh throughout the research period, while yield data was collected from the crop database in the Statistical Year Book Bangladesh-2018.
Socio-economic data
The socio-economic data from the period of 1980 to 2011, including GDP, population, total power of agricultural machines, effective irrigation area, disaster area, and plowed area, are collected from the statistical yearbooks of the Census (2011) of Bangladesh.
2.3 Calculation of CWDI
This study uses the crop water deficit index (CWDI) to detect drought during the Aman rice growing season. The Agricultural Water Deficit Index (CWDI) is an effective method for determining deficits in crop water and it's frequently used in agricultural drought analyses (Yang et al., 2017; Jhajharia et al., 2012; Gao et al., 2019). The CWDI is determined by subtracting the crop water requirement from the water demand at the very same phase (Gao et al., 2019; Hu et al., 2021a). Eq. (1) shows how to calculate the CWDI for a discrete 10-day period:
$${CWDI}_{i}=\left\{\begin{array}{c}\frac{{ET}_{i}-{P}_{i}}{{ET}_{i}}\times 100\% {ET}_{i}\ge {P}_{i}\\ 0 {ET}_{i}<{P}_{i}\end{array}\right. \left(1\right)$$
Where CWDIi is the total water deficit index for the ith 10-day period (%), ETi is the rice collective water need for the ith 10-day period (mm), and Pi is the cumulative precipitation for the ith 10-day period (mm), which is equivalent to crop watersupply. Five consecutive 10-day weighted CWDIs make up the total CWDI. The CWDI is determined by the quantity of rainfall and the volume of water required.The CWDI has several benefits, including the ability to properly show drought from the perspective of Aman rice development stages and the ability to be used in humid areas.
As a result, the noncumulative CWDI was chosen to expose the combined effects of plants and agricultural water demand is influenced by meteorological conditions.
ETi was computed with the crop coefficient-ET0 method by the given equation:
$${ ET}_{i}={K}_{c}\times {ET}_{0} \left(2\right)$$
The crop coefficient is Kc, while the reference evapotranspiration is ETo.
The FAO P-M (Penman-Monteith) equation, given by the FAO (Food and Agriculture Organization), was used to calculate ETo (Allen et al. 1998).
3
Where ET0 is the reference evapotranspiration (mm day-1), is the vapor pressure curve (k Pa° C-1), Rn is the net radiation at the crop surface (MJ m2 day-1), G is the soil heat flux density (MJ m2 day-1), is the psychrometric factor, which should be 0 according to the FAO standard, T denotes the daily average atmospheric temperature at a height of 2 meters (°C), u (kPa).
Based on the CWDI values, north-western part of Bangladesh can be divided into three sub-regions (Fig. 1). The drought classification criteria based on the CWDI of Aman rice are shown in Table S1.
2.4 Drought characteristic index calculation
Drought frequency
The drought onset throughout the crop growth period is expressed as drought frequency (DF). Using Eq. 4, F can be obtained.
DF = (n / N) × 100% (4)
Where N is the total year numbers with meteorological data, and n is the number of drought years.
Drought station ratio
The drought station ratio (Pj), which is used to assess drought-affected areas, is calculated as follow:
5
The overall number of climatic sites is M, the amount of drought sites is m, and the year is j..
Yield reduction rate
Eq. 6 is used to calculate the rate of rice yield reduction.
R = (CY / TY)× 100% (6)
where The yield decrease rate is R, the climatic yield is CY, and the time trend yield is TY. CY is a variable that represents the annual variation in wheat production due to climatic circumstances, and it may be expressed using Eq. 7.
CY i = AYi–TYi (7)
Where CY is for climatic yield, AY stands for actual wheat yield per unit area, TY stands for temporal trend yield, which can be estimated using the 5-year rolling average actual yield, and I stands for the year. Positive (negative) CY values indicate increases (decreases) in wheat production due to favorable (unfavorable) environmental conditions..
2.5 Drought disaster risk assessment model
Eq. 8 was used to compute the drought disaster risk index (DDRI) for the Aman rice growth period in western Bangladesh:
DDRI = (DH×DHw) + (DS×DSw) + (ED×EDw) – (PL×PLw) (8)
Here, DS is the sensitivity of Aman rice to drought; DH is the drought hazard, which is the frequency and intensity of drought; ED is the degree of exposure of Aman rice to drought; and PL is the regional production level, which is used to assess drought disaster prevention and reduction capabilities. The normalized indices are all DS, DH, PL, and ED values. The weights assigned to the DH, DS, ED, and PL are DHw, DSw, EDw, and PLw, respectively.
DH was determined with drought intensity and frequency using the methods of Zhang et al. (2016) and Islam et al. (2017). DS was reflected using the slope of the regression function between yield reduction rate and drought intensity. DH and DS were determined month by month, taking into account the crop's sensitivity to dryness throughout distinct growth periods. The ratio of Aman rice grown area to total cultivated area was used to illustrate ED. Per capita GDP, effective irrigated cultivated area, and total power of agricultural machinery were used to calculate PL. Standard treatments of DH, DS, ED, and PL were conducted, and the weights of DH, DS, ED, and PL were given using the analytic hierarchy process.
2.6 Spatial distribution analysis
In general, bayesian kriging, standard kriging, and inverse distance weighted kriging were used to generate geographic variability of AGD risk at various crop growth phases (IDW) (Islam et al. 2019; Zinat et al. 2020). When compared to alternative interpolation methods, the IDW model was chosen for spatial variability analysis for this study due to its acceptance and estimation accuracy (Islam et al. 2021).