Land Use / Land Cover Assessment over Time Using a New Weighted Environmental Index (WEI) Based on an Object-Oriented Model and GIS Data

: For the ﬁrst time, this paper introduces and describes a new Weighted Environmental Index (WEI) based on object-oriented models and GIS data. The index has been designed to integrate all the available information from extensive and detailed GIS databases. After the conceptual deﬁnition of the index has been justiﬁed, two applications for the regional and local scales of the WEI are shown. The applications analyze the evolution over time of the environmental value from land-use change for two di ﬀ erent case studies in Spain: the Valencian Region and the L’Alcora municipality. Data have been obtained from the Spanish Land Occupation Information System (SIOSE) public database and integrate GIS information about land use / land cover on an extensive, high-detailed scale. Results demonstrate the application of the WEI to real case studies and the importance of integrating statistical analysis of WEI evolution over time to arrive at a better understanding of the socio-economic and environmental processes that induce land-use change.


Introduction
Current efforts to define and establish environmental indicators stem from early debates about sustainability, popularized by the "Brundtland Report" and "Agenda 21" [1,2]. In recent years, environmental indicators have become a fundamental tool in environmental assessment at detailed, local, regional and national levels [3]. These environmental indicators significantly influence environmental management and the formulation of environmental policies [4], as well as monitoring and evaluation processes [5,6]. Environmental indicators have become important because they provide a sign that conveys a complex message, potentially from numerous sources, in a simplified and useful way [7]. Therefore, there is a growing need to establish appropriate environmental indicators based on truthful and verifiable information [8,9].
The use of indicators to emphasize the relevance of environmental data has many advantages. It is easier for scientists to understand and agree on the existing purpose of a particular monitoring program, and clients clearly understand what they are paying for and why [10]. It also helps to understand interactions between different groups of measurements, such as pollutant inputs, concentrations in sentinel organisms and biological effects [11]. However, it is a great challenge to determine which of the numerous measures of environmental systems characterize the entire system but are simple enough to be effectively and efficiently monitored and modeled [12].
Current and future work on environmental indicators should consider the following aspects [13]: (i) indicators are the product of numerous measuring processes that oversimplify environmental trends, while ignoring important social and political factors produced by the indicators; (ii) the establishment Composite indicators or indices are a mathematical combination of a set of simple indicators that summarize a multidimensional concept in a simple or unidimensional index based on an underlying conceptual model. They can be quantitative or qualitative, depending on the analyst's requirements [38].
Similarly, in [39], compound indicators are defined as mathematical combinations of simple indicators that do not have a common unit of measurement.
The number of composite indicators grows every year, and they are applied in different areas of interest, since they have the ability to explain complex concepts [39]. The European Union and the Organization for Economic Cooperation and Development (OECD) are pioneering organizations in the development of initiatives related to these concepts in different fields of study (innovation/technology, society, globalization, environment, economy, etc.), generating a collection of documentation that serves as a starting point for its study [38]. Organizations such as the United Nations and the European Commission have developed highly interesting composite indicators [40][41][42][43][44][45][46].
The increasing number of these tools is a clear symptom of their political importance and their operational relevance in decision-making [60][61][62].
However, there are certain limitations of composite indicators that must be known in order to improve their design and avoid possible criticism about their construction [20,39]. Table 1 lists the main advantages and disadvantages of composite indicators [38,39]. Table 1. Main advantages and disadvantages of the composite indicators.

Advantages Disadvantages
• Can summarize complex, multi-dimensional realities with a view to support decisionmakers. • Are easier to interpret than a battery of many separate indicators.

•
Can assess progress of countries over time.

•
Reduce the visible size of a set of indicators without dropping the underlying information base.

•
Make it possible to include more information within the existing size limit. • Place issues of country performance and progress at the center of the policy arena. Some of the best-known composite indicators or environmental indices are: the Environmental Performance Index [9,63], Environmental Vulnerability Index [64], Living Planet Index [65] and Ecological Footprint [66].
Most composite environmental indicators measure multidimensional concepts about a group of countries for later comparison between them (rankings), being difficult to find environmental indices on smaller scales.

Environmental Indicators Based on Land Use
Land use is one of the main causes of the transformations of the terrestrial ecosystem [67,68], and these changes affect climate [69], biodiversity [70] and landscape ecology [71,72]. Furthermore, it is also a key factor in the context of policy and reporting schemes [73].
Due to these different approaches, the terminology related to land use is diverse. For example, two basic definitions, land cover and land use, are often mixed up or used synonymously. "Land cover" refers to the physical material on the surface of the land, while "land use" most often refers to the functional dimension and describes how the area is used for urban, agricultural, forestry and other purposes.
Two main approaches can be found in the literature that address environmental indicators based on land use [74]. The first has to do with the amount of land occupied and/or transformed for a certain activity. The second is based on models that quantify impacts in terms of variation in soil properties.
The amount of land occupied and/or transformed is multiplied by characterization factors that reflect changes in soil properties for each type of land occupation and transformation. Soil quality can be measured using a single indicator, such as soil organic matter [75], soil organic carbon [76], soil erosion [77], or multiple indicators (for example, using LANCA ® models (Land Use Indicator Value Calculation in Life Cycle Assessment) [78] and LUCI (Land Utilization and Capability Indicator) [79].
But the difficulty of having the necessary data to apply one or another indicator on many occasions means that each country implements its own method of evaluating land use [49,80]. Continuously and depending on the information available, new indicators and indices based on land use are appearing [81,82].
Recent studies show the importance of using updated data in order to study the territory with environmental indexes and indicators related to land cover/use or landscape ecology. Practically all of these investigations are supported by the use of official databases of the area to be analyzed.
For example, in [83], the landscape composition and configuration changes of an area in southern Ecuador were evaluated with data from the United States Geological Survey (USGS). To perform a spatiotemporal analysis of land-use and land-cover change in the Atlantic forests of Brazil [84], obtained available data from a free-access dataset developed by a consortium of Brazilian and international research institutes, universities, private organizations, and NGOs aiming to generate national coverage of land cover/use information. In China, national databases are also used to study environmental aspects based on land use/land cover changes [85]. The continuous acquisition of images by orbital satellites also provides enough information to be able to evaluate land cover change and habitat configuration from the visual comparison of the generated maps [86][87][88].

Materials and Methods
Environmental indices and indicators based on GIS analysis of land use/land cover are usually based on qualitative values of different parameters assigned to the different plots into which the study area has been divided [68,83,84,86,88]. In these situations, the number of soil categories is usually small (between five to eleven), and it is easy to make comparisons between corresponding maps from different years. However, these indices do not allow for the performing quantitative assessments of environmental quality changes, beyond analysis of land uses per unit area. More elaborated environmental indices quantify environmental value based on land use changes. This is the case in [81], which studies the variation in the degree of anthropization through analysis of a single factor called the Relative Integrated Anthropization Index (INRA). This index includes five different categories of soil anthropization, with relative values changing between 0 and 1.
Subsequent studies have modified the INRA index to include more categories. In [89], 27 additional subcategories are added, providing a greater degree of detail in the definition of the INRA index. The relative value of anthropization assigned to each of these new categories is obtained as a result of applying multiple-criteria decision analysis, obtained by expert judgement.
Recent studies focus on the application of fuzzy logic methodologies for the evaluation of land use. The overall suitability assessment of land units is based on the definition of weighting factors of the relevant characteristics. In these methodologies, the choice of weight values is of critical importance. These weights are usually decided on the basis of expert knowledge considering local advice, experimental data or previous land evaluation methods [90].
This paper introduces a new environmental index named WEI (Weighted Environmental Index), based on land use analysis techniques, that allows all the information obtained from official public databases to be integrated on a detailed scale. The WEI definition and its application to a practical case is shown below.

Description of Mapping Techniques Using GIS
Land occupation mapping is thematic mapping that represents two distinct but interrelated components: the occupation of land surface according to its biophysical properties, called land cover, and the characterization of the territory according to its socioeconomic dedication, called land use. Therefore, mapping of land uses and covers involves the natural and socioeconomic factors used in a given space and time [88]. Its main objective is the planning and monitoring of resources, such as the changes that affect natural cover caused by land use. This process is generally driven by natural phenomena and anthropogenic activities, which affect the natural ecosystem.
Geographical information systems (GIS) are useful tools for the elaboration of land use and cover mapping. They offer an important advantage by integrating different information technologies such as remote sensing (RS) and global navigation satellite systems (GNSS). This integration allows for a more efficient study of changes and the quantification of the dynamics of land use [91]. The National Geographic Institute of Spain (IGN) coordinates, in parallel, two projects aimed at structuring information on land use and coverage: Corine Land Cover and SIOSE.

Corine Land Cover
The Corine Land Cover project, established by the European Union and coordinated by the European Environment Agency, is aimed at homogenizing information on all of Europe, to facilitate the performance of territorial analysis, the state of the environment and natural resources, and the establishment of European policies. This information, integrated into an international geographic information system, is structured as a hierarchical data model of 44 classes at level 3 and 58 classes at level 4. It is defined on a 1:100,000 spatial scale and has a minimum polygonal surface of 25 hectares. The first product was obtained in 1990 and later updates have been made in 2000, 2006, 2012 and 2018 [92,93].

SIOSE
The need for more detailed information on land use on a national scale led to the development of the Spanish Land Occupation Information System (SIOSE) in 2005. This system has been structured as an object-oriented conceptual data model, with 40 simple classes and 46 compound classes; its spatial scale is 1:25,000. Its work unit is polygons, with a minimum mappable surface of 2 hectares for agricultural, forest and natural areas; 1 ha for urban areas and bodies of water; and 0.5 hectares for crops. The conceptual model includes two super classes: land use and cover. The coverage can be of a simple type when it is unique within the polygon, or composite when the polygon includes two or more types of simple or composite coverage in turn. However, land use refers to the type of socioeconomic activity and does not necessarily correspond to a physical aspect. For example, a forest cover can have one more use types of use (recreational or/and economic [93]). In this way, the SIOSE model does not describe a single coverage for each polygon but can assign one or multiple simple or composite coverages for a single polygon through its attributes and occupancy percentages. Consequently, SIOSE offers more detailed information geared towards user needs [92] (Figure 1).

Figure 1.
Example of classification using a hierarchical data modeling system and the Spanish Land Occupation Information System (SIOSE) system.

The Weighted Environmental Index (WEI). Conceptual scheme.
This work introduces the Weighted Environmental Index (WEI), a new index for the analysis of environmental value based on land use.
This new index has been defined to fulfill the following characteristics: 1. WEI must integrate all the characteristics of indices that vary continuously in space. 2. WEI values should be justified in a simple way from pre-established classifications of land use. 3. It must be able to be used to carry out land use assessments based on information integrated into geographic information systems (GIS). 4. It must be able to be used both in general studies carried out on a large scale and in detailed studies that use cartography obtained by very high-resolution GIS techniques. 5. Its application in the same geographical area at different times should allow for trend analysis to determine the impact of correction measures that are implemented through territorial, urban or environmental planning tools.
The process of determining environmental index values for each land use (EIj) has been carried out, taking into account the joint consideration of the following five evaluation factors (Fi): F1: Anthropic or natural nature of activity developed in soil. F2: Water consumption associated with land use. F3: Soil degradation (use of chemicals). F4: Environmental sustainability of land use (stability of the ecosystem). F5: Landscape value of activity carried out in the analyzed area.
The determination of values of the evaluation factors for each land use is carried out individually so that a quantitative value is assigned for each factor Fi and land use j, in such a way that: For each one of the land use categories included in the SIOSE legend, the corresponding environmental index (EIj) has been obtained as the value of the weighted average of each of the values

The Weighted Environmental Index (WEI). Conceptual scheme.
This work introduces the Weighted Environmental Index (WEI), a new index for the analysis of environmental value based on land use.
This new index has been defined to fulfill the following characteristics: 1. WEI must integrate all the characteristics of indices that vary continuously in space.

2.
WEI values should be justified in a simple way from pre-established classifications of land use.

3.
It must be able to be used to carry out land use assessments based on information integrated into geographic information systems (GIS).

4.
It must be able to be used both in general studies carried out on a large scale and in detailed studies that use cartography obtained by very high-resolution GIS techniques.

5.
Its application in the same geographical area at different times should allow for trend analysis to determine the impact of correction measures that are implemented through territorial, urban or environmental planning tools.
The process of determining environmental index values for each land use (EI j ) has been carried out, taking into account the joint consideration of the following five evaluation factors (F i ): F 1 : Anthropic or natural nature of activity developed in soil. F 2 : Water consumption associated with land use. F 3 : Soil degradation (use of chemicals). F 4 : Environmental sustainability of land use (stability of the ecosystem). F 5 : Landscape value of activity carried out in the analyzed area.
The determination of values of the evaluation factors for each land use is carried out individually so that a quantitative value is assigned for each factor F i and land use j, in such a way that: For each one of the land use categories included in the SIOSE legend, the corresponding environmental index (EI j ) has been obtained as the value of the weighted average of each of the values assigned to each of the previous factors (F i ), considering the corresponding weights (α i ), as shown in Equations (2) and (3): where: EI j : environmental index of land use j (0 ≤ E ij ≤ 100) α i : assigned weights to factor i F i : evaluation factor i ncat: land use categories The application of environmental index values (EI j ) is carried out on a discretization in irregular polygons of variable surface that together constitute the entire area under study: where: A total : total area of study A k : area of polygon k npol: total number of polygons in the discretization Therefore, once the values of the environmental indices corresponding to each land use have been established, the weighted environmental index of a certain polygon (WEI k ) is determined based on the values of the environmental index of each land use included inside the polygon, considering as weights the proportion of the area assigned to each land use with respect to the total area of the polygon, as shown in Equation (6). where: WEI k : weighted environmental index of polygon k. EI j : environmental index of land use j. A jk : area assigned to land use j inside the polygon k. β jk : land use weighting factor j in polygon k.
n jk : number of land uses (j) inside polygon k.
The value of the weighted environmental index obtained by Equation (6) adopts values that vary between 0 and 100, so that values close to 0 indicate a very low environmental value, while values close to 100 indicate a high environmental value. This is in accordance with the five evaluation factors (F i ) considered in the definition of the environmental index for each type of land use.
Thus, WEI k values are determined from EI j values, which depend on the values assigned to the evaluation factors (F i ) and the weights associated with each factor (α i ). Therefore, the value of the WEI index depends on the corresponding values of the evaluation factors (F i ) and their corresponding weights (α i ). The values of F i and α i should be decided by the modeler on the basis of expert knowledge considering local advice. Table 2 shows the values of the weighted environmental index for each land use (WEI k ) included in the SIOSE legend as a result of the linear combination of the five evaluation factors (F i ) considered in the definition of the index. The values of each environmental factor are the ones that have been used for demonstration purposes in the two case studies shown in this paper (Valencia Region and L'Alcora municipality), which are described in detail below. In both case studies, equal values of the weights associated with each factor (α i = 0.2) have been considered.
The values shown in Table 2 can be modified or adapted by the user in each case. The user is responsible for justifying the values of F i and α i , for which the existence of particular conditions in the area under study that could modify the proposed values must be taken into account. These values have been designed so final results in terms of the WEI index allow the evolution of environmental value on a certain region to be studied, prioritizing the natural uses of the soil with low water consumption, low soil degradation, high sustainability of the ecosystem and high landscape value, following a Multi-Criteria Decision Analysis (MCDA) technique. MCDA is currently used to establish the value of environmental indicators [94][95][96][97]. When analyzing a territory, it is necessary to take into account that sustainability assessment is a multi-criteria decision process that comprises of economic, social, and environmental practice [98]. The purpose of MCDA is to compare and rank alternative options and to evaluate environmental consequences according to the criteria established [99]. One of its greatest strengths is the possibility of using the criteria with their own dimensions. The main weakness of MCDA is the subjectivity of the weighting step that is needed to value the different criteria [100].   Commercial and offices  20  20  20  20  20  20  TCH  Hotels  20  20  20  20  20  20  TPR  Recreational park  20  20  20  20  20  20  TCG  Camping  20  40  40  50  50  40  EAI  Institutional administrative  20  20  20  20  20  20  ESN  Medical and sanitary  20  20  20  20  20  20  ECM  Cementery  20  20  20  20  20  20  EDU  Education  20  20  20  20  20  20  EPN  Penitentiary  20  20  20  20  20  20  ERG  Religious  20  20  20  20  20  20  ECL  Cultural  20  20  20  20  20  20  EDP  Sport  25  15  20  20  20  20  ECG  Golf course  40  10  70  50  80  50  EPU  Urban park  60  65  70  80  75   The application of the WEI to each land use considered by the SIOSE legend allows a classification to be established based on the discrimination by ranges shown in Table 3:

Results and Discussion
Two applications of the WEI for regional analysis (Valencian Region) and for detailed analysis at a municipal level (L'Alcora) are shown below. Data used for these studies included the SIOSE information available for the Valencia Community in 2005, 2009 and 2015 downloaded from the Spanish National Geographic Institute (IGN) platform with geodetic reference system ETR89 and transverse Mercator universal projection geographic system (UTM) in time zone 30. This information was based on the photointerpretation of SPOT5 images, orthophotos from the National Plan for Aerial Orthophotography (PNOA), IGN official cartography databases and information provided by the Autonomous Community on land uses at different scales (SIOSE2015).
Using spatial analysis data techniques, a report was generated from the results obtained from a query to the SIOSE data model in which information on the percentage of occupation, the surface in hectares and type of ground cover was obtained for each polygon mapped the Valencian Region in 2005, 2009 and 2015. This information was the basis for the analysis and generation of the WEI, and later it was linked to the GIS through a polygon identifier for spatial representation of the evolution of this index on the three considered dates. The values of the evaluation factors (Fi) and the weights (α i ) are the same in both cases. The evaluation factors used in both analyses are shown in Table 2, and they were selected by expert judgment considering equal weight to every evaluation factor.

Large Scale Analysis: Valencian Region (2005-2015)
The Valencian Region is one of the 17 autonomous communities into which Spain is divided. It is made up of three provinces (Alicante, Castellón and Valencia), with its capital in the city of Valencia. It is located to the east of the Iberian Peninsula, on the coast of the Mediterranean Sea. It has an area of 23,255 km 2 and a population of 5,003,769 inhabitants, according to the National Statistical Institute (INE). Its economy is based mainly on the service sector (70%), followed by industry, construction and finally agriculture. The latter has significantly lost its importance in the last five decades [101]. To analyze the evolution of environmental status at a regional level, the new WEI has been applied to the Valencian Region for which an extensive set of data obtained from the SIOSE database was available. These data referred to the land use on the whole territory for the 2005-2015 period. Following the methodology explained in Section 3.2, the WEI for every polygon has been computed.  Results in Figure 2 show that the lowest environmental value (WEI < 50) is found in the surroundings of Valencia city, located at the center of the eastern coast of the study area. Higher WEI values (WEI > 80) are almost exclusively found on the non-altered mountains located inland where human actions are not relevant.
The application of the WEI to the Valencian Region demonstrates that WEI can be successfully used to analyze environmental status at a regional level if enough accurate data are available.
Additionally, WEI allows us to analyze the evolution over time of environmental value. Figure  3 shows the WEI differences map between 2005 and 2015 for the whole Valencian Region. Results in Figure 2 show that the lowest environmental value (WEI < 50) is found in the surroundings of Valencia city, located at the center of the eastern coast of the study area. Higher WEI values (WEI > 80) are almost exclusively found on the non-altered mountains located inland where human actions are not relevant.
The application of the WEI to the Valencian Region demonstrates that WEI can be successfully used to analyze environmental status at a regional level if enough accurate data are available.
Additionally, WEI allows us to analyze the evolution over time of environmental value. Figure 3 shows the WEI differences map between 2005 and 2015 for the whole Valencian Region. In Figure 3, green values represent areas for which the WEI has increased during this ten-year period, while red values represent areas where the WEI has decreased. As a result of the analysis, a large area located inland at a western region from Valencia city where WEI values had decreased was detected. A specific analysis was performed to investigate the reasons for this sudden change, leading to the conclusion that the WEI had been affected by a change in the criteria used by SIOSE to map and define the land use of polygons of this area. Regarding the accuracy of the data, WEI acts as a tool for checking the land use databases provided by SIOSE. The irregularities detected are due to (i) recoding of land use, (ii) grouping of polygons or (iii) errors in the database for specific years.
It has been precisely when applying the WEI to the Valencian Region that inconsistencies in the numerical values of the SIOSE database, which were unnoticed before, have been detected.
Furthermore, the spatial distribution of WEI differences allows us to easily identify the position of areas that have improved or worsened their environmental value during the ten-year period. A statistical analysis can be performed, comparing the WEI value inside each polygon in which the area is discretized.
In addition to the results and cartography shown, which are both of great value in the visualization of the results, statistical analysis of the temporal evolution of the WEI provides very valuable results from the point of view of land use management and the impact of the policies implemented in the territory. Due to the extremely large size of the Valencian Region, an example of this statistical analysis is shown below, applying the WEI to the municipality of L'Alcora. In Figure 3, green values represent areas for which the WEI has increased during this ten-year period, while red values represent areas where the WEI has decreased. As a result of the analysis, a large area located inland at a western region from Valencia city where WEI values had decreased was detected. A specific analysis was performed to investigate the reasons for this sudden change, leading to the conclusion that the WEI had been affected by a change in the criteria used by SIOSE to map and define the land use of polygons of this area. Regarding the accuracy of the data, WEI acts as a tool for checking the land use databases provided by SIOSE. The irregularities detected are due to (i) recoding of land use, (ii) grouping of polygons or (iii) errors in the database for specific years.
It has been precisely when applying the WEI to the Valencian Region that inconsistencies in the numerical values of the SIOSE database, which were unnoticed before, have been detected.
Furthermore, the spatial distribution of WEI differences allows us to easily identify the position of areas that have improved or worsened their environmental value during the ten-year period. A statistical analysis can be performed, comparing the WEI value inside each polygon in which the area is discretized.
In addition to the results and cartography shown, which are both of great value in the visualization of the results, statistical analysis of the temporal evolution of the WEI provides very valuable results from the point of view of land use management and the impact of the policies implemented in the territory. Due to the extremely large size of the Valencian Region, an example of this statistical analysis is shown below, applying the WEI to the municipality of L'Alcora.

Municipal Scale Analysis: L'Alcora municipality (2005-2015)
Using the same values of the evaluation factors (F i ) and weights (α i ), the WEI has been applied to analyze land use evolution over time at a local level in the municipality of L'Alcora (Spain). The municipality of L'Alcora is located within the province of Castellón with an area of 95.26 km 2 . It has a population of 10,405 inhabitants, and its main economic sector is the ceramic industry [102].
The municipality of L'Alcora has been selected to verify the suitability of using the WEI at a municipal level, as its socioeconomic structure includes a large number of land uses associated with the territorial distribution of urban use, industrial use, agricultural use and forestry use. In this way, in the municipality of L'Alcora, it is of great interest to carry out analysis of the evolution of land use to assess the relationships between the growth (or maintenance) of industrial use, economic development and employment of the area and the status and environmental value of the territory.
Similar to how it was done in Section 4.1 and following the methodology explained in Section 3.2, the WEI for every polygon has been computed. Figure 4 shows the spatial distribution of the WEI for 2005, 2009 and 2015. The spatial distributions of the WEI allow us to easily identify the position of the areas of highest and lowest environmental value within the municipality.

Municipal Scale Analysis: L'Alcora municipality (2005-2015)
Using the same values of the evaluation factors (Fi) and weights (α ), the WEI has been applied to analyze land use evolution over time at a local level in the municipality of L'Alcora (Spain). The municipality of L'Alcora is located within the province of Castellón with an area of 95.26 km 2 . It has a population of 10,405 inhabitants, and its main economic sector is the ceramic industry [102].
The municipality of L'Alcora has been selected to verify the suitability of using the WEI at a municipal level, as its socioeconomic structure includes a large number of land uses associated with the territorial distribution of urban use, industrial use, agricultural use and forestry use. In this way, in the municipality of L'Alcora, it is of great interest to carry out analysis of the evolution of land use to assess the relationships between the growth (or maintenance) of industrial use, economic development and employment of the area and the status and environmental value of the territory.
Similar to how it was done in Section 4.1 and following the methodology explained in Section 3.2, the WEI for every polygon has been computed. Figure 4 shows the spatial distribution of the WEI for 2005, 2009 and 2015. The spatial distributions of the WEI allow us to easily identify the position of the areas of highest and lowest environmental value within the municipality.  Roughly speaking, low environmental value areas (WEI < 40) correspond to urban use and roads, while high environmental value areas (WEI > 80) correspond to forest land. It has been observed that low-WEI areas tend to concentrate along the main road that crosses the municipality.
Unlike the case of the Valencian Region shown in Section 4.1, the adequate size of the study area allows us to perform a detailed statistical analysis that is now useful to better understand the evolution over time of the WEI inside the L'Alcora municipality.
Evolution over time of WEI values can be analyzed by comparing the values of the WEI for every polygon. Results are shown in Figure 5, which shows the WEI differences map between 2005 and 2015 for the municipality of L'Alcora. Though no changes are observed along the main road, a decrease in WEI values has been detected in agricultural and forest land disseminated throughout the territory. This fact can be numerically objectivized by computing the average value of the WEI inside the study area as the weighted average of the WEI inside each polygon, considering each polygon's area as weights.
Sustainability 2020, 12, x FOR PEER REVIEW  14 of 23 Roughly speaking, low environmental value areas (WEI < 40) correspond to urban use and roads, while high environmental value areas (WEI > 80) correspond to forest land. It has been observed that low-WEI areas tend to concentrate along the main road that crosses the municipality.
Unlike the case of the Valencian Region shown in Section 4.1, the adequate size of the study area allows us to perform a detailed statistical analysis that is now useful to better understand the evolution over time of the WEI inside the L'Alcora municipality.
Evolution over time of WEI values can be analyzed by comparing the values of the WEI for every polygon. Results are shown in Figure 5, which shows the WEI differences map between 2005 and 2015 for the municipality of L'Alcora. Though no changes are observed along the main road, a decrease in WEI values has been detected in agricultural and forest land disseminated throughout the territory. This fact can be numerically objectivized by computing the average value of the WEI inside the study area as the weighted average of the WEI inside each polygon, considering each polygon's area as weights.      Table 4 shows the values of the deciles of the WEI distribution functions for each one of the years under study and the summary of the basic statistics. Figure 7 graphically shows this information as the correspondent cumulative distribution functions (CDFs).   Table 4 shows the values of the deciles of the WEI distribution functions for each one of the years under study and the summary of the basic statistics. Figure 7 graphically shows this information as the correspondent cumulative distribution functions (CDFs). As no dramatic changes in land use have been observed during the period 2005-2015, the shape of the CDFs is similar for the three different dates. As expected, the shape of the CDFs shows a trend in high WEI values, which is in accordance with the high Average WEI value obtained before, which was higher than 73 for every year. It has been observed, though, that the largest area for low WEI values was found in 2009. This fact is also seen when comparing the shape of the CDFs, computing the differences between the correspondent deciles for each year. Table 5 and Figure 8 show the details of these calculations and their graphical representation, respectively. As no dramatic changes in land use have been observed during the period 2005-2015, the shape of the CDFs is similar for the three different dates. As expected, the shape of the CDFs shows a trend in high WEI values, which is in accordance with the high Average WEI value obtained before, which was higher than 73 for every year. It has been observed, though, that the largest area for low WEI values was found in 2009. This fact is also seen when comparing the shape of the CDFs, computing the differences between the correspondent deciles for each year. Table 5 and Figure 8 show the details of these calculations and their graphical representation, respectively.   Analysis of the evolution over time of the deciles of the CDFs is shown in Figure 8 and lead to interesting results. Figure 8a shows the CDF's decile values. The maximum value for each year (2005,2009 and 2015) has always been obtained for the WEI class [70,80[, but only in 2015, this class reached values higher than 40%. To fully understand the evolution over time of environmental value inside the L'Alcora municipality using the WEI, Figure 8b,d must be analyzed carefully. These figures show the incremental analysis of the CDF's decile evolution through the analyzed time period, allowing us to obtain specific results for the L'Alcora municipality. Figure 8b shows that between 2005 and 2009, a significant loss of environmental value was observed. The lowest WEI class [10,20[ increased its area by 1.81% due to the loss of area of other low-WEI classes. As shown in Figure 8c,   Analysis of the evolution over time of the deciles of the CDFs is shown in Figure 8 and lead to interesting results. Figure 8a shows the CDF's decile values. The maximum value for each year (2005,2009 and 2015) has always been obtained for the WEI class [70,80[, but only in 2015, this class reached values higher than 40%. To fully understand the evolution over time of environmental value inside the L'Alcora municipality using the WEI, Figure 8b,d must be analyzed carefully. These figures show the incremental analysis of the CDF's decile evolution through the analyzed time period, allowing us to obtain specific results for the L'Alcora municipality. Figure 8b shows that between 2005 and 2009, a significant loss of environmental value was observed. The lowest WEI class [10,20[ increased its area by 1.81% due to the loss of area of other low-WEI classes. As shown in Figure 8c, between 2009 and 2015 a significant loss (−2.77%) of a high-WEI class [80,90[ was changed into a lower class [70,80[ at almost exactly the same rate (2.73%). Additionally, a significant loss of the highest WEI class [90,100] (−0.90%) was observed together with the increase (1.35%) of a lower WEI class [60,70[. A very significant result is shown in Figure 8d, which shows the comparison between the CDF's deciles in the 2005−2015 period. The largest difference in the whole analysis (−3.82%) was found in the loss of the highest WEI classes [80,90[+[90,100], which were turned into lower WEI classes, finally leading to a decrease in the average WEI, which has been described above. A very significant result is shown in Figure 8d, which shows the comparison between the CDF's deciles in the 2005−2015 period. The largest difference in the whole analysis (−3.82%) was found in the loss of the highest WEI classes [80,90[+[90,100], which were turned into lower WEI classes, finally leading to a decrease in the average WEI, which has been described above.

Conclusions
A new Weighted Environmental Index (WEI) based on the SIOSE object-oriented model and using GIS data has been introduced for analyzing environmental status through the evolution of land use over time. The versatility of the WEI is based on the fact that the user can define or modify the values of the evaluation factors (F i ) in order to adapt them to the case study under analysis. The methodology is completed by defining each specific land use weighting factor (β jk ), providing great versatility for analyzing land use/land cover change over time.
A demonstration of the application of the WEI to two different case studies (at regional and municipal levels) has been shown. The application of the WEI to these two case studies has demonstrated that the WEI is a powerful tool for analyzing land use change over time and has two major advantages over other environmental indexes. Firstly, the WEI is built based on periodically updated official data, so it avoids subjectivity. Secondly, the WEI can be applied to the analysis of land use change at different scales, and its application allows for the performing of local, regional or even national analyzes and comparisons. Additionally, the WEI is a flexible tool that covers a whole range of situations since it is able to analyze land use evolution over time based on the SIOSE object-oriented model and GIS data.
The application of WEI allows for the performing of detailed statistical analyses, leading to key conclusions about land use changes inside the study area and their environmental implications, quantifying and analyzing trends of environmental quality. The WEI index is based on the definitions of the values of evaluation factors (F i ) and their corresponding weights (α i ). The values of F i and α i should be decided by the modeler on the basis of expert knowledge considering local advice, following the methodology used by previous and simpler environmental indices identified in scientific literature. Strategically selecting and justifying the appropriate values of evaluation factors and weights allows us to use the WEI both for overall and fast screening or for precise and in-depth evaluation purposes [103].