Near-term regional climate change over Bangladesh

Bangladesh stands out as a climate change hot spot due to its unique geography, climate, high population density, and limited adaptation capacity. Mounting evidence suggests that the country is already suffering from the effects of climate change which may get worse without aggressive action. Here, we use an ensemble of high-resolution (10 km) regional climate model simulations to project near-term change in climate extremes, mainly heat waves and intense rainfall, for the period (2021–2050). Near-term climate projections represent a valuable input for designing sound adaptation policies. Our climate projections suggest that heatwaves will become more frequent and severe in Bangladesh under the business-as-usual scenario (RCP8.5). In particular, extremes of wet-bulb temperature (a temperature and humidity metric important in evaluating humid heat stress) in the western part of Bangladesh including Bogra, Ishurdi, and Jessore are likely to exceed the extreme danger threshold (according to U.S. National Weather Service criterion), which has rarely been observed in the current climate. The return periods of extreme heat waves are also significantly shortened across the country. In addition, country-averaged rainfall is projected to increase by about 6% during the summer months, with the largest increases (above 10%) in the eastern mountainous areas, such as Sylhet and Chittagong. Meanwhile, insignificant changes in extreme rainfall are simulated. Our results suggest that Bangladesh is particularly susceptible to climate extremes in the near future, in the form of extreme heat waves over the western part of the country.


Introduction
Scientific evidence has shown that rising anthropogenic greenhouse gas (GHG) concentrations due to human activity are unequivocally accelerating global warming (IPCC 2013). As global warming steadily progresses since the industrial revolution, weather-related extremes, like heat waves and intense rainfall, are increasing in likelihood, threatening human life and causing socio-economic damage worldwide (Meehl and Tebaldi 2004;Kharin et al. 2013;Sillmann et al. 2013). Developing and financially challenged countries are especially vulnerable to the impacts of climate change due to their limited resources and low adaptive capacity, and some of these countries are already facing a climate crisis (Huq 2001;Kreft et al. 2017). Therefore, it is critical to project the impacts of near-term climate change in the most vulnerable regions, such as Bangladesh (Huq 2001;Kreft et al. 2017), in order to inform sound adaptation strategies.
Bangladesh, one of the world's most impoverished countries, stands out as a climate change hot spot because of its weak infrastructure, geographic conditions and dense population (more than 1100 inhabitants per square kilometer shown in Fig. 1; Huq 2001;World Bank 2013;Kreft et al. 2017). The climate in Bangladesh is typically humid due to its proximity to the Bay of Bengal, and hence is prone to the impacts of extreme weather events including heat stress and heavy rainfall (Khan et al. 2000;Shahid 2010Shahid , 2011Shahid et al. 2016). Since the country belongs to the Asian monsoon regime, predominant southerly flows during summer monsoon season can further worsen the risk of heat stress by transporting warm and moist air masses from the Bay of Bengal to Bangladesh (Shahid et al. 2016). Given the aforementioned factors, the country is particularly vulnerable to the impacts of climate change, and the risks associated with these impacts will likely be exacerbated as climate change becomes more pronounced (Huq 2001).
As exemplified by recent climate records, heat and rainfall extremes are already deadly in Bangladesh. According to modern record , record-shattering heat waves impacted South Asian countries (e.g., India, Pakistan and Bangladesh) along the Indus and Ganges river valleys (Im et al. 2017a;Raymond et al. 2020). Recently, a maximum wet-bulb temperature above 31 °C was recorded in Bangladesh, which is extremely dangerous for human health (Im et al. 2017aKang et al. 2019b;Sherwood and Huber 2010;Raymond et al. 2020). The 2015 South Asian heat wave was estimated to have caused 4,000 deaths as a consequence of heat-related diseases (UN ESCAP 2015;Ghumman and Horney 2016;Zommers et al. 2016). In addition to temperature extremes, devastating floods resulted in a death toll over 1,000 and affected an estimated 40 million people over South Asia in 2017 (Gettleman 2017). Bangladesh is frequently threatened by devastating cyclones with intense rainfall every three years on average (MoEF 2009). Without serious mitigation measures, Bangladesh will face an increasing risk of more severe impacts in the coming decades.
Based on future climate projections using General Circulation Models (GCMs), southwest Asian countries are likely to be exposed to deadly weather extremes including heat waves and extreme heavy rainfall in the future (OECD 2003;May 2004;Sillmann et al. 2013;Kirtman et al. 2013;Coffel et al. 2018). Since the vulnerability to climate extremes is largely dependent on geographical location and topographical characteristics, GCMs with coarse resolution and inadequate representation of physical processes have limitations in simulating the fine scale climate extremes (Seneviratne  Flato et al. 2013;Im et al. 2017b). Dynamically downscaled simulations using a regional climate model facilitate in-depth analysis of extremes over Bangladesh. Most previous regional climate modeling studies were confined to the far future period (Gu et al. 2012;Im et al. 2017a) close to the end of the century (2071-2100), and have analyzed the likelihood of extremes based on maximum temperature (e.g., Gu et al. 2012).
In our earlier study by Im et al. (2017a), we projected changes in wet-bulb temperature over south Asia towards the end of the century using a regional climate model at 25-km resolution. In this study, we investigate the changes in extreme heat stress and rainfall in response to the business-as-usual near-term climate change scenario based on ensemble of projections using the high-resolution (10 km) MIT Regional Climate Model (MRCM), allowing a more detailed representation of topography, coastlines, physical processes, and extreme climatic events.

Regional climate model description
The MRCM is a three-dimensional hydrostatic numerical model, which was developed at MIT to understand how global climate and land-use changes may impact society (e.g., Im et al. 2017a). The general structure and dynamical core of the MRCM originate from the Abdus Salam International Centre for Theoretical Physics (ICTP) Regional Climate Model version 3 (RegCM3; Pal et al. 2007). Characteristically, the Integrated Biosphere Simulator (IBIS; Foley et al. 1996) is coupled to the MRCM to obtain better representation of land surface processes (Winter et al. 2009). In addition to coupling with the IBIS, the MRCM notably employs several improvements of model physics, such as a new surface albedo assignment (Marcella 2012;Marcella and Eltahir 2012), an irrigation scheme (Marcella 2012;Marcella and Eltahir 2014), a new convective cloud, and rainfall autoconversion schemes (Gianotti 2012;Gianotti andEltahir 2014a, 2014b), and a modified boundary layer height and boundary layer cloud scheme (Gianotti 2012). The MRCM has been extensively evaluated in its ability to simulate the fine-scale structure of climatic variables over various domains, such as North America (Winter et al. 2009), West Africa (Im et al. 2014), Southwest Asia (Pal and Eltahir 2016), South Asia (Im et al. 2017a), and the Maritime Continent . Pal and Eltahir (2016) have improved the capability of the MRCM through prescribing the satellite-derived albedo and adjusting the associated parameters. They showed that the upgraded version of MRCM captured the key observed climate features over the desert and semi-desert region, particularly Southwest Asia. Im et al. (2017a) showed that the MRCM could suitably reproduce the spatial pattern and magnitude of temperature fields (e.g., dry-bulb and wet-bulb temperatures) over the South Asia region. Given the satisfactory performance of the model in this recent study over South Asia, the MRCM is applied in this study to simulate regional climate over Bangladesh.

Optimization of the model parameters
Despite various updates and improvements of the MRCM, the latest version of the MRCM with default settings still has limitations in reproducing region-specific climate information. Thus, it is necessary to optimize the model performance for the current domain before projecting future climate. Several parameters in the cloud microphysics and convection schemes are updated bringing them close to the values recommended in recent version of RegCM version 4 (Giorgi et al. 2016): (1) The auto-conversion rate is 0.50 × 10 -3 s −1 ; (2) accretion rate is 6 m 3 kg −1 s −1 ; (3) raindrop evaporation rate is 1.0 × 10 -5 (kgm −2 s −1 ) −1/2 s −1 ; (4) maximum negative parcel temperature perturbation is 0.9 K; and (5) approach rate to quasi-equilibrium is 0.2 kgm −2 K −1 s −1 . Furthermore, we adjusted the parameters of cloud liquid water based on the study of : (6) climatological cloud liquid water in continental and maritime environments are set to 1.0 g m −3 and 0.4 g m −3 , respectively; (7) Threshold of cloud liquid water is 1.1 g m −3 in continent, and 0.45 g m −3 in maritime.

Experimental design
The simulation domain covers South Asia centered in Bangladesh with a horizontal resolution of 10 km on a Lambert conformal projection (Fig. 1). It consists of 84 × 81 × 18 grid points along the zonal, meridional, and vertical (sigma coordinate) directions, respectively. To minimize interference from the elevated topography located in the vicinity of the lateral boundary zone, we choose the slightly larger domain which fully covers Bangladesh as well as the surrounding mountainous area. We use this setup for evaluation runs and climate simulations.
Before applying the regional climate modeling system to project future climate, it is important to assess the capability of the model in reproducing the present climate conditions in terms of the spatial distribution as well as the seasonal cycle of climatic variables (i.e., temperature and precipitation) to identify possible systematic biases, present in the model. In this regard, we conducted an evaluation run (i.e., known as perfect boundary simulation) using the MRCM driven by 6-hourly European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim) with a resolution of 1.5° (Uppala et al. 2008). The 1° × 1° weekly NOAA Optimum Interpolation Sea Surface Temperature (OISST) v2 data (Reynolds 2002) is prescribed as the lower boundary condition of the MRCM. The simulation period of evaluation run is from 1982 to 2012, because NOAA OISST data (i.e., sea surface boundary condition) is available from 1982, and at least a 30-year simulation is required for climate study.
To test the necessity of 10-km high-resolution simulation over Bangladesh where the topography is almost flat with an elevation of less than 50 m, we additionally conducted low-resolution (30 km) MRCM simulation for the 10-year period (1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992) using the same model configuration (e.g., physics scheme, boundary condition from the ERA-Interim reanalysis). By comparison with the low-resolution result, 10-km MRCM simulation can provide added value by better representing spatial distribution of temperature and precipitation particularly along the coastal and eastern mountainous regions (Suppl. Figure 1). This result is consistent with recent studies that demonstrated the advantages of high-resolution simulations (Ciarlo et al., 2020;Qiu et al., 2020).
Through a rigorous screening process, Im et al. (2017a) recommended three Coupled Model Inter-comparison Project Phase 5 (CMIP5) GCMs [Community Climate System Model Version 4 (CCSM), Australian Community Climate and Earth System Simulator Version 1.0 (ACCESS), and Max-Planck-Institution Earth System Model running on Medium Resolution grid (MPI)] to use as lateral boundary conditions of the MRCM climate simulation over South Asia. They found that the MRCM simulations driven by these GCMs are able to capture local climate characteristics (e.g., seasonal mean climatology and extreme temperatures) of South Asia. In accordance with this recommendation, the same three GCMs are selected in this study as driving boundary conditions for the MRCM to project future climate change and its impacts. The general information describing these CMIP5 models is summarized in Table 1. We performed three reference simulations for the period 1975-2005 with the historical greenhouse gas (GHG) concentration and conducted three time-slice simulations under the businessas-usual anthropogenic emissions scenario, namely, RCP8.5, for the period 2020-2050. The first year of each simulation is regarded as the spin up period, and thus are not considered in the analysis. We assumed that the skill of the model under current and future climate conditions is the same (Suh et al. 2016).

Bias correction
All simulations by a regional climate model suffer from systematic bias arising from lateral boundary conditions and their imperfect dynamics and physical processes (Liang et al. 2008;Ehret et al. 2012). In order to reduce the aforementioned bias, we applied the Parametric Quantile Mapping (PQM) technique. This method is commonly used in climate change studies to correct for the structural biases in the rainfall and temperature fields although several problematic aspects are possible: parametric quantile mapping may artificially corrupt the modeled climate change signals (Cannon et al., 2015;Switanek et al., 2017). High-resolution Terraclimate data (horizontal resolution is about 4 km; Abatzoglou et al. 2018) is used to correct for biases in drybulb temperature and precipitation in the model output, and ERA-interim is used to correct the bias of wet-bulb temperature for the entire domain at the monthly time scale. Also, weather observations at 7 major cities ( Fig. 1; Dhaka, Jessore, Bogra, Sylhet, Ishurdi, Chittagong, and Cox's Bazar) in Bangladesh are applied to correct the mean biases of MRCM simulations at daily time scale. To calculate the daily maximum wet-bulb temperature at stations, daily minimum relative humidity is estimated using the following formulation (Allen et al. 1998).
where RH min is the daily minimum relative humidity, T max the daily maximum temperature, e a the actual vapor pressure of air, e s (T max ) saturated vapor pressure of air at maximum temperature. For bias correction, we specifically applied the method following Piani et al. (2010) and Kang et al. (2019aKang et al. ( , 2019b observed climatic variables (i.e., rainfall, dry-bulb and wetbulb temperatures) for each month with the assumption that temperatures and rainfall distributions are well approximated by the gaussian and gamma distributions, respectively. The modeled cumulative frequency distribution for the variable of interest, simulated for the past, is then corrected to fit the observed distribution. The same correction factors are then applied to correct the bias in the simulated future variable considered.

Evaluation of "perfect boundary" simulation
We assume that the MRCM has sufficient skill for use in this study to project climate change over Bangladesh based on comparison with temperature and rainfall observations (Figs. 2 and 3). Figure 2 shows the observed and simulated seasonal climatology of temperature over Bangladesh during the reference period . The MRCM adequately reproduces the mean surface temperature distribution over the target region in all seasons, as compared with the corresponding observations. For pre-monsoon season (March, April, May (MAM)), the mean temperature shows a distinct gradient from Southwest to Northeast, warmest in Southwest coastal region and coolest near the Northeast mountainous area. Meanwhile, a south to north gradient is evident across the region in winter. The observed temperatures are high from monsoon to post-monsoon seasons, decreasing markedly as winter approaches. All these features are well captured by the MRCM, although a cold bias exists in certain seasons. During post-monsoon season, the model shows a noticeable cold bias of up to 2-3℃ over most of the region. This systematic discrepancy has also been reported by previous studies (e.g., Rahman et al. 2012). The MRCM reasonably reproduces the seasonal variation of rainfall over Bangladesh as compared to TerraClimate observation. Figure 3 presents seasonal mean rainfall derived from TerraClimate observation and MRCM simulation. The observed rainfall exhibits a west-east gradient, and it peaks in monsoon season with more than 60% of the total annual rainfall in Bangladesh. According to the geographical setting, high rainfall is observed along the high mountainous regions, and relatively dry conditions are observed in the low-lying plain regions. These features are well captured by high-resolution MRCM simulation, although it suffers from some systematic biases across the region. The MRCM tends to underestimate mean rainfall, especially in the low elevation regions during monsoon and post-monsoon seasons, meanwhile it shows a large overestimation of mean rainfall along the northeastern parts of Bangladesh for premonsoon and monsoon seasons. The magnitude of the bias is comparatively higher in the mountainous regions. Some of these biases are expected due to the difficulty representing dynamics over mountainous regions with sharp elevation gradients, in regional climate models .

Evaluation of reference simulation
As discussed in the previous sections, the MRCM has been tested in several climate change studies. The next steps in testing the reliability of future projections is to identify systematic errors apparent in historical climate simulations (i.e., present-day climate simulations) and to remove the structural biases using the statistical bias correction method (i.e. PQM method; see Sect. 2.4). It is assumed that model bias does not change with time.
Here, we evaluate the reliability of the MRCM ensemble simulations. Figure 4 illustrates the spatial distribution of mean dry-bulb (TA) and wet-bulb (TW) temperatures during June-September (JJAS) derived from observation and model biases. The MRCM ensemble reasonably captures the spatial patterns of dry-and wet-bulb temperatures during monsoon season, although it overall underestimates the temperatures showing a noticeable cold bias with amounts greater than 0.6 ℃ in the southwestern low-lying region (Fig. 4b, e). In addition, MRCM simulations driven by the three different models consistently overestimate the temperatures along the southeastern border of the country. We also show that structural model biases are effectively decreased through the bias-correction step with the PQM method. As expected, the bias-corrected MRCM ensemble shows better agreement with observed climatological patterns in terms of spatial distribution and magnitude.
The spatial distributions of mean rainfall during the whole year and monsoon season are shown in Fig. 5. As discussed in the previous section, the observed rainfall is characterized by a pronounced gradient between low rainfall in the low-elevation plain and high rainfall in mountainous areas (Fig. 5a, d). The ensemble mean of the MRCM simulations tends to underestimate rainfall in the low-lying coastal region. Moreover, it noticeably overestimates rainfall over northeastern regions. The spatial pattern of rainfall bias from the reference simulation is similar to that from the perfect boundary forcing simulation driven by ERA-Interim (Fig. 3), which can thus be interpreted as the MRCM's systematic bias. These systematic biases are in line with the results of previous studies and can be effectively eliminated through the PQM method (Kang et al. 2019a(Kang et al. , 2019b. The bias-corrected rainfall derived from the MRCM ensemble agrees well with observation, with bias less than 1 mm/day. The above results, in terms of systematic temperature and rainfall biases, provide an overview of the MRCM performance on

Projected changes in near-term temperature
Based on the MRCM ensemble projection, heat stress over Bangladesh is expected to worsen in the near future as climate change progresses. Figure 6 shows the changes in dryand wet-bulb temperatures under the business-as-usual scenario. The ensemble mean of the MRCM simulations shows a generalized warming of mean temperature with magnitude above 1.0 ℃ across the country during monsoon season. The largest increase is expected in the western part of the country including Dhaka, Bogra, Ishurdi, and Jessore (Fig. 6c). Along with rising temperature, significant increases in wetbulb temperature (increase by average of about 1.1 ℃) are also identified in the overall region, although the spatial pattern of change in wet-bulb temperature differs from that for the dry-bulb temperature. That is, the increase in wet-bulb temperature is spatially homogeneous, with no major difference depending on the region, except for some coastal and mountainous regions (e.g., Cox's Bazar; Fig. 6f). To quantitatively examine the projected risk of heat stress, we applied the U.S. National Weather Service (USNWS) heat index risk categories, including the danger and extreme danger thresholds which approximately correspond to a wet-bulb temperature of 27 ℃ and 31.2 ℃, respectively, under the specific climatic condition of Bangladesh (Suppl. Figure 2; Kang and Eltahir, 2018;Kang et al., 2019b). This approach, defining the thresholds of heat stress using the impacts of temperature and humidity, is in line with Nissan et al. (2017). They found significant effects of day-and nighttime temperatures and humidity in shaping heat-related deaths and proposed a heat-wave definition for Bangladesh that could inform the development of a heat early warning system. The mean wetbulb temperature is expected to exceed the USNWS danger threshold of 27 °C, primarily in the interior region, which implies that dangerous heat stress will become normal in most low-lying plains over Bangladesh in next few decades.
Consistent with the increases in mean dry-and wet-bulb temperatures, we also project the amplification of heat stress seasonally for the near-term future. Figure 7 shows the projected changes in annual cycle of maximum wet-bulb temperature for seven major cities. In a present-day climate, the daily maximum wet-bulb temperature is approaching or has already exceeded the USNWS danger threshold (27 ℃) over most stations during monsoon season. As the century progresses, more warming is projected for the monsoon season as compared to the other seasons. Especially, increases in mean maximum wet-bulb temperature during monsoon season are expected to be well above the danger threshold across all target stations under the business-as-usual scenario. Of particular note, the period of suffering from the heat stress is expected to be longer than those we observed in the current climate. Western Bangladesh cities, including Jessore, Bogra, and Ishurdi, where heat-related risk is already high, will experience more extreme heat waves in the near-term future.
To qualitatively assess the change in extreme heat waves events, we also analyze the probability distribution function of daily maximum wet-bulb temperature (Fig. 8). For the historical period, the daily maximum wet-bulb temperature is almost never observed more than the 31.2 °C extreme danger threshold except for Jessore. However, under the business-as-usual scenario, the extremes of daily maximum wet-bulb temperature are projected to exceed the extreme danger threshold (31.2 °C) for many urban population centers, including Dhaka, Jessore, Bogra, Sylhet, Ishurdi, Chittagong, and Cox's Bazar, by the middle of the century. Particularly, western Bangladesh cities, including Jessore, Bogra and Ishurdi, are projected to experience extreme heat stress episodes with a maximum wet-bulb temperature exceeding 32 °C and approaching 33 °C, which is considered extremely dangerous for most humans. While the general shape of the distribution of maximum wet-bulb temperature in the seven target cities does not change dramatically, the RCP8.5 projection shows the distributions shift to higher temperatures relative to the present climate with higher frequency of threshold exceedance.
Analyses of return period (known as recurrence interval) preliminarily suggest that the enhanced GHG forcing under the RCP8.5 scenario will increase the likelihood of the reoccurrence of extreme heat stress. Figure 9 shows return periods for annual maxima of wet-bulb temperature at several stations derived from bias-corrected values from all MRCM ensemble members. Results indicate significant decreases in the return periods across the region with the RCP8.5 scenario. For the current climate condition, the return periods of the dangerous heat stress events (maximum wet-bulb temperature is above USNWS extreme danger threshold, 31.2 °C) have a recurrent period of longer than 50 years. However, these extremes are likely to become a 1-in-3 to 1-in-20 year event in the 7 major cities by the middle of the twenty-first century. In particular, the impact of near-term climate change would lead to a stronger projected decrease in return period in the western regions, in particular, Bogra, Ishurdi, and Jessore. These cities can be considered hotspots because they are likely to experience more frequent extremes of wet-bulb temperature (Fig. 8), with the greatest reduction in return period of extreme heat at less than 10 years. Figure 10 shows the projected change in mean rainfall during annual and monsoon seasons. Based on bias-corrected MRCM ensemble projections, annual mean precipitation generally increases by about 4% in Bangladesh in a warmer climate. Particularly, a small but significant increase in annual mean rainfall is projected mostly in Northern region (e.g., Bogra and Sylhet) of Bangladesh with a magnitude of more than 6% as compared to the reference period (Fig. 10c). In addition, our ensemble projection shows a significant increase of monsoon rainfall, by about 6%, over the target The 30-year mean seasonal cycle of TWmax for the seven cities (a Dhaka, b Jessore, c Bogra, d Sylhet, e Ishurdi, f Chittagong, g Cox's Bazar) over Bangladesh derived from station observation (1984-2013; black), and ensemble mean of MRCM simulations under historical (1976MRCM simulations under historical ( -2005blue) and RCP8.5 scenarios (2021RCP8.5 scenarios ( -2050red  , and RCP8.5 (2021-2050) scenarios, respectively. Vertical dashed red line indicates the extreme danger (31.2 ℃) threshold defined by the U.S. National Weather Service region, with the largest increases (above 10%) along the eastern mountainous areas (e.g., Sylhet and Chittagong). The increase in annual rainfall is mostly a result of the monsoon rainfall increases. The consistent rainfall increases among the simulations are also notable in some areas, such as Western, Northeastern, and Southeastern regions of Bangladesh.

Projected changes in near-term rainfall
Using our high-resolution projections to assess the potential contribution of future climate change on the seasonal cycle of rainfall at the local scale, we find that the Eastern cities, Sylhet and Chittagong, show general increases in monsoon season rainfall. Figure 11 presents the projected changes in rainfall seasonality for major cities. In general, rainfall reaches an annual maximum in the monsoon season, peaking in June/July and then gradually decreases until the winter season (Figs. 3 and 11). Particularly, rainfall in three sub-regions, Sylhet, Chittagong, and Cox's Bazar, shows strong seasonal variation with a maximum of up to 25 mm/day in monsoon season. Whereas there is relatively weaker seasonality in the western cities including Dhaka, Jessore, Bogra, and Ishurdi. In the future projection, increases in rainfall during the monsoon season appear in most target cities, though small. Eastern cities, such as Sylhet and Chittagong will likely experience increased rainfall of around 14%, during monsoon season in the coming few decades. The consistency among all simulations driven by the three different GCMs gives high confidence in the results. The change in the incidence of extreme rainfall events is further examined through return period analysis (Fig. 12). Based on the MRCM projections with bias correction, the return period of annual maximum precipitation does not change significantly in response to the RCP8.5 GHG forcing despite the increase in mean precipitation over the region. Most cities show slight decreases in return periods, which are, however, not significant except for Cox's Bazar. These weak changes reinforce the insignificant impact of near-term climate change on the rainfall extremes. Coloured dots c, f indicate the future changes at 7 meteorological stations (Dhaka, Jessore, Bogra, Sylhet, Ishurdi, Chittagong, and Cox's Bazar). Stippling indicates agreement by three MRCM simulations on the sign of the change. Area-averaged value over land are given on the bottom left corner of each plot (unit: mm/day) Fig. 11 The 30-year mean seasonal cycle of rainfall for the seven cities (a Dhaka, b Jessore, c Bogra, d Sylhet, e Ishurdi, f Chittagong, g Cox's Bazar) over Bangladesh derived from station observation (1979-2008; black), and ensemble mean of MRCM simulations under historical (1976-2005; blue)

Summary and conclusions
Many parts of the world experience climate extremes, such as heat waves and heavy rainfalls. Among the countries, Bangladesh is considered one of most vulnerable countries to the adverse impacts of climate change. It is generally agreed that the projected increase in temperature and humidity will exacerbate the risk of extreme climate events in the country over the next few decades (e.g., Im et al. 2017a). However, there is lack of understanding of the direction, magnitude and spatial pattern of projected climatic changes at regional to local scales.
In this study, we used an ensemble of high-resolution regional climate model simulations driven by three CMIP5 GCMs to project near-term changes in temperature and rainfall extremes across Bangladesh under the business-asusual scenario. We showed that the MRCM ensemble with statistical bias correction is capable of capturing regionspecific climate information over the target domain.
Under the business-as-usual scenario we show that heatwaves will become more frequent and severe in Bangladesh. It is also projected that the western part of Bangladesh, including Bogra, Ishurdi, and Jessore, is particularly susceptible to extreme heat stress with the maximum wetbulb temperature exceeding the extreme danger threshold (31.2 ℃ defined by USNWS), which was rarely observed over this area in the current climate. Furthermore, the enhanced GHG forcing under the RCP8.5 scenario will increase the likelihood of the reoccurrence of extreme danger level heat extremes. However, changes in rainfall, unlike temperature changes, are small across the region, with a slight increasing tendency. Area-averaged annual mean and monsoon season rainfall in Bangladesh are expected to slightly increase by about 4% and 6%, respectively, with the maximum centered in eastern mountainous regions during the summer months. Furthermore, the return period of annual maximum precipitation will experience small changes in Cox's Bazar, however the same variable is likely to remain close to its present climate value in most regions.
Our study further illustrates that Bangladesh faces climate change associated risks beyond the risks in coastal areas. It is generally agreed that coastal regions of Bangladesh are particularly susceptible to the impacts of sea level rise, the extreme cyclones (Ali 1999;OECD 2003;Dasgupta et al. 2014). However, using an ensemble of projections, we also show that inland locations are likely to be threatened by deadly heat waves, especially in localized hotspots in western parts of Bangladesh. Therefore, populations moving inland due to coastal climate related disasters will still be faced with extreme heatwaves away from the coast. This study provides a crucial glimpse of the near future climate, with projections based on high-resolution MRCM simulation that provide significant information to help inform the development of sound adaptation strategy.