Geospatial analysis of forest fragmentation in Uttara Kannada District,India

Ramachandra T V,Bharath Setturuand Subash Chandran

Geospatial analysis of forest fragmentation in Uttara Kannada District,India

Ramachandra T V1,2,3*,Bharath Setturu1and Subash Chandran1

Background:Landscapes consist of heterogeneous interacting dynamic elements with complex ecological, economic and cultural attributes.These complex interactions help in the sustenance of natural resources through bio-geochemical and hydrological cycling.The ecosystem functions are altered with changes in the landscape structure.Fragmentation of large contiguous forests to small and isolated forest patches either by natural phenomena or anthropogenic activities leads to drastic changes in forest patch sizes,shape,connectivity and internal heterogeneity,which restrict the movement leading to inbreeding among Meta populations with extirpation of species.

Methods:Landscape dynamics are assessed through land use analysis by way of remote sensing data acquired at different time periods.Forest fragmentation is assessed at the pixel level through computation of two indicators, i.e.,Pf(the ratio of pixels that are forested to the total non-water pixels in the window)and Pff(the proportion of all adjacent(cardinal directions only)pixel pairs that include at least one forest pixel,for which both pixels are forested).

Results:Uttara Kannada District has the distinction of having the highest forest cover in Karnataka State,India.This region has been experiencing changes in its forest cover and consequent alterations in functional abilities of its ecosystem.Temporal land use analyses show the trend of deforestation,evident from the reduction of evergreensemi evergreen forest cover from 57.31%(1979)to 32.08%(2013)Forest fragmentation at the landscape level shows a decline of interior forests 64.42%(1979)to 25.62%(2013)and transition of non-forest categories such as crop land,plantations and built-up areas,amounting now to 47.29%.PCA prioritized geophysical and socio variables responsible for changes in the landscape structure at local levels.

Conclusion:Terrestrial forest ecosystems in Uttara Kannada District of Central Western Ghats have been experiencing threats due to deforestation with land use changes and fragmentation of contiguous forests,as is evident from the decline of interior forests and consequent increases in patch,transitional,edge and perforated forests.Interior or intact forest cover in this ecologically fragile region is now 25.62%.Considering the accelerating rates of forest fragmentation in recent times,the focus should be on reforestation and regeneration of natural vegetation to sustain food and water security and the livelihood of local populations.This requires innovation with holistic approaches in the management of forests by involving all local stakeholders to minimize the encroachment of forests,and improvements in regeneration.

Land use dynamics,Forest fragmentation,Interior forests,Biodiversity,Geophysical variables

Background

Any landscape is a mosaic of heterogeneous interacting dynamic elements,i.e.,manifestations of natural and anthropogenic processes.The structure of a landscape (size,shape and configuration)affects its functional aspects such as bio-geo chemical cycling and hydrologic regimes.The interactions among the landscape elements result in the flow of nutrients,minerals and energy, which contribute to the functioning of the landscape. Forest ecosystems constitute a key component of the global carbon cycle that account for over two-thirds of net primary production on land through photosynthesis converting solar energy into biomass(Roy et al.2001; MEA 2005;Ramachandra et al.2013).Forest ecosystems offer timber and non-timber forest products(NTFP), such as medicinal resources,fuel wood and as well provide recreational values(Kindstrand et al.2008). They aid as biodiversity repositories(Li et al.2009),restrain soil erosion(Nandy et al.2011),prevent landslides given that tree roots bind soil,regulate air humidity, temperature and mitigate global warming(Cabral et al. 2010)by absorbing 30%of fossil fuel CO2emissions (Pan et al.2011).The goods and services provided by forested landscapes are vital to the socioeconomic development of human populations(DeFries et al.2004)and their survival(Ramachandra et al.2013).At a large-scale land use,more recent land cover changes(LULC)are altering the ecosystem structure,affecting the goods and services of the ecosystem.These disturbances have resulted in fragmentation of forests with a mosaic of naturalpatchessurroundedbyotherlanduses (Ramachandra and Kumar 2011).A host of anthropogenic activities,such as tree logging,conversion of forest land to agriculture,intense agricultural practices,forest fire and unplanned infrastructural development have contributed to the disruption of the contiguity of forests in predominantly natural landscapes(Buskirk et al.2000; Boogaert et al.2004).An alteration in forest structure through fragmentation of forests has affected its functional abilities,as is evident from the decline in water yield,carbon sequestration potential and biodiversity. (Diaz et al.2006;Ramachandra and Kumar 2011).

Fragmentation processes involve alterations in the structure and composition of native forests through the division of contiguous forest into smaller noncontiguous fragments with a sharp increase in edges (Riitters et al.2002;Broadbent et al.2008).This will have detrimental effects such as disruption in bio-geo chemical cycling,nutrient and water cycling,ecological processes(Fahrig 2003;Holway 2005),easier access and further land use changes(Holway 2005; Bennett and Saunders 2010;Ramachandra et al.2012a). The edge effect may even destroy large trees within 300 m of the forest edge which are then replaced by densely spaced short-lived pioneers(Laurance 1999),resulting in the decline of forest biomass(Harper et al.2005).The negative impact of edges include alterations in plant and animal community compositions with subsequent changes in diversity(Cagnolo et al.2006),seed dispersion,pollination,predation,fire susceptibility,altered microclimate and increased carbon emissions(Laurance et al.2002). Consequences of edge effects(Carolina 1995;Asner et al. 2006)also include abiotic changes such as in its micro climate(Kampichler et al.2012)and water availability. Figure 1 explains the process of fragmentation with drivers.The earlier stages of forest fragmentation lead to changes in forest composition due to the transition of evergreen forests to a mix of semi evergreen and scrub forest patches leading to alterations in micro climate and habitat(Ramachandra et al.2013).During the succession stage,edges will become more prominent with higher light availability and loss of soil moisture.This will directly affect seed dispersal,seed viability,species distribution and propagation of invasivespecies(RamachandraandKumar2011). These cascading effects have the potential to disrupt seriously many basic ecological processes over large areas.The changes in forest structure will result in abrupteffectsanddisturbingtheproductivityof goodsandservicesfromthenaturalecosystem. Analyses of the extent of forest fragmentation will provide insights of complex dynamic interactions, which help in adopting appropriate location specificconservationmeasurestomitigatedisruptionsin ecological processes.

Quantification of forest fragmentation

Land use(LU)changes driven by anthropogenic activities alter the structure of a landscape,which adversely affect the functional aspects of an ecosystem.Land use patterns are the collective result of interactions among local geophysical indicators,such as elevation,slope and rainfall and agro-climatic indicators,demographic variation,market forces and related development policies (Munroe et al.2004;Ameztegui et al.2010).Numerous studies focusing on deforestation at the landscape level have explored the spatial patterns and interactions among the geophysical elements of the landscape and its dynamics(Nelson and Geoghegan 2002;Alix-Garcia et al.2005;Echeverria et al.2008;Ramachandra and Kumar 2011;Ramachandra et al.2013).These changes, measured at temporal scales,help in monitoring ecosystems and aid in the implementation of location specific mitigation measures.Spatial data acquired remotely through space-borne sensors(remote sensing data)at regular intervals help in assessing the temporal changes in spatial patterns(Ramachandra et al.2014a).Remote sensing data(RS)with geographic information systems (GIS)have made significant contributions to the examination of spatial-temporal patterns and processes of forest ecosystems(Nandy et al.2011;Ramachandra et al. 2012a)in criteria based decision-making and selection of the optimal alternative.The availability of remote sensing data,with improvements in resolutions(spatial, spectral and temporal resolutions),enables the creation of land use and land cover maps(Chen et al.2015),as well as that of innovative analytical techniques and helps to monitor changes in cost-effective ways(Bharath et al. 2012a).Changes in forested landscapes of the Western Ghats have been more sudden since the middle of the last century due to the impetus of industrialization policies as a consequence of globalization.Quantification of forest fragmentation in an ecologically fragile region,such as the Western Ghats will help in formulating appropriate mitigation measures towards the conservation of biodiversity.

Objectives

The objectives of this assessment of land use dynamics and forest fragmentation are to:

1.Understand the prevailing forest cover dynamics from 1979 to 2013 and evaluate the causal factors contributing to forest changes;

2.Explore the spatial-temporal patterns of forest fragmentation and quantification of the extent of fragmentation;

3.Evaluate the role of geophysical variables in forest cover changes and assess the effects of such variables on the patterns of forest loss at a temporal scale and

4.Suggest management options to mitigate forest loss and fragmentation to restore and sustain the ecosystem.

Method

Study area

Uttara Kannada(13.92 to 15.53°North and 74.09 to 75.09°East)is the fifth largest district of Karnataka State in India with abundant natural resources,perennial (water flows throughout the year)rivers,abundant flora, fauna and a long coastal line(about 140 km).The district consists of three distinct agro-climatic zones covering 11 taluks(local administrative divisions),i.e.,the coastal taluks Karwar,Ankola,Kumta,Honnavar and Bhatkal,the Sahyadri Interior taluks Supa,Yellapura, Sirsi and Siddapur and two eastern plains(connected to Deccan plateau)Haliyal and Mundgod.The total population of the district is 1,437,169(as per 2011 census) with a density of 140 persons per km2.The district covers about 15,055 km of road length,ranging from national and state highways to village roads.The four forest types(Fig.2)are dry deciduous,moist deciduous, evergreen to semi evergreen and coastal mixed deciduous forests,the result of its varied geographical features. Recent changes in the forest cover and associated habitat fragmentations are due to unplanned developmental activities.The implementation of hydroelectric and nuclear power projects with the consequent submergence of large tracts of evergreen forests and croplands has displaced thousands of families affecting their livelihood as well as the sustainability of natural resources.This necessitates a better understanding of forest dynamics withdisturbanceregimes,includingtheextentof fragmentations.

Quantification of spatial-temporal forest changes and extent of fragmentation

Land use(LU)changes in Uttara Kannada District were assessed using temporal remote sensing(RS)data with ancillary data and information obtained from the field. Figure 3 outlines the method adopted for the analysis of spatial patterns of forest changes.RS data used in the study are Landsat MSS(1979),TM(1989,1999),Landsat ETM+(2013)(downloaded from the archive-http://glovis. usgs.gov/)and Google Earth(http://earth.google.com). Slope maps were generated by using ASTER DEM(30 m) (http://gdem.ersdac.jspacesystems.or.jp/).Rainfall data was procuredfromtheDirectorateofEconomicsand Statistics,Government of Karnataka(http://des.kar.nic.in). Population data was procured from the Directorate ofCensus Operations,Karnataka(http://censuskarnataka. gov.in).Ancillary spatial data include cadastral revenue maps(1:6000),the Survey of India(SOI)topographic maps(1:50,000 and 1:250,000 scales)were procured from the Survey of India(http://surveyofindia.gov.in),etc. Ground control points(GCP’s)used to geo-rectify remote sensing data were digitized from the topographic maps and from field observations.Various forest cover types to classify RS data from the 1980’s were obtained from the vegetation map of South India at a scale of 1:250,000 (Pascal 1986).Training data required for classification and validation of the latest remote sensing data were compiled from field observations using a pre-calibrated hand-held GPS(Global Positioning System).GPS calibration was carried out from known ground coordinates(benchmarks, intersection of roads and other landmarks)as per the most commonly adopted protocol(Yuan 2008).FCC(false color composites)helped in identifying heterogeneous land use patches.Select land use patches(uniformly distributed in the study region covering at least 10%of the study region),were digitized and loaded to the GPS.Attribute information are collected(land use type,vegetation details)for these patches(training polygons)from the field.

Land use analysis and detection of change

Changes in land use during the last four decades were assessed using temporal RS data.This helped in determining the extent and causes of transitions in the landscape.Spatial data analyses involved(i)pre-processing, (ii)vegetation cover and(iii)land use analyses.Preprocessing involved geo-referencing,rectification and cropping of data pertaining to the study region.Georeferencing was carried out through ground control points collected from the field using pre calibrated GPS and from known points,such as road intersections,collected from geo-referenced topographic maps published by the Survey of India.The Landsat data of 1979,with a spatial resolution of 57.5 m×57.5 m,was resampled to 30 m comparable to the 1989–2013 data,which are 30 m×30 m(nominal resolution).The Landsat ETM+ bands of 2013 were corrected for the SLC-off(Scan Line Corrector failed)by using image enhancement techniques and nearest-neighbor interpolation.

Land use analysis involved(i)the generation of False Color Composites(FCC)of remote sensing data (bands–green,red and NIR),which helped in locating heterogeneous patches in the landscape,(ii)selection of training polygons by covering 15%of the study area(polygons are uniformly distributed over the entire study area),(iii)loading the co-ordinates of these training polygons into GPS,(vi)the collection of corresponding attribute data(land use types)for these polygons from the field,where the GPS helped in locating respective training polygons in the field,(iv) supplementing this information with Google Earth.In the end,60%of the training data was used for classification,with the balance going for validation or accuracy assessment(Ramachandra et al.2012c).The land use analysis was performed using a supervised classification technique based on the Gaussian maximum likelihood(GML)algorithm with training data (collected from the field using GPS).This evaluated quantitatively the variance and covariance of spectral response patterns of land uses based on a GML estimator(Atkinson and Lewis 2000;Ramachandra et al. 2012c).Land use classifications using temporal data was carried out through the open source program GRASS-Geographical Resources Analysis Support System(http://ces.iisc.ernet.in/grass).Accuracy assessments of the classified information have been performed(Lillesand and Keifer 1987;Liu et al.2007)to evaluate the performance of classifiers through the computation of an error matrix,kappa(κ)statistics and overall accuracies.

Forest fragmentation model implementation

Fragmentation of forests at the pixel level are estimated through the computation of Pf(the ratio of the number of pixels that are forested to the total number of nonwater pixels in the window)and Pff(the proportion of all adjacent(in cardinal directions)pixel pairs that include at least one forest pixel,for which both pixels are forested)as given in equations 1 and 2(Riitters et al. 2002;Kuèas et al.2011;Ramachandra and Kumar 2011).

Pffestimates the conditional probability that,given a pixel of forest its neighbor is also forest based,the proportion of all adjacent(cardinal directions)pixel pairs.Pfand Pffwere computed through a moving window of 5×5 pixels(Fig.4),given that the results of the model are scale-dependent and threshold dependent.A kernel of size 5×5 is selected for spatial data of 30 m, based on earlier work(Riitters et al.2002;Wickham et al.2007;Kuèasetal.2011;Prasadetal.2009; Ramachandra and Kumar 2011),which maintain a fair representation of the proportion(Pf)of pixels and also interior forest at an appropriate level(Riitters et al. 2000;2002;Wickham et al.2007;Kuèas et al.2011). Riitters et al.(2004)demonstrated the appropriateness of 5×5 kernels.Kernel sizes smaller than 5×5 have an effect of decreasing the average inter patch distance,indicating less fragmentation,even though disintegration of interior forests might be evident(Boogaert et al.2004; Lindenmayer et al.2008).Similarly,an increase of kernel size decreases the core area and this process may transform small core areas to form a discontinuous landmass (Ostapowicz et al.2008;Kuèas et al.2011).Depending on the indices Pfand Pff,spatial extent of forest fragmentation were mapped and details are presented in Table 1.Water bodies or river courses are considered non-fragmenting features,for these constitute natural corridors in a forested landscape,while anthropogenic landscape elements(such as buildings,roads,agricultural field and barren land)are drivers of forest fragmentation.

Principal Component Analysis(PCA)

Unplanned urbanization,geography,industrialization, government policies,economic reforms and population growth are the major factors driving landscape changes (Gong et al.2013).Socio-economic and bio-geophysical variables elucidate the role of anthropogenic forces in forest transitions,altering landscape structures and their composition(Ramachandra et al.2014a).These geophysical variables and socio-economic factors aid as drivers of land use changes in the landscape(Nelson et al.1999; Barbier 2001;Timar 2011;Cho et al.2015).

Table 1 Fragmentation components and their description

The primary objective of geospatial statistical analysis is to quantify the correlation between socio-economic and bio-geophysical variables with the fragmentation of forests over a set of non-water pixels in a landscape through multiple linear regression(MLR).But,MLR explains the output data as a weighted sum of individual correlations with the assumption that any individual variable/feature in the input set is linearly independent. A high level of interdependence among input variables imply that forest fragmentation correlations do not have an independent biological interpretation.To reduce redundancy in multivariate data,a decomposition of eigen values is used in Principal Component Analysis(PCA). It helps in noise reduction and also in prioritizing variables,responsible for the variation in spatial landscape processes(Morris et al.2009).Prioritizing these variablesbased on their role in land use changes will help to evolve management strategies and sustain landscape elements with their interactions(Todd and Kerr 2009).

PCA explores multivariate patterns of all forest patch types versus spatial processes based on correlation matrices and frames clusters and relates spatial heterogeneity among variables correlated with one another and responsive to multiple principal components(PC).These derived components are linear function of the original data set.The components derived are ordered by decreases in variance,i.e.,PC1 will have the largest variance among n components and PCn will have the smallest variance.After the extraction of the PCA eigenvalues,Scree plot is generated to prioritize the number of significant principal components(Jackson 1993).Finally,a series of ordinations were created to interpret the PCs visually and determine their fragmentation patterns.The sum of the eigenvalues is equal to the variance of the original data set,which preserves the original variation.The table of factor loadings explains the contribution of each variable to the derived components.A component is oriented towards that variable which has the maximum loading on it.

PCA with geo-physical variables aid in understanding the causal factors of environmental vulnerability.The positively correlated spatial parameters with geophysical and socio-economic variables serve as key descriptors of land use transitions(Salvati et al.2008).PCA analyzes the variance of variables and reorganizes it into a new set of uncorrelated independent components(principal components)equal to the number of original variables as linear combinations of the measured variables(Swan and Sandilands 1995).These combinations are based on weights(eigenvectors)and the loading for each item/ variable is the correlation between components,which serve to demarcate clusters of similar patterns(Colson et al.2011).PCA aided in prioritizing bio-geophysical and socio variables that act as agents of changes in vegetation cover.The variables used in PCA were normalized (Abdi and Williams 2010;Bell et al.2015)by computing Z-scores(Normalized(X)).

Results

Temporal land use analyses(Fig.5)show a decrease in the evergreen forest cover from 57.31%(1979)to 32.08%(2013).Enhanced agricultural activities are evident from the increase in agricultural land use(Table 2) from 10.02(1979)to 14.13%(2013)and areas under human habitation have increased during the last four decades from 0.95%(1979)to 3.07%(2013).In addition to these,various on-going unplanned developmental projects have contributed to the decline in forest cover. The increase in plantations of exotic species such as Acacia auriculiformis,Casuarina equisetifolia,various Eucalyptus spp.and Tectona grandis serve mainly to meet the demand by forest-based industries.These plantations now(2013)constitute 12.04%(Acacia auriculiformis,Casuarina equisetifolia,various Eucalyptus spp) and 6.60%(Tectona grandis)respectively in the district. Degradation of vegetation cover in the coastal zone had triggered a series of landslides in 2009 at 21 locations along a national highway,i.e.,NH-17,leading to the loss of property and human life(Ramachandra et al.2012b). The dry deciduous forest cover is has declined considerably from 2.83 to 0.96%(during 1973 to 2013)and is found mainly in the northeastern part of the district in Mundgod and Haliyal taluks.Accuracy of the classification ranges from 87 to 92%indicating consistent classification results.

Land use changes during the past four decades were estimated using temporal remote sensing data from 1979 to 2013 and are shown in Fig.5a,b,c,d.Naturalforest categories include moist deciduous,evergreen to semi evergreen,scrub/grass and dry deciduous forest types.Category-wiseforestcoverisextractedand depicted in Fig.6a-d,which helped in assessing the impacts of unplanned developmental activities.Figure 6a illustrates the natural forest cover in 1979(represented in green color)while Fig.6d shows forest cover in 2013. Table 3 provides the spatial extent of forest cover loss during 1979 to 2013.The natural forest cover of 780,778 ha in 1979 was reduced to 733,692 ha by 1989 due to a large scale hydro-electric project at Supa,as well as by the conversion of native forests to monoculture plantations for forest based industries and cultivation.The forest changes during this ten year period was moderate,i.e.,a 6.03%reduction,whereas rapid losses were taken place during the period from 1989 to 1999 (forest cover 63.93%)due to implementation of a large number of infrastructure projects,a series of hydroelectric projects,the construction of the Tattihalla and Bommanhalli reservoirs and other associated developments. The major disturbances noted in the forests of rural areas are due to unauthorized conversion of forest land into agriculture(encroachments),the diversion of forest lands for other purposes and logging of wood by forest based industries.Diversion of forestland under various schemes accounted for the loss of 66,443 ha(KFD 2010). The current area of forest cover is 542,475 ha(52.71%) and forest losses amounted to about 35,021 ha(17.54%) due to unauthorized land conversions and implementation of projects such as the Kodasalli reservoir,the Kadra dam,the Kaiga Nuclear Power plant,the Project Sea Bird,the Gerusoppa dam and others.

Table 2 Temporal land use changes from 1979 to 2013 and accuracy assessment

Table 3 Forest cover change from 1979 to 2013

Fragmentation of forests at the landscape level was assessed in order to understand the spatial-temporal patterns in forest degradation using temporal land use (1979 to 2013)information.Figures 7a,b,c,d and 8 depict the temporal pattern of the fragmentation process in Uttara Kannada District.Table 4 lists the spatial extent of various types of fragments(interior,perforated, edge,transitional and patch forests).Figure 7a depicts the domination of interior forests(64.42%)in the district.Edge forests(7.32%)are located along linear corridors,such as roads,rivers and boundary pixels of large forest patches.Patch forests are mainly located at the interfaces of forests,intermixed with agriculture and urban classes over small portions.The unscientific forest exploitation by the industrial sector peaked with the impetus of forest based industries during the period from 1960 to the 1980’s leading to selective felling of trees in the evergreen forests(Gadgil and Chandran 1989).This has created canopy gaps and the spread of invasive exotic species,adversely affecting faunal species.Mining activities in the district leave significant ecological,economic and social footprints much beyond the physical boundaries of mines by disrupting continuous forest patches(Ramachandra et al.2014b).

Figure 7b shows a decline in the area of interior forests from 64.42%(1979)to 53.33%(1989),with increase in edge forests(12%).The major activities during this period were industrialization,infrastructure development,intensified agriculture,manganese mining,a ferromanganese plant,a paper mill and plantations.The provision of forest resources to industries at highly subsidized rates and permission to polluting industries in the ecologically sensitive regions have contributed to the decline of forests and contamination of natural resources.Unplanned developmental activities such as a series of large scale power projects,manganese mining,a ferromanganese plant,a paper mill and several irrigation projects have led to the retreat of forests,with their degradation evident in the form of barren hilltops.The mismanagement of Kans(‘sacred forests’protected by local communities)and reserve forests also aggravated the situation towards the loss of interior,contiguous forests(Chandran 1989).Figure 7c shows the fragmentation status for the year 1999;the region lost a major portion of its interior forest and reached 40.74% from 53.33%(1989)with the increase in edge forests to 16.35%.Drivers of these changes are the implementation of a series of hydroelectric projects,the construction of national routes NH-17,NH-63,NH-204, the Konkan railway line and other infrastructure projects.Figure 7d illustrates the status of forests(in 2013)with 25.62%interior forests and 17.48%of edge forests,as well as the loss of connectivity between interior forest patches.These interior forests exist now only in the form of protected areas-sanctuaries,protected areas,sacred groves or Kans.The areaundernon-forestshasincreasedfrom36.07 (1999)to 47.3%(2013)with an increase of edge and perforated patches.

Table 4 Temporal changes in forest fragmentation at landscape level from 1979 to 2013

PCA aided in prioritizing bio-geophysical and socioeconomic variables and explained the variation in spatial forest cover.Components PC1 and PC2(for 1979 data) explain a large proportion of the variation(79.67%) (Tables 5 and 6).Axis 1 represents the greater loading of variables such as evergreen to semi evergreen forest, dry deciduous forest,moist deciduous forest,slope, rainfall,and scrub/grass land.Axis 2 represents elevation,population density and other physical factors indicating that Supa,Siddapur and Ankola taluks with greater slopes had higher forest cover.Elevation as a variable had a partial effect on scrub/grass, while moist deciduous forest types in the Sirsi and Yellapurataluksformedclusters.Coastaltaluks(Karwar,Kumta,Honnavar)suggest the effect of the variables slope and rainfall.Mundgod represents the effect of dry a deciduous and moist deciduous cover Fig.9.

Table 5 Component loadings&variance for the year 1979

Together,PC1,PC2 and PC3(based on 2013 data) explain a large proportion of the variation(88.04%)presented in Tables 7 and 8.Axis 1 indicates the higher loadings of the variables evergreen,rainfall,slope,scrub/ grass and dry deciduous forests.Axis 2 represents elevation,moist deciduous and population density.Supa and Ankola taluks have evergreen forests cover on slopes and edge effects are observed towards the periphery due to road network.Forest transitions in the Ankola, Kumta,Honnavar and Karwar(coastal)taluks have lost major forest cover towards the west side due to agriculture and aquaculture activities.The Bhatkal taluk has experienced loss of evergreen cover over the 1979 to 2013 period at its lower slopes due to intense activities of agriculture,plantations and other anthropogenic activities. Haliyal and Mundgod taluks have lost deciduous forest cover due to agricultural land expansions with the construction of reservoirs.Our analyses substantiate the role of geophysical variables in the changes of vegetation cover. However,further analyses are required in order to account for the nonlinear relationships among variables,since land transformations are often associated with nonlinear dimensions of human demand,government policies and socio-economic factors(Gong et al.2011)Fig.10.

Table 6 Component variances for the year 1979

Discussion

Forest clearing due to anthropogenic activities has been a major ecological problem(Laurance 1999;Etter et al. 2006),affecting biodiversity.Fragmentations of forests have resulted in habitat destruction and changes in the dispersal and migration processes(Armenteras et al. 2003;Etter et al.2006;Eldegard et al.2015).Shrinkages in animal habitat have led to inbreeding pressure,resulting in the extirpation of species,highlighting the intimate relationship between species and habitat.Edges contain communities different from interior forests due to an altered climate with higher light availability,loss of soil moisture,increased incursion of predators and competitors.As a consequence of fragmentation,changes in microclimatic near edges have favored the establishment of alien species such as Lantana camara,Chromolaena odorata and other species.This has caused a decline in the of native species,particularly in forests highly fragmented by transmission lines such as in the Haliyal and Mundgod taluks and at the lower slopes of Supa and Sirsi taluks.Forest patches shelter rare endemic species in human dominated landscapes.Studies have reported their role in pollination,maintenance of different life cycle phases of species,diversity and seed dispersal by harboring honey bees,small mammals,avifauna survival and many others(Bodin et al.2006;Page et al.2010). The Uttara Kannada District with relic forests(sacred forests/groves)and highly productive landscapes conserve local biodiversity and offer important ecological services,as well as improvement in the livelihood of local communities(Ray and Ramachandra 2010).These intact forests provide shelter for wild fauna and also benefit village communities with an array of forests goods and services such as hydrological functions,fuel wood and timber.

Forest fragmentation analyses provide vital insights to the potential effects of human disturbances(Hobbs and Yates 2003),through spatial descriptors of landscapes besides social and economic factors.The analysis of spatial patterns of forest changes should aid in formulating appropriate management strategies to conserve thesethreatened ecosystems.Analyses of forest changes are possible due to the availability of temporal remote sensing data.However,the spatial resolution of remote sensing data plays an important role in the analysis of spatial landscape patterns.Kernel sizes and resolution are intimately related in fragmentation analysis.Smaller kernel sizes have the effect of decreasing the average inter patch distance by considering smaller patches in the neighborhood,which also reduces the edge effect.Fragmentation analysis is also affected by up-scaling the resolution of remote sensing data as changes in spatial heterogeneity on a micro scale may not be detected using coarse spatial resolutions(Kitron et al.2006). Bharath et al.(2012b)has elucidated the dependency of the role of spatial resolution on measuring landscape structures in order to understand the different landscape patterns.The results reveal that landscape metrics based on area or patch cover are sensitive to spatial resolution. Fragmentation indices are also explicit in sensitivity across various spatial resolutions of remote sensing data. Gupta et al.(2000)highlights the problem of upscaling data by comparing correlations among reflectance and areas of land use categories.Their results indicate the loss of land use information of about 40%in proportion where upscaling was performed.

Table 7 Component loadings for the year 2013

Kernel sizes and resolution are intimately related in fragmentation analysis.The kernel or scale size can be consistent with respect to the pixel resolution and an increase of kernel size does not change the input data but increases the width of the non-interior classes,at the expense of interior forests,which maintain their overall proportion and the shapes of their features.Ostapowicz et al.(2008)explains the relationship between the resolution and kernel size by assessing and monitoring the structure of landscape patterns from multi-scale landcover maps.Riitters et al.(2004)demonstrated the appropriateness of 5×5 kernel at a nominal resolution of 30 m.The scale effects on different forest composition and configuration were appraised by sensitivity analysis with various combinations of pixel size and kernel size parameters and comparing frequencies of pattern classes in the entire forest area under study.

Reforestation and afforestation practices have been helpful in partially addressing the negative consequences of forest loss,through carbon sequestration,erosion control and non-consumptive use of forest products. However,the introduction of exotic species wouldimpact the native forest patches in the neighborhood (Shigesada and Kawasaki 1997;Ramachandra et al.2013). The introduction of Acacia auriculiformis,Tectona grandis and other exotic species in the forested regions and grasslands of this district has appalling effectsonbiodiversityduetohabitatdestruction, decline in water resources and unavailability of food especially for grazing mammals,which in turn has affected the prey stock of wild carnivores(Rao et al. 2012).Many commercial plantations have come up in the valleys by removing the natural vegetation,even in places with ecologically important ecosystems such as the Myristica swamps(Chandran and Mesta 2001). Inthiscontext,thecurrentconservationstrategy needs to focus on the local regeneration of natural forests and maintain continuity of forests,which helps in sustaining the livelihood of dependent populations. This requires motivation,conviction and commitment amongmajorstakeholders,i.e.,theforestfringe dwellers as well as forest officials.Restoration of forests with native species at watershed levels will help in mitigating the impact of forest fragmentations and improve hydrological services and biodiversity.The existing village grazing lands needs to be demarcated and managed by involving local stakeholders.This would help in mitigating grazing impacts in natural forests and also improves the prospects of forest regeneration.Prohibition of clear felling in the intact primeval forests would aid in preserving the structure of this ecosystem while enhancing its functional aspects.Joint management of forests by involving all stakeholders-local communities and others,would help in curtailing illegal logging,encroachments,wildlife protection and sustainable management of forests.

Table 8 Component variances for the year 2013

Conclusions

The Western Ghats,the repository of diverse biological organisms is one among 35 global hotspots of biodiversity.This region has been experiencing large-scale land cover changes with the fragmentation of primeval forests.Current research is attempting to quantify spatialtemporal patterns of land use dynamics and fragmentation of forests using temporal remote sensing data.The analyses would be useful in evolving appropriate forest management strategies to mitigate impacts of adverse land use dynamics.Temporal land use analyses show a decrease in the evergreen forest cover from 57.31% (1979)to 32.08%(2013).Forest fragmentation analysis based on remote sensing data of the 1979–2013 period helpedinassessingthespatialpatternsofforests changes under patch,transitional,edge,perforated and interior cover.The study region,as an ecologically fragile area,is now left with only 25.62%of interior forests and the spatial extent of non-forests is 47.3%(2013), which highlights the need to restore forests.The district has 18.5%of its area under monoculture plantations in the Haliyal and Mundgod taluks.Conservation planning of forest ecosystems needs to be holistic at watershed levels involving all stakeholders.Restoration of forests with native species would enhance hydrological services and biodiversity.By active participation in forest restoration initiatives and micro-level planning,stakeholders of the Western Ghats are likely to gain as promoters and guardians of biodiversity and hydrology.Rendering such service would help in mitigating global climatic change and serve the cause of forest ecosystems in global biodiversity hotspots.

Competing interests

The authors declare that they have no competing interests.

Authors’contribution

All authors read and approved the final manuscript.BS carried our remote sensing data analysis,MDS participated in the discussion and field investigations and TVR worked on the manuscript,data interpretation and language editing.

Acknowledgements

We are grateful for the sustained financial and infrastructure support to ecological research in the Western Ghats to the following institutions:(i)the Ministry of the Environment,Forests and Climate Change,Government of India,(ii)the NRDMS division of the Ministry of Science and Technology

(DST),Government of India,(iii)the Karnataka Biodiversity Board,Government of Karnataka,(iv)the Western Ghats Task Force,Government of Karnataka and(v) the Indian Institute of Science.

Author details

1Energy&Wetlands Research Group,Center for Ecological Sciences[CES], Indian Institute of Science,http://ces.iisc.ernet.in/energy.2Centre for Sustainable Technologies(astra),Bangalore,Karnataka 560 012,India.3Centre for infrastructure,Sustainable Transportation and Urban Planning[CiSTUP], Indian Institute of Science,Bangalore,Karnataka 560 012,India.

References

Abdi H,Williams LJ(2010)Normalizing data.In:Salkind NJ,Dougherty DM,Frey B (eds)Encyclopedia of research design.Sage,Thousand Oaks,pp 935-938

Alix-Garcia J,Janvry AD,Sadoulet E(2005)A tale of two communities:explaining deforestation in Mexico.World Dev 33(2):219-235

Ameztegui A,Brotons L,Coll L(2010)Land-use changes as major drivers of mountain pine(Pinus uncinata Ram.)expansion in the Pyrenees.Global Ecol Biogeogr 19(5):632-641

Armenteras D,Gast F,Villareal H(2003)Andean forest fragmentation and the representativeness of protected natural areas in the eastern Andes.Colombia Biol Conserv 113:245-256

Asner GP,Broadbent EN,Oliveira PJC,Keller M,Knapp DE,Silva JN(2006)Condition and fate of logged forests in the Brazilian Amazon.PNAS 103:12947-12950

Atkinson PM,Lewis P(2000)Geostatistical classification for remote sensing:an introduction.Comp Geosci 26(4):361-371

Barbier EB(2001)The economics of tropical deforestation and land use:an introduction to the special issue.Land Econ 77(2):155-171

Bell AR,Caviglia-Harris JL,Cak AD(2015)Characterizing land-use change over space and time:Applying principal components analysis in the Brazilian Legal Amazon.J Land Use Sci 10(1):19-37

Bennett A,Saunders D(2010)Habitat fragmentation and landscape change.In:Sodhi NS,Ehrlich PR(eds)Conservation biology for all.Oxford University Press,UK

Bharath S,Bharath HA,Rajan KS,Ramachandra TV(2012a)Cost effective mapping,monitoring and visualisation of spatial patterns of urbanisation using FOSS.FOSS 4G:10/2012

Bharath HA,Bharath S,Durgappa SD,Ramachandra TV(2012b)Empirical patterns of the influence of spatial resolution of remote sensing data on landscape metrics.Intl J Eng Res Appl(IJERA)2(3):767-775

Bodin R,Tengo M,Norman A,Lundberg J,Elmqvist T(2006)The value of small size:loss of forest patches and ecological thresholds in southern Madagascar. Ecol Appl 16(2):440-451

Boogaert J,Ceulemans R,Salvador-van ED(2004)Decision tree algorithm for detection of spatial processes in landscape transformation.Environ Manage 33:62-73

Broadbent E,Asner GP,Keller M,Knapp D,Oliveira P,Silva J(2008)Forest fragmentation and edge effects from deforestation and selective logging in the Brazilian Amazon.Biol Conserv 140:142-155

Buskirk S,William H,Federick R,Smith W,Richard Knight L(2000)An Overview of Forest Fragmentation in the Southern Rocky Mountains 3-14.In:Richard LK, Frederick SW,Steven BW,William RH,William BBL(ed)Forest Fragmentation in the Southern Rocky Mountains,University Press of Colorado,Colorado

Cabral AIR,Vasconcelos MJ,Oom D,Sardinha R(2010)Spatial dynamics and quantification of deforestation in the central-plateau woodlands of Angola (1990-2009).Appl Geogr 31:1185-1193

Cagnolo L,Cabido M,Valladares G(2006)Plant species richness in the Chaco Serrano woodland from central Argentina:ecological traits and habitat fragmentation effects.Biol Conserv 132:510-519

Carolina M(1995)Edge effects in fragmented forests:implications for conservation.Trends Ecol Evol 10(2):58-62

Chandran MDS(1989)Environmental impact of forest based industries on the evergreen forests of Uttara Kannada District:A case study.Department of Ecology and Environment,Government of Karnataka

Chandran MDS,Mesta DK(2001)In:Uma Shanker R,Ganeshaia KN,Bawa KS(eds) On the conservation of the myristica swamps of the Western Ghats.IBH Publishing Co,Oxford,pp 1-29

Chen J,Chen J,Liao A,Cao X,Chen L,Chen X,He C,Han SG,Peng M,Tong LX, Zhang W(2015)China's global land cover mapping at 30 m resolution: A POK-based operational approach.ISPRS J Photogramm Remote Sens 103:7-27

Cho SJ,McCarl BA,Wu X(2015)Climate change adaptation via US land use transitions: A spatial econometric analysis.2015 annual meeting,January 31-February 3,2015. Southern Agricultural Economics Association,Atlanta,No.196684

Colson F,Bogaert J,Ceulemans R(2011)Fragmentation in the legal Amazon, Brazil:can landscape metrics indicate agricultural policy differences?Ecol Indicat 11:1467-1471

DeFries R,Foley J,Asner G(2004)Land-use choices:balancing human needs and ecosystem function.Front Ecol Environ 2(5):249-257

Diaz S,Fargione J,Chapin F,Tilman D(2006)Biodiversity loss threatens human well-being.PLoS Biol 4(8):277,1300-1306

Echeverria C,Coomes DA,Hall M,Newton AC(2008)Spatially explicit models to analyze forest loss and fragmentation between 1976 and 2020 in southern Chile.Ecol Model 212(3):439-449

Eldegard K,Totland ?,Moe SR(2015)Edge effects on plant communities along power line clearings.J Appl Ecol.doi:10.1111/1365-2664.12460

Etter A,Mcalpine C,Wilson K,Phinn S,Possingham H(2006)Regional patterns of agricultural land use and deforestation in Colombia.Agric Ecosyst Environ 114:369-386

Fahrig L(2003)Effects of habitat fragmentation on biodiversity.Ann Rev Ecol Evol Syst 34:487-515

Gong C,Chen J,Yu S(2011)Spatiotemporal dynamics of urban forest conversion through model urbanization in Shenzhen,China.Intl J Remote Sens 32(24):9071-9092

Gong C,Yu S,Joesting H,Chen J(2013)Determining socioeconomic drivers of urban forest fragmentation with historical remote sensing images.Landsc Urban Plan 117:57-65,ISSN 0169-2046,http://dx.doi.org/10.1016/j. landurbplan.2013.04.009

Gupta RK,Prasad TS,Krishna Rao PV,Bala Manikavelu PM(2000)Problems in upscaling of high resolution remote sensing data to coarse spatial resolution over land surface.Adv Space Res 26(7):1111-1121,ISSN 0273-1177, http://dx.doi.org/10.1016/S0273-1177(99)01127-8

Harper KA,MacDonald SE,Burton PJ,Chen J,Brosofske KD,Saunders SC,Euskirchen ES,Roberts D,Jaiteh MS,Esseen PA(2005)Edge influence on forest structure and composition in fragmented landscapes.Conserv Biol 19:768-782

Hobbs RJ,Yates CJ(2003)Impacts of ecosystem fragmentation on plant populations:generalising the idiosyncratic.Aust J Bot 51:471-488

Holway D(2005)Edge effects of an invasive species across a natural ecological boundary.Biol Conserv 121:561-567

Jackson DA(1993)Stopping rules in principal components analysis:a comparison of heuristical and statistical approaches.Ecology 74:2204-2214

Kampichler C,van Turnhout CAM,Devictor V,van der Jeugd HP(2012)Largescale changes in community composition:Determining land use and climate change signals.PLoS One 7(4):e35272.doi:10.1371/journal.pone.0035272

Karnataka Forest Department(KFD)Internal Evaluation Report,2009-10,Evaluation of Plan and Non-Plan works implemented during 2004-05,2005-06,2006-07& 2008-09,“Restoration of Degraded Forests,Biodiversity Management, Community and Tribal development”.Aranya Bhavan,Bangalore-560003.

Kindstrand C,Norman J,Boman M,Mattsson L(2008)Attitude towards various forest functions:a comparison between private forest owners and forest officers.Scand J Forest Res 23:133-136

Kitron U,Clennon JA,Cecere MC,Gürtler RE,King CH,Vazquez-Prokopec G(2006) Upscale or downscale:applications of fine scale remotely sensed data to Chagas disease in Argentina and schistosomiasis in Kenya.Geospat Health 1(1):49-58

Kuèas A,Trakimas G,Balèiauskas L,Vaitkus G(2011)Multi-scale analysis of forest fragmentation in Lithuania.Baltic Forest 17(1):128-135

Laurance WF(1999)Reflections on the tropical deforestation crisis.Biol Conserv 91:109-117

Laurance WF,Lovejoy TE,Vasconcelos HL,Bruna EM,Didham RK,Stouffer PC, Gascon C,Bierregaard RO,Laurance SG,Sampaio E(2002)Ecosystem decay of Amazonian forest fragments:a 22-year investigation.Conserv Biol 16:605-618

Li M,Huang C,Zhu Z,Wen W,Xu D,Liu A(2009)Use of remote sensing coupled with a vegetation change tracker model to assess rates of forest change and fragmentation in Missippi,USA.Int J Remote Sens 30:6559-6574

Lillesand TM,Keifer RW(1987)Remote sensing and image interpretation.John Willey and Sons,New York

Lindenmayer D,Hobbs J,Montague-Drake R,Bennett A,Burgman M,Cale P et al (2008)A checklist for ecological management of landscapes for conservation. Ecol Lett 11:78-91

Liu C,Frazier P,Kumar L(2007)Comparative assessment of the measures of thematic classification accuracy.Remote Sens Environ 107:606-616

Millennium Ecosystem Assessment(MEA)(2005)Ecosystems and human wellbeing:Desertification synthesis.World Resources Institute,Washington DC

Morris AJ,Donovan JJ,Strager M(2009)Geospatial analysis of climatic and geomorphic interactions influencing stream discharge,Appalachian Mountains, USA.Environ Model Assess 14:73-84.doi:10.1007/s10666-008-9145-7

Munroe DK,Southworth J,Tucker CM(2004)Modeling spatially and temporally complex land-cover change:the case of western Honduras.Professional Geographer 56:544-559

Nandy S,Kushwaha SPS,Dadhwal VK(2011)Forest degradation assessment in the upper catchment of the river Tons using remote sensing and GIS.Ecol Ind 11:509-513

Nelson GC,Geoghegan J(2002)Deforestation and land use change:sparse data environments.Agric Econ 27(3):201-216

Nelson GC,Harris V,Stone SW(1999)Spatial econometric analysis and project evaluation:modeling land use change in the Darien.Inter-American Development Bank,Environment Division,Sustainable Development Department,Washington,DC,33

Ostapowicz K,Vogt P,Riitters KH,Kozak J,Estreguil C(2008)Landscape.Ecology 23(9):1107-1117,http://dx.doi.org/10.1007/s10980-008-9271-2

Page NV,Qureshi Q,Rawat GS,Kushalappa CG(2010)Plant diversity in sacred forest fragments of Western Ghats:a comparative study of four life forms. Plant Ecol 206(2):237-250

Pan Y,Birdsey RA,Fang J,Houghton R,Kauppi PE,Kurz WA,Phillips OL, Shvidenko A,Lewis SL(2011)A large and persistent carbon sink in the world's forests.Science 333(6045):988-993,1201609

Pascal J-P(1986)Explanatory booklet of vegetation maps 1,2 and 3.Institut Fran?ais de Pondichéry,India.88 pp.http://hal.archives-ouvertes.fr/hal-00444285

Prasad PRC,Sringeswara AN,Reddy CS,Nagabhatla N,Rajan KS,Giriraj A,Dutt CBS(2009)Assessment of forest fragmentation and species diversity in North Andaman Islands(India):a geospatial approach.Int J Ecol Dev 14(F09):33-46

Ramachandra TV,Kumar U(2011)Characterisation of landscape with forest fragmentation dynamics.J Geogr Inf Syst 3(03):242

Ramachandra TV,Krishnadas G,Setturu B,Kumar U(2012a)Regional bioenergy planning for sustainability in Himachal Pradesh,India.J Energy,Environ Carbon Credits 2(1):13-49

Ramachandra TV,Subash Chandran MD,Joshi NV,Pallav J,Uttam K,Aithal BH, Prakash M,Rao GR,Vishnu M(2012b)Landslide susceptible zone mapping in Uttara Kannada Central Western Ghats,ENVIS Technical Report 28.Energy& Wetlands Research Group,Centre for Ecological Sciences,Indian Institute of Science,Bangalore,560012

Ramachandra TV,Bharath H,Aithal,Durgappa,Sanna A(2012c)Insights to urban dynamics through landscape spatial pattern analysis.Int J Appl Earth Obs Geoinf 18:329-343

Ramachandra TV,Subash Chandran MD,Bhat SP,Aithal BH,Rao GR,Vishnu M (2013)Status of forest in Shimoga district,Karnataka,Sahyadri conservation series 23,ENVIS Technical Report:53.ENVIS-Environmental Information System,Centre for Ecological Sciences,Indian Institute of Science,Bangalore

Ramachandra TV,Bharath S,Bharath A(2014a)Spatio-temporal dynamics along the terrain gradient of diverse landscape.J Environ Eng Landsc Manag 22(1):50-63

Ramachandra TV,Subash Chandran MD,Bharath S,Rao GR,Vishnu DM(2014b) Reclamation of mine regions at Bisgod:Approaches and challenges,Sahyadri conservation series 43,ENVIS Technical Report 80.CES,Indian Institute of Science,Bangalore,560012

Rao GR,Subash Chandran MD,Ramachandra TV(2012)Grasslands of Anshi-Dandeli Tiger reserve,Sahyadri conservation series 12,ENVIS Technical Report:36

Ray R,Ramachandra TV(2010)Small sacred grove in local landscape:Are they really worthy for conservation?Curr Sci 98:1078-1080

Riitters KJ,O'Neill R,Wickham JD,Jones B,Smith E(2000)Global-scale patterns of forest fragmentation.Conserv Ecol 4(2):3

Riitters K,Wickham J,O'Neill R,Jones K,Smith E,Coulston J,Wade T,Smith J(2002) Fragmentation of the continental United States forests.Ecosystems 5:815-822

Riitters KH,Wickham JD,Coulston JW(2004)A preliminary assessment of Montréal process indicators of forest fragmentation for the United States. Environ Monit Assess 91:257-276

Roy J,Saugier B,Mooney HA(2001)Terrestrial global productivity.Academic, London

Salvati L,Zitti M,Ceccarelli T(2008)Integrating economic and environmental indicators in the assessment of desertification risk:a case study.Appl Ecol Environ Res 6:129-138

Shigesada N,Kawasaki K(1997)Biological invasion:theory and practice.Oxford University Press,Oxford

Swan ARH,Sandilands M(1995)Introduction to geological data analysis. Blackwell Science,Oxford,p 446

Timar L(2011)Rural Land Use and Land Tenure in New Zealand.Motu Working Paper,11,Moto Economic and Policy Research,Wellington,http://motuwww.motu.org.nz/wpapers/11_13.pdf Accessed 14Sept 2015

Todd M,Kerr S(2009)How does changing land cover and land use in New Zealand relate to land use capability and slope?Motu Working Paper 09-17, Moto Economic and Policy Research,Wellington

Wickham JD,Riitters KH,Wade TG,Coulston JW(2007)Temporal change in forest fragmentation at multiple scales.Landscape Ecol 22:481-489

Yuan X(2008)A novel method of systematic error compensation for a position and orientation system.Prog Nat Sci 18:953-963

*Correspondence:cestvr@ces.iisc.ernet.in

1Energy&Wetlands Research Group,Center for Ecological Sciences[CES], Indian Institute of Science,http://ces.iisc.ernet.in/energy

2Centre for Sustainable Technologies(astra),Bangalore,Karnataka 560 012, India

Full list of author information is available at the end of the article

?2016 T V et al.Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/),which permits unrestricted use,distribution,and reproduction in any

medium,provided you give appropriate credit to the original author(s)and the source,provide a link to the Creative Commons license,and indicate if changes were made.

7 January 2015 Accepted:14 April 2016