Results and Discussion
Results and discussion
Assessment of historical LU transitions and simulation
Temporal LU analyses were carried out using RS data of 1985, 2005, and 2019 through a supervised classifier based on the Gaussian maximum likelihood algorithm. Fig. 6 depicts LU in Karnataka, indicating that the Western Ghats districts have a higher forest cover. The state has witnessed large-scale transitions with the decline of forests from 21% (1985) to 15% (2019). LU transitions gained impetus with globalization and consequent industrialization, urbanization, and infrastructure development during the post-1990s. Unplanned developmental activities have resulted in the loss of forest cover in the Western Ghats districts of Belgaum, Uttara Kannada, Dakshina Kannada, Shimoga, Kodagu, and Chikmagalur. Large tracts of forests were lost due to developmental activities such as constructing dams and reservoirs, land conversion for built-up areas, creating SEZs and townships. The abrupt LU changes are evident from the decline of productive agricultural lands near cities such as Bengaluru, Mysore, Hubli-Dharwad, and Shimoga. The forest cover is restricted to major conservation reserves such as protected areas, national parks, and wildlife sanctuaries. The category-wise LU dynamics are presented in Table 4, highlighting an increase in the built-up cover from 0.47% to 3% and horticulture area from 8.8% (1985) to 11.1% (2019). The decline in the spatial extent of agricultural lands, forests, and lakes highlights the need for sustainable LU policies to arrest deforestation and abrupt land conversions.
Table 4. Ecosystem Extent – Karnataka state (extent in km2 and percentages) – based on temporal RS data analyses
LU Categories |
Built-up |
Cropland |
Horticulture |
Fallow land |
Evergreen Forest |
Moist Deciduous |
Dry Deciduous |
Scrub \Grass lands |
Water |
Total Area |
|
Year |
Units |
||||||||||
1985 |
(km2) |
904 |
128468 |
16790 |
1678 |
14293 |
10960 |
7622 |
6733 |
4344 |
191791 |
% |
0.5 |
67 |
8.8 |
0.9 |
7.5 |
5.7 |
4 |
3.5 |
2.3 |
100 |
|
2005 |
(km2) |
2666 |
127196 |
20209 |
1185 |
12445 |
9900 |
7410 |
5604 |
5177 |
191791 |
% |
1.4 |
66.3 |
10.5 |
0.6 |
6.5 |
5.2 |
3.9 |
2.9 |
2.7 |
21.13 |
|
2019 |
(km2) |
5748 |
127962 |
21325 |
2854 |
10888 |
7892 |
4281 |
4907 |
5934 |
191791 |
% |
3 |
66.7 |
11.1 |
1.5 |
5.7 |
4.1 |
2.2 |
2.6 |
3.1 |
17.68 |
|
LU Changes during 1985 to 2019 |
|||||||||||
1985 |
(km2) |
904 |
128468 |
16790 |
1678 |
14293 |
10960 |
7622 |
6733 |
4344 |
191791 |
2019 |
(km2) |
5748 |
127962 |
21325 |
2854 |
10888 |
7892 |
4281 |
4907 |
5934 |
191791 |
Net change of extent (during 1985 to 2019) |
|||||||||||
Extent |
(km2) |
4844 |
-505 |
4536 |
1175 |
-3405 |
-3068 |
-3341 |
-1826 |
1590 |
0 |
% |
535.8 |
-0.4 |
27 |
70 |
-23.8 |
-28 |
-43.8 |
-27.1 |
36.6 |
The natural forests show a decline, evident from the decrease in evergreen forests from 7.5% (1985) to 5.7% (2019), moist deciduous forests from 5.7% (1985) to 4.1% (2019), and dry deciduous forests from 4.0% (1985) to 2.2 % (2019). An accuracy assessment of the classified RS data was done through the computation of category-wise accuracies and Kappa statistics. Table 5 lists the category-wise accuracies, which indicate that the overall accuracy ranges from 86.84% (1985), 90.08% (2005) to 89.39% (2019), and Kappa ranges from 0.85 (1985) to 0.88 (2019). The LU transitions across ecosystems for Karnataka are provided in Table 6. New urban agglomerations were noticed across cities and major towns such as Bengaluru, Mangalore, Hubli, Hassan, and Mysuru. Large-scale monoculture plantations of eucalyptus, rubber, acacia, teak, and areca nut have increased and now cover 12% of the state. These abrupt changes result in an imbalance of ecosystem services, affecting the hydrologic regime and availability of natural resources. Tier-1 cities such as Bengaluru and Mangalore, and tier-2 cities such as Mysuru, Hubli-Dharwad, and Belgaum are experiencing loss of agricultural areas in the sub-urban regions, with new layouts and satellite towns.
Table 5. Accuracy assessment
Year |
1985 |
2005 |
2019 |
|||
Category |
PA (%) |
UA (%) |
PA (%) |
UA (%) |
PA (%) |
UA (%) |
Built-up |
93.97 |
95.19 |
92.58 |
88.01 |
98.99 |
96.27 |
Cropland |
94.01 |
94.71 |
97.85 |
100.00 |
93.97 |
93.15 |
Horticulture |
63.98 |
69.40 |
90.34 |
81.38 |
93.13 |
88.27 |
Fallow land |
84.58 |
76.09 |
98.14 |
100.00 |
91.34 |
90.88 |
Evergreen Forest |
91.77 |
47.95 |
88.53 |
52.43 |
79.46 |
100.00 |
Moist Deciduous Forest |
94.89 |
75.42 |
67.05 |
86.23 |
66.95 |
83.99 |
Dry Deciduous Forest |
69.93 |
69.70 |
77.77 |
79.76 |
93.80 |
87.63 |
Scrub_Grass lands |
66.95 |
64.54 |
80.24 |
92.82 |
39.21 |
81.23 |
Water |
91.45 |
79.95 |
91.60 |
100.00 |
97.62 |
86.39 |
Overall Accuracy (%) |
86.84 |
90.08 |
89.39 |
|||
Kappa |
0.85 |
0.87 |
0.88 |
Table 6. Transitions across LU categories during 1985 to 2019 – Karnataka state (extent in km2 and percentages)
2019 LU Categories 1985 |
Units |
Built-up |
Cropland |
Horticulture |
Fallow land |
Evergreen Forest |
Moist deciduous |
Dry Deciduous Forest |
Scrub / Grass lands |
Water |
Total (Opening Stock) 1985 |
Built-up |
(km2) |
859 |
20 |
15 |
2 |
1 |
3 |
1 |
2 |
2 |
904 |
% |
95.0 |
2.2 |
1.7 |
0.2 |
0.2 |
0.4 |
0.1 |
0.2 |
0.3 |
||
Cropland |
(km2) |
3169 |
114150 |
4050 |
4529 |
69 |
119 |
117 |
254 |
2008 |
128468 |
% |
2.5 |
88.9 |
3.2 |
3.5 |
0.1 |
0.1 |
0.1 |
0.2 |
1.6 |
||
Horticulture |
(km2) |
775 |
3344 |
11661 |
185 |
243 |
258 |
93 |
108 |
122 |
16790 |
% |
4.6 |
19.9 |
69.5 |
1.1 |
1.4 |
1.5 |
0.6 |
0.6 |
0.7 |
||
Fallow land |
(km2) |
50 |
843 |
32 |
710 |
1 |
2 |
12 |
13 |
15 |
1678 |
% |
3.0 |
50.2 |
1.9 |
42.3 |
0.1 |
0.1 |
0.7 |
0.8 |
0.9 |
||
Evergreen Forest |
(km2) |
175 |
416 |
2372 |
151 |
9097 |
1402 |
211 |
286 |
182 |
14293 |
% |
1.2 |
2.9 |
16.6 |
1.1 |
63.6 |
9.8 |
1.5 |
2.0 |
1.3 |
||
Moist Deciduous Forest |
(km2) |
190 |
1973 |
1648 |
485 |
388 |
5581 |
470 |
164 |
61 |
10960 |
% |
1.7 |
18.0 |
15.0 |
4.4 |
0.4 |
54.1 |
4.3 |
1.5 |
0.6 |
||
Dry Deciduous Forest |
(km2) |
85 |
3374 |
779 |
458 |
68 |
419 |
2306 |
87 |
47 |
7622 |
% |
1.1 |
44.3 |
10.2 |
6.0 |
0.9 |
5.5 |
30.3 |
1.1 |
0.6 |
||
Scrub/Grass lands |
(km2) |
327 |
3056 |
701 |
438 |
136 |
113 |
72 |
1787 |
103 |
6733 |
% |
4.9 |
45.4 |
10.4 |
6.5 |
2.0 |
1.7 |
1.1 |
26.5 |
1.5 |
||
Water |
(km2) |
93 |
734 |
114 |
35 |
14 |
16 |
7 |
9 |
3321 |
4344 |
% |
2.1 |
16.9 |
2.6 |
0.8 |
0.3 |
0.4 |
0.2 |
0.2 |
76.5 |
||
Closing Stock, 2019 |
(km2) |
5725 |
127910 |
21371 |
6994 |
10018 |
7914 |
3288 |
2710 |
5862 |
191791 |
% |
3.0 |
66.7 |
11.1 |
3.6 |
5.2 |
4.1 |
1.7 |
1.4 |
3.1 |
100 |
Table 7 summarises the ecosystem extent account for LULC in Karnataka state. The LU transition from 1985 to 2009 was accounted to understand the area of change and probability of change in each LU type from time t1 to time t2 using the Markovian process. The accuracy for the simulated 2019 LU map was assessed by accounting for the consistency of land quantity and spatial position similarity at the pixel level (Table 8). Kappa statistics were computed to evaluate the consistency of land quantity simulated in comparison with the actual LU. The similarity of the spatial position of simulated LU with actual with respect to a certain pixel was also assessed. Fig. 7 illustrates the simulated LUs of 2019 with area details. The simulated LU depicts the expansion of existing built-up cover in peri-urban regions, through LU conversion from agriculture, based on the transition accounted for from 2005 (as per CA neighbourhood influence). The accuracy of the model was found to be 95%; the simulated LU and actual LU appear visually consistent, but spatial differences exist.
Table 7. Ecosystem extent account for LULC in Karnataka state, India
Level-1 |
Level-2 |
Karnataka |
||||
Opening Stock 1985 |
Additions to Stock |
Reduction in Stock |
Closing Stock 2019 |
Net change (in%) during 1985 to 2019 |
||
Built-up land |
Built-up |
904 |
4866 |
45 |
5725 |
533.1 |
Urban |
||||||
Sub-Total 1 |
904 |
4866 |
45 |
5725 |
533.1 |
|
Agricultural Land |
Horticulture |
16790 |
9711 |
5129 |
21371 |
27.3 |
Cropland |
128468 |
13760 |
14317 |
127910 |
-0.4 |
|
Fallow Land |
1678 |
6284 |
968 |
6994 |
316.7 |
|
Sub-Total 2 |
146936 |
29754 |
20414 |
156275 |
6.4 |
|
Forests |
Evergreen/Semi-Evergreen |
14293 |
921 |
5196 |
10018 |
-29.9 |
Moist Deciduous |
10960 |
2333 |
5379 |
7914 |
-27.8 |
|
Dry Deciduous |
7622 |
981 |
5316 |
3288 |
-56.9 |
|
Scrub Forest |
6733 |
922 |
4946 |
2710 |
-59.8 |
|
Forest Plantation |
||||||
Swamp/Mangroves |
||||||
Sub-Total 4 |
39607 |
5158 |
20836 |
23929 |
-39.6 |
|
Grass / Grazing |
Grass / Grazing |
|||||
Sub-Total 5 |
||||||
Wetlands / Water bodies |
Coastal Wetland |
|||||
River/stream/canals |
||||||
Waterbodies |
4344 |
2541 |
1023 |
5862 |
35.0 |
|
Sub-Total 6 |
4344 |
2541 |
1023 |
5862 |
35.0 |
|
Grand Total (Sq. Km) |
191791 |
42319 |
42319 |
191791 |
Table 8. Accuracy of the simulation
Information of Allocation |
No [n] |
Medium [m] |
Perfect [p] |
Kappa |
Perfect [P(x)] |
P(n)=0.3706 |
P(m)= 0.9541 |
P(p)=1 |
Kno=0.9003 |
Perfect Stratum [K(x)] |
K(n)=0.3706 |
K(m)=0.9541 |
K(p)=1 |
Klocation=0.9273 |
Medium Grid [M(x)] |
M(n)=0.3559 |
M(m)=0.9102 |
M(p)=0.9009 |
KlocationStrata=0.9273 |
Medium Stratum [H(x)] |
H(n)=0.1 |
H(m)=0.3502 |
H(p)=0.3598 |
Kstandard=0.8619 |
No [N(x)] |
N(n)=0.1 |
N(m)=0.3502 |
N(p)=0.3598 |
Fig. 6. LU changes from 1985 to 2019 in Karnataka
Fig. 7. Simulated land uses in 2019
Scenario-based prediction and validation
The distance maps for various urban growth factors (agents) were prepared (Fig. 8), including slope maps, distance to road, proposed railway projects, and industries. Five different constraints maps were prepared for the chosen scenarios. The site suitability maps were generated to define the suitability of each pixel for transition to any LU type. Each pixel in the suitability maps has a value ranging from 0 to 255, with 0 representing unsuitable and 255 representing the highly suitable area for a particular LU type (Pontius and Malanson 2005).
The distance influence of each agent on specific LUs was evaluated; this helped in effectively deriving the weights, unlike traditional weighing. Fig. 4 shows the influence of factors such as roads and industries on the built-up area. Industries have influence up to 2000 m from urban areas and moderate influence up to 10,000 m; the influence gradually decreases after that. Roads have shown significant influence up to 500 m, with a reduction in influence after 1000 m. Fig. 9 shows the factors that influence forest cover LU, which helped in understanding the distance influence of each factor. Industries show high conversion of forests to other LU up to 7000 m, and after that, the influence reduced; roads showed stronger influence up to 2000 m.
The AHP was used to assign weights and evaluated its consistency by computing the consistency ratio. Scenario-wise weights were assigned for each factor and are given in Annexure 1 (Tables a–e). A consistency ratio of <0.1 was achieved for each scenario, and the eigenvector weights for individual factors considered for the scenario are listed in Tables a-e (Annexure). The site suitability of LU transition was accounted for using MCE and was provided as an input for CA-based prediction. The CA has been used for computing projected LUs under five different scenarios based on varied inputs (Fig. 10).
BAU can be considered as a base case that helps in comparing or differentiating with alternative growth strategies. It depicts the likely increase in built-up area to 11.5% (2033) from 3% in 2019 (Table 9, Fig. 11). The forest cover and the agricultural regions would lose significant tract in case of likely expansion of built-up cover. The dry and moist deciduous forest cover has witnessed a higher change with a loss of 6.5% forest cover. The existing infrastructure and built-up cover will transform deciduous forests and agriculture areas along highways across the State. The cities such as Bengaluru, Mangalore, and Dharwad (tier-1) depicted compact growth, whereas Mysuru, Belgaum, Hassan, and Tumkur (tier-2) indicated peri-urban development.
Table 9. Likely LU in 2033 for various scenarios
Land use categories |
BAU_2033 |
ALT_2033 |
RFP_2033 |
AF_2033 |
SDP_2033 |
|||||
Area (km2) |
% |
Area (km2) |
% |
Area (km2) |
% |
Area (km2) |
% |
Area (km2) |
% |
|
Built-up |
22,015 |
11.5 |
28,100 |
14.7 |
21,703 |
11.2 |
21,462 |
11.2 |
20,720 |
10.8 |
Cropland |
1,17,646 |
61.3 |
1,13,309 |
59.1 |
1,18,946 |
62.0 |
1,18,586 |
61.8 |
1,16,672 |
60.8 |
Plantation |
21,777 |
11.4 |
20,667 |
10.8 |
21,984 |
11.5 |
21,692 |
11.3 |
21,521 |
11.2 |
Horticulture |
6,688 |
3.5 |
6,688 |
3.5 |
6,851 |
3.6 |
6,688 |
3.5 |
6,851 |
3.6 |
Evergreen Forest |
9,209 |
4.8 |
9,152 |
4.8 |
9,328 |
4.9 |
9,591 |
5.0 |
9,725 |
5.1 |
Moist Deciduous Forest |
5,774 |
3.0 |
5,721 |
2.9 |
5,414 |
2.8 |
5,752 |
3.0 |
6,670 |
3.5 |
Dry Deciduous |
1,524 |
0.8 |
1,512 |
0.8 |
1,246 |
0.7 |
1,569 |
0.8 |
2,369 |
1.2 |
Scrub_Grass |
1,160 |
0.6 |
1,144 |
0.6 |
915 |
0.5 |
979 |
0.5 |
1,967 |
1.0 |
Water |
5,997 |
3.1 |
5,497 |
2.8 |
5,404 |
2.8 |
5,471 |
2.9 |
5,298 |
2.8 |
Total |
1,91,791 |
Upgrading urban infrastructure (roads, commercial establishments, information networks, industrial units, and sports/recreation centers) will result in a built-up cover of 15% compared to BAU. The ALT scenario will increase paved surfaces by 2033 due to the development of towns and compact urban centers (consisting of high-density residential and industrial layouts). The BAU and ALT scenario predictions indicate a significant reduction in forest area and agricultural LU. The ALT scenario shows that forest cover would likely be only 9% compared to 15% in 2019 (Table 9, Fig. 11). The agricultural area showed a loss of 2%, which might affect food availability in the region. The regions where dry crops are grown will witness higher transformation, which signifies the external driver’s role in LU alterations. The ALT scenario indicates higher LU conversion due to transportation and other projects with compact, peri-urban, and mixed LU developments; this necessitates effective LU policies.
The RFP scenario visualizes protecting reserve forests without any alterations for unplanned developmental activities and expansion of earlier LUs. The regions covered under reserve forest protection do not favour development and other LU conversions. The districts that are part of the Western Ghats region might not affect the BAU and ALT scenarios. The regions with dry deciduous forest cover and not falling under reserve forest area protection will experience abrupt growth due to transportation corridors and urban clusters. Unlike other scenarios, urban areas of 2019 are considered as a factor influencing LU changes. The RFP scenario shows a likely built-up cover of 11%, retaining a good amount of evergreen forest cover of 5% (Table 9). As depicted in Fig. 11, the AF scenario shows likely improvement in forest cover due to plantation activities and protection of degraded forests from uncontrolled LU conversion. Here, the forest cover is 9.3%, a slight increase as compared to BAU and ALT. The built-up cover will limit existing cities and towns by over 11%. The new plantation activities will compensate for the abrupt LU conversion, especially in moist, dry deciduous cover regions. Districts such as Bidar, Kolar, and Koppala did not show any improvement due to the availability of the least forest land area.
The SDP scenario signifies stringent policy implementation initiatives with minor disturbances, and the simulated LU using this scenario is depicted in Fig. 11. The forest cover will remain 11% compared to 17% in 2019 (Table 9), the least probable loss among all the scenarios. The existing built-up cover will remain in the tier-1 and tier-2 cities; the least increase was shown in peri-urban areas. Districts such as Bidar, Gulbarga, Bagalkot, and Kolar might experience higher LU conversion rates due to existing infrastructure and be devoid of any forest policy. Overall, the RFP, AF, and SDP scenarios favour the protection of evergreen and moist deciduous forest categories. The BAU and ALT scenarios show a merging of urban clusters, indicating that smaller clusters of paved areas are no longer available and tend to agglomerate with bigger ones to form urban centers, satellite towns, and dispersed growth at the periphery of tier-2 city limits. There might be new constructions in the rural hinterland, resulting in loss of evergreen, moist deciduous cover and agricultural areas. Categories such as plantation, water, and open field had similar values across all scenarios. The scenarios will directly help decision-makers and planners to mitigate environmental impact with strict growth rules; otherwise, the conversion of primary, secondary forest/ moist/dry deciduous forest into urban areas will be favoured (BAU and ALT) can induce ecological imbalances. Another important observation in terms of statistics was that slope steepness restriction played an important role in restricting urban settlements from further infiltrating into forested areas. Collectively, sustainable development seems feasible with scenarios 3–5, without altering the existing edges of agriculture, forest, water bodies, and other protected areas or conservation zones.
The limitations of the current study lie in the factors considered, which vary based on developmental issues and regional dynamics. The model results may vary depending on the assumptions, modeling strategies, computational approaches, and data. The driving factors are generally region-specific, and the likely LU changes are expressed based on local variations in geographic and socio-economic conditions. The modeled scenarios applied at the micro (e.g., cities) level might show variations in results than the entire study region (state). The challenges to be addressed in LU transition models are: (a) lack of framework for systematic assessment of LULC changes and hence, its validation considering the decision-making process and authenticity of interactions between various agents in a system; (b) resolution of data, availability, and trade-offs between different resolutions. LU change is dependent on the land management policy framed by local, regional, national, and global agencies. Therefore, there is a need to integrate these factors along with human-environment interactions to achieve desirable results.
Fig. 8. Distance influence on LU change-built-up
Fig. 9. Distance influence on LU change-forest
Fig. 10. Projected LU under BAU, ALT, RFP, AF, and SDP scenarios
Fig. 11. Projected LU (2033) of policy scenarios
The modeled scenarios aid as a base for the decision-making process, to balance LU changes associated with developmental activities and conservation. The abrupt changes induced by developmental activities such as creating industries, upgradation of infrastructure, and establishing new linear corridors in regions dominated by forests will enhance deforestation and degradation. The unscientific LU changes in forest landscapes result in ecosystem imbalance and increased carbon emissions in the atmosphere (Ramachandra and Bharath 2019a). Landscapes with increased paved surfaces tend to have considerably greater air temperature than their peri-urban/rural surroundings, altered microclimatic conditions, increased anthropogenic heat release, and increased land surface temperature, which further contributes to global warming (Yao et al., 2017). The buildings and road surfaces act as high heat storages, exposing humans to health-threatening heat. The increase in paved surfaces due to uncontrolled urbanization in Karnataka state can result in an urban heat island effect due to the conversion of latent heat flux into sensible heat flux, thereby threatening human well-being (Ramachandra and Uttam, 2009). Other consequences of unregulated paved surfaces in a landscape are: (a) increase in the energy consumption for cooling with enhanced land surface temperature (Ramachandra et al., 2017b), (ii) escalation in the carbon footprint (Ramachandra and Shwetmala, 2009), (iii) reduction in water availability (Ramachandra et al., 2020), (iv) alteration in the seasonal rainfall pattern (Buyantuyev and Wu, 2012), (v) increase in flood instances (Kumar et al., 2021), (vi) effect on air quality (Feizizadeh and Blaschke, 2013; Fuladlu and Altan, 2021), (vii) phenological changes (Allen et al., 2015), (ix) impact on biodiversity (Ramachandra et al., 2018), (x) reduction in the net primary productivity of vegetation (Jackson and Baker 2010; Bharath et al. 2013; Alavipanah et al. 2015), and finally (xi) vegetation die-off (Breshears et al. 2005; Zhou et al. 2016). In this regard, the Karnataka state authorities should focus on sustainable developmental planning with judicious resource usage and existing vegetation cover improvement. Banning highly-polluting industries, illegal mining, and unscientific land conversion, and regulating urbanization should be prioritized to mitigate impending climate change.