Next Article in Journal
Applying the Evaluation of Cultural Ecosystem Services in Landscape Architecture Design: Challenges and Opportunities
Previous Article in Journal
Relationship between Physical Environment Satisfaction, Neighborhood Satisfaction, and Quality of Life in Gyeonggi, Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Symmetry of Statistical Methods Applied for the Identification of Agricultural Areas

by
Justyna Wójcik-Leń
1 and
Przemysław Leń
2,*
1
Department of Environment and Agriculture, District Authority Office in Brzozów, 36-200 Brzozów, Poland
2
Department of Environmental Engineering and Geodesy, University of Life Sciences in Lublin, 20-950 Lublin, Poland
*
Author to whom correspondence should be addressed.
Submission received: 13 May 2021 / Revised: 15 June 2021 / Accepted: 19 June 2021 / Published: 23 June 2021
(This article belongs to the Section Land Planning and Landscape Architecture)

Abstract

:
The main priorities of the common agricultural policies of the European Union (EU) are improvement of the quality of life in rural areas for their inhabitants as well as the optimum utilisation of rural resources. The most efficient tools to improve the management conditions and utilise the potential of land are land consolidation works aimed at creating more favourable management conditions in agriculture and forestry through improving the territorial structure of farms, forests and forestland; the reasonable configuration of land, aligning the limits of real properties with the system of irrigation; and drainage facilities, roads and terrain. The development of agriculture in Poland and its production capacity are considerably differentiated in terms of space. At present, Poland has agricultural areas which, in many respects, have a chance of competing with agriculture in the other member states of the European Union. However, in some areas, agricultural production run by private farms owned by individuals is on the verge of falling below the limit of profitability or falls below the limit of profitability. Currently, Poland lacks tools (strategies) allowing identification of land for intensive agricultural production as well as information about agricultural land that should be developed for non-agricultural purposes. Therefore, it is necessary to develop a methodology for identifying similar areas using available tools that can facilitate reliable identification of the areas relating to the indicated factors. Taxonomic methods can be used for clustering purposes. The study materials are data derived from real property register databases referring to one of the districts (poviats) situated in east-central Poland. As a final result, a method of clustering villages according to similar land-use categories was developed. It was created using two independent statistical methods: Ward’s method and the complete-linkage method. The highest consistency was observed in two groups of identified types of areas sharing very similar characteristics. A high index of similarity of both methods—the so-called Rand index—testified to the reliability of the results of calculations. The results of clustering corresponded to a large extent to actual features defining the use of land in the analysed villages as well as the terrain relief.

1. Introduction

The development of agriculture and its production capabilities, when considered globally, are characterised by high variability in many countries of the world. This situation is due to the processes of long-term transformations in politics and agricultural management in areas with different social, demographic and economic situations. Issues related to agricultural crops have been an area of research and interest among many scientists all over the world due to agricultural production capabilities and the need for protecting such land. As shown by studies carried out by many authors, the area of arable land throughout the world is decreasing [1,2,3,4,5,6,7,8,9]. Research reveals that, despite the continuous growth of the world’s population, the supply of land for the needs of the population is limited on all continents. The rapid growth of the world’s population reflects the dynamic development of civilisation [10,11], creating an unprecedented requirement for lands other than crops or forestland [12]. According to the studies by [13,14,15], crop and forest resources have been almost completely exhausted all over the world. Thus far, the use of land and the changes in the use of land have been widely investigated, mainly in the context of their effects on the landscape [16,17,18,19,20,21,22,23,24], soil protection [25,26,27,28], climate change [29,30,31,32,33,34], policies used at a local level aimed at protecting crops and forestland [35,36,37] as well as the land’s comprehensive reconstruction through land consolidation works [38,39,40,41,42,43].
Studies indicate [44,45,46,47] that there are areas where private farms run agricultural production on the verge of or below the level of profitability. It can be seen that such areas are the most severely affected by all changes and system transformations and are often referred to as “problem areas”. A main factor contributing to the formation of problem areas includes the unreasonable utilisation of natural resources, which intensifies erosive degradation and soil acidity as well as the depletion of soil organic matter. Other hazards to the environment and agriculture are the concentration of industrial production, the locations of landfill sites and dust emissions contributing to local pollution of agricultural soils [48]. In addition, these are grounds upon which agricultural production is inefficient and onerous: “(...) agricultural production in respective regions is determined by natural conditions that almost “automatically” lead to underdevelopment of the whole agricultural infrastructure and culture, which results in the region’s backwardness” [49]. Therefore, attention should be paid to areas with limited production potential, lower income per capita and delayed economic development. These areas are increasingly exposed to marginalisation and exclusion from the list of areas with potential for development. Agriculture in such areas is doomed to fail. However, to this end, it is necessary to develop a comprehensive, reliable methodology using statistical tools that will allow identifying areas that are similar in terms of land-use categories because the tools provide a possibility of determining the variety of land covered by the study.
In studies regarding spatial and economic phenomena occurring in larger rural areas, villages are clustered into larger typological units. These areas can be confined within the administrative boundaries of a village or can be larger typological areas referring to, for instance, soil classes. This is associated with both the differentiation and similarity of rural areas. Identification of the most similar areas allows analysing and capturing their spatial differentiation [50]. In Poland [51,52,53,54,55,56,57] and in many other countries of the world [58,59,60,61,62], available methods in many fields of science are used for delimiting areas. Areas with similar land-use categories have been delimited using recognised taxonomic tools and in particular, agglomerative, hierarchical clustering procedures. Hierarchical methods are among the most frequently used taxonomic procedures. They lead to the identification of a full hierarchy of clusters with a monotonically increasing similarity index. The clusters of higher orders include separable clusters of lower orders. There are agglomerative (bottom-up) and divisive (top-down) clustering procedures. The agglomerative methods presume that each unit is initially a separate cluster and then, in a sequential manner, the numbers of existing clusters are reduced by agglomerating them into higher order groups. The proceeding ends when one cluster is obtained that includes all units of the set. Divisive methods use a reverse algorithm. The analysed set of units is initially treated as one group, and the classification gradually increases the number of clusters by dividing the currently existing groups until one-element groups are obtained. An advantage of hierarchical methods is the possibility of presenting the results of classification in a compact graphic form by means of a dendrogram illustrating subsequent links between groups of higher orders. The most popular agglomerative, hierarchical clustering method is Ward’s procedure, determined as a squared Euclidean distance to distances between objects. The idea behind Ward’s method is to minimise intragroup variability (maximise intergroup variability) at each clustering stage [63]. This was the reference method used in this study.
In this paper, another clustering method was used as a control—the complete-linkage approach—in which the underlying element was a city block distance that, as a percentage structure, generally complemented the index of similarity between two categories by 100%. The results of research were presented as a dendrogram illustrating accurately the whole clustering process and the mean output variables in each cluster and, in addition, indicators of mean values for clusters calculated as the ratio between the mean value for a specific cluster and the total mean value. Highly consistent results between two statistical methods can be a guarantee of well-made decisions allowing the reliable identification of picked objects. One method does not provide a sufficient guarantee in this case. In the cluster analysis of objects constituting a study sample, a dendrogram can be deemed an estimator of the hierarchical structure of the whole population only when the analysed objects are described by essential features [64].
The clustering covered 44 villages of the Brzozów district, according to 24 indicators describing land-use categories. The study area was not a random choice because the authors had already carried out studies in this area, covering the dynamics of changes in the use of land in the years 1872–2008 and a preliminary land use analysis in 2010. The calculations were based on unprocessed values of land use indicators: arable land, orchards, permanent meadows, permanent pastures, built-up agricultural land, pond bottoms, ditch bottoms, agricultural land with tree stands and shrubs, wasteland, forests, other land with tree stands and shrubs, housing grounds, industrial grounds, other built-up grounds, land under building development, leisure grounds, surface mining grounds, roads, other transport grounds, grounds for the construction of roads, grounds for the construction of roads, water courses, still waters, ecological areas and other various grounds because this way, higher weights could be naturally assigned to types of land extending over relatively larger areas. The final result of the work was the development of clusters of villages that were similar in terms of land-use categories in the area of 44 villages of the Brzozów district.

2. Materials

The study materials were actual field data derived from real property register databases. There were created using two independent statistical methods: Ward’s method and the complete-linkage method. The use of these methods resulted in a combination with outcome conformity, which constituted sufficient study material.
The study material was a complex of 44 villages forming part of 6 communes in the Brzozów district situated in southeastern Poland in the Subcarpathian region (Figure 1).
The study area covered 44 precincts with an area totalling 53,941.00 ha, divided into 144,610 administrative plots owned by individual farmers. The study made use of 24 indicators describing land-use categories. Table 1 presents the distribution of respective indicators in the whole analysed group.
The rural areas of the Brzozów district were characterised by high variability in their categories of land use a (% forests), b (% arable land), c (% fields), d (% pastures). For instance, the share of arable land ranged from 2.66% to 71.15%, and that of forestland from 6.49% to 96.85%. Most indicators showed very strong right-sided asymmetry, which meant that for most areas, the indicator had low or average values, but for a number of areas, the values were high or very high. For example, the share of permanent meadows was highly asymmetrical (A = 2.54) due to the fact that in spite of the mean level being low (not more than 2.80% in every second commune), the maximum value was 33.94%. The spatial distribution of the four factors (uses of land) with the highest percentage shares within the study area is illustrated in Figure 2.

3. Methods

The research comprised calculations and comparison of the results of two statistical approaches: Ward’s method and complete-linkage using the STATISTICA PLUS program developed by StatSoft Polska. Using Ward’s method (based on data in Table 1), it was possible to identify clusters of rural areas in the Brzozów district, according to land-use categories. The computational algorithm was based on differences (a so-called distance matrix) between the analysed rural areas determined as squared Euclidean distances. The optimum number of clusters was obtained by cutting the arms of the dendrogram at points at which they became longer and the distances between clusters were significantly larger, which made it possible to create five clusters of villages differing in the proposed features.
In complete-linkage clustering, five groups were selected to facilitate comparing them with the results obtained using Ward’s method. However, the division into four groups is more clearly marked on the chart—when the cutting line is moved slightly to the right, groups A and B or six groups are combined—when the cutting line is moved slightly to the left, C is divided. The accurate courses of the clustering process are presented as dendrograms (Figure 3), and Figure 4 illustrates their spatial distribution.
Both methods selected for comparison were slightly different. Ward’s method analysis based on the squared Euclidean distance matrix contributed to increasing the significance of specific types of land (forestland, arable land, meadows etc.) in the results of classification. For complete-linkage, the variation in less-significant types of land could have had a relatively bigger share in the clustering results. Of course, the types of land for which the share was a fraction of a percent were insignificant for both methods used.

4. Results

Table 2 shows the mean values for each land use indicator within each of the five groups. For each type of land, the colour green means a relatively low share of the group, and red corresponds to a higher share. In addition, the value of test probability p evaluating the statistical significance of variations between the identified groups of rural areas in terms of the share of a specific type of land is given.
The description of the group is complemented by Table 3, showing mean group indicators that illustrate the group-to-total-population mean ratio. For instance, the mean share of arable land in the areas from group A was almost identical to the general mean (0.97), for areas from group E it was, on average, four times lower (0.25) and in group D, it was almost two times higher (1.74) than the mean for the whole population.
Complete-linkage clustering is presented in Table 4 and—in terms of spatial distribution—in Figure 3 and Figure 4. Having analysed the mean values for respective types of land in the identified groups, it could be stated that group E featured the highest percentage of forestland and pond bottoms. These were mainly leisure grounds with multiple ponds and waterholes, and most of the land had maintained its natural landscape. Group C showed the decidedly highest share of arable land in comparison to other groups with a relatively small share of forestland. Villages from Group C were typically agricultural areas with good quality soils and wide-scale agricultural management. Group D had a relatively high share of arable land, the highest share of meadows and the lowest of forestland. These were areas with varied terrain relief. Group A featured a high share of pastures and quite a high share of forestland, while most types of land had an average share in group B, except that it clearly had the highest share of the relative area of agricultural land with tree stands and shrubs.
The descriptive analysis of the identified groups can be supplemented by a compilation of group indicators, providing a relatively easy way of characterising the distinctive elements of each group (Table 5).
In order to verify the methods used the index of similarity of both identified divisions—the so-called Rand’s index—was determined. Rand’s index is calculated based on the following concept—all pairs of objects are analysed, and those falling into the same group in both divisions being compared or falling into different groups in both divisions are taken into account. The number of object pairs divided by the number of all possible pairs is Rand’s index value. This leads to the formula:
R = Z / ( n · ( n 1 ) 2 )
where
R is Rand’s index,
N is the number of objects, and
Z is the number of object pairs linked in both classifications.
Surveys showed that its value was 0.67, which meant that the similarity of both identified divisions was slightly above average, as this index could assume values from 0 to 1.

5. Discussion and Conclusions

The studies made it possible to present the results of clustering the rural areas in the Brzozów district with respect to land-use categories, using Ward’s method and the complete-linkage method. In the first method, differences (a so-called distance matrix) between the analysed rural areas were determined as a squared Euclidean distance. The accurate course of the clustering process was presented as a dendrogram (Figure 3). The division into five groups identified as A–E was quite clear (in addition, the division into three groups of A and B and C and D as well as E was very clear), whereas the division using the complete-linkage method led to the identification of five groups that were more similar in terms of count than in Ward’s method. Verification of the methods using an index of similarity of both identified divisions—the so-called Rand’s index—showed that the value of similarity was 0.67 (range 0–1).
Studies carried out based on two independent groupings showed that the areas most similar in terms of land-use categories were E/E and C/D, where the index of similarity—Rand’s index—was 1, with 100% village coverage (E/E) and 0.71% (C/D). The area of the E/E group comprised villages in the highlands situated in the eastern part of the district (Hroszówka, Ulucz, Wołodź). These villages were characterised by the highest share of forestland and pond bottoms. Other uses of land constituted very small percentages. These villages have a historical background because in 1947, during “Operation Vistula”, the villagers were expelled, and now there are few residents in that area. The territory maintained its natural landscape with little signs of human interference. Villages in groups C (Stara Wieś, Buków, Haczów, Jasionów, Trześniów, Wzdów, Jasienica Rosielna) and D (Stara Wieś, Buków, Haczów, Jasionów, Trześniów, Wzdów) were situated in plain and hilly terrain with a high share of arable land and meadows and the highest share of pastures. The identified similar areas featured the lowest share of forestland. The identified area comprised plains featuring very good conditions for agricultural production.
The results of clustering to a large extent correspond to actual features regarding the use of land in the analysed villages as well as terrain relief, as corroborated by previous studies by Leń [65,66] in the study area. The statistical methods used for clustering of similar areas provided satisfactory results in terms of assessment of the actual degree of similarity of the spatial distribution of land-use categories in the study area. High consistency between the results of each combination of pairs of the selected statistical methods could guarantee the right decisions in the development strategy for the specific area. A strategy based on the calculation results derived from a single method is not a sufficient guarantee in this case. The credibility of results is confirmed by the conformity of calculations with the actual needs of the region. The highest consistency was observed between the results of clustering for groups E/E and C/D, which was corroborated by the high Rand’s index of similarity. Considering the purpose of the study, it should be emphasised that two independent statistical methods were used, and the calculated index of similarity in village clustering provides a reliable source of data on clustered areas with very similar features.
It should be noted that funds allocated to the consolidation of land are certainly insufficient in countries like Poland which—due to its historical background—is now facing huge problems in agriculture in connection with the defective spatial structure of rural areas, especially in the eastern, south-eastern, southern and central parts of Poland. The proposed methods allow identifying areas featuring the best quality of agricultural land where research should be continued to identify more specifically the areas where land consolidation is advisable. Such a methodology will be used by the authors in further detailed studies concerning the options of clustering agricultural problem areas to ensure the possibility of their management during land consolidation works.

6. Conclusions

In order to identify rural areas featuring low soil productivity for planning land surveying works in the process of land consolidation, a methodology was developed based on data fully reflected in actual features that was related to the quality of land in the surveyed villages. The results of this study could be helpful in identifying areas for the process of land consolidation aimed at, among other objectives, aligning the chances for development and preservation of the agricultural nature of areas with unfavourable natural and landscape conditions. It should be noted that the designed algorithm could be one of the elements taken into account in creating strategies at the local, district and regional level because it also allows identifying areas with the most favourable soil class, which does not fall under the scope of this study. Familiarity with the spatial variability of soil quality classes within the area of the commune, district and voivodeship has an impact on planning the development of agricultural areas and the effective allocation of money to rural spatial restructuring. This element also provides important information for the process of consolidation works programming at the regional governmental level. An important advantage of the developed method is that it is versatile and can be applied in virtually any conditions.

Author Contributions

Conceptualization, J.W.-L.; methodology, J.W.-L.; software, J.W.-L.; validation, P.L. and J.W.-L., formal analysis, J.W.-L.; investigation, J.W.-L.; resources, J.W.-L.; data curation, J.W.-L.; writing—original draft preparation, J.W.-L.; writing—review and editing, J.W.-L.; visualization, P.L.; supervision, J.W.-L.; project administration, J.W.-L. Both authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gurtin, M.E.; MacCamy, R.C. Non-linear age-dependent population dynamics. Arch. Ration. Mech. Anal. 1974, 54, 281–300. [Google Scholar] [CrossRef]
  2. Guckenheimer, J.; Oster, G.; Ipaktchi, A. The dynamics of density dependent population models. J. Math. Biol. 1977, 4, 101–147. [Google Scholar] [CrossRef]
  3. Ehrlich, P.R.; Ehrlich, A.H. The Population Explosion; Simon and Shuster: New York, NY, USA, 1990. [Google Scholar]
  4. Ranta, E.; Kaitala, V.; Lindström, J.; Linden, H. Synchrony in population dynamics. Proc. R. Soc. Lond. B Biol. Sci. 1995, 262, 113–118. [Google Scholar]
  5. Bjørnstad, O.N.; Ims, R.A.; Lambin, X. Spatial population dynamics: Analyzing patterns and processes of population synchrony. Trends Ecol. Evol. 1999, 14, 427–432. [Google Scholar] [CrossRef]
  6. Gignoux, C.R.; Henn, B.M.; Mountain, J.L. Rapid, global demographic expansions after the origins of agriculture. Proc. Natl. Acad. Sci. USA 2011, 108, 6044–6049. [Google Scholar] [CrossRef] [Green Version]
  7. Linnér, H.; Messing, I. Agricultural land needs protection. Acta Agric. Scand. Sect. B Soil Plant Sci. 2012, 62, 706–710. [Google Scholar] [CrossRef]
  8. Kurowska, K.; Rudnicki, R. Changes in Land Use in Poland–Comparative Study of Period 2002–2010. 2017. Available online: http://conf.rd.asu.lt/index.php/rd/article/view/155 (accessed on 10 May 2021).
  9. Mondal, P.P.; Zhang, Y. Research Progress on Changes in Land Use and Land Cover in the Western Himalayas (India) and Effects on Ecosystem Services. Sustainability 2018, 10, 4504. [Google Scholar] [CrossRef] [Green Version]
  10. Parry, J.A.; Ganaie, S.A.; Bhat, M.S. GIS based land suitability analysis using AHP model for urban services planning in Srinagar and Jammu urban centers of J&K, India. J. Urban Manag. 2018, 7, 46–56. [Google Scholar] [CrossRef]
  11. Zambon, I.; Rontos, K.; Serra, P.; Colantoni, A.; Salvati, L. Population Dynamics in Southern Europe: A Local-Scale Analysis, 1961–2011. Sustainability 2019, 11, 109. [Google Scholar] [CrossRef] [Green Version]
  12. Bradbury, I.; Kirkby, R.; Guanbao, S. Development and environment: The case of rural industrialization and small-town growth in China. Ambio 1996, 25, 204–209. [Google Scholar]
  13. Xiu-Bin, Z.H.A.O.; Yu-Luan, L.I. Forest transition, agricultural land marginalization and ecological restoration. China Popul. Resour. Environ. 2011, 21, 91–95. [Google Scholar]
  14. Munsi, M.; Areendran, G.; Joshi, P.K. Modeling spatio-temporal change patterns of forest cover: A case study from the Himalayan foothills (India). Reg. Environ. Chang. 2012, 12, 619–632. [Google Scholar] [CrossRef]
  15. Kanianska, R.; Kizeková, M.; Nováček, J.; Zeman, M. Land-use and land-cover changes in rural areas during different political systems: A case study of Slovakia from 1782 to 2006. Land Use Policy 2014, 36, 554–566. [Google Scholar] [CrossRef]
  16. Turner, M.G.; Ruscher, C.L. Changes in landscape patterns in Georgia, USA. Landsc. Ecol. 1988, 1, 241–251. [Google Scholar] [CrossRef]
  17. Moreira, F.; Godinho-Ferreira, P.; Rego, F.C.; Bunting, S. Landscape changes and breeding bird assemblages in northwestern Portugal: The role of fire. Landsc. Ecol. 2001, 16, 175–187. [Google Scholar] [CrossRef]
  18. Kozak, J. Forest Cover Change in the Western Carpathians in the Past 180 Years. Mt. Res. Dev. 2003, 23, 369–375. [Google Scholar] [CrossRef] [Green Version]
  19. Junhui, L.; Jixi, G. Changes of land use and landscape pattern in the boundary change areas in farming-pastoral ecotone of northern China. Trans. Chin. Soc. Agric. Eng. 2008, 11. [Google Scholar] [CrossRef]
  20. 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. Geoinf. 2015, 38, 1–14. [Google Scholar] [CrossRef]
  21. Clavero, M.; Villero, D.; Brotons, L. Climate Change or Land Use Dynamics: Do We Know What Climate Change Indicators Indicate? PLoS ONE 2011, 6, e18581. [Google Scholar] [CrossRef] [Green Version]
  22. Parcerisas, L.; Marull, J.; Pino, J.; Tello, E.; Coll, F.; Basnou, C. Land use changes, landscape ecology and their socioeconomic driving forces in the Spanish Mediterranean coast (El Maresme County, 1850–2005). Environ. Sci. Policy 2012, 23, 120–132. [Google Scholar] [CrossRef]
  23. Meneses, B.M.; Reis, E.; Pereira, S.; Vale, M.J.; Reis, R. Understanding Driving Forces and Implications Associated with the Land Use and Land Cover Changes in Portugal. Sustainability 2017, 9, 351. [Google Scholar] [CrossRef] [Green Version]
  24. Janus, J.; Taszakowski, J. Spatial differentiation of indicators presenting selected barriers in the productivity of agricultural areas: A regional approach to setting land consolidation priorities. Ecol. Indic. 2018, 93, 718–729. [Google Scholar] [CrossRef]
  25. Gzyl, J. Soil protection in Central and Eastern Europe. J. Geochem. Explor. 1999, 66, 333–337. [Google Scholar] [CrossRef]
  26. Zuazo, V.H.D.; Pleguezuelo, C.R.R.; Flanagan, D.; Tejero, I.G.; Fernández, J.L.M. Sustainable Land Use and Agricultural Soil; Alternative Farming Systems, Biotechnology, Drought Stress and Ecological Fertilisation; Springer: Berlin/Heidelberg, Germany, 2010; pp. 107–192. [Google Scholar]
  27. Louwagie, G.; Gay, S.H.; Sammeth, F.; Ratinger, T. The potential of European Union policies to address soil degradation in agriculture. Land Degrad. Dev. 2010, 22, 5–17. [Google Scholar] [CrossRef]
  28. Arredondo-Ruiz, F.; García-Montero, L.; Valverde-Asenjo, I.; Menta, C. Soil-Quality Indicators for Forest Management. Quant. Tech. Particip. For. Manag. 2016, 179–240. [Google Scholar] [CrossRef]
  29. Gibbard, S.; Caldeira, K.; Bala, G.; Phillips, T.J.; Wickett, M. Climate effects of global land cover change. Geophys. Res. Lett. 2005, 32. [Google Scholar] [CrossRef] [Green Version]
  30. Bala, G.; Caldeira, K.; Wickett, M.; Phillips, T.J.; Lobell, D.; Delire, C.; Mirin, A. Combined climate and carbon-cycle effects of large-scale deforestation. Proc. Natl. Acad. Sci. USA 2007, 104, 6550–6555. [Google Scholar] [CrossRef] [Green Version]
  31. Paeth, H.; Born, K.; Girmes, R.; Podzun, R.; Jacob, D. Regional Climate Change in Tropical and Northern Africa due to Greenhouse Forcing and Land Use Changes. J. Clim. 2009, 22, 114–132. [Google Scholar] [CrossRef] [Green Version]
  32. Fall, S.; Niyogi, D.; Gluhovsky, A.; Pielke Sr, R.A.; Kalnay, E.; Rochon, G. Impacts of land use land cover on temperature trends over the continental United States: Assessment using the North American Regional Reanalysis. Int. J. Climatol. 2010, 30, 1980–1993. [Google Scholar] [CrossRef] [Green Version]
  33. Pielke Sr, R.A.; Pitman, A.; Niyogi, D.; Mahmood, R.; McAlpine, C.; Hossain, F. Land use/land cover changes and climate: Modeling analysis and observational evidence. Wiley Interdiscip. Rev. Clim. Chang. 2011, 2, 828–850. [Google Scholar] [CrossRef]
  34. Gaitanis, A.; Kalogeropoulos, K.; Detsis, V.; Chalkias, C. Monitoring 60 Years of Land Cover Change in the Marathon Area, Greece. Land 2015, 4, 337–354. [Google Scholar] [CrossRef] [Green Version]
  35. Ożóg, K. The use of unmanned aerial vehicles for the assessment of land boundaries accuracy. J. Water Land Dev. 2020, 45, 94–99. [Google Scholar]
  36. Kurowska, K.; Kryszk, H.; Marks-Bielska, R.; Mika, M.; Leń, P. Conversion of agricultural and forest land to other purposes in the context of land protection: Evidence from Polish experience. Land Use Policy 2020, 95, 104614. [Google Scholar] [CrossRef]
  37. Janus, J.; Markuszewska, I. Land consolidation—A great need to improve effectiveness. A case study from Poland. Land Use Policy 2017, 65, 143–153. [Google Scholar] [CrossRef]
  38. Leń, P. An algorithm for selecting groups of factors for prioritization of land consolidation in rural areas. Comput. Electron. Agric. 2018, 144, 216–221. [Google Scholar] [CrossRef]
  39. Stańczuk-Gałwiaczek, M.; Sobolewska-Mikulska, K.; Ritzema, H.; Van Loon-Steensma, J.M. Integration of water management and land consolidation in rural areas to adapt to climate change: Experiences from Poland and the Netherlands. Land Use Policy 2018, 77, 498–511. [Google Scholar] [CrossRef]
  40. Janus, J.; Markuszewska, I. Forty years later: Assessment of the long-lasting effectiveness of land consolidation projects. Land Use Policy 2019, 83, 22–31. [Google Scholar] [CrossRef]
  41. Wójcik-Leń, J.; Leń, P.; Mika, M.; Kryszk, H.; Kotlarz, P. Studies regarding correct selection of statistical methods for the needs of increasing the efficiency of identification of land for consolidation—A case study in Poland. Land Use Policy 2019, 87, 104064. [Google Scholar] [CrossRef]
  42. Muchová, Z. Assessment of land ownership fragmentation by multiple criteria. Surv. Rev. 2017, 51, 265–272. [Google Scholar] [CrossRef]
  43. Janus, J. A new approach to calculating distances to parcels: A way to increase the accuracy of farm efficiency analyses and the assessment of land consolidation projects. Comput. Electron. Agric. 2020, 175, 105512. [Google Scholar] [CrossRef]
  44. Wójcik-Leń, J.; Sobolewska-Mikulska, K. Specific features of development of selected agricultural problematic areas in the land consolidation process. J. Water Land Dev. 2017, 34, 249–258. [Google Scholar] [CrossRef] [Green Version]
  45. Wójcik-Leń, J.; Sobolewska-Mikulska, K. Issues related to marginal lands with reference to selected agricultural problematic areas. J. Water Land Dev. 2017, 35, 265–273. [Google Scholar] [CrossRef]
  46. Wójcik-Leń, J.; Leń, P.; Sobolewska-Mikulska, K. The proposed algorithm for identifying agricultural problem areas for the needs of their reasonable management under land consolidation works. Comput. Electron. Agric. 2018, 152, 333–339. [Google Scholar] [CrossRef]
  47. Wójcik-Leń, J.; Sobolewska-Mikulska, K.; Sajnóg, N.; Leń, P. The idea of rational management of problematic agricultural areas in the course of land consolidation. Land Use Policy 2018, 78, 36–45. [Google Scholar] [CrossRef]
  48. Jadczyszyn, J. Regional Differentiation of Problem Areas of Agriculture in Poland; Instrukcja Upowszechnieniowa; IUNG-PIB: Puławy, Poland, 2009; 80p, Available online: http://opr.iung.pulawy.pl/publikacje/IU-163.pdf (accessed on 15 April 2021).
  49. Bis, K. Agricultural problem areas with particular consideration of marginal soils. In Problem Areas in Agriculture in Poland with Particular Consideration of the Lublin Region; Jedut, R., Ed.; UMCS-PTG: Lublin, Poland, 1990; pp. 11–17. [Google Scholar]
  50. Noga, K. Topology of Mountain Villages for the Needs of Rural Management Works; Land Surveying, Scientific Papers; University of Agriculture: Kraków, Poland, 1979; Volume 6, pp. 17–29. [Google Scholar]
  51. Liczkowski, J. Delimitation of agricultural intensity regions using the average differences method. Agric. Econ. 1972, 2, 28–37. [Google Scholar]
  52. Fajferek, A. Economic Region and Methods of Regional Analysis; Dissertations Series; ZN WSE: Kraków, Poland, 1964; Volume 6, pp. 44–59. [Google Scholar]
  53. Steczkowski, J. Principles and Methods of Identifying Agricultural Production Regions; PWRiL: Warsaw, Poland, 1966. [Google Scholar]
  54. Domański, S. Typological procedure in economic and agricultural surveys. Geogr. Rev. 1969, XVI Pt 1, 55–68. [Google Scholar]
  55. Król, B. Taxonomic Method of Smallest Differences Applied in Classification of Farms; Scientific Papers of the School of Agriculture in Kraków, Series Economics; University of Agriculture: Punjab, Pakistan, 1969; Part 5. [Google Scholar]
  56. Noga, K. Method of Grouping the Spatial Structure of Mountain Villages. Land Surv. 1988, 10, 221. [Google Scholar]
  57. Stręk, Ż.; Noga, K. Method of Delimiting the Spatial Structure of Villages for the Purposes of Land Consolidation and Exchange. Remote Sens. 2019, 11, 1268. [Google Scholar] [CrossRef] [Green Version]
  58. Shelia, V.; Hoogenboom, G. A new approach to clustering soil profile data using the modified distance matrix. Comput. Electron. Agric. 2020, 176, 105631. [Google Scholar] [CrossRef]
  59. Abu Bakar, M.A.; Ariff, N.M.; Jemain, A.A.; Nadzir, M.S.M. Cluster Analysis of Hourly Rainfalls Using Storm Indices in Peninsular Malaysia. J. Hydrol. Eng. 2020, 25, 05020011. [Google Scholar] [CrossRef]
  60. David, J.; De Pessemier, T.; Dekoninck, L.; De Coensel, B.; Joseph, W.; Botteldooren, D.; Martens, L. Detection of road pavement quality using statistical clustering methods. J. Intell. Inf. Syst. 2019, 54, 483–499. [Google Scholar] [CrossRef]
  61. Ward, J.H. Hierarchical grouping to optimize an objective function. J. Am. Stat. Assoc. 1963, 58, 236–244. [Google Scholar] [CrossRef]
  62. Marek, T. Analiza Skupień w Badaniach Empirycznych: Metody SAH [Cluster Analysis in Empirical Research SAHN Methods]; PWN: Warszawa, Poland, 1989; p. 121. [Google Scholar]
  63. Kurowska, K.; Kryszk, H.; Marks-Bielska, R.; Kietlinska, E. Spatial Analysis of Afforestation in Poland under Rural Development Programme 2007–2013. 2014. Available online: https://llufb.llu.lv/conference/Research-for-Rural-Development/2014/LatviaResearchRuralDevel20th_volume2-14-21.pdf (accessed on 7 April 2021).
  64. Kurowska, K.; Kryszk, H. Profitability of the farmland afforestation within the Rural Development Programme. Sylwan 2017, 161, 1035–1045. [Google Scholar]
  65. Leń, P. Kierunki Zmian W Strukturze Użytków Gruntowych W Latach 1872–2008 w Powiecie Brzozów [Areas of Change in the Pattern of Ground Water Use in the Period 1872–2008 in Rural District of Brzozów]. 2009. Infrastructure and Ecology of Rural Areas, Nr 4/2009. pp. 167–175. Available online: http://infraeco.pl/pl/art/a_15494.htm?plik=577 (accessed on 7 April 2021).
  66. Leń, P. Podział Przestrzeni Rolniczej Powiatu Brzozowskiego Pod Względem Wartości Produkcyjnej Gruntów Ornych Oraz Użytków Zielonych [Break Down of County Agricultural Space Brzozowski in Terms of Production Value of Cropland and Grassland]. 2010. Infrastructure and Ecology of Rural Areas, Nr 12/2010. pp. 37–44. Available online: http://infraeco.pl/pl/art/a_16166.htm?plik=928 (accessed on 7 April 2021).
Figure 1. Map of location of the analysed area, in the region and in Europe.
Figure 1. Map of location of the analysed area, in the region and in Europe.
Land 10 00664 g001
Figure 2. (ad) The spatial distribution of the percentage shares of four uses of land.
Figure 2. (ad) The spatial distribution of the percentage shares of four uses of land.
Land 10 00664 g002
Figure 3. Division of the study area using two statistical methods.
Figure 3. Division of the study area using two statistical methods.
Land 10 00664 g003
Figure 4. Division of the study area by means of the methods used.
Figure 4. Division of the study area by means of the methods used.
Land 10 00664 g004
Table 1. The distribution of land-use categories in the Brzozów district.
Table 1. The distribution of land-use categories in the Brzozów district.
No.Indicators x ¯ MeSminmaxA
1arable land38.1237.6414.82.6671.150.15
2orchards0.140.030.2301.082.42
3fields5.392.87.170.0233.942.54
4pastures10.911.095.220.129.430.79
5built-up agricultural land2.582.341.2406.290.5
6pond bottoms0.2900.7702.752.56
7ditch bottoms0.150.040.2701.12.16
8agricultural land with tree stands and shrubs3.782.44.430.0718.321.75
9wasteland0.160.110.1600.781.79
10forests32.8930.4518.216.4996.851.52
11other land with tree stands and shrubs0.0600.1500.783.62
12housing grounds0.760.38105.062.45
13industrial grounds0.070.010.1901.215.39
14other built-up grounds0.250.20.180.010.640.69
15land under building development0.030.020.0400.171.59
16leisure grounds0.110.050.2301.394.55
17surface mining grounds0.0400.1801.115.56
18roads2.822.871.010.225.420.22
19other transport grounds00000.025.79
20grounds for the construction of roads000.0100.043.99
21water courses1.340.641.750.019.12.78
22still waters0.0700.301.925.86
23ecological areas0.0100.0900.596.55
24various grounds0.0300.1500.955.62
x ¯ mean, Me—median, s—standard deviation, min.—minimum value, max.—maximum value, A—asymmetry.
Table 2. Mean values in village clusters identified using Ward’s method.
Table 2. Mean values in village clusters identified using Ward’s method.
Type of LandMean Values in Groupsp
ABCDE
arable land37.1522.7242.3766.499.480.0000 ***
orchards0.170.010.040.250.110.1595
permanent meadows3.801.5317.592.020.970.0011 **
permanent pastures11.6115.0111.827.373.290.0054 **
built-up agricultural land2.652.383.212.950.170.0213 *
pond bottoms0.290.690.030.000.830.4231
ditch bottoms0.110.040.540.010.020.0024 **
agricultural land with tree stands and shrubs5.392.542.840.220.180.0017 **
wasteland0.160.140.170.220.000.0726
forests32.7050.7814.2215.4183.310.0000 ***
other land with tree stands and shrubs0.070.020.070.000.000.4593
housing grounds0.820.201.580.250.000.0004 ***
industrial grounds0.080.040.080.020.010.2120
other built-up grounds0.270.120.360.210.060.0284 *
land under building development0.040.000.070.020.000.0074 **
leisure grounds0.080.070.370.030.000.0630
surface mining grounds0.050.100.010.000.000.3456
roads3.021.963.193.330.580.0036 **
other transport grounds0.000.000.000.000.000.7989
grounds for the construction of roads0.000.000.010.010.000.5901
water courses1.401.641.231.170.970.5501
still waters0.110.010.030.010.000.2356
ecological areas0.030.000.000.000.000.8168
various grounds0.000.000.190.000.020.0293 *
Source: Own elaboration. p: value of test probability calculated using the Kruskal–Wallis test (statistically significant differences, significance levels 0.05, 0.01 and 0.001 marked with *, ** and ***).
Table 3. Values of indicators for identified clusters.
Table 3. Values of indicators for identified clusters.
Type of LandIndicators of Mean Values for Clusters
ABCDE
arable land0.970.601.111.740.25
orchards1.220.060.271.800.77
permanent meadows0.700.283.260.370.18
permanent pastures1.061.381.080.680.30
built-up agricultural land1.030.921.241.140.07
pond bottoms1.012.380.100.002.85
ditch bottoms0.710.293.480.090.13
agricultural land with tree stands and shrubs1.420.670.750.060.05
wasteland1.040.851.081.420.03
forests0.991.540.430.472.53
other land with tree stands and shrubs1.330.301.300.090.00
housing grounds1.070.262.080.320.00
industrial grounds1.250.601.230.290.12
other built-up grounds1.090.481.420.830.26
land under building development1.080.051.990.550.08
leisure grounds0.690.593.290.270.00
surface mining grounds1.282.490.280.000.00
roads1.070.701.131.180.21
other transport grounds1.370.001.390.000.00
grounds for the construction of roads0.490.002.792.460.00
water courses1.051.230.920.880.72
still waters1.610.090.420.100.00
ecological areas1.760.000.000.000.00
various grounds0.110.005.660.000.53
Source: own elaboration.
Table 4. Mean values in village clusters identified using the complete-linkage method.
Table 4. Mean values in village clusters identified using the complete-linkage method.
Type of LandMean Values in Clustersp
ABCDE
arable land29.8933.7451.1641.149.480.0000 ***
orchards0.080.080.260.040.110.2726
permanent meadows2.793.993.7219.500.970.0015 **
permanent pastures14.845.559.6911.653.290.0002 ***
built-up agricultural land2.582.022.993.080.170.0039 **
pond bottoms0.340.670.140.030.830.8867
ditch bottoms0.080.120.090.610.020.0091 **
agricultural land with tree stands and shrubs4.5914.691.253.050.180.0000 ***
wasteland0.100.270.220.150.000.0236 *
forests39.5531.5724.7113.7183.310.0000 ***
other land with tree stands and shrubs0.060.000.070.080.000.3226
housing grounds0.550.700.781.670.000.0031 **
industrial grounds0.030.020.120.050.010.4700
other built-up grounds0.180.480.270.310.060.0135 *
land under building development0.020.030.040.070.000.0155 *
leisure grounds0.090.100.070.340.000.0959
surface mining grounds0.030.010.080.000.000.5921
roads2.892.963.082.960.580.0629
other transport grounds0.000.000.000.000.000.0462 *
grounds for the construction of roads0.000.000.000.000.000.8273
water courses1.292.861.081.310.970.4555
still waters0.010.150.130.010.000.3341
ecological areas0.000.000.040.000.000.9304
various grounds0.010.000.000.220.020.1246
Source: Own elaboration. p: value of test probability calculated using the Kruskal–Wallis test (statistically significant differences, significance levels 0.05, 0.01 and 0.001 marked with *, ** and ***).
Table 5. Values of indicators for identified clusters.
Table 5. Values of indicators for identified clusters.
Type of LandIndicators of Mean Values for Clusters
ABCDE
arable land0.780.891.341.080.25
orchards0.580.531.820.280.77
permanent meadows0.520.740.693.620.18
permanent pastures1.360.510.891.070.30
built-up agricultural land1.000.781.161.190.07
pond bottoms1.192.320.490.102.85
ditch bottoms0.490.750.603.950.13
agricultural land with tree stands and shrubs1.213.890.330.810.05
wasteland0.641.691.360.970.03
forests1.200.960.750.422.53
other land with tree stands and shrubs1.020.001.261.430.00
housing grounds0.720.911.022.190.00
industrial grounds0.520.361.850.810.12
other built-up grounds0.711.931.091.240.26
land under building development0.560.791.212.130.08
leisure grounds0.800.890.633.060.00
surface mining grounds0.690.262.040.000.00
roads1.031.051.091.050.21
other transport grounds0.000.992.500.000.00
grounds for the construction of roads0.800.001.920.220.00
water courses0.972.130.810.980.72
still waters0.122.272.000.190.00
ecological areas0.240.002.520.000.00
various grounds0.170.020.026.590.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wójcik-Leń, J.; Leń, P. Evaluation of the Symmetry of Statistical Methods Applied for the Identification of Agricultural Areas. Land 2021, 10, 664. https://0-doi-org.brum.beds.ac.uk/10.3390/land10070664

AMA Style

Wójcik-Leń J, Leń P. Evaluation of the Symmetry of Statistical Methods Applied for the Identification of Agricultural Areas. Land. 2021; 10(7):664. https://0-doi-org.brum.beds.ac.uk/10.3390/land10070664

Chicago/Turabian Style

Wójcik-Leń, Justyna, and Przemysław Leń. 2021. "Evaluation of the Symmetry of Statistical Methods Applied for the Identification of Agricultural Areas" Land 10, no. 7: 664. https://0-doi-org.brum.beds.ac.uk/10.3390/land10070664

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