Dynamics and drivers of land use and land cover changes in Bangladesh

Bangladesh has undergone dramatic land use and land cover changes (LULCC) in recent years, but no quantitative analysis of LULCC drivers at the national scale exists so far. Here, we quantified the drivers of major LULCC in combination with biophysical and socioeconomic observations at the sub-district level. We used Landsat satellite data to interpret LULCC from 2000 to 2010 and employed a Global Surface Water Dataset to account for the influences of water seasonality. The results suggest that major LULCC in Bangladesh occur between agricultural land and waterbodies and between forest and shrubland. Exclusion of seasonal waterbodies can improve the accuracy of our LULCC results and driver analysis. Although the gross gain and loss of agricultural land are large on the local scale, the net change (gross gain minus gross loss) at a country scale is almost negligible. Climate dynamics and extreme events and changes in urban and rural households were driving the changes from forest to shrubland in the southeast region. The conversion from agricultural land to standing waterbodies in the southwest region was mainly driven by urban household dynamics, population growth, distance to cities and major roads, and precipitation dynamics. This study, which is the first effort accounting for water seasonality and quantifying biophysical and socioeconomic drivers of LULCC at the national scale, provides a perspective on overall LULCC and underlying drivers over a decadal time scale and national spatial scale and can serve as a scientific basis for developing land policies in Bangladesh.


Introduction
Land use and land cover changes (LULCC) is the key topic in global change studies since they can alter regional and global climate through changing biophysical, biogeochemical, and biogeographical characteristics of the Earth system (Jain et al. 2013;Robinson et al. 2013;Xu et al. 2016). Understanding LULCC dynamics and drivers can help to Communicated by Jasper van Vliet Electronic supplementary material The online version of this article (https://doi.org/10.1007/s10113-020-01650-5) contains supplementary material, which is available to authorized users. better understand the LULCC processes and mechanisms to develop models and land policies for a country. LULCC is driven by anthropogenic activities, such as socioeconomic development. It is also influenced by biophysical factors, such as climate, terrain, and soil quality (Dewan et al. 2012). Globally, about 60% of land changes are associated with direct human activities and the rest with indirect drivers, such as climate change during the period 1982(Song et al. 2018. Specifically, 27% of global forest loss can be attributed to commodity production over the period 2001 to 2015 (Curtis et al. 2018). At the regional scale of Southeast Asia, it is reported that forest loss is related to the cropland expansion (Zeng et al. 2018a, b). However, the quantitative analyses on the LULCC drivers at the national scale are still limited, especially in developing countries, which have experienced greater LULCC in recent years (Dewan and Yamaguchi 2009;Sloan and Sayer 2015).
Bangladesh has experienced LULCC over the years due to population and economy growth, infrastructure expansion (Islam and Hassan 2011), and climate change (Rahman and Manprasert 2006). The agricultural land in Bangladesh is threatened by soil salinity, productivity loss (especially in the coastal areas) (Hossain 2015), and climate events, such as floods , which may devastate vegetation and man-made facilities and therefore cause LULCC at the local scale. In addition, despite the efforts made by local government and international programs, forest areas continue to decrease (Chowdhury and Koike 2010;Hasan et al. 2017;Rasul et al. 2004;Reddy et al. 2016) and may disappear in the next 30-40 years or earlier (Chowdhury and Koike 2010). Extensive water resources (in the form of ponds, natural depressions, lakes, canals, rivers, estuaries, and coastal areas) have contributed to the aquaculture industry expansion in recent decades (Gias 2005;Shamsuzzaman et al. 2017). The aquaculture in Bangladesh can be broadly divided into two types, freshwater aquaculture (mainly comprised of pond fish farming) and coastal aquaculture (mainly shrimp farming) (Gias 2005). The area of fish ponds in Bangladesh has increased from 2,655 km 2 in 2001-2002 to 3,700 km 2 in 2013-2014, while the area of coastal shrimp farms has nearly doubled from 1,414 km 2 in 2001-2002 to 2,753 km 2 in (Fisheries Resources Survey System 2003 due to increased demand of shrimps in the international market (Inderbitzin et al. 2010). In addition, the conversion of traditional paddy culture land to shrimp culture ponds is a wellestablished practice in the southwest coastal area of Bangladesh (Ali 2006;Khan et al. 2015).
With these observed LULCC in Bangladesh, a number of studies have focused on understanding the dynamics and drivers of LULCC and their impacts (Ali 2006;Dewan and Yamaguchi 2009;Hossain 2015;Hossain et al. 2001;Reddy et al. 2016;Soil Resource Development Institute 2013;Yesmin et al. 2014). However, most of these local and regional-scale studies cannot be used to quantify the relations between LULCC and the driving factors for the entire country, which has 8 divisions, 64 districts (zila), and 484 sub-districts (upazila) with different socioeconomic conditions (one zila consists of approximately 8 upazilas, and the average size of an upazila is~300 km 2 ). Thus, the objective of our study is to improve our understanding of the dynamics and drivers of LULCC at the national scale, which is of great importance for managing future challenges in LULCC. We especially address the following key questions regarding LULCC in Bangladesh: (1) What are the major LULC conversions between 2000 and 2010? (2) What are the major biophysical and socioeconomic drivers of these conversions?
To address these questions, we (1) quantify LULCC at national scale using a wall-to-wall analysis of 30 m resolution Landsat imageries at a decadal time interval (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010), (2) compile the spatial biophysical drivers from global gridded datasets and socioeconomic drivers from national census data at sub-district level, (3) apply logistic regression to study the relationships between major LULCC and the drivers, and (4) evaluate and reinforce our findings on spatial drivers using the synthesis of case studies that incorporate field knowledge of the dynamics and drivers of LULCC in Bangladesh.

Materials and methods
The overall approach to study the dynamics and drivers of LULCC is broken down into six steps (Fig. 1). We describe these steps in detail in the following sections.

Quantifying the land use and land cover change
We use the USGS Landsat 5 Surface Reflectance Tier 1 cloudfree images at 30 m spatial resolution taken from Google Earth Engine (Masek et al. 2006). This dataset is produced from Landsat 5 TM L1TP data (level-1 precision and terrain corrected product) and is atmospherically corrected using LEDAPS (Schmidt et al. 2013). We use the scenes of winter months (October-February) to minimize the effects of cloud and seasonal variation. We use the images in adjacent years of 2000 and 2010 to ensure the image quality. Path/row and date information of the used scenes are shown in Table S1.
We apply a geographic object-based image analysis (GEOBIA) classification technique to extract the LULC information from Landsat satellite scenes using eCognition (Trimble Inc., Sunnyvale, CA, USA, 2016) Uddin et al. 2015). The GEOBIA method first identifies adjacent pixels in satellite images with similar spatial and spectral characteristics and partitions the images into distinct segments (Hay and Castilla 2008) and then classifies these segments into different LULC types. We classify Bangladesh into eight LULC types-agricultural land, shrubland, bare land, forest, flowing waterbodies, standing waterbodies, settlement, and mangrove. Our study only focuses on the permanent LULCC. Therefore, we define the standing and flowing waterbodies as areas covered by water throughout the entire year (i.e., permanent water). A detailed description of the LULC classes can be found in Table S2. Our estimated agricultural land area does not subdivide the areas into croplands, plantation and grass/pasturelands, because of the lack of ground observed information.
We use a multi-resolution segmentation method in eCognition. We first identify the scale, shape, and compactness parameters appearing in this method from previous studies Uddin et al. 2015) and then tune and test these parameters to obtain a satisfied segmentation result. We use the following parameters for segmentation: scale, 25; shape, 0.1; and compactness, 0.5. We use highresolution satellite information (such as QuickBird and IKONOS available within Google Earth Pro) to select reference segments (at least 10 for each LULC classes) as training datasets. Overall, 75% of training datasets are randomly selected for developing the classification rules, and the rest are used for validation (Uddin et al. 2015). The classification rules are developed by analyzing the spectral bands, the geometry of the segments, and the spatial contextual parameters. The detailed parameters include TM5 spectral bands, normalized difference vegetation index (NDVI), landwater mask (LWM), soil-adjusted vegetation index (SAVI), brightness, altitude, aspect, and slope (from SRTMGL1 v003 product on Google Earth Engine). The detailed calculation methods of these indices can be found in Gilani et al. (2015).
After the rule-based classification, we use a GIS layer of the rivers and streams (US Geological Survey 2001) to separate flowing and standing waterbodies. This separation only applies to waterbodies and has no effect on rice paddy, which has been classified as a part of agricultural land during the rule-based classification. The waterbodies, which are overlaid with (and close to) the rivers and streams, are classified as flowing waterbodies. Other waterbodies are classified as standing waterbodies. The separation of flowing and standing waterbodies helps to identify the dynamics of lands used for aquaculture, which is an important economic sector in Bangladesh. We consider the impacts from water seasonality by excluding seasonal waterbodies in this study using the Global Surface Water Dataset (Pekel et al. 2016) (see Supplementary Text S1 for detailed procedures and Fig. S1 for the water seasonality maps). Water seasonality discriminates between permanent and seasonal waterbodies, which are defined as waterbodies that are underwater throughout the year and less than 12 months of the year (Pekel et al. 2016). Inclusion of seasonal waterbodies may hinder the driver analysis of waterbody changes. For example, only permanent standing waterbodies can be used for aquaculture; lack of water seasonality information may affect the driver analysis. Therefore, our study focuses only on permanent LULCC. We use the confusion matrix method to assess the classification accuracy (Congalton 1991;Foody 2002) (Supplementary Text  S2 and Tables S3 and S4).

Field observations and Google Earth Pro images
Step 2 Step 3 Step 4 Binary LULCC data in hotspot regions Step 5 Step 1 Step 6

Results
Step 1 Step 2 Step 3 Step 4 Step 5 Step 6 Identifying the hotspot zones The inclusion of low-value regions may dilute the importance of key drivers (Xu et al. 2019). To exclude the regions with smaller areas of LULCC in our driver analysis, we use Hot Spot Analysis (Getis-Ord Gi*) tool in ArcMap 10.4 (ESRI Inc., Redlands CA, USA, 2016) to identify the aggregated areas of major LULCC. The Gi* statistic, which measures zscores of spatially clustered values to identify hotspot regions (Ord and Getis 1995), is performed within the context of neighboring features (upazilas in our case) inside a specified distance band. We use the fixed distance band, which means all neighboring features inside a fixed distance band (Euclidean distance that ensures every feature has at least one neighbor) are all incorporated with equal weight in the computations for the target upazila. Therefore, we do not specify any weight matrix file in this analysis. We consider upazilas with z-scores greater than 3 (99% confidence level) as the hotspot regions.

Proxy data of the biophysical and socioeconomic drivers
Biophysical proxy data Considering the significant impacts of climate change on LULCC in Bangladesh (Chowdhury and Ndiaye 2017), we include seasonal mean temperature and precipitation from 2000 to 2010 as proxy data of the biophysical drivers. We calculate the average values, increasing rate, and standard deviation for temperature and precipitation, for whole year, monsoon months (June-September), and postmonsoon months (October-November) to represent their inter-annual trends based on the CRU NCEP Reanalysis data (Jones and Harris 2014). The increasing rates of precipitation and temperature are calculated using the linear regressions over the years. We use elevation information from the SRTMGL1 v003 product to indicate the topographical features (Yang et al. 2011) and include annual mean values and standard deviation of soil moisture as driver proxies (Le Quéré et al. 2016). A detailed list of proxy data for biophysical drivers is provided in Table S5.

Socioeconomic proxy data
Our socioeconomic database includes 13 variables for 484 sub-districts (  (Sloan 2015). We also include the Euclidean distances to major cities and major highways and other major roads as two socioeconomic proxies based on the spatial data shown in Fig. S2.
Addressing the multi-collinearity of driver proxy data We first calculate the matrix of correlation coefficients for the driver proxy data (Table S6 and S7) to identify a higher level of multi-collinearity among the proxy data for all changes. Then we use the principal component analysis (PCA) method to address the multi-collinearity prior to the logistic regression (Du et al. 2014), because the PCA method generates a set of orthogonal variables from original driver proxies while keeping all driver proxies, which helps us to better interpret the LULCC driver results. We include the components that have ≥ 85% cumulative proportions of variance in the logistic regression (Du et al. 2014). After modeling the relationship between major LULCC and components ("Modeling the drivers using logistic regression" section), we convert the regression coefficients of the PCs back to the original driver proxies (detailed procedures can be found in Supplementary Text S3).

Modeling the drivers using logistic regression
We use the LULCC data in the hotspot zones to generate the binary dependent variable (1 for LULCC activity occurred and 0 for no LULCC occurred) for the LULCC driver analysis (see Supplementary Text S3 for detailed explanation). We keep the original resolutions of the socioeconomic (3 00 km 2 on national average) and biophysical proxy data (0.5°× 0.5°), because downscaling the proxy data to finer resolutions will change the mapping unit of the original data, which may induce the modifiable areal unit problem (MAUP) that causes bias for our logistic regression.
We analyze the drivers of two major LULCC, forest to shrubland and agriculture to standing waterbodies. These activities have relatively large areas and are affected by both biophysical and socioeconomic factors. The 10-year interval conversion data from forest to shrubland also include the changes caused by shifting cultivation (locally referred to as Jhum). This practice involves clearing of forests and shrubs by slashing and burning the crop cultivation, followed by a fallow period (Hossain 2011). The cycling period of shifting cultivation has reduced from 15-25 years in the 1960s to 2-4 years in 1990s (Hossain 2011). A short fallow time does not allow the natural regeneration of forest; when next cultivation starts, the fallow lands are mostly shrubland. Thapa and Rasul (2006) also report an increase of the shrubland area due to the short fallow time. Thus, we treat the change from forest to shrubland over the 10-year interval as "permanent" deforestation activity. We have not analyzed the drivers of the changes between flowing waterbodies, barren land, and other LULC types, because these changes are mainly due to migration of river channel (Bristow 1987) rather than our driver proxies. We have also estimated the drivers of other three important LULCC activities: shrubland to agricultural land, shrubland to forest, and agricultural land to settlement. We use logistic regression to model the drivers of the LULCC and Wald Chi-Square test to evaluate the statistical significance of the logistic regression results. Detailed data processing and model algorithms can be found in Supplementary Text S3.

Synthesis of case studies
We have synthesized existing case studies on the drivers of LULCC in Bangladesh to (1) get a better understanding of the LULCC dynamics, (2) identify potential drivers of LULCC, (3) identify the gaps in the existing studies (e.g., lack of understanding of driving factors at a sub-district level), and (4) collect evidences to complement and evaluate the results of our study (Table S8).  (Table 1 and Fig. S3). The spatial patterns of these LULCs are similar in both years. Agricultural land is widely distributed in both years. Most forests are located in the eastern part of the Chittagong Division (Fig. 2). Flowing waterbodies are mainly consisted of Padma (Ganges), Jamuna (Brahmaputra), and Meghna Rivers, while the standing waterbodies are mainly ponds in the Khulna Division, which is in the southwest of the country.

Land use and land cover change
There are significant exchanges between agricultural land and other land cover classes (e.g., shrubland, barren land, flowing and standing waterbodies, and forest) as shown in Fig. 3 and Table S9. The roughly similar areas of gross gain and loss in agricultural land cause the negligible net change (Fig. 3). Similarly, mangrove also has a small net change. Areas increased from 2000 to 2010 for shrubland by 21.76%, barren land by 41.06%, settlement by 32.52%, and standing waterbodies by 8.63% (Table 1) and decreased for forest by 8.90%, flowing waterbodies by 2.88%. The forest is mainly replaced by shrubland, which has increased by 694 km 2 . The decreasing rate of the forest is the greatest among all LULC types (Table 1), indicating intensive deforestation activities. The net increased rate of barren land of 40.90% is the highest among all LULC types. The sources of increased barren land are mainly from agricultural land and flowing waterbodies ( Fig. 3 and Table S9). The standing waterbody area has increased by 8.58%, because a large area of agricultural land is converted to standing waterbodies. Our hotspot regional analyses show that the conversions from agriculture to standing waterbodies have mainly occurred in the southwest region of the country (Khulna Division) and forest to shrubland conversion is concentrated mainly in Chittagong Division, which is in the southeast region of the country (Fig. 2). The hotspot regions for the changes from shrubland to agricultural land and forest are overlaid in Chittagong Division, while the changes from agricultural land to the settlement have mainly occurred in the northern part of Dhaka City and eastern Chittagong Division (Fig. S4) (Table S12 and S13). We find a similar estimation of LULCC between our data and the ESA CCI data, but greatly different from the GlobeLand30 data. We discuss the comparison results in detail in Supplementary Text S5.

Major drivers
We use the standardized coefficients of driver proxies from the logistic regressions to indicate the impacts of different drivers (Fig. 4). The standardized coefficients refer to how many standard deviations the LULCC area will change per standard deviation increase in the driver proxy data. We have standardized the driver proxy data; therefore, the standardized coefficients allow comparisons of the relative impacts of driver proxies measured on different scales. The relative importance is determined by the absolute mean estimates across coefficients. A positive relationship means that higher values of a driver proxy correspond to greater LULCC areas, whereas a negative relationship refers to higher values of a driver proxy which correspond to smaller LULCC areas. In this section, we only describe the results of the driver analysis. The interpretation and discussion of these results are described in the "Interpretation and explanation of the drivers" section.

Conversion from forest to shrubland
There are seven major biophysical drivers for the conversion from forest to shrubland ( Fig. 4a and Table S14). Increasing rate and standard deviation of temperature are the two most important drivers, and both have positive impacts (the positive relationship between drivers and LULCC) ( Fig. 4a and Table S14). The increasing rate of temperature in monsoon months also has a positive impact on this change. In contrast, four precipitationrelated drivers, mean and standard deviation of precipitation both annually and in monsoon months, are negatively correlated to the changes (the negative relationship between drivers and LULCC). Socioeconomic drivers, including increasing rates of rural and urban household sizes and rural household numbers, are negatively correlated to this change.

Conversion from agricultural land to standing waterbodies
There are eight socioeconomic drivers and two biophysical drivers ( Fig. 4b and Table S15) out of the ten most important drivers for the conversion from agricultural land to standing waterbodies. Population, urban and rural household numbers, and increasing rate of rural household size have negative impacts. The increasing rate of urban household size is positively associated with this change. Meanwhile, the distance to major cities and the connection ratio of electricity (percentage of total households that have access to electricity) are negatively associated with the change, but the distance to major highways and other major roads is positively correlated to the change. The biophysical drivers, distance to rivers, and increasing rate of precipitation in monsoon months all have negative impacts on the changes from agricultural land to standing waterbodies.

Discussion
Interpretation and explanation of the drivers

Conversion from forest to shrubland
Our results show that the forest proportion of Bangladesh in 2010 is 7.37%, which is relatively small compared with other SSEA countries (Xu et al. 2019). Forest and its changes (specifically to shrubland) are mainly located in the Chittagong Division in southeast Bangladesh, which have also been reported in many other studies (Chowdhury and Koike 2010;Hasan et al. 2017;Rasul et al. 2004;Reddy et al. 2016). Biophysical drivers, mainly temperature and precipitation dynamics (Fig. 4a), greatly controlled this change. Increasing but unstable temperatures (higher values of increasing rate and standard deviation of temperature), especially in monsoon months, as well as lower and relatively stable precipitation conditi ons (l ower val ues of m ean preci pi tation and increasing rate of precipitation in monsoon months) facilitate forest loss (Fig. 4a). This type of condition might indicate extreme climate events, such as droughts, which may result in forest loss (Chowdhury and Ndiaye 2017). Meanwhile, lower and stable precipitation conditions may decrease forest Fig. 3 Gross gains, gross losses, and net changes in land use and land cover at a national scale (km 2 ) over the period 2000-2010.
The areas of the gains and losses from one land use/cover type (the y-axis) to the other land use/cover types between 2000 and 2010 are showed in different colors (except black). Black bars demonstrate the net changes (gross gainsgross losses) of different land use/ cover types (the y-axis) productivity  and may lead the forest degradation. Biswas and Choudhury (2007) suggest that socioeconomic factors, such as development of infrastructure and industry, caused deforestation in Bangladesh. However, based on our results we see that the dynamics in the urban and rural household are more important for this type of change. The negative relationships (Fig. 4b) illustrate that the conversion from forest to shrubland, which mainly occurred in regions with slower increasing rates of household size and number. Based on our study, such regions consist mainly of rural regions in the Chittagong Division where the majority of forests are concentrated (Fig. 2). In such areas, local people earn their livelihood by using resources and fuels from forests.

Conversion from agricultural land to standing waterbodies
The results of this study and synthesis analysis demonstrate that the aquaculture industry played an important role in driving the changes in standing waterbodies, which mainly include open waterbodies, such as haors and beels (natural depressions), and closed waterbodies, like baors (oxbow lakes) and ponds. The low-lying areas (haors), (e.g., northeast Sylhet and Mymensingh Divisions) are inundated in rainy seasons. When the dry season comes, water in haors recedes, leaving only more depressed lands under water, called beels (Banglapedia 2015). Our analysis suggests that the standing waterbodies in the northeast region are mainly seasonal water due to flooding during monsoon seasons. We used the Global Surface Water Dataset (Pekel et al. 2016) to exclude the seasonal standing waterbodies. As a result, the permanent standing waterbodies are mainly observed in a very limited area, particularly in Khulna Division in the southwest of the country (Fig. 2), in forms of baors and ponds, which are mainly used for aquaculture (Islam 2010). The most identifiable limitation of the global surface water data is the geographic and temporal discontinuities of the Landsat archive data (Pekel et al. 2016). There are some regions or periods lacking satellite observations, which may affect the discrimination of the seasonal and permanent water (Fig. S1). However, our analysis only accounts for permanent water. Therefore, the effect of missing data on our analysis results is minimum, because the regions or periods with no observations have been treated as seasonal water as well. All permanent water according to the global surface water data has continuous Landsat observations across the whole year.
The development of aquaculture is happening at the cost of agricultural land. This phenomenon is also observed by other satellite studies such as Abdullah et al. (2019). Our driver analysis indicates that the changes from agricultural land to standing waterbodies occur mainly in less populous rural regions with smaller household size, which are close to major Fig. 4 Ten most prominent drivers for the land use and land cover changes of a forest to shrub land and b agricultural land to standing waterbody from 2000 to 2010. Blue dots show the standardized coefficients, error bars show the 95% confidence interval. Standardized coefficients demonstrate the relative importance of each driver cities. Compared with urban areas, these regions have relatively less access to infrastructures, such as electricity and major highways and roads (note that these regions still have enough infrastructure accesses than remote rural areas that enable preservation and transportation of aquaculture products). These regions are the most suitable areas for aquaculture development. On one hand, in urban regions, there are limited lands for aquaculture development and land expenses are higher. On the other hand, in remote rural areas far away from major cities, the preservation and transportation costs of aquaculture products are higher. In these less populous rural regions close to major cities, the relatively lower land expenses and preservation and transportation costs increase the economic profitability of aquaculture. This will lead more local families to convert their agricultural lands to standing waterbodies (Ahmed and Diana 2015;Ali 2006;Islam and Tabeta 2016). The regions close to major cities have higher domestic demands for aquaculture products. Meanwhile, the international demands for aquaculture products also significantly increased in past few years: the shrimp export (particularly to the USA and European Union) became the second largest export industry in Bangladesh in recent years (Ahmed and Diana 2015), and its value increased from 240 million US dollars in 2000-2001to 459 million in 2015(Fisheries Resources Survey System 2016. Both domestic and international demands have incentivized farmers to change their agricultural land to standing waterbodies for aquaculture products. Water availability directly affects aquaculture. The regions close to rivers (smaller value of distance to rivers) have more available water, which can facilitate the aquaculture ponds development. Both of our LULCC analyses (Fig. 2) and case studies (Abdullah et al. 2019;Islam 2010) suggest that shrimp farming is mainly located in the water-rich coastal regions of southwest Bangladesh. Changes in climate, such as the increasing rate of precipitation in monsoon months in our results, have negative impacts on agricultural production in the southwest of the country Huq et al. 2015), particularly the rice cultivation (Ali 2006). Literature research suggests that farmers might opt for aquaculture to avoid possible risks to agricultural production due to climate change, because Bangladesh has suitable agro-climatic conditions, adequate water resources, and cheap labor force (Paul and Vogl 2011) and the profit of aquaculture is higher than cultivation (e.g., shrimp farming is about 12 times more profitable than rice cultivation in Bangladesh (Ali 2006)).
We have also studied the drivers of the changes from shrubland to agricultural land, shrubland to forest, and agricultural land to settlement. The changes from shrubland to agricultural land are mainly driven by biophysical factors such as increasing rate of precipitation and mean soil moisture. Socioeconomic drivers, such as increasing rate of rural household size and increasing rate of population, are dominant in driving the change from shrubland to forest. For the changes from agricultural land to settlement, both socioeconomic (e.g., distance to major highways and other major roads and urban household number) and biophysical drivers (e.g., distance to rivers and elevation) are important (Fig. S5). We have discussed the drivers of these three LULCC activities in detail in Supplementary Text S6.

Limitation and caveats
While our analysis of spatial drivers of LULCC for Bangladesh is important and our analysis results are corroborated with other studies, there are some limitations related to the data and methods used. First, our analysis for the dynamics of LULCC captures only the decadal changes and can mask within-decade variations including intermediary land uses. However, wall-to-wall analysis of Landsat scenes at much finer temporal resolution is laborious and beyond the scope of this study. Second, our Landsat analysis detects changes in land type only when the magnitude of modification is large enough to cause a shift from one land cover category to another (e.g., forest to shrubland). However, in reality, forest loss and regrowth are gradual and cause subtle modifications to land cover. Nevertheless, our analysis bias is likely to be minimal because our statistical estimation weighs each observation (at grid cell) by the magnitude of land change; thus, small changes have less influence in our model analysis. Third, quality of the socioeconomic data in some regions may be poor due to misreporting, human errors in computerization, quality of village/town boundaries, or unavailability of data due to separation and union of different upazilas boundaries. Yet this is the most consistent and detailed data we can access. Finally, our analysis does not extend beyond 2010 due to lack of socioeconomic data. However, we will extend our land cover conversion estimates and its driver as soon as the data are available.

Conclusions
Through a comprehensive methodology, this study successfully reveals the dominant biophysical and socioeconomic drivers of major LULCC activities in Bangladesh. Climate dynamics and extremes are critical to the conversion from forest to shrubland, while population and accessibility to infrastructure are controlling factors of the conversion from agricultural land to standing waterbodies. This study is unique, because it advances a general understanding of the dynamics and drivers of LULCC at the national scale of Bangladesh, which will guide effective national-level planning and policymaking. Particularly, we have accounted for the impacts of water seasonality on LULCC analysis; such effort has not been reported in previous studies. Our study adds an additional quantitative dimension by providing the explanation of drivers of LULCC. In addition, the synthesis of case studies provides a more generalized understanding of the drivers of LULCC and reinforced the findings of our spatial determinants of LULCC in Bangladesh over the time period of 2000-2010.
Data availability Our satellite LULC data for Bangladesh for years 2000 and 2010 can be downloaded for free from the ISAM website (http:// climate.atmos.uiuc.edu/Bangladesh/). We are sharing the data at 30 m spatial resolution. Our geospatial sub-district level socioeconomic database (covering 2001 and 2011) is also available at the same website.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.