The Land Cover Naturalness Index (LCN) describes the state and trend of land use/land cover (LULC) within the basin, according to the amount of human-induced transformation present. A basin in its undisturbed state, with intact forests and wetlands, generally maintains a sufficient quantity and quality of water to support indigenous flora and fauna. Naturalness exists on a gradient from completely natural to completely artificial or human dominated (Angermeier, 2000). Human conversion of lands and waterways are associated with increases in pollutant loads (non-point source from agriculture, point-source from urban and industrial), changes to infiltration and runoff regimes, and losses of regulating services (flood mitigation, erosion prevention, water purification). The Naturalness Index is, therefore, a proxy indicator for the degree to which these naturally-occurring functions are preserved within the basin. It is based on similar efforts to categorize and quantify this gradient over landscapes, such as Machado’s (2004) index of naturalness and the hemeroby hvb scale (Sukopp, 2004; Paracchini and Capitani, 2011; Walz and Stein 2014). More detailed investigations of specific LULC types, such as forest conversion to agriculture, may be warranted as a secondary step and can be calculated using the same data.

Attributes
Scale of calculation: Sub-basin; aggregate to single value per basin
Reference of indicator (if any): Based on index of naturalness methods described in Machado (2004)
Input required:
  1. Land use land cover data (raster or shapefile) for one year[1] , preferably, the most recent year for which land use land cover maps exist.

  2. Degree of Naturalness classification table

  3. Sub-basin shapefile

Suggested sources of ‘minimum’ data to enable calculation: ESA CCI land cover (2000, 2005, 2010, 300m resolution); Normalized Difference Vegetation Index (NDVI)

Data Preparation:

  1. Review and revise naturalness weights:

The Degree of Naturalness classification table (Deg_of_N.csv) below contains descriptions of LULC types as well as cultural practices (e.g., irrigation) that correspond to “naturalness” weighted on a 0-100 gradient. Sub-classifications are suggested based on three factors:

  1. Management of the water cycle: manually altering the flow and/or use of water to maintain a particular land-use type

  2. Pollution: chemical and physical pollutants entering the local water cycle due to human practices, such as fertilizer and pesticide use, and increased soil runoff from croplands as well as urban runoff and point-source wastewater loads from urban and industrial lands

  3. Vegetation characteristics: degree of native vegetation and permanence of vegetative cover

    The proposed weighting includes ranges of values to help highlight transitions from “natural” to “transformed” systems, i.e., from forests and wetlands to cultivated lands or from cultivated lands to urban areas (see Table below). It is strongly recommended that the default weights in the classification table be reviewed and, based on expert judgment, adjusted to be compatible with local conditions. For example, in some regions, flooded rice paddies may be considered to have a higher degree of naturalness than other irrigated crops, due to their ability to mimic some aspects of wetlands (which they may have replaced). In this case, a different classification and higher relative weight may be appropriate. Similarly, local or region-specific land use datasets may include highly detailed and differentiated classes of land use that will require expert judgment on their relative weight.

Table. Proposed “naturalness” characteristics and weights

Degree of naturalness Management of water cycle Pollution emissions Vegetation characteristics Examples Weight
Natural and semi-natural None None Native Forest (primary and secondary); lakes (natural) and wetlands; native grasslands; native shrublands 100
Cultural assisted system Low Low Mixed, high diversity Mosaic native vegetation (>50%, vegetation cover <50%) 70
Low Low Mixed, moderate diversity Mosaic cropland (>50%, natural vegetation <50%) 60
Transformed system Low Low Permanent cover with atypical species Permanent pasture land; agroforestry; tree crops 50
Low to Moderate Moderate to High Seasonal cover with atypical species Non-irrigated arable land 40
High Moderate to High Seasonal cover with atypical species Permanently irrigated arable land 30
Completely artificial High Moderate to High Sparse cover with grass Urban park space; low-density suburban areas; barren land 10
High High None Urban commercial areas; mining areas 0

Calculation in FHI Toolbox:

1. GIS processing

Intersect buffer polygon with land-cover data to generate table of land cover.

2. Classification (User-input)

For each land-cover type identified, assign weight based on decision-matrix above

3. Aggregation

Scores within a sub-basin can be calculated as:

$$ LCN = \left( \ \frac{\sum_{i = 1}^{m}{N_{i}\text{LCN}_{i}}}{\text{TN}} \right)*100 $$

where, TN is the total number of raster cells in buffer polygon of m type, Ni is the number of cells of land cover type ith, and LCNi is the appropriate score from the table directly above for that land-cover.