Next Article in Journal
Evaluation of Freshwater Provisioning for Different Ecosystem Services in the Upper Mississippi River Basin: Current Status and Drivers
Next Article in Special Issue
Water Level Loggers as a Low-Cost Tool for Monitoring of Stormwater Control Measures
Previous Article in Journal
Isotopic Characterization of River Waters and Water Source Identification in an Inland River, Central Asia
Previous Article in Special Issue
Impacts of Salinity on Saint-Augustin Lake, Canada: Remediation Measures at Watershed Scale
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Probability of Surface Artificialization in Zêzere Watershed (Portugal) Using Environmental Data

1
Centre for Geographical Studies, Institute of Geography and Spatial Planning, Universidade de Lisboa, Edif. IGOT, Rua Branca Edmée Marques, 1600-276 Lisboa, Portugal
2
General Directorate for Territorial Development (DGT), Rua da Artilharia Um, 107, 1099-052 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Water 2016, 8(7), 289; https://doi.org/10.3390/w8070289
Submission received: 5 April 2016 / Revised: 7 June 2016 / Accepted: 7 July 2016 / Published: 13 July 2016
(This article belongs to the Special Issue Watershed Protection and Management)

Abstract

:
The land use and land cover (LUC) of the Zêzere watershed (Portugal) have undergone major changes in recent decades, with the increase of artificial surfaces. This trend is quantified in some studies, but the probability of the increase of this type of LUC, nor the places where the next transitions or land use/cover changes (LUCC) for artificial surfaces will have high probability of occurrence has not yet been assessed. This research presents an evaluation of these two aspects, by means of bivariate statistical models (fuzzy logic and information value) and environmental data. The artificialization probability by sectors within the same watershed is also evaluated, to further understand which areas will require greater attention, taking into account the environmental conditions favorable to the occurrence of this process and bearing in mind the conditions under which this process took place in the past. The results obtained using these models were assessed independently, through curves of success, noting that the modeling through the fuzzy gamma presents slightly better efficiency in determining the probability of artificialization surfaces in the study area. The area with the highest probability of artificialization is mostly located in the SW of this watershed, but high probabilities are also present in the upstream sector, being those areas that require further preventive measures once they have influence on the water quality and quantity in the main reservoirs of this watershed.

Graphical Abstract

1. The Artificialization of Surfaces and Their Assessment

Urban growth has been evaluated in different territories around the world, considering the spatiotemporal aspect, as well as the factors that induced it [1,2,3,4,5,6]. Other studies have emerged in the context of land use and land cover changes (LUCC) to assess where the soils are being converted into artificial surfaces [7,8,9,10,11], especially when there is loss of soils essential for the development of agricultural practices, or the conversion of forest areas. These processes have multiple socio-economic and environmental impacts [12,13,14,15,16]. In this sense, the determination of driving forces (socio-economic, environmental, or other) that are the cause of LUCC is fundamental to understand the factors that induced them and for the creation of measures aimed at the sustainability of land use [15,17]. Some studies demonstrating the links between driving forces and LUCC have appeared recently, also quantifying how much each factor represents to the observed changes [13,18,19,20,21,22].
In Portugal, studies addressed the impacts of LUCC in the environment, in particular the major transitions of land use and land cover (LUC) resulting from deforestation and their contribution to emissions and removals of CO2 [23], or in water bodies [17,24,25].
The evaluation of disturbances caused by LUCC in water availability (quantity and quality) are also essential for a sound land use planning [17], particularly in places where there was already degradation of this natural resource by the occurrence of certain harmful events, e.g., forest fires [26].
Urban growth can be considered a negative factor on water availability, particularly where this causes water stress [17]. In this context it is desirable to minimize the artificialization in the neighborhood of important water reservoirs used to supply the populations and their activities, namely within the influence area of water catchment towers.
The artificialization of surfaces has been quantified in the Portuguese territory based on cartography produced in recent years (LUC maps of Portugal, Landyn maps, and CORINE Land Cover—CLC). Data consistency increased a great deal with the most recently reviewed and published datasets [8,27,28].
However, so far there are no estimates of the possible increase in this type of LUC (in area), or the soils which will most likely be occupied or suffer LUCC for this type of LUC. In hydrographic basins the increase of this type of LUC has negative impacts, especially in the increase of surface runoff and also in the increase of physical and chemical substances of anthropic origin that are drained to the water bodies, causing their degradation [17].
The evaluation of urban growth has been based on methodologies using Geographic Information Systems (GIS) and remote sensing tools for the assessment of LUCC [1,3,5,29]. These tools for the collection, processing, and analysis of geographic information (GI) allow more detailed analysis of LUCC and the monitoring of natural or anthropogenic processes [30] that occur in the territories (e.g., artificialization surfaces).
For the evaluation of the LUCC, different methodologies have been used, especially the cellular automata model and artificial neural networks, with good results in the estimation and representation of land cover dynamics [6,31,32,33,34]. In the estimation of future LUCC, the probabilistic cellular automata-Markov model has been used in several studies [6,35,36,37,38].
Many studies have determined the probability of occurrence of a phenomenon in the future (e.g., landslides, risk of forest fire, floods, etc.) [39,40,41,42] based on the conditions observed in the past that caused certain phenomena in a certain place. These methodologies have a strong statistical component, differing only in the integration of the respective variables and method of calculation (e.g., bivariate methods and logistic regression). These factors can provide different results and, consequently, different interpretations. Thus, the methods used for the validation of results are essential to understand what is the best method for modeling a given phenomenon [42].
Some methodologies for validating results have emerged, including the modeling of the probability of occurrence of a natural phenomenon with verified occurrences and, thus, allowing the verification of the final results (through the intersection of GI) if the same occurrences fall on the areas with the highest probability of occurrence. Another option is to use part of the dependent variable for modeling and another part to validate the results obtained (random partition). These methods allow the development of success or prediction curves [43] and receiver operating characteristic (ROC) [44,45], and to quantify their robustness for modeling.

2. Main Objectives

This research aims to test two bivariate statistical models for the identification of areas with higher probability of surface artificialization in the Zêzere watershed (Portugal), bearing in mind the increase verified in the last two decades.
The information used is essentially environment-related data, once we considered in this study that socio-economic conditions remain constant, so this information was not part of the model.
The validation of results will take place with information of LUC, namely information about artificial surfaces for two different times: 1990 and 2012. The information of the first year will be used to determine the areas with highest probability of artificialization and the information of the last year will serve to validate the results. We also considered an intermediate date (2000) for verification of the probability of surface artificialization, obtained with data from 1990, being that this data is complementary for the assessment of the validation performed with 2012 data.
Another goal of this study is to determine the differences between the two models in different locations (sectors) within the watershed, in order to quantify and understand the differences between all sectors. Special attention is taken in assessing which are the most important independent variables for determining the probability of artificialization for this territory.

3. Research Area

The study area is the Zêzere watershed (Portugal) (Figure 1). This watershed has an area of 502,278.4 ha, and includes one of the main drinking water reservoirs (Castelo de Bode) in Continental Portugal. The surface artificialization in 1990 was approximately 1.2% of the total area of the watershed, but in the last two decades it has increased approximately 0.5% (according to data from CLC of 1990 and 2012), especially in the vicinity of water bodies. This factor can induce water stress within this watershed, namely decreasing water quality and increasing water treatment costs [17,46].
The watershed was divided into three sectors (Figure 1): two of them based on the location of the main reservoirs (A and B) and one based on the sub-basins without reservoirs (C). The objective of this division areas is to assess the differences in the probability of artificialization: initially within the watershed, considered as a whole, with the selected statistical models, and later to evaluate and differentiate, spatially, this probability between the different watershed sectors and their behavior when compared to the whole watershed results.
In this watershed, since 1990, the artificial surfaces areas increased, particularly in the sub-class discontinuous urban fabric (Table 1). This increase in dispersed artificial surfaces carries greater economic and financial burden, associated in particular with the construction of basic infrastructures (roads, sewage systems, drinking water, power supply, etc.). Here there are also records of environmental disturbances, particularly at the level of domestic waste and water treatment, confirming a deficiency in this area.
In the watershed upstream sector (A) it is located the Estrela mountain (Serra da Estrela). Here the hillsides with slopes exceeding 20% are one of the constraints to urbanization and, thus, the location of major cities and towns is primarily in locations with lower altitude and reduced slope. The Cabril dam is located in this sector of the watershed; an important infrastructure for public water supply.
Sector B comprises the Castelo de Bode reservoir, this area being mainly occupied by scrub and/or herbaceous vegetation associations, resulting mainly from large forest transitions that occurred between 1990 and 2012. This sector also comprises hillsides with slopes exceeding 20%.
The geomorphology of the downstream sector (Sector C) is totally different from the upstream sectors, i.e., the relief is more flat, and was one of the factors that facilitated the urban settlement (e.g., Tomar City). This sector is also characterized by other factors, such as soil characteristics, favorable to the development of agricultural practices, high sun exposure, higher temperature, proximity to water courses, and proximity to railways connecting the area with larger urban areas, like Lisbon, among others. Nevertheless, this is the most relevant sector within the Zêzere watershed once it includes the most relevant water catchment within the Portuguese drinking water supply infrastructure.
The artificial surfaces increase at this watershed occurred in areas very close to the major urban centers within the watershed (Tomar, Fátima, and Covilhã). The analysis of the geographical dispersion of the artificial surfaces (from data of the CLC) found that the new areas have emerged as the expansion of existing ones, particularly on the periphery of the main larger conurbations (Figure 2). It was also found that the largest increase in artificial surfaces occurred until the year 2000, within less than 5 km of the urban centers previously referred.

4. Data, Tools, and Methods

4.1. Statistical Methods and Validation

The areas with higher probability of artificialization were determined using fuzzy logic [47] and information value [48] methods. This last methodology was adapted to the study of surface artificialization probability. In these methods it is assumed that a particular phenomenon that occurred in a given territory has a probability to happen in the future under the same conditions under which it occurred in the past [42,49,50,51].
For the determination of information value, first we calculated a priori and conditional probabilities. This methodology was presented in Meneses and Zêzere et al. [42,52]. After calculating these probabilities, we used the following equation to determine the information value [48]:
I i = L n C P j i P p
In Equation (1), Ii is the information value; CPji is the conditional probability of surface artificialization on class i of the thematic map j; and Pp is the a priori probability of occurrence of surface artificialization.
In assessing the final probability of surfaces artificialization, i.e., integrate information values of all independent variables, we considered the following equation:
I j = i = 0 n X i j I i
In Equation (2), Ij is the total information value of pixel j, Ii is the information value of each pixel of each independent variable, n is the number of variables, Xij assumes the value 1 or 0 depending on the presence or not of the variable in the field unit.
The fuzzy logic methodology, developed by Zadeh [53], admits the variation between 0 and 1 (or 0 and 100%) of an existing element in a given set, this being expressed by a fuzzy membership value. According to Bonham-Carter [47], the assignment of values of fuzzy membership to every variable is typically made from the subjective evaluation (expert opinion) of their importance in the model, so it is considered an heuristic method. However, in this study, the fuzzy membership values were assigned objectively to each class of independent variables, i.e., its importance was calculated in proportion to each calculated conditional probability. The maximum value of all of the independent variables was determined and, from this, for each class, the respective fuzzy membership value was calculated as the result of dividing the respective conditional probability by the maximum value found earlier.
Among the various operators for combining fuzzy membership values, we used the fuzzy gamma operator, because this combines two operators [47]: sum and algebraic product. The Fuzzy gamma operator uses the following equation [47]:
F u z z y   G a m m a = ( 1 i = 1 n ( 1 μ i ) ) y ( i = 1 n μ i ) 1 y
In Equation (3), µi is the fuzzy association values (i = 1, 2, 3, …, n) for the variables 1, 2, 3, …, n; n corresponds to the number of variables considered, and y the parameter set by the operator.
The geographic information of the artificial surfaces considered in this study is the result of the CLC for the years 1990 and 2012 (level 2), comprising approximately 5978 and 8588 ha, respectively. The GI of the first year (dependent variable) was used for modeling with the methods presented earlier; and the GI for the last year was used for validation of the results, considering in this procedure only the artificial surfaces that have emerged between the two years (approximately 2610 ha). By applying this procedure, we want to know if the artificial surfaces of 2012 coincide with the areas with the highest probability of artificialization obtained by the models previously described.
The process used for the validation of results included the elaboration of success curves and the measurement of the area under the curve (AUC), according to the methodology described by Meneses [42] and Tehrany et al. [43]. This method of validation enables the evaluation of the robustness of the models presented for the determination of artificial surfaces probabilities.
In the analysis of the results the sectors A, B, and C (Figure 1) were distinguished in order to determine the possible differences between the results of the two models on the probability of artificialization surface within the same watershed.
The importance of each independent variable in the process of artificialization for each sector was also determined, in order to understand what the spatial influence of each predisposition factor in the development of this process is. In this procedure the accountability (AI) and reliability (RI) indexes [42,54] were determined using Equations (4) and (5). AI explains how various classes of predisposition factors are relevant in the analysis because they contain artificial surfaces, while RI depends on the average density of artificial surfaces in classes of predisposition factors most relevant to the development of this process.
A I = i = 1 n k N 100
R I = i = 1 n k i = 1 n y 100
In Equations (4) and (5), A I is the accountability index; R I is the reliability index; k is the area of artificial surfaces in classes with values of conditioned probabilities superior to a priori probabilities; N is the total area of artificial surfaces; y the area of each class of independent variable with conditioned probability above the a priori probability.

4.2. Data Collection and Tools

The GI themes considered in this research are diversified (Figure 3). Cartography of the Portuguese Environment Agency (soil maps, insolation, humidity, temperature, and precipitation), available online, was used. The LULC considered is the CLC data produced by the Portuguese General Directorate for Territorial Development. Due to the similarity between the spatial geological units and the soil types, we opted to use only this last variable.
Variables derived from the Digital Elevation Model (DEM) from the GMES RDA project (EU‑DEM) made available by the European Environment Agency, were also used, in particular the slope, aspect, and hillshade. These variables derived from the DEM also calculated the topographic wetness index—TWI [55].
In addition to this GI, the distance to water bodies and watercourses (DWBW) was also obtained, to assess the influence of their location on the surface artificialization that occurred over the last two decades, since the vicinity of the main reservoirs in the Zêzere watershed [46] was subject to an increase in housing developments.
All collected GI was harmonized and manipulated in the Geographic Information Systems (GIS), using the software ArcGIS 10.2 and ILWIS 3.4. The GI of the variables in vector structure was converted to raster (pixel 10 × 10 m) for carrying out the various processing procedures. The selection of this resolution was the result of several geoinformation conversion tests, where different modeling resolutions were tested in this research, and it was found that the results of the adopted resolution are similar to those obtained using data with higher resolution. After several vector-to-raster conversions, the resolution considered is the one that allows better results in the determination of artificialization probabilities, since the surface loses a lot of information in the process of generalization that occurs during data conversion, in particular in the conversion of LUC, where small buildings were not considered in the modeling due to the pixel size (e.g., 20 m × 20 m = 400 m2), and some of the buildings that appeared during the period under evaluation are also not included in the validation process.
For the intersection of the variables presented in Figure 3 we used only ArcGIS 10.2. This software was also used for determining the areas with greater probability of artificialization. On the fuzzification of fuzzy membership values we used the Spatial Data Modeler (ArcSDM), module added to ArcGIS 10.2.
In order to verify that the artificialization surface probabilities obtained from the year 1990 (CLC data) have sequences 10 years later, we also calculated the probabilities for the year 2000. With these results we intended to check the correspondence between the two years (similar odds for the same class of independent variables) so as to enhance the use of the information of 2012 in the context of the models adopted.

5. Results and Analyses

5.1. Conditioned Probabilities of Artificialization Surfaces in the Zêzere Watershed

Each class of independent variables considered in this study presents a different probability of surface artificialization (Table 2). The flat surfaces are those that have a higher probability (higher conditional probabilities) in the variable Aspect, a fact confirmed also in the variable Slope (greater in weak slopes); for hillshade, with records varying between 0 and 300, the most influential is the class 150–200; also with high probability are the surfaces with higher humidity, insolation, and temperature. The distance to water bodies and watercourses is also one factor that conditions the surface artificialization, where the 4–5 km buffer in DWBW has the highest probability (CPji AS 1990 in Table 2), due to the total artificial surface included in it. However, the artificial surface is high along water courses (up to 1 km vicinity), revealing a pattern in the preference for these surfaces for the location of infrastructures (housing, industrial complexes, etc.). The Topographic Wetness Index (TWI) reveals that the surfaces with lower values are the ones that have less probability of artificialization.
For calculating the information value, the a priori probability (probability to find artificial surfaces) should be also calculated, which is 0.012. This value results from the division of the total area of artificial surfaces by the total area of each independent variable.
The relationship (coefficient of determination—R2) between the conditional probabilities obtained from information of artificial surfaces of 1990/2000 is high for most independent variables (Figure 4), with the exception of DWBW, where it was observed a significant increase of artificial surfaces in more distant areas to water bodies and watercourses in the year 2000, given the area occupied by artificial surfaces in 1990. The R2 between the conditional probabilities of 1990/2012 is also high, with some variables with lower R2 compared to 1990/2000, such as the Insulation and TWI, but the DWBW highlights with a higher R2, indicating the greater probability of artificialization in the vicinity of water bodies.
Considering the conditional probabilities of 1990/2000, 2000/2012, and 1990/2012 of all independent variables it was verified that a high correspondence between the results was obtained (R2 = 0.994, 0.999 and 0.992, respectively).

5.2. Spatial Variation of the Probability of Artificialization Surfaces in the Zêzere Watershed

The artificialization surface probability is higher in the downstream sector of the watershed, as shown by the maps resulting from application of statistical methods described above (Figure 5).
The map with the results of information value shows the sectors upstream and downstream of the watershed, especially this last one, with a high probability of artificialization surface.
The map with the results of the fuzzy gamma operator shows the highest probability in the vicinity of the main water lines in the southwest of the Nabão River, but also shows that the central areas in the upstream sector have some probability of artificialization surface.
The success curve characterizes the quality of a forecast system by describing the system’s ability to correctly anticipate the development or non-development of a predefined event [56,57]. In this case, the predefined event is the artificialization of the surfaces in the Zêzere watershed, and the probability of artificialization obtained was validated with the artificial surfaces that have emerged up to 2012. During the validation of results it was found that both models have similar robustness for the determination of artificial surfaces probabilities, as noted through with success rate curves (Figure 6). With around 35% of total area of the watershed classified in descending order of probability of artificialization, about 80% of the artificial surfaces of 2012 are validated, and with about 50% validated, 90% of the total artificialized area (Figure 6). However, modeling using a fuzzy gamma operator presented slightly better results, with AUC of 79.4%, compared to the AUC of 79.1% obtained by the information value.

5.3. Artificialization Surfaces by the Sectors of the Zêzere Watershed and the Relevance of the Environmental Predisposing Factors to Artificialization Process

The surface artificialization probability is highly variable, spatially, depending on the conditions existing in the territory. The statistical results presented in Table 3 demonstrate high variability of artificialization probability between the sectors defined in the watershed, showing areas A and B with the minimum values, and sector C with maximum probability of artificialization.
In order to assess the artificialization probability it is also important to know which environmental predisposition factors are the most important. Thus, the accountability and reliability index for each sector of the Zêzere watershed were determined, where the three most relevant predisposition factors in the analysis are highlighted because they contain artificial surfaces in the year 1990 (bold values in Table 4).
The slope turned out to be a very important factor for the artificialization in the three sectors, mainly because of the high area of artificial surfaces in the smaller slope classes. Hillshade is also an important predisposing factor, particularly in sectors B and C. Once sector A includes areas of higher altitude and with a more rugged relief, the humidity and the insolation are the most important factors for the artificialization in this sector; however, the average density of artificial surfaces in each precipitation class is also important considering the set of predisposition factors.
In sector B the temperature has an important influence on surface artificialization, but an analysis of the medium density of artificial surfaces for each independent variable in this sector shows the importance of the type of soils, in particular the Calcic luvisols, where the main drinking water reserves are located.
The C sector soil is also important, considering the high artificialization that occurred in certain classes, but the insolation is an important factor in the analysis of the average density of artificial surfaces by class.

6. Discussion

The probability of surface artificialization obtained by the information value and fuzzy gamma methods have many similarities in the study area.
These results indicate a good performance of these methods in modeling surface artificialization probabilities, i.e., the artificial surfaces used for modeling (1990) made it possible to differentiate the territory with different probability of artificialization, confirming that a part of these surfaces with a high probability the artificialization occurred (a fact confirmed by overlapping them with the artificial surfaces of 2012). However, it was also found that the new artificial surfaces are located mostly in the periphery of the existing ones, in particular of larger artificial surfaces, which justifies the good performance of the models, i.e., if the expansion occurred mainly from the larger urban centers, it is also more likely that the expansion areas exhibit the same characteristics as the areas of the original urban centers.
The robustness of the methods used has been tested in different studies for determining the susceptibility or probability of occurrence of natural phenomena (e.g., landslides, forest fires, etc.), where the best results have been achieved by the information value model [42]; however, in this research the fuzzy gamma model presented the best performance (Figure 6), although by a marginal difference in relation to the results of the information value model.
The two maps (Figure 5) show that the central sector (B) of the watershed is where there is less probability of artificialization (soils occupied mostly by forest, scrub, and herbaceous vegetation). The natural conditions of this sector are less favorable to human occupation, in particular due to the slope (>25%). This conjugation of less favorable factors to surface artificialization has been referred in some other studies, e.g., Druga and Faltan [58].
The spatial distribution of artificial surfaces in the Zêzere watershed is very uneven throughout the territory and its location is conditioned mostly by the same morphology, this factor being the most important on the distribution of artificial surfaces [58].
The change in the area of artificial surfaces in the Zêzere watershed was evaluated at different times (1990, 2000, 2006, and 2012), which identified the consequences of an increase in the quality of surface water [17,59]. Taking into account the importance of this natural resource and the interference of artificial surfaces in its quality, it is essential to know which areas are those with the greatest artificialization probability in order to avoid new construction, especially in the vicinity of water bodies. These areas are currently experiencing an increase in demand due to the scenic context and watersports, among others [17].
However, in this research it was not proved that the distance to water bodies and watercourses is an important variable in determining the probability of artificial surfaces, because urbanization is still in the process of development in the vicinity of water bodies and infrastructure currently located in these areas is scattered. Additionally, the spatial resolution that characterizes CLC cartography (25 ha) induces limitations to this analysis that require further work. This fact demonstrates the importance of knowledge about the properties of the available GI datasets, and their influence in the results presented in this research that have to be taken into account.
The urban development is also influenced by the location of the main housing clusters [5], which also influence the surface artificialization process, particularly on the periphery of these clusters, where there have been some cases of new housing and infrastructure construction (e.g., roads, railways, industrial complexes, and support equipment, etc.) [3,6]. As a result of this research we highlight the areas with the greatest probability of artificialization (e.g., the SW of sector C), taking into account that under the same environmental conditions which created artificial surfaces in the past, these areas might be artificialized in the future, based on the concept of uniformitarianism implicit in the methodologies used [42,60,61]. Since human intervention was not included in the modeling procedures, and is one of the main agents of LUCC [17,59], some discussion here is justifiable. On the one hand, the application of these variables is impossible, as there are not enough data to demonstrate certain conditions in the past to the watershed under study (e.g., socio-economic power and conditions of the families, infrastructure, or urbanization index, search for housing, or construction of certain infrastructures). On the other hand, there is a very high uncertainty about these human conditions for large periods (next decades), so this approach can only be carried out assuming different kinds of scenarios for the future.
The influence of these man-made factors in the artificialization surface, namely the urbanization in the proximity of road networks, can form the basis of the explanation of the results that are not explained by the models used in this research. However, the option to use only environmental data is due to the fact that the variables used do not present large variations in relatively short periods, such as those considered in this investigation. However, it is acknowledged that this information resulting from anthropic actions is important for determining the areas that will, in the future, become artificial surfaces, but this factor depends only on human actions, such as building new roads or other infrastructure essential to the location of people and goods. Thus, these anthropic actions are encompassed in the process of artificialization and the resulting data are considered artificial surfaces in CLC data.

7. Conclusions

Surfaces artificialization in the Zêzere watershed is more likely to occur in the downstream sector (sub-basin of Nabão River), and this is the place with the highest density of artificial surfaces at present.
The determination of the areas with highest probability of artificialization in the Zêzere watershed, using only environmental data, showed good results, a fact confirmed in the validation (through the curves of success) of the results obtained by the methods of information value and fuzzy gamma.
By comparing the artificialization probability of the sectors delimited in this watershed, it was observed that there is similarity between the results of the two methods obtained for the whole area of the watershed. However, the results differ among sectors, with the highest probability of artificialization in sector C. This spatial differentiation is essential for decision-making in land use planning; in particular, for determining the possible interferences resulting from artificialization in the vicinity of important water bodies to the public water supply. Once there are favorable environmental conditions for this process to occur in these areas, these surfaces present some probability of artificialization. Yet, in this context, it was noted that there are conditions for the increase of the artificialization in the upstream area of the watershed, but the development of this process at this location can be negative in the maintenance of downstream water bodies (water stress). This process leads to important impacts on water quality, mainly due to urban growth in a disorderly manner, and characterized by deficient sewage network systems in the vicinity of the reservoirs.
The spatial differentiation of the artificialization probability will allow better decisions (preventive or reactive) in the territory. In our case it is a watershed with major water bodies in Continental Portugal, and better decisions will contribute to minimize possible consequences resulting from transitions of other types of LUC to artificial surfaces. However, this assessment must include other factors, such as socio-economic conditions [62], restrictions, or obligations set out in territorial and sectorial management plans and programs, along with the chosen strategy to monitor legal framework implementation.

Acknowledgments

This work was supported by a grant from Institute of Geography and Spatial Planning and Universidade de Lisboa (BD2015), and eENVplus Project (CIP-ICT-PSP–grant No. 325232). The work of Eusébio Reis is supported by the FCT project FORLAND (PTDC/ATPGEO/1660/2014).

Author Contributions

Bruno M. Meneses worked on the methodology for LULC mapping, performed the experiment work (probability of surfaces artificialization) and statistical analysis, produced the tables and figures, and wrote the paper. Eusébio Reis, Maria J. Vale and Rui Reis read and made improvements to the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bagan, H.; Yamagata, Y. Remote Sensing of Environment Landsat analysis of urban growth: How Tokyo became the world’s largest megacity during the last 40 years. Remote Sens. Environ. 2012, 127, 210–222. [Google Scholar] [CrossRef]
  2. Berdegué, J.A.; Carriazo, F.; Jara, B.; Modrego, F.; Soloaga, I. Cities, Territories, and Inclusive Growth: Unraveling Urban–Rural Linkages in Chile, Colombia, and Mexico. World Dev. 2015, 73, 56–71. [Google Scholar] [CrossRef]
  3. Hegazy, I.R.; Kaloop, M.R. Monitoring urban growth and land use change detection with GIS and remote sensing techniques in Daqahlia governorate Egypt. Int. J. Sustain. Built Environ. 2015, 4, 117–124. [Google Scholar] [CrossRef]
  4. Kourtit, K.; Nijkamp, P.; Reid, N. The new urban world: Challenges and policy. Appl. Geogr. 2014, 49, 1–3. [Google Scholar] [CrossRef]
  5. Shafizadeh-Moghadam, H.; Helbich, M. Spatiotemporal variability of urban growth factors : A global and local perspective on the megacity of Mumbai. Int. J. Appl. Earth Obs. Geoinform. 2015, 35, 187–198. [Google Scholar] [CrossRef]
  6. Shafizadeh-Moghadam, H.; Helbich, M. Spatiotemporal urbanization processes in the megacity of Mumbai, India: A Markov chains-cellular automata urban growth model. Appl. Geogr. 2013, 40, 140–149. [Google Scholar] [CrossRef]
  7. Direção Geral do Território (DGT). Land Use and Land Cover Evolutions in Continental Portugal. Work to Support Reporting of Emissions and Carbon Sequestration in the Sector Use and Land Use Changes. Kyoto Protocol and United Nations Framework Convention on Climate Changes; DGT: Lisbon, Portugal, 2014. [Google Scholar]
  8. Meneses, B.M.; Vale, M.J.; Reis, R.; Marrecas, P.; Barreira, E. Metodologias para a avaliação do uso e ocupação do solo em diferentes épocas. In Uso e Ocupação do Solo em Portugal Continental: Avaliação e Cenários Futuros. Projeto LANDYN; Direção-Geral do Território: Lisboa, Portugal, 2014; pp. 5–16. [Google Scholar]
  9. Quan, B.; Chen, J.-F.; Qiu, H.-L.; Romkens, M.J.M.; Yang, X.-Q.; Jiang, S.-F.; Li, B.-C. Spatial-temporal pattern and driving forces of land use changes in Xiamen. Pedosphere 2006, 16, 477–488. [Google Scholar] [CrossRef]
  10. Rodriguez-Galiano, V.; Chica-Olmo, M. Land cover change analysis of a Mediterranean area in Spain using different sources of data: Multi-seasonal Landsat images, land surface temperature, digital terrain models and texture. Appl. Geogr. 2012, 35, 208–218. [Google Scholar] [CrossRef]
  11. Teixeira, Z.; Teixeira, H.; Marques, J.C. Systematic processes of land use/land cover change to identify relevant driving forces: Implications on water quality. Sci. Total Environ. 2014, 471, 1320–1335. [Google Scholar] [CrossRef] [PubMed]
  12. Fiquepron, J.; Garcia, S.; Stenger, A. Land use impact on water quality: Valuing forest services in terms of the water supply sector. J. Environ. Manag. 2013, 126, 113–121. [Google Scholar] [CrossRef] [PubMed]
  13. Long, H.; Tang, G.; Li, X.; Heilig, G.K. Socio-economic driving forces of land-use change in Kunshan, the Yangtze River Delta economic area of China. J. Environ. Manag. 2007, 83, 351–364. [Google Scholar] [CrossRef] [PubMed]
  14. Nunes, A.N.; de Almeida, A.C.; Coelho, C.O.A. Impacts of land use and cover type on runoff and soil erosion in a marginal area of Portugal. Appl. Geogr. 2011, 31, 687–699. [Google Scholar] [CrossRef]
  15. Ribeiro, M.; Vale, M.J.; Reis, R. Identificação das principais forças motrizes: Abordagem quantitativa. In Uso e Ocupação do Solo em Portugal Continental: Avaliação e Cenários Futuros. Projeto LANDYN; Direção-Geral do Território: Lisboa, Portugal, 2014; pp. 81–88. [Google Scholar]
  16. Schoene, D.H.F.; Bernier, P.Y. Adapting forestry and forests to climate change: A challenge to change the paradigm. For. Policy Econ. 2012, 24, 12–19. [Google Scholar] [CrossRef]
  17. Vale, M.J.; Meneses, B.M.; Reis, R.; Saraiva, R.A.; Villegas, J.; Cruz, M. Water Quality Impact Assessment of Land Use and Land Cover Changes. A dynamic IT model for territorial integrated management. In Proceedings of the GEOProcessing 2015: The Seventh International Conference on Advanced Geographic Information Systems, Applications, and Services, Lisbon, Portugal, 22–27 February 2015; pp. 38–42.
  18. Bičı́k, I.; Jeleček, L.; Štěpánek, V. Land-use changes and their social driving forces in Czechia in the 19th and 20th centuries. Land Use Policy 2001, 18, 65–73. [Google Scholar] [CrossRef]
  19. Regos, A.; Ninyerola, M.; Moré, G.; Pons, X. Linking land cover dynamics with driving forces in mountain landscape of the Northwestern Iberian Peninsula. Int. J. Appl. Earth Obs. Geoinform. 2015, 38, 1–14. [Google Scholar] [CrossRef]
  20. Liu, C.; Ma, X. Analysis to driving forces of land use change in Lu’an mining area. Trans. Nonferrous Met. Soc. China 2011, 21, 727–732. [Google Scholar] [CrossRef]
  21. Teixeira, Z.; Teixeira, H.; Marques, J.C. Systematic processes of land use/land cover change to identify relevant driving forces: Implications on water quality. Sci. Total Environ. 2014, 470–471, 1320–1335. [Google Scholar] [CrossRef] [PubMed]
  22. Ascione, M.; Campanella, L.; Cherubini, F.; Ulgiati, S. Environmental driving forces of urban growth and development. Landsc. Urban Plan. 2009, 93, 238–249. [Google Scholar] [CrossRef]
  23. Melo, J.L.B. Landyn: Land Use and Land Cover Change in Continental Portugal: Characterization, Driving Forces and Future Scenarios. Technical Report for the Landyn Project. Final Report: CO2 Emissions and Removals from the LULUCF; Direção Geral do Território: Lisbon, Portugal, 2014. [Google Scholar]
  24. Meneses, B.M. The Influence of Forest Fire on Water Quality of São Domingos Stream Located in the Western Region of Portugal; High Institute of Agronomy of University of Lisbon: Lisbon, Portugal, 2013. [Google Scholar]
  25. Vale, M.J.; Saraiva, R.A. GI and Watershed Management: Discussing GI Relevance Towards Monitoring EU Legal Framework. In Proceedings of the INSPIRE Conference 2012, Istanbul, Turkey, 23–27 June 2012.
  26. Meneses, B.M. Nuno Cortez Effect of a forest fire on physicochemical properties of the water of São Domingos Stream (Western Region of Portugal). J. Waters Resour. 2015, 36, 1–10. [Google Scholar]
  27. Vale, M.J.; Reis, R.; Meneses, B.M. A caraterização do uso e ocupação do solo de Portugal Continental. In Uso e Ocupação do Solo em Portugal Continental: Avaliação e Cenários Futuros. Projeto LANDYN; Direção-Geral do Território: Lisboa, Portugal, 2014; pp. 1–5. [Google Scholar]
  28. Meneses, B.M.; Saraiva, R.; Reis, R.; Vale, M.J. A importância da avaliação das transições de uso e ocupação do solo na compreensão da evolução da paisagem. In Espaços de Fronteira, Territórios de Esperança: Das Vulnerabilidades às Dinâmicas de Desenvolvimento; Jacinto, R., Ed.; Iberografias, Centro de Estudos Ibéricos: Guarda, Portugal, 2014; pp. 79–92. [Google Scholar]
  29. Srivastava, P.K.; Han, D.; Rico-Ramirez, M.A.; Bray, M.; Islam, T. Selection of classification techniques for land use/land cover change investigation. Adv. Space Res. 2012, 50, 1250–1265. [Google Scholar] [CrossRef]
  30. Paz, J.M.; Sánchez, J.; Visconti, F. Combined use of GIS and environmental indicators for assessment of chemical, physical and biological soil degradation in a Spanish Mediterranean region. J. Environ. Manag. 2006, 79, 150–162. [Google Scholar] [CrossRef] [PubMed]
  31. Verstegen, J.A.; Karssenberg, D.; van der Hilst, F.; Faaij, A.P.C. Identifying a land use change cellular automaton by Bayesian data assimilation. Environ. Model. Softw. 2014, 53, 121–136. [Google Scholar] [CrossRef]
  32. Basse, R.M.; Omrani, H.; Charif, O.; Gerber, P.; Bódis, K. Land use changes modelling using advanced methods: Cellular automata and artificial neural networks. The spatial and explicit representation of land cover dynamics at the cross-border region scale. Appl. Geogr. 2014, 53, 160–171. [Google Scholar] [CrossRef]
  33. García-Frapolli, E.; Ayala-Orozco, B.; Bonilla-Moheno, M.; Espadas-Manrique, C.; Ramos-Fernández, G. Biodiversity conservation, traditional agriculture and ecotourism: Land cover/land use change projections for a natural protected area in the northeastern Yucatan Peninsula, Mexico. Landsc. Urban Plan. 2007, 83, 137–153. [Google Scholar] [CrossRef]
  34. Mitsova, D.; Shuster, W.; Wang, X.H. A cellular automata model of land cover change to integrate urban growth with open space conservation. Landsc. Urban Plan. 2011, 99, 141–153. [Google Scholar] [CrossRef]
  35. Kamusoko, C.; Aniya, M.; Adi, B.; Manjoro, M. Rural sustainability under threat in Zimbabwe—Simulation of future land use/cover changes in the Bindura district based on the Markov-cellular automata model. Appl. Geogr. 2009, 29, 435–447. [Google Scholar] [CrossRef]
  36. Guan, D.; Li, H.; Inohae, T.; Su, W.; Nagaie, T.; Hokao, K. Modeling urban land use change by the integration of cellular automaton and Markov model. Ecol. Model. 2011, 222, 3761–3772. [Google Scholar] [CrossRef]
  37. König, H.J.; Uthes, S.; Schuler, J.; Zhen, L.; Purushothaman, S.; Suarma, U.; Sghaier, M.; Makokha, S.; Helming, K.; Sieber, S.; et al. Regional impact assessment of land use scenarios in developing countries using the FoPIA approach: Findings from five case studies. J. Environ. Manag. 2013, 127, S56–S64. [Google Scholar] [CrossRef] [PubMed]
  38. Chen, X.; Yu, S.X.; Zhang, Y.P. Evaluation of Spatiotemporal Dynamics of Simulated Land Use/Cover in China Using a Probabilistic Cellular Automata-Markov Model. Pedosphere 2013, 23, 243–255. [Google Scholar] [CrossRef]
  39. Santos, P.P.; Tavares, A.O.; Zêzere, J.L. Risk analysis for local management from hydro-geomorphologic disaster databases. Environ. Sci. Policy 2014, 40, 85–100. [Google Scholar] [CrossRef]
  40. Zêzere, J.L.; Vaz, T.; Pereira, S.; Oliveira, S.C.; Marques, R.; Garcia, R.A.C. Rainfall thresholds for landslide activity in Portugal: A state of the art. Environ. Earth Sci. 2014, 73, 2917–2936. [Google Scholar] [CrossRef]
  41. Pereira, S.; Zêzere, J.L.; Quaresma, I.D.; Bateira, C. Landslide incidence in the North of Portugal: Analysis of a historical landslide database based on press releases and technical reports. Geomorphology 2014, 214, 514–525. [Google Scholar] [CrossRef]
  42. Meneses, B.M. Susceptibility and Risk of Landslides in Tarouca County. Master’s Thesis, Institute of Geography and Spatial Planning, University of Lisbon, Lisbon, Portugal, 2011. [Google Scholar]
  43. Tehrany, M.S.; Pradhan, B.; Jebur, M.N. Spatial prediction of flood susceptible areas using rule based decision tree (DT) and a novel ensemble bivariate and multivariate statistical models in GIS. J. Hydrol. 2013, 504, 69–79. [Google Scholar] [CrossRef] [Green Version]
  44. Conforti, M.; Pascale, S.; Robustelli, G.; Sdao, F. Evaluation of prediction capability of the artificial neural networks for mapping landslide susceptibility in the Turbolo River catchment (northern Calabria, Italy). Catena 2014, 113, 236–250. [Google Scholar] [CrossRef]
  45. Brown, C.D.; Davis, H.T. Receiver operating characteristics curves and related decision measures: A tutorial. Chemom. Intell. Lab. Syst. 2006, 80, 24–38. [Google Scholar] [CrossRef]
  46. Vale, M.J. Web Based Colaboratory to Support Water Resources Management and Landuse Integrated Planning; Instituto Superior de Estatística e Gestão de Informação, Universidade Nova de Lisboa: Lisbon, Portugal, 2002. [Google Scholar]
  47. Bonham-Carter, G. Geographic Information Systems for Geoscientists: Modeling with GIS; Elsevier Science Inc. (Pergamon Press): New York, NY, USA, 1994. [Google Scholar]
  48. Yin, K.; Yan, T. Statistical Preditiction Model for Slope Instability of Metamorphosed Rocks. In Proceedings of the Fifth International Symposium in Landslides; Bonnard, C., Ed.; Balkema: Rotterdam, The Netherlands, 1988; pp. 1269–1272. [Google Scholar]
  49. Varnes, D. Slope movement types and processes. In Landslides: Analysis and Control (Special Report—Transportation Research Board, National Research Council; 176); Schuster, R., Krizek, R., Eds.; Transportation Research Board: Washington, DC, USA, 1978; pp. 11–28. [Google Scholar]
  50. Guzzetti, F. Risk-Advanced Weather Forecast System to Advise on Risk Events and Management. Review and Selection of Optional Geological Models Related to Spatial Information Available; IRPI CNR (INTERREG III B CADSES; Europion Union): Perugia, Italy, 2005. [Google Scholar]
  51. Eeckhaut, M.; Reichenbach, P.; Guzzetti, F.; Rossi, M.; Poesen, J. Combined landslide inventory and susceptibility assessment based on different mapping units: An example from the Flemish Ardennes, Belgium. Nat. Hazards Earth Syst. Sci. 2009, 9, 507–521. [Google Scholar] [CrossRef]
  52. Zêzere, J.; Faleh, A.; Sadiki, A.; Garcia, R.; Oliveira, S. Modelação da Susceptibilidade a Deslizamentos na Bacia do Oued Sra (Rif central, Marrocos) com Métodos Estatísticos Bi-variados. In Publicações da Associação Portuguesa de Geomorfólogos, Vol. VI, APGEOM; APGEOM: Braga, Portugal, 2009; pp. 203–210. [Google Scholar]
  53. Zadeh, L.A. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef]
  54. Blahut, J.; Westen, C.; Sterlacchini, S. Analysis of landslide inventories for accurate prediction of debris-flow source areas. Geomorphology 2010, 119, 36–51. [Google Scholar] [CrossRef]
  55. Beven, K.J.; Kirkby, M.J. A physically based variable contributing model of basin hydrology. Hydrol. Sci. Bull. 1979, 24, 43–69. [Google Scholar] [CrossRef]
  56. Yesilnacar, E.; Topal, T. Landslide susceptibility mapping: A comparison of logistic regression and neural networks methods in a medium scale study, Hendek region (Turkey). Eng. Geol. 2005, 79, 251–266. [Google Scholar] [CrossRef]
  57. Chung, C.J.; Fabbri, A.G. Predicting landslides for risk analysis—Spatial models tested by a cross-validation technique. Geomorphology 2008, 94, 438–452. [Google Scholar] [CrossRef]
  58. Druga, M.; Faltan, V. Influences of environmental drivers on land cover structure and its long term changes: A case study of the villages of Malachov and Podkonice in Slovakia. Morav. Geogr. Rep. 2014, 22, 29–41. [Google Scholar]
  59. Meneses, B.M.; Reis, R.; Vale, M.J.; Saraiva, R. Land use and land cover changes in Zêzere watershed (Portugal)—Water quality implications. Sci. Total Environ. 2015, 527–528, 439–447. [Google Scholar] [CrossRef] [PubMed]
  60. Erwin, D.H. Evolutionary uniformitarianism. Dev. Biol. 2011, 357, 27–34. [Google Scholar] [CrossRef] [PubMed]
  61. Baker, V.R. Uniformitarianism, earth system science, and geology. Anthropocene 2014, 5, 76–79. [Google Scholar] [CrossRef]
  62. Dumanski, J.; Pieri, C. Land quality indicators: Research plan. Agric. Ecosyst. Environ. 2000, 81, 93–102. [Google Scholar] [CrossRef]
Figure 1. Zêzere watershed (study area).
Figure 1. Zêzere watershed (study area).
Water 08 00289 g001
Figure 2. Dispersion of new artificial surfaces. Relation between the distance of centroids of these surfaces (2000 and 2012) to principal urban centers (Tomar, Fátima and Covilhã).
Figure 2. Dispersion of new artificial surfaces. Relation between the distance of centroids of these surfaces (2000 and 2012) to principal urban centers (Tomar, Fátima and Covilhã).
Water 08 00289 g002
Figure 3. Independent variables used in modeling artificial surfaces and LUC (CLC) in the years 1990 and 2012.
Figure 3. Independent variables used in modeling artificial surfaces and LUC (CLC) in the years 1990 and 2012.
Water 08 00289 g003
Figure 4. R2 resulting between conditioned probabilities of artificial surfaces of 1990/2000, 2000/2012, and the total period (1990/2012).
Figure 4. R2 resulting between conditioned probabilities of artificial surfaces of 1990/2000, 2000/2012, and the total period (1990/2012).
Water 08 00289 g004
Figure 5. Artificialization surface probability determined by information value and fuzzy gamma, with artificial surfaces from 1990 (the classification of this map follows the method of natural breakdown).
Figure 5. Artificialization surface probability determined by information value and fuzzy gamma, with artificial surfaces from 1990 (the classification of this map follows the method of natural breakdown).
Water 08 00289 g005
Figure 6. Success rate curve for the statistical models used in the determination of artificial surface probabilities.
Figure 6. Success rate curve for the statistical models used in the determination of artificial surface probabilities.
Water 08 00289 g006
Table 1. Areas occupied by artificial surfaces (%) in the Zêzere watershed, from 1990 to 2012.
Table 1. Areas occupied by artificial surfaces (%) in the Zêzere watershed, from 1990 to 2012.
LUC1990200020062012
Continuous urban fabric0.020.020.020.02
Discontinuous urban fabric1.011.131.151.32
Industrial or commercial units0.060.160.250.24
Road and rail networks and associated land0.000.000.010.01
Airports0.010.010.010.03
Mineral extraction sites0.030.050.060.06
Dump sites0.010.010.020.03
Construction sites0.040.050.010.00
Sport and leisure facilities0.010.010.010.01
Table 2. Area of artificial surfaces (AS), conditional probability (CPji), information value (IV), and fuzzy membership value (FM) for each class of the independent variables, in the Zêzere watershed.
Table 2. Area of artificial surfaces (AS), conditional probability (CPji), information value (IV), and fuzzy membership value (FM) for each class of the independent variables, in the Zêzere watershed.
VariablesClassesTotal Area (ha)AS 1990 (ha)AS 2000 (ha)AS 2012 (ha)CPji (AS 1990)CPji (AS 2000)CPji (AS 2012)IV (AS 1990)FM (AS 1990)
AspectNorthwest51,549.15338.96409.33482.170.0070.0080.009−0.5930.007
West56,574.22411.9489.83549.450.0070.0090.010−0.4910.008
Southwest61,918.88705.13848.98967.340.0110.0140.016−0.0440.013
South57,996.57884.721050.811209.060.0150.0180.0210.2480.017
Southeast61,927.05786.6953.711138.820.0130.0150.0180.0650.014
East54,624.12642.44787.59961.530.0120.0140.018−0.0120.013
Northeast48,357.44525.06657.68791.940.0110.0140.016−0.0920.012
North42,475.03300.68379.59458.830.0070.0090.011−0.520.008
Flat66,855.891382.431682.282028.420.0210.0250.0300.5520.023
Hillshade250–300221.623.303.153.150.0150.0140.0140.2240.016
200–25089,250.44259.45314.87352.120.0030.0040.004−14100.003
150–200306,686.85193.936328.975620.0170.0210.0250.3530.019
100–15085,814.72475.65559.81614.410.0060.0070.007−0.7640.006
50–10018,890.3943.4150.5453.170.0020.0030.003−16450.003
0–501414.382.182.532.710.0020.0020.002−20440.002
Humidity (%)75–8059,071.721432.731775.392442.130.0240.0300.0410.7120.027
70–75249,233.663745.744436.335012.90.0150.0180.0200.2330.017
65–70171,005.31745.06993.691051.270.0040.0060.006−10050.005
<6522,967.6654.3954.3981.260.0020.0020.004−16150.003
Insolation (Hours)2700–28004373.774.9583.4683.460.0170.0190.0190.3650.019
2600–2700149,922.622085.832603.913409.020.0140.0170.0230.1560.015
2500–2600123,480.1919.68986.321051.590.0070.0080.009−0.4690.008
2400–2500117,329.611419.261629.111747.500.0120.0140.0150.0160.013
2300–2400105,098.181478.21957.002295.990.0140.0190.0220.1670.015
2200–23002074.14000000−0.469 *0
Precipitation (mm)2400–28005012.8135.0635.0635.060.0070.0070.007−0.5320.008
2000–240010,286.57132.26132.26131.620.0130.0130.0130.0770.014
1600–200020,541.19430.63506.08531.680.0210.0250.0260.5660.023
1400–160066,379.08550.95769.49916.670.0080.0120.014−0.360.009
1200–1400132,983.261623.232009.962402.330.0120.0150.0180.0250.013
1000–120090,577.8945.341115.631311.520.010.0120.014−0.1310.011
800–1000130,193.41226.41474.781642.620.0090.0110.013−0.2340.01
700–80039,288.9677.81804.971048.720.0170.0200.0270.3710.019
600–7007015.34356.24411.57567.340.0510.0590.0811,4510.056
Slope (%)40–45175000000−0.967 *0
35–40874.97000000−0.967 *0
30–352670.41.351.351.350.0010.0010.001−31590.001
25–307100.7733.5533.5533.550.0050.0050.005−0.9240.005
20–2516,322.4361.6273.1672.730.0040.0040.004−11480.004
15–2037,573.66264.08276.15290.430.0070.0070.008−0.5270.008
10–1579,761.58360.91421.9439.290.0050.0050.006−0.9670.005
5–10156,451.831267.971596.841863.890.0080.0100.012−0.3840.009
0–5201,347.713988.444856.855886.320.0200.0240.0290.5090.022
SoilHumic Cambisols110,548.72708.37800.85819.070.0060.0070.007−0.6190.007
Rankers12,992.1531.4831.4831.480.0020.0020.002−15920.003
Dystric Cambisols48,878.67581.17829.41067.760.0120.0170.022−0.0010.013
Dystric Fluvisols3017.16.396.396.390.0020.0020.002−17260.002
Eutric Lithosol187,479.931475.361776.122140.970.0080.0090.011−0.4140.009
Calcic Cambisols8275.4977.7486.9886.980.0090.0110.011−0.2370.01
Calcic Luvisols34,032.091300.161509.651793.440.0380.0440.05311660.042
Hortic Luvisols40,991.55320.14376.3421.690.0080.0090.010−0.4210.009
Calcic-chromic Cambisols13,019.96169.54215.53215.550.0130.0170.0170.090.014
Eutric Cambisols30,585.641213.241512.21792.560.040.0490.05912040.044
Chromic Cambisols6065.9380.23100.8197.570.0130.0170.0330.1060.015
Hortic Podzols6390.7913.813.7713.770.0020.0020.002−17070.002
Eutric Fluvisols0.330.30.330.330.9091.0001.00043361000
Temperature (°C)16.0–17.583,699.312360.492918.13851.40.0280.0350.0460.8630.031
15.0–16.060,518.29819.04918.61929.710.0140.0150.0150.1290.015
12.5–15.0135,430.82977.721250.241320.630.0070.0090.010−0.50.008
10.0–12.5135,296.25814.51994.851122.50.0060.0070.008−0.6820.007
7.5–10.072,973.84935.961094.531279.850.0130.0150.0180.0750.014
<7.514,359.8470.283.4783.470.0050.0060.006−0.890.005
DWBW (km)5–6135.231.2611.5211.520.0090.0850.085−0.2450.01
4–52872.22140.47226.48236.480.0490.0790.08214130.054
3–424,827.85458.91563.72676.750.0180.0230.0270.440.02
2–380,622.991053.361258.361404.170.0130.0160.0170.0930.014
1–2154,963.461443.751803.562327.760.0090.0120.015−0.2450.01
0–1238,856.62880.173396.163930.880.0120.0140.0160.0130.013
TWI>25910.2412.5210.4810.430.0140.0120.0110.1450.015
20–252612.1743.2247.553.520.0170.0180.0200.330.018
15–2018,209.96297.73348.24401.760.0160.0190.0220.3180.018
10–15174,288.63317.834030.444772.980.0190.0230.0270.470.021
5–10306,257.382306.622818.863342.930.0080.0090.011−0.4570.008
Note: * These values correspond to the smaller IV observed in the variable under analysis.
Table 3. Statistical description of information value (IV) and fuzzy gamma (FG) by sectors (S) of the Zêzere watershed (significance level p < 0.05).
Table 3. Statistical description of information value (IV) and fuzzy gamma (FG) by sectors (S) of the Zêzere watershed (significance level p < 0.05).
DescriptionIVFG
S (A)S (B)S (C)S (A)S (B)S (C)
Min.−8.89−8.89−4.49000.0014
Max.3.344.2510.410.00510.00670.0082
Mean−2.58−1.671.880.00200.00230.0041
Std. Dev.1.841.981.860.00070.00080.0011
Table 4. Accountability ( A I ) and reliability ( R I ) indexes by sectors (S) of the Zêzere watershed.
Table 4. Accountability ( A I ) and reliability ( R I ) indexes by sectors (S) of the Zêzere watershed.
VariablesS (A)S (B)S (C)
A I R I A I R I A I R I
Aspect59.91.054.20.954.53.8
Hillshade72.51.191.11.094.63.4
Humidity98.51.567.50.944.85.0
Insolation84.41.565.01.336.86.5
Precipitation44.12.248.01.056.45.7
Slope77.80.969.41.287.23.4
Soil30.21.269.33.581.33.9
Temperature47.31.686.70.874.73.6
DWBW26.91.354.70.867.93.9
TWI45.01.064.21.270.44.0

Share and Cite

MDPI and ACS Style

Meneses, B.M.; Reis, E.; Vale, M.J.; Reis, R. Modeling the Probability of Surface Artificialization in Zêzere Watershed (Portugal) Using Environmental Data. Water 2016, 8, 289. https://doi.org/10.3390/w8070289

AMA Style

Meneses BM, Reis E, Vale MJ, Reis R. Modeling the Probability of Surface Artificialization in Zêzere Watershed (Portugal) Using Environmental Data. Water. 2016; 8(7):289. https://doi.org/10.3390/w8070289

Chicago/Turabian Style

Meneses, Bruno M., Eusébio Reis, Maria J. Vale, and Rui Reis. 2016. "Modeling the Probability of Surface Artificialization in Zêzere Watershed (Portugal) Using Environmental Data" Water 8, no. 7: 289. https://doi.org/10.3390/w8070289

APA Style

Meneses, B. M., Reis, E., Vale, M. J., & Reis, R. (2016). Modeling the Probability of Surface Artificialization in Zêzere Watershed (Portugal) Using Environmental Data. Water, 8(7), 289. https://doi.org/10.3390/w8070289

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop