MAPPING OF LAND UNITS AND LAND CAPABILITY CLASSIFICATION IN PORTUGAL: THE CASE OF THE MUNICIPALITY OF LOURINHÃ

Land units can be an advantageous tool in landscape planning, in land suitability and land capability classification. it is vital to deepen the knowledge on their delimitation procedures, particularly in areas where available data is scarce and presents a certain complexity. The municipality of Lourinhã, based on sedimentary rocks, was chosen as a case study. Land units were mapped using lithological units (1:50  000), landform types and soil unit (1:25  000) information. several inadequacies were observed both in the soil and soil capability mapsi. There was evidence of a strong information deficiency on the basic characteristics of soil units and the presence of rock outcrops was not revealed in such maps. from the lithological units and landforms, 33 basic physiographic units were delimited, and their association with soil units resulted in an equal number of land units, afterwards genera lized. The delimited land units have proven to be suitable for a land capability classification at the municipality scale, revealing great potential to identify the nature and severity of the major land (and soil) degradation risks. also, they prove to be useful for improvement of the soil capability classification and to standardize the definition of the national agricultural reserve by using identical delimitation criteria.


i. intrODUCtiOn
Land units are an essential tool for several landscape planning applications as they express a "tract of land that is ecologically homogeneous at the scale level concerned" (Zonneveld, 1989, p. 68). Land units have also been used as a mapping tool and provide a basis for land suitability and land capability assessment (Zonneveld, 1989, p. 68;OeH, 2012). recent studies demonstrate that land units can be used as a main concept in the resolution of landscape planning issues. for instance, Blasi, Zavattero, Marignani, smiraglia, Copiz, rosati, and Del Vico (2008) used the land unit approach to design an ecological network of sites in the rome Province (Lazio, italy), in order to fulfill the main ecological needs of species, community and ecosystems. also, Capotorti, Guida, siervo, smiraglia, and Blasi (2012) used land units as a tool to describe the ecological classification of land in italy and show how these units can act as reliable frameworks to implement national conservation strategies. in addition, the land capability mapping in tasmania was developed using land units, that is, homogeneous units in terms of climate, relief, lithology and soils (Grose, 1999) in Portugal, the land unit concept has been only partially developed in soil surveys associated with the evaluation of land suitability. in particular, agroconsultores and Coba (1991) in the elaboration of the Portugal northeast soil map delimited natural regions corresponding to large landscape units combining physiographic features (that is, landform types), climate, vegetation and land use, thus creating homogenous zones approaching the land unit concept. Later, agroconsultores and Geometral (1995) also produced the entre--Douro and Minho soil and land suitability maps at 1:100 000 scale, delimiting large morpho -climatic units, relatively homogenous as to climate, lithology, geomorphology, topographical conditions, vegetation and land use. Meanwhile, abreu, Pinto, & Oliveira (2004) defined landscape units, at the national level, as areas where the landscape presents a specific pattern associated with a particular character, defined by geomorphology, lithology, soils, land use, farm dimension, human settlement pattern, climate, proximity of the sea and presence of important structures and infrastructures. although the delimitation procedure of landscape units (abreu et al., 2004) is not quite clear, the concept adjusts itself fairly to the land definition as expressed by faO (2007).
Despite the scarce applications of the land unit concept in Portugal, the legislation regulating the national agricultural reserve (Reserva Agricola Nacional, ran) recognizes the need for the delimitation of homogenous areas defining land units and pointing out that "the land units that present a high or moderate suitability for agriculture", are to be integrated in the ran, after being classified "according to the classification methods of the land suitability classification recommended by the food and agriculture Organization (faO)" (articles 8 th and 6 th of Decree Law nº 73/2009, 31 st March, respectively). as large areas of Portugal still requires classification, such law foresees that, when it is alternatively possible to use the faO classification, it is also possible to affiliate into ran "the areas with soil capability classes a, B and Ch", or "alluvial or colluvial areas" or else areas where these classes are largely represented. Besides the absence of the classification rec-106 ommended by faO (1976,2007) in the most part of the country, the soil capability map shows deficiencies, namely in the littoral central region, where soil capability units were not adequately framed with the general system of the Official soil service (srOa, 1973;iDrHa & sPCs, 2005). it should be emphasized that "soil capability map" and soil "capability classification" are a translation of the Portuguese terms carta de capacidade de uso do solo and classificação da capacidade de uso do solo (srOa, 1972), which approximately correspond to the land capability classification (Klingebiel & Montgomery, 1961). although the aim was to develop a classification based on the land capability classification (by the United states soil Conservation service), the soil capability classification (in Portugal) did not consider the climate and was based on less detailed soil information compared to that used in the Usa (srOa, 1972).
since the land unit concept is rarely used in Portugal and is considered in the legislation that regulates ran areas, it is therefore crucial to implement the methodology regarding land units mapping and to evaluate whether ran areas can be assessed, especially where the information is scarce and presents great complexity. in this context, a study was developed aiming to (i) analyze the land unit concept applicability, (ii) develop a methodology in areas where the existing data is considered insufficient, (iii) evaluate how much land units are in accordance with the existing soil capability units, and (iv) assess whether the land unit delimitation can lead to a unifying methodology for the ran delimitation. as a case study, the municipality of Lourinhã (an area based on sedimentary rocks) was chosen, since it is used intensively for agricultural purposes.
ii. MetHODOLOGY

Study area
The study was carried out in the municipality of Lourinhã located in Portugal's West Coast, part of the Lisbon district ( fig. 1). it comprises an area of about 147km 2 with a population of ca. 25 000 inhabitants. The majority of the territory is used for temporary crops ( fig. 1); thence the landscape is profoundly influenced by agricultural activity ( fig. 1), especially horticulture. The area exhibits mild temperatures, high relative humidity, and low frost risk in winter. according to the climate based on the closest weather stations and/or pluviometric stations (reis & Gonçalves, 1981), for the period of 1951/1981, the mean annual precipitation (MaP) in Lourinhã varied from 600 to 700mm, since it was 567mm in Cabo Carvoeiro, 608mm in Caldas da rainha and 780mm in torres Vedras. The daily temperature range is relatively narrow, given the small difference between the mean maximum and the mean minimum temperature, either in the warmest month or in the coldest month. in Cabo Carvoeiro, the closest station to the ocean, the range is 4.5ºC in the warmest and 5.3ºC in the coldest month; in Caldas da rainha the range is more pronounced: 8.5ºC in the warmest month and 7.7ºC in the coldest month, probably because of its greater distance from the ocean. The monthly mean relative humidity of the air is high and above 77% in Caldas da rainha, and 81% in Cabo Carvoeiro.

Reorganization of data and land unit delimitation
The procedure for land unit delimitation (see flowchart of fig. 2) was adjusted to the available data for the study area, respecting the fundamental concepts introduced by Zonneveld (1995). after collecting basic available data, field observations were carried out. following the field observations, it was possible to aggregate and dimension the geological formations and landform types. Geological formations were aggregated into lithological units, the alluvial and colluvial areas being taken from the soil map of Portugal, since its scale was more appropriate. Landform analysis was performed according to the Hammond (1954Hammond ( , 1964aHammond ( , 1964b methodology and aggregated into landform types. soil information from the Portugal soil map (using the Portugal soil Classification, PsC) (srOa, 1973;iHDra, 1999) was inferred to the World reference Base for soil resources (WrB, 2006) resulting in eight soil -mapping units. Lithological units and landform types were intersected resulting in basic physiographic units; to each of these units, information regarding soil units, as a qualitative attribute, and rock outcrops were added. The land cover map COs 2007 Level 2 (1:25 000) (Direcção -Geral do território, 2007) and the soil capability map (1:25 000) (DGaDr) were also used to support the analysis of land units. Therefore, the obtained land units were afterwards grouped into generalized land units.

Lithological units
Geological information was gathered from the geological map of Portugal (1:50 000) -sheets 26 -D (Caldas da rainha - Zbyszweski & almeida, 1960aZbyszweski & almeida, , 1960b, 26 -C (Peniche -frança, Zbyszewski, & almeida, 1960a;frança, Zbyszewski, & almeida, 1960b), 30 -a (Lourinhã - Manuppella, antunes, Pais, ramalho, & rey, 1999a;Manuppella, antunes, Pais, ramalho, & rey, 1999b) and 30 -B (Bombarral -Zbyszewski, 1966a, 1966b. The sheets were appended, in order to obtain a geologic map of the study area. These sheets show a distinct legend and content, thus a lithological map was derived from a simplification and aggregation of the existing 22 geological formations, aiming to reach the closest uniformity to lithological units. Lithological units were delineated through the similarity described on top of the lithological formations following the criteria of its topographical position, texture of the detrital or clastic sedimentary rocks and moreover the nature or composition of the rocks (calcareous or not calcareous). for the lithological unit delimitation information from the Portugal soil map regarding the alluvial and colluvial areas was also used. This was due to the fact that the alluvial areas in the geological map did not coincide with those in the soil map. firstly, because the geological map and the soil map have different scales; secondly, because some alluvial areas of the geological map are colluvial areas according to the soil map. Hence, it was preferable to use the alluvial and colluvial areas of the soil map for the land unit delimitation. finally, some information on the characteristics of the soil units (e.g. texture) was also jointly used with the grain size of rocks regarding different geological formations.

Landform units
a digital terrain model (DtM) was prepared by projecting data from shuttle radar topography Mission (srtM) into the etrs89 -tM06 reference system, using bilinear resampling to a 30 meter pixel resolution. This DtM was used to identify the landform types based on the Hammond (1954Hammond ( , 1964aHammond ( , 1964b) methodology, following the procedures specified by Morgan and Lesh (2005). The analysis was implemented in arcGis using esri modelbuilder, as implemented by Mendes (2010). to run the model it was necessary to define a parameter of the focal statistics tool, the pixel radius circular window. That is, the pixel radius that determined the number of cells within a certain circular neighbourhood surrounding each pixel; this tool allows to calculate the percentages needed in the estimation of the various parameters (Morgan & Lesh, 2005). in this case, a 15 pixel radius was used because it was considered to be the most appropriate for the scale of this study.
according to the procedures specified by Morgan and Lesh (2005), the landform analysis is divided in four sub -models, three of them are related with the parameters set by Hammond (1954Hammond ( , 1964aHammond ( , 1964b) as means to determine three different characteristics which are determinant in land form analysis -near -level areas, local relief and profile type -and the fourth sub -model reclassifies the results of the first three sub--models based on the local interpretation according to the field observations. firstly, the percentage of near -level areas is considered -corresponding to areas presenting slope values below 8% -and later a first landform type will be described as "flat to nearly flat sloping", which consists of landforms predominantly dominated by near--level areas. Later, the local relief and profile type sub -models will assist in the delimitation of different landform types as, for instance, the "gently sloping" areas, which can be further divided into "gently sloping lowlands" and "gently sloping uplands", differentiating the surfaces according to its concavity/convexity (Mendes, 2010). finally, the last sub -model will consist of a reclassification of the different areas obtained by the combination of the previously mentioned attributes (near -level areas, local relief and profile type), in order to obtain a final map of aggregated landform types. The geological formations were used to support the landforms delimitation.

Soils
information regarding soil was taken from the Portuguese soil map (1:25 000)sheets 337, 338, 349, 350, 361 e 362 (Direcção -Geral da agricultura e Desenvolvimento rural, DGaDr). The information of the Portugal soil map (srOa, 1973;iHera, 1999) was extrapolated to the WrB (2006) system as shown in table i. in the first column are represented the identification numbers of which the general soil mapping unit obtained as a result. some soil mapping units are marked with ( a ) as they may include small areas of Leptosols and eventually Luvisols and Vertisols; with ( B ) in general can be qualified as calcic; and with ( C ) may include small areas of Luvisols and Vertisols. also, for each of the soil families an indication is given of the various soils phases as (a) agropedic phase; (d) shallow phase; (h) hydromorphic phase; (i) flooded phase; (p) stony phase. it should be emphasized that the delimitation of soil map units in the Portuguese soil map was not supported by the previous delimitation of basic physiographic units (landforms). also, soil units identified in the study area generally reflect their relationship with the local soil parent material, and are not framed in the report of the geological map.   several difficulties emerged with the analysis of soil data. The scale (1:25 000) suggests a detailed soil map, but it presents great spatial variability related to the soil parent material, as designations such as "fine sandstones, clays and claystones" or "fine and coarser interstratified sandstones" are used. Many soil units presented in the soil map of Portugal are distinguished generally by the color, which is not necessarily related to different physical and chemical characteristics. several soil units considered in the soil map of Portugal are absent in the PsC Handbook (srOa, 1973), lacking therefore a general morphological description and a typified characterization (such as thickness and effective soil depth, horizons, texture, concentration of calcium carbonate, complex exchange characteristics...) other than their designation (iHera, 1999). although characteristics of some soil units in the soil map are shown in the PsC Handbook (srOa, 1973), morphological and analytical characterization details for the same units occurring in the area under study are not known. Given these difficulties in establishing reliable inferences, it was only possible to consider general aspects of the relation between the soil parent material and the different soil units presented in the soil map. There were two options: either considering the lithological information from the soil map or the general lithological units obtained from the geological map. The authors opted for the second one, in order to acquire a result that could be referenced to the geological information available but also because the designations of the soil parent material suggest a great heterogeneity and would bring even more subjectivity to the results. Consequently, the criteria used to define the soil -mapping units (reference soil Groups sensu WrB, 2006) were mostly associated with the possible soil profile differentiation, the calcium carbonate occurrence, the nature of limestone rocks and the topographical position. The different soil types were transferred from one system to the other according with the correlations of diagnostic soil horizons established by sousa, Madeira, & Monteiro (2004) and also taking into account the studies carried out by Gomes (1978), Madeira (1986) and Marques (2010). it is noteworthy that the soil information available is not sufficient for an accurate soil classification at first level within the WrB (2006) system.

Land units
in the present case study, the concept of basic physiographic units, a concept already developed in some soil maps of Portugal (agroconsultores & Coba, 1991;agroconsultores & Geometral, 1995), was used to support the delimitation of land units. These units resulted from the intersection of lithological units and landform types; after the intersection it was also necessary to eliminate the small sliver polygons that resulted from the intersection. in this case, the authors opted to eliminate all polygons with an area under 11ha and associate them with the biggest adjacent area. such polygons were eliminated because their shape/area did not make sense for practical terms, being interpreted as "noise".
The basic physiographic units were used to frame the information regarding soils and rock outcrops in order to delimitate the land units. Moreover, the predominant soil mapping units within each land unit were also considered, as well as the presence/absence of rock outcrops Therefore, information regarding soils was only framed on the basic physiographic unit polygons based on the most relevant soil types inside each polygon. information regarding rock outcrops was framed using the soil capability map and confirmed by field observations. after framing the information concerning the soil attribute and presence or absence of rock outcrops with the basic physiographic units, land units were obtained as a result. Thus, land units were described by geology, landforms, soils and presence/absence of rock outcrops. additionally, the resulting land units had to be generalized in order to compare them with the soil capability map. Thus, the texture and hardness of the lithological materials, topographical position and presence of rock outcrops were used as criteria for land unit assessment.
iii. resULts  The a (24.9ha) and D (38.1ha) lithological units are defined by its location and are associated with sandy beach and dunes, respectively. The e lithological unit corresponds to Lourinhã's oldest geological formation (Dagorda marls and dolomites), covering an area of 109.2ha.

Lithological units
Both the f and K lithological units are located in the Cesareda plateau and are distinguished from each other by the nature of their lithological materials; the former (860.6ha) includes calcareous formations in thick benches, while the latter (752.2ha) comprises marls and sandy limestones.
Unit G is the largest unit of the area (9,240.5ha) and includes materials corresponding to sandstones (variable grain size), marls, and some formations composed by clays and siltstones, and gravels.
Unit i (1,835.1ha) marks the transition between the G and K lithological units. This unit is distinct from the former G presenting an apparently homogeneous texture composed by quartzitic sandstones, and from unit K by the absence of calcareous materials.
in the southern part of the Lourinhã municipality area, the M lithological unit includes another calcareous formation (Vimeiro limestones), comprising an area of 69.7ha. finally, the n and J units, located in the typhonic valley of Bolhos (Manuppella et al., 1999b), distinguished themselves by the consolidation and maybe grain size of the lithological formations: unit n (206.1ha) corresponds to sandy deposits, whereas unit J (125.5ha) is associated with sandstones.

Landforms
The aggregated landform types resulting from the landform analysis model are described in table iii and represented in figure 4. The landform description was given according to the prevalent slope gradient class. The landform class flat to nearly flat sloping, occupying 5,020.6ha, includes alluvial and colluvial areas as well as its surrounding areas; also, it includes the Cesareda plateau, where several colluvial/alluvial areas occur, and the area corresponding to the typhonic valley of "Bolhos" with sandy materials.
The landform class gently sloping, comprises the largest area of the municipality (5,319.8ha) and is related to the top areas corresponding to the ridgelines, where most of the settlement occurs.
The strongly sloping landform class is primarily associated with the fluvial incision shaping the valleys; it is largely associated with the Rio Grande watershed and occupies about 3,890.9ha.
finally, the landform class steep slopes, has a smaller expression comprising 485.3ha of the municipality area. This class describes the areas along the coastline in transition to the cliffs and the cliffs themselves, and also areas sectioned by geological faults, as Cruz da Columbeira and Vila Verde faults (Manuppella et al., 1999b) that slice up the i lithological unit. fig. 4 -Distribution of the major landforms.

Soil -mapping units
The distribution of resulting eight soil -mapping units is shown in figure 5. Unit 1, mostly composed by Regosols and Cambisols (table iV) is the largest soil -mapping unit in the Lourinhã municipality, occupying 9, 394.5ha; small areas of Leptosols and eventually Luvisols and Vertisols may also occur. The cartographic information regarding these soil units suggests that the land use system has a strong influence on soil characteristics, as agropedic phases are frequent; in addition, mostly in the northern area of the municipality (Cesareda plateau), stony phases occur along with bedrock outcrops (fig. 9). The soils in this unit, generally present a coarse to medium or even fine texture, although local variations may occur since the parent material varies from fine sandstones, clays and mudstones, to coarser sandstones and related rocks (iHera, 1999); also, erosion and/or landslide deposits and soil creep may also influence the soils of this unit, as it can be realized by some families' designations resulting from colluvial deposition (iHera, 1999) and by references regarding landslide occurrences (Cipriano, 2001).
fig. 5 -Distribution of the general soil -mapping units.

Fig. 5 -Distribuição das unidades cartográficas gerais de solos.
Mapping of land units and land capability classification in Portugal: the case of the municipality of Lourinhã Unit 2 covers 2,758.8ha and may include Regosols, Calcisols and Leptosols mostly developed on compact calcareous rocks; the Regosols and Leptosols may be qualified as calcic (sensu WrB, 2006) in some areas. Unit 3 (129.2ha) including Calcisols and Regosols (mostly calcic) is largely developed over non -compact calcareous materials, and is mostly located in the northern part of the municipality.
Unit 4 includes Fluvisols and occurs in alluvial areas (as designated in both geological and soil maps), whereas unit 5 (associated with Regosols and Fluvisols) occurs in colluvial areas (as designated in the soil map), which are mostly designated as alluvial areas in the geological map; these soil -mapping units together comprise 1,454.7ha.

Land units
from the intersection between the landforms and lithological units, 33 basic physiographic units were delineated. as mentioned, after framing the soil attribute and the presence of rock outcrops to the basic physiographic units, 33 land units were defined. These resulting land units were then grouped in 11 general land units (table V, fig. 6). although most of these units are explained by their topographical position, the texture of the lithological materials and presence of rock outcrops was also used as criteria for their aggregation.
Land unit 1 represents the alluvial and colluvial areas mostly located at lower altitudes with flat slopes, corresponding mainly to Fluvisols (and eventually Regosols), and comprising an area of 1,473.8ha; field observations indicate that this area is not intensively used for agriculture, possibly due to flooding events and to hydromorphic soil conditions affecting plant rooting.
Unit 2 (27.4ha) is associated with the beach area, with sandy materials and flat slopes, whereas the land unit 3 is associated with Arenosols (over sand dunes) and occupies 36.9ha.
Unit 4 (2,528.1ha) occurs in the areas with gentle slopes adjacent to the main alluvial areas that are located on the slope base, and soils mostly correspond to Regosols (in part calcic) and Cambisols; drainage limitations may occur in this unit since some hydromorphic phases exist. fig. 6 -Distribution of the generalized land units. Unit 5 (3,360.4ha) represents the slope areas bordering the alluvial and colluvial areas and presents great grain size variability, corresponding mostly to Regosols, Cambisols and Leptosols, at a low extent; in this unit, erosion evidence was observed in some areas.
Land unit 6, where most of the settlement is concentrated, is the largest unit of the area (5,021.1ha) and describes the top sloping areas that include the ridgelines, Regosols and Cambisols being dominant, and, calcic Regosols and Leptosols occur at a lower extent; this unit, as mentioned for unit 5, presents great grain size variability associated with the parent materials corresponding to sandstones. Unit 7 is mainly located along the coast representing the cliffs and areas sectioned by geological faults and it describes exceptional situations of steep slopes, where Regosols, Cambisols and Leptosols may occur. it comprises an area of 387ha.
Unit 8 (1,160.6ha) is located in the Cesareda plateau and represents the areas with gentle slopes where limestone formations occur and the presence of rock outcrops (not clear by the information reported in the soil map; fig. 5 and table iV) is a main landscape feature ( fig. 7); soils mostly correspond to Regosols and Cambisols, but Calcisols and calcic Regosols may also occur in small areas.  The land units 9 and 10 are also located in the Cesareda plateau, likewise they have a significant presence of rock outcrops. Unit 9 (375.1ha), associated with Calcic Regosols, Regosols, Cambisols, Calcisols and Leptosols, describes the sloping areas, while unit 10 (29.2ha), associated with Regosols (partially calcic), Cambisols and Leptosols (partially calcic), describes areas with steep slopes.
finally, the land unit 11 (316.9ha) corresponds to a gentle sloping landform and is associated with a mosaic of Arenosols, Regosols and Podzols developed on the sandy deposits located in the "Bolhos" typhonic valley.

iV. DisCUssiOn
The applicability of the land unit concept involved a certain complexity given the type and quality of the available biophysical data for the area under study. firstly, difficulties arose from the scale and legend heterogeneity of the geological map, and strong uncertainty may exist regarding the top of the lithological formations as a result of the sedimentary nature of geological formations and also due to the fact that their description is primarily related with the age or deposition environments rather than lithological differentiation; that is, a strong heterogeneity of soil parent material is expected. secondly, the land unit delimitation process endorsed the recognition of strong deficiency of information on basic characteristics of soil units and therefore a limiting capacity to establish their classification; also, the rockiness (boulders and outcrops of bedrock material) was not revealed in the soil map and was confounded by the stoniness (coarse fragments). although the presence of bedrock outcrops (as a soil characteristic) is indicated in the soil capability map, it did not overcome the deficiency (and inadequacy) of information provided by the soil map (1:25 000) ( fig. 8, table i). The reorganization and improvement of information following the land unit assessment creates conditions for a deeper understanding of the land characteristics and qualities which may support the land capability and land suitability classification, enabling an adaptation of methodologies already used in other regions for the whole country (agroconsultores & Coba, 1991;agroconsultores & Geometral, 1995;agroconsultires & Geometral, 2004). That is, the land unit approach may be a useful tool to support land planning decisions (at different scales) and legislation (ran) using identical criteria.
Our study indicates some similarity between the land units and the distribution of some soil capability classes and subclasses (figs. 6 and 8). for example, land unit 1 largely coincides with areas of the a soil capability class and some areas of the Bh and Bs subclasses (figs 6 and 8); also, Ds and es subclasses largely coincides with the land units 8 and 9. However, our study also suggests that the landscape discrimination by the land unit approach is more detailed than that corresponding to the distribution of soil capability classes and subclasses, as shown in the soil capability map ( fig. 8). for example, the Ce soil capability subclass comprising 6,751.7ha (about 45.9% of the total area; table Vi) corresponds to the largest area of the Lourinhã Municipality, and includes mostly areas of land units 4 and 6 and areas of land unit 5 (figs. 6 and 8). This suggests that the land unit delimitation, following an improved framework of permanent biophysical features of the land, may facilitate a deeper land capability assessment (soil capability classes). similar considerations could be developed for the De soil capability subclass.
according to the soil Capability Handbook of Portugal (srOa, 1972), the subclass Ce corresponds to areas with severe limitations regarding erosion risks, with slopes up to 25% and limited effective soil depth (> 25cm), and may be used for occasional cropping. nevertheless, the areas of the subclass Ce (soil capability map; fig. 8) mostly correspond to those of land units 4, 6 and 5, where the dominant slope gradient classes are 0 -5%, 5 -8% and 8 -12%, respectively ( fig. 6 and table iii). although some soils in such subclass areas may show an effective depth lower than 50cm, most of them have an effective depth greater than 50cm (srOa, 1973), and therefore could be considered in the Be capability subclass (slope up to 15% and effective soil depth > 35cm); also, most areas of those land units could be included in the a2 land suitability established by agroconsultores & Geometral (1995) for Entre Douro e Minho for areas where slopes may range up to 15% and the effective soil depth is at least 50cm. indeed, in field observations it was observed that most of the Ce subclass areas are intensely used for agriculture, as opposed to what is considered for the Ce soil capability subclass (srOa, 1972). as long as agricultural management practices are implemented to control the erosion processes, such areas may not show biophysical constraints to be classified as Ce. similar considerations could be outlined for most of the areas classified in De soil capability subclass (figs. 6 and 8). according to the soil capability map, only 14.4% of the area (6% of class a and 8.4% of class B) is well suited for intensive cropping and grazing activities; however, most of this land is not used for intensive cropping, which may be associated with flooding susceptibility and drainage limitations in some areas, as suggested by the hydromorphic and flooding phases reported in the soil map (table i). Therefore, considering that the subclass Ch only occupy 0.2% of the municipality area (table Vi), only 14.6% of the area would be, in principle, considered within the national agricultural reserve (according to the Ministerial Order nº 73/2009, 31 st March). as mentioned, large areas of the Lourinhã municipality, especially those corresponding to land units 4, 6 and 5, could be positioned in the B capability class (mostly Be subclass), or in a2 suitability class (moderate suitability for agricultural crops) established for The Entre -Douro e Minho region (agroconsultores & Geometral, 19995). Consequently, following the same Ministerial Order nº 73/2009, 31 st March, that recognizes that "are part of the ran the land units that present a high or moderate suitability for agriculture", it is concluded that ran areas could reach a high proportion of the municipality area, as the proportion of those land units reach about 74% (table V). This large discrepancy is clearly associated with the fact that the soil capability units were not delimitated within a framework of land units, with most of the area considered to have steep slopes (Ce subclass). in short, the land unit delimitation may allow an improved framework of the soil capability classes and subclasses as well as an improved outline to identify and delimitate ran areas.
Land physical limitations in the Lourinhã municipality are mostly associated with erosion, rockiness, drainage/flooding and climate. Meanwhile, the area may be largely prone to erosion susceptibility, as about 72.5% (table Vi) of the area is included in erosion subclasses of different classes. although the apparent contradiction with aforementioned considerations regarding land capability classification, it should be emphasized that ero-sion (sheet and rill erosion) is common and pronounced even in areas with gentle slopes (as land unit 6 with a prevalent slope gradient between 5 -8%), as recognized by field observations ( fig. 9).
fig. 9 -soil erosion observed during the field observations. Colour figure available online. Besides the soil erodibility and the slope angle, soil erosion is also associated with the removal of hedgerows, stonewalls and other fences (delimitating land parcels) for easier crop management, resulting in an increment of slope length and uniformity, and with reduction of soil protection ( fig. 10). That is, despite the slope and soil erodibility effects, erosion susceptibilities in the municipality of Lourinhã are strongly related to the great changes in the land use system and the absence of required conservation measures. also, landslides and soil creep should be also taken into account for land capability purposes, as Cipriano (2001, p. 111) reported that "the nature of the clay soils may be the cause of the population's decay, as prolonged rains create landslides that drag everything with it".
fig. 10 -Coastal areas with hedgerows (a) and innermost areas with large crop fields (B).
Colour figure available online. Physical limitations related to rockiness are mostly associated with the limestone formations (land units 8 and 9) in the Cesareda plateau (ribeiro, 1940; fig. 7), and may have a strong impact on trafficability. Likewise, limitations regarding flooding and drainage must be considered, especially in those areas that apparently did not present constraints for agricultural purposes (as the soil capability map states). Despite drought associated with uneven rainfall distribution, limitations regarding climate may be also strongly associated with the effect of wind in exposed areas, as is suggested by the landscape compartment in areas close to the coast ( fig. 10), which in turn may contribute to minimize soil erosion.
V. COnCLUsiOns from the present study, it becomes evident that land unit mapping can potentially provide greater information than the soil capability map. The land unit mapping is an essential tool to assist in landscape planning, by assisting data collection and by improving the information available in soil and soil capability maps, as it may aid in the reorganization, requalification and homogenization of available information. even in areas with great complexity, as the area under study, land unit mapping enables the outlining of unified methodologies, with identical criteria, to assess the capability and suitability of the land. The land unit mapping may enable the development of a unified system to more accurately delineate management tools (such as the ran) for the whole country, meeting the general requirements of the existing legislation. Thereby it is an essential support for the landscape planning decisions. as land unit mapping is a top -down process that requires a careful data analysis, further research should be developed in larger areas to clarify its potential benefits regarding land planning in Portugal. The land unit delimitation is an essential tool to improve land limitations assessment and support land capability and suitability classifications. aCKnOWLeDGeMents authors gratefully acknowledge isabel roussado for helping handling data and Paulo Marques for fieldwork assistance. The support of Prof. ana Carla Madeira for english style improvement is also acknowledged.