A Land-Use Change Model to Support Land-Use Planning in the Mekong Delta (MEKOLUC)

: Agricultural land-use changes pose challenges for land managers in terms of ensuring the implementation of local land-use plans. This paper aims to build a land-use change model named MEKOLUC (Mekong Delta land-use change) for simulating land-use changes under the impacts of socioeconomic factors (proﬁtability of land-use types, societal impacts on neighborhoods) and environmental factors (soil, salinity, persistence of salinity). The salinity diffusion map was generated using GAMA software and employed Markov cellular automata to depict the spread of salinity under the inﬂuence of dike and sluice gate system operations. The land-use decision-making process was based on multi-criteria selection of the main factors, which were land suitability, land convertibility, density of land use in the neighborhood and proﬁtability of land-use patterns. The input data for the case study were historical land-use maps from 2005, 2010 and 2015 of Soc Trang, a coastal province in the Mekong Delta. The model was calibrated using a land-use map from 2010 (with kappa = 0.86) and was veriﬁed with land-use maps from 2015 and 2020 with deviations from 0 to 19%. The simulated results showed that shrimp–rice farming areas have been shrinking, even though these are recommended as sustainable farming systems. Inversely, intensive rice crops tended to change to rice–vegetable crops, vegetable crops or perennial fruit trees, which are projected to be well adapted to climate and salinity intrusion by 2030. This case study shows that the developed model is an essential tool for helping land managers and farmers build land-use plans.


Introduction
The Mekong Delta of Vietnam comprises 13 provinces, which are located in the lower basin of the Mekong River. This is a major agricultural production area in Vietnam, combining many ecological zones: fresh, salty and brackish water. In this region, the production of rice, fruit and aquatic products takes place; the production in this area accounts for 55.7% of Vietnam's total rice production and more than 69% of the country's aquaculture production [1].
The Mekong Delta is expected to be strongly affected by climate change and a rise in sea levels [2]. Therefore, understanding land-use change in order to serve agricultural land-use planning plays an important role in the development of the region. Land-use plans are established every ten years, with a mid-term land-use inventory produced every five years to adjust the plan. However, the land use patterns observed at the end of the planning periods are often very different from those that are planned [3]. One of the main reasons for this is the lack of tools that help planners to predict land-use changes and test land-use solutions under the impact of different driving factors. as the provincial or regional scale, because of the question of spatial data size. However, cellular automata allowed users to define cell sizes that were appropriate to the scale of the maps. Thus, the cellular automata approach, with multi-criteria analysis for understanding land-use changes influenced by driving factors, was suitable for large-scale areas in the low lands and coastal areas.
The objectives of this study were as follows: (1) to build a simulation model based on cellular automata and multi-criteria analysis approaches that takes into account socioeconomic, salinity intrusion and environmental factors; (2) to use the model to understand future land-use changes as a basis for agricultural land-use planning; (3) to simulate land-use adaptation scenarios when environmental changes cause sea levels to rise. The model was designed based on the cellular automata modeling approach in the GAMA simulation platform [35]. This is a software that provides many powerful tools for large-scale, equation-based cellular automata and agent-based modeling, especially the ability to integrate complex GIS data with an intuitive design interface [36].

Materials and Methods
In this work, a new model named MEKOLUC 1 (Mekong Delta land-use change) was developed to study land-use changes and assess land-use strategies in order to support land-use planning. The conceptual model of MEKOLUC is shown in Figure 1. The input of the model includes a salinity raster map, dike and sluice gates, the profits of land-use types and the ability to convert between a land-use-type matrix and a land-use raster map. There were two main dynamics in each simulation step (one step was one year): salinity diffusion and land-use decision. The salinity diffusion was based on Markov-cellular automata with the saline neighbor rule, where a cell's salinity conversion was based on the proportion of neighborhood saline cells and limited by the dike and sluice gate systems. In each step, the salinity diffusion process updated the salinity level of the cells in the salinity layer. The salinity level was then combined with the salinity duration, soil type and acid sulfate type to create a land mapping unit that was used for calculating the land suitability of each land-use type. The second process was used to select the appropriate land-use type for each cell based on the multi-criteria evaluation among the land-use candidates.
Land 2021, 10, x FOR PEER REVIEW 3 many factors, with a detailed description of the current status of each land plot. How this level of detail was not effective when it was applied on a larger scale, such a provincial or regional scale, because of the question of spatial data size. However, ce automata allowed users to define cell sizes that were appropriate to the scale of the m Thus, the cellular automata approach, with multi-criteria analysis for understanding use changes influenced by driving factors, was suitable for large-scale areas in th lands and coastal areas. The objectives of this study were as follows: (1) to build a simulation model bas cellular automata and multi-criteria analysis approaches that takes into account soci nomic, salinity intrusion and environmental factors; (2) to use the model to under future land-use changes as a basis for agricultural land-use planning; (3) to simulate use adaptation scenarios when environmental changes cause sea levels to rise.

Overview of the Model
The model was designed based on the cellular automata modeling approach GAMA simulation platform [35]. This is a software that provides many powerful too large-scale, equation-based cellular automata and agent-based modeling, especial ability to integrate complex GIS data with an intuitive design interface [36].
In this work, a new model named MEKOLUC 1 (Mekong Delta land-use change developed to study land-use changes and assess land-use strategies in order to su land-use planning. The conceptual model of MEKOLUC is shown in Figure 1. The of the model includes a salinity raster map, dike and sluice gates, the profits of lan types and the ability to convert between a land-use-type matrix and a land-use raster There were two main dynamics in each simulation step (one step was one year): sa diffusion and land-use decision. The salinity diffusion was based on Markov-cellul tomata with the saline neighbor rule, where a cell's salinity conversion was based o proportion of neighborhood saline cells and limited by the dike and sluice gate sys In each step, the salinity diffusion process updated the salinity level of the cells i salinity layer. The salinity level was then combined with the salinity duration, soi and acid sulfate type to create a land mapping unit that was used for calculating the suitability of each land-use type. The second process was used to select the appro land-use type for each cell based on the multi-criteria evaluation among the land-us didates.   The diagram in Figure 2 shows the main objects in the model: the world, cells, land-use types, dike-protected regions, land units and the land-suitability dataset. The main objects are the cells. Each cell belongs to a land unit. The cells belong to the area inside or outside the dike. The land suitability of the cells is determined by the composition of land-use types and land units. The land unit provides information on soil type, salinity level and time of saline intrusion. Land units are associated with land-use patterns to provide a level of land suitability. Each cell also has a list of neighboring cells. Each land-use type provides the benefit and implementation costs, along with a map that measures how convenient it is to convert from one land use to another. The diagram in Figure 2 shows the main objects in the model: the world, cells, landuse types, dike-protected regions, land units and the land-suitability dataset. The main objects are the cells. Each cell belongs to a land unit. The cells belong to the area inside or outside the dike. The land suitability of the cells is determined by the composition of landuse types and land units. The land unit provides information on soil type, salinity level and time of saline intrusion. Land units are associated with land-use patterns to provide a level of land suitability. Each cell also has a list of neighboring cells. Each land-use type provides the benefit and implementation costs, along with a map that measures how convenient it is to convert from one land use to another. Object diagram of MEKOLUC. * means many in one-to-many relationship, one object in an entity can be associated with one or more objects in another entity.

Multi-Criteria Land-Use Decision
The main land-use types (LUTs) selected are the dominant agricultural land-use types in the Mekong Delta, which include land growing rice, rice-vegetable rotations, vegetables, fruit trees, rice-shrimp rotations and shrimp.
Based on previous studies [6,9,34], the selection of land-use type for each cell is based on multi-criteria decision making. There are four selected criteria affecting land-use change that were taken into account.
• Land suitability: This criterion represents the suitability of the land's environmental conditions (soil and water) to a particular land use. As environmental conditions change, the adaptation level also changes. Land suitability assessments were carried out in accordance with the FAO soil assessment process (FAO, 1976) [8]. The adaptability of the soil was based on the characteristics of the soil, water, topography and land cover, and was classified into four levels: S1-highly suitable; S2-suitable; S3-slightly suitable; and N-not suitable. The four suitability levels were standardized from 0 (not suitable) to 1 (most suitable).
• Changeable: This criterion measures how convenient it is to convert from the current land use to another (Table 1). In farming systems, the transition from one type to another depends on the conditions in which new farming types can be established. In this study, the ease of conversion was determined based on the similarity between the cultivation types. Convertibility was coded on a Likert scale with five levels (0, 0.25, 0.5, 0.75, 1), where 0 was non-changeable and 1 was easily changeable. Object diagram of MEKOLUC. * means many in one-to-many relationship, one object in an entity can be associated with one or more objects in another entity.

Multi-Criteria Land-Use Decision
The main land-use types (LUTs) selected are the dominant agricultural land-use types in the Mekong Delta, which include land growing rice, rice-vegetable rotations, vegetables, fruit trees, rice-shrimp rotations and shrimp.
Based on previous studies [6,9,34], the selection of land-use type for each cell is based on multi-criteria decision making. There are four selected criteria affecting land-use change that were taken into account.
• Land suitability: This criterion represents the suitability of the land's environmental conditions (soil and water) to a particular land use. As environmental conditions change, the adaptation level also changes. Land suitability assessments were carried out in accordance with the FAO soil assessment process (FAO, 1976) [8]. The adaptability of the soil was based on the characteristics of the soil, water, topography and land cover, and was classified into four levels: S1-highly suitable; S2-suitable; S3-slightly suitable; and N-not suitable. The four suitability levels were standardized from 0 (not suitable) to 1 (most suitable).
• Changeable: This criterion measures how convenient it is to convert from the current land use to another (Table 1). In farming systems, the transition from one type to another depends on the conditions in which new farming types can be established. In this study, the ease of conversion was determined based on the similarity between the cultivation types. Convertibility was coded on a Likert scale with five levels (0, 0.25, 0.5, 0.75, 1), where 0 was non-changeable and 1 was easily changeable. • Density: The density of a LUT candidate in the neighborhood affects the selection of LUT. The higher the percentage is of a LUT candidate in neighboring cells, the greater the influence it has on the selection.
• Profit: The profitability of a LUT is one of the factors that leads to changes. This factor expresses the willingness of farmers to increase their income.
At each simulation step, each cell first evaluated the factors, calculated the convertibility index for each candidate using a weighted mean equation of the four criteria (Equation (1)) and then selected the optimized LUT among the candidates. where: • LUT i : current land use of the cell.
• LUT j : candidate land-use type; LUT j is one of the land-use types in the land-use map.
• C D , C LS , C P , C C : the four criteria-namely, density, land suitability, profit and changeability, respectively.
• W D , W LS , W P , W C : weights of criteria C D , C LS , C P , and C C , respectively.  The LUT scale values of neighboring cells were calculated based on the 8-cell Moore topology for each land-use type. For each cell, the value of each land-use type was calculated according to the proportion of the neighboring cell.

Case Study
Soc Trang is a coastal province in the Mekong Delta ( Figure 3) with an area of 3311.87 km 2 and approximately 1.3 million inhabitants. This province is located in the estuary area of the Hau river, a tributary of the Mekong river. The northeast region of the province is bordered by the Hau River and the southeast region is bordered by the sea. This region is also characterized by dike systems that prevent saline intrusion. The province is located in different specific agro-ecological zones of the delta: the fresh water region in the north (composed of the Ke Sach, Chau Thanh and My Tu districts), the coastal region protected by dike systems (composed of Long Phu, Tran De and Cu Lao Dung), the brackish water region (with 6 months of salt water) in the southwest (My Xuyen) and the saline water region in the south (Vinh Chau). With many ecological zones, the main land-use types of the region are very diverse and include rice cultivation, rice-vegetable rotations, rice-shrimp rotations, shrimp, fruit trees, vegetables and mangroves.
Land 2021, 10, x FOR PEER REVIEW 6 of 16 the region are very diverse and include rice cultivation, rice-vegetable rotations, riceshrimp rotations, shrimp, fruit trees, vegetables and mangroves. Soc Trang is also one of the provinces in the Mekong Delta that has been strongly impacted by saline intrusion and drought caused by extreme weather. In 2016 and 2020, saline water penetrated deeply into the river mouth of the Hau River, causing freshwater areas to become saline and drought in the dike regions, which seriously affected crop cultivation in the dry season of the province [38,39]. These irregular phenomena show that it is necessary to study and project land-use changes in order to support land-use planning towards adapting to environmental changes.

Socioeconomic and Environmental Data
The  Soc Trang is also one of the provinces in the Mekong Delta that has been strongly impacted by saline intrusion and drought caused by extreme weather. In 2016 and 2020, saline water penetrated deeply into the river mouth of the Hau River, causing freshwater areas to become saline and drought in the dike regions, which seriously affected crop cultivation in the dry season of the province [38,39]. These irregular phenomena show that it is necessary to study and project land-use changes in order to support land-use planning towards adapting to environmental changes. For economic data, profitability of land-use types are collected based on household surveys. These surveys showed that the profit of shrimp was the highest 2 (389 million VND/ha), followed by fruit trees (294 million VND/ha), vegetables (180 million VND/ha) and rice-vegetable rotations (64 million VND/ha). Rice gave the lowest profit at 34 million VND/ha.  For economic data, profitability of land-use types are collected based on household surveys. These surveys showed that the profit of shrimp was the highest 2 (389 million VND/ha), followed by fruit trees (294 million VND/ha), vegetables (180 million VND/ha) and rice-vegetable rotations (64 million VND/ha). Rice gave the lowest profit at 34 million VND/ha.

Socioeconomic and Environmental Data
To study the change in land-use status up to 2030 under the impact of climate change (a rise in sea level), this study uses additional input data for the model from the saline intrusion scenario map due to the sea-level rise (SLR) under normal conditions by 2030 [41].

Calibration and Validation Methods
We used the genetic algorithm (installed in GAMA) to find the best combination of parameters-i.e., the one that maximizes kappa between the simulated data and the observed data (land use in 2010). Each parameter value could vary from 0.1 to 1, with a step of 0.1. To study the change in land-use status up to 2030 under the impact of climate change (a rise in sea level), this study uses additional input data for the model from the saline intrusion scenario map due to the sea-level rise (SLR) under normal conditions by 2030 [41].

Calibration and Validation Methods
We used the genetic algorithm (installed in GAMA) to find the best combination of parameters-i.e., the one that maximizes kappa between the simulated data and the observed data (land use in 2010). Each parameter value could vary from 0.1 to 1, with a step of 0.1.
The fuzzy-kappa [42] action installed in GAMA was used to calculate the precision of the simulation, including the root mean square error (RMSE) value that was used to validate the simulated area with the historical land-use statistic. where Area sim is the simulated area of each land-use type; Area statistic is the statistical landuse area; and T is the number of simulated land-use types.

Historical Land-Use Changes
To analyze past transitions between land-use types, the areas of land-use types in the period 2005-2020 were analyzed ( Figure 5).
The fuzzy-kappa [42] action installed in GAMA was used to calculate the preci of the simulation, including the root mean square error (RMSE) value that was use validate the simulated area with the historical land-use statistic.
where Areasim is the simulated area of each land-use type; Areastatistic is the statistical la use area; and T is the number of simulated land-use types.

Historical Land-Use Changes
To analyze past transitions between land-use types, the areas of land-use types in period 2005-2020 were analyzed ( Figure 5).  Figure 5 shows the trend of shifting away from rice production (including r shrimp rotations) to other land uses. In the first period (2005-2010), the total land area shifted away from rice production represented more than 13,800 ha. At the same time area used for shrimp aquaculture increased from 45,054 ha to 54,519 ha and the area u for vegetables and perennial trees slightly increased. Young et al. (2002) [47] showed in early 2000 the market price of rice was near or below its production cost, which plained why a huge number of farmers shifted their land use away from rice. In the n period (2010-2015), the area used for the production of rice slightly increased to aro 149,000 ha, then decreased again in the period 2015-2020. After the initial increase in 2 the aquaculture area slightly increased over the next 10 years because the regions in the dikes were well protected and hence suitable for rice production. In this period o years, the forested area decreased from more than 12,000 ha to 9820 ha.
Historical data have revealed a difficult problem for land managers: while farm tend to convert their land use to intensive rice production in the dike system or transi from rice-shrimp rotations to shrimp farming in the outer dike areas, the areas out  Figure 5 shows the trend of shifting away from rice production (including rice-shrimp rotations) to other land uses. In the first period (2005-2010), the total land area that shifted away from rice production represented more than 13,800 ha. At the same time, the area used for shrimp aquaculture increased from 45,054 ha to 54,519 ha and the area used for vegetables and perennial trees slightly increased. Young et al. (2002) [47] showed that in early 2000 the market price of rice was near or below its production cost, which explained why a huge number of farmers shifted their land use away from rice. In the next period (2010-2015), the area used for the production of rice slightly increased to around 149,000 ha, then decreased again in the period 2015-2020. After the initial increase in 2010, the aquaculture area slightly increased over the next 10 years because the regions inside the dikes were well protected and hence suitable for rice production. In this period of 15 years, the forested area decreased from more than 12,000 ha to 9820 ha.
Historical data have revealed a difficult problem for land managers: while farmers tend to convert their land use to intensive rice production in the dike system or transition from rice-shrimp rotations to shrimp farming in the outer dike areas, the areas outside the dikes were converted to shrimp production, which raised the question of how to build infrastructure plans and land-use strategies for sustainable development.

Salinity Diffusion
Agricultural activities in Soc Trang province are heavily influenced by local policies on the management of dikes and sluice gates, which prevent saline intrusion. The salinity transmission mechanism was built using a Markov chain based on salinity maps from 2005 and 2015, where the salinity data were divided into five levels: <4‰; 4-8‰; 8-12‰; 12-20‰; and greater than 20‰. The mechanism of salinity propagation was constructed according to a Markov-cellular automata for the areas outside the dikes. Figure 6 shows the results of salinity diffusion, where the map on the left shows the simulated results of salinity in 2014 (Figure 6b) compared to the actual salinity map in 2014 (Figure 6a). In the coastal areas of the Tran De district, the salinity changed due to the expanded dike boundary region surrounding the shrimp farming areas. For the rest of the areas, the saline area, at 4-8‰, was opened up to the dike and the 8-12‰ salinity zone did not change during this period. This phenomenon showed that under normal climate conditions, the dikes mainly affected the saline intrusion in this region. Therefore, the dike management policy in the current situation affects the water environment and will affect future land use. the dikes were converted to shrimp production, which raised the question of how to build infrastructure plans and land-use strategies for sustainable development.

Salinity Diffusion
Agricultural activities in Soc Trang province are heavily influenced by local policies on the management of dikes and sluice gates, which prevent saline intrusion. The salinity transmission mechanism was built using a Markov chain based on salinity maps from 2005 and 2015, where the salinity data were divided into five levels: <4‰; 4-8‰; 8-12‰; 12-20‰; and greater than 20‰. The mechanism of salinity propagation was constructed according to a Markov-cellular automata for the areas outside the dikes. Figure 6 shows the results of salinity diffusion, where the map on the left shows the simulated results of salinity in 2014 (Figure 6b) compared to the actual salinity map in 2014 (Figure 6a). In the coastal areas of the Tran De district, the salinity changed due to the expanded dike boundary region surrounding the shrimp farming areas. For the rest of the areas, the saline area, at 4-8‰, was opened up to the dike and the 8-12‰ salinity zone did not change during this period. This phenomenon showed that under normal climate conditions, the dikes mainly affected the saline intrusion in this region. Therefore, the dike management policy in the current situation affects the water environment and will affect future land use.

Calibration
In this case study, the MEKOLUC model was calibrated with a land-use map of Soc Trang province from 2010. The calibration process determined a set of weights (WD, WLS, WP, WC), where the simulated land-use map fitted the observed map from 2010 well. The calibrated results showed that the set of parameters for the criteria were set to 0.6, 0.7, 0.8 and 0.2 according to the weight of land suitability, changeability, density and profitability of the use types. The best kappa coefficient was 0.86.

Validation
The simulated results were then compared to the 2015 and 2020 land-use inventory datasets. Two indicators were taken into account: the difference between the simulated area and the statistical area in 2015 and 2020, and the fuzzy kappa values. Table 2 shows the deviations of the simulated values in 2015 and 2020. The deviation between the two

Calibration
In this case study, the MEKOLUC model was calibrated with a land-use map of Soc Trang province from 2010. The calibration process determined a set of weights (W D , W LS , W P , W C ), where the simulated land-use map fitted the observed map from 2010 well. The calibrated results showed that the set of parameters for the criteria were set to 0.6, 0.7, 0.8 and 0.2 according to the weight of land suitability, changeability, density and profitability of the use types. The best kappa coefficient was 0.86.

Validation
The simulated results were then compared to the 2015 and 2020 land-use inventory datasets. Two indicators were taken into account: the difference between the simulated area and the statistical area in 2015 and 2020, and the fuzzy kappa values. Table 2 shows the deviations of the simulated values in 2015 and 2020. The deviation between the two datasets showed small values of 0% to 19% for rice and fruit trees. The RSME values were 2087.16 ha and 2379.54 ha (per total of 240,034 ha) for 2015 and 2020, respectively. Since 2015 and 2020, in the statistics data, rice was not classified in detail (as rice or a rice-shrimp rotation), the rice area in Table 2 includes rice, rice-vegetable rotations and rice-shrimp rotations. The map in Figure 7 shows a comparison between the simulated land-use map and the observed map of Soc Trang province in 2015. The obtained fuzzy kappa coefficient was 0.88. The simulated map showed that there was a big gap in the estuary area of the Long Phu district. This area was simulated as a rice-vegetable rotation, but in reality it was still a rice field. Nevertheless, the results also showed that the area of rice land decreased and the areas used for the production of vegetables and aquatic products both increased.
2087.16 ha and 2379.54 ha (per total of 240,034 ha) for 2015 and 2020, respectively. Since 2015 and 2020, in the statistics data, rice was not classified in detail (as rice or a rice-shrimp rotation), the rice area in Table 2 includes rice, rice-vegetable rotations and rice-shrimp rotations. The map in Figure 7 shows a comparison between the simulated land-use map and the observed map of Soc Trang province in 2015. The obtained fuzzy kappa coefficient was 0.88. The simulated map showed that there was a big gap in the estuary area of the Long Phu district. This area was simulated as a rice-vegetable rotation, but in reality it was still a rice field. Nevertheless, the results also showed that the area of rice land decreased and the areas used for the production of vegetables and aquatic products both increased.

Land Use in 2030 under the Conditions of Climate Change
As the model was verified with the observed land-use map, the model was used to estimate land use up until 2030 under the impact of climate change (in the case of a rise in sea level). In this scenario, the predicted saline intrusion (following a rise in sea level) in 2030 was used (Figure 8). In the case of Soc Trang province, it was assumed that the dike system and the level of operation of the dike lines were kept constant.

Land Use in 2030 under the Conditions of Climate Change
As the model was verified with the observed land-use map, the model was used to estimate land use up until 2030 under the impact of climate change (in the case of a rise in sea level). In this scenario, the predicted saline intrusion (following a rise in sea level) in 2030 was used (Figure 8). In the case of Soc Trang province, it was assumed that the dike system and the level of operation of the dike lines were kept constant.
The results of the simulated land-use change up until 2030 are shown in Figure 9. The main trend showed there would be a slow conversion from rice and rice-shrimp rotations to other types. Rice-vegetable rotations, fruit trees and vegetable areas are predicted to increase until 2025, but then in the next 5-year period these land-use types are predicted to decrease. All of the decreasing areas are shown to be shifted to shrimp. The results of the simulated land-use change up until 2030 are shown in Figure 9. The main trend showed there would be a slow conversion from rice and rice-shrimp rotations to other types. Rice-vegetable rotations, fruit trees and vegetable areas are predicted to increase until 2025, but then in the next 5-year period these land-use types are predicted to decrease. All of the decreasing areas are shown to be shifted to shrimp.
The reason for the slow conversion is the lack of water for irrigation in the dry season. When salinity increases, even saline water is protected by the dike systems, which cause a lack of water for irrigation in the dry season. This situation affects the adaptation of plants. This shows that the government needs to invest in infrastructure to improve irrigation systems, even when the farming system uses less water for irrigation.

Discussion
In this section, we discuss the points that arose when using the model throughout the application case.
The advantage of this study is that it developed a model of land-use change that is The reason for the slow conversion is the lack of water for irrigation in the dry season. When salinity increases, even saline water is protected by the dike systems, which cause a lack of water for irrigation in the dry season. This situation affects the adaptation of plants. This shows that the government needs to invest in infrastructure to improve irrigation systems, even when the farming system uses less water for irrigation.

Discussion
In this section, we discuss the points that arose when using the model throughout the application case.
The advantage of this study is that it developed a model of land-use change that is adapted to the characteristics of the Mekong Delta, which will be affected by saline intrusion on a large scale compared with the highland region [33,48]. The model used cellular automata with a multi-criteria method that has been promoted in many similar studies [18,22], but with a new approach of coupling salinity diffusion with land-use changes. In our study, the multi-criteria assessment took into account the main driving factors of land-use changes, including economic (profitability), social (neighborhood and land use changeability) and environmental (soil, salinity, land suitability) factors. This helps planners to understand the impact of these factors on the lowlands of developing countries. Of these criteria, land suitability is an important criterion for choosing the type of land use where there are salt water intrusion dynamics; in this case, land suitability was re-evaluated during simulation. The multi-criteria assessment process then used these updated environmental data to increase the accuracy of the model.
In terms of accuracy, the model automatically calibrated the weight of the criteria experiment based on the genetic algorithm and batch experiment provided by the GAMA platform. These were convenient and easy to apply because of the existing historical land-use maps, as well as being more precise, compared with using the AHP approach for calculating the weights of the criteria [22,49].
The verification results with the land-use maps of Soc Trang province from 2015 and 2020 showed that the maximum deviation was 19% for fruit trees. The over-estimated area in this case meant that transforming land use to fruit trees was, in reality, slower than expected, and the model was not a good simulation for this type of land-use. There may have been some other factors, such as capital, affecting fruit trees land use.
The simulation results showed that the agricultural land area tended to shrink and be converted to aquaculture. Specifically, the area of rice-shrimp rotations, which is considered to be a form of environmentally sustainable farming [50,51], shrank due to the profit-driven behavior of the people. The model allowed planners to forecast the likelihood that, when environmental factors change, people's need to increase their income will be increased. The trend of transformation seems inevitable, which requires the government to have solutions for minimizing the damage caused by inexperienced people when converting to a new type of farming, such as investment capital for farmers.
Regarding input data, the land-use data sourced from the land-use map in the case study were based on the results of the land inventory, which is a land management requirement in Vietnam. However, the land-use types could easily be replaced by ones deduced from satellite images in similar application cases. In addition, the system of sluices and dikes, used to prevent salinity, helped to protect the rice growing area. This could be modified to test land-use solutions. Thus, planners could test policies related to environmental data by changing the dike-protection boundaries. As local management policies change, the choice of land use will change according to economic and environmental data.
During the extreme dry season [52][53][54], the source of fresh water on the canals was exhausted, and high temperatures and a lack of precipitation, combined with highly saline water, affected agricultural activities. However, the current model did not assess the risks of transitioning. These points should be explored and, supplemented by the possibility of a model-scale risk assessment during land-use change. Thus, it is necessary to extend the properties of the land cell to add more attributes concerning the capital and the level of farming technique of households to identify the risks involved, as well as quantify the required capital and technical support in the case of land-use conversion.
The model was designed to be completely open with the source code being published, allowing additional research to take place to complete the model.

Conclusions
In this study, a new land-use change model named MEKOLUC was proposed based on the coupling of salinity diffusion and land-use dynamics using cellular automata. Land-use decisions for each cell were based on multi-criteria decision making using four criteria: land suitability, land-use convertibility, land-use type in the neighborhood and the profitability of land-use patterns.
The model was applied to Soc Trang province in the Mekong Delta, Vietnam. In the case study, the model was verified using the 2015 and 2020 land-use statistical data of Soc Trang province, with RMSE values of 2087.16 ha and 2379.54 ha, respectively, and the fuzzy kappa value in 2015 was 0.88.
The simulated results of the land-use change scenario in 2030, according to the sealevel rise scenario, showed that the shrimp-rice rotation areas will decrease and aquaculture areas will increase to the maximum available size for the area outside the dike. In the area protected by the dike, people will gradually convert from rice cultivation to vegetable and fruit tree cultivation to adapt to climate change. This scenario is an example of using MEKOLUC to help managers and planners understand the potential for land-use change and test planning solutions to select land-use that will be most suitable for the environmental conditions.
Regarding the criteria, the model was predefined with four basic criteria, which could be improved with additional criteria to improve the precision of the simulation model. Additional criteria, such as the capital capacity and agricultural skills of households in the region, could be added to improve the accuracy of the model. In addition, climatic factors affecting land use were not included; it is therefore necessary to integrate climate factors related to climate change, such as temperature and precipitation, into further studies to account for the dynamic interaction between climate change and land use change. This will help residents and managers understand adaptation plans and test appropriate plans for each study area.
The model was designed for the Mekong Delta, a coastal region; hence, the model can be applied in delta areas with similar conditions. In addition, the model involves fully open-source code that operates on many different operating systems, so it could be extended to be applied to many land-use studies.