Evaluating the Contribution of Trees outside Forests and Small Open Areas to the Italian Landscape Diversification during the Last Decades

Land use by humans strongly alters the landscape mosaic, either by reducing or increasing its heterogeneity. One of the most recent and widespread land use changes in Europe has been the spontaneous reforestation of marginal agricultural lands. These primarily affected small landscape patches, such as trees outside forests (TOF) and small open areas (SOA), often represent the most diversifying features of landscape’ structures. Nevertheless, only small-scale studies can be found in the literature and thus it remains a relatively unexplored issue. Integrating inventory and cartographic approaches, this work assesses changes in abundance, coverage, and average size of small patches in Italy between 1990 and 2013. Main results showed an overall increase in number and coverage of small patches during the reference period. The average patch size remains unaltered for TOF but decreases significantly for SOA, due to trees encroachment and canopy cover increasing in forests. Our findings confirm the important changes in Mediterranean land mosaics and contribute to a better understanding of current conditions and recent trends regarding TOF and SOA. The integrated approach has proven to be helpful for the large-scale assessment of small patches dynamics, representing a viable monitoring tool to encourage the inclusion of small patches in landscape policy and planning.


Introduction
Natural processes and human activities might change landscape patterns over time.This will affect the configuration and distribution of the spatial elements of the land mosaic and their contribution to ecological processes [1].The spatial elements of a land mosaic contribute to its heterogeneity and fragmentation, and have positive effects on ecosystem functionality and resilience [2,3].Patches-relative homogeneous areas that differ from their surroundings-are among the spatial elements of the land mosaic whose geometry and spatial configuration are important for ecosystem functionality [4].Indeed, the proportion of patch types and their spatial configuration within the land mosaic have a strong impact on ecological processes and, in turn, on the functioning of socio-economic and ecological systems [5,6].
Recently, humans have driven landscape dynamics much more than natural evolution [7,8].Since the Industrial Revolution, intensive farming systems and urban land expansion at a global scale have increasingly altered the structures and functions of the existing landscapes at an unprecedented rate [9].The anthropogenic influence is even more marked in landscapes characterized by higher levels of heterogeneity and built by centuries of human activities and modifications, such as those in the Mediterranean area [10][11][12].During the last decades, the Mediterranean landscapes have indeed experienced a general decrease of agricultural lands, mainly due to their abandonment, followed by tree encroachment and subsequent woodland expansion [13,14], particularly in mountain and less profitable marginal lands [15][16][17].Oppositely, the agricultural intensification and the expansion of urban areas, which usually follow a high scattered and fragmented spatial pattern, urban sprawl, and urban sprinkling [18,19], are occurring in lowlands and plains [20][21][22].These phenomena are strictly interlinked and can alter the sustainability of such systems, especially in the medium-long term [22].
Tree encroachment can have both positive and negative implications in terms of ecosystems' functionality and their capacity to provide ecosystem services (e.g., [23][24][25][26][27][28][29]).The risk associated with landscape homogenization [30] is particularly relevant if combined with urban expansion, leading to the loss of fertile soils and urban green spaces [31], as well as agricultural intensification in more profitable lands [32,33].The negative impacts of such land cover changes are particularly exacerbated when occurring at the expense of small landscape patches, due to their limited extension [27,30,32,34].Examples of small patches are single trees, small woods, and hedgerows outside the forested matrix, collectively identified as trees outside forests (TOF) [35], and the small open areas (SOA) enclosed by forests.TOF and SOA are intended here as homogeneous areas covering between 0.05 and 0.5 ha, which is the threshold set by FAO to discriminate forestlands [36].
The relevance of TOF as landscapes' features, particularly the hedgerows, has emerged since the 70s [37].During the last decades, the systems for their inventorying have been increasingly explored and reviewed.TOF contributes to human wellbeing at both the global and local scale [38], through the conservation of biodiversity [35] and the mitigation of climate change [39].TOF also offers a wide range of provisioning [38], regulating [40,41], and cultural services (ecotourism, recreation, and education; [38]).Conversely, to our knowledge, poor attention has been given to the systematic inventory of SOA in forested landscapes, as demonstrated by the fact that a clear and unambiguous definition is still missing in the literature.The spatial configuration of SOA has a high variability and varies depending on the successional stage.For example, SOA in early-stage secondary forests (i.e., areas under active trees encroachment phenomena) are usually settled close to the forest edges and characterized by sparse trees in low-density stands, whereas SOA in mature stands are enclosed within forest stands and characterized by more evident and clearly defined borders [42].These differences result in the ecological functionality of these patches, and therefore in the ecosystem services delivered [42][43][44].
The large-scale assessment and mapping of landscape dynamics from a socio-ecological perspective is a prerequisite for sustainable planning [8].However, assessing the dynamics of small forest patches in diversified landscapes is limited by the currently adopted minimum thresholds for forest cover and area in many inventories [45].Therefore, assessing the spatial configuration of small patches in fragmented landscapes and the relative influence on biodiversity and ecosystem services provision remains an open issue [46,47].Recently, different methodological approaches and techniques for small patches monitoring were developed and implemented using both inventory and cartographic approaches [48][49][50][51][52][53].The inventory approaches are convenient in terms of costs and are usually adopted to estimate the total coverage of small patches and additional quantitative attributes, such as trees characteristics and biodiversity.Especially thanks to the availability of very high resolution images, as well as the progresses in proximal and remote sensing techniques [52], the cartographic approach actually represents a reliable option to assess small patches.Nevertheless, current data and studies are often limited to small-scale and synchronic studies (i.e., for Italy, [53,54]), mainly due to the lack of historical data over large scales and high computational costs [51].Accordingly, though some isolated attempts to integrate small patches within large scale monitoring systems were recently made (e.g., for TOF, [39]), further efforts should be oriented towards developing time and cost-effective approaches to analyze their consistency and characteristics (e.g., shape and average size) as well as changes over time.Furthermore, particular emphasis should be given to statistical accuracy and reliability of such information to strengthen their supporting role in decision making contexts [55].
According to these issues, the present study aims to assess coverage, abundance, and average size of small patches; i.e., TOF and SOA, in 1990 and 2013 as well as their dynamics in this time-span, using Italy as a case study.Accordingly, a design-based approach integrating inventory and cartographic data was implemented with the threefold objective to (i) reduce costs and time for monitoring activities, (ii) provide statistically-sound estimates, and (iii) increase the informative power of pure inventory approaches.In this way, the present study provides useful insights to support landscape policy and planning, and conservation strategies, and to facilitate the inclusion of small patches in existing large-scale monitoring systems.

Study Area
The study area overlaps the entire national surface of Italy (about 300,000 km 2 ) with high ecological and socio-economic variability.The mountain chains, namely Alps and Apennines, dominate the northern and the central part of the country, respectively, have been widely characterized by outmigration, land use abandonment, and forest expansion in the last decades [56].During the same time-span, the low hills, plains, and coastlines experienced agricultural intensification and the growth of urban areas.Currently, agricultural areas cover almost 56%, forests 32%, and settlements 7.3% of the national territory.Urban areas are mainly located along the coasts and in flat areas where the majority of population is settled [21].As briefly introduced, changes in the socio-economic asset of the country (i.e., outmigration towards urban areas, industrialization, etc.) had strongly influenced recent land use and land cover changed dynamics nationwide [57].Hence, on the one hand, depopulation of mountain areas resulted in the abandonment of traditional land uses, thus facilitating forest expansion [58].On the other hand, urban growth on croplands increasingly affects low hills and plains [31].

The Classification System and Inventory Framework
The Italian Land Use Inventory (IUTI) has been adopted to implement the following methodology in order to analyze small landscape patches and their dynamics during the 1990-2013 time-span.IUTI has been promoted and implemented by The Italian Ministry of Environment and Protection of Land and Sea as a key instrument for the National Registry for forest carbon sinks [59].According to its tessellated stratified sampling scheme, the whole Italian territory was covered by a network of 25-ha squares and one sample point was randomly located within each square [60].For this work, a sub-sampling strategy was applied in order to further reduce cost and times for photointerpretation [60].Accordingly, about 13,000 IUTI sampling points were classified based on high-resolution imagery in 1990 and 2013, and according to the classification system proposed by the Intergovernmental Panel on Climate Change (IPCC) [61].Moreover, a minimum dimensional threshold of 0.5 ha was used to discriminate the land use categories [61].The coarsest classification adopted six land use classes: forests, cropland, grassland, wetland, settlements, and other lands.
In addition to the basic land use classification, during the photointerpretation phase, sampling points falling within homogeneous patches less than 0.5 ha in size and with a land use different from the surrounding area were marked as "enclosed patches".Particularly, we distinguished between: (i) small patches with trees within land use classes other than forests, namely TOF, and (ii) small patches without trees within forest, namely SOA.In general terms, the latter can be principally associated with: (i) human activities (i.e., clearings; [44]); (ii) natural disturbances (i.e., treefall causing tree canopy gaps [42]); (iii) the tree crown architecture typical of certain forest types due to the so-called "canopy shyness" phenomenon [52]; (iv) the traditional extensive grazing activities (e.g., the so-called "maggenghi" meadows in the Alps; [62]); (v) the canopy cover gaps due to the ongoing encroaching process in new forest stands [25].Taking into account forest dynamics and forest ecosystem functions and related services (e.g., biodiversity conservation), in this study we divided SOA into two main sub-types, namely enclosed small open areas (encld-SOA) and encroaching small open areas (encrg-SOA).The former mainly corresponds to gaps associated with traditional grazing activities, while the latter to gaps in open recent secondary forest stands.See Figure 1 for some examples of encld-SOA and encrg-SOA.SOA due to forest harvesting or natural treefall are not considered in the present study because of their temporary nature (i.e., short-term natural evolution processes).
Forests 2018, 9, x FOR PEER REVIEW 4 of 13 "canopy shyness" phenomenon [52]; (iv) the traditional extensive grazing activities (e.g., the so-called "maggenghi" meadows in the Alps; [62]); (v) the canopy cover gaps due to the ongoing encroaching process in new forest stands [25].Taking into account forest dynamics and forest ecosystem functions and related services (e.g., biodiversity conservation), in this study we divided SOA into two main sub-types, namely enclosed small open areas (encld-SOA) and encroaching small open areas (encrg-SOA).The former mainly corresponds to gaps associated with traditional grazing activities, while the latter to gaps in open recent secondary forest stands.See Figure 1 for some examples of encld-SOA and encrg-SOA.SOA due to forest harvesting or natural treefall are not considered in the present study because of their temporary nature (i.e., short-term natural evolution processes).In order to apply the methodological approach based on the integration of inventory and cartographic data proposed by [63,64] for assessing green spaces in urban areas, each polygon surrounding a sampling point classified as TOF or SOA was drawn for both inventory occasions.This integrated approach allows to estimate not only the total coverage, but also the abundance and average size of small patches, which can offer a valuable contribution for better understanding and characterize the historical dynamics as well as their possible implications on landscape structure and functionality.

Trees outside Forests and Small Open Areas Estimation
Total coverage of small patches (A), abundance (N) and average size (a) at both inventory occasions were estimated, applying the design-based estimation approach developed by [48].Let Q be the extent of the area covered by the n square grid cells fully overlapping the national territory under the IUTI sampling scheme, the estimate of A is given by: where: where nu is the number of sample points classified as small patches (or their sub-types).The variance of A can be estimated as: Let S and aj be, respectively, the set of small patches (or their sub-types) selected by the n sampling points and the size of the j th patch.Whether aj values are negligible with respect to Q, the estimate of N is given by: In order to apply the methodological approach based on the integration of inventory and cartographic data proposed by [63,64] for assessing green spaces in urban areas, each polygon surrounding a sampling point classified as TOF or SOA was drawn for both inventory occasions.This integrated approach allows to estimate not only the total coverage, but also the abundance and average size of small patches, which can offer a valuable contribution for better understanding and characterize the historical dynamics as well as their possible implications on landscape structure and functionality.

Trees outside Forests and Small Open Areas Estimation
Total coverage of small patches (A), abundance (N) and average size (a) at both inventory occasions were estimated, applying the design-based estimation approach developed by [48].Let Q be the extent of the area covered by the n square grid cells fully overlapping the national territory under the IUTI sampling scheme, the estimate of A is given by: where: where n u is the number of sample points classified as small patches (or their sub-types).The variance of A can be estimated as: Let S and a j be, respectively, the set of small patches (or their sub-types) selected by the n sampling points and the size of the j th patch.Whether a j values are negligible with respect to Q, the estimate of N is given by: with estimated variance equal to: (5) Accordingly, the estimate of a is given by A/N, i.e., with estimated variance equal to: Considering the standard normal distribution of the variable, a statistical test was performed to assess whether a real variation occurs between estimates at 1990 and 2013.Accordingly, in the present study we considered variation as (i) extremely significant for p value < 0.01, (ii) very significant for p value < 0.05, (iii) significant for p value < 0.1, and (iv) not significant for p value > 0.1.

Results
The main findings of the study reveal that the overall number and coverage of small patches, considering both TOF and SOA, increased during the 1990-2013 period.On the contrary, their average size slightly decreased, even if results in this case show different trends between TOF and SOA.Indeed, the average size of TOF essentially remained stable, while that of SOA decreased during the observed time-span (Table 1), passing from 0.16 to 0.07 ha (−56% than the average size at 1990).Even though the increase in the number of patches is higher for SOA (+181.6%)than for TOF (+17.1%), the opposite trend is visible for their total coverage (+18.0%and +27.4%, respectively).The estimates of standard errors are satisfactory for all features at both reference times, ranging from 1.3% to 12.5%.
Table 1.Number (N), coverage (A, in ha) and average size (a, in ha) of small patches (SOA and TOF) at 1990 and 2013 (±their estimated relative standard errors, in percentage).Significance levels (not significant for p value > 0.1, significant for p value ≤ 0.1, very significant for p value ≤ 0.05 and extremely significant for p value ≤ 0.01) of changes between 1990 and 2013 (∆) are also reported in the summary tables with the relative symbol (n.s., *, **, ***, respectively).To better understand the dynamics of small patches, their changes were also compared with those of the two land matrices in which they can occur (i.e., forestlands for SOA, and non-forest lands except water and other lands without vegetation, for TOF), using data reported by [65] as a reference for land use estimates at 2013 in Italy.During the 1990-2013 period, the relative coverage of both TOF and SOA with respect to the surrounding land use matrix increased (+0.5% and +0.7%, respectively).In fact, the relative coverage of TOF passed from 1.6% to 2.1% while that of SOA passed from 5.5% to 6.2% during the observed period.However, it is worth to highlight that while forests increased by about 550,000 ha, the non-forest land uses decreased by about 580,000 ha during this time-span.
The overall increase of SOA, in terms of both abundance and coverage, is almost solely due to the dynamics of encld-SOA.In fact, as shown in Table 2, encld-SOA contribute to 93% of the overall increase in number of SOA, and even more, they are able to even compensate for the negative trend concerning encrg-SOA from 1990 to 2013.Table 2. Number (N), coverage (A, in ha) and average size (a, in ha) of SOA (encld-SOA and encrg-SOA) at 1990 and 2013 (±their estimated relative standard errors, in percentage).Significance levels (not significant for p value > 0.1, significant for p value ≤ 0.1, very significant for p value ≤ 0.05 and extremely significant for p value ≤ 0.01) of changes between 1990 and 2013 (∆) are also reported in the summary tables with the relative symbol (n.s., *, **, ***, respectively).Accordingly, the relative contribution of encld-SOA respect to SOA as a whole, passed from 64.8% to 74.3% in terms of coverage, and from 63.4% to 82.6% in terms of abundance, during the observed period.The average size decreased for both encld-SOA and encrg-SOA, but with major emphasis in the former, for which the 2013 value is 0.06 ha.Even in this case, estimates of standard errors are satisfactory, being in any case below 12%.
Relative transitions between the two SOA sub-types (i.e., encld-SOA and encrg-SOA) and with respect to forest and non-forest matrices (i.e., land use change occurred during the time-span towards a different land use class), were also analyzed (Table 3).It is worth to highlight that values in rows refer to flows from 1990 to 2013 (loss), while values in columns to flows from 2013 to 1990 (gain).Overall, the results confirm a general increase of SOA coverage, with encld-SOA losing 7.5% mostly due to their shift toward other land uses (i.e., forests), and encrg-SOA losing about 13.3% and mainly converted into encld-SOA (74.4%).Contrarily, encld-SOA show a higher overall increase than encrg-SOA (24.2% and 10.5%, respectively).This gain predominantly occurs at the expense of other land uses for encrg-SOA, while amounts to almost 59% for encld-SOA with the remaining 41% occurring at the expense of encrg-SOA.

Discussion
The results reveal that the number and relative coverage of small patches generally increased across the entire Italian territory from 1990 to 2013.In particular, our results show that TOF have increasingly contributed to the diversification of semi-natural (i.e., agricultural lands) and artificial land uses.On the other hand, the decrease in average size of SOA and their strong relation with successional stages, as demonstrated by the transition from encrg-SOA to encld-SOA and from the decrease of the latter within established stands, indicate that their contribution might slightly decrease in the medium-long term due to forest expansion (landscape homogenization).
These findings extend those of several studies carried out at regional and local scales.However, the comparison with such studies is not an easy task, because of the different thresholds usually fixed to distinguish small from large patches.In addition, forests are usually aggregated, and enclosed small open patches are included in the open forest class, as a mixture of agricultural and non-agricultural lands, and thus they are difficult to distinguish from aerial photographs [66].The same occurs at the global scale, where spontaneous reforestation led to an increase and a decrease of forest area and open patch sizes and numbers, respectively, but with some inconsistencies probably related to the lack of distinction between large and small patches [29].Usually, the reduction of the total size of open lands is interpreted as negative for ecological connectivity in Europe, because a large part of the plant and animal biodiversity is associated with open habitats of human origin.However, species with small home ranges and a short dispersal distance are more influenced by changes at the small patch level, while only species with larger home ranges and longer dispersal distances will be affected by the loss of larger patches.
The relative coverage of TOF with respect to the national surface increased from 1.1% at 1990 to 1.4% at 2013, with an average patch size almost stable around 0.1 ha.These results are in line with those obtained by [53] in Molise Region, where TOF cover 1.7% of the regional territory, having an average patch size of 0.14 ha at 2008.Our results are slightly lower than those obtained by [54] in Central Italy, where TOF cover almost 2.3% of the territory and with an average patch size of 0.21 ha, but with large variability among the four regions considered.However, as explicitly reported by [54], the analyses were performed on four administrative regions where TOF coverage is particularly high, thus justifying these differences.If our estimates are compared with those reported by the 2017 Great Britain's National Forest Inventory [67], which used the same classification system adopted by the present study, Italy has less TOF relative coverage than GB (1.4% and 3.2%, respectively), but a larger forest area [68].In spite of the recent debate around TOF, only few studies have been conducted so far and limited to regional scales [69].However, if compared with those by [51] at the regional scale in the Czech Republic during the 1953-2014 period, our results show similar dynamics about the overall increase of TOF coverage.Consistently with [51]' results, we found that TOF both disappeared because it was incorporated in recently expanded forest lands, and appeared in agricultural lands.Both processes have been driven by the abandonment of agricultural practices, which in turn originated in forest expansion and TOF formation and their enlargement in other land uses other than forest.Indeed, we observed that the increase of TOF in the 1990-2013 period inversely relates with the loss of agricultural lands (+27.4% and −8.5%, respectively; [65]).This trend can be interpreted as an important extension of the green infrastructure network, through which the ecological connectivity at lower scales is further strengthened (i.e., landscape to regional) [70].The mitigation of threats and pressures coming from the surrounding land uses and processes (e.g., urbanization and agricultural intensification) offered by TOF's network should be considered, and the TOF's functions evaluated according to different landscape patterns [54,71].In our case, the average patch size of TOF remained almost stable, despite a general increase in already existing TOF due to the small size of newly formed TOF.
The raising relative contribution in terms of coverage for both TOF and SOA has probably different characteristics whether taking into account the opposite trends in the land use matrices in which they occur.In fact, the relative contribution of TOF mainly increases because of the decrease of the land use matrix during the same period, while the increase in the relative contribution of SOA is exclusively due to the their own coverage increase.Moreover, results show that non-forested landscape are becoming more diverse thanks to the increase of small patches, while the decrease in average size of SOA and their coverage increase mostly related to encld-SOA, thus mostly coming from trees encroachment processes, could led to a progressive homogenization and simplification of forested landscapes in the medium-long term period.
Results show a reduction in average patch size of SOA, particularly regarding encld-SOA; (−62% with respect to the average patch size at 1990).This trend indicates the closure of these patches due to forest expansion and tree densification as a consequence of the progressive abandonment of extensive grazing and forest management practices [72].We observed a rate of reduction similar to what has been detected at local scales [30].Similarly, other authors have detected larger differences in the response of metrics over time at smaller spatial extents, and a stronger reduction in small patch size in completely abandoned forestlands, rather than in low-intensively managed watersheds [28].This is confirmed by the in-depth analysis of the IUTI sampling points classified as SOA in 1990 and 2013, showing that 45.5% of them experienced a change in the reference period (i.e., new SOA, disappearing or shifting from one sub-type to another).This is also consistent with the evident dynamics in Mediterranean forests, triggered by land use and climate change, especially regarding the modifications of species composition and spatial structure [73].Considering the key role played by forest structure on productivity [74], biodiversity [75], and the resistance/resilience [76] of forest ecosystems, these dynamics deserve particular attention from forest managers.At a broader scale (i.e., landscape), forest management and planning should valorize the role of small forest patches as contributing to ecological connectivity (in terms of e.g., species corridors, refugia, seed dispersers, and pollinators; [77,78]).
Looking at our results, the low increase in the number of patches of encrg-SOA compared to encld-SOA, and their relatively unchanged coverage overtime, might be related to the progressive slowdown of forest expansion processes in Italy during the most recent years [65].Considering the origin of encrg-SOA (usually associated with recent secondary forests; [79]), future active forest management should be considered as a valuable option towards the long-term maintenance of their functionality and stability [80].Results also suggest that recently formed encld-SOA are mainly represented by residual canopy gaps in the successional phases.Therefore, if a further decrease in forest harvesting and extensive grazing coupled with a progressive decrease in forest expansion in the future is hypothesized, a general negative trend in SOA in the medium to long term will be expected, as well as negative implications on their ecological roles.This is particularly true if considering that these small patches are prone to the rapid invasion of fast growing and light demanding species [80].Furthermore, it is worth noting that only 41.2% of SOA belongs to permanent encld-SOA, probably having the highest ecological value if it is considered that they effectively contribute to the heterogeneity of mature stands [42].Accordingly, from a management perspective, management activities oriented towards their maintenance should be set as a conservation priority.With this regard, a future scenario of increase of the harvesting rate could be also possible, consistently with the sustainable mobilization of wood [81].In this case, competent authorities should resist the temptation to develop a bimodal forest landscape pattern (frequent timber harvest, including clearcut, and no timber harvest) and to simplify the forest plan requirements, because this will result in the loss of small patches [82].

Conclusions
In the present work, we demonstrated that small patches do not represent a negligible component, at least in quantitative terms, of heterogeneous landscapes.TOF represent essential elements of the "green infrastructure" network (e.g., [83]), and SOA contribute to forest heterogeneity, biodiversity conservation, and ecosystem services' provision.The observed results underline that the management of secondary forests may reverse the impact of land use change (LUC), and give an opportunity to maintain ecosystem functionality and resilience in diversified landscapes.
Open patches play also a key role for the conservation of EU habitats [77].To bridge the gap between forest planning and habitat conservation, enclosed patches may be included in the ordinary silvicultural plans, and the (larger) open patches may be a focus of rural farming plans.A support is therefore needed for those restoration actions aiming at levelling out habitat gains and losses through the conservation of small patches.The current scenario of wood accumulation might conduct the forest landscape of Italy to a bifurcated model of intensive logging against total conservation, undermining the historical land mosaic patterns, namely of small patches, which are a key-habitat feature of small home-range species, as well as non-matching with stakeholders' expectations (e.g., in Switzerland; [66]).
From a methodological perspective, the integrated approach (i.e., coupling inventory and cartographic data) we propose can be considered particularly valuable for large-scale assessment of small landscape patches and related temporal and spatial dynamics.Moreover, the proposed approach seems to be cost-and time-effective in updating processes, as well as for obtaining accurate and reliable data.In addition, it enables a cross-comparison between small patches dynamics and surrounding landscape patterns.Considering the ecological and economic relevance of small forest patches, a deeper understanding of their spatial arrangement may increase their "edge-effect", and enhance the ecosystem services' flow, thus counterbalancing the disservices coming from the agricultural intensification [84].Further work should deal with the inventory of temporal and spatial changes in other patch metrics like shape, which might be equally important than size and density in assessing habitat quality and landscape ecosystem services.In conclusion, sustainable land management and conservation planning should (i) incorporate small forest patches and associated dynamics, into current intervention plans; (ii) balance biological conservation with other landscape services (e.g., small vs. large heaths; [85]); (iii) include effective landscape restoration actions [25]; and (iv) consider local knowledge, in combination with remote-sensing techniques, when assessing and mapping the historical landscape changes (e.g., participatory mapping; [45]).As at the basis of sustainability, the consideration of small elements in landscape planning indeed contributes to further understand the impacts of landscape complexity and the related slow processes of change for local wellbeing, otherwise neglected when looking at broader scales.

Table 3 .
Cross-tabulation matrix of the flows from one class to another during the 1990-2013 period, expressed in relative terms with respect to the total surface of SOA.Entries along the diagonal refer to persistence from 1990 to 2013.