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.