Understanding Driving Forces and Implications Associated with the Land Use and Land Cover Changes in Portugal

Understanding the processes of land use and land cover changes (LUCC) and the associated driving forces is important for achieving sustainable development. This paper presents the LUCC in Portugal at the regional level (NUTS II) from 1995 to 2010 and discusses the main driving forces and implications associated with these LUCC. The main objectives of this work are: (a) to quantify the land use and land cover (LUC) types (level I of LUC cartography) by NUT II in Portugal for the years 1995, 2007 and 2010; (b) to assess the spatio-temporal LUCC; and (c) to identify and discuss the main driving forces of LUCC and corresponding implications based on correlations and Principal Components Analysis. The results revealed large regional and temporal LUCC and further highlighted the different and sometimes opposite time trends between neighboring regions. By associating driving forces to LUCC, different influences at the regional level were observed, namely LUCC into agriculture land derived from the construction of dams (Alentejo region), or the conversion of coniferous forest into eucalypt forest (Centre region) associated with increased gross value added (GVA) and employment in industry and forestry. Temporal differentiation was also observed, particularly in the settlements that expanded between 1995 and 2007 due to the construction of large infrastructures (e.g., highways, industrial complexes, or buildings), which is reflected on employment in industry and construction and respective GVA. However, certain LUCC have implications, particularly in energy consumption, for which different behavior between regions can be highlighted in this analysis, but also on land-use sustainability.


Introduction
The land use and land cover (LUC) of territories is changing and these transitions have influences in the landscape [1][2][3][4].The evaluation of these changes, both spatially and temporally, is increasingly important in the context of the sustainable use of the territories [5].The evaluations of land use intensity and also the land use and land cover changes (LUCC) with environmental, economic and social impacts are particularly important in this context [6][7][8].
Methods for the preparation of land cover mapping have evolved, in particular, the methods of data retrieval by remote sensing (supervised or unsupervised analyses) that fostered the proliferation of studies about LUCC in different territories and scales [1,[9][10][11][12][13][14].
Deforestation is currently a major problem due to its large effect on the reduction of air and water quality and climate changes [23,24].These LUCC have indirect influence on the environment and contribute to the reduction of the quality of life, and further affect all living beings that depend on specific life conditions in the territories subject to these LUCC.In some countries, forest land has been turned into agriculture land [25][26][27][28].In order to maximize the agricultural production, some soils are being intensively exploited through the use of machinery and agrochemicals [29][30][31][32].However, the intense use of these soils generates negative impacts that contribute to the physiochemical and biological degradation of them [33], some of which are currently undergoing desertification process.
Another current issue is the high artificialization of soils [12], mainly explained by the urban expansion process [34][35][36][37][38][39] that include the construction of roads, telecommunications infrastructures or urban facilities [4].Soil sealing reduces the infiltration of rain water into the soil [40,41] and consequently the runoff increases, often causing floods, loss of human lives and material damage.
There are several research studies involving LUCC assessments that aim to identify the driving forces (e.g., [3,22,42]) and studies about the impacts on the environment, economy and society, and its influence on the development of the territories [11,[43][44][45].
The driving forces of land changes have been categorized into different types, such as economic, cultural, social, political, and, in some cases, may include two or more categories simultaneously [13,[46][47][48].However, the changing patterns can be the result of a complex interaction between different driving forces at different scales of action and LUC [49].The proven relationships between the LUC and economic variables have been analyzed and explained in detail [50,51].It is also important to consider the temporal component of the driving forces in this analysis [46] in order to understand if these have the same influence on LUCC over time, but also to figure out if these temporal variations are similar (among regions).
Landscape changes that occurred in the last decades in Portugal have been evaluated in some parts of the territory and in different contexts, such as the LUCC estimation in the LANDYN Project [52], in the reporting of emissions and carbon sequestration in the LUC sector and in governmental reports [4], and other publications revealing clear evidence of LUCC [52][53][54][55][56].These studies show that the great landscape transitions that occurred in the last decades in Portugal, derived largely from LUCC along the coastal areas with the increase of land artificialization and reduction of forest land due to forest fires and anthropic actions [28].
The evaluation of LUCC is crucial for landscape changes assessments and evaluation of territorial dynamics, as well as for the implementation of new investments, policy development (sustainability development) [37] and spatial planning actions among others.However, the assessment of landscape changes that quantifies the full spatial variations to allow regional differentiation, as well as the temporal variations to identify the main driving forces, has not yet been made in Portugal.
The identification of driving forces associated with LUCC has been done in different locations and in different contexts in Portugal, e.g., LANDYN Project, detecting the relationship between LUCC and the resulting impacts on natural resources, such as water [3,57].On the other hand, these studies do not address the identification of the relationships between LUCC and socio-economic disturbances in land-use planning.This topic has already been studied by other researchers albeit in other study areas [19,35,58,59].
In this context, the current research has the following objectives: (a) to quantify the LUC types (level I of LUC cartography) by NUT II in Portugal for the years 1995, 2007 and 2010; (b) to assess the spatio-temporal changes (LUCC); and (c) to identify and discuss the main driving forces of LUCC and corresponding implications (socio economic and environment factors).

Study Area
The study area is Portugal (only Continental Portugal, 88,962.50km 2 ), a European country divided into five NUTS II (Nomenclature of Territorial Units for Statistics, level II): North, Centre, Lisbon, Alentejo and Algarve.
The relief in Portugal is quite irregular, characterized by deep incised valleys surrounded by mountains in the North and by lower less rugged relief in the South.The climate is strongly influenced by the Atlantic Ocean and the Mediterranean Sea, performing the transition between the Mediterranean and the Atlantic climatic conditions.The rainfall regime is characterized by high spatial and inter-seasonal variability.The mean annual precipitation (MAP) ranges from less than 500 mm, in the south and northeast, to more than 2000 mm, in the northwest.The spatial variation of MAP reflects the influence of latitude, elevation and distance to the ocean.In fact, MAP tends to increase with increasing latitude, elevation and proximity to the Atlantic.Summer months (June, July and August) are particularly dry and the rainfall concentrates, mainly, in the period lasting from October to March.
The population is concentrated along the coast, especially in the metropolitan areas of Lisbon and Oporto.In recent years, some cities in the countryside have expanded towards the periphery (e.g., Viseu in Centre region), but rural areas still prevail.Some urban macrocephaly can be found in the study area which resulted from the concentration of activities and population that once migrated towards the large urban centers and settled there (rural exodus) [60].
In Portugal, there are large spatial contrasts in terms of LUC (Figure 1), within the NUTS II coinciding with large spatial variations of LUC for the years 1995, 2007 and 2010 (Table 1).The relief in Portugal is quite irregular, characterized by deep incised valleys surrounded by mountains in the North and by lower less rugged relief in the South.The climate is strongly influenced by the Atlantic Ocean and the Mediterranean Sea, performing the transition between the Mediterranean and the Atlantic climatic conditions.The rainfall regime is characterized by high spatial and inter-seasonal variability.The mean annual precipitation (MAP) ranges from less than 500 mm, in the south and northeast, to more than 2000 mm, in the northwest.The spatial variation of MAP reflects the influence of latitude, elevation and distance to the ocean.In fact, MAP tends to increase with increasing latitude, elevation and proximity to the Atlantic.Summer months (June, July and August) are particularly dry and the rainfall concentrates, mainly, in the period lasting from October to March.
The population is concentrated along the coast, especially in the metropolitan areas of Lisbon and Oporto.In recent years, some cities in the countryside have expanded towards the periphery (e.g., Viseu in Centre region), but rural areas still prevail.Some urban macrocephaly can be found in the study area which resulted from the concentration of activities and population that once migrated towards the large urban centers and settled there (rural exodus) [60].
In Portugal, there are large spatial contrasts in terms of LUC (Figure 1), within the NUTS II coinciding with large spatial variations of LUC for the years 1995, 2007 and 2010 (Table 1).Forest land and others LUC occupy more than 65% of the LUC of the country showing different spatial distributions.The Northern region is dominated by pinewoods (Pinus pinaster), rain-fed crops and shrub land, whereas, in the Centre region, these types of occupation are also dominant but there is a higher percentage of land occupation with Eucalyptus, which differentiates the forest species in comparison to the Northern region.The LUC of the Lisbon Region is dominated by rain-fed crops and settlements.In Alentejo, the forest is totally differentiated from the previous regions with predominance of evergreen oaks (Quercus suber and Quercus rotundifolia), rain-fed crops, and even with high percentage of grassland.Shrub land is predominant in Algarve and this region is also characterized by the dominance of Quercus suber in forest areas (Figure 2).Forest land and others LUC occupy more than 65% of the LUC of the country showing different spatial distributions.The Northern region is dominated by pinewoods (Pinus pinaster), rain-fed crops and shrub land, whereas, in the Centre region, these types of occupation are also dominant but there is a higher percentage of land occupation with Eucalyptus, which differentiates the forest species in comparison to the Northern region.The LUC of the Lisbon Region is dominated by rain-fed crops and settlements.In Alentejo, the forest is totally differentiated from the previous regions with predominance of evergreen oaks (Quercus suber and Quercus rotundifolia), rain-fed crops, and even with high percentage of grassland.Shrub land is predominant in Algarve and this region is also characterized by the dominance of Quercus suber in forest areas (Figure 2).

Data and Methods
Figure 3 presents the methodological scheme used in the LUC cartography, LUC type variations, spatio-temporal LUCC, driving forces and implications assessments.

LUC Cartography
The LUC data used in the evaluation of the LUCC of Portugal are derived from Land Cover Maps of Portugal for the years 1995, 2007 and 2010 produced by the General Directorate for Territorial Development (DGT), using a nomenclature (hierarchical) defined to support Kyoto reporting of emissions and carbon sequestration in Portugal [4].The LUC data have the following characteristics: scale 1/25,000; 1 ha of minimum mapping unit; vectorial data model (polygons); and minimum distance of 20 m between lines.These are the most recent LUC data with higher cartographic detail for the complete study area.The cartographic information of 2007 was obtained from photointerpretation of orthophotos (0.5 × 0.5 m), a process aided by IRS and AWIFS satellite imagery and geographic information of cadastral surveys (agriculture and forestry) done at DGT.All results were validated (geometrically and thematically) in order to obtain data with high accuracy and quality [61].The DGT used for the LLUC cartography of 2010 the same methodology used for the cartography of 2007.The LUC cartography of 1995, also produced by DGT, was obtained using the vectorial data of LUC boundaries of 2007 that were updated to the 1995 LUC, based on orthophotos and satellite images of this year.
In the present work, the LUC cartography for the year 1995 is only available with the nomenclature used in the Kyoto report.For this reason, and in order to allow LUC comparison between different periods, the legend of the LUC maps of 2007 and 2010 follows the nomenclature used in the Kyoto report.

Data and Methods
Figure 3 presents the methodological scheme used in the LUC cartography, LUC type variations, spatio-temporal LUCC, driving forces and implications assessments.

Data and Methods
Figure 3 presents the methodological scheme used in the LUC cartography, LUC type variations, spatio-temporal LUCC, driving forces and implications assessments.

LUC Cartography
The LUC data used in the evaluation of the LUCC of Portugal are derived from Land Cover Maps of Portugal for the years 1995, 2007 and 2010 produced by the General Directorate for Territorial Development (DGT), using a nomenclature (hierarchical) defined to support Kyoto reporting of emissions and carbon sequestration in Portugal [4].The LUC data have the following characteristics: scale 1/25,000; 1 ha of minimum mapping unit; vectorial data model (polygons); and minimum distance of 20 m between lines.These are the most recent LUC data with higher cartographic detail for the complete study area.The cartographic information of 2007 was obtained from photointerpretation of orthophotos (0.5 × 0.5 m), a process aided by IRS and AWIFS satellite imagery and geographic information of cadastral surveys (agriculture and forestry) done at DGT.All results were validated (geometrically and thematically) in order to obtain data with high accuracy and quality [61].The DGT used for the LLUC cartography of 2010 the same methodology used for the cartography of 2007.The LUC cartography of 1995, also produced by DGT, was obtained using the vectorial data of LUC boundaries of 2007 that were updated to the 1995 LUC, based on orthophotos and satellite images of this year.
In the present work, the LUC cartography for the year 1995 is only available with the nomenclature used in the Kyoto report.For this reason, and in order to allow LUC comparison between different periods, the legend of the LUC maps of 2007 and 2010 follows the nomenclature used in the Kyoto report.

LUC Cartography
The LUC data used in the evaluation of the LUCC of Portugal are derived from Land Cover Maps of Portugal for the years 1995, 2007 and 2010 produced by the General Directorate for Territorial Development (DGT), using a nomenclature (hierarchical) defined to support Kyoto reporting of emissions and carbon sequestration in Portugal [4].The LUC data have the following characteristics: scale 1/25,000; 1 ha of minimum mapping unit; vectorial data model (polygons); and minimum distance of 20 m between lines.These are the most recent LUC data with higher cartographic detail for the complete study area.The cartographic information of 2007 was obtained from photointerpretation of orthophotos (0.5 × 0.5 m), a process aided by IRS and AWIFS satellite imagery and geographic information of cadastral surveys (agriculture and forestry) done at DGT.All results were validated (geometrically and thematically) in order to obtain data with high accuracy and quality [61].The DGT used for the LLUC cartography of 2010 the same methodology used for the cartography of 2007.The LUC cartography of 1995, also produced by DGT, was obtained using the vectorial data of LUC boundaries of 2007 that were updated to the 1995 LUC, based on orthophotos and satellite images of this year.
In the present work, the LUC cartography for the year 1995 is only available with the nomenclature used in the Kyoto report.For this reason, and in order to allow LUC comparison between different periods, the legend of the LUC maps of 2007 and 2010 follows the nomenclature used in the Kyoto report.
LUC maps were created for the years 1995, 2007 and 2010, followed by the assessment of losses and gains of area by LUC type between the different years.LUC types include the level I, corresponding to the principal LUC (forest land and others, cropland, grassland, wetlands, settlements and other land), disaggregated into level II (Table 1).
LUC cartography of Portugal with similar characteristics was available only for those years that have conditioned the choice of LUCC analysis periods: the first period comprises 12 years, and the last three years.This can influence the presented results, especially when comparisons are made between the two periods, regarding LUCC, driving forces and implications.We have opted for the in-depth analysis of the relative variations calculated for each period.

LUC Type Variations
The variance between main groups of LUC types (corresponding at level I of LUC cartography) that differ from the average is obtained by the following coefficient of variation: where L cv is the coefficient of area variation of the LUC group i in a time series, σ is the standard deviation area of the classes (level II) that integrated the group of LUC (level I) and A the mean area of these classes (level II).
The measure of relative dispersion, L cv , is used to evaluate and compare different distributions of LUC areas, i.e., the variability of the areas relative to the mean.Small values of L cv are related to more homogeneous data sets (in this case, areas of the LUC types-level II-that integrated the principal LUC class corresponding to level I).Temporal changes in area by LUC type were obtained using the following equations [44]: where A k is the absolute variation area of the determined LUC type for the period k, A t1 is the area in the initial period, A t2 is the area in the final period and A Kr refers to the relative variation area of determined LUC type for the period k.
The quantification of LUC is an innovation to the DGT data because there is a lack of comparisons at the national level (level I) made with the LUC nomenclature using the Kyoto report [4].Furthermore, the coefficient of area variation (L cv ) is original as there are no comparisons between the LUC groups for the three periods under analysis.

Spatio-Temporal LUCC
To further understand the dynamics of transition between the various types of LUC, the thematic maps (Figure 2) (vector data) of the different years were crossed.The analysis of LUCC is supported by transfer matrices [62,63]: where A ij refers to the changed area from the i LUC type for period k to the j LUC type for the period k + 1, and n is the number of LUC types.
Using the transfer matrix the kappa coefficient is also calculated [64] to assess the spatial agreement of all the LUC products of the different years, using the following equation: where K is the Kappa coefficient, N is the total area, and P pi and P li are the areas in the row and rank, respectively, of a given LUC type.
The results of the variations by LUC type (areas) were correlated among the NUTS II.This procedure is intended to find relationships between the LUCC of different periods obtained from the data previously submitted and to determine if there are regional differences in the LUCC in the study area.

Driving Forces and Implications
Statistical data (socioeconomic and environmental) collected from the databases available on the web of the National Statistics Institute of Portugal (INE), National System of Hydrological Resources (SNIRH) and Pordata (Table S1 in Supplementary Materials) for the correspondent years of LUC data were used for the analysis of the driving forces and implications associated to LUCC.This analysis is primarily based on the identification of correlations between all collected data (see Supplementary Materials, Table S1) with the previously determined LUCC.Only the variables presented in Table 2 (driving forces and implications) were considered in the study because these are the ones that have the higher correlation with the LUCC (calculated previously) determined for different regions of Portugal.
The information about subsidies and agro-environmental measures was also considered when selecting the variables driving forces and implications.However, these data are not available for the entire period under analysis.A multivariate exploratory technique, the Principal Components Analysis (PCA), was performed using the software Statistica 7 (Stat Soft.Inc., Tulsa, USA) to develop the exploratory analysis of socio-economic and environmental standardized data with the LUCC.When performing this analysis for each type of LUC, only the variables that had a high correlation with LUCC were selected.Redundant variables were excluded from this analysis.
In recent decades, there were significant changes in areas occupied by water bodies [28], especially in the Alentejo, due to the construction of the Alqueva dam.Therefore, the volumes stored by watershed (VSW) were also considered, to find relationships with changes in area of irrigated crops or other LUC that need water.The VSW were analyzed separately with the LUCC, because the Lisbon region (NUTS) does not have any reservoirs or water storage.The watersheds were grouped according to their location in each NUTS and the average of all watersheds was considered for two years (year of LUC cartography and the previous one).

Land Use/Cover Types Variations
The LUC suffered wide changes and transitions between 1995 and 2010.According to Figure 4, the settlements distribution show that, on average, the deviations from the mean value of the areas of all the LUC classes reached approximately 14.5% of the settlements.The wetlands distribution reaches 10%, decreasing to 1% in the last distribution (forest land).It should be noted that intermediate changes in the classes that integrated each principal group of LUC (level I) can occur during this analysis, i.e., certain LUC classes (level II) can reduce their area and increase it later or vice versa, which requires a more detailed analysis of LUCC.These results are innovative and important to understand the LUCC in this territory.
A multivariate exploratory technique, the Principal Components Analysis (PCA), was performed using the software Statistica 7 (Stat Soft.Inc., Tulsa, USA) to develop the exploratory analysis of socio-economic and environmental standardized data with the LUCC.When performing this analysis for each type of LUC, only the variables that had a high correlation with LUCC were selected.Redundant variables were excluded from this analysis.
In recent decades, there were significant changes in areas occupied by water bodies [28], especially in the Alentejo, due to the construction of the Alqueva dam.Therefore, the volumes stored by watershed (VSW) were also considered, to find relationships with changes in area of irrigated crops or other LUC that need water.The VSW were analyzed separately with the LUCC, because the Lisbon region (NUTS) does not have any reservoirs or water storage.The watersheds were grouped according to their location in each NUTS and the average of all watersheds was considered for two years (year of LUC cartography and the previous one).

Land Use/Cover Types Variations
The LUC suffered wide changes and transitions between 1995 and 2010.According to Figure 4, the settlements distribution show that, on average, the deviations from the mean value of the areas of all the LUC classes reached approximately 14.5% of the settlements.The wetlands distribution reaches 10%, decreasing to 1% in the last distribution (forest land).It should be noted that intermediate changes in the classes that integrated each principal group of LUC (level I) can occur during this analysis, i.e., certain LUC classes (level II) can reduce their area and increase it later or vice versa, which requires a more detailed analysis of LUCC.These results are innovative and important to understand the LUCC in this territory.The class of settlements increased during the overall period under analysis owing, largely, to the urban sprawl of metropolitan areas (Lisbon and Oporto).However, the heaviest urbanization of soils occurred along the coastal areas [65], mainly in the North, Lisbon and Algarve.Major investments in road construction throughout the country were made during this period.This generalized artificialization of land is reflected in the results of the settlements coefficients of variation (Lcv) shown in Figure 4 (each group corresponding to level I of LUC, see bold legend in LUC maps).
The wetlands class also increased between 1995 and 2010, which is explained by a large increase in the construction of new dams or reservoirs for the use of surface water for energy production and for the supply of the populations and their activities (agricultural and industrial).The infrastructure with the largest contribution to the increase in the area occupied by wetlands was the Alqueva dam (Alentejo) that comprises a reservoir of about 250 km 2 .However, in some cases, new infrastructures that were built around these new water bodies resulted in landscape fragmentation.This process of The class of settlements increased during the overall period under analysis owing, largely, to the urban sprawl of metropolitan areas (Lisbon and Oporto).However, the heaviest urbanization of soils occurred along the coastal areas [65], mainly in the North, Lisbon and Algarve.Major investments in road construction throughout the country were made during this period.This generalized artificialization of land is reflected in the results of the settlements coefficients of variation (L cv ) shown in Figure 4 (each group corresponding to level I of LUC, see bold legend in LUC maps).
The wetlands class also increased between 1995 and 2010, which is explained by a large increase in the construction of new dams or reservoirs for the use of surface water for energy production and for the supply of the populations and their activities (agricultural and industrial).The infrastructure with the largest contribution to the increase in the area occupied by wetlands was the Alqueva dam (Alentejo) that comprises a reservoir of about 250 km 2 .However, in some cases, new infrastructures that were built around these new water bodies resulted in landscape fragmentation.This process of landscape fragmentation in wetlands has been investigated in different locations, because of its implications on the degradation of its ecological functions [64].
The other land class presented a L cv very close to the wetlands (Figure 4).This class includes shrub land and other types of LUC.However, the first type of occupation can evolve into forest or can be lost due to the occurrence of forest fires, which are very frequent in Portugal [24,66,67].This class ranges approximately 9% in the L cv .The smallest coefficients of variation of LUC were observed in cropland and forest areas.L cv results show that large variations of LUC occurred in the settlements group, in relation to the area occupied by these types of LUC in the year 1995.
The LUC transitions shown in Tables S2-S4 presented in the Supplementary Materials, registered a Kappa coefficient above 81% for the different periods under analysis (Table 3) and a high percentages of overall unchanged LUC areas.In the complete period between 1995 and 2010 overall changes reached 16.6%.

Spatio-Temporal Transitions
In Portugal, there are large spatial contrasts in terms of LUC (Figure 2).For the years 1995, 2007 and 2010 (Table 1), the largest spatial variations of LUC are visible at the NUTS II level of analysis.
LUCC analysis was differentiated in absolute (ha) and relative (%) terms in order to understand the landscape changes of each Portuguese region (Figure 5) between 1995 and 2010.Regarding the absolute variations (Figure 5A), the most important change is the loss of forest in the Centre of Portugal.The forest land area in relative terms does not present significant variations in the country, although this result does not reflect the internal transitions that occurred in this main LUC type.
Sustainability 2017, 9, 351 9 of 21 landscape fragmentation in wetlands has been investigated in different locations, because of its implications on the degradation of its ecological functions [64].
The other land class presented a Lcv very close to the wetlands (Figure 4).This class includes shrub land and other types of LUC.However, the first type of occupation can evolve into forest or can be lost due to the occurrence of forest fires, which are very frequent in Portugal [24,66,67].This class ranges approximately 9% in the Lcv.The smallest coefficients of variation of LUC were observed in cropland and forest areas.Lcv results show that large variations of LUC occurred in the settlements group, in relation to the area occupied by these types of LUC in the year 1995.
The LUC transitions shown in Tables S2-S4 presented in the Supplementary Materials, registered a Kappa coefficient above 81% for the different periods under analysis (Table 3) and a high percentages of overall unchanged LUC areas.In the complete period between 1995 and 2010 overall changes reached 16.6%.

Spatio-Temporal Transitions
In Portugal, there are large spatial contrasts in terms of LUC (Figure 2).For the years 1995, 2007 and 2010 (Table 1), the largest spatial variations of LUC are visible at the NUTS II level of analysis.
LUCC analysis was differentiated in absolute (ha) and relative (%) terms in order to understand the landscape changes of each Portuguese region (Figure 5) between 1995 and 2010.Regarding the absolute variations (Figure 5A), the most important change is the loss of forest in the Centre of Portugal.The forest land area in relative terms does not present significant variations in the country, although this result does not reflect the internal transitions that occurred in this main LUC type.A sharp reduction of Pinus pinaster and shrub land can also be observed in Central.The reduction of Pinus pinaster forests and Pinus pinea is due to the transition to Eucalyptus, shrub land, settlements, other broadleaves and rain-fed crops is the most important (see Table S4 in the Supplementary Materials).This LUC type has been affected by forest fires [24] fueled by the presence of a high number of resinous trees, forest maintenance and lack of shrub land cleaning.Forest fires cause major changes in the landscape, have destroyed a large part of the Portuguese forests [24,68,69] and the recovery of burned areas (pine forest) is very slow when the natural vegetation is regenerated without human intervention.The LUCC in Portugal are not unique within the larger picture of European LUCC, for example the reduction of pine forest in Catalonia (NE Spain) also caused by several forest fires [70].
Figure 5A shows that in the North and Centre regions there has been a large increase in the area occupied by Eucalyptus, settlements, and a reduction in the area of rain-fed crops.The significant A sharp reduction of Pinus pinaster and shrub land can also be observed in Central.The reduction of Pinus pinaster forests and Pinus pinea is due to the transition to Eucalyptus, shrub land, settlements, other broadleaves and rain-fed crops is the most important (see Table S4 in the Supplementary Materials).This LUC type has been affected by forest fires [24] fueled by the presence of a high number of resinous trees, forest maintenance and lack of shrub land cleaning.Forest fires cause major changes in the landscape, have destroyed a large part of the Portuguese forests [24,68,69] and the recovery of burned areas (pine forest) is very slow when the natural vegetation is regenerated without human intervention.The LUCC in Portugal are not unique within the larger picture of European LUCC, for example the reduction of pine forest in Catalonia (NE Spain) also caused by several forest fires [70].
Figure 5A shows that in the North and Centre regions there has been a large increase in the area occupied by Eucalyptus, settlements, and a reduction in the area of rain-fed crops.The significant transition of Pinus pinaster into Eucalyptus is explained by the investments in this type of forest to meet the demand of raw material for the cellulose industry [71].Eucalyptus is a fast-growing tree that allows several cuts in shorter time-periods, thus ensuring a higher yield for the producer.Furthermore, the high incidence of pinewood nematode vector in the region of Alentejo (coast) and in the Centre region was responsible for heavy losses of pine trees [72].
In relative terms, the variations of the different LUC types by NUTS II between 1995 and 2010 are very different, as shown in Figure 5B.In the North and in the Algarve, there was a high increase in the areas occupied by Pinus Pinea (175% and 85%, respectively), while in the Centre there was an increase in the areas taken by other coniferous (70%).
In the Lisbon region, there was mostly an increase in settlements when the absolute changes are observed (Figure 5A).Indeed the increase of settlements has taken place all across the country albeit more pronounced in the Alentejo and Algarve in relative terms (Figure 5B), owing largely to the construction of new resorts and other infrastructures to cater for the increase in tourism in recent decades [73].This land artificialization was the result of the marked transition that occurred on soils occupied by rain-fed crops, irrigated crops, Pinus pinaster and shrub land.However, between 2007 and 2010 the transition from irrigated crops was lower when compared with other types of LUC previously mentioned, indicating an increase in the use of agricultural land that had suffered from abandonment and disinvestment during the last two decades [74].The increase of the urban areas was also observed in other Mediterranean regions [5].
In relative terms the Lisbon region registered a rise of 170% in the area dedicated to rice crops between 1995 and 2010, and to a lesser degree in the soil occupied by other quercus.
In the Alentejo region the absolute changes in LUC show an increase of Quercus suber and other coniferous, and wetlands, followed by a high reduction of area occupied by rain-fed crops and grassland (Figure 5A); the dryland was converted into irrigated land due to the greater availability of water resulting from the construction of the Alqueva dam [28,52].In relative terms, the Alentejo shows an increase of 60% of wetlands areas due to the area occupied by the reservoir of the Alqueva dam, which spurred the intensification of irrigated crops, vineyards and olive trees.These LUCC are thus considered as water-dependents [13] highlighting the availability of water in quality and quantity for the LUC transitions.The expansion of water bodies provided an increase of water storage thus allowing the conversion of certain types of LUC to irrigated areas, and partially justifying the increase in the area of these types of LUC in the total period under review.Some examples of the expansion of agricultural areas in lowland Bolivia [26] and the increasing area dedicated to olive trees in the Mediterranean region [75] can be found in the literature.
Absolute changes of LUC in the Algarve region were characterized by an increase in the area occupied by Pinus pinea, settlements and by the reduction of soil occupied by shrub land.The same trend is observed in the relative changes, including the reduction of vineyards and grassland areas (Figure 5B).
When LUCC are analyzed per time variations (∆1995-2007 and ∆2007-2010) different results of correlation coefficients (r) and coefficients of determination (r 2 ) were obtained (Table 4).The Northern Region showed an increase in area for several types of LUC during the first period; however, a large proportion of these showed a decline in the second period, which accounts for the negative correlation shown in Table 4.The remaining NUTS II show a positive correlation of LUCC for the two periods, with the higher correlation in the relative changes observed in Lisbon.These results show that the LUCC tendency is very similar in the two periods in Lisbon.Although the majority of the relative changes in the remaining NUTS have positive correlations, they are not very high.It is worth noting that the analysis of the results should consider that the time periods do not have the same length which can condition the results of the LUCC dynamics.
Table 5 presents all possible correlations of the LUCC between each NUTS II.It is possible to observe higher correlations between LUCC of the NUTS II between the two dates (bold values in Table 5).However, we would like to stress the proximity effect of the NUTS II in LUCC with the closest NUTS showing positive correlation.This fact is explained by the distribution of different types of LUC throughout Portugal, where there is a high variation of area with LUC type between regions.The results show that the NUTS North and Algarve have the highest positive correlations in two periods for the relative changes, reflecting mainly the increase of Pinus pinea and settlements, and the reduced area occupied by grassland.These last results indicate similar trends in relation to variations of these LUC classes.

Driving Forces and Implications of the LUCC
The LUCC vary considerably over time in the NUTS II, but it is necessary to understand the driving forces in the LUCC in each region and possible resulting implications.The obtained correlations between socio-economic and environmental data available for the three years under review (see the variables in Table 2) and the LUCC (Table 6) show that variations of settlements correlate strongly with the employment (according to INE classification includes the following activities: agriculture, forestry, fishing, industry, construction and water), volume of exports, gross value added (industry, construction and water activities); resident population and roads.The previous factors are the driving forces behind the increase on artificial surfaces.
The changes in resident population per NUTS are positively correlated with the increase in settlements, owing largely to the construction of infrastructures (housing, commercial and industrial) during the last decade of the 20th century [28].The LUCC for settlements has some implications, in particular, the increase in energy consumption by domestic and industrial consumers that showed high positive correlations (bold values in Table 6).The forest land changes have different correlations with the variables considered in different subtypes of LUC covered by this class (Table 6).Employment in agriculture, forestry and fishing is the most important driving force in forest land changes, especially changes in areas occupied by Pinus pinaster, Eucalyptus, others quercus, broadleaves, coniferous and irrigated crops.In the case of Eucalyptus, the area variation was high, especially in the Centre (location of major pulp mills); this increase was reflected on the conversion of forested areas into Eucalyptus forest which brought about an increase in electricity consumption associated with the industrial sector, which is an implication derived from this type of LUCC.
LUCC of rain-fed and rice crops presented a high positive correlation with cereals production.The first type of LUC shows high correlation with the gross value added in the agriculture sector.In the case of variations of irrigated crops, there was a high positive correlation with the gross value added in agriculture sector, volume of exports, employment in agriculture and olive trees production.Although the LUC cartography of the irrigated crops is different from that of the olive groves, the statistical data include all the irrigated productions including the irrigated olive production from Alentejo.In recent years, the increase in number of irrigated plantations and new olive groves was due to the construction of a new dam and other infrastructures for water storage.These findings are attested in the case of irrigated olive groves in an extensive cultivation, especially in the Alentejo region, contributing to the increase of the gross value added in the agriculture, forestry and fishing.
It is noteworthy that, regarding the cropland class, the changes in the areas of vineyards correlated positively with wine production, volume of exports and employment in agriculture and industry.However, these LUCC also indicate the increase in electric energy consumption and in the number of electricity consumers in agriculture and industry.
The variation of the volumes stored by watershed is also clear in LUCC, especially in cultures that require water, as well as in the areas of expansion of the main water bodies in the total period under review.As for the variations of environmental expenditure in protection of biodiversity and landscape, this exhibit negative correlation with the LUCC which are potentially explained by fragmentation of LUC, indicating higher LUCC area mean lesser environmental expenditure or vice versa.
The per capita gross domestic product and the per capita purchasing power variables have very weak correlations with the LUCC.The per capita purchasing power presents a high negative correlation with the rain-fed and olive crops.These results can be explained by the increase and efficiency of the cereals and olive productions (increase of production volume in smaller areas) but also by the supply and demand of the food products markets and their derivatives from the types of LUC.
As previously explained, the LUCC in Portugal are associated with different driving forces, resulting in different implications at the regional level.These results were obtained using the Principal Component Analysis (PCA) performed with the LUCC in five NUTS II and the respective variables of driving forces and implications of these LUCC (Figure 6).These results showed distinct groupings among regions and also some asymmetry between driving forces and implications (Figure 6A,B).For example, the driving forces with LUCC of the Alentejo and Algarve regions in factor 1 (Figure 6A) present the higher positive values, while in the case of the implications (Figure 6B) the results are the opposite.According to the obtained results, factor 1 can be related with the spatial distribution of the regions and factor 2 identifies, in general terms, the LUCC associated with the driving forces as well as the implications.
Loadings of factor 1 of the driving forces analysis (Figure 6A) highlight the PCGDP and PCPP variables (0.47 and 0.59, respectively).The temporal analysis of these variables showed an increasing trend over time and also an effect of spatial continuity between the NUTS (increasing values from the North region to Lisbon and decreasing to the Algarve).These observations, in particular the effect of spatial continuity between NUTS, are also evident for the loadings of factor 1 of the implications analysis (Figure 6B), highlighting the high positive loadings for the CEEI and ECI variables (0.9 and 0.91, respectively).
The analysis of the loadings for the factor 2 of the driving forces analysis (Figure 6A) reflected the positive influence of the GVAA and OTP variables (0.60 and 0.68, respectively) representing the main driving forces for the analyzed LUCC.In the case of factor 2, the implications analysis (Figure 6B) highlights the CEEA positively (0.73).This variable constitutes one of the main implications of the observed LUCC.Among the considered LUCC and driving forces (Figure 6A), there is a clear distinction between the northern (North, Centre and Lisbon) and southern (Alentejo and Algarve) regions of Portugal.In fact, the landscape changes can be the result of the influence of regionally distinct factors (job offer, agricultural productivity, and water availability), culminating in the fragmentation of the territory.This fragmentation is also referred to in EEA reports [76] and Ribeiro et al. [77].fragmentation of the territory.This fragmentation is also referred to in EEA reports [76] and Ribeiro et al. [77].When performing the analysis of LUCC and its implications (Figure 6B), the clusters formed by regions are not so well defined (greater proximity or linkage distance) as in the case of driving forces, but there is still an approximation between them, providing a continuity from north to south or vice versa (Figure 6A,B).This means that implications associated with LUCC are observed in a spatially threaded process, with less demarcated transitions between regions when compared with the driving forces.These implications associated with environmental repercussions of LUCC have been mentioned by Chen et al. [78], and can include interference in energy sustainability in order to produce more energy (need for more dams, wind and photovoltaic farms and supply networks) to assure the energy demands of consumers (industry, agriculture and population).However, more energy supply also creates new possibilities to generate new LUCC.

Discussion and Conclusions
In Portugal, the LUCC are not static in time and space, because LUC that characterizes these landscapes is very dynamic along the total period under analysis.This LUCC dynamics are also referred to in other studies carried out in other territories (e.g., [17,20,63,79]).
In this research, LUCC exhibits different regional trends that are conditioned by the latitude factor, which was proved by the higher correlation among regions with higher spatial neighborhood.Different regional trends are also explained by different types of dominant LUC that characterize these regions at a certain time.The results show that the LUCC are differentiated for each NUTS II, considering the relative variation of area for each LUC type.
Crosschecking the area of LUCC with socio-economic and environmental indicators allowed the identification of the factors that have the greatest influence on LUCC, and consequently on landscape changes and ecological pressures (forest land degradation).However, the LUCC have different variations in the considered periods, in particular, the measures, policies and investments When performing the analysis of LUCC and its implications (Figure 6B), the clusters formed by regions are not so well defined (greater proximity or linkage distance) as in the case of driving forces, but there is still an approximation between them, providing a continuity from north to south or vice versa (Figure 6A,B).This means that implications associated with LUCC are observed in a spatially threaded process, with less demarcated transitions between regions when compared with the driving forces.These implications associated with environmental repercussions of LUCC have been mentioned by Chen et al. [78], and can include interference in energy sustainability in order to produce more energy (need for more dams, wind and photovoltaic farms and supply networks) to assure the energy demands of consumers (industry, agriculture and population).However, more energy supply also creates new possibilities to generate new LUCC.

Discussion and Conclusions
In Portugal, the LUCC are not static in time and space, because LUC that characterizes these landscapes is very dynamic along the total period under analysis.This LUCC dynamics are also referred to in other studies carried out in other territories (e.g., [17,20,63,79]).
In this research, LUCC exhibits different regional trends that are conditioned by the latitude factor, which was proved by the higher correlation among regions with higher spatial neighborhood.Different regional trends are also explained by different types of dominant LUC that characterize these regions at a certain time.The results show that the LUCC are differentiated for each NUTS II, considering the relative variation of area for each LUC type.
Crosschecking the area of LUCC with socio-economic and environmental indicators allowed the identification of the factors that have the greatest influence on LUCC, and consequently on landscape changes and ecological pressures (forest land degradation).However, the LUCC have different variations in the considered periods, in particular, the measures, policies and investments applied in each region.Thus, it was found that the landscape of Portugal suffered some changes between the analyzed periods, but this spatial differences are derived from certain LUCC in different regions, especially the resinous forest transition to other types of land cover (e.g., North and Centre) as well as the artificialization of soils (important in the Algarve region).The latter is marked essentially by the urban growth, that characterizes the complex urban system of the Iberian Peninsula [86,87].These new land cover types exert great pressures on the territory's ecosystems and they can cause problems in their sustainability, in particular in the consequent degradation of more fertile soils and other natural resources (e.g., drinking water) [36,57].
The temporal analysis of the driving forces results (Figure 6A) shows that the regions maintain proximity in most types of the LUC analyzed.However, there are some exceptions such as the settlements with greater linkage distance between 1995 observations and those of 2007 and 2010.Portugal went into economic recession in the last decade and this factor may have contributed to a slowdown in construction, which was particularly marked in the Centre and North regions where a higher percentage of settlements is located.According to Auch et al. [46] during times of recession new development typically declines and thus the rate of LUCC decreases.
By associating driving forces to LUCC different influences at the regional level were observed, namely LUCC into agriculture land derived from the construction of dams (Alentejo region), or the conversion of coniferous forest into eucalypt forest (Centre region) associated with increased gross value added (GVA) and employment in industry and forestry.Temporal differentiation of LUCC driving forces was also observed, particularly in the settlement expansion between 1995 and 2007 due to the construction of large infrastructures (e.g., highways, industrial complexes, or buildings), which is reflected on employment in industry and construction and respective GVA.However, certain LUCC have implications, particularly in energy consumption, for which different behavior between regions can be highlighted in this analysis.
It was also observed that certain LUCC have multiple driving forces of economic and social nature, like for instance the LUCC for settlements associated with employment (agriculture, forestry and industry), agricultural production (e.g., wine) and gross value added.The increase in irrigated crops was driven by the construction of dams, but also presents a strong correlation with the employment in the agricultural sector and the volume of exports.The implications associated with the analyzed LUCC highlights the consumption of electric energy.The increase of vineyards has been reflected in larger volume of exports to meet the increased demand for Portuguese wine products in recent years, according to the Portuguese Wine Institute (IVV).Because of these reasons, the soil conversion of these plantations was promoted and its expansion took place in different regions.
The driving forces of LUCC are very diverse depending on the territory.For example, in Cimandiri and Cibuni Watersheds (Java Island), Kelarestaghi and Jeloudar [79] highlighted the rainfall, soil type, slope, population, population density, and distance to urban areas as the most important variables for the LUCC ; in northern Iran, in addition to the variables above mentioned, Aroengbinang [88] also highlighted the distance from drainage network; and, in Alt Emporda county (northeast of Spain), Serra et al. [42] referred to economic factors (e.g., agricultural subsidies) to promote irrigate crops versus dry crops, or social factors (e.g., senior agrarian holders) to justify the cultivated areas reduction.In Portugal, where the forest and agricultural land predominates, these LUCC driving forces also are important, i.e., the allocation of subsidies (to forest and agriculture), agricultural abandonment due to social and economic factors (e.g., low wages, emigration of the population, poor agricultural and forestry production, reduced competitiveness in Portuguese and international markets, etc.), the construction of large infrastructures such as dams, among other factors.In Ribeiro et al. [77] already referred to some of these factors as driving forces boosters of LUCC in Portugal (regional level), especially in agricultural land.
The landscape is the identity of a territory and it contributes to human well-being promoting important cultural, ecological, environmental and social functions [89].In this context, landscape should be the target of constant evaluations due to political, socioeconomic and territorial decisions that are taken at the macro scale (e.g., distribution of social funds and agricultural development strategies and guidelines for the regional territory) with direct implications in LUC and natural ecosystems.Portugal has benefited from European Union (EU) funds for agricultural, forestry, transportation infrastructure, among others, allowing the development of large projects that culminated in LUCC in very large areas, especially in the Alentejo region [90,91].In this sense, the results presented are important to estimate the potential implication of economic decisions to the economical (corporations and state) and political strategy to develop profitable economic sectors (e.g., Eucalyptus forest, vineyards, olive trees, rice crops).These decisions have inevitable consequences on the employment, electricity and energy consumption, volume of exports and gross value added.In addition, studies about the driving forces that control LUCC and possible implications on society, economy, environment and territory are important to the implementation of new land management guidelines and develop more effective strategies for ecosystems protection [92].
This research also shows that the LUCC projections are a challenging task because they involve several dynamic factors that can affect the LUCC in short periods.In addition, the example of the Alqueva dam construction caused deep LUCC in a short period, which would not have been predicted in a projection made with the 1995 LUC data.
The obtained results are also important to the LUC sustainability of the Portuguese territory, because the forest ecosystems have been continuously degraded by forest fires and human activities (e.g., disforestation).The forest surrounding the principal urban areas has diminished in each decade and the ecological pressure in these areas is high.

Figure 3 .
Figure 3. Methodological scheme used in the LUCC and driving forces and implications assessment.

Figure 3 .
Figure 3. Methodological scheme used in the LUCC and driving forces and implications assessment.

Figure 3 .
Figure 3. Methodological scheme used in the LUCC and driving forces and implications assessment.

Figure 6 .
Figure 6.Principal components results of the LUCC classification and driving forces (A) and implications (B) variables by NUTS II for the complete period.N, North; C, Centre; L, Lisbon; AL, Alentejo; and AG, Algarve; and the numbers after the letter correspond to the year: 95, 1995; 07, 2007; and 10, 2010.

Figure 6 .
Figure 6.Principal components results of the LUCC classification and driving forces (A) and implications (B) variables by NUTS II for the complete period.N, North; C, Centre; L, Lisbon; AL, Alentejo; and AG, Algarve; and the numbers after the letter correspond to the year: 95, 1995; 07, 2007; and 10, 2010.

2. Study Area The
study area is Portugal (only Continental Portugal, 88,962.50km 2 ), a European country divided into five NUTS II (Nomenclature of Territorial Units for Statistics, level II): North, Centre, Lisbon, Alentejo and Algarve.

Table 1 .
Area (%) of Portugal by LUC types in NUTS II for the years 1995, 2007 and 2010.

Table 2 .
Socioeconomic and environmental statistical data considered in analysis.

Table 5 .
Correlation coefficients between LUCC area variations (absolute and relative changes) between all NUTS II combinations (significance level p < 0.05).

Table 6 .
Correlation coefficients average values between the LUCC areas of the NUTS and the driving forces and implications variables for the complete period (significance level p < 0.05).