What Is the Influence of the Planning Framework on the Land Use Change Trajectories? Photointerpretation Analysis in the 1958–2011 Period for a Medium/small Sized City

Medium/small sized cities create a polycentric urban system representing the backbone of their territory, characterized by profound changes on land use. The objective of this study is to evaluate the spatial and temporal dynamics of land use in a medium/small city in Portugal during the period 1958–2011, and the relation with the planning framework. The assessment involves land use changes recognition in the period 1958–2011, calculation of the stability grade indicator, and the losses and gains between classes. The rate of artificialization and its relation with the planning framework was also evaluated. The results for the city of Viseu showed a main decrease in annual crops and arable lands, with an increase of the continuous and discontinuous urban fabric. The changes are systematic transitions, marked by planning framework and its typology, objectives and scale enforcement. A plan addressing the city encouraged land use changes in the fringe, while a municipal Master Plan determined sequent transformation in all areas. The land use changes and the artificialization processes enabled the identification of three periods where planning framework was forthright. The study underlines the importance of a planning framework for medium/small sized cities for urban sprawl and artificialization processes.


Introduction
Land changes in urban areas, especially the conversion of cropland and forest land to urban uses, is one of the most important forms of global environmental changes [1].Cities of small and medium scale are frequently characterized by profound changes in land use, a heterogeneous occupation, and a high level of infrastructure and equipment, where hazardous processes, natural and technological, have cumulative presence with economically and socially relevant damages and losses [2].The relevance of these cities is also pointed out by some authors [3][4][5][6], presenting a complex social distribution and organization, with the increase of territorial vulnerability.However, as stated by the European Territorial Strategy and the Territorial Agenda of European Union [7,8], these cities create a polycentric urban system, where their interdependency forms are especially important for rural regions and represent the backbone of their territory because of the ability to generate jobs and permit the dissemination of certain urban functions [9].
Different studies consider land use and the urban expansion as a combined outcome of societal needs, supported by planning preparation actions [10][11][12][13].Understanding long-term dynamics of the landscape and transitions of land use is fundamental to dimension the past and future effects and for the recognition of present landscape structure [14][15][16][17][18].
The use of satellite imagery and aerial photographs for interpretation has made it possible to integrate sets of images from several observation periods obtained from different sources and using different scales.The reconstitution of the transformation path using those images, in accordance with uniform classification principles and field validation capabilities, has been explored in other studies to reveal the dynamics of land-use transformation [20,[24][25][26][27].
The analysis of changes on land use and occupation in the territory, with special focus in the processes from small and medium cities, allows the evaluation of consistency in urban expansion scenarios, the vulnerability of ecosystems, and the maintenance of the distinctive morphological and socio-cultural features [1,4,21,26,[28][29][30].
Several authors present indicators to monitor urban expansion and temporal dynamics which are major triggering factors in land use change processes [31][32][33].Different research groups point sets of indicators responsible for the characteristics of landscapes and land use change [34][35][36][37], namely related with the planning framework [13,[38][39][40].
The main objective of this study is to evaluate the spatial and temporal dynamics of land use over time in an inland medium/small city of Portugal during the period 1958-2011, and the relation with the planning framework.The research questions we want to answer are: What are the changes in land use in the period 1958-2011, using five moments of analysis?What are the stability grade and the dynamics of losses and gains between classes of land use?What is the zoning model for land use change and artificialization process?What is the relation between land use trajectories and the planning framework?

Study Area
The study area is located in the Central Region of Portugal and it is a part of the municipal area of Viseu (Figure 1).The county of Viseu is a district capital and centre of the NUT III (Dão-Lafões), which presents one of the highest population consolidations in the inland of Portugal, polarizing several neighbouring municipalities.In morphological terms, the area is surrounded by mountains where the regional hypsometric values can reach 1043 m about 17 km away from the city centre.The study area is characterized by hills sometimes marked by river valleys, with altitudes values ranging from 300 to 600 m, approximately.
The geology of the area consists essentially of hercinic granite rocks, with small expressions of metasedimentary rocks and sedimentary and alluvial deposits preserved in tectonic depressions and related to the large rivers that cross the region [41].
The region has a Mediterranean climate with a well-defined dry season, although short (July and August).The average annual temperatures vary between 10 °C and 15 °C and the annual precipitation ranges from 400 to 2400 mm, with unstable winter atmospheric conditions, due to the influence of mountains, with long lasting periods of concentrated rainfall [42,43].
The land occupation is a further evidence of the heterogeneity of the municipality of Viseu.It is dominated by forest areas, essentially wild pine followed by agricultural areas, uncultivated areas, and artificial areas.The most distant parishes of the city centre have a higher forest area [42].
In biogeographical terms, the municipality of Viseu is located in the Eurosiberian Region.More precisely, the study area represents the Galician-Portuguese Sector, specifically in the Miniense Subsector, that is characterized by the existence of seasoned oak of Rusco aculeate-Quercetum roboris quercetosum suberis and furze (Ulex europaeus and Ulex micranthus) [44].
Regarding human occupation, the municipality of Viseu has 99,274 inhabitants according to the 2011 Census [45]. Figure 2 shows the evolution, in the period 1950-2011, of the resident population in the municipal area and in the parishes that are within the study area.It presents a small variation by the year 1970, with a sequent population increase (24%), where the growth of the municipal population is essentially in the study area.Figure 2 presents some municipal demographic and socio-economic indicators which demonstrate the progressive number of the resident population and the crescent number of family households.The population increase between 1970 and 1990 was influenced by the social integration of returned Portuguese, both from the former colonies and emigrants in Europe.The population and household increase after 1990 expresses the general improvement of well-being after joining the European community and concentration on the city from residents from surrounding rural municipalities.These socio-economic dynamics are expressed by the sharp increase in the number of electricity consumers, accompanying the increase of the population active in the tertiary sector at the same time reducing the importance of the primary sector.Industrialization, in the seventies, and the tertiary activities in the eighties and nineties, were influenced by road infrastructure construction, with the correspondent accessibility improvement.
Those aspects reinforced the increment of centrality and attractiveness of the town, which is expressed by a spatial urban growth, and achieving the characteristics of a medium-sized city in the national context [47], and a small-sized city according the EC definition [48].
To evaluate the land use and occupation changes, and the progressive urban expansion and infrastructure, which reflect the planning framework, a study area was selected, presented in Figure 1.It is centred in the City Council of Viseu, neighbour of the XV ancient town, and has a circular shape with 4 km radius, with 50.26 km 2 of area.The circle represents or partially covers nine parishes of the municipal territory, including the urban core and the peri-urban areas.
As mentioned by J. Almeida [47], the urban model of the city of Viseu is defined by a dominant and compact centre and develops in "stellar" form, growing the city quickly to the outside along the pre-existent road network, which provides the appearance of peri-urban areas of some size, later connected by concentric pathways [48,49].For S. Almeida [49] the urban fabric development has been conditioned by urban structures, namely the encircling roads, but also the by collective equipment's location.
In the last 50 years two former master plans influenced urban expansion, the infrastructure network, and the constrained areas protection.A city urbanization plan was approved in 1952 and a Municipal master plan was published in 1995, which is currently being replaced.Between these two instruments there were partial plans, such as the Hidroprojecto Urban Plan (1971) and Macroplan Urban Plan (1983) that has never been approved, but with informal application for new urbanization areas, or the Polis Program (2000)-Urban Rehabilitation Program and Environmental Improvement of Cities.In Figure 3 some of the cartographic extracts of these plans are represented, emphasizing that by 1983 the plans were partial, focusing mainly in the urban core of the city, and in 1995 acquired the municipal level, which has been complemented by plans further detailed (layout plans, urbanization plans, and a rehabilitation plan), and also the 2013 Municipal Master Plan under approval.

Land Cover Data
The analysis of the land use and occupation change was supported by a set of five photographic mosaics from different periods, between 1958 and 2011, as described in Table 1.According to the initial characteristics of the images, they were digitalized, georeferenced, processed and analyzed using the ArcGIS 10.2 software (ESRI ® ).To adjust the different imagery scale three procedures were carried out (single operator, minimum unit of area, and backward analysis).Due to differing contrast levels, the 1958, 1974, and 1985 images had to be cross-referenced with other topographic representations.In order to integrate images from different sources consistently, a supervised analysis was performed using fieldwork as a control.A Datum_73_Hayford_Gauss_IGeoE was used as the geographic coordinate system in the process of georeferencing.In another phase, based on the categories of land use and occupation of the Corine Land Cover [50], a classification system has been developed and identified by classes organized according to numerical codes (Table 2) [51,52].Then, the classification of polygons of land use from each set of images was performed and incorporated into the geographic information system generated [17,18,53] wherein all polygons were classified using a minimum unit of analysis with 1 hectare and supported on a 1:10,000 scale.

Spatial and Temporal Changes in Land Use and Annual Rate of Artificialization
In order to evaluate the dynamics of land use and occupation a map to map comparison method was used involving the set of successive images, with cross-references to define the transitions of land use and land cover and also the measurement of the stability grade (SG) in each time period.
The maps created for each year were used to produce matrices for the periods 1958-1974, 1974-1985, 1985-1995, 1995-2011, and 1958-2011.In the Table 3, the calculation algorithms for the different indicators considered are presented.

Designation Equation Legend
Change in land cover CHi = (pi.
CHi-change in land cover for the period in column i pi.-the column total for the grid cells p.i-row total for the grid cells in the same category i Conversion areas corresponding to gains and losses for a given year AS-evaluate the speed of urbanization in the study area TA n + i -total study area to be calculated at the time point i + n UA n + i and UA i -surfaces with urban fabric and infrastrutures in the target unit at time i + n and i n-number of years in each time period The Equation (1), used in other studies [13,36,54] was applied in order to express the land use transition dynamics between two periods.
The conversion areas corresponding to gains and losses for a given year, relative to the year under comparison were calculated in relation to the total land-use type using Equation (2).
To express the total area that has not experienced a possible transition to a different category of land use and occupation, the indicator stability grade (SG) was defined as shown in Equation ( 3).
An annual rate of artificialization (AS) was calculated using Equation ( 4), expressing the unceasing increase in urbanization, namely related with peri-urbanization.This indicator combines urban transformation processes that express soil sealing and vegetation cover change, from the continuous and discontinuous urban fabric, industrial or commercial units, road and rail networks, mineral extraction areas, dump site, construction sites, green urban areas, and sport and leisure facilities.
In order to understand if the major changes in land use and occupation in the study area occur in the central area or in the periphery, and according to other studies [22,55,56], it was decided to divide the area into 4 different sections with increasing distances from the center point (city hall), as shown in Figure 4.For this purpose, four buffers were defined with distances of 1, 2, 3, and 4 km respectively, relative to city council point.As a result, a central circle and three rings with increasing distances from the central point were originated.The focus of this analysis was to establish a model of land use change dynamics according with the spatial planning framework which presents different scales and spatial enforcement.Thereafter the stability grade was calculated for the four sections in the five different periods of analysis.
To analyze the relationship between land use change and planning processes in the city of Viseu, different instruments of planning were compiled, between 1952 and 2008.For each of them there was an analysis of the objectives, the type, and area of enforcement.
For these plans, a content analysis of the planning documents and associated reports was made, supported by methods [57,58].The special protection legal regimes which present constraints to planning and land use transformations were also identified.
Table 4 presents the planning framework for the period of analysis, and also indicates the previous plans and the recent master plan.As can be observed, during the period of analysis the first initial plans had a central city enforcement, and were responsible for creating new residential areas, infrastructures and an encircling road network.These plans were supported by local plans for urban occupation proposal.The plans with municipal characteristics and enforcement were approved in 1995, and supported by a regulatory agricultural and ecological protection regime for non-urban areas.Yet, in addition to these plans, some urbanization plans, layout plans and a rehabilitation plan emerged.

Temporal Patterns for Land Use Change in the Study Area
Using the methods described a systematic analysis and classification of land use and occupation was produced for the five moments (1958,1974,1985,1995,2011) in the study area.
Table 5 shows the classes applied to the use and occupation of land, the number of polygons of each class and the respective area and percentage.Twenty six different classes of land use and occupation have been identified in the circle area under the analysis (Figure 5).The images only presented five classes with episodic representation (mineral extraction areas, dump site, construction sites, burnt areas, and water bodies).As observable, over the analysis period there was a steady increase in the urban area, namely on the discontinuous urban fabric.This process was still associated with increasing complex cultivation patterns.Since 1985 there has been a considerable decrease in areas with annual crops associated with permanent crop and pasture use.The display also shows the maintenance of forest areas, with constant values above 30% of the study area.The cartographic outputs present in Figure 6 show differentiated changes in the different sectors of the area with a variety of trajectories.The analysis indicates some systematic transitions involving an increase in urban areas and a decrease in agricultural areas, as well as the transformation of forest types and the episodic occurrence of wildfires.The analysis also reveals a continuous increase in the number of polygons between 1958 and 2011, which demonstrates the fragmentation of the landscape, representing a fast expansion process with new corridors [4] with impacts on natural systems and landscape homogeneity [59,60].In the study area, and having as focal point the city center-center of the considered circle, two different transformations can be observed.Firstly, there is an increase of the urban areas with a continuous fabric, discontinuous fabric, industrial or commercial units, road and rail networks, green urban areas, and sport and leisure facilities, as well as complex cultivation patterns.Secondly, a systematic decrease in the annual crops associated with permanent crops, and a routing on the forest occupation.
The analysis of land use change shows some temporal transitional processes (e.g., forest and semi-natural areas, forest and burnt areas) representing casuistic losses of forest space.
The results also show the growth of pattern complexity and reveal the general increase of the artificial areas from the central part of the study area to the periphery, contrasting with a continuous decrease of the agricultural areas.These results illustrate the peri-urbanization and rururbanization processes with the increase of classes of urban fabric and complex cultivation patterns.

Dynamic Analysis of Changes in Land Use in the Study Area
The transformation matrices for 1958-1974, 1974-1985, 1985-1995, 1995-2011, and 1958-2011 allowed for a detailed study of the dynamics of land use and occupation in five periods of analysis.For each period of analysis a transformation matrix was generated, and the respective stability grade calculated.
The results obtained in the matrices allowed the creation of a scheme (Figure 7), illustrating the dynamics of land use in the different periods and highlighting the visual trajectory of transformation between classes.The scale used to achieve the results presented in the figure, is comprised of three ranges (0.2 to 0.5 km 2 ; 0.5 to 1 km 2 , and > 1km 2 ) that represent the higher or lower amount of area transferred between classes from period to period.The results show both the dominant dynamics and the differentiated processes for the various periods analyzed, highlighting specific gains and losses within the overall transformation and indicates a continuous trajectory of urbanization for the overall period , namely supported by the annual crops.The scheme also demonstrated a dynamic involving agricultural areas, with a loss in the annual crops and pastures and a gain in the complex cultivation patterns and agriculture with natural spaces and agro-forestry areas.The period between 1958-1974 shows an increase in the forest, with consequent evolution to woods and degraded forest areas, with the exception of the casuistic transition associated with the burnt areas.
In order to express the percentage of conversions in relation to the total land cover type, an internal transition matrix was created, for the period 1958-2011.Table 6 detailed the gains and losses for the period under analysis (1958-2011) for the 26 classes of land use and occupation.This matrix was obtained by cross-referencing the maps of 1958 and 2011 and the results show the units of change for each land use class and indicate the (%) rate of change.The results in general show an increase in the urban fabric and infra-structures associated, supported by the arable lands, fruit trees and olive groves, pastures, annual crops and agriculture with natural spaces and agro-forestry areas.In the agricultural areas, it is possible to verify an increase in the complex cultivation patterns (360.22%),opposed to the rest of the classes that present losses.Gains in broad-leaved, mixed forest and woods should also be noted.The results also reflect the decrease in coniferous forest, degraded forest areas, and sparsely vegetated areas.The matrix demonstrates the low stability grade (SG) for the period 1958-2011 (20.61%) reflecting the changes and dynamics between the different land use classes.This means that 79.39% of the land cover has experienced changes, showing intensive land use dynamics during the study period.The final results show the general gains and losses for the discontinuous urban fabric from annual crops and coniferous forest (positive values), generating continuous urban fabric, industrial and commercial units and green urban areas (negative values).The conversion areas involving agriculture with natural spaces and agro-forestry areas display substantial losses to discontinuous urban fabric, complex cultivation patterns, and mixed forest.

Model for Land Use Change and Artificialization Processes
Expressing the land use change dynamics in the study area, Table 7 presents the aggregation of five classes of land use, scaling up the classes determined over time the process of artificialization.The results point, for each period, to the gain of artificial areas with the percentage that remains artificial and the one added from other previous land use.This expresses soil sealing, vegetation cover change, and antrophic actions with soil consummation for urban fabric, industrial or commercial units, road and rail networks, mineral extraction areas, dump sites, construction sites, green urban areas and sport and leisure facilities.As can be observed in Figure 6 and Table 7, the biggest transformations happened from the agricultural and shrub vegetation areas.The forest conversion to artificial area is successive, with a major rate in the period 1995-2010.When analyzing the evolution for the center zone and the three rings it is noted that the artificialization of agricultural areas decreases from the center to the periphery, the decrease from 1985 in center area and on ring 1 being especially evident.The artificialization of forest areas occurred rapidly since 1974 in the center area, verifying that there is a progressive use on the different rings, but more noticeably from 1995 onwards.
The use of shrub areas and sparse vegetation for artificialization processes has a tendency of gradual use over time, more evident in the center than at the periphery, but also expressing temporal and spatial peak usage.
This process of artificialization in general involves a sequence of land use change and occupation from: (1) agricultural areas to forest; (2) forest to shrub or sparsely vegetation areas; (3) shrub or sparse vegetation areas to artificial areas, is recognized as common dynamic in the peri-urbanization processes.
The results also emphasize the changes to artificialized areas with an urban fabric and infrastructures associated, namely during the periods of 1985-1995 and 1995-2011.In order to realize the major changes in the land use and occupation and the relation with the spatial variation from the central area or in the periphery, the partial stability grade (SG) was calculated for the 4 sections considered, in the five different periods of analysis.
Looking to the stability grade study into the four sectors, the dynamics of land use also reveal distinctive frames, as we can see in Figure 8, which shows the results of the stability grade of each sector in the different periods of time analyzed.The Figure 8 analysis shows a similar trend in the periods 1958-1974 and 1974-1985 for stability grade, where the SG decreases in the following order: central sector, ring 3, ring 1 and ring 2. This means that larger changes occur in the intermediate sectors (ring 1 and ring 2) in relation to the proximity to the city center, but also expressing the city planning framework.This made a tight transformation process for the excluded plan areas in the fringe of the central area possible, and potentiated the peri-urbanization and rururbanization processes.In the period 1985-1995 the trend of SG is different, presenting a continuous decrease from ring 3 to the central sector, showing that the dynamic of land use is more intensive in the urban core of the study area with the influence of the Hidroprojecto Urban Plan (1971) and Macroplan Urban Plan (1983) that increased the artificialization.In the last period (1995-2011), with a municipal master plan that includes all the sectors in the planning framework, the trend is the inverse of the period 1985-1995, with a continuous decrease from the center to the ring 3.In this period the high value of the SG in the central sector (63.06%) is tangible compared with the rings 1, 2 and 3, which means that larger changes occurred in the peripheral zones of urban core.The stabilization of the urban core and consequent intensive dynamics of land use in peripheral zones are in accordance with the directives of the Master Plan approved in 1995.
When analyzing the overall trend for the period 1958-2011, it can be observed that the periurbanization and rururbanization processes that occurred outside the boundaries of the planning framework trace the whole evolution of land use and occupation, corresponding to both ring 1 and ring 2 low stability grades (18.60% and 17.22%, respectively).

Land Use Trajectories and Planning Framework
The focus of this analysis is to deepen the model of land use change dynamics according to the spatial planning framework which presents different scales and spatial enforcement.
The graphic present in Figure 9 also show a relationship between the planning process and the rate of artificialization.The results show that the informal partial Plans (Hidroprojecto Urban Plan and Macroplan Urban Plan) influence the increase of the artificialization rate, as they respond to pressures to the occupation of the territory for new residential or service areas.These plans are directly related to the peri-urbanization processes and the definition of new road infrastructures.With the approval of the Master Plan in 1995 there has been a decrease in the artificialization rate.This plan fits all land use and occupation processes for the whole county and is a regulatory plan supported by two restriction regimes (agriculture-RAN and ecological-REN), leading to a decrease in the annual rate of artificialization, despite the increase in continuous and discontinuous urban areas.
The increasing distances from the center point (city hall), as shown in Figure 4, and the different planning framework typology, objectives, and enforcement in the study area produced during time different dynamics on land use change.
A zoning model for land use change and the artificialization processes can be summarised as: (a) In the period 1958-1985 the major transformation occurred in the fringe of the central area (ring 1 and ring 2) which are outer areas for the central urban planning enforcement.The agricultural areas are the principal source of the artificialization process, with the subsequent use of the forest areas; (b) For the period 1985-1995 the stability grade is directly related to the distance for the city council central point, which reflects the influence of the non-formalised central urban plan, and the encircling road consolidation.This period also reflects the intensification of the artificialization processes namely using the shrub vegetation areas and sparsely vegetated areas, which represent abandoned areas or areas with reserve urban qualification.
The period 1995-2011 expresses the influence of a master plan with municipal enforcement.The degree of stability intensely downloaded to the central area now operates in reverse to the different rings.So, there is a major change towards the peripheral ring, which expresses a process of municipal urban sprawl with periurbanization and rururbanization processes.Observing the artificialization areas, the gains are mainly from the agricultural areas, forest areas, shrub vegetation areas, and sparsely vegetated areas.
These processes namely represent a deep transformation on rings 1 and 2, supported by partial plans and the road infrastructure design and construction, which often outweighs the planning processes by economic or political decision propose.Those artificialized areas present considerable gains from annual crops, arable lands, and pastures, defining those classes as one of the most important in the trajectories of land use and occupation in the study area, but associated with systematic transitions involving the forest, shrub vegetation, and sparsely vegetated areas.The casuistic transformation with land use classes resulting from forest fires, the appearance of mineral extraction areas, and waste disposal areas do not reflect the planning framework, in opposition to the modification on the stream waters margins related with legal constrains or rehabilitation plan.

Discussion and Conclusions
The five cartographic outputs for changes in land use and occupation underline the main decrease in annual crops and arable lands, the increase in artificial areas, namely the continuous and discontinuous urban fabric, the industrial and commercial units and the road network as well as the complex cultivation patterns and the forest routing involving different typologies.
The data representation demonstrates that urban occupation has increased gradually since 1958.The increase of artificial areas since the 70s results in the follow-up, at first, of industrial growth and afterwards, tertiary sector growth with the overall improvement of well-being [47].
The transformation of the agricultural land to urban areas is evident in this study case, as well in other contexts [5,61,62], and reflecting similar changes in Portugal [13,63].In the area, an increase in landscape fragmentation, accompanied by increased complex cultivation patterns, translates the evolution of the rural areas [64,65] with peri-urbanization and rururbanization processes.
In general terms, there has been a continuous increase in the urban area (from 5.7% to 27.9% of the area), and the stability of the area occupied by agriculture with natural spaces and agro-forestry areas (around 6%).It was also possible to observe a decrease in the annual crops (from 27.5% to 5.9% of the area) and in the coniferous forests (from 16.1% to 5.3% of the area), especially underlined in the 1995 and 2011 images.This represents the general transition of the rural areas, with loss of agricultural areas and forest degradation, and can be clearly associated with small and medium sized city urban expansion.
The results for the study area confirm processes that can be classified as systematic transitions, contrasting with some land use changes marked by casuistic trajectories in the case of land use classes resulting from forest fires or the appearance of mineral extraction areas and waste disposal areas.
In a broad sense there has been a continuous increase in the artificial areas, with a general transition of the rural areas that can be clearly associated with the small and medium sized cities urban expansion.For this study area, but also referred to other urban contexts [13,56,66,67] a general sequence of land use change and occupation from: (1) agricultural areas to forest; (2) forest to shrub or sparsely vegetated areas; (3) shrub or sparsely vegetated areas to artificial areas, is recognized as common dynamic in the artificialization processes.
This study stressed the importance of monitoring the stability grade to identify the sequences of major changes in land use.This indicator points to a low level of stabilization, 20.61%, which reflects the intensive dynamics of the changing classes in the area.The temporal dynamics between the 26 classes identified also reveals distinct transformation phases with loss and gain conversions, justifiable by territorial planning framework [47], which is also noted in other geographical contexts [13,63,68].Some of the systematic transitions are clearly marked by the planning process, involving partial Plans (city urban, urbanization or layout plans), Master Plans, and Regulatory Protection Regimes.
This does point out that a plan that addresses only the city encourages the land use change in the fringe, while the municipal Master Plan determines the transformation processes in all areas, revealing systematic transformations.
The results demonstrated that the transformation processes within the study area are supported by a history of urban occupation, including the discontinuous fabric and encircling road network as well as the loss of agricultural activities to the emergence of new housing residences and equipments.The results have also shown the urban planning framework as a driver on land use and occupation changes, differencing and marking the territory and its systematic transformations, according to the type, objectives and scale, as pointed by [40].
The Figure 10 summarized the two forces that act on the land use change: the territorial dynamics and the planning framework.
The zoning model approach for land use change and the artificialization processes enables the identification of three periods where the planning framework was forthright.The results show different dynamics on land use change, to the central areas and the fringe, directly related with planning scope and the scale enforcement.This conclusion highlighted the importance of the planning framework for medium/small sized cities on processes of urban sprawl and the periurbanization and rururbanization processes.The study also underlined the relevance of road infrastructure design and construction on zoning land use changes and enabled territorial dynamics with systematic transformations.According to Antrop [4] the evolution can translate an expanding changing pattern for the period 1958-1985, an exploding changing for the period 1985-1995 corresponding to a several urbanization and layout plans on the central area and the close rings, and a broad ring development for the period 1995-2011 promoted by the municipal Master Plan.
The results obtained show temporal and spatial processes of land use and occupation changes that affect the medium and small sized cities as recognized from different authors [1,3,24,[69][70][71][72], and underlining the national cycles of evolution on the processes of urbanization and infrastructures construction that spread the artificialization in different ways through space [31,40,66,73,74].
These achievements support the discussion about the influence of the national planning system in relation to the spatial inequalities and as drivers on the land use change and management processes [75,76].
Land use change studies of medium/small cities and there specific patterns and processes of artificialization are needed.Planning the territory in a grounded and creative way requires a deep knowledge about the systematic and casuistic processes of land use and occupation and the different forcers, stressing the goals pointed out by different authors [33,77,78].
Giving consistent attention to historical planning framework ascendance on the land use change processes could support new strategies and regulatory regimes to manage and maintain the landscape heterogeneity and the multifunctional land use, avoiding conflicts and promoting recognition from actors.
Reshaping or rebuilding the medium/small sized cities is a constant challenge for planner and local inhabitants, balancing new functionalities and infrastructures with diversity values and ecological corridors where the land use dynamics are drivers to create new functions and offer a basis for planning and decision making.

Figure 1 .
Figure 1.Location of the study area.

Figure 2 .
Figure 2. Demographic and socio-economic indicators: (a) municipal and study area resident population; (b) number of family households; (c) number of electricity consumers; (d) employed population by activity sectors, according the PORDATA database [46].
percentage by type j in the total conversion of category row i pj,i and pi,j-express the individual entry in the change matrix Stability grade SG = [(C ix A 2…5 − C ix A 1 )/TA] x 100 (3) SG-indicator of the consistency and represent the global stability of the classes of land use in the year A1 for the year A 2…5 C ix A 2…5 -sum of the areas of the different classes (ix) at the subsequent time point C ix A 1 -sum of the areas of the different classes (ix) at time point 1 TA-total area studied (km 2 ) Annual rate of artificialization of surfaces

Figure 4 .
Figure 4. Subdivision of the study area in four different sections to calculate the stability grade in each one.

Figure 5 .
Figure 5. Representations of subclasses of land use from 1958 to 2011 in the study area.

Figure 6 .
Figure 6.Land use and occupation in the area studied, based on image classifications by year.

Figure 7 .
Figure 7. Dynamic transformations of land use for the study area and the stability grade (SG) for each period.CUF: Continuous urban fabric; DUF: Discontinuous urban fabric; ICU: Industrial or commercial units; RRN: Road and rail networks; MEA: Mineral extraction areas; DS: Dump site; CS: Construction sites; GUA: Green urban areas; SLF: Sport and leisure facilities; AL: Arable land; V: Vineyards; FT/OL: Fruit trees and olive groves; P: Pastures; AC: Annual crops associated with permanent crops; CCP: Complex cultivation patterns; ANS/AF: Agriculture with natural spaces and agro-forestry areas ; BF: Broad-leaved forest; CF: Coniferous forest; MF: Mixed forest; W: Woods; DFA: Degraded forest areas; BR: Bare rock; SVA: Sparsely vegetated areas; BA: Burnt areas; WC: Water courses; WB: Water bodies.

Figure 8 .
Figure 8. Stability grade for the 4 sections of the study area in the periods of analysis.

Figure 9 .
Figure 9. Annual rate of the artificialization of surfaces (AS) and the respective mean for the period 1958-2011 and the relation with the planning framework.

Figure 10 .
Figure 10.Forcers acting on land use and occupation changes, promoting the artificialization processes.

Table 1 .
Data source characteristics.

Table 5 .
Number of patches, areas, and percentages for each land use and occupation in the study area, according to the Corine Land Cover classification.

Table 6 .
Transformations in land use and occupation in the study area (%).

Table 7 .
Transformations in land use and occupation in the study area (%).