Temporal Indices of Landscape Change: A Proposal to Measure Variations in the Conservation Status of Vegetation at Fine Resolution

: The goal of this study was to propose scientiﬁc and objective indices capable of measuring the changes that occur in the conservation status of the vegetation of a particular area over a period of time. To this end, phytosociologically-based vegetation cartography at a detailed scale was used, carried out at two di ﬀ erent times, and the distance from the climax stage of the territory was calculated for each time. Three temporal indices of landscape change are proposed: Conservation Status Variation Index (ConSVI), Conservation Status Variation Velocity Index ( ConSVVe) and Change Ratio (ChanRat) . These enable the intensity, velocity, and percentage of change to be measured, and to determine whether this change is progressive or regressive—in other words, whether it is approaching or receding from the climax. To test the proposal, it was applied to a territory in Northwest Spain. The proposed indices are universally applicable to any territory and are the ﬁrst of their kind to operate at a detailed scale with a phytosociological basis. They also enable an objective measurement to be made of the landscape change that has occurred, meaning that they have immense practical utility in studies of managing and planning territorial resources.


Introduction
Vegetation changes as a result of ecological processes that act at different temporal and spatial scales [1]. It is necessary to understand vegetation changes in order to manage natural habitats. Toward this end, the classification of vegetation aims to summarize the variation of vegetation using a limited number of abstract entities [1]. Natural plant communities must therefore be differentiated from plantations, crops, and other non-natural formations. Such formations have benefited in recent years from the development of the human economy, with an attendant reduction in native forests. The natural resources management applied over the course of recent decades has not managed to stop the disappearance of natural vegetation. However, rural areas affected by depopulation and the reduction of livestock and agricultural impacts could improve the state of conservation of native forests.
How to measure the conservation status of vegetation is a question that has been approached in various scenarios [2], some of them based on the study of phytosociological communities at the landscape level [3][4][5] and others requiring a more complete and detailed knowledge of plant composition. In accordance with other research [6], ecological studies need to be contextualized within the frame of a vegetation typology, since the vegetation type being assessed may considerably affect the results of the studies.
The present study proposes to measure the variations in the conservation status of a territory's vegetation based on a phytosociological typology. Phytosociology is a branch of vegetation science that studies, among other concepts, current plant communities at the spatial size of vegetation patches [7]. Current phytosociology takes a statistical approach that aims to characterize vegetation types by the combined information from many different plots. Symphytosociology is the division of phytosociology that studies plant communities in their different states and how they replace each other according to their state of maturity in a specific homogeneous area [8,9]. The typological unit in symphytosociology is the sigmetum or vegetation series. Vegetation series includes the vegetation type representative of the series head and the serial communities that precede or replace it [10][11][12]. Thus, conservation status is viewed from a dynamic-serial perspective, as a measure of the distance to the hypothetical situation in which the entire territory is occupied by the series head or climax plant communities concerned.
The landscape level is needed to reflect vegetation health status in real time. The scale of the work and the typology of the markings that are used to divide up the landscape are conditioned by this level. From this perspective, the landscape can be regarded as a spatial configuration of markings of relevant dimensions for the phenomenon in question [13]. When the scale and perception of the analysis require the classification of the landscape mosaic based on vegetation markings, the use of excessively simplistic systems may lead to a loss of information, undermining the description of the place being studied. An example of this would be the decision to unify different scrubland communities within a single type when classifying the Orocantabrian landscape. Even communities that are close to each other in physiognomic and syntaxonomic terms may comprise totally distinct ecosystems reflecting different local climate conditions, playing a distinct serial role and harboring a different ecological complex [14]. The level of detail the researcher decides to use, depending on the goals set for the study, should account for the loss of information caused when choosing the types of landscape. When asking oneself what to differentiate and what to aggregate, one needs to be aware of the loss entailed by one's choice; there will be a need to base the decision on knowledge of the landscape's components-something that nearly always requires an appropriate floristic and vegetational foundation, without which the results may be erroneous.
From the 1970s to the present day, many efforts have been made to link satellite imagery and vegetation cover [15]. Specifically, near the study area, attempts have been made to study the recent dynamics of the landscape through aerial photographs [16]. Based on the previous attempts to detect changes in the vegetation cover, various indices have been proposed, most of them based on remote sensing data for land-cover change analysis, using measurements derived from infrared radiation and Landsat images [15,[17][18][19][20]. Some of these indices are widely used, such as the normalized difference vegetation index (NDVI) [21] and the enhanced vegetation index (EVI) [22]. These approaches are based on vegetation changes emerging from the study of satellite and orthophotographic images, which calculate an index or ratio derived from a series of numerical parameters. The potential of these vegetation indices has been validated by studies that analyze crop yields [23]. In general, such indices are applied to large territories and on not a fine scale. However, few studies propose indices that address changes in vegetation based on its state of maturity. For example, recently published research assesses change in vegetation in terms of the growth or shrinkage of plant cover, but does not take into account the naturalness of the change [24,25]. Structural changes in vegetation have also been addressed in recent works that analyzed this parameter on the basis of its homogeneity [25] or some very broad vegetation classes [26].
In order to be able to determine vegetation change, this study proposes indices that enable the conservation status of a territory to be calculated at a certain period of time and to compare it with an earlier period. Toward this end, the authors rely on symphytosociological concepts, according to which the optimal conservation status of a territory arises when it arrives at climax in all of its points [27]. The proposed indices enable researchers to ascertain and quantify how this conservation status of the territory varies in its entirety over time. It should thus be possible to determine whether the change in vegetation of a specific area is progressive, regressive, or neutral. The proposed approximation could be complementary to other indices, providing a type of information that may prove crucial in decision-making. This proposal could be very useful for determining the correct measures required in a territory and for interpreting the results obtained from the analysis of satellite images.

Materials and Methods
First, having selected the working area on which the various vegetation indices are going to be applied, the level of detail and required scale were established. Vegetation patches were defined as physiognomically homogeneous units with a shared tonality and texture that could be differentiated from the surroundings with the help of GIS (Geographic Information Systems). In order to be able to carry out a comparison at the landscape level, we used orthophotographs and georeferenced aerial photographs from different periods of time.
All operations involving the use of GIS were performed with ArcGIS ® software [28]. An important part of planning the work consisted of establishing the units or types of vegetation common to both moments, such that they enable the plant classification of the landscape and enable analysis of the changes that have taken place. These types of vegetation correspond to the grouping of communities that appear together on a detailed map, and were established as a new unit of work (Supplementary Materials, Annex 1).
Having established the types of vegetation, in order to determine the changes in the plant cover, the Potentiality Distance Index (PDI) put forward by Penas et al. [4] was used. This index helped us to interpret the state of conservation of the territory at the landscape scale, understanding it as a measure of the distance to the hypothetical situation in which the entire territory is occupied by the corresponding climax communities. According to the latter, the PDI for a given territory is calculated as the sum of the Potentiality Distance values for each type of vegetation (DI i ) multiplied by the relative area occupied by that type of vegetation (Ω i ) compared to the total surface area in question (Ω T ): where DI = 1 − [(3P − NI)/3n] and P: ordinal position that the vegetation type occupies in its respective vegetation series relative to the head of series; n: defined number of permissible stages resulting from the succession process; NI: Naturalness Index, proposed by Penas et al. [4] which assigns a weighting to the anthropic influence in the area where the vegetation type in question is located [29].
The PDI ranges of variation proposed by [4], as set out in Table 1, were reviewed and modified by the authors of the present study. A temporal dimension was incorporated into the PDI by studying the same territory at two different moments in time. Three new indices are proposed, which will jointly be referred to as where ConSVI is the Conservation Status Variation Index; PDI f is the Final Potentiality Distance Index; PDI i is the Initial Potentiality Distance Index. This index is the difference in the Potentiality Distance Index (PDI) value of the territory at the end of the passage of time in question (PDI f ) and the territory's PDI value at the outset of this passage of time (PDI i ).
In order to apply this index, a vegetation map is generated for each point in time of the comparison. Positive ConSVI scores indicate the recovery or dominance of a progressive direction in the plant succession-in other words, that there has been a narrowing of the distance to the respective heads of series in the dynamic processes. Therefore, an increase in the conservation status of the territory at the landscape scale is deemed to have occurred. Negative ConSVI scores on the other hand reveal a widening distance or regressive direction of the succession relative to the respective climax communities; therefore, a reduction in the conservation status of the territory is deemed to have taken place at the landscape scale. The absolute value of the ConSVI indicates the intensity of the progressive or regressive forces over the passage of time in question. This index can vary between −1 and +1, and seven categories were established as set out in Table 2. N is the number of decades (number of years of the time interval/10). This is deemed to be the variation in the conservation status of a territory over a unit of time. ConSVVe enables the conservation status variation (ConSVI) to be calculated in relative terms depending on the duration of the temporal series, and for the variations in the conservation status that have occurred in different zones to be compared relative to the potentiality distance at the landscape level Sustainability 2019, 11, 5887 5 of 16 over a certain period of time. The authors established a decade as the minimum temporal unit for dynamic-serial studies. It was considered that whereas regressive anthropic processes may take place over short periods of time (plantations, crops, clearance operations, opening mines, urbanization, etc.), dynamic-serial natural processes occur much more slowly, and their detectability was therefore estimated at ten years. ConSVVe retains the positive/negative sign of its respective ConSVI, entailing that a negative velocity indicates regression in the conservation status of a territory, and a positive velocity indicates an improvement.
This index ranges between −1 and +1. Higher positive values represent a faster recovery velocity towards the climactic stage. When the index has negative values, the higher its absolute values the greater the velocity of degradation will be-in other words the distance from the climax ( Table 3). Table 3. Relationship between the Conservation Status Variation Velocity Index (ConSVVe) scores and the variation in the territory's conservation status over a unit of time.

ConSVVe Scores Variation in Conservation Status Over a Unit of Time
Slow degradation 0 Neutral 0 to 0. 25 Slow recovery 0.25 to 0.5 Moderate recovery 0.5 to 1 Rapid recovery This index is derived by calculating the variation observed in the conservation status over the time period being studied relative to the initial conservation status. In other words, it derives from the values of ConSVVe and PDI i . It represents the ratio of the change observed per decade elapsed to the initial conservation status, and is expressed as a percentage. Thus, in two scenarios where the variation velocity is the same, the landscape with the lower initial conservation status (i.e., smaller PDI i value) will have the higher change ratio. The percentage of change in the landscape will be greater if the initial situation reflects more degraded successional stages. This ratio emphasizes the degree of recovery or alteration based on landscapes with a predominance of initial successional stages that are far removed from the hypothetical climactic situation. Negative ChanRat scores indicate regressive alterations, whereas positive scores indicate recovery or progressive changes towards the climax.
These three indices enable researchers to assess the progressive or regressive changes occurring in a territory in an objective manner and to quantify them in relation to the conservation status (increase, decrease, maintenance) based on knowledge of the types of vegetation present in the territory and their dynamic-serial significance in the landscape.

An Applied Example of Temporal Change Indices
In order to validate these indices, a territory located in the Cantabrian Mountains (León, Spain) was selected. It is made up of the sub-basin of the Dueñas river and has a surface area of 55.85 km 2 ( Figure 1). The district exhibits the typical topography of a mid-mountain landscape, comprising materials from the Paleozoic age, mainly highly fissured limestones and sandstones, pertaining to the Central Carboniferous Basin. The valley, which is approximately 11 km in length, is bounded by ridges and peaks with an altitude of around 2000 m. The average altitude is 1400 m, with a pronounced average slope of 21 • . From a biogeographical perspective it belongs to the Orocantabrian subprovince (European Atlantic Province, Atlantic-Central European Subregion, Eurosiberian Region). The predominant bioclimatic threshold is the supratemperate hyperhumid, where woods of beech and Sustainability 2019, 11, 5887 6 of 16 oak (Quercus petraea (Mattuschka) Liebl., Q. pyrenaica Willd.) flourish, followed by the orotemperate hyperhumid, which approximately occupies altitudes superior to 1600 m and where creeping junipers predominate. There are four small hamlets in the district, with a combined population of 118 inhabitants. The traditional manner of farming the mountain has molded a pastoral landscape that is now undergoing transformation owing to new uses. The uniform environmental and socioeconomic conditions enable the pertinent anthropic variables to be analyzed within the framework of the changes that have occurred in the territory [14].
Sustainability 2019, 11, x FOR PEER REVIEW 6 of 17 orotemperate hyperhumid, which approximately occupies altitudes superior to 1600 m and where creeping junipers predominate. There are four small hamlets in the district, with a combined population of 118 inhabitants. The traditional manner of farming the mountain has molded a pastoral landscape that is now undergoing transformation owing to new uses. The uniform environmental and socioeconomic conditions enable the pertinent anthropic variables to be analyzed within the framework of the changes that have occurred in the territory [14]. In order to ascertain the current vegetation, the authors established the UTM (Universal Transverse Mercator) grid of 1 km 2 as the working scale to delimit the patches, which led to the issuing of information layers at a fixed on-screen scale of 1:3000. Once the current vegetation was established and the types of vegetation ascertained, cartography was carried out at a scale of 1:65,000 and at two different times, 1957 and 2008. These years were selected because they corresponded to the oldest and most modern years of flight available, so this time period is a proof of concept. The DI for the vegetation types was calculated, as was the PDIi for the landscape in 1957 and the PDIf for the landscape in 2008. Using these figures, the authors applied the indices set out above to ascertain the changes in the territory's conservation status occurring over the course of the 51 years in question. In order to detect natural dynamic-serial processes (slower in time than those related to human action), the oldest series of aerial photographs available was selected for this study. In order to ascertain the current vegetation, the authors established the UTM (Universal Transverse Mercator) grid of 1 km 2 as the working scale to delimit the patches, which led to the issuing of information layers at a fixed on-screen scale of 1:3000. Once the current vegetation was established and the types of vegetation ascertained, cartography was carried out at a scale of 1:65,000 and at two different times, 1957 and 2008. These years were selected because they corresponded to the oldest and most modern years of flight available, so this time period is a proof of concept. The DI for the vegetation types was calculated, as was the PDI i for the landscape in 1957 and the PDI f for the landscape in 2008. Using these figures, the authors applied the indices set out above to ascertain the changes in the territory's conservation status occurring over the course of the 51 years in question. In order to detect natural dynamic-serial processes (slower in time than those related to human action), the oldest series of aerial photographs available was selected for this study. The fundamental tools used to carry out the cartography of the vegetation were orthophotography (for 2008) and georeferenced aerial photography (for 1957) (Figure 2). The orthophotographs derive from the summer photogrammetric flight carried out within the framework of the National Aerial Orthophotographic Plan (PNOA), in .ecw (enhanced compression wavelet) format. Each orthophotograph corresponds to a raster of 14,240 columns and 10,040 rows, with a cell size of 0.25 × 0.25 and a pixel depth of 8 Bits, which encompasses approximately 3.28 × 2.51 km 2 . The aerial photographs used to analyze the earlier landscape derive from digitalized stills from a flight in 1957, not ortho-corrected, saved in Geotiff format, with 6780 × 6071 pixels, 600 ppp, and 8 Bits of depth, held in the cartographic archives of the Duero Hydrographic Confederation [30]. The fundamental tools used to carry out the cartography of the vegetation were orthophotography (for 2008) and georeferenced aerial photography (for 1957) (Figure 2). The orthophotographs derive from the summer photogrammetric flight carried out within the framework of the National Aerial Orthophotographic Plan (PNOA), in .ecw (enhanced compression wavelet) format. Each orthophotograph corresponds to a raster of 14,240 columns and 10,040 rows, with a cell size of 0.25 × 0.25 and a pixel depth of 8 Bits, which encompasses approximately 3.28 × 2.51 km 2 . The aerial photographs used to analyze the earlier landscape derive from digitalized stills from a flight in 1957, not ortho-corrected, saved in Geotiff format, with 6780 × 6071 pixels, 600 ppp, and 8 Bits of depth, held in the cartographic archives of the Duero Hydrographic Confederation [30].
The 1957 stills in Geotiff format exhibited serious errors of calibration, some as great as ±300 m. To remedy this, the authors carried out a new prior referencing for each still that reduced the error to a maximum of ±100 m.   The 1957 stills in Geotiff format exhibited serious errors of calibration, some as great as ±300 m. To remedy this, the authors carried out a new prior referencing for each still that reduced the error to a maximum of ±100 m.
Each type of vegetation is defined by a physiognomic structure, a dynamic-serial context and a single naturalness index (NI), combined with height, bioclimatic, and geological data. This provided the foundation for applying the new indices using two maps of the area's vegetation, showing the landscape in 1957 and 2008 ( Figure 3). The fundamental tools used to carry out the cartography of the vegetation were orthophotography (for 2008) and georeferenced aerial photography (for 1957) (Figure 2). The orthophotographs derive from the summer photogrammetric flight carried out within the framework of the National Aerial Orthophotographic Plan (PNOA), in .ecw (enhanced compression wavelet) format. Each orthophotograph corresponds to a raster of 14,240 columns and 10,040 rows, with a cell size of 0.25 × 0.25 and a pixel depth of 8 Bits, which encompasses approximately 3.28 × 2.51 km 2 . The aerial photographs used to analyze the earlier landscape derive from digitalized stills from a flight in 1957, not ortho-corrected, saved in Geotiff format, with 6780 × 6071 pixels, 600 ppp, and 8 Bits of depth, held in the cartographic archives of the Duero Hydrographic Confederation [30].
The 1957 stills in Geotiff format exhibited serious errors of calibration, some as great as ±300 m. To remedy this, the authors carried out a new prior referencing for each still that reduced the error to a maximum of ±100 m.    A final map of the variation in the conservation status of the vegetation was made. To do this, a vectorial information layer (corresponding to a continuous surface formed by the union of different polygons defined by a type of vegetation in 1957 and 2008) was converted to raster format. The pixel size was 20 m 2 and contained the DI value of each type of vegetation. The difference between the final and initial states at pixel level was calculated. The subtraction of the final and initial raster layers allowed another layer to be obtained that reflects the evolution of the vegetation. The final map was reclassified into three classes: no change, progressive direction (advancement to climax), or regressive direction. Table 4 shows the values calculated for the DIs of each of the territory's vegetation types, the surface area covered by these in both 1957 and 2008, the PDI i value for 1957, and the PDI f for 2008. Twenty-five vegetation types were identified in the territory. The PDI i value was 0.564. This value suggests a landscape moderately removed from the climax communities and a good conservation status. The PDI f value was 0.621. This value also suggests a landscape moderately removed from the climax communities and therefore that the overall conservation status of the territory is good. Although both values fall within the variation range of the Sustainability 2019, 11, 5887 9 of 16 index, the higher value obtained for 2008 suggests a better state of plant conservation in the territory compared to 1957.

Calculation of the Temporal Indices of Landscape Change: ConSVI, ConSVVe, and ChanRat
The intensity and velocity of the changes at the landscape level over the period of time in question are shown in the following results (Table 5). Table 5. Scores for the indices calculated for the study area: initial indices (PDI i and PDI f ) and temporal change indices (ConSVI, ConSVVe, and ChanRat).

Indices
Abbreviation Score The conservation status variation index (ConSVI) obtained a score of 0.057, which falls within the "moderately positive" range ( Table 2). This indicates a recovery (albeit small) or a predominance of the progressive direction in the plant succession. Put another way, over the interval concerned, the distance to the respective heads of series in the dynamic processes narrowed. The authors therefore conclude that there has been a small enhancement in the degree of conservation of the territory at the landscape level between 1957 and 2008.

Initial indices
The score for the conservation status variation velocity (ConSVVe) was 0.011, which falls within the range of "slow recovery" over 10 years (Table 3).
The rate of change (ChanRat) obtained was close to 2%; this reflects the proportion represented by the observed change per decade relative to the initial conservation status. The rate of change was positive, albeit small, which indicates a change in a progressive direction from a scenario that was not very distant from the climax. The rate of change will enable future comparisons to be made with other territories.
The results of the three indices are clearly linked to each other. The positive values of the three indices indicate that the territory is recovering towards climax: ConSVI shows a moderate intensity of this progressive evolution, ConSVVe indicates that this recovery is slow, and ChanRat points to a small change as the territory was not much altered in its initial state. Although the first index can be applied alone, the others detail and complete the information on the changes that have occurred in the vegetation. All of them are very useful and complement each other.

Changes in the Structure of the Vegetation in the Sample Territory
The results of analyzing the cartography of the vegetation in the sample territory at the two points in time are shown in Figure 3. Annex 1 shows the correspondence between types of vegetation and phytosociological communities. In 1957, the landscape was dominated by the following types of vegetation: supratemperate perennial grasslands (39.30%), supratemperate chamaephytic scrubland (heaths and gorse) (10.74%), mature woodlands (7.84%), immature woodlands (6.75%), and a mosaic of supratemperate nanophanerophytic scrubland and perennial grasslands (broom and dandelion pastures) (6.13%). Grouping these together in large physiognomic units (grassland/scrubland/woodland), the overall pattern has the proportions 2.7/1.5/1. In 2008, the landscape was dominated by three types of vegetation: supratemperate perennial grassland (24.74%), mature woodland (21.25%), supratemperate chamaephytic scrubland (13.21%), supratemperate nanophanerophytic scrubland (broom and hawthorn) (10.87%), and orotemperate basophilous creeping juniper (7.71%). Grouping these together in large physiognomic units (grassland/scrubland/woodland), the overall pattern is characterized by the proportions 1/1.1/1. The percentages were calculated by measuring the surface areas occupied by each type of vegetation in the vectorial layer generated by the ArcGIS software.
A final map of variation in the conservation status of vegetation in the studied time period is shown in Figure 4. The land cover change data are shown in Table 6. In the entire territory studied, the percentage of vegetation cover in which no changes occurred was 39.62%; 35.85% of the territory advanced towards a climax state (progressive direction), while 24.53% moved away from the climax state (regressive direction). percentages were calculated by measuring the surface areas occupied by each type of vegetation in the vectorial layer generated by the ArcGIS software. A final map of variation in the conservation status of vegetation in the studied time period is shown in Figure 4. The land cover change data are shown in Table 6. In the entire territory studied, the percentage of vegetation cover in which no changes occurred was 39.62%; 35.85% of the territory advanced towards a climax state (progressive direction), while 24.53% moved away from the climax state (regressive direction).

Discussion
The indices proposed here, based on working units known as vegetation types, which entail enhanced knowledge of the vegetation, provide more precise and comprehensive information than that offered by other vegetation indices based on Landsat images. The latter have proved themselves to be useful in terms of recording the changes in plant coverage, but they do so in absolute terms of losses and gains, without taking into account the naturalness of such change [24]. One advantage of the proposed indices is that they can be applied to any territory. They may even be applied to arid regions, where the study of the vegetation by means of reflectance techniques may incur difficulties owing to the minimal wooded area occupied by the vegetation in these ecosystems, preventing substantial changes in the plant coverage from being detected [31]. The solutions to this problem currently involve the use of unmanned aerial vehicles (UAVs), as well as the establishment of vegetation types that require exhaustive knowledge of the plant communities [32].
The integration of distinct possible plant communities into a single dynamic-structural unit simplifies and facilitates the process of characterizing the landscape. In the present case study, it was especially useful with the stills from 1957, in which the task of identification at the community level would not be easy despite the availability of the current cartography. The poor image definition and the major changes in use that had taken place in a single location would render the characterization of the markings visible in earlier stills virtually impossible at the level of plant association. The vegetation types created should enable their identification and delimitation in the monochrome aerial photographs while also attaining the goals proposed: determining the structure and composition of the landscape, quantifying the findings recorded for each type of vegetation, and modelling any changes occurring over the time period in question. The resulting vegetation typology is therefore one of the most important results presented in this research. The result is not shown as a simple checklist listing vegetation types in a series of classes. This result enables the communities to be grouped in evaluable and comparable sets according to their naturalness, which helps the nature of the change to be understood based on very complete information. In addition, the results are reflected in a map that enables the most vulnerable areas, as well as those that have undergone negative changes, to be identified relatively quickly.
The results obtained from the temporal indices calculated for the case study suggest a change in the plant landscape of the territory in the direction of a higher maturity status, as is evident in Figure 4. The basis of this change is explained and accounted for by analyzing and interpreting the vegetation cartography. The most significant changes detected between 1957 and 2008, and which are likely to be relevant to the result obtained in the indices, were as follows: 1. A significant expansion of bushy formations at the expense of a reduction in the perennial grasslands, owing to successional processes (transition towards scrubland).
An expansion process involving bushy formations was observed, closer in the plant succession process to climax communities, in detriment to the area occupied by grasslands. Thus, the tracts of broom notably increased their area of occupation between 1957 and 2008. Heaths suffered a reduction, which is explained by the replacement of part of the land occupied by this type of vegetation in 1957 by forest plantations. The grasslands were reduced almost by half owing to an intense successional process, as well as replacement by cultivated forestry species. Sporadically, based on the evolution of the grassland, a type of vegetation not evident in 1957 emerged: the orophilous community of Orocantabrian oak Quercus orocantabrica Rivas Mart. et al. Despite being a species that hybridizes and forms part of woodland communities, the formations at these levels exhibit a dense bushy structure, similar to those of the scrublands of Genista obtusiramea Spach and Citysus oromediterraneus (G. López and C.E. Jarvis) Rivas Mart. et al., and their successional significance is similar.

Intense process of forestation in the supratemperate vegetation and the depths of the valley.
This change proceeded from earlier successional phases (perennial grasslands and scrublands) that transitioned into mature woodlands during the temporal period in question. Moreover, the socioeconomic changes causing the abandonment of crops and pasture meadows, which are no longer exploited by humans, determined the transition of these regions towards woodland communities. The extent of woodland formations increased almost twofold. In addition to increasing quantitatively, they have also increased qualitatively: whereas in 1957 they accounted for 53.72% of climatophilous phanerophytic formations, in 2008 they reached 88.58%.
These changes were fundamentally due to a demographic decline that involved a reduction in agricultural activity. The disappearance of dryland crops, shepherd-led livestock, certain management practices (e.g., irrigation of pasture, burning of gorse), and the replacement of small livestock (sheep and goats) with large animals (horses and cows) were the most important changes in the factors influencing the landscape. In general terms, the phenomenon of plants recolonizing a territory as a consequence of crop farming and other traditional practices being abandoned is one that has been observed in various parts of Europe [19,33,34].
The changes in the woodland vegetation are the ones that contributed most to the change detected by the proposed temporal indices. They involved a significant increase in the surface area occupied by climax communities in the territory (intense forestation) as well as a significant increase in the surface area covered by bushy stages close to climax, brooms being a case in point (transition to scrubland). In general terms these changes reflect a recovery in the vegetation, but also a homogenization of the landscape. This phenomenon has also been reported in areas of the Pyrenees, where the homogenization has been linked to the expansion of mature woodlands and with processes prior to heterogenization [19]. These authors link both phenomena to the elimination of preceding crops, which paves the way to the colonization of pioneering species and communities related to the natural succession and, subsequently, to the stability of the mature woodland. If the trends that have been detected persist, in 100 years the landscape will be less diverse, albeit more mature from the perspective of the serial dynamic. The homogenization of the vegetation cover related to the change in the uses of the territory is a phenomenon that has been shown in the Cantabrian mountain range at different scales and that has been related to the loss of ecosystem services [35]. This apparent contradiction between intervention and non-intervention in the natural world, management and conservation, diversity, and maturity, is a subject open to debate and of great importance to ecologists, botanists, politicians, and society at large [36,37]. Therefore, although a progression of vegetation towards climax stage is desirable, maintaining a degree of heterogeneity in the landscape mountain systems would be beneficial for the maintenance of certain ecosystem services [37]. The structure of the current landscape in the area in question, which is balanced and diverse, stems from a varied scenario dominated by pasturelands and tends, over the short term, towards a more homogeneous scenario dominated by woodlands.
An increase in the recorded surface area of vegetation in a territory is a phenomenon that has already been observed in other regions, such as China [20] and southwestern South America [38]. The fact that in some territories deforestation has stalled or even that regeneration of natural vegetation has taken place has, to a large extent, been driven by conservationist environmental policies, as has been postulated in China [20] and Argentina [38]. Although rapid economic growth and unchecked population growth lead to the degradation of vegetation and the environment-as has occurred in various parts of Mexico [18] and eastern Paraguay [38]-human activities are negatively correlated with elevation and with slope [39]. In addition, the degree of naturalness of the vegetation and the altitude have been positively related [2]. These phenomena account for the fact that such activities are concentrated in less-mountainous areas, triggering migration away from those rural areas with relatively harsh topography, as in the case of the area studied here.
The use of temporal indices enables trends of change to be detected objectively and globally for a certain territory, enabling subsequent detailed analysis to be carried out into which plant communities are shrinking or increasing, what losses of biodiversity these transformations entail, and how the use of the territory by humans is intrinsically linked to all these processes.

Conclusions
The temporal indices of landscape change being proposed enable the changes that befall the vegetation in a territory to be detected and quantified in an objective way.
These indices need to be applied on a cartographic basis established in accordance with the methodology of modern dynamic-serial phytosociology.
They constitute an efficacious and objective tool for detecting the trend of change in a landscape, the velocity at which it is taking place and the percentage of change that will occur over specified units of time.
The indices can be applied in all territories as long as the vegetation (its composition, structure and functioning) is known at a level of detail such as that provided by the phytosociological method.
Furthermore, the precision of the information provided by these indices means that, apart from its intrinsic value, it is of great utility for supplementing other indices that analyze changes in plant coverage, such as those based on the reflectance of the territory.
In terms of the data obtained when the temporal indices were applied to a sample territory, the ConsVI score (0.057) showed a positive variation, suggesting a predominance within the territory of the progressive direction of succession, which the authors interpret on the basis of their analysis of the vegetation as an increase in bushy plants and the nanophanerophytic scrublands close to climax. The velocity of change shown by the ConSVVe score (0.011) reflected slow but positive change, in other words slow recovery of the vegetation towards more mature and stable stages in the plant succession. Lastly, the rate of change indicated by ChanRat was 1.981 (almost 2%) over 10 years, which provides stronger evidence for this slow but continuous progress towards climax stages of vegetation at the landscape scale.