Increase in the Value of Agricultural Parcels—Modelling and Simulation of the Effects of Land Consolidation Project

: In the theory and practice of valuation, it is commonly accepted that the key feature determining the value of agricultural land is its location, both general and in a speciﬁc part (zone) of a village. The model approach used in the present study can provide the answer to the question of how to maximize the value of agricultural land as part of a conducted arrangement, agricultural works. The study used data on the market sale of agricultural parcels in 10 Polish municipalities. Each parcel was described using a set of features (parameters) that were key to its value and entered into a database. Using the database, two statistical models were built: a multiple regression analysis model (MRA) and an artiﬁcial neural network model (ANN). The studies conducted have shown that changes in such features as surface area, shape, and access to a public road were accompanied by signiﬁcant changes in the market values of parcels. Another important observation was that potential decreases in the value of agricultural parcels as a result of changes in their surface areas were offset (where it was reasonable) by the elimination of their excessive elongation and providing them with an access to a public road. Based on the ﬁndings, it has been concluded that change in land value should be considered one of the effects of executed land consolidation projects (LCP), during which the parameters of agricultural parcels are subject to the biggest changes.


Introduction
Areas used for agricultural purposes play a key socio-economic and cultural role in the development of virtually any country. Agricultural areas are also essential environmentally and ecologically [1]. According to Eurostat data in the European Union, utilized agricultural areas cover 173 million ha. In terms of the size of agricultural areas, Poland comes fifth, behind France, Spain, Great Britain, and Germany [2].
Of particular importance for many of the above-mentioned functions of rural areas is the layout of property borders, often equated with division into parcels. From the perspective of the operation of farms, the layout of parcels and their number, surface area, and shape largely determine the profitability of agricultural activity [3,4], although the number of factors impacting the profitability is naturally much bigger. Excessive land fragmentation increases the costs of farms' operation by decreasing yields and increasing labor and transportation costs [5][6][7][8][9]. In the long run, these factors have a negative impact on the incomes and living conditions of the population, increasing negative migration trends from rural areas to cities [10] and leading over time to an increasing phenomenon of permanent abandonment of land cultivation [11][12][13]. Both these phenomena are not only negative but also hard to reverse. Land acquires a real value when it is sold or used for securing a credit. Estimation of the value of land related to its productivity or market value is also necessary to determine the amount of numerous taxes. It is also needed when securing mortgage credits, exchanging farmland, or determining compensations for the acquisition of farmland to carry out public investment projects. What affects the value and market prices of parcels located in rural areas?
Parcels used for agricultural purposes are sold on the market relatively less often than parcels located in investment areas, in particular those designated for building. Theory suggests [14][15][16][17][18] that when determining the value of any property or its parts, the specificity of the subject of valuation should be taken into account in each case. In the case of agricultural parcels, one should certainly focus on those characteristics of land that are significant for agricultural production (impacting yield and income levels) and those factors that have a significant impact on the value of land as viewed by the owners. It should be stressed that much of the trade in agricultural land takes place between neighbors. It takes place between land owners who live in one locality and similarly perceive the local conditions affecting the assessment of the land's agricultural suitability in a given village.
The most important characteristics determining the value of agricultural parcels include ( Figure 1) location in a specific part (zone) of a village, proximity of other agricultural areas and buildings, soil usability (soil classification), variety and spatial layout of agricultural areas, surface area and shape of a parcel, topography, accessibility, as well as availability of structures and equipment used for agricultural production (land improvements or irrigation system). Among other significant factors are obstacles to cultivation (e.g., flintiness, permanent elements of technical infrastructure, etc.) and agricultural condition [4,6,[16][17][18][19].
Agriculture 2021, 11, x FOR PEER REVIEW 2 of 15 of permanent abandonment of land cultivation [11][12][13]. Both these phenomena are not only negative but also hard to reverse. Land acquires a real value when it is sold or used for securing a credit. Estimation of the value of land related to its productivity or market value is also necessary to determine the amount of numerous taxes. It is also needed when securing mortgage credits, exchanging farmland, or determining compensations for the acquisition of farmland to carry out public investment projects. What affects the value and market prices of parcels located in rural areas? Parcels used for agricultural purposes are sold on the market relatively less often than parcels located in investment areas, in particular those designated for building. Theory suggests [14][15][16][17][18] that when determining the value of any property or its parts, the specificity of the subject of valuation should be taken into account in each case. In the case of agricultural parcels, one should certainly focus on those characteristics of land that are significant for agricultural production (impacting yield and income levels) and those factors that have a significant impact on the value of land as viewed by the owners. It should be stressed that much of the trade in agricultural land takes place between neighbors. It takes place between land owners who live in one locality and similarly perceive the local conditions affecting the assessment of the land's agricultural suitability in a given village.
The most important characteristics determining the value of agricultural parcels include ( Figure 1) location in a specific part (zone) of a village, proximity of other agricultural areas and buildings, soil usability (soil classification), variety and spatial layout of agricultural areas, surface area and shape of a parcel, topography, accessibility, as well as availability of structures and equipment used for agricultural production (land improvements or irrigation system). Among other significant factors are obstacles to cultivation (e.g., flintiness, permanent elements of technical infrastructure, etc.) and agricultural condition [4,6,[16][17][18][19]. Some of the above-listed factors (surface area, shape, and location) are also key in determining the values of land fragmentation parameters (LF). Thus, we can posit the statement that LF parameters impact the value of land or how it is assessed by current or potential land owners. If we accept the relationship between parcel fragmentation parameters and land value as obvious, we can pose the question: How much will the Some of the above-listed factors (surface area, shape, and location) are also key in determining the values of land fragmentation parameters (LF). Thus, we can posit the statement that LF parameters impact the value of land or how it is assessed by current or potential land owners. If we accept the relationship between parcel fragmentation parameters and land value as obvious, we can pose the question: How much will the changes in these parameters affect changes in the value of land? This way of addressing the problem is in line with the findings of numerous existing studies into the impact of land fragmentation on production profitability in the different branches of agriculture or general conditions of agriculture [20][21][22][23].

Materials and Methods
In order to run a simulation of the impact of changes in selected parameters of parcels on the value of agricultural land, data on the market sale of agricultural properties in 10 municipalities located in six Polish regions ( Figure 2) were used: Little Poland (Szczurowa and Jordanów municipalities), Lublin (Łopiennik Górny municipality), Silesia (Lelów, Pilica and Mykanów municipalities), Lower Silesia (Mściwojów and Gaworzyce municipalities), Mazovian (Strzegowo municipality), and Opole (Opole municipality). changes in these parameters affect changes in the value of land? This way of addressing the problem is in line with the findings of numerous existing studies into the impact of land fragmentation on production profitability in the different branches of agriculture or general conditions of agriculture [20][21][22][23].

Materials and Methods
In order to run a simulation of the impact of changes in selected parameters of parcels on the value of agricultural land, data on the market sale of agricultural properties in 10 municipalities located in six Polish regions ( Figure 2) were used: Little Poland (Szczurowa and Jordanów municipalities), Lublin (Łopiennik Górny municipality), Silesia (Lelów, Pilica and Mykanów municipalities), Lower Silesia (Mściwojów and Gaworzyce municipalities), Mazovian (Strzegowo municipality), and Opole (Opole municipality). Based on the information obtained from relevant county offices maintaining registers of property prices and values, unbuilt agricultural parcels sold on the market over the period 2018-2019 (292 transactions in total) were selected. Next, a database was created with each transaction described using information about the surface area and structure of the agricultural areas, the price, date, and parties involved in the transaction, as well as the location (voivodeship, tax district, municipality, cadastral region, and parcel number). Upon verification of the initial list of transactions using maps and aerial photographs, built-up parcels and parcels located within residential areas or in areas designated for building were rejected. Only those parcels that were used for agricultural purposes as confirmed by aerial photographs were left for further analysis (270 transactions). Each transaction was described with an additional set of characteristics derived from cadastral maps and orthophotomaps. Thus, the initial data sets were enriched with geometry information: length, width, and perimeter of each sold parcel. Moreover, the parcels were assessed in terms of their shape, detailed location, vicinity, access to a public road, obstacles to usage, level of agricultural condition, and distance from the nearest built-up area. These are the aspects that based on both literature review and the authors' experience [6,19] could be considered as having a potentially significant impact on prices and values.
Two models were used to assess the impact of market characteristics on the value of parcels: the multiple regression analysis model (MRA) and the artificial neural network model (ANN). The MRA model captures the relationships between parcel characteristics and their market value. Unfortunately, due to its specificity, the model does not capture Based on the information obtained from relevant county offices maintaining registers of property prices and values, unbuilt agricultural parcels sold on the market over the period 2018-2019 (292 transactions in total) were selected. Next, a database was created with each transaction described using information about the surface area and structure of the agricultural areas, the price, date, and parties involved in the transaction, as well as the location (voivodeship, tax district, municipality, cadastral region, and parcel number). Upon verification of the initial list of transactions using maps and aerial photographs, builtup parcels and parcels located within residential areas or in areas designated for building were rejected. Only those parcels that were used for agricultural purposes as confirmed by aerial photographs were left for further analysis (270 transactions). Each transaction was described with an additional set of characteristics derived from cadastral maps and orthophotomaps. Thus, the initial data sets were enriched with geometry information: length, width, and perimeter of each sold parcel. Moreover, the parcels were assessed in terms of their shape, detailed location, vicinity, access to a public road, obstacles to usage, level of agricultural condition, and distance from the nearest built-up area. These are the aspects that based on both literature review and the authors' experience [6,19] could be considered as having a potentially significant impact on prices and values.
Two models were used to assess the impact of market characteristics on the value of parcels: the multiple regression analysis model (MRA) and the artificial neural network model (ANN). The MRA model captures the relationships between parcel characteristics and their market value. Unfortunately, due to its specificity, the model does not capture non-linear relationships, which also appear between the analyzed feature and the variables that characterize it. In order to capture such relationships, ANNs were used to complement the results of regression. ANNs are a tool that can be used wherever non-linear and non-monotonic relationships occur [24,25].
The list of predictors was established based on the theoretical premises of agricultural land valuation and initial results of market analyses conducted by filtering observations by variants of the characteristics, generating plots, and analysis of variance. Initially, these were the 10 features of the property: tax region, general location, detailed location, soil quality, share of arable land, obstacles to usage, vicinity, surface area, parcel elongation, and access to a public road ( Table 1). The authors followed the principle of modeling simplicity [26], the generalized form of which is known under the name of Ockham's razor: in explaining phenomena one should seek simplicity, choosing those explanations that are based on the least number of assumptions and concepts. With reference to statistical models, this principle can be understood as follows: if several different models explain a phenomenon in question to a similar extent, one should choose the model that is conceptually simpler. In this context, in the case of regression, preference was given to a linear model, with the necessary minimum number of predictors ensured by a stepwise procedure. The multiple regression analysis model (MRA) is conceptually simpler and easier to interpret than non-linear regression models. Numerous references to the practical application of such a model in real estate market analyses can be found in Polish and international literature [27][28][29][30][31][32][33][34][35][36][37][38][39][40][41][42][43]. In recent decades, such models were used to determine real estate value, and they provided satisfactory results. There are known MRA models that made it possible, thanks to the stepwise procedure, to find out which factors impact agricultural properties and to what extent [35,44].
In the case of the neural network, preference was given to its simplest architectures, with the necessary minimum number of predictors ensured by analysis of network sensitivity. At the same time, the authors, in accordance with the aim of the study, bore in mind the need to take into account those market characteristics that undergo significant changes as a result of the implementation of a land consolidation project.
Based on initial analyses, it was found that constructing a useful regression model for the full set of observations (270 transactions) was impossible. As far as surface area was concerned, it showed a non-monotonic impact on the market value of parcels: the smallest and the largest parcels were the most expensive. Consequently, this feature was not taken into consideration in the linear regression model. In addition, in the stepwise procedure, the following factors were eliminated from the MRA model: tax region and soil quality. Despite numerous attempts and deletion of outlier observations, these factors remained insignificant for MRA.
Apart from the verification of the individual predictors, it was also necessary to conduct an analysis of residuals and to reject a certain percentage of observations for which the models were unable to predict the values correctly. The analysis of residuals was conducted as part of the construction of the MRA model. The set of observations selected as part of this analysis was used to determine the forms of both the models.
The results of the analyses conducted as part of the construction of the MRA model were used in the ANN model. The latter provided more details to and supplemented the mathematical relationships presented by the regression equation. Surface area was expressed by means of five classes based on the assumption that change in surface area by one unit does not cause such marked changes in prices as change in surface area by a class (e.g., transition from the category of mid-sized parcels to the category of large parcels). Apart from surface area, the model also takes into account soil quality and location in a tax region.
The MRA and ANN models were developed using Statistica 13 with Data Mining module. With regards to the MRA model, the option of multiple forward selection stepwise regression was used, whereas with regards to the ANN model there is the option of ANNs.  Obstacles to usage A discrete quantitative feature assuming three states: 1-large (two or more of the following aspects appear at the same time: flintiness, wetness of the area, visible soil variability, bushes and coppices, wasteland, proximity of forests, sharp angles of parcel borders) 2-moderate (one of the above) 3-lack Vicinity A discrete quantitative feature assuming two states: 0-not directly bordering a built-up area 1-directly bordering a built-up area Surface area * A qualitative feature assuming five variants: up to 0.5 ha (very small parcels) from 0.5 to 1 ha (small parcels) from 1 to 2 ha (mid-sized parcels) from 2 to 3 ha (large parcels) above 3 ha (very large parcels) Parcel elongation A continuous quantitative feature assuming values in the range from 1.0 to 114.6, calculated as a quotient of the parcel's sides (longer/shorter) Access to a public road A discrete quantitative feature assuming two states: 0-no access to a road 1-access to a road * Features used in the ANN model; in the MRA model they were insignificant. ** Wb was calculated as the quotient of the area, which is the basis for calculating the agricultural tax (conversion hectare) and the physical area (physical hectare).

Linear Regression
The linear mathematical model of multiple regression specifies the linear relationship between the dependent variable (Y) and the independent variables X 1 , X 2 , and X i as outlined in the formula below: where: Y-dependent variable β 0 -intercept β 1 , β 2 , β k -independent variables' coefficients X 1 , X 2 , X i -independent variables ξ-random component The parameters of the regression model were estimated by using the least-squares method.
In accordance with the methodological assumptions, the method of forward selection stepwise regression was used, which enables inclusion into the model of only those predictors that explain the dependent variable in a statistically significant way. The analysis of residuals was conducted using the model of multiple linear regression. Deletion of outlier observations brought positive results in subsequent regression models: coefficient of determination R2 changed as follows:  (Table 2). Statistical verification: global F-test and Student's t-tests for each used explanatory variable (Table 2) also yielded a positive result. The model also respected the lower limit of the acceptable number of observations for regression. This number should not be less than 3 × (number of predictors + 1) [45]. The model predicted price volatility with 68% accuracy.

Artificial Neural Networks
Neural networks (ANN) represent an analytical technique modeled on the cognitive system and neurological functions of the brain. First references to neural networks date back to the 1950s, when an attempt was made to present a computing model of the internal processes of the human brain [46]. However, due to limited computing power and insufficient data, the potential of neural networks has not been fully exploited until recent decades [47]. Currently, neural networks are among the most popular tools of artificial intelligence [48]. The neural network used in the present study acquired the ability to predict the values of agricultural parcels based on the predictors fed at its input, after a training process was performed based on gathered market data (203 observations selected as part of the analysis of residuals). The size and structure of the network was determined by the complexity of the phenomenon under study. Although network layers can theoretically consist of any number of neurons, the principle of modeling simplicity had to be followed here as well. A more complex network analyzes more input variables, but it also requires more observations [49]. This is because, generally, each used predictor corresponds to one input neuron.
The connections between a network's neurons are assigned randomly selected weights when the training process is started. As successive examples are presented to the network, these weights are iteratively adjusted so that the differences between the network's predictions and the actual values could be the smallest [26,50].
For the purpose of constructing the network, the data were divided into three sub-sets in the following proportions: -training set: 70% -test set: 15% -validation set: 15% The division of data into subsets reflects how a neural network works. A training set is used to modify randomly selected weights of the connections between neurons. A test set is designed to monitor the training process on a continuous basis (assessment of the network's approximation capability). A validation set is used to assess the network after the training process is completed (generalization capability).
Of several dozens of networks tested concurrently over 200 epochs, the best results were achieved by multilayer perceptron (MLP). This network model is among the most commonly used [25]. It consists of interconnected units, called perceptrons or nodes, which are arranged in layers (Figure 3). The first layer is input data (or independent variables), the last layer is the output layer, and the layers in between are called hidden layers [51].
The model had one hidden layer numbering 11 neurons. The network was trained using the Broyden-Fletcher-Goldfarb-Shanno optimization algorithm (BFGS) as the criterion for seeking the minimum of the error function and the sum of squares (SoS) as its preferred function. The Sum of Squares (SoS) function typically has the formula: where: For hidden neurons, the activation function was the hyperbolic tangent, whereas for the output neuron, it was the linear function. These were the recommended settings for multilayer perceptrons solving regression problems using the sum of squares as an error function.

Results
The next stages of the computation process determined the values predicted by the mathematical form of the multiple regression analysis model (MRA) in different configurations of parcel elongations and access to a road. Then, predictions of the artificial neural network model (ANN) were computed by entering at its inputs data in different configurations of parcel elongations, classes of parcel surface areas, and road access. The results of the studies were presented separately for the model of multiple linear regression and the neural network model. The model had one hidden layer numbering 11 neurons. The network was trained using the Broyden-Fletcher-Goldfarb-Shanno optimization algorithm (BFGS) as the criterion for seeking the minimum of the error function and the sum of squares (SoS) as its preferred function. The Sum of Squares (SoS) function typically has the formula: where:

Yi-predicted values, ti-actual values, n-number of observation
For hidden neurons, the activation function was the hyperbolic tangent, whereas for the output neuron, it was the linear function. These were the recommended settings for multilayer perceptrons solving regression problems using the sum of squares as an error function.

Results
The next stages of the computation process determined the values predicted by the mathematical form of the multiple regression analysis model (MRA) in different configurations of parcel elongations and access to a road. Then, predictions of the artificial neural network model (ANN) were computed by entering at its inputs data in different configurations of parcel elongations, classes of parcel surface areas, and road access. The

Multiple Linear Regression Model
The regression model constructed in accordance with the above-described methodology showed a satisfactory explanatory capacity of prices. The coefficient of determination R2 indicated that residual variation (resulting from the impact of unmodeled factors) was 32% of the total variation of the observed features. The model passed a global test, as evidenced by the value of the Fisher-Snedecor statistic (F = 61). The coefficients for explanatory variables were estimated with low standard errors and passed Student t-tests. Moreover, the MRA model, based on seven selected predictors, passed a substantive verifi-cation: regression coefficient signs indicated relationships consistent with the theory (Table 2).
Based on the data from Table 2, the mathematical form of the multiple linear regression model explaining how the values of agricultural parcels were established can be written as follows (Formula (1)): Y = −20, 811 + 9251 * X 1 + 6160 * X 2 − 56 * X 3 + 3413 * X 4 + 3590 * X 5 + 4948 * X 6 + 2497 * X 7 With a parcel elongation (feature X 3 ), there was a negative regression coefficient. This confirmed an important fact from the perspective of consolidation assessment, namely: a decrease in excessive elongation of agricultural parcels led to an increase in their market value. The other features, when assuming higher values, contributed to an increase in the modeled market prices. Based on standardized regression coefficients (Table 2, Column 2), it was found out that location factors and the share of arable land were most significant in  Table 1, the multiple regression model did not take into account location in a tax region, a parcel's surface area (an important feature from the perspective of the conducted analyses), and soil quality. In the case of linear regression, these features were statistically insignificant even after analysis of residuals was conducted and outlier observations were removed.

Artificial Neural Network Model (ANN)
Another approach to estimating changes in the values of parcels resulting from the implementation of consolidation projects was the artificial neural network model. In the case of this method, the role of the individual predictors can be assessed based on a sensitivity analysis [52]. This is a network's functionality that enables evaluation of the contribution of the individual predictors to explain how the prices of parcels were established (Table 3). Of particular importance were a parcel's access to a public road and location factors (general location-in a given municipality, in the village, and in a given tax region). Other important factors included soil quality, obstacles to usage, and a parcel's vicinity. A less important role as compared to the other predictors was played, according to the network, by surface area, share of arable land, and a parcel elongation. For the neural network model, the correlation levels between the actual and predicted values were as follows: 0.89-for the training set, 0.82-for the test set, and 0.80-for the validation set. The principle is that the result for the validation set should not deviate significantly from the result for the test set. Then, it can be acknowledged that the network generalizes the data structure well and will work well with completely new data.

Simulation of the Change in the Values of Agricultural Parcels with a Change in Their Elongation, Surface Area, and Access to a Public Road
Once the final forms of both the models were determined and verified, a simulation was run of the changes in the value of 1 ha of agricultural land as a result of adjusting the parcel's surface area, elongation, and providing it with a public road access. These are the features that change as a result of consolidation [53]. In the case of the multiple linear regression model, potential effects of consolidation projects were assessed for different configurations of parcels' elongation and road access. In the case of the artificial neural network (ANN), data in various configurations of parcels' elongation, size classes, and road access were fed at the input. The results of these simulations were (Table 4) presented separately for the multiple linear regression model and the artificial neural network model. The computations were carried out assuming that other important factors with a potential impact on market prices were constant (100% share of agricultural land, peripheral location, no buildings in the vicinity, no obstacles to usage, and moderate quality of soil). A sample result of such estimation was presented for the municipality Gaworzyce (Lower Silesia voivodeship, Poland). With respect to other market characteristics, it was assumed that a hypothetical parcel is characterized by a 100% share of arable land, peripheral location, no buildings in the vicinity, no obstacles to usage, and moderate soil quality. The simulations run using both the models presented a positive impact of the change in agricultural parcels' parameters on the changes in parcels' values. Both the regression model and neural network model showed that the value of 1 ha of agricultural land increased significantly when a parcel was provided with a road access and its excessive elongation was reduced. In the analyzed sample of the municipality Gaworzyce, the value of 1 ha of agricultural land, once the parcel's elongation was corrected and the parcel was provided with an access to a public road, increased, according to the regression model, by about PLN 10.5 thousand (i.e.: a-b). The neural network predicted the same effect at the level of PLN 24.5 thousand (for large parcels, i.e., c-d) to PLN 28 thousand (for very small parcels, i.e., e-f). In addition to that, the neural network model revealed a non-monotonic impact of the changes in a parcel's size on the value of consolidated agricultural land. The biggest value was assigned by the neural network to very small parcels. Small parcels (in various variants of elongation and access to a road) were cheaper than very small parcels by PLN 3.3 thousand/ha (i.e., f-g) to PLN 4.3 thousand/ha (i.e., e-h). The neural network model indicated that the values of very big parcels were comparable to those of the smallest parcels. An interesting conclusion resulting from this simulation for the practical implementation of consolidation projects in the municipality Gaworzyce was as follows: to maximize the increase of the value of consolidated land, possibly the largest parcels should be designed. As Table 4 shows, even if consolidation projects implemented in the municipality Gaworzyce resulted in the elimination of very small parcels (up to 0.5 ha), which were replaced by small parcels (from 0.5 to 1 ha), middle-sized parcels (from 1 to 2 ha), or large parcels (from 2 to 3 ha), their lower value would be offset by the elimination of their excessive elongation and providing them with an access to public roads.

Discussion
Both the approaches used for assessing the impact of the changes in parcels' features on their value can be used to assess this phenomenon separately for each of the analyzed municipalities. The regression model did not take into account a parcel's size. Yet, it showed clearly that as the elongation of parcels increases, their values decrease. A road access was also significant here.
A fuller assessment of the effects of the changes in parcels' parameters was provided by the neural network model. Because this approach can show non-linear and non-monotonic relationships, it was possible to analyze the impact of the factor of a parcel's size. It is worth stressing that the ANN model assessed this aspect individually for each local market-i.e., the area of a given municipality.
In the case of the municipality Gaworzyce, the smallest parcels (up to 0.5 ha) showed higher unit prices, with similar trends observed across all the analyzed municipalities. However, based on the model, the largest parcels (>3 ha) had almost the same value on the local markets as the smallest ones. This relationship between the value of a parcel and its size can be explained by the fact that small-sized parcels are concentrated in areas where the construction of single-family houses is often likely in the future, even though the current regulations do not provide for such a possibility. The high price (per unit area) of large parcels is connected to their much better agricultural suitability, both with respect to current economic activity and in the context of greater potential interest in the purchase of such parcels by those wishing to expand their farms through purchase or lease of land.
Similar conclusions were drawn in the past by other scholars, who pointed out that a property's size impacts the attractiveness of parcels and the possibility of using them for specific purposes [54]. In agriculture, a parcel's size determines the amount of labor. Large parcels with regular shapes make it easier to carry out agrotechnical operations. Surface area is one of the most important features taken into account when estimating the value of agricultural land [55]. According to [56], the prices of larger parcels are often lower, whereas the prices of smaller parcels per unit area are higher.
The ANN model showed that in the municipality Gaworzyce (as well as in the other municipalities where market analysis was conducted), parcels with a model elongation and a road access across all the size classes had a relatively higher value (similar relationship was presented by the MRA model). The impact of the shape of parcels on the effects of agricultural works, as well as on the prices of agricultural land, is highlighted in numerous academic works [57,58]. A parcel's elongation, which is determined by the ratio of the parcel's width to its length, impacts the duration of the work in the fields and the loss in yields [59]. Similarly, the effectiveness of the work in the fields is impacted by the shape of the agricultural property. According to [60], the shape of a parcel should be based on a statistical similarity to a rectangle. Parcels with regular shapes make it easier to carry out agrotechnical operations. According to [16], the shape and slope of a parcel may significantly determine the prices of agricultural land.
The sensitivity analysis showed that the neural network assigned considerable importance to the issue of access to a public road. This observation was confirmed in the simulations that were run (Table 4). In each size class of parcels and with each variant of elongation, the differences in parcels' values in that respect were significant.
The results confirm the principle already indicated in the literature [61] that in the case of little diversity in a given cadastral region in terms of climatic, water, and soil conditions, the impact of location on the prices of agricultural properties increases. In such areas, the prices are most determined by distance-related location factors, i.e., distance from settlement development, distance from the municipality center, and distance from a hard paved road, as well as morphometric characteristics, the size and shape of a parcel [62].

Conclusions
Based on the simulations, it can be concluded that a change in value is significantly determined by the range of changes in parcel features such as size, shape, and access to a road. The impact of changes in surface area on the value of 1 ha of agricultural land is non-monotonic: the most expensive were the smallest and the largest parcels. Potential decreases in the value of agricultural parcels as a result of changes in their surface areas are successfully offset by elimination of the excessive elongation of parcels and providing an access to a public road.
Although the results were presented for only one municipality, it is worth highlighting the significant importance of general location-this feature had a large impact on the value of agricultural land, according to MRA and ANN models. The created models enable assessment of the impact of the implementation of a consolidation project on the value of agricultural land in w given region. Identifying universal rules seems difficult due to wide regional variations of the factors impacting the prices of agricultural properties and their changeability over a longer period of time. However, an attempt at such a universal approach can and should be a subject of future research using a range of parameters that adjust the model of changes in the value of land to local conditions.
Having said that, change in the value of land should be considered one of the effects of implemented land consolidation projects, during which the parameters of agricultural land are subject to the greatest changes.