1. Introduction
Land consolidation (LC) is recognized as a process that is important and effective for sustainable development and spatial planning of rural areas but is also expensive and long-lasting [
1].
The importance of land consolidation and its positive effects on agricultural land effective use and rural revitalization are thoroughly researched. Even though the positive effects of land consolidation have been well known for centuries, it is still not provided on a satisfactory level, and their positive effects are not exhausted yet. According to Sayilan [
2], even though the LC works in Turkey were initiated in 1961, sufficient success could not be achieved. According to research by Yin et al. [
3], based on 92 peer-reviewed articles, “rural land consolidation had played more positive than no/negative roles in facilitating rural revitalization, with cases reporting generally positive outcomes accounting for 74%.” The main rural issues recognized, such as land fragmentation, eco-environmental destruction, industrial lag and rural hollowing, could be alleviated by proper utilization of LC [
4]. The paper by Ying et al. [
5] stated that “scant attention has been paid to the relationship between population hollowing and rural sustainability”, and one conclusion stated that the difference of rural hollowing in different villages was significant and that per capita of the village was one of three variables which affect the household directly. Land fragmentation has numerous negative effects on agricultural production, which cause a decrease in yield, an increase in treatment duration, fuel cost and time required to access the parcels, and consequently decreases in agricultural profitability [
6]. Bearing in mind these results, it is obvious that LC could play a significant role in the alleviation of rural hollowing, also containing significant possibilities for rural revitalization.
Starting from the assumption that the LC process is long-lasting and expensive, it could be stated that it is not possible to provide enough financial and other necessary resources (experts, technical and logistic) in order to realize all required LC processes at the same time and on rural areas which need revitalization. This fact causes priority determination and choosing appropriate models for proper decision making. In order to avoid the subjective decision, different optimization methods are utilized in the issue of prioritization of land consolidation [
7,
8,
9,
10,
11].
Since the complexity of land consolidation caused by numerous possible criteria and limited resources, as well as the imperfections of optimization methods, it is questionable if there exists an optimization method that could result in the best solution. The optimization method for decision making could be based on the different approaches and could have similarities and differences [
12], consequently resulting in different ranks of alternatives. In order to avoid the uncertainties in different optimization methods, it is possible to utilize a combination of them to find the optimal solution. Tomić et al. [
13] utilized PROMETHEE I, PROMETHEE II, TOPSIS, WSM-EW and WSM-AHP methods in order to calculate the land consolidation suitability index. Concluding that “the ranking results was more dependent on indicator weights than the choice of multi criteria analysis (MCA)”, they recommend the utilization of the simplest WSM method and that the decision-makers review and approve indicators’ weighting (preferences). This conclusion immediately implicates that the ranking is significantly dependent on the decision-makers’ ability to choose the proper weighting for each utilized method.
In this research, we start from the assumption that chosen weights are burdened with errors and that it is not possible to determine their true values and that, consequently, resulted ranking is also burdened by random errors. This characteristic of optimization methods allows the statistical analysis of obtained results, and it is the main issue of this research. The rank analysis was provided by introducing 8 and 10 criteria, and obtained results were compared.
The analysis of CM ranking in this research is based on finding the average ranking of each municipality and the standard deviation of ranks. The average ranking and the standard deviation define the span in which the rank could belong with the same probability. In the case of limited resources, it is possible to determine the number of CM which could be consolidated. The candidates for land consolidation are ranked better or equal to the number of land consolidations defined in advance.
2. Materials and Methods
The research area covers the territory of the municipality of Bela Crkva, located in the South Banat district of the Autonomous Province of Vojvodina in the Republic of Serbia. The research included the areas unavailable for construction development of fifteen cadastral municipalities with a total area of 33,630 hectares, 68,094 parcels and 10,968 future land consolidation participants. The material used for quantifying the criterion of each alternative was obtained from the database of the Real Estate Cadastre (SKN) of Bela Crkva, the Secretariat for Agriculture of the Municipality of Bela Crkva and the Statistical Office of the Republic of Serbia.
The main criteria originated from the main goals and tasks of LC are ranked from f1 to f8. The additional criteria f9 and f10 were introduced in order to research their influence on cadastral municipalities (CM) ranking related to the main goals of LC. This research was provided firstly by all 10 criteria and after that with the 8 basic criteria.
This database is, in fact, the basis of the Republic Geodetic Institute for the Municipality of Bela Crkva, and it is located in the Bela Crkva Real Estate Cadastre Department (RCD). This database contains information about all properties on the territory of the municipality. From it are taken the data on the areas unavailable for construction development (
Table 1):
(f1) The share of agricultural land in the total area of the land consolidation area;
(f2) Average surface area of cadastral parcels for each land consolidation area;
(f3) Average number of cadastral parcels per participant for each land consolidation area;
(f4) Average surface area of the participants’ property for each land consolidation area.
This database is the local base of the Municipality of Bela Crkva, and it contains data on agricultural land and estates. The data on the future land consolidation projects are taken from this base, i.e., (
Table 2):
The Statistical Office of the Republic of Serbia (RZS) is a special professional organization in the state administration system in the Republic of Serbia. In the system of official statistics of the Republic of Serbia, according to the provisions of the Law on Official Statistics, the RZS performs activities within the system of the Republic of Serbia as the principal participant. The RZS, among other things, carries out activities related to economic production, the implementation of the census, conducting household research, carrying out research relating to the economy and agriculture, and the introduction and management of statistical registers, with the exception of certain financial sector research.
From this database, the data on active agricultural population (criterion f
7) in the future land consolidation areas are collected (
Table 3).
The total costs of land consolidation and cost per hectare (criterion f
8) are explicated in EUR and given in
Table 4.
The integral assessment method is based on the utilization of 6 multi-criteria analysis methods (AHP, VIKOR, COPRAS, ELECTRE, TOPSIS and SAW) and analysis of the differences in obtained rankings for each method. Analytic Hierarchy Process (AHP) is based on the structure of hierarchy in decision problems [
14]. The hierarchy consists of three levels: on the top level is the goal of the decision; at the second level, the criteria are situated and followed by alternatives on the third level. According to the same literature [
14]: “Hierarchical decomposition of complex systems appears to a basic device used by the human mind to cope with diversity.“ Another explication of the AHP reads, “The Analytic Hierarchy Process (AHP) is a theory of measurement through pairwise comparisons and relies on the judgements of experts to derive priority scales” [
15]. AHP method is applied in different land consolidation projects expressing consent for solving problems in the land consolidation reallocation process (where the satisfaction of participants was 91.5% with the solution based on the AHP model compared with 62.7% satisfaction with the solution based on classical interview-based model) [
16] and regional land consolidation [
17]. Although the AHP method is developing still criticism exists, the possibility of the appearance of the Rank Reversal phenomenon, variability of set criteria as well as its construction in accordance with reasons and simplicity led to the statement that “almost all applications of AHP are potentially flawed” [
18].
Similarities and differences between VIKOR (VIKOR is the acronym in the Serbian language which means “Multicriteria Optimization and Compromise Solution”) and TOPSIS (Technique for Order of Preference by Similarity to Ideal Solution) multiple criteria decision-making methods are based on an aggregating function, which represents “closeness to the ideal”, where in the VIKOR method the linear normalization, while in the TOPSIS method vector normalization is used to eliminate the units of criterion function [
12]. The modified VIKOR method was utilized in solving a complex problem in the land consolidation process which included “Production dimension”, “Environmental dimension” and “Life dimension” [
19]. The correctness of the TOPSIS method is also examined for the possibility of increasing the efficiency of identification of land for consolidation, and results showed that “the lowest conformity of results was achieved for a combination of Hellwig’s method and TOPSIS” [
1]. This result indicates that combinations of methods should be chosen very carefully and proven by other methods.
The essential principle of the COPRAS (Complex PRoportional ASsessment of alternatives) method is the possibility of combining values of all indicators, and it uses classical normalization and assumes direct and proportional dependence of priority and utility degree of studying alternatives [
11]. The COPRAS method and its modifications are utilized in different kinds of problems [
20], including land consolidation [
21].
The SAW (Simple Additive Weighting) is the oldest and most widely known practically used method [
22]. The SAW method clearly shows the main concept of multi-criteria evaluation methods integrating criteria values and weights into a single magnitude. The SAW or WSM (Weighted Sum Model) could be applied effectively in many different situations during land consolidation project realization where a lot of alternatives could appear (the exchange of land parcels between participants, public facility allocation, land use suitability analysis, etc.) [
23]. The multi-criteria decision analysis method ELECTRE (ELimination Et Choice Translating REality) “has a good potential to solve multi-objective problems when compensation among criteria is not allowed” [
24]. The ELECTRE method consists of eight steps [
24]: normalization of the decision matrix, multiplying the columns of the decision matrix by the associated weights, determining the concordance and discordance set, calculating the concordance matrix, calculating the discordance matrix, determining the concordance dominance matrix, determining the discordance dominance matrix, determining the aggregate dominance matrix and eliminating the less favorable alternatives.
The development of multi-criteria analysis through the development of multi-criteria decision-making methods showed is twofold: on the one side, it broadens the models involved, and on the other side, it broadens the mathematical models involved. This development is caused by the fact that existing methods (1) are not perfect and (2) should be improved in order to cover the different domains of practice. For example, the AHP methods combined with fuzzy sets results in FAHP [
25], while the modified VIKOR method combined Fermatean hesitant fuzzy sets resulted in novel multi-attribute decision-making approach [
26]. In the paper by Ozdagoglu et al. [
27], the PIPRECIA [
28] method is integrated with COPRAS for solving transportation problems. In the paper by Zahid et al. [
29], the ELECTRE method is combined with fuzzy set theory in order to “empower it to encompass the abstinence part of vague information in addition to satisfaction and dissatisfaction in two dimensional scenarios”. Deng and Chen [
30] modified the TOPSIS method and combined it with intuitionistic fuzzy set (IFS) in order to describe the uncertainty better, while Zhang and Dai [
31] introduced novel TOPSIS method decision-theoretic rough fuzzy sets. All of those multiple-criteria decision-making improvements result with additional expert knowledge needed in order to utilize them properly. In this paper, we propose a simple approach: to use existing models with respect to their imperfections which is explicated by using statistical methods for analysis.
The general summary of the given methods is that they are applicable to land consolidation problems, but no one could be used as an ultimate solution, i.e., every single method is at risk of shortcoming influence on final results.
In this research, we introduce the assumption that every obtained rank by every utilized method is burdened by random error, and it could be explicated in the following way:
where
- -
–right (ideal) rank of -th cadastral municipality
- -
–rank of -th cadastral municipality obtained by -th method and
- -
–random error in ranking of -th cadastral municipality by -th method
In this case, the “ideal” rank of each cadastral municipality could be determined (estimated) as a mean value:
where
denotes the mean value of ranks for
-th cadastral municipality and
the number of utilized methods.
The standard deviation of ranks is calculated as follows:
The interval in which could belong the rank of
-th cadastral municipality than is determined as follows:
where
is the quantile of Student’s distribution,
degrees of freedom and
is level of significance. In this research we adopt
and
. Consequently the value for the quantile of the Student’s distribution is
.
The final interval of ranking
is determined as follows:
4. Discussion
In this research, the multi-criteria analysis approach was utilized in order to determine the rank of prioritization of cadastral municipalities for realizing land consolidation projects. The basic assumption for providing the analysis is that the resources for land consolidations are limited, and consequently, an objective approach is required. However, the imperfections of multi-criteria analysis methods take a real risk that the wrong decisions could be taken. To avoid this possibility or to reduce risk as much as possible in this research, the idea of random errors in ranking results is introduced. Every utilized method is considered a random process, and consequently, the results of ranking by each utilized method are treated as a result burdened by random error. This allowed the statistical analysis, i.e., by defining the interval in which the rank of each cadastral municipality belongs with equal probability. In our study, the probability is adopted at the level of 95% and of 5 degrees of freedom. Those assumptions resulted in the different positions of priority as well as the different possible intervals where each cadastral municipality belongs. To test the sensitivity of the method, the analysis was provided by calculating ranks with ten and eight criteria.
In the case that financial resources were limited to the six cadastral municipalities that could be land consolidated, it is obvious that six cadastral municipalities could be eliminated in both ten and eight criteria cases; four cadastral municipalities can be adopted with probability of 95%, and five municipalities are candidates for further analysis for final decision making. Final decision making could be provided by using additional criteria or by using obtained statistical parameters (lowest average rank, best rank within intervals, or smallest standard deviation, etc.).
The proposed method reduces the subjective or arbitrary decisions in the process of resource distribution, respecting the imperfections of multi-criteria analysis methods.