Next Article in Journal
The Impact of Mosaic Land Use and Land Cover on the Quality of River Waters (Case Study: Lubelskie Province, E Poland)
Next Article in Special Issue
Review of Soil Salinity and Sodicity Challenges to Crop Production in the Lowland Irrigated Areas of Ethiopia and Its Management Strategies
Previous Article in Journal
Detectability of the Critically Endangered Araucaria angustifolia Tree Using Worldview-2 Images, Google Earth Engine and UAV-LiDAR
Previous Article in Special Issue
Integrating Climate Change Adaptation and Mitigation into Land Use Optimization: A Case Study in Huailai County, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimizing Land Use and Land Cover Allocation for Flood Mitigation Using Land Use Change and Hydrological Models with Goal Programming, Chaiyaphum, Thailand

by
Athiwat Phinyoyang
and
Suwit Ongsomwang
*
School of Geoinformatics, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
*
Author to whom correspondence should be addressed.
Land 2021, 10(12), 1317; https://doi.org/10.3390/land10121317
Submission received: 19 October 2021 / Revised: 26 November 2021 / Accepted: 27 November 2021 / Published: 30 November 2021
(This article belongs to the Special Issue Land Use Optimisation)

Abstract

:
Floods represent one of the most severe natural disasters threatening the development of human society worldwide, including in Thailand. In recent decades, Chaiyaphum province has experienced a problem with flooding almost every year. In particular, the flood in 2010 caused property damage of 495 million Baht, more than 322,000 persons were affected, and approximately 1046.4 km2 of productive agricultural area was affected. Therefore, this study examined how to optimize land use and land cover allocation for flood mitigation using land use change and hydrological models with optimization methods. This research aimed to allocate land use and land cover (LULC) to minimize the surface for flood mitigation in Mueang Chaiyaphum district, Chaiyaphum province, Thailand. The research methodology consisted of six stages: data collection and preparation, LULC classification, LULC prediction, surface runoff estimation, the optimization of LULC allocation for flood mitigation and mapping, and economic and ecosystem service value evaluation and change. According to the results of the optimization and mapping of suitable LULC allocation to minimize surface runoff for flood mitigation in dry, normal, and wet years using goal programming and the CLUE-S model, the suitable LULC allocation for flood mitigation in 2049 under a normal year could provide the highest future economic value and gain. In the meantime, the suitable LULC allocation for flood mitigation in 2049 under a drought year could provide the highest ecosystem service value and gain. Nevertheless, considering future economic and ecosystem service values and changes with surface runoff reduction, the most suitable LULC allocation for flood mitigation is a normal year. Consequently, it can be concluded that the derived results of this study can be used as primary information for flood mitigation project implementation. Additionally, the presented conceptual framework and research workflows can be used as a guideline for government agencies to examine other flood-prone areas for flood mitigation in Thailand.

1. Introduction

Floods represent one of the most severe natural disasters threatening the development of human society worldwide, including in Thailand. They cause enormous losses to economies, societies, and ecological environments [1], and the flood-related damage to agriculture and other related activities impacts a country’s economy and development [2].
In general, the primary cause of flooding is heavy rainfall [3]. However, many other causes are also due to human activities, such as land degradation; deforestation of catchment areas, urban growth, and increased population along riverbanks [4,5,6]; poor land use planning, zoning, and control of flood plain development; poor drainage, particularly in cities; and insufficient management of discharge from river reservoirs [7].
In the last two decades, Chaiyaphum province has experienced a problem with flooding almost every year, causing a loss of lives, as well as economic losses, asset or housing losses, inundated farmlands, and decreased crop productivity for people who live in this area. In particular, the flood in 2010 caused property damage of 495 million Baht. More than 322,000 persons were affected, at least seven persons lost their lives, and approximately 1046.4 km2 of productive agricultural area was affected [8].
Due to the risk of large-scale damage to public and private property in Chaiyaphum province, the Royal Thai Government has allocated a significant budget to mitigate flood effects using structural measures, such as channel modification, bank protection, dikes, and reservoir development. However, the problems persist and are becoming exacerbated [9]. However, it is difficult to fundamentally mitigate flood damage using only flood prevention facilities [10]. A comprehensive flood control measure should consider land use and land cover change and optimum land use allocation.
In general, LULC strongly influences flood risk and affects the probability of floods and their consequences in several ways [11]. LULC change can affect the hydrological characteristics of a river basin through the influence of land use on runoff generation processes [12,13]. This study chose the SCS-CN method, which represents a distributed hydrologic model, to estimate the time series surface runoff according to LULC changes in the study period (2001–2019). These changes may alter the quantity of surface/subsurface runoff generation, river flooding regimes, and the extent [14]. Thus, defining optimal strategies for appropriate flood management, especially LULC management, is very important and necessary [15] for flood mitigation in Chaiyaphum province.
Land use optimization is one of the proper solutions for soil and water conservation at the watershed level. It can help decision-makers determine the best scenario of various land use alternatives without sacrificing the economic value obtained from the available land use [16,17]. Land use arrangement can be optimized using a programming model to increase land use earnings and reduce environmental impacts, especially surface runoff [16]. The essence of management science, manifested in modeling and programming techniques, is considered an essential tool for optimally allocating rare resources to gain the most benefits [18].
In recent decades, the new programming methods that have been developed can be employed under conflicting conditions of the goals and limited resources for decision-makers. In natural resource management, there are many optimization techniques. Some approaches such as linear programming (LP), goal programming (GP), and weighted goal programming (WGP) are widely employed in land use optimization at the watershed level [11]. For instance, Yeo et al. applied LP to optimize land use to peak discharge minimization at the Old Woman Cheek watershed, Ohio State, USA [19]. Owji et al. applied LP for land use optimization in the Jajrood watershed, Iran, to reduce surface runoff and sediment yield [20]. Likewise, Aldea et al. used GP for forest management in the Pinar Grande Forest, Spain [21]. Further, Gonfa and Kumar applied LP and GP for optimum land use to minimize soil erosion and maximize the net benefit in Ethiopia’s Mojo watershed [22], and Al-Zahrani et al. developed GP for optimizing water resources in Riyadh, Saudi Arabia [23]. Similarly, Sokouti and Nikkami applied LP to optimize land use patterns to reduce soil erosion in West Azerbaijan province, Iran [24]. WGP has been applied to optimize LULC allocation for surface runoff and sediment load minimization at Bayg watershed [11]. Moreover, LP has been used to maximize cropland allocation in Abaro Kebele, Ethiopia [25]. Recently, Han et al. applied LP to optimize the land use structure for carbon emission reduction in Shenzhen, China [26].
Nevertheless, the integration of the optimization technique (GP), advanced land use change modeling (CLUE-S model), and the distributed hydrological model (SCS-CN model) to minimize surface runoff for flood mitigation does not exist in Thailand. Therefore, a novel classification method, random forests, was first applied to classify LULC data in 2001, 2010, and 2019. Then, the classified data were further used to predict a time series LULC between 2001 and 2019 using the CLUE-S model for time series surface runoff estimation using the SCS-CN model. After this, goal programming was applied to minimize surface runoff for flood mitigation based on the surface runoff coefficient value of each LULC type in 2029, 2039, and 2049 in dry, normal, and wet years. Finally, economic and ecosystem service value change between the existing LULC data in 2019 and the suitable LULC allocation in dry, normal, and wet years was examined in terms of gain and loss using the present value (PV) model and the simple benefit transfer method.
The specific objectives of the study were (1) to classify LULC data in 2001, 2010, and 2019 using the random forest classifier, (2) to predict LULC change in two periods (2002–2009 and 2011–2018) using the CLLUE-S model, (3) to estimate surface runoff between 2001 and 2019, (4) to optimize and map LULC allocation for flood mitigation under three rainfall conditions, and (5) to evaluate economic and ecosystem service values and change for the most suitable LULC allocation for flood mitigation.

2. Study Area

The study area was the Second Part of the Lam Nam Chi watershed, Chaiyaphum province, Thailand, under the Chi River basin, covering approximately 3794 km2. As mentioned earlier, the selected study area covers the flood-prone area in Mueang Chaiyaphum district, Chaiyaphum province. The topography of the area is generally characterized by rolling hilly terrain and flat areas. The elevation ranges from 162 m above the mean sea level (MSL) in the lower part of the watershed to approximately 1034 m above MSL in the upper part of the watershed (Figure 1). The study area consists of nine soil groups: clay, clay loam, loam, loamy sand, sandy loam, sandy clay loam, silty clay, silty loam, and silty clay loam [27]. Meanwhile, the top three dominant land use types in 2015 were paddy fields (43.47%), cassava (12.69%), and forest land (12.48%) [28].

3. Materials and Methods

The research methodology consisted of data collection and preparation and five significant components, which were (1) LULC classification, (2) LULC prediction, (3) surface runoff estimation, (4) optimization of LULC allocation for flood mitigation and mapping, and (5) economic and ecosystem service value evaluation and change (Figure 2). Details of each stage are separately described in the following sections.

3.1. Data Collection and Preparation

The required input data for data analysis included GIS data, remote sensing data, and primary and secondary data, which were collected and prepared, as summarized in Table 1.

3.2. LULC Classification

Landsat imageries in 2001, 2010, and 2019 were downloaded from the USGS website (www.earthexplore.usgs.gov, accessed on 22 November 2021) for LULC classification using the RF classifier of the EnMap-Box software. In practice, the training areas of each LULC type in a specific year were separately prepared to extract multiple decision trees for LULC classification. Spectral reflectance data (visible, NIR, and SWIR bands), additional spectral bands, and elevation were applied to classify the LULC types. The spectral bands that enhance particular features for LULC classification include the Normalized Difference Vegetation Index (NDVI) to represent vegetation features [29], the Modified Normalized Difference Wetness Index (MNDWI) to signify the moisture regime [30], and the Normalized Difference Built-up Index (NDBI) to indicate built-up areas [31]. Likewise, elevation is directly related to the spatial distribution of LULC type, e.g., paddy fields are generally situated in the floodplain, while forests are primarily located in mountainous areas.
In this study, the modified land use classification of the LDD consisted of (1) urban and built-up areas, (2) paddy fields, (3) sugarcane, (4) cassava, (5) other field crops, (6) para rubber, (7) perennial trees and orchards, (8) forest land, (9) waterbodies, (10) rangeland, (11) marshes and swamps, and (12) unused land.
After classification, the LULC maps in 2001, 2010, and 2019 were assessed for thematic accuracy (overall accuracy and Kappa hat coefficient) based on the reference data from color orthophotograph in 2000–2001, very high spatial resolution imageries from Google Earth in 2010, and field surveys in 2020, respectively. This study estimated the number of sample sizes for thematic accuracy assessment based on multinomial distribution with a stratified random sampling scheme, as suggested by [32].

3.3. LULC Prediction

The CLUE-S model was selected to predict LULC data in two periods, 2002–2009 and 2011–2018, for filling the gap of LULC data between three classified LULC data in 2001, 2010, and 2019. As a result, time-series LULC data between 2001 and 2019 will be available for annual surface runoff estimation in this study.

3.3.1. Optimal Local Driving Factors on Land Use Change Identification for LULC Prediction

The land use change model CLUE-S was selected to predict the LULC data in two periods, 2002–2009 and 2011–2018, for data analysis. The local driving factors on land use change for LULC prediction were identified by comparing the predicted LULC map in 2019 with the classified LULC map in 2019. The basic parameters of the CLUE-S model, which include (1) elasticity value, (2) LULC conversion matrix, and (3) land requirement of each LULC type in 2019, were firstly extracted based on the final LULC map in 2001 and 2010 using the Markov Chain model. At the same time, the selected three driving factor categories on LULC change, including physical, socioeconomic, and proximity data, which were reviewed from the previous studies of many researchers [33,34,35,36,37,38,39,40,41,42,43,44,45], were examined multicollinearity, and significant driving factors for LULC identified by allocating using binomial logistic regression analysis, as follows:
Log ( P i 1 P i ) = β 0 + β 1 X 1 , i + β 2 X 2 , i + β n X n , i
where Pi is the probability of a grid cell for the considered land use type on location i, and the Xs are the location factors. The coefficients (β) were estimated through logistic regression using the actual land use pattern as the dependent variable [46].
After this, the predicted LULC map in 2019 was compared to the classified map in 2019 using a wall-to-wall thematic accuracy assessment with overall accuracy and Kappa hat coefficient. If the overall accuracy and Kappa hat coefficient were equal to or more than 80%, then the derived significant driving factors by binomial logistic regression analysis were chosen as the optimal local driving factors on land use change for LULC prediction using the CLUE-S model.

3.3.2. LULC Prediction of Two Periods: 2002–2009 and 2011–2018

The optimal local driving factors on land use change for LULC prediction, namely, elasticity value, LULC conversion matrix, and land requirement of each LULC type in two time periods (2002–2009 and 2011–2018), which were extracted using the Markov Chain model based on the corresponding LULC data in 2001, 2010, and 2019, were applied to predict the LULC data in the two periods using the CLUE-S model.

3.4. Surface Runoff Estimation

This study estimated the time series surface runoff between 2001 and 2019 based on the classified and predicted LULC data, soil series, and rainfall data using the SCS-CN method with suitable AMC via ESRI ArcGIS software.
In practice, the required input data included LULC, soil series, rainfall, and hydrologic soil group data, prepared and operated for surface runoff estimation using the SCS-CN method in raster format with a cell size of 30 m in a raster-based GIS environment. The surface runoff depth in each cell was semi-automatically generated based on runoff curve numbers (CNs) according to hydrologic soil group–land cover complex using the Model Builder of ArcGIS, as shown in Figure 3.
The soil series data, which were applied to classify soil texture classes based on the percentage of sand, silt, and clay of each soil series unit, are displayed in Figure 4a. Meanwhile, the hydrologic soil group (HSG), which presents potential runoff in the study area according to soil texture classes, is shown in Figure 4b.
In this study, two significant steps implemented under this component include (1) suitable AMC for surface runoff estimation using the SCS-CN method and (2) surface runoff estimation between 2011 and 2019.

3.4.1. Suitable AMC for Surface Runoff Estimation Using the SCS-CN Method

The suitable AMC condition for surface runoff estimation using the SCS-CN method was examined based on the classified and predicted LULC data between 2001 and 2010 with three different CN values of three different AMCs, as suggested by [47] using the following equations:
CN I = 4 . 2 CN II 10     0 . 058 CN II
CN III = 23 CN II 10 + 0 . 13 CN II
where CN I is the runoff curve number value of each LULC type of AMC-I, CN II is the runoff curve number value of each LULC type of AMC-II, and CN III is the runoff curve number value of each LULC type of AMC-III (Table 2).
In practice, the three runoff curve numbers of the hydrologic soil group of each LULC type were separately applied to estimate the potential maximum storage under three different AMCs using Equation (4).
S = 25.4 1000 CN   10
where CN is the runoff curve number of the hydrologic soil group (HSG)–land cover complex.
The calculated potential maximum storage was further applied to estimate the surface runoff depth of the three different AMCs [48,49], as follows:
Q = ( P     0.2 S ) 2 ( P + 0.8 S )
where Q is the surface runoff depth (mm), P is the annual rainfall (mm), and S is the potential maximum storage.
Then, the estimated surface runoff depth of the three different AMCs from 2001 to 2010 was converted into surface runoff volume using Equation (6).
Surface   runoff   volume = Surface   runoff   depth 1000   ×   cell   size
Later, they were used to identify the suitable AMC using model performance scale, including Nash and Sutcliffe’s coefficient of efficiency (NSE), the coefficient of determination (R2), and the percent of bias (PBIAS) (Equations (7)–(9)), as suggested by [50] (Table 3).
NSE = 1 [ i n ( Q simi Q obsi ) 2 i n ( Q obsi Q avg ) 2 ]
where n is the number of years, Q simi is the simulated surface runoff, Q obsi is the observed surface runoff, and Q avg is the average observed surface runoff over the simulation period. The values for E can vary from –∞ to 1, with 1 indicating a perfect fit.
R 2 = { i = 1 n ( Q obsi Q obsavg ) ( Q simi Q simavg ) [ i = 1 n ( Q obsi Q obsav g ) 2 i = 1 n ( Q simi Q simavg ) 2 ] 0 . 5 } 2
where Q obsi is the observed surface runoff in year i , Q simi is the simulated surface runoff in year i , Q obsavg is the average of the observed surface runoff over the calibration or validation period, Q simavg is the average of the simulated surface runoff over the validation period, i is the year, and n is the total count of data pairs.
PBIAS = [ i = 1 n ( Y i obs Y i sim ) × ( 100 ) i = 1 n Y i obs ]
where Y i obs is the observed surface runoff in time step i , and Y i sim is the simulated surface runoff in time step i .
In this study, the observed runoff data between 2001 and 2010 from the hydrological station at E.21, E.23, and E.6C of the RID were used to calculate an average NSE, R2, and PBIAS for suitable AMC identification (see the location of the station in Figure 1).

3.4.2. Surface Runoff Estimation between 2011 and 2019

The runoff curve number of the suitable AMC was applied to estimate the surface runoff between 2011 and 2019 based on the classified and predicted LULC data in the same period. The estimated surface runoff data between 2011 and 2019 were further examined for model validation based on the observed runoff data in the same period from the same gauges of the RID using NSE, R2, and PBIAS.
The time series surface runoff estimation (between 2001 and 2019) from two steps will be further applied to extract the average runoff coefficient of each LULC type of three rainfall conditions for optimizing LULC allocation for flood mitigation in the next component.

3.5. Optimization of LULC Allocation for Flood Mitigation and Mapping

Goal programming was first applied to optimize LULC allocation for flood mitigation under dry, normal, and wet years. Then, the suitable LULC allocation data of dry, normal, and wet years were mapped using the CLUE-S model.

3.5.1. SPI Calculation for the Rainfall Condition Identification

Annual rainfall data between 1987 and 2019 were first used to calculate the 12-month SPI values. Then, they were reclassified into three rainfall conditions according to the SPI drought classification of [51] as follows:
  • If any year had a 12-month SPI value less than or equal to –0.50, the annual rainfall was categorized as a dry year;
  • If any year had a 12-month SPI value between –0.49 and 0.49, the annual rainfall was categorized as a normal year;
  • If any year had a 12-month SPI value more than or equal to 0.50, the annual rainfall was categorized as a wet year.
After that, the average surface runoff coefficient for each LULC type of three rainfall conditions (dry, normal, and wet years) was extracted from the time series surface runoff and LULC data using zonal statistical analysis for optimization of LULC allocation to minimize surface runoff for flood mitigation.

3.5.2. Optimization of LULC Allocation to Minimize Surface Runoff for Flood Mitigation

This study assigned the constraint sets for optimizing LULC allocation in 2029, 2039, and 2049 using the Markov Chain model based on the historical LULC development between 2010 and 2019. The changing area of each LULC type was considered according to the derived transitional change area from the Markov Chain model. Then, the derived average runoff coefficient and constraints of the objective function were applied to optimize LULC allocation to minimize surface runoff for flood mitigation under dry, normal, and wet years using goal programming with “What’s Best!” as an extension program in an MS Excel environment.
The goal programming model, working as the surface runoff minimization function, can be expressed as the following equations:
Minimize surface runoff:
Min ( z ) = i = 1 n C i X i
Subject to constraints:
i = 1 n X i = A
i = 1 n X i   A
i = 1 n X i   A
X i     0
where Z is the total annual surface runoff of the study area (m3/year), C i is the average surface runoff coefficient in each land use type (m3/km2/year), X i is the area of land use class i (km2), n is the number of land use classes, and A is the total area of land use classes (km2).
Finally, the suitable LULC allocation to minimize surface runoff for flood mitigation under dry, normal, and wet years was separately identified based on surface runoff reduction compared to the actual data in 2019.

3.5.3. Mapping of Suitable LULC Allocation for Flood Mitigation

The suitable LULC allocation to minimize surface runoff for flood mitigation under dry, normal, and wet years was separately mapped using the CLUE-S model. The required input data for LULC mapping included LULC data in 2019, local driving factors on land use change, elasticity value, LULC conversion matrix, and suitable LULC allocation data of three rainfall conditions as the land requirement.

3.6. Economic and Ecosystem Service Value Evaluation and Change

Economic and ecosystem service value evaluation and change were separately evaluated using the present value (PV) model and the simple benefit transfer method based on LULC data in 2019 and suitable LULC allocation data to minimize surface runoff for flood mitigation under dry, normal, and wet years in terms of gain and loss.

3.6.1. Economic Value Evaluation and Change

For economic value evaluation, the LULC data in 2019 and the suitable LULC allocation data to minimize surface runoff under dry, normal, and wet years were first calculated in terms of economic values using Equation (15), as suggested by [52]. Then, the future economic value changes between the LULC data in 2019 and the suitable LULC allocation data to minimize surface runoff under dry, normal, and wet years were compared using the image algebra change detection algorithm for gain and loss.
PV = FV [ 100 100 + IR ] Y
where PV is the present value, FV is the future value, IR is the interest rate in percent, and Y is the number of years from the present, counting from zero.

3.6.2. Ecosystem Service Value Evaluation and Change

In general, ecosystem services represent a dynamic field in current scientific research, linking ecological, economic, and social aspects, demanding practical applications and methodologies at different spatial scales, and maintaining environmental management and decision-making processes [53,54,55].
This study applied a simple benefit transfer method [53] to evaluate the ecosystem service values of LULC data in 2019 and suitable LULC allocation data to minimize surface runoff under dry, normal, and wet years using Equation (16) with the coefficient value for different LULC types (Table 4). Then, the ecosystem service value change between the LULC data in 2019 and the suitable LULC allocation data to minimize surface runoff under dry, normal, and wet years was compared using the image algebra change detection algorithm for gain and loss.
ESV =   ( A k   ×   VC k )
where ESV denotes the total value of the ecosystem service, while A k and VC k represent the area and value coefficient for proxy LULC type k , respectively.
Finally, the future economic and ecosystem service values and their changes in suitable LULC allocation for flood mitigation by each rainfall condition (dry, normal, and wet years) were compared to justify the most suitable LULC allocation for flood mitigation among three rainfall conditions.

4. Results and Discussion

4.1. Land Use and Land Cover Classification and Change Detection

The results of the LULC classification in 2001, 2010, and 2019 using the random forests (RF) classifier based on Landsat imageries and supplementary data, including NDVI, MNDWI, NDBI, and DEM, are reported in Table 5 and displayed in Figure 5.
The results reveal that in 2001, the top three most dominant LULC types were paddy fields (61.79%), forest land (16.66%), and cassava (14.05%). On the contrary, the top three least dominant LULC types were other field crops (0.06%), marshes and swamps (0.30%), and para rubber (0.44%), which are randomly distributed in the study area. Meanwhile, the top three most dominant in 2019 were paddy fields (54.58%), cassava (16.59%), and forest land (15.94%). Conversely, the top three least dominant LULC types in 2001 were other field crops (0.14%), marshes and swamps (0.79%), and para rubber (0.88%). Recently, the top three most dominant LULC types in 2019 were paddy fields (53.03%), cassava (12.91%), and forest land (12.68%). In contrast, the top three least dominant LULC types were other field crops (0.16%), marshes and swamps (0.73%), and waterbodies (1.40%). The temporal change in LULC type coverage in 2001, 2010, and 2019 is presented in Figure 6.
Furthermore, a simple comparison of the LULC type coverage with the changing area, annual change rate, and percentage of in two short periods (2001–2010 and 2010–2019) and in the long period (2001–2019) are reported in Table 6, Table 7 and Table 8. The annual change rate of the short- and long-term periods plays a significant role in land requirement estimation for LULC prediction using the CLUE-S model.
Moreover, by considering the performance of the RF classifier for LULC classification, the overall accuracy and Kappa hat coefficient of the classified LULC maps in 2001 were 89.88% and 84.88%, in 2010 were 90.71% and 87.03%, and in 2019 were 91.37% and 88.26%. The overall accuracy of more than 85% of the three maps provides acceptable results, as suggested by [57]. Likewise, the Kappa hat coefficient of the three maps was more than 80%, representing strong agreement or accuracy between the classified map and the ground reference information [58]. The derived overall accuracy and Kappa hat coefficient values of the classified LULC maps are comparable to those of other researchers who classified LULC maps based on Landsat data with the RF classifier [44,59,60,61,62,63,64,65]

4.2. Driving Factor Identification for LULC Change

The results of the binary logistics regression analysis for identifying LULC type location preference after the multicollinearity test are reported in Table 9. The most common vital driving factor for all LULC type changes was elevation. Meanwhile, the most important driving factors for field crops (sugarcane, cassava, and other field crops) included elevation, slope, and annual rainfall. The specific driving factors for each LULC type preference from binary logistics regression are further applied to allocate LULC type for predicting LULC change under the CLUE-S model.
Moreover, the derived area under the curve (AUC) values for each LULC type allocation varied from 0.6185 to 0.9824. This finding suggests a fair-to-excellent fit between the predicted and real LULC transition [66,67,68].

4.3. Local Parameter for LULC Prediction Using the CLUE-S Model

After a wall-to-wall accuracy assessment, by comparison, for the predicted LULC in 2019 with the classified LULC in 2019, the overall accuracy was 86.95%, and the Kappa hat coefficient was 80.72%. Thus, the significant driving factors by binomial logistic regression analysis were chosen as the optimal local driving factors on land use change for LULC prediction using the CLUE-S model in this study (see Table 9).
In the meantime, the optimal conversion matrix for each LULC type possibly changed between 2010 and 2019 (Table 10) and the land use type resistance (elasticity) values based on the transition probability matrix of LULC changed between 2010 and 2019 from the Markov Chain model (Table 11), identified as a local parameter for LULC prediction. Thus, the elasticity value of the urban and built-up areas, paddy fields, sugarcane, cassava, other field crops, para rubber, perennial trees and orchards, forest land, waterbodies, rangeland, marshes and swamps, and unused land were 1.00, 0.93, 0.93, 0.65, 0.99, 0.80, 0.99, 0.80, 0.91, 0.89, 0.39, and 0.96, respectively.

4.4. LULC Prediction between 2002 and 2009

LULC prediction data between 2002 and 2009, which were simultaneously allocated based on the conversion matrix of LULC change (Table 10), elasticity value (Table 11), annual land demand (Table 12), and driving factors on LULC change for specific LULC type location preference (Table 9), are reported in Table 13 and displayed in Figure 7.
As a result, the LULC prediction between 2002 and 2009 was dictated by the historical LULC development between 2001 and 2010. The results indicate that the most increasing LULC types were cassava and sugarcane, with an increasing annual change rate of 10.71 and 10.24 km2 per year, while the paddy field was the most decreasing LULC type with a decreasing annual change rate of 30.42 km2 per year.

4.5. LULC Prediction between 2011 and 2018

The derived optimum local parameter of the CLUE-S model with annual land demand between 2010 and 2019 based on the annual change rate of the Markov Chain model (Table 14) was simultaneously combined with driving factors for LULC change for a specific LULC type location preference (Table 9) to predict LULC data, as shown by the results in Table 15 and Figure 8.
As with the previous results, the LULC prediction between 2011 and 2018 was enforced by the historical LULC development between 2010 and 2019. Again, the results revealed that sugarcane was the most increasing LULC type, with an increasing annual rate of 17.04 km2 per year, while cassava and forest land were the most decreasing LULC types, with decreasing annual change rates of 15.49 and 13.71 km2 per year.
In summary, the proportional area of the time series LULC data between 2001 and 2019 based on RF classification and CLUE-S prediction is displayed in Figure 9.

4.6. Surface Runoff Estimation between 2001 and 2010

The surface runoff estimation between 2001 and 2010 was estimated with three runoff curve number (CN) values assigned based on three different AMC conditions (see Table 3) for suitable AMC identification in the study area. In principle, the AMC is an indicator of watershed wetness and availability of soil storage. The annual surface runoff volume of the three AMC conditions with rainfall data between 2001 and 2010 is presented in Table 16 and Figure 10.
According to the results, the accumulated surface runoff volume under three AMC conditions was relatively different, because the CN values of these conditions were different (Table 3), while annual rainfall data were the same. The patterns of surface runoff volume under the three AMC conditions and annual rainfall between 2001 and 2010 were similar. The higher the annual rainfall, the higher the surface runoff. This finding was confirmed by simple linear regression analysis, as shown in Figure 11. The surface runoff volume under the three AMC conditions positively correlated with annual rainfall, with R2 values of 0.8503, 0.8511, and 0.8513, respectively. These coefficient values show a strong relationship between annual rainfall and surface runoff, according to [50].
In addition, the results of model performance using NSE, R2, and PBIAS for identifying the suitable AMC for surface runoff estimation based on observed and simulated data between 2001 and 2010 (Table 17) are reported in Table 18.
As a result, in Table 18, under the AMC-I condition, the NSE value varied from 0.67 to 0.95, while the R2 value varied from 0.71 to 0.96 (Figure 12). The PBIAS values ranged from –1.63% for overestimation to 2.62% for underestimation. Meanwhile, under the AMC-II condition, the NSE value varied from 0.86 to 0.94, while the R2 value varied from 0.96 to 0.97 (Figure 13). The PBIAS values ranged from –4.58% for overestimation to 1.21% for underestimation. In the meantime, under the AMC-III condition, the NSE value varied from 0.32 to 0.96, while the R2 value varies from 0.96 to 0.97 (Figure 14). The PBIAS values ranged from –8.63% to –0.11% for overestimation.
These findings indicate that the NSE values under the three AMCs conditions can provide a perfect fit for surface runoff estimation. Likewise, the R2 values indicate a very high correlation between the observed and estimated surface runoff for all AMC conditions. Similarly, the PBIAS values show a relatively different overestimation and underestimation among three AMC conditions less than ±10%. These results show a very good fit between the observed and estimated surface runoff for the surface runoff estimation using the SCS-CN method, as suggested by [50] (see Table 3).
Furthermore, the NSE, R2, and PBIAS of all of the AMCs were compared to identify the suitable AMC for surface runoff estimation between 2011 and 2019, as shown in Table 19. As a result, AMC-II provided better overall average statistics measurements than the other AMCs. Thus, the AMC-II condition was chosen as the suitable AMC for surface runoff estimation in the second period (2011–2019).

4.7. Surface Runoff Estimation between 2011 and 2019

The LULC and annual rainfall data between 2011 and 2019 as dynamic input with the hydrologic soil group were applied to estimate the surface runoff with the suitable AMC-II using the Model Builder of ArcGIS (See Figure 3). The annual surface runoff volume with rainfall data between 2011 and 2019 is reported in Table 20 and Figure 15.
As a result, the surface runoff volume ranged from approximately 1004 million m3 in 2014 to approximately 6142 million m3 in 2011. The pattern of surface runoff volume and annual rainfall between 2011 and 2019 was similar. As mentioned in the previous section, the surface runoff volume positively correlated with annual rainfall, with an R2 of 0.8511. This finding is consistent with previous studies [69,70].
Furthermore, the results of model performance using NSE, R2, and PBIAS for validating surface runoff estimation based on observed and simulated data between 2011 and 2019 (Table 21) are reported in Table 22.
As a result, in Table 22, the NSE value from three stations varied from 0.87 to 0.91, with an average of 0.8933. The R2 value varied from 0.87 to 0.94, with an average of 0.9033 (Figure 16). The PBIAS values ranged from –5.71% for overestimating to 3.66% for underestimating.
According to a statistical report of model performance for surface runoff estimation between 2011 and 2019, the derived NSE and R2 values were more than 0.87 and the PBIAS value was less than ±10. These results show a very good fit for the surface runoff estimation with a very high correlation between the observed and estimated surface runoff, as suggested by [50] (see Table 3). Thus, it can be concluded that the surface runoff estimation using the SCS-CN method in the current study can be validated with acceptable results.

4.8. SPI Calculation for Rainfall Condition Identification

The available historical rainfall data records from 33 years (1987–2019) of the Chaiyaphum meteorological station (Figure 17) were applied to calculate the 12-month SPI index, as shown in Figure 18. The derived SPI values were classified into seven drought types and grouped into three rainfall conditions: dry, normal, and wet years, as summarized in Table 23.
As a result, the dry years were found to be 2001, 2003, 2004, 2014, 2015, 2018, and 2019. Meanwhile, normal years occurred in 2002, 2005, 2006, 2012, and 2016, and wet years in 2007, 2008, 2009, 2010, 2011, 2013, and 2017.

4.9. Optimization of LULC Allocation for Flood Mitigation

Goal programming of multi-objective decision analysis (MODA) was applied to allocate the optimum LULC to minimize surface runoff for flood mitigation based on the average surface runoff coefficient from LULC under three rainfall conditions. In this study, an average of each LULC type was separately extracted from the time series surface runoff between 2001 and 2019 for dry, normal, and wet years, as presented in Table 24.
At the same time, the constraints of goal programming subjected to the change of LULC area in the study area were assigned for optimizing the LULC allocation for flood mitigation based on the historical LULC development between 2010 and 2019. In this study, a 10-year period was chosen to predict LULC data in 2029, 2039, and 2049 based on a period of input data (2010–2019) for transitional area prediction by the Markov Chain model. Accordingly, the changing area of each LULC type was categorized into two groups: Decreased area (paddy fields, cassava, forest land, waterbodies, rangeland, and marshes and swamps) and increased area (urban and built-up areas, sugarcane, other field crops, para rubber, perennial trees and orchards, and unused land) according to the derived transitional change area from the Markov Chain model. Details of the LULC change to minimize surface runoff in 2029, 2039, and 2049 as constraints of goal programming are summarized in Table 25.
Based on the linearity of the objective function and constraints, the objective functions of the surface runoff minimization problem for optimizing LULC allocation in dry, normal, and wet years were formulated as shown in Equations (17)–(19).
Min ( Z ) = 0 . 45 X 1 + 0 . 44 X 2 + 0 . 38 X 3 + 0 . 45 X 4 + 0 . 39 X 5 + 0 . 37 X 6 + 0 . 38 X 7 + 0 . 28 X 8 + 0 . 49 X 9 + 0.42 X 10 + 0.54 X 11 + 0.46 X 12
Min ( Z ) = 1.03 X 1 + 1.01 X 2 + 0.90 X 3 + 1.00 X 4 + 0.88 X 5 + 0.83 X 6 + 0.84 X 7 + 0.76 X 8 + 1.04 X 9 + 0.92 X 10 + 1.14 X 11 + 1.02 X 12
Min ( Z ) = 1.46 X 1 + 1.44 X 2 + 1.30 X 3 + 1.41 X 4 + 1.31 X 5 + 1.30 X 6 + 1.31 X 7 + 1.17 X 8 + 1.51 X 9 + 1.36 X 10 + 1.59 X 11 + 1.41 X 12
The constraints for optimizing LULC allocation in 2029, 2039, and 2049 to minimize runoff under three rainfall conditions using objective functions (Equations (17)–(19)) are summarized in Table 26. Additionally, two mandatory constraints, as described by Equations (11) and (14) in Section 3.5.2, were set under the three rainfall conditions as follows:
The first mandatory constraint, related to the area of all land use classes, must be equal to the allowable area of 3794.22 km2, as shown in Equation (20).
X 1 + X 2 + X 3 + X 4 + X 5 + X 6 + X 7 + X 8 + X 9 + X 10 + X 11 + X 12 = 3794 . 22
The second mandatory constraint is related to the non-negative variable. The area of each land use class should be more than or equal to 0 km2, as shown in Equation (21).
X 1 + X 2 + X 3 + X 4 + X 5 + X 6 + X 7 + X 8 + X 9 + X 10 + X 11 + X 12     0
The objective functions for dry, normal, and wet years (Equations (17)–(19)) were then transformed into a goal programming form as follows:
f k ( x ) + d k   d k   + = a k
The results of the optimized LULC allocation to minimize surface runoff for flood mitigation in 2029, 2039, and 2049 under the three rainfall conditions (dry, normal, and wet years) are presented in Table 27, Table 28 and Table 29.
In dry years, the significant increasing LULC type after optimizing LULC allocation in 2029, 2039, and 2049 were sugarcane, with approximately 80, 148, and 210 km2. In contrast, the significantly decreasing LULC type was cassava, with approximately 82, 133, and 200 km2. Additionally, the annual surface runoff decreased in 2029, 2039, and 2049 by 12.95, 20.78, and 27.89 million m3, or approximately 0.82%, 1.32%, and 1.77% of total estimated surface runoff in 2019 (See Table 27).
Meanwhile, the most increased LULC area under the normal year was also sugarcane, with an area of approximately 79, 145, and 205 km2. In contrast, the significant decreasing LULC types were paddy fields, with approximately 65, 133, and 200 km2, and sugarcane, with an area of approximately 82, 130, and 160 km2. Additionally, the annual surface runoff decreased in 2029, 2039, and 2049 by 21.34, 37.59, and 52.31 million m3, or 0.59%, 1.03%, and 1.43% of the total estimated surface runoff in 2019 (see Table 28).
In the meantime, the most increasing LULC type in wet years was sugarcane, with an area of approximately 115, 213, and 293 km2, respectively. Conversely, the significant decreasing LULC types were paddy fields, with approximately 65, 133, and 200 km2, and sugarcane, with approximately 82, 131, and 160 km2. Additionally, the annual surface runoff decreased in 2029, 2039, and 2049 by 18.52, 33.06, and 46.73 million m3, or approximately 0.35%, 0.63%, and 0.89% from the total estimated surface runoff in 2019 (see Table 29).
According to these results, the optimized LULC allocation data in 2049 of the dry, normal and wet years were suitable for flood mitigation according to surface runoff reduction. The surface runoff under dry, normal and wet years in 2049 reduced by approximately 27.89, 52.31, and 46.73 million m3 (Figure 19).
Moreover, notable reductions in annual surface runoff under the three rainfall conditions occurred in paddy fields and cassava to sugarcane, para rubber, and perennial trees and orchards after optimization, which changed the hydrological properties. Paddy and cassava fields provided higher runoff coefficients than sugarcane, rubber, and perennial trees and orchards (see Table 27, Table 28 and Table 29). This finding is consistent with previous studies [19,20,71].
The deviation in annual surface runoff after minimization by goal programming is presented in Table 30. As a result, the percentage deviation from the goal varied from –1.77% to –0.35%, which is acceptable.

4.10. Mapping of LULC Allocation for Flood Mitigation

LULC data in 2019 and the derived optimum local parameter of the CLUE-S model were applied to map the suitable LULC allocation to minimize surface runoff for flood mitigation in 2049 of the three rainfall conditions. The conversion matrix for each possible LULC type change in 2049 was based on transitional LULC change between 2010 and 2019, as shown in Table 31. Meanwhile, the elasticity values were assigned based on the transition probability matrix of LULC change between 2019 and 2049 from the Markov Chain model, as shown in Table 32. Thus, the elasticity values of the urban and built-up areas, paddy fields, sugarcane, cassava, other field crops, para rubber, perennial trees and orchards, forest land, waterbodies, rangeland, marshes and swamps, and unused land for LULC prediction in 2049 were 1.00, 0.83, 0.83, 0.31, 0.97, 0.52, 0.96, 0.70, 0.76, 0.71, 0.76, and 0.87, respectively.
In addition, the optimized LULC allocation in 2049 in dry, normal, and wet years (see Table 27, Table 28 and Table 29) was applied to calculate the annual land requirement between 2019 and 2049 for LULC prediction.
The spatial distribution of the LULC allocation maps in 2049 for dry, normal, and wet years is displayed in Figure 20. Meanwhile, the area and percentage of LULC type in 2049 for dry, normal, and wet years are summarized in Table 33. The derived map of suitable LULC allocation to minimize surface runoff for flood mitigation of dry, normal, and wet years was further applied to evaluate economic and ecosystem service values and change.

4.11. Economic Value Evaluation and Change

The actual LULC data in 2019 and suitable LULC data in 2049 for flood mitigation in dry, normal, and wet years (Table 34) were applied to evaluate future economic values using the PV model (Table 35), and the results are shown in Table 36 and Figure 21.
According to results, the top three dominant LULC types of actual LULC data in 2019 were paddy fields, cassava, and forest land, while the top three dominant LULC types of suitable LULC allocation data in 2049 for flood mitigation in dry, normal, and wet years were paddy field, sugarcane, and forest land. Moreover, the top three highest future values were forest land, perennial trees and orchards, and para rubber (Table 35).
As a result, in Table 36, the suitable LULC allocation data for flood mitigation in normal years provided the highest economic value, approximately 149,211 million Baht. In contrast, the actual LULC allocation in 2019 provided the lowest economic value, approximately 144,068 million Baht. The contribution of the economic value by LULC type from the actual LULC and suitable LULC allocation data for flood mitigation from three rainfall conditions are displayed in Figure 22.
According to the data in Figure 22, the top three dominant LULC types of actual LULC in 2019, including paddy fields, cassava, and forest land, will provide future economic value in 2049, about 89% of the total value. Meanwhile, the top three dominant LULC types of suitable LULC allocation for flood mitigation under dry, normal, and wet years, including paddy fields, sugarcane, and forest land, deliver future economic value in 2049 of approximately 85%, 86%, and 88% of the total value, respectively.
Moreover, the future economic value of forest land from actual LULC and suitable LULC allocation data for flood mitigation (dry, normal, and wet years) contributed the highest values compared to the other LULC types because the present economic value, approximately 25,000,000 Baht per km2, or the future economic value, approximately 165,359,154 Baht per km2, was exceptionally high when compared to other LULC types (see Table 35).
Furthermore, the results of the future economic value change by comparing the values of actual LULC data and each suitable LULC allocation for flood mitigation (dry, normal, and wet years) in terms of gain (+sign) and loss (-sign) are reported in Table 37, Table 38 and Table 39 and spatially displayed in Figure 23.
According to change detection, suitable LULC allocation for flood mitigation in 2049 in normal years gained the highest future economic value of approximately 4322 million Baht (Table 38). On the contrary, suitable LULC allocation for flood mitigation in 2049 in wet years gained the highest future economic value of approximately 3124 million Baht (Table 39). These results show the consequence of LULC allocation for flood mitigation in 2049 using goal programming on future economic value, because the future economic value depends on the areas of LULC type and their values.

4.12. Ecosystem Service Value Evaluation and Change

The actual LULC data in 2019 and the suitable LULC data in 2049 for flood mitigation under dry, normal, and wet years (Table 40) were applied to evaluate the ecosystem service values using a simple benefit transfer method according to a coefficient value of LULC type (see Table 4), and the results are shown in Table 41 and Figure 24.
As a result, the ESV of actual LULC in 2019 and the suitable LULC allocation for flood mitigation in 2049 in dry, normal, and wet years are slightly different. The suitable LULC allocation for flood mitigation in 2049 under the dry year provided the highest ESV, approximately 15,481 million Baht, while the suitable LULC allocation for flood mitigation in 2049 in wet years delivered the lowest ESV, approximately 15,237 million Baht. The contributions of ESV of each LULC type from actual LULC in 2019 and the suitable LULC allocation for flood mitigation in 2049 under the three rainfall conditions are compared in Figure 25.
As shown in Figure 25, paddy fields provided the highest ESV value from the actual LULC and the suitable LULC allocation data for flood mitigation in 2049 under the three rainfall conditions, because paddy fields were the most dominant type in the area. Additionally, the top three dominant LULC types, namely, paddy fields, field crops, and forest land, delivered ESVs of approximately 86%, 88%, 89%, and 88% of the total value.
Furthermore, the results of the ESV change upon comparison between values of actual LULC data and each suitable LULC allocation for flood mitigation (dry, normal, and wet years) in terms of gain (+sign) and loss (-sign) are reported in Table 42, Table 43 and Table 44 and spatially displayed in Figure 26.
According to the results in Table 42, Table 43 and Table 44, the suitable LULC allocation for flood mitigation in 2049 under dry years gained the highest ESV of approximately 207 million Baht. On the contrary, the suitable LULC allocation for flood mitigation in 2049 in wet years resulted in losses in ESV of approximately 38 million Baht. These findings indicate that the ESV of suitable LULC allocation for flood mitigation in 2049 in wet years is lower than the actual LULC in 2019 (see Table 45). Similar to the change in future economic value, these results show the consequence of LULC allocation for flood mitigation in 2049 using goal programming on ESV, because the ESV depends on areas of LULC type and their coefficient values. Moreover, it can be observed that forest land, waterbodies, and marshes and swamps provide a gain in ESV (+sign), while urban and built-up areas, paddy fields, field crops, rangeland, and unused land results in a loss of ESV (-sign) for all three rainfall conditions. These findings indicate that the ecosystem service value was dictated by the coefficient value of each LULC type.
In summary, the suitable LULC allocation for flood mitigation in 2049 in normal years provides the highest value for future economic value evaluation and the highest gain value compared to actual LULC in 2019. In the meantime, the suitable LULC allocation for flood mitigation in 2049 in dry years provided the highest value for ecosystem service evaluation and the highest gain the ESV by comparing it with actual LULC in 2019. Meanwhile, the suitable LULC allocation for flood mitigation in 2049 in normal years can reduce the highest surface runoff by approximately 52 million m3 compared to the actual LULC in 2019 (see Table 45).
Consequently, it can be concluded that the most suitable LULC allocation for flood mitigation in 2049 in Chaiyaphum district, Chaiyaphum province under the Second Part of the Lam Nam Chi watershed, based on the future economic value and ecosystem service value evaluation, is a suitable LULC allocation for flood mitigation in 2049 in the normal year scenario. This information can be used as primary data for supporting project implementation.

5. Conclusions

This study applied the supervised method to classify LULC data in 2001, 2010, and 2019 based on Landsat 5-TM and Landsat 8-OLI with supplementary data, including NDVI, MNDWI, NDBI, and DEM, using an RF classifier under EnMap BOX software. The derived thematic accuracy of the LULC maps showed an overall accuracy and Kappa hat coefficient between classified LULC maps and ground reference data in 2001, 2010, and 2019 of 89.88% and 84.88%, 90.71% and 87.03%, and 91.37% and 88.26%, respectively. Later, the classified LULC data in 2001, 2010, and 2019 were further applied to predict the LULC change in two periods, 2002–2009 and 2011–2018, using the CLUE-S model. The significant driving factors of LULC change for specific LULC type location preferences included elevation, slope, annual rainfall, average income per capita at the sub-district level, population density at the sub-district level, distance to the road network, distance to a stream, and distance to the existing urban area. As a result, the LULC prediction of both periods was dictated by the historical LULC development between 2001 and 2010 and 2010 and 2019, respectively. Then, time series surface runoff data between 2001 and 2019 were estimated using the SCS-CN method under a GIS raster-based environment. The process worked on spatial variation of land use, hydrologic soil group, and rainfall data. In this study, a suitable AMC condition was first examined and validated for time series surface runoff estimation between 2001 and 2010. Then, a suitable AMC condition was further chosen to estimate the time series surface runoff between 2011 and 2019.
After this, goal programming was applied to minimize the surface runoff for flood mitigation based on the surface runoff coefficient value of LULC types in dry, normal, and wet years for 2029, 2039, and 2049. Accordingly, the surface runoff could be reduced under all three rainfall conditions, and suitable LULC allocation for flood mitigation in dry, normal, and wet years was in 2049. The suitable LULC allocation for flood mitigation in 2049 of the normal year provided the highest value and gain for future economic value compared to the actual LULC in 2019. Meanwhile, the suitable LULC allocation for flood mitigation in 2049 of dry years provided the highest value and gain for ecosystem services compared to the actual LULC in 2019. Nonetheless, considering the future economic and ecosystem service values and changes in surface runoff reduction, the most suitable LULC allocation for flood mitigation in 2049 was normal years.
In conclusion, the derived results of this study can be used as primary information for flood mitigation project implementation in Chaiyaphum province. Likewise, the presented conceptual framework and research workflows can be used as a guideline for government agencies to examine flood-prone areas for flood mitigation in Thailand.
However, to apply the proposed method in other areas, we recommend that the CN value of the AMC-II condition, as the identified suitable value in this study, can directly apply to estimate time series surface runoff. In addition, rainfall conditions identification using SPI can be ignored to increase the number of years for calculating the average runoff coefficient value of each LULC type. This value plays a vital role in minimizing surface runoff for flood mitigation using goal programming.

Author Contributions

Conceptualization, S.O. and A.P.; methodology, S.O. and A.P.; software, A.P.; validation, S.O. and A.P.; formal analysis, S.O. and A.P.; investigation, S.O. and A.P.; data curation, A.P.; writing—original draft preparation, A.P.; writing—review and editing, S.O.; visualization, A.P.; supervision, S.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Suranaree University of Technology for supporting the facilities to undertake this research. Additionally, the special thanks from the authors go to anonymous reviewers for their valuable comments and suggestions that improve our manuscript from various perspectives.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, D.; Xie, P.; Dong, X.; Hu, B.; Ji, L.; Li, Y.; Peng, T.; Ma, H.; Wang, K.; Xu, S. Improvement of the SWAT model for event-based flood simulation on a sub-daily timescale. Hydrol. Earth Syst. Sci. Discuss. 2018, 22, 5001–5019. [Google Scholar] [CrossRef] [Green Version]
  2. Jothityangkoon, C.; Maskong, H.; Sangthong, P.; Kosa, P. Development processes of a master plan for flood protection and mitigation in a community area: A case study of Roi Et province. KKU Eng. J. 2015, 42, 287–291. [Google Scholar]
  3. Kuntiyawichai, K.; Sri-Amporn, W.; Wongsasri, S.; Chindaprasirt, P. Anticipating of potential climate and land use change impacts on floods: A case study of the lower Nam Phong river basin. Water 2020, 12, 1158. [Google Scholar] [CrossRef] [Green Version]
  4. Mbow, C.; Diop, A.; Diaw, A.T.; Niang, C.I. Urban sprawl development and flooding at Yeumbeul suburb (Dakar-Senegal). Afr. J. Environ. Sci. Technol. 2008, 2, 75–88. [Google Scholar]
  5. Prasad, A.S.; Pandey, B.W.; Leimgruber, W.; Kunwar, R.M. Mountain hazard susceptibility and livelihood security in the upper catchment area of the river Beas, Kullu Valley, Himachal Pradesh, India. Geoenvironmental Disasters 2016, 3, 1–17. [Google Scholar] [CrossRef] [Green Version]
  6. Billi, P.; Alemu, Y.T.; Ciampalini, R. Increased frequency of flash floods in Dire Dawa, Ethiopia: Change in rainfall intensity or human impact? Nat. Hazards 2015, 76, 1373–1394. [Google Scholar] [CrossRef]
  7. Danumah, J.H.; Odai, S.N.; Saley, B.M.; Szarzynski, J.; Thiel, M.; Kwaku, A.; Kouame, F.K.; Akpa, L.Y. Flood risk assessment and mapping in Abidjan district using multi-criteria analysis (AHP) model and geoinformation techniques, (cote d’ivoire). Geoenvironmental Disasters 2016, 3, 1–13. [Google Scholar] [CrossRef] [Green Version]
  8. Department of Disaster Prevention and Mitigation. Report of Damage from a Flooding Situation; Ministry of Interior: Bangkok, Thailand, 2019.
  9. Sriwongsitanon, N. Flood forecasting system development for the upper Ping River basin. Kasetsart J. Nat. Sci. 2010, 44, 717–731. [Google Scholar]
  10. Banba, M. Influences of regional development on land use of Nagara Basin and flood risk control. In Proceedings of the 3rd European Conference on Flood Risk Management, Lyon, France, 17–21 October 2016; pp. 1–7. [Google Scholar]
  11. Tajbakhsh, S.M.; Memarian, H.; Kheyrkhah, A. A GIS-based integrative approach for land use optimization in a semi-arid watershed. Glob. J. Environ. Sci. Manag. 2018, 4, 31–46. [Google Scholar]
  12. Garg, V.; Nikam, B.R.; Thakur, P.K.; Aggarwal, S.P.; Gupta, P.K.; Srivastav, S.K. Human-induced land use land cover change and its impact on hydrology. HydroResearch 2019, 1, 48–56. [Google Scholar] [CrossRef]
  13. Leta, M.K.; Demissie, T.A.; Tränckner, J. Hydrological responses of watershed to historical and future land use land cover change dynamics of Nashe watershed, Ethiopia. Water 2021, 13, 2372. [Google Scholar] [CrossRef]
  14. Kuntiyawichai, K. Interactions between Land Use and Flood Management in the Chi River Basin; Wageningen University: Wageningen, The Netherlands, 2012. [Google Scholar]
  15. Tingsanchali, T.; Karim, F. Flood-hazard assessment and risk-based zoning of a tropical flood plain: Case study of the Yom River, Thailand. Hydrol. Sci. J. 2010, 55, 145–161. [Google Scholar] [CrossRef]
  16. Riedel, C. Optimizing land use planning for mountainous regions using LP and GIS towards sustainability. India. J. Soil Conserv. 2003, 34, 121–124. [Google Scholar]
  17. Sadeghi, S.H.R.; Jalili, K.; Nikkami, D. Land use optimization in watershed scale. Land Use Policy 2009, 26, 186–193. [Google Scholar] [CrossRef]
  18. Nikkami, D.; Elektorowicz, M.; Mehuys, G.R. Optimizing the management of soil erosion. Water Pollut. Res. J. Can. 2002, 37, 577–586. [Google Scholar] [CrossRef] [Green Version]
  19. Yeo, I.-Y.; Gordon, S.I.; Guldmann, J.-M. Optimizing patterns of land use to reduce peak runoff flow and nonpoint source pollution with an integrated hydrological and land use model. Earth Interact. 2004, 8, 1–20. [Google Scholar] [CrossRef]
  20. Owji, M.R.; Nikkami, D.; Mahdian, M.H.; Mahmoudi, S. Minimizing surface runoff by optimizing land use management. World Appl. Sci. J. 2012, 20, 170–176. [Google Scholar]
  21. Aldea, J.; Martínez-Peña, F.; Romero, C.; Diaz-Balteiro, L. Participatory goal programming in forest management: An application integrating several ecosystem services. Forests 2014, 5, 3352–3371. [Google Scholar] [CrossRef] [Green Version]
  22. Gonfa, Z.B.; Kumar, D. Optimal land use planning in mojo watershed with multi-objective linear programming. Am. Int. J. Res. Hum. Arts Soc. Sci. 2015, 13, 10–17. [Google Scholar]
  23. Al-Zahrani, M.; Musa, A.; Chowdhury, S. Multi-objective optimization model for water resource management: A case study for Riyadh, Saudi Arabia. Environ. Dev. Sustain. 2016, 18, 777–798. [Google Scholar] [CrossRef]
  24. Sokouti, R.; Nikkami, D. Optimizing land use pattern to reduce soil erosion. Eurasian J. Soil Sci. 2017, 6, 75–83. [Google Scholar] [CrossRef]
  25. Mellaku, M.T.; Reynolds, T.W.; Woldeamanuel, T. Linear programming-based cropland allocation to enhance performance of smallholder crop production: A pilot study in Abaro Kebele, Ethiopia. Resources 2018, 7, 76. [Google Scholar] [CrossRef] [Green Version]
  26. Han, D.; Qiao, R.; Ma, X. Optimization of land-use structure based on the trade-off between carbon emission targets and economic development in Shenzhen, China. Sustainability 2019, 11, 11. [Google Scholar] [CrossRef] [Green Version]
  27. Land Development Department. Soil Series of Thailand; Ministry of Agriculture and Cooperatives: Bangkok, Thailand, 2019.
  28. Land Development Department. Land use data of Thailand; Land Development Department, Ministry of Agriculture and Cooperatives: Bangkok, Thailand, 2016.
  29. Rouse, J.; Haas, R.H.; Schell, J.A.; Deering, D. Monitoring vegetation systems in the great plains with ERTS. In Proceedings of the Third Earth Resources Technology Satellite-1 Symposium- Volume I: Technical Presentations. NASA SP-351, Greenbelt, MD, USA, 1 January 1973; pp. 3010–3017. [Google Scholar]
  30. Xu, H. A new index for delineating built-up land features in satellite imagery. Int. J. Remote Sens. 2008, 29, 4269–4276. [Google Scholar] [CrossRef]
  31. Zha, Y.; Gao, J.; Ni, S. Use of normalized difference built-up index in automatically mapping urban areas from TM imagery. Int. J. Remote Sens. 2003, 24, 583–594. [Google Scholar] [CrossRef]
  32. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data: Principles and Practices, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  33. Trisurat, Y.; Alkemade, R.; Verburg, P.H. Projecting land-use change and its consequences for biodiversity in northern Thailand. Environ. Manag. 2010, 45, 626–639. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Han, H.; Yang, C.; Song, J. Scenario simulation and the prediction of land use and land cover change in Beijing, China. Sustainability 2015, 7, 4260–4279. [Google Scholar] [CrossRef] [Green Version]
  35. Ongsomwang, S.; Iamchuen, N. Integration of geospatial models for optimum land use allocation in three different scenarios. Suranaree J. Sci. Technol. 2015, 22, 377–396. [Google Scholar]
  36. Zheng, H.W.; Shen, G.Q.; Wang, H.; Hong, J. Simulating land use change in urban renewal areas: A case study in Hong Kong. Habitat Int. 2015, 46, 23–34. [Google Scholar] [CrossRef] [Green Version]
  37. Xu, X.; Du, Z.; Zhang, H. Integrating the system dynamic and cellular automata models to predict land use and land cover change. Int. J. Appl. Earth Obs. Geoinf. 2016, 52, 568–579. [Google Scholar] [CrossRef]
  38. Gao, C.; Zhou, P.; Jia, P.; Liu, Z.; Wei, L.; Tian, H. Spatial driving forces of dominant land use/land cover transformations in the Dongjiang River watershed, Southern China. Environ. Monit. Assess. 2016, 188, 84. [Google Scholar] [CrossRef]
  39. Ongsomwang, S.; Boonchoo, K. Integration of geospatial models for the allocation of deforestation hotspots and forest protection units. Suranaree J. Sci. Technol. 2016, 23, 283–307. [Google Scholar]
  40. Li, X.; Wang, Y.; Li, J.; Lei, B. Physical and socioeconomic driving forces of land-use and land-cover changes: A case study of Wuhan City, China. Discrete Dyn. Nat. Soc. 2016, 2016, 1–11. [Google Scholar] [CrossRef] [Green Version]
  41. Phompila, C.; Lewis, M.; Ostendorf, B.; Clarke, K. Forest cover changes in Lao tropical forests: Physical and socio-economic factors are the most important drivers. Land 2017, 6, 23. [Google Scholar] [CrossRef]
  42. Arowolo, A.O.; Deng, X. Land use/land cover change and statistical modelling of cultivated land change drivers in Nigeria. Reg. Environ. Change 2018, 18, 247–259. [Google Scholar] [CrossRef]
  43. Palchowdhuri, Y.; Roy, P.S. Driver based statistical model for simulating land use/land cover change in Indus river basin, India. Remote Sens. Land 2018, 2, 15–30. [Google Scholar] [CrossRef] [Green Version]
  44. Ongsomwang, S.; Pattanakiat, S.; Srisuwan, A. Impact of land use and land cover change on ecosystem service values: A case study of Khon Kaen City, Thailand. Environ. Nat. Resour. J. 2019, 17, 43–58. [Google Scholar] [CrossRef] [Green Version]
  45. Nguyen, H.H.; Dargusch, P.; Moss, P.; Aziz, A.A. Land-use change and socio-ecological drivers of wetland conversion in Ha Tien Plain, Mekong Delta, Vietnam. Land Use Policy 2017, 64, 101–113. [Google Scholar] [CrossRef]
  46. Verburg, P.H.; Lesschen, J.-P. Practical: Explorative Modeling of Future Land Use for the Randstad Region of the Netherlands; Wageningen University: Wageningen, The Netherlands, 2014. [Google Scholar]
  47. Chow, V.T.; Maidment, D.R.; Mays, L.W. Applied Hydrology; McGraw-Hill: New York, NY, USA, 1988. [Google Scholar]
  48. United States Department of Agriculture. Urban Hydrology for Small Watersheds, Title 210-VI-TR-55, 2nd ed.; Conservation Engineering Division, National Resources Conservation Service, United States Department of Agriculture: Washington, DC, USA, 1986.
  49. Weng, Q. Remote Sensing and GIS Integration Theories, Methods, and Applications; McGraw-Hill: New York, NY, USA, 2010. [Google Scholar]
  50. Me, W.; Abell, J.M.; Hamilton, D.P. Effects of hydrologic conditions on SWAT model performance and parameter sensitivity for a small, mixed land use catchment in New Zealand. Hydrol. Earth Syst. Sci. 2015, 19, 4127–4147. [Google Scholar] [CrossRef] [Green Version]
  51. Liu, Y.; Zhou, Y.; Ju, W.; Wang, S.; Wu, X.; He, M.; Zhu, G. Impacts of droughts on carbon sequestration by China’ s terrestrial ecosystems from 2000 to 2011. Biogeosciences 2014, 11, 2583–2599. [Google Scholar] [CrossRef] [Green Version]
  52. Rossiter, D.G. Lecture Notes: “Land Evaluation”, Part 4: Economic Land Evaluation; College of Agriculture & Life Sciences, Department of Soil, Crop, & Atmospheric Sciences: Cornell University: Ithaca, NY, USA, 1994. [Google Scholar]
  53. Costanza, R.; d’Arge, R.; de Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, R.V.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  54. Millennium Ecosystem Assessment. In Ecosystems and Human Well-Being: Synthesis; Island Press: Washington, DC, USA, 2005.
  55. TEEB. The Economics of Ecosystems and Biodiversity Ecological and Economic Foundations; Kumar, P., Ed.; Earthscan Publications: London, UK, 2010. [Google Scholar]
  56. Mamat, A.; Halik, Ü.; Rouzi, A. Variations of ecosystem service value in response to land-use change in the Kashgar Region, Northwest China. Sustainability 2018, 10, 200. [Google Scholar] [CrossRef] [Green Version]
  57. Anderson, J.R.; Hardy, E.E.; Roach, J.T.; Witmer, R.E. A Land Use and Land Cover Classification System for Use with Remote Sensor Data; 964; United States Government Printing Office: Washington, DC, USA, 1976; pp. 1–28. [Google Scholar]
  58. Fitzpatrick-Lins, K. Comparison of sampling procedures and data analysis for a land-use and land-cover map. Photogramm. Eng. Remote Sens. 1981, 47, 343–351. [Google Scholar]
  59. Na, X.; Zhang, S.; Li, X.; Yu, H.; Liu, C. Improved land cover mapping using random forests combined with landsat thematic mapper imagery and ancillary geographic data. Photogramm. Eng. Remote Sens. 2010, 76, 833–840. [Google Scholar] [CrossRef]
  60. Rodriguez-Galiano, V.F.; Ghimire, B.; Rogan, J.; Chica-Olmo, M.; Rigol-Sanchez, J.P. An assessment of the effectiveness of a random forest classifier for land-cover classification. ISPRS J. Photogramm. Remote Sens. 2012, 67, 93–104. [Google Scholar] [CrossRef]
  61. Gartzia, M.; Alados, C.L.; Pe’rez-Cabello, F.; Bueno, C.G. Improving the accuracy of vegetation classifications in mountainous areas: A case study in the Spanish Central Pyrenees. Mt. Res. Dev. 2013, 33, 63–74. [Google Scholar] [CrossRef] [Green Version]
  62. Jhonnerie, R.; Siregar, V.P.; Nababan, B.; Prasetyo, L.B.; Wouthuyzen, S. Random Forest Classification for mangrove land cover mapping using Landsat 5 TM and Alos Palsar imageries. Procedia Environ. Sci. 2015, 24, 215–221. [Google Scholar] [CrossRef] [Green Version]
  63. Eisavi, V.; Homayouni, S.; Yazdi, A.M.; Alimohammadi, A. Land cover mapping based on random forest classification of multitemporal spectral and thermal images. Environ. Monit. Assess. 2015, 187, 291. [Google Scholar] [CrossRef] [PubMed]
  64. Kulkarni, A.; Lowe, B. Random Forest algorithm for land cover classification. Int. J. Recent Innov. Trends Comput. Commun. 2016, 4, 58–63. [Google Scholar]
  65. Pareeth, S.; Karimi, P.; Shafiei, M.; De Fraiture, C. Mapping agricultural landuse patterns from time series of Landsat 8 using random forest based hierarchical approach. Remote Sens. 2019, 11, 601. [Google Scholar] [CrossRef] [Green Version]
  66. Pontius, R.G.; Schneider, L.C. Land-cover change model validation by an ROC method for the Ipswich watershed, Massachusetts, USA. Agric. Ecosyst. Environ. 2001, 85, 239–248. [Google Scholar] [CrossRef]
  67. Vilar del Hoyo, L.; Martín Isabel, M.P.; Martínez Vega, F.J. Logistic regression models for human-caused wildfire risk estimation: Analysing the effect of the spatial accuracy in fire occurrence data. Eur. J. For. Res. 2011, 130, 983–996. [Google Scholar] [CrossRef]
  68. Liang, X.; Liu, X.; Chen, G.; Leng, J.; Wen, Y.; Chen, G. Coupling fuzzy clustering and cellular automata based on local maxima of development potential to model urban emergence and expansion in economic development zones. Int. J. Geogr. Inf. Syst. 2020, 34, 1930–1952. [Google Scholar] [CrossRef]
  69. Canqiang, Z.; Wenhua, L.; Biao, Z.; Moucheng, L. Water yield of Xitiaoxi River basin based on InVEST modeling. J. Resour. Ecol. 2012, 3, 50–54. [Google Scholar] [CrossRef]
  70. Kasei, R.A.; Ampadu, B.; Sapanbil, G.S. Relationship between rainfall-runoff on the White Volta River at Pwalugu of the Volta Basin in Ghana. Environ. Earth Sci. 2013, 3, 92–99. [Google Scholar]
  71. Hargreaves, J.J.; Hobbs, B.F. Optimal selection of priority development areas considering tradeoffs between hydrology and development configuration. Environ. Model. Assess. 2009, 14, 289–302. [Google Scholar] [CrossRef]
  72. Office of Agricultural Economics. Agricultural Statistics of Thailand 2019. Available online: https://www.oae.go.th/assets/portals/1/files/jounal/2563/yearbook62edit.pdf (accessed on 19 April 2021).
  73. Office of the Cane and Sugar Board. Annual Report of Sugarcane Production in 2019. Available online: http://www.ocsb.go.th/upload/journal/fileupload/923-1854.pdf (accessed on 19 April 2021).
  74. Rubber Authority of Thailand. Thailand Rubber Price in 2019. Available online: https://www.raot.co.th/rubber2012/rubberprice_yr.php (accessed on 19 April 2021).
  75. Wittawatchutikul, P.; Jirasuktaveekul, W. The Technique for Court Witness: A Case of Compensation Claim due to Deforestation, Technical Report No 9/2548; National Park, Wildlife and Plant Conservation Department: Bangkok, Thailand, 2005.
  76. Bank for Agriculture and Agricultural Cooperatives. Minimum Retail Rate in 2019. Available online: https://www.baac.or.th/th/content-rate.php?content_group=9&content_group_sub=2&inside=1 (accessed on 19 April 2021).
Figure 1. Terrain characteristics of the study area with runoff gauge stations.
Figure 1. Terrain characteristics of the study area with runoff gauge stations.
Land 10 01317 g001
Figure 2. Workflow of the research methodology.
Figure 2. Workflow of the research methodology.
Land 10 01317 g002
Figure 3. Schematic diagram of Model Builder for surface runoff estimation.
Figure 3. Schematic diagram of Model Builder for surface runoff estimation.
Land 10 01317 g003
Figure 4. Spatial distribution of (a) soil series and (b) hydrologic soil group.
Figure 4. Spatial distribution of (a) soil series and (b) hydrologic soil group.
Land 10 01317 g004
Figure 5. Spatial distribution of the LULC classification in: (a) 2001, (b) 2010, and (c) 2019.
Figure 5. Spatial distribution of the LULC classification in: (a) 2001, (b) 2010, and (c) 2019.
Land 10 01317 g005
Figure 6. Multitemporal change of LULC type coverage in 2001, 2010, and 2019.
Figure 6. Multitemporal change of LULC type coverage in 2001, 2010, and 2019.
Land 10 01317 g006
Figure 7. Spatial distribution of predicted LULC data between 2002 and 2009.
Figure 7. Spatial distribution of predicted LULC data between 2002 and 2009.
Land 10 01317 g007
Figure 8. Spatial distribution of the predicted LULC data between 2011 and 2018.
Figure 8. Spatial distribution of the predicted LULC data between 2011 and 2018.
Land 10 01317 g008aLand 10 01317 g008b
Figure 9. Proportion area of LULC type of LULC data between 2001 and 2019.
Figure 9. Proportion area of LULC type of LULC data between 2001 and 2019.
Land 10 01317 g009
Figure 10. Patterns of surface runoff volume of the three AMC conditions and annual rainfall data between 2001 and 2010.
Figure 10. Patterns of surface runoff volume of the three AMC conditions and annual rainfall data between 2001 and 2010.
Land 10 01317 g010
Figure 11. Relationship between surface runoff volume and annual rainfall: (a) AMC-I, (b) AMC-II, and (c) AMC-III.
Figure 11. Relationship between surface runoff volume and annual rainfall: (a) AMC-I, (b) AMC-II, and (c) AMC-III.
Land 10 01317 g011
Figure 12. Relationship between the observed and estimated runoff between 2001 and 2010 under the AMC-I condition at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Figure 12. Relationship between the observed and estimated runoff between 2001 and 2010 under the AMC-I condition at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Land 10 01317 g012
Figure 13. Relationship between the observed and estimated runoff between 2001 and 2010 under the AMC-II condition at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Figure 13. Relationship between the observed and estimated runoff between 2001 and 2010 under the AMC-II condition at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Land 10 01317 g013
Figure 14. Relationship between the observed and estimated runoff between 2001 and 2010 under the AMC-III condition at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Figure 14. Relationship between the observed and estimated runoff between 2001 and 2010 under the AMC-III condition at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Land 10 01317 g014
Figure 15. Pattern of surface runoff volume and annual rainfall of AMC-II between 2011 and 2019.
Figure 15. Pattern of surface runoff volume and annual rainfall of AMC-II between 2011 and 2019.
Land 10 01317 g015
Figure 16. Relationship between the observed and estimated runoff between 2011 and 2019 at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Figure 16. Relationship between the observed and estimated runoff between 2011 and 2019 at the three stations: (a) E.21 station, (b) E.23 station, and (c) E.6C station.
Land 10 01317 g016
Figure 17. Annual rainfall of Chaiyaphum meteorological station (1987–2019).
Figure 17. Annual rainfall of Chaiyaphum meteorological station (1987–2019).
Land 10 01317 g017
Figure 18. The 12-month SPI values of Chaiyaphum meteorological station (1987 and 2019).
Figure 18. The 12-month SPI values of Chaiyaphum meteorological station (1987 and 2019).
Land 10 01317 g018
Figure 19. Surface runoff reduction in 2029, 2039, and 2049 in dry, normal, and wet years.
Figure 19. Surface runoff reduction in 2029, 2039, and 2049 in dry, normal, and wet years.
Land 10 01317 g019
Figure 20. Spatial distribution of predicted LULC data in 2049: (a) Dry years, (b) normal years, and (c) wet years.
Figure 20. Spatial distribution of predicted LULC data in 2049: (a) Dry years, (b) normal years, and (c) wet years.
Land 10 01317 g020
Figure 21. Spatial distribution of economic value in 2049: (a) Actual LULC 2019, (b) dry years, (c) normal years, and (d) wet years.
Figure 21. Spatial distribution of economic value in 2049: (a) Actual LULC 2019, (b) dry years, (c) normal years, and (d) wet years.
Land 10 01317 g021
Figure 22. Contribution of the future economic value of LULC type of actual LULC and suitable LULC allocation for flood mitigation (dry, normal, and wet years).
Figure 22. Contribution of the future economic value of LULC type of actual LULC and suitable LULC allocation for flood mitigation (dry, normal, and wet years).
Land 10 01317 g022
Figure 23. Gain and loss of future economic value of suitable LULC allocation for flood mitigation in 2049: (a) Dry years, (b) normal years, and (c) wet years.
Figure 23. Gain and loss of future economic value of suitable LULC allocation for flood mitigation in 2049: (a) Dry years, (b) normal years, and (c) wet years.
Land 10 01317 g023
Figure 24. Spatial distribution of the ecosystem service value: (a) Actual LULC 2019, (b) dry years, (c) normal years, and (d) wet years.
Figure 24. Spatial distribution of the ecosystem service value: (a) Actual LULC 2019, (b) dry years, (c) normal years, and (d) wet years.
Land 10 01317 g024
Figure 25. Contribution of the ecosystem service value of each LULC type of actual LULC and suitable LULC allocation for flood mitigation (dry, normal, and wet years).
Figure 25. Contribution of the ecosystem service value of each LULC type of actual LULC and suitable LULC allocation for flood mitigation (dry, normal, and wet years).
Land 10 01317 g025
Figure 26. Gain and loss of ESV of suitable LULC allocation for flood mitigation in 2049: (a) Dry years, (b) normal years, and (c) wet years.
Figure 26. Gain and loss of ESV of suitable LULC allocation for flood mitigation in 2049: (a) Dry years, (b) normal years, and (c) wet years.
Land 10 01317 g026
Table 1. List of data collection and preparation for data analysis in this study.
Table 1. List of data collection and preparation for data analysis in this study.
DataData CollectionData PreparationSource
PrimaryGround reference--
SecondaryRunoff-RID
Annual rainfallInterpolationTMD
Socioeconomic dataPopulation densityDOPA
Income per capitaNESDC
Remote SensingLandsat 5 TM: Path 128 Row 49, 6 January 2001
Landsat 5 TM: Path 129 Row 49, 14 February 2001
Landsat 5 TM: Path 128 Row 49, 16 February 2010
Landsat 5 TM: Path 129 Row 49, 23 February 2010
Landsat 8 OLI: Path 128 Row 49, 24 January 2019
Landsat 8 OLI: Path 129 Row 49, 31 January 2019
-USGS
Satellite image from Google Earth in 2010-Google
Color orthophotograph-RTSD
GISAdministrative boundary-DEQP
Soil (soil series)RecodeLDD
Watershed boundary-RID
ElevationExtract from DEMSRTM
SlopeExtract from DEMSRTM
Road networkBufferingMOT, DEQP
StreamBufferingRTSD
Urban areaBufferingLULC data
Note: USGS, United States Geological Survey; RTSD, Royal Thai Survey Department; DEQP, Department of Environmental Quality Promotion; TMD, Thai Meteorological Department; RID, Royal Irrigation Department; LDD, Land Development Department; NESDC, Office of the National Economic and Social Development Council; SRTM, Shuttle Radar Topography Mission; MOT, Ministry of Transport; DOPA, Department of Provincial Administration.
Table 2. Runoff curve number under AMC-I, -II, and -III.
Table 2. Runoff curve number under AMC-I, -II, and -III.
LULC TypesCN Value of AMC-ICN Value of AMC-IICN Value of AMC-III
ABCDABCDABCD
Urban and built-up area77.2682.8586.8188.868992949594.996.3697.397.76
Cassava and other field crops51.9264.1675.4980.947281889185.5490.7594.495.88
Sugarcane 28.7548.3261.2468.84969798468.8583.6689.6492.35
Paddy fields43.8257.0868.875.496576848881.0387.9392.3594.4
Rangeland28.7548.3261.2468.84969798468.8583.6689.6492.35
Para rubber and perennial trees and orchards24.0643.8257.0865.684365768263.4481.0387.9391.29
Forest land15.2533.9249.4958.443055707749.6473.7684.2988.51
Waterbodies and marshes and swamps95.3795.3795.3795.379898989899.1299.1299.1299.12
Unused land58.4472.0780.9486.817786919488.5193.3995.8897.3
Table 3. Model performance scale.
Table 3. Model performance scale.
Statistics MeasurementPerformance Ratings
UnsatisfactorySatisfactoryGoodVery Good
NSE<0.50.5–0.650.65–0.750.75–1
R2<0.50.5–0.60.6–0.70.7–1
PBIAS>2515–2510–15<10
Table 4. LULC type and coefficient value for ESV evaluation.
Table 4. LULC type and coefficient value for ESV evaluation.
No.LULC Classification for RFLULC Classification for ESV 1Coefficient Values (USD/ha/year) 1
1Urban and built-up areasConstruction land12.7
2Paddy fieldsCultivated land1032.3
3Field cropsCultivated land1032.3
4Para rubberForest land1949.0
5Perennial trees and orchardsForest land1949.0
6Forest landForest land1949.0
7WaterbodiesWaterbodies6873.7
8RangelandRangeland808.6
9WetlandWetland9368.7
10Miscellaneous landUnused land96.3
1 Modified from [56].
Table 5. Area and percentage of LULC data in 2001, 2010, and 2019.
Table 5. Area and percentage of LULC data in 2001, 2010, and 2019.
No.LULC TypeLULC 2001LULC 2010LULC 2019
Area (km2)PercentageArea (km2)PercentageArea (km2)Percentage
1Urban and built-up areas46.171.2253.211.465.841.74
2Paddy fields2344.3961.792070.7154.582012.1653.03
3Sugarcane61.251.61153.524.05306.858.09
4Cassava532.9514.05629.3316.59489.9112.91
5Other field crops2.090.065.190.146.190.16
6Para rubber16.560.4430.050.7997.032.56
7Perennial trees and orchards55.761.4750.211.3288.952.34
8Forest land63216.66604.715.94481.312.68
9Waterbodies36.810.9757.461.5153.31.4
10Rangeland26.030.6972.111.971.651.89
11Marshes and swamps11.640.3133.40.8827.730.73
12Unused land28.570.7534.330.993.322.46
Total3794.22100.003794.22100.003794.22100.00
Table 6. Simple LULC change detection between 2001 and 2010.
Table 6. Simple LULC change detection between 2001 and 2010.
LULCLULC Type (Area, km2)
URPASUCAPCPRPOFOWARAMAUL
In 200146.172344.3961.25532.952.0916.5655.76632.0036.8126.0311.6428.57
In 201053.212070.71153.52629.335.1930.0550.21604.7057.4672.1133.4034.33
Change area7.04–273.6892.2796.383.1113.49–5.54–27.3020.6546.0821.765.76
Annual change rate0.78–30.4110.2510.710.351.50–0.62–3.032.295.122.420.64
Percentage of change0.19–7.212.432.540.080.36–0.15–0.720.541.210.570.15
Table 7. Simple LULC change detection between 2010 and 2019.
Table 7. Simple LULC change detection between 2010 and 2019.
LULCLULC Type (Area, km2)
URPASUCAPCPRPOFOWARAMAUL
In 201053.212070.71153.52629.335.1930.0550.21604.7057.4672.1133.4034.33
In 201965.842012.16306.85489.916.1997.0388.95481.3053.3071.6527.7393.32
Change area12.63–58.55153.33–139.411.0066.9938.73–123.41–4.16–0.47–5.6758.99
Annual change rate1.40–6.5117.04–15.490.117.444.30–13.71–0.46–0.05–0.636.55
Percentage of change0.33–1.544.04–3.670.031.771.02–3.25–0.11–0.01–0.151.55
Table 8. Simple LULC change detection between 2001 and 2019.
Table 8. Simple LULC change detection between 2001 and 2019.
LULCLULC Type (Area, km2)
URPASUCAPCPRPOFOWARAMAUL
In 200146.172344.3961.25532.952.0916.5655.7663236.8126.0311.6428.57
In 201965.842012.16306.85489.916.1997.0388.95481.353.371.6527.7393.32
Change area19.67–332.23245.6–43.044.180.4733.19–150.716.4945.6216.0964.75
Annual change rate1.09–18.4613.64–2.390.234.471.84–8.370.922.530.893.60
Percentage of change0.52–8.766.47–1.130.112.120.87–3.970.431.200.421.71
Table 9. Multiple linear equations of each LULC type location preference and area under the curve value from logistic regression analysis.
Table 9. Multiple linear equations of each LULC type location preference and area under the curve value from logistic regression analysis.
Driving FactorsLULC Type
URPDSUCAFCPROPFOWARAMAUL
Constant0.093010.883194.6560–38.4766207.978424.1870–38.483187.3484–4.6862–80.01185.0436–8.6344
X1n. s.–0.0947–0.00010.00420.01010.0083n. s.0.00320.0068–0.0043–0.04670.0093
X2n. s.n. s.–0.1943–0.1030–0.2144–0.0730–0.13140.1886–0.6715n. s.n. s.–0.0290
X3n. s.0.0130–0.12640.0477–0.2850–0.03980.0454–0.1167n. s.0.0991n. s.n. s.
X4n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.
X50.0008–0.0018–0.0117–0.0148n. s.–0.0186–0.0043–0.0113n. s.–0.0060–0.0059n. s.
X6–0.0064n. s.n. s.n. s.n. s.n. s.n. s.n. s.0.0014n. s.n. s.n. s.
X7n. s.n. s.n. s.0.0010n. s.–0.0014n. s.n. s.–0.0037n. s.–0.0053n. s.
X8–0.0182n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.n. s.
AUC0.98240.97630.76770.79420.93550.94400.61860.92890.83180.74450.86280.8693
Note: Elevation (X1), slope (X2), annual rainfall (X3), average income per capita at the sub-district level (X4), population density at the sub-district level (X5), distance to road network (X6), distance to stream (X7), and distance to the existing urban area (X8).
Table 10. Conversion matrix of possible changes between 2010 and 2019.
Table 10. Conversion matrix of possible changes between 2010 and 2019.
LULC TypesLULC Type Possible Change in 2019
URPASUCAFCPRPOFOWARAMAUL
LULC in 2010Urban and built-up areas (UR)100000000000
Paddy fields (PA)011100100110
Sugarcane (SU)111100100000
Cassava (CA)111101100101
Other field crops (FC)001110000000
Para rubber (PR)000111000001
Perennial trees and orchards (PO)001000100000
Forest land (FO)100001110101
Waterbodies (WA)101000001010
Rangeland (RA)010000100100
Marshes and swamps (MA)111100100010
Unused land (UL)110000000001
Table 11. Elasticity of LULC change for LULC prediction between 2010 and 2019.
Table 11. Elasticity of LULC change for LULC prediction between 2010 and 2019.
LULC TypesLULC Type Possible Change in 2019
URPASUCAFCPRPOFOWARAMAUL
LULC in 2010Urban and built-up areas (UR)1.00-----------
Paddy fields (PA)-0.930.030.02--0.01---0.01-
Sugarcane (SU)0.010.010.930.05--0.01-----
Cassava (CA)0.010.060.160.65-0.070.01----0.04
Other field crops (FC)----0.99-0.01-----
Para rubber (PR)---0.110.010.800.03----0.05
Perennial trees and orchards (PO)--0.010.01--0.99-----
Forest land (FO)-0.02-0.05-0.050.010.80---0.06
Waterbodies (WA)-0.040.020.01----0.91-0.01-
Rangeland (RA)-0.08-0.02--0.01--0.89--
Marshes and swamps (MA)0.010.550.03---0.02---0.39-
Unused land (UL)0.010.03---------0.96
Table 12. Annual land requirement of each LULC type between 2001 and 2010.
Table 12. Annual land requirement of each LULC type between 2001 and 2010.
YearLULC Type (Area, km2)Total
URPASUCAFCPRPOFOWARAMAUL
200146.172344.3961.25532.952.0916.5655.76632.0036.8126.0311.6428.573794.22
200246.952314.0171.53543.602.4518.0355.15628.9839.1231.1614.0329.193794.22
200347.712283.5981.77554.342.7919.5454.54625.9541.4036.2816.4629.843794.22
200448.502253.1891.96565.053.1321.0953.90622.9043.7041.3818.9430.473794.22
200549.302222.70102.27575.783.5022.5453.27619.8846.0046.5121.3031.163794.22
200650.052192.30112.51586.503.8324.1052.66616.8348.2851.6323.7831.743794.22
200750.822161.94122.76597.164.1725.5752.07613.8150.5756.7726.1732.403794.22
200851.612131.48133.03607.874.5127.0851.47610.7752.8661.9028.6033.033794.22
200952.412101.12143.27618.634.8728.5350.84607.7455.1666.9930.9733.693794.22
201053.212070.70153.51629.335.1930.0550.21604.7057.4672.1133.4034.333794.22
Annual rate0.79–30.4210.2410.710.331.52–0.63–3.032.295.122.440.64
Table 13. Area of predicted LULC type between 2002 and 2009.
Table 13. Area of predicted LULC type between 2002 and 2009.
LULC TypesArea of Predicted LULC Type (km2)
20022003200420052006200720082009
Urban and built-up areas (UR)46.9547.7148.549.350.0550.8251.6152.41
Paddy fields (PD)2314.012283.592253.182222.702192.302161.942131.482101.12
Sugarcane (SU)71.5381.7791.96102.27112.51122.76133.03143.27
Cassava (CA)543.6554.34565.05575.78586.5597.16607.87618.63
Other field crops (FC)2.452.793.133.53.834.174.514.87
Para rubber (PR)18.0319.5421.0922.5424.125.5727.0828.53
Perennial trees and orchards (PO)55.1554.5453.953.2752.6652.0751.4750.84
Forest land (FO)628.98625.95622.9619.88616.83613.81610.77607.74
Waterbodies (WA)39.1241.443.74648.2850.5752.8655.16
Rangeland (RA)31.1636.2841.3846.5151.6356.7761.966.99
Marshes and swamps (MA)14.0316.4618.9421.323.7826.1728.630.97
Unused land (UL)29.1929.8430.4731.1631.7432.433.0333.69
Total3794.223794.223794.223794.223794.223794.223794.223794.22
Table 14. Annual land requirement of each LULC type between 2010 and 2019.
Table 14. Annual land requirement of each LULC type between 2010 and 2019.
YearLULC Type (Area, km2)Total
URPASUCAFCPRPOFOWARAMAUL
201053.212070.71153.52629.335.1930.0550.21604.7057.4672.1133.4034.333794.22
201154.582064.15170.56613.865.3237.5154.49591.0257.0072.0932.7540.883794.22
201255.992057.68187.61598.355.4344.9358.80577.3056.5472.0332.1347.423794.22
201357.462051.17204.62582.865.5352.3963.17563.5356.0671.9831.4853.983794.22
201458.812044.70221.65567.395.6559.8367.42549.8755.6171.8530.9060.553794.22
201560.272038.20238.68551.895.7467.2671.78536.1455.1571.8130.2167.073794.22
201661.652031.64255.75536.395.8774.7176.06522.4354.7071.8129.5973.613794.22
201763.062025.15272.80520.905.9782.1480.37508.7454.2571.7128.9580.183794.22
201864.412018.69289.83505.426.0889.5984.63495.0253.7671.7128.3386.753794.22
201965.842012.16306.85489.916.1997.0388.95481.3053.3071.6527.7393.323794.22
Annual rate1.40–6.5117.04–15.490.117.444.30–13.71–0.46–0.05–0.636.55
Table 15. Area of predicted LULC type between 2011 and 2018.
Table 15. Area of predicted LULC type between 2011 and 2018.
LULC TypesArea of Predicted LULC Type (km2)
20112012201320142015201620172018
Urban and built-up areas (UR)54.5855.9957.4658.8160.2761.6563.0664.41
Paddy fields (PD)2064.152057.682051.172044.702038.202031.642025.152018.69
Sugarcane (SU)170.56187.61204.62221.65238.68255.75272.8289.83
Cassava (CA)613.86598.35582.86567.39551.89536.39520.9505.42
Other field crops (FC)5.325.435.535.655.745.875.976.08
Para rubber (PR)37.5144.9352.3959.8367.2674.7182.1489.59
Perennial trees and orchards (PO)54.4958.863.1767.4271.7876.0680.3784.63
Forest land (FO)591.02577.3563.53549.87536.14522.43508.74495.02
Waterbodies (WA)5756.5456.0655.6155.1554.754.2553.76
Rangeland (RA)72.0972.0371.9871.8571.8171.8171.7171.71
Marshes and swamps (MA)32.7532.1331.4830.930.2129.5928.9528.33
Unused land (UL)40.8847.4253.9860.5567.0773.6180.1886.75
Total3794.223794.223794.223794.223794.223794.223794.223794.22
Table 16. Annual surface runoff volume of the three AMC conditions and rainfall data between 2001 and 2010.
Table 16. Annual surface runoff volume of the three AMC conditions and rainfall data between 2001 and 2010.
YearSurface Runoff Volume (million m3)Annual Rainfall * (mm)
AMC-IAMC-IIAMC-III
20011178.41939.701537.88939.70
20024057.741191.604652.401191.60
20031444.46900.801836.21900.80
20041335.34915.401716.51915.40
20051754.451039.002185.221039.00
20064229.131196.004847.851196.00
20073234.721342.903795.551342.90
20085869.541695.206558.701695.20
20094189.381502.104804.031502.10
20105761.811506.306437.311506.30
* Annual rainfall data at Chaiyaphum meteorological station, located at the center of the study area.
Table 17. Annual surface runoff volume of the three AMC conditions and rainfall data between 2001 and 2010.
Table 17. Annual surface runoff volume of the three AMC conditions and rainfall data between 2001 and 2010.
YearE.21 Hydrological StationE.23 Hydrological StationE.6C Hydrological Station
QobsAMC-IAMC-IIAMC-IIIQobsAMC-IAMC-IIAMC-IIIQobsAMC-IAMC-IIAMC-III
QsimQsimQsimQsimQsimQsimQsimQsimQsim
2001633.4601.82695.9725.99373.1292.79354.69382.8736.6425.8245.1961.25
20021954.201820.051959.272001.761500.70929.331032.161074.94166.38113.11169.23204.97
2003656.7650.03746.73777.44575.6316.3379.93408.7450.6321.637.6451.78
2004549.9838.14945978.53542.4409.53482.1514.1341.2131.1653.4171.19
2005854.71077.711195.421232.06653.8530612.45648.284.3354.3890.48116.27
20062230.502439.222590.812636.541630.701246.651364.561412.1973.3735.2261.2281.35
20071914.301450.861581.281621.40845.5724.36819.84859.7585.4247.5882.78107.02
20083494.503173.603334.973383.211932.201625.721754.191804.98180.45146.87214.7256.12
20092180.202181.082328.332372.981560.001121.681234.521280.50151.03122.82184.17222.61
20103188.903033.463189.973237.172174.401570.831698.551749.15206.15181.51255.45299.02
Average1765.731726.601856.771896.711178.84876.72973.31013.54107.5678.01119.43147.16
Table 18. Statistical data of model performance for suitable AMC identification and model validation.
Table 18. Statistical data of model performance for suitable AMC identification and model validation.
AMCYearE.21 StationE.23 StationE.6C Station
NSER2PBIASNSER2PBIASNSER2PBIAS
AMC-II20010.950.960.180.670.960.680.710.961.01
20020.764.854.95
20030.042.202.70
2004–1.631.130.93
2005–1.261.052.78
2006–1.183.263.55
20072.621.033.52
20081.822.603.12
2009–0.013.722.62
20100.885.122.29
AMC-II20010.940.96–0.350.820.960.160.850.97–0.79
2002–0.033.97–0.27
2003–0.511.661.21
2004–2.240.51–1.13
2005–1.930.35–0.57
2006–2.042.261.13
20071.890.220.25
20080.901.51–3.18
2009–0.842.76–3.08
2010–0.014.04–4.58
AMC-III20010.940.96–0.520.860.96–0.080.320.97–2.29
2002–0.273.61–3.59
2003–0.681.42–0.11
2004–2.430.24–2.79
2005–2.140.05–2.97
2006–2.301.85–0.74
20071.66–0.12–2.01
20080.631.08–7.03
2009–1.092.37–6.66
2010–0.273.61–8.63
Table 19. Comparison of the average statistics measurement for suitable AMC examination.
Table 19. Comparison of the average statistics measurement for suitable AMC examination.
AMCAverage Statistics Measurement
NSER2PBIAS
AMC-I0.77670.96001.8443
AMC-II0.87000.96330.0423
AMC-III0.70670.9633–1.0067
Table 20. Annual surface runoff volume and rainfall data between 2011 and 2019.
Table 20. Annual surface runoff volume and rainfall data between 2011 and 2019.
YearSurface Runoff Volume (million m3)Annual Rainfall (mm)
20116142.431428.30
20123583.031087.20
20134200.481333.30
20141003.60793.50
20151233.80919.30
20163475.161044.20
20175433.331281.80
20182588.57809.00
20191445.54752.60
Table 21. Comparison of the observed (Qobs) and simulated (Qsim) surface runoff between 2011 and 2019 of the three stations.
Table 21. Comparison of the observed (Qobs) and simulated (Qsim) surface runoff between 2011 and 2019 of the three stations.
E.21 StationE.23 StationE.63 Station
YearQobsQsimQobsQsimQobsQsim
20112943.402919.232283.301954.57175.78158.26
20121989.101694.14677.20882.97112.15151.54
20131884.802058.571646.701492.0170.5979.15
2014339.30600.38330.40299.1719.8723.05
2015440.80547.04418.70272.6227.3253.76
20161563.301490.671058.30771.6687.2675.27
20172732.702668.971823.201563.06182.55166.29
2018923.601089.07665.80573.7664.4878.83
201960.10795.5084.20412.5538.2559.31
Average1430.791540.40998.64913.6086.4793.94
Note: Observed (Qobs) and simulated (Qsim) surface runoff in a million m3.
Table 22. Statistical model performance data for the surface runoff estimation between 2011 and 2019 at the three stations.
Table 22. Statistical model performance data for the surface runoff estimation between 2011 and 2019 at the three stations.
YearE.21 StationE.23 StationE.6C Station
NSER2PBIASNSER2PBIASNSER2PBIAS
20110.910.940.190.900.943.660.870.902.25
2012 2.29 –2.29 –5.06
2013–1.351.72–1.10
2014–2.030.35–0.41
2015–0.831.63–3.40
20160.563.191.54
20170.492.892.09
2018–1.281.02–1.84
2019–5.71–3.65–2.71
Table 23. The 12-month SPI values, drought classification, and rainfall conditions of Chaiyaphum meteorological station (1987–2019).
Table 23. The 12-month SPI values, drought classification, and rainfall conditions of Chaiyaphum meteorological station (1987–2019).
YearSPIDrought Classification 1Rainfall ConditionsYearSPIDrought Classification 1Rainfall Conditions
19870.77Mild wetWet year2004–0.88Mild droughtDry year
1988–0.21Near normalNormal year2005–0.32Near normalNormal year
1989–0.58Mild droughtDry year20060.32Near normalNormal year
1990–0.02Near normalNormal year20070.87Mild wetWet year
1991–0.35Near normalNormal year20082.06Extreme wetWet year
1992–0.57Mild droughtDry year20091.43Moderate wetWet year
1993–0.88Mild droughtDry year20101.45Moderate wetWet year
1994–0.12Near normalNormal year20111.18Moderate wetWet year
19951.09Moderate wetWet year2012–0.12Near normalNormal year
19960.36Mild wetWet year20130.84Mild wetWet year
1997–1.65Severe droughtDry year2014–1.48Moderate droughtDry year
1998–0.19Near normalNormal year2015–0.86Mild droughtDry year
1999–0.14Near normalNormal year2016–0.30Near normalNormal year
20002.14Extreme wetWet year20170.65Mild wetWet year
2001–0.76Mild droughtDry year2018–1.40Moderate droughtDry year
20020.31Near normalNormal year2019–1.70Severe droughtDry year
2003–0.95Mild droughtDry year
Note: 1 SPI drought classification by [51].
Table 24. Runoff coefficient and its average in dry years.
Table 24. Runoff coefficient and its average in dry years.
No.LULC TypeAverage Surface Runoff Coefficient (million m3/km2)
Dry YearNormal YearWet Year
1Urban and built-up areas0.451.031.46
2Paddy fields0.441.011.44
3Sugarcane0.380.901.30
4Cassava0.451.001.41
5Other field crops0.390.881.31
6Para rubber0.370.831.30
7Perennial trees and orchards0.380.841.31
8Forest land0.280.761.17
9Waterbodies0.491.041.51
10Rangeland0.420.921.36
11Marshes and swamps0.541.141.59
12Unused land0.461.021.41
Table 25. Existing and predicted areas of LULC in 2029, 2039, and 2049 for constraint setting.
Table 25. Existing and predicted areas of LULC in 2029, 2039, and 2049 for constraint setting.
LULC TypesExisting Area (km2)LULC in 2029 (km2)LULC in 2039 (km2)LULC in 2049 (km2)
Predicted AreaRemarkPredicted AreaRemarkPredicted AreaRemark
Urban and built-up areas (X1)65.8477.380.17% increase88.310.34% increase99.030.50% increase
Paddy fields (X2)2012.161947.120.03% decrease1879.300.06% decrease1812.160.09% decrease
Sugarcane (X3)306.85424.950.38% increase520.310.69% increase599.360.95% increase
Cassava (X4)489.91408.260.16% decrease358.650.26% decrease330.320.32% decrease
Other field crops (X5)6.197.720.27% increase9.45250.52% increase11.260.81% increase
Para rubber (X6)97.03134.430.38% increase152.020.56% increase164.410.69% increase
Perennial trees and orchards (X7)88.95125.360.40% increase160.120.80% increase193.131.17% increase
Forest land (X8)481.30382.630.20% decrease304.440.36% decrease242.230.49% decrease
Waterbodies (X9)53.3049.520.07% decrease47.010.13% decrease42.780.19% decrease
Rangeland (X10)71.6570.540.01% decrease69.910.02% decrease66.050.07% decrease
Marshes and swamps (X11)27.7325.130.09% decrease24.660.11% decrease21.610.07% decrease
Unused land (X12)93.31141.200.51% increase180.070.93% increase212.001.27% increase
Total3794.223794.22 3794.22 3794.22
Table 26. Summary of the constraint setting for each LULC type in 2029, 2039, and 2049.
Table 26. Summary of the constraint setting for each LULC type in 2029, 2039, and 2049.
LULC TypesConstraints Setting for Each LULC Type
Constrain in 2029Constrain in 2039Constrain in 2049
Urban and built-up areas (X1)65.84 ≤ X1 ≤ 77.3865.84 ≤ X1 ≤ 88.3165.84 ≤ X1 ≤ 99.03
Paddy fields (X2)2012.16 ≥ X2 ≥ 1947.122012.16 ≥ X2 ≥ 1879.302012.16 ≥ X2 ≥ 1812.16
Sugarcane (X3)306.85 ≤ X3 ≤ 424.95306.85 ≤ X3 ≤ 520.31306.85 ≤ X3 ≤ 599.36
Cassava (X4)489.91 ≥ X4 ≥ 408.26489.91 ≥ X4 ≥ 358.65489.91 ≥ X4 ≥ 330.32
Other field crops (X5)6.19 ≤ X5 ≤ 7.726.19 ≤ X5 ≤ 9.456.19 ≤ X5 ≤ 11.26
Para rubber (X6)97.03 ≤ X6 ≤ 134.4397.03 ≤ X6 ≤ 152.0297.03 ≤ X6 ≤ 164.41
Perennial trees and orchards (X7)88.95 ≤ X7 ≤ 125.3688.95 ≤ X7 ≤ 160.1288.95 ≤ X7 ≤ 193.13
Forest land (X8)481.30 ≥ X8 ≥ 382.63481.30 ≥ X8 ≥ 304.44481.30 ≥ X8 ≥ 242.23
Waterbodies (X9)53.30 ≥ X9 ≥ 49.5253.30 ≥ X9 ≥ 47.0153.30 ≥ X9 ≥ 42.78
Rangeland (X10)71.65 ≥ X10 ≥ 70.5471.65 ≥ X10 ≥ 69.9171.65 ≥ X10 ≥ 66.05
Marshes and swamps (X11)27.73 ≥ X11 ≥ 25.1327.73 ≥ X11 ≥ 24.6627.73 ≥ X11 ≥ 21.61
Unused land (X12)93.31 ≤ X12 ≤ 141.2093.31 ≤ X12 ≤ 180.0793.31 ≤ X12 ≤ 211.89
Table 27. Optimization of LULC allocation to minimize surface runoff in dry years.
Table 27. Optimization of LULC allocation to minimize surface runoff in dry years.
LULC Types2019202920392049
Area of LULC Surface Runoff Allocated LULC Surface Runoff Allocated LULC Surface Runoff Allocated LULC Surface Runoff
Urban and built-up areas65.8429.8465.8429.8465.8429.8465.8429.84
Paddy fields2012.16942.321947.12856.541879.30826.711812.16797.17
Sugarcane 306.8570.96387.24148.87454.94174.90517.13198.81
Cassava 489.91242.89408.26184.27359.65162.33330.32149.09
Other field crops6.191.756.192.396.192.396.192.39
Para rubber97.0319.65134.4349.86152.0256.39164.4160.99
Perennial trees and orchards88.9525.76125.3647.40160.1260.54193.1373.02
Forest land481.30159.21481.30136.03481.30136.03481.30136.03
Waterbodies53.3023.8149.5224.2947.0123.0642.7820.99
Rangeland71.6523.2370.5429.3669.9129.1066.0527.49
Marshes and swamps27.7312.6925.1313.5924.6613.3421.6111.69
Unused land93.3125.8593.3142.5793.3142.5793.3142.57
Total3794.221577.963794.221565.013794.221557.183794.221550.07
Note: Area of LULC in kilometers squared, surface runoff in million meters cubed, and allocated LULC in kilometers squared.
Table 28. Optimization of LULC allocation to minimize surface runoff in normal years.
Table 28. Optimization of LULC allocation to minimize surface runoff in normal years.
LULC Types2019202920392049
Area of LULC Surface RunoffAllocated LULC Surface RunoffAllocated LULC Surface RunoffAllocated LULC Surface Runoff
Urban and built-up areas65.8467.5865.8467.5865.8467.5865.8467.58
Paddy fields2012.162038.951947.121973.051879.301904.331812.161836.29
Sugarcane 306.85277.30385.71348.56451.67408.18512.06462.76
Cassava 489.91489.07408.26407.56359.65359.03330.32329.75
Other field crops6.195.437.726.789.458.3011.269.88
Para rubber97.0380.36134.43111.33152.02125.89164.41136.16
Perennial trees and orchards88.9574.28125.36104.68160.12133.70193.13161.27
Forest land481.30366.03481.30366.03481.30366.03481.30366.03
Waterbodies53.3055.2349.5251.3047.0148.7042.7844.33
Rangeland71.6565.9570.5464.9369.9164.3566.0560.80
Marshes and swamps27.7331.5625.1328.6024.6628.0621.6124.59
Unused land93.3195.3393.3195.3393.3195.3393.3195.33
Total3794.223647.073794.223625.723794.223609.483794.223594.76
Note: Area of LULC in kilometers squared, surface runoff in million meters cubed, and allocated LULC in kilometers squared.
Table 29. Optimization of LULC allocation to minimize surface runoff in wet years.
Table 29. Optimization of LULC allocation to minimize surface runoff in wet years.
LULC Types2019202920392049
Area of LULCSurface RunoffAllocated LULC Surface Runoff Allocated LULC Surface Runoff Allocated LULC Surface Runoff
Urban and built-up areas65.8495.8265.8495.8265.8495.8265.8495.82
Paddy fields2012.162889.201947.122795.811879.302698.431812.162602.02
Sugarcane 306.85400.38422.14552.77520.31678.91599.36782.05
Cassava 489.91688.86408.26574.05358.65505.70330.32464.46
Other field crops6.198.096.198.099.4512.3611.2614.71
Para rubber97.03126.18134.43174.82152.02197.69164.41213.81
Perennial trees and orchards88.95116.4288.95116.4290.48119.73103.83138.51
Forest land481.30563.35481.30563.35481.30563.35481.30563.35
Waterbodies53.3080.2450.0274.5548.0170.7743.7864.41
Rangeland71.6597.1270.5495.6069.9194.7666.0589.52
Marshes and swamps27.7344.1026.1339.9725.6639.2122.6134.37
Unused land93.31131.7493.31131.7493.31131.7493.31131.74
Total3794.225241.523794.225223.013794.225208.463794.225194.79
Note: Area of LULC in kilometers squared, surface runoff in million meters cubed, and allocated LULC in kilometers squared.
Table 30. Deviation of annual surface runoff after minimization by goal programming.
Table 30. Deviation of annual surface runoff after minimization by goal programming.
ItemsSurface Runoff Minimization (million m3)
Dry YearsNormal YearsWet Years
202920392049202920392049202920392049
Goal of annual surface runoff (million m3)1577.961577.961577.963647.073647.073647.075241.525241.525241.52
Annual surface runoff after optimization (million m3)1565.011557.181550.073625.723609.483594.765223.015208.465194.79
Deviation from goal (million m3)–12.95–20.78–27.89–21.34–37.59–52.31–18.52–33.06–46.73
Deviation from goal (%)–0.82–1.32–1.77–0.59–1.03–1.43–0.35–0.63–0.89
Table 31. Conversion matrix of the possible change in 2029, 2039, and 2049.
Table 31. Conversion matrix of the possible change in 2029, 2039, and 2049.
LULC TypesLULC Type Possible Change in 2029, 2039, and 2049
URPASUCAFCPRPOFOWARAMAUL
LULC in 2019Urban and built-up areas (UR)100000000000
Paddy fields (PA)011100000110
Sugarcane (SU)101101100000
Cassava (CA)111101100101
Other field crops (FC)000010000000
Para rubber (PR)001111000001
Perennial trees and orchards (PO)001100100000
Forest land (FO)000001110101
Waterbodies (WA)000000001010
Rangeland (RA)010100100100
Marshes and swamps (MA)011000100010
Unused land (UL)100000000001
Table 32. Elasticity of LULC change for LULC prediction between 2019 and 2049.
Table 32. Elasticity of LULC change for LULC prediction between 2019 and 2049.
LULC TypessLULC Type Possible Change in 2049
URPASUCAFCPRPOFOWARAMAUL
LULC in 2019Urban and built-up areas (UR)1.00-----------
Paddy fields (PA)-0.830.080.04--0.03---0.01-
Sugarcane (SU)0.020.030.830.09-0.010.02----0.01
Cassava (CA)0.030.120.310.31-0.110.03--0.01-0.08
Other field crops (FC)----0.97-0.02-----
Para rubber (PR)0.010.020.050.180.020.520.07----0.13
Perennial trees and orchards (PO)-0.010.020.01--0.96-----
Forest land (FO)0.010.010.030.04-0.100.040.70-0.01-0.04
Waterbodies (WA)0.010.100.070.03--0.01-0.76-0.01-
Rangeland (RA)-0.210.020.03--0.02--0.71--
Marshes and swamps (MA)0.010.060.080.03--0.05---0.77-
Unused land (UL)0.020.080.010.01--0.01----0.87
Table 33. Area and percentage of optimized LULC allocation to mitigate flood in 2049 in dry, normal, and wet years.
Table 33. Area and percentage of optimized LULC allocation to mitigate flood in 2049 in dry, normal, and wet years.
LULC TypeDry YearNormal YearWet Year
Area (km2)PercentageArea (km2)PercentageArea (km2)Percentage
Urban and built-up areas65.841.7465.841.7465.841.74
Paddy fields1812.1447.761812.1447.761812.1847.76
Sugarcane517.1213.63517.1213.63599.3415.80
Cassava330.338.71330.338.71330.348.71
Other field crops6.190.166.190.1611.250.30
Para rubber164.424.33164.424.33164.414.33
Perennial trees and orchards193.175.09193.175.09103.882.74
Forest land481.2612.68481.2612.68481.312.69
Waterbodies42.771.1342.771.1343.781.15
Rangeland66.071.7466.071.7466.011.74
Marshes and swamps21.570.5721.570.5722.560.59
Unused land93.32.4693.32.4693.282.46
Total3794.22100.003794.22100.003794.22100.00
Table 34. Areas of actual LULC in 2019 and suitable LULC allocation for flood mitigation in 2049 in dry, normal, and wet years.
Table 34. Areas of actual LULC in 2019 and suitable LULC allocation for flood mitigation in 2049 in dry, normal, and wet years.
No.LULC TypeActual LULC 2019Suitable LULC Allocation in 2049 (km2)
Dry YearsNormal YearsWet Years
1Urban and built-up areas65.8465.8465.8465.84
2Paddy fields2012.161812.141812.131812.18
3Sugarcane306.85517.12512.07599.34
4Cassava489.91330.33330.33330.34
5Other field crops6.196.1911.2711.25
6Para rubber97.03164.42164.41164.41
7Perennial trees and orchards88.95193.17193.15103.88
8Forest land481.30481.26481.30481.30
9Waterbodies53.3042.7742.8043.78
10Rangeland71.6566.0766.0466.01
11Marshes and swamps27.7321.5721.5822.56
12Unused land93.3293.3093.2893.28
Total3794.223794.223794.223794.22
Table 35. Present and future economic value of agricultural and forest land.
Table 35. Present and future economic value of agricultural and forest land.
LULC TypePrice
(Baht/ton)
Yield
(ton/km2)
Present Value
(Baht/ km2)
Discount Rate 5
(%)
Future Value in 2049 Using PV Model
(Baht/km2)
Paddy field 113,287.75218.752,906,695.316.5019,225,947.12
Sugarcane 2900.004468.754,021,875.006.5026,602,153.91
Cassava 11430.002240.633,204,093.756.5021,193,049.28
Other field crops 18092.50415.633,363,445.316.5022,247,058.87
Para rubber 343,685.83131.255,733,765.636.5037,925,225.34
Perennial trees and orchard 125,600.00247.246,329,440.006.5041,865,233.77
Forest land 4--25,000,000.006.50165,359,154.08
Source: 1 [72], 2 [73], 3 [74], 4 [75], and 5 discount rate was based on the minimum retail rate of Bank for [76].
Table 36. Economic value by LULC types of actual LULC 2019 and suitable LULC allocation for flood mitigation in 2049.
Table 36. Economic value by LULC types of actual LULC 2019 and suitable LULC allocation for flood mitigation in 2049.
LULC TypeEconomic Value in 2049 (Baht)
Actual LULC in 2019Suitable LULC Allocation for Flood Mitigation in 2049
Dry YearNormal YearWet Year
Paddy fields38,654.2034,629.2034,822.3034,826.09
Sugarcane8147.4213,678.8713,614.0315,920.81
Cassava10,347.306783.216979.086968.65
Other field crops137.81136.47200.36200.68
Para rubber3675.106145.566235.436222.15
Perennial trees and orchards3717.469041.538055.824404.83
Forest land79,388.7178,164.7979,303.8979,315.05
Total in million Baht144,068.00148,579.64149,210.92147,858.27
Table 37. Future economic value change in 2049 between actual LULC in 2019 and suitable LULC allocation in dry years.
Table 37. Future economic value change in 2049 between actual LULC in 2019 and suitable LULC allocation in dry years.
LULC TypeSuitable LULC Allocation in 2049 of Dry YearsTotal
PASUCAFCPRPOFO
Actual LULC in 2019Paddy fields (PA)01507.3028.210.0218.11602.051656.233811.92
Sugarcane (SU)–188.270–105.11–0.9914.05173.00790.00682.69
Cassava (CA)–28.21207.9600.28813.501983.531249.234226.29
Other field crops (FC)–0.010.91–0.11012.190.42104.33117.73
Para rubber (PR)0–2.68–25.43–8.9902.451167.211132.55
Perennial trees and orchards (PO)–193.53–123.84–86.03–0.35–5.100237.52–171.32
Forest land (FO)–131.26–2240.87–630.19–95.31–2467.68–305.090–5870.41
Total–541.27–651.22–818.66–105.34–1614.922456.365204.513929.45
Table 38. Future economic value change in 2049 between actual LULC in 2019 and suitable LULC allocation in normal years.
Table 38. Future economic value change in 2049 between actual LULC in 2019 and suitable LULC allocation in normal years.
LULC TypeSuitable LULC Allocation in 2049 of Normal YearsTotal
PASUCAFCPRPOFO
Actual LULC in 2019Paddy fields (PA)01453.4828.900.0238.03563.561656.233740.22
Sugarcane (SU)–190.750–108.03–1.0621.53151.87790.00663.57
Cassava (CA)–29.13224.3900.59975.321546.541252.863970.57
Other field crops (FC)–0.010.94–0.11012.050.41104.46117.74
Para rubber (PR)0–2.69–26.55–41.8502.451206.431137.79
Perennial trees and orchards (PO)–196.66–123.25–88.76–0.64–6.050237.74–177.62
Forest land (FO)–152.56–2237.12–645.50–142.20–1666.34–286.860–5130.59
Total–569.12–684.25–840.06–185.14–625.451977.965247.724321.67
Table 39. Future economic value change in 2049 between actual LULC in 2019 and suitable LULC allocation in wet years.
Table 39. Future economic value change in 2049 between actual LULC in 2019 and suitable LULC allocation in wet years.
LULC TypeSuitable LULC Allocation in 2049 of the Wet YearTotal
PASUCAFCPRPOFO
Actual LULC in 2019Paddy fields (PA)01525.9527.550.0258.57323.421656.233591.74
Sugarcane (SU)–190.820–105.43–1.0619.2797.94790.00609.90
Cassava (CA)–29.16565.9900.59971.75245.811252.863007.84
Other field crops (FC)-0.010.95–0.12011.920.35104.46117.55
Para rubber (PR)0–2.69–31.35–41.8902.451206.891133.40
Perennial trees and orchards (PO)–196.68–141.72–90.29–0.64–6.210237.74–197.81
Forest land (FO)–152.83–2360.51–653.55–143.61–1569.08–259.190–5138.76
Total–569.50–412.03–853.19–186.59–513.78410.785248.173123.86
Table 40. Area of each LULC type for ESV evaluation of actual LULC and suitable LULC allocation for flood mitigation under dry, normal, and wet years.
Table 40. Area of each LULC type for ESV evaluation of actual LULC and suitable LULC allocation for flood mitigation under dry, normal, and wet years.
No.ESV-LULC TypeActual LULC 2019Suitable LULC Allocation in 2049
Dry YearsNormal YearsWet Years
1Urban and built-up areas65.8465.8465.8465.84
2Paddy fields2012.161812.161812.161812.16
3Field crop802.95853.64853.64940.94
4Forest land667.28838.84838.84749.54
5Waterbodies53.342.7842.7843.78
6Rangeland71.6566.0566.0566.05
7Marshes and swamps27.7321.6121.6122.61
8Unused land93.3293.3193.3193.31
Total3794.223794.223794.223794.22
Table 41. Ecosystem service value by ESV-LULC types of actual LULC 2019 and suitable LULC allocation for flood mitigation in 2049.
Table 41. Ecosystem service value by ESV-LULC types of actual LULC 2019 and suitable LULC allocation for flood mitigation in 2049.
ESV-LULC TypeEcosystem Service Value (Baht)
Actual LULC in 2019Suitable LULC Allocation for Flood Mitigation in 2049
Dry YearsNormal YearsWet Years
Urban and built-up areas (UR)2.602.622.622.62
Paddy fields (PA)6478.735804.115836.475837.11
Field crops (FC)2580.202708.132739.313017.20
Forest land (FO)4050.695175.685088.764556.47
Waterbodies (WA)1143.071019.07923.23951.98
Rangeland (RA)180.65176.46167.47163.55
Marshes and swamps (MA)810.57567.00642.80679.99
Unused land (UL)27.9227.9427.9427.94
Total in million Baht15,274.4215,481.0015,428.6115,236.86
Table 42. Ecosystem service value change between actual LULC in 2019 and suitable LULC allocation in 2049 in dry years.
Table 42. Ecosystem service value change between actual LULC in 2019 and suitable LULC allocation in 2049 in dry years.
LULC TypeSuitable LULC Allocation in 2049 of Dry YearsTotal
URPAFCFOWARAMAUL
Actual LULC in 2019Urban and built-up areas (UR)033.0714.0314.064.880.172.050.0868.34
Paddy fields (PA)–36.3600111.3040.04–10.79119.19–1.86221.52
Field crops (FC)–14.6900495.1416.80–3.065.11–18.86480.44
Forest land (FO)–9.60–27.03–102.47013.13–3.090.65–91.86–220.28
Waterbodies (WA)–4.93–46.67–15.46–20.490–0.1439.86–0.82–48.65
Rangeland (RA)–0.168.884.3012.490.1900.75–0.0226.42
Marshes and swamps (MA)–1.97–184.81–6.07–249.41–0.91–1.6400–444.81
Unused land (UL)–0.071.288.82112.930.630.0100123.60
Total–67.79–215.28–96.84476.0174.76–18.53167.60–113.35206.58
Table 43. Ecosystem service value change between the actual LULC in 2019 and the suitable LULC allocation in 2049 in normal years.
Table 43. Ecosystem service value change between the actual LULC in 2019 and the suitable LULC allocation in 2049 in normal years.
LULC TypesSuitable LULC Allocation in 2049 of Normal YearsTotal
URPAFCFOWARAMAUL
Actual LULC in 2019Urban and built-up areas (UR)033.2114.2513.474.680.152.180.0868.02
Paddy fields (PA)–36.3600109.4931.33–10.30119.54–1.86211.83
Field crops (FC)–14.6900460.3016.74–2.724.87–18.86445.63
Forest land (FO)–9.60–27.85–110.14013.13–2.750.65–91.86–228.41
Waterbodies (WA)–4.93–46.67–16.38–19.960–0.1470.47–0.82–18.44
Rangeland (RA)–0.168.984.2720.190.1500.77–0.0234.19
Marshes and swamps (MA)–1.97–184.81–54.01–238.61–0.69–1.6400–481.73
Unused land (UL)–0.071.329.31111.900.630.0100123.09
Total–67.79–215.81–152.70456.7765.98–17.39198.47–113.35154.18
Table 44. Ecosystem service value change between the actual LULC in 2019 and the suitable LULC allocation in 2049 in wet years.
Table 44. Ecosystem service value change between the actual LULC in 2019 and the suitable LULC allocation in 2049 in wet years.
LULC TypesSuitable LULC Allocation in 2049 of Wet YearsTotal
URPAFCFOWARAMAUL
Actual LULC in 2019Urban and built-up areas (UR)033.2215.5910.914.760.142.100.0866.81
Paddy fields (PA)–36.360082.2734.28–10.04130.26–1.86198.55
Field crops (FC)–14.6900268.9316.79–2.643.65–18.86253.17
Forest land (FO)–9.60–27.85–117.37013.13–2.620.85–91.86–235.33
Waterbodies (WA)–4.93–46.67–19.20–17.570–0.1461.42–0.82–27.92
Rangeland (RA)–0.168.994.6821.780.1500.72–0.0236.14
Marshes and swamps (MA)–1.97–184.81–214.16–48.07–0.76–1.6400–451.40
Unused land (UL)–0.071.3210.00110.520.630.0100122.41
Total–67.79–215.81–320.47428.7868.99–16.92199.01–113.35–37.56
Table 45. Future economic and ESV value evaluation and change and reduction in surface runoff of each suitable LULC allocation for flood mitigation by comparison to the baseline information of LULC data in 2019.
Table 45. Future economic and ESV value evaluation and change and reduction in surface runoff of each suitable LULC allocation for flood mitigation by comparison to the baseline information of LULC data in 2019.
ItemSuitable LULC Allocation for Flood Mitigation in 2049
Dry YearsNormal YearsWet Years
Future economic value (million Baht)148,579.64149,210.92147,858.27
Gain or loss by economic value (million Baht)3929.454321.673123.86
Ecosystem service value (million Baht)15,481.0015,428.6115,236.86
Gain or loss by ESV (million Baht)206.58154.18–37.56
Runoff reduction (million m3)27.8952.3146.73
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Phinyoyang, A.; Ongsomwang, S. Optimizing Land Use and Land Cover Allocation for Flood Mitigation Using Land Use Change and Hydrological Models with Goal Programming, Chaiyaphum, Thailand. Land 2021, 10, 1317. https://doi.org/10.3390/land10121317

AMA Style

Phinyoyang A, Ongsomwang S. Optimizing Land Use and Land Cover Allocation for Flood Mitigation Using Land Use Change and Hydrological Models with Goal Programming, Chaiyaphum, Thailand. Land. 2021; 10(12):1317. https://doi.org/10.3390/land10121317

Chicago/Turabian Style

Phinyoyang, Athiwat, and Suwit Ongsomwang. 2021. "Optimizing Land Use and Land Cover Allocation for Flood Mitigation Using Land Use Change and Hydrological Models with Goal Programming, Chaiyaphum, Thailand" Land 10, no. 12: 1317. https://doi.org/10.3390/land10121317

APA Style

Phinyoyang, A., & Ongsomwang, S. (2021). Optimizing Land Use and Land Cover Allocation for Flood Mitigation Using Land Use Change and Hydrological Models with Goal Programming, Chaiyaphum, Thailand. Land, 10(12), 1317. https://doi.org/10.3390/land10121317

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