Journal of Geographical Sciences >
Simulating land use for protecting food crop areas in northeast Thailand using GIS and Dyna-CLUE
Author: K. Sakayarote, E-mail: kanda839@hotmail.com
Received date: 2018-08-04
Accepted date: 2018-12-20
Online published: 2019-04-19
Copyright
Land use in the northeast region of Thailand has changed dramatically in the past two decades. These changes are mainly due to the government policies, which launched a scheme to promote rubber plantation during 2003-2013 targeting to solve the problem of poverty in the region. At least 50,000 ha of paddy fields were found to be converted to other land use types between 2002 and 2012. This study was conducted in Nong Khai and Bueng Kan province of northeast Thailand, where massive rubber plantation is going on prompting significant amount of land-use change, with the objective of investigating how land-use changes will affect on food availability in future. We analyzed land-use changes of the past and simulated future land uses using GIS and Landsat Thematic Mapper Data. The most obvious change was the decrease in paddy field and an increase in rubber plantation. This eventually leads to decreased paddy production affecting food supply of farm households. The land use projections for 2032 were done for three scenarios using Dyna-CLUE model. Unlike business as usual scenario, which will further decrease the paddy area, other scenarios with different land use policies if implemented will help protect paddy areas and thus achieving higher food production locally. The lack of implementation of proper spatial policies will lead to a further loss of paddy areas at macro level. The smallholder farmers may be highly vulnerable to land use-change and experience significant food crop losses, food insecurity and income loss when they change the land to rubber and there is market failure.
Key words: land-use change; rubber plantation; scenarios; Dyna-CLUE; northeast Thailand
SAKAYAROTE Kanda , SHRESTHA Rajendra P. . Simulating land use for protecting food crop areas in northeast Thailand using GIS and Dyna-CLUE[J]. Journal of Geographical Sciences, 2019 , 29(5) : 803 -817 . DOI: 10.1007/s11442-019-1629-7
Figure 1 Location map of the study area showing major land uses in 2012 |
Table 1 Change of rubber area and paddy area in Thailand from 2007 to 2016 |
Year | Para rubber area (ha) | Paddy area (ha) | ||
---|---|---|---|---|
Of the country | Of the study area | Of the country | Of the study area | |
2007 | 2,888,800 | 85,043 | 9,181,760 | 143,183 |
2008 | 3,009,440 | 102,052 | 9,187,520 | 141,119 |
2009 | 3,140,000 | 96,259 | 9,199,520 | 111,788 |
2010 | 3,242,240 | 134,126 | 10,331,840 | 199,671 |
2011 | 3,386,400 | 136,310 | 10,448,640 | 180,584 |
2012 | 3,597,120 | 136,432 | 10,392,160 | 179,838 |
2013 | 3,711,040 | 162,506 | 9,932,800 | 156,297 |
2014 | 3,738,999 | 162,666 | 9,726,560 | 153,909 |
2015 | 3,733,108 | 162,706 | 9,290,080 | 152,849 |
2016 | 3,734,400 | 166,934 | 9,348,320 | 156,716 |
Source: OAE, Agricultural Statistics of Thailand 2007-2016 |
Table 2 Land suitability for rice and rubber |
Land characteristics | Limiting values for rice | Limiting values for rubber | ||||||
---|---|---|---|---|---|---|---|---|
S1 | S2 | S3 | N | S1 | S2 | S3 | N | |
Annual mean temperature | >24 | 21-24 | 18-21 | <18 | 26-28 | 23-25, 29-34 | 20-22 | >34, <20 |
(℃) | ||||||||
Elevation (m) | 0-600 | 600-1200 | 1200-1800 | >1800 | < 250 | 250-400 | 400-500 | >500 |
75% probability rainfall (mm) | >1300 | 900-1300 | 500-900 | <500 | 1500-2000 | 2000-3000 | 1500-2000 | 2000-3000 |
Soil drainage class | Poorly drained | Imperfectly drained | Moderately well drained | Excessively drained | Well to moderately well | Moderately well to somewhat poorly | Poorly | Very poorly |
Soil texture | C, ZC, ZCL, L | SC, SCL, ZL, Z | SL | S, LS | SSiL, SiL | SCL, CL | LS | - |
Soil depth (cm) | >80 | 60-80 | 40-60 | <40 | >150 | 100-150 | 50-100 | <50 |
pH of flooded soil | 7-Jun | 5 6, 7-8 | 4.5-5, 8-8.5 | <4.5 | 5.0-7.3 | 7.3-8.0 | 3.5-4.0 | >8.0, <3.5 |
ECe (mS cm-1) | <3 | 3-5 | 5-7 | >7 | >10 | 3-10 | <3 | - |
Slope angle (degrees) | <1 | 1-2 | 2-6 | >6 | <20 | 20-40 | 40-60 | >60 |
Sources: LDD (1996) |
Table 3 Land demand characteristics of three scenarios |
Land use type | Percentage change from the current trend | ||
---|---|---|---|
Scenario 1: Business as usual (A continuation of the land transformation rates for the past ten year is assumed) | Scenario 2: Policy to protect high potential land for paddy cultivation (2032) | Scenario 3: Policy to protect high and moderate potential land for paddy cultivation (2032) | |
Field crops Forest Paddy Perennial crops Water Rubber Orchard Urban Others | -2.69 -1.31 -24.09 5.71 -2.19 18.72 0.03 0.15 5.68 | -4.16 1.59 -17.46 2.84 -2.22 14.23 -0.33 0.60 4.91 | -2.80 -2.62 -11.87 3.76 -1.88 9.96 0.08 0.16 5.21 |
Spatial policy | No restriction on any land use class | Retention of natural forest and high potential land for paddy field (S1) | Restriction on high (S1) and moderate (S2) potential land for paddy field, forest, and unsuitable land for rubber |
Table 4 Location factors used for future land use simulation |
Factors | Variable scale | Minimum | Maximum |
---|---|---|---|
Elevation (m) | Continuous | 158 | 245 |
Distance to main road (km) | Continuous | 0 | 6.08 |
Distance to water resource (km) | Continuous | 0 | 6.9 |
Land with property right | Binary | 0= non-property right | 1= property right |
Distance to district (km) | Continuous | 30.9 | 41 |
Slope (%) | Continuous | 0 | 21.39 |
Land suitability for rice | Binary | 0= not suitable for rice | 1= suitable for rice |
Land suitability for rubber | Binary | 0= not suitable for rubber | 1= suitable for rubber |
Figure 2 Land uses of the study area in different years (1993 to 2012) |
Table 5 Land use transition matrix of the study area from 2002 to 2012 (ha) |
Land use type | Field crops | Forest | Paddy field | Perennial crops | Water | Rubber | Orchard | Urban area | Other | Total (2002) |
---|---|---|---|---|---|---|---|---|---|---|
Field crops | 6089 | 4676 | 4163 | 1448 | 1221 | 8542 | 1374 | 1177 | 2,580 | 31,269 |
Forest | 1466 | 38,151 | 12,988 | 2069 | 714 | 6479 | 647 | 1396 | 3,551 | 67,462 |
Paddy field | 7533 | 9605 | 153,439 | 16,149 | 10,191 | 53,919 | 919 | 10,755 | 22,650 | 285,160 |
Perennial crops | 1919 | 2548 | 6449 | 8788 | 839 | 17,417 | 120 | 887 | 2,908 | 41,875 |
Water | 851 | 537 | 7425 | 1176 | 8452 | 1786 | 143 | 1185 | 5,840 | 27,395 |
Rubber | 6289 | 2192 | 25,915 | 13,747 | 2933 | 91,684 | 485 | 4479 | 12,164 | 159,888 |
Orchard | 314 | 251 | 2301 | 332 | 283 | 641 | 145 | 475 | 604 | 5,347 |
Urban area | - | - | - | - | - | - | - | 11,352 | 839 | 12,193 |
Other | 9725 | 7641 | 22,172 | 7089 | 5414 | 20,707 | 303 | 3639 | 17,309 | 93,998 |
Total (Year 2012) | 34,187 | 65,600 | 234,852 | 50,798 | 30,047 | 201,174 | 4136 | 35,345 | 68,445 | 724,590 |
Table 6 Rubber cultivation in different land suitability levels for rubber in the study area |
Land suitability for rubber | Total study area (ha) | Rubber cultivation under different land suitability levels (ha) | |||
---|---|---|---|---|---|
2002 | 2006 | 2009 | 2012 | ||
Highly suitable | 31,666 | 1669 | 3735 | 4356 | 6019 |
Moderately suitable | 49,068 | 6484 | 14,825 | 21,456 | 26,497 |
Marginally suitable | 353,072 | 64,623 | 87,246 | 106,691 | 119,868 |
Unsuitable area and forests | 290,785 | 7075 | 9476 | 3720 | 4778 |
Total | 724,590 | 79,850 | 115,282 | 136,223 | 157,164 |
Table 7 Factors influencing land use type |
Variable | Field crops | Forest | Paddy | Perennial crops | Water | Rubber | Orchard | Urban | Other |
---|---|---|---|---|---|---|---|---|---|
Constant | -0.839 | -2.210 | 0.581 | -3.034 | -3.710 | -1.700 | -5.063 | -3.027 | -3.198 |
β values | |||||||||
Elevation | -0.012 | 0.008 | n.s. | 0.008 | -0.003 | 0.004 | -5.021 | n.s. | -0.003 |
Distance to main road | n.s. | -0.002 | 0.001 | n.s. | -0.001 | 0.008 | n.s. | 0.002 | -0.001 |
Distance to water resource | 0.001 | -0.002 | -0.078 | -0.004 | -0.002 | n.s. | -0.001 | 0.043 | n.s. |
Hill shade | -0.006 | -0.801 | -0.01 | 0.003 | 0.003 | n.s. | 0.005 | n.s. | 0.004 |
Land property rights | 0.597 | -0.266 | -0.10 | 0.357 | -0.397 | 0.415 | 0.466 | 0.289 | -0.296 |
Distance to town | n.s. | -0.002 | 0.01 | -0.001 | 0.001 | -0.003 | 0.002 | 0.005 | 0.013 |
Slope | -0.0014 | -0.348 | -0.06 | n.s. | -0.002 | 0.176 | 0.002 | -0.095 | 0.098 |
Land suitability for rice | -0.003 | 0.192 | 0.01 | -0.288 | -1.17 | -0.232 | n.s. | -0.004 | -0.003 |
Land suitability for rubber | 0.216 | 0.485 | -0.176 | -0.154 | -1.89 | 0.171 | 0.495 | -0.218 | 0.002 |
Note: n.s. = not significant at 0.01 level |
Figure 3 Future land use maps of the study area under three scenarios |
Table 8 Area coverage under projected land use types by 2032 in the study area |
Land use | Scenario 1: Business as usual (2032) | Scenario 2: Policy to protect high land potential for paddy cultivation (2032) | Scenario 3: Policy to protect high and moderate land potential for paddy cultivation (2032) | |||
---|---|---|---|---|---|---|
ha | % | ha | % | ha | % | |
Field crops | 22,236 | 3.07 | 22,810 | 3.15 | 32,642 | 4.50 |
Forest | 45,531 | 6.28 | 83,084 | 11.47 | 52,649 | 7.27 |
Paddy | 225,785 | 31.16 | 275,452 | 38.01 | 315,979 | 43.61 |
Perennial crops | 41,532 | 5.73 | 31,834 | 4.39 | 38,509 | 5.31 |
Water | 31,502 | 4.35 | 31,636 | 4.37 | 34,078 | 4.70 |
Rubber | 255,118 | 35.21 | 183,539 | 25.33 | 152,577 | 21.06 |
Orchard | 1,813 | 0.25 | 1,779 | 0.25 | 4,749 | 0.66 |
Urban | 41,036. | 5.66 | 35,964 | 4.96 | 32,771 | 4.52 |
Other | 60,038 | 8.29 | 58,491 | 8.07 | 60,635 | 8.37 |
Total | 724,590 | 100 | 724,590 | 100 | 724,590 | 100 |
Table 9 Household rice production and income |
Variable | Average | Conversion (Change of paddy to rubber) (n=100) | Non-conversion (No change) (n=111) | T-value |
---|---|---|---|---|
1. Rice production (ton/household /year) | 3.14 | 1.83 (0.27) | 4.54 (0.37) | ** |
2. Rice productivity (ton/ha/year) | 1.77 | 1.74 (0.39) | 1.81 (0.25) | * |
3. Rice purchase for household consumption (ton/household/year) | 0.34 | 1.02 (0.79) | 0.02 (0.013) | ** |
4. Rice expenditure (USD102/household/year) | 1.57 | 2.26 (1.69) | 0.03 (0.027) | ** |
5. Household rice consumption (ton/ household/year) | 0.81 | 0.74 (0.25) | 0.86 (0.31) | * |
6. Rice consumption (kg/person) | 229 | 224 (111) | 238 (172) | * |
7. Household member (person) | 4.02 | 3.59 (0.91) | 4.45 (1.02) | * |
8. Net income of the household (103USD/year) | 7.58 | 8.89 (4.66) | 8.91 (4.89) | - |
9. Average yield of rubber (ton/ha/year) | 1.17 | 1.14 (0.24) | 1.18 (0.33) | * |
10. Income from rubber production (USD/ha/year)1 | 1,495 | 1477 (341) | 1491 (376) | * |
11. Income from rice production (USD/ha/year) 1 | 713 | 708 (245) | 722 (159) | * |
*Significant at P<0.05; **Significant at P<0.01. Figure in parentheses is standard deviation; 1 USD = approx. 35 Thai baht 1Rubber price (Cup lump) = 1.28 USD/kg; Rice price (Jasmine rice) = 0.40 USD/kg (in 2012) |
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
/
〈 | 〉 |