Investigating the Effect of Tractor’s Tire Parameters on Soil Compaction Using Statistical and Adaptive Neuro-Fuzzy Inference System (ANFIS) Methods

: Many factors contribute to soil compaction. One of these factors is the pressure applied by tires and tillage tools. The aim of this study was to study soil compaction under two sizes of tractor tire, considering the effect of tire pressure and trafﬁc on different depths of soil. Additionally, to predict soil density under the tire, an adaptive neuro-fuzzy inference system (ANFIS) was used. An ITM70 tractor equipped with a lister was used. Standard cylindrical cores were used and soil samples were taken at four depths of the soil inside the tire tracks. Tests were conducted based on a randomized complete-block design with three replications. We tested two types of narrow and normal tire using three inﬂation pressures, at trafﬁc levels of 1, 3 and 5 passes and four depths of 10, 20, 30 and 40 cm. A grid partition structure and four types of membership function, namely triangular, trapezoid, Gaussian and General bell were used to model soil compaction. Analysis of variance showed that tire size was signiﬁcant on soil density change, and also, the binary effect of tire size on depth and trafﬁc were signiﬁcant at 1%. The main effects of tire pressure, trafﬁc and depth were signiﬁcant on soil compaction at 1% level of signiﬁcance for both tire types. The inputs of the ANFIS model included tire type, depth of soil, number of tire passes and tire inﬂation pressure. To evaluate the performance of the model, the relative error ( ε ) and the coefﬁcient of explanation (R 2 ) were used, which were 1.05 and 0.9949, respectively. It was found that the narrow tire was more effective on soil compaction such that the narrow tire signiﬁcantly increased soil density in the surface and subsurface layers.


Introduction
Soil is a resource with a potentially rapid degradation rate and a very slow regeneration and formation process [1].Thus, sustainable soil use is the only solution to global problems such as food security, energy and water demands, climate change and biodiversity [2].In order to increase productivity in the agricultural sector, farms should be a large size, in which productivity can be increased by the use of more efficient machines [3].For many years, the agricultural sector has tended to increase the size and weight of the tractor, which increases the risk of severe soil compaction [4].
One of the causes of soil compaction in mechanized agriculture is the passage of agricultural implements and tools.Soil compaction as a form of reducing the pore space between soil particles can significantly reduce soil production capacity [5].According to data from seven countries in Europe and North America, a 14% drop in crop yield was reported in the first year after machine traffic [6].Soil compaction is a form of physical degradation of soil that causes compacted soil particles to be closer to each other and the disappearance of voids in the soil, reducing porosity and soil permeability.This manifests as an increase in soil bulk density [7].One method of assessing soil compaction is to measure its bulk density.When soil compaction occurs, the bulk density of the soil increases as the mass remains constant and the volume decreases.Density reduction decreases hydraulic conductivity and thus causes surface runoff and soil erosion by water.This induces forced currents in the soil pores, thus facilitating the transport of absorbable colloidal nutrients and pesticides to deeper horizons and groundwater, thereby reducing root growth and causing a loss of nitrogen and production of greenhouse gases through nitrogenization by anaerobic processes.As a result, soil compaction is one of the main causes of environmental and agricultural problems, and as a result causes significant economic damage to society and agriculture [8].
Increasing the load on the wheels of agricultural machinery is the most important cause of soil compaction created in deeper soil layers.The topsoil of arable lands is loosened every year by plowing operations, however, the subsoil often remains intact, which means that whenever the maximum load capacity is exceeded, additional soil compaction will occur deep in the soil.In addition, these effects are lasting, and the results of natural and artificial soil loosening are disappointing [9].Excessive movement of agricultural machinery on the field and plowing at a constant depth gradually causes soil compaction, and this compaction restricts plant root growth and reduces crop yield [10].To remove this hard layer created in the subsoil, which limits the longitudinal growth of roots and proper soil drainage, deep plowing is required, which greatly increases the required energy [11].Another solution is to use energy-efficient combined (active-passive or passive-passive) tillage implements in place of conventional ones to reduce soil compaction, labor, and fuel cost as well as to save time in preparing the seedbed [12].Another aspect affecting the protection of soil structure from compaction is the evolution of tires [13].In this article we compare the effects of five tractor tires introduced between the 1970s and 2017 on soil stress and soil structure.Soil physical properties before and after one pass of each of the five tires were measured on undisturbed soil.The results of stress measurements showed the lowest mean normal stress at all depths with newer tires, both below the center of the tire and at the edge, even though the effect decreased with increasing depth.The results indicate that tire design can play an important role in reducing soil stress for a given wheel load.Tire design might then also help to reduce the compaction risk for larger wheel loads and other soil conditions which were tested in the study.
Maintaining proper ratios between solid, liquid and gas phases is of particular importance, because increasing the solid phase (i.e., reducing the liquid and gas phases) increases the mechanical strength of the soil and consequently reduces the growth and development of plant roots inside the soil.Reducing the development of root growth causes less water and nitrogen to be available to the plant.When the machine enters the field, the solid phase of the soil increases by a certain volume and the gas phase decreases.Soil density is one of the indicators of the destruction of the physical structure of the soil, which is defined as an increase in soil bulk density (solid phase mass to total volume), or a decrease in porosity [14].Wheel load, tire type, and inflation pressure can increase soil bulk density and play an important role in soil compaction [15].Almost all tires significantly increase soil compaction along the wheel track, while only some increase soil compaction near the track [15].As the distance from the wheel path increases, a general decrease in soil density occurs, especially in the subsoil.Many researchers have reported that working with low-pressure tires can significantly reduce soil compaction and increase crop yield [16,17], while high tire pressure increases soil compaction [18].Tire pressure control is an important factor in the control of and in soil compaction.A new tire pressure control system in forest machines was developed using a PressurePro solution.It was found that the installation of an automatic tire pressure control system leads to compaction reduction by 20%.Additionally, reduced tire pressure and its automatic control contribute to a minimal reduction in humus content and soil compaction over time [19].Another researchers have investigated the effect of different tire inflation pressures of 100 and 200 kPa on selected physical soil properties by conducting field tests.The undisturbed samples were collected both in and between the tracks at depths of 0 to 0.5 m.The results indicated that fewer negative changes were found in the variant with a lower inflation pressure for all of the soil properties.Change in soil physical properties caused by the passage of the tractor were statistically significant for both tire pressures only at depths ranging from 0 to 0.1 m [20].The effect of three types of wide, low-inflation pressure tires with similar dimensions on mean normal stress on the soil profile was also investigated.The results revealed a very limited effect of tire construction on mean normal stress.No differences were measured beneath the center, and differences at +0.3 m were found only at 0.2 m depth for the tires at the rear axle.The effect of low tire inflation pressure was limited to the upper parts of the soil profile for beneath the center of the tire, which was significant at depths of 0.2 m and 0.4 m.It was found that in order to reduce soil stress, tire design should allow for a large contact area and low inflation pressure [21].
Machine traffic plays an important role in soil compaction.Experiments have been carried out to study the effect of multiple passes of tillage operation on soil compaction beneath the tillage depth with three different tillage systems.The authors observed that active tillage machinery generally caused more compaction compared to passive tillage machinery for the same number of passes.However, active machinery required fewer vehicular passes [22].In plowed soils, the first pass can dramatically increase soil compaction and cause soil degradation, while in firm soil and deep areas soil compaction occurs over time and with large numbers of tire passes [23].Research has also shown that if the traffic exceeds 10 times, the benefits of using light tractors are lost and soil compaction is greatly increased.These studies have suggested that heavy tractors with less traffic should be used as much as possible because high traffic leads to compaction in the subsoil [24,25].
One method of assessing soil compaction is to measure its bulk density.When soil compaction occurs, the bulk density of the soil increases as the mass remains constant and the volume decreases.In this method, after determining soil bulk density at different depths and, based on the soil texture, measured data are compared with the standard values to determine the level of soil compaction.For example, soil compaction is considered significant for loamy clay soil when bulk density reaches more than 1.5 g/cm 3 , and for sandy and loamy soil more than 1.6 g/cm 3 [26].
In recent years, the artificial intelligence approach has shown that it has great power to model and predict complex systems, thus it can be an alternative method to physical models that do not respond well to a large number of input variables [26].Because agricultural systems such as soil are so complex, most researchers have focused on artificial intelligence to model different components of agricultural systems.The fuzzy neural adaptive inference system (ANFIS) is a multi-layered adaptive network consisting of the main elements and functions of fuzzy logic systems developed by Jang [27].Each of the fuzzy systems and artificial neural networks (ANNs) has advantages and disadvantages.The fuzzy system is able to use human language and can use human experiences and experts while it is not able to learn [28].In other words, the fuzzy system cannot be trained using observational data and it should be noted that this method does not provide good results in unpredictable conditions [29].However, neural networks are self-learning using data sets.At the same time, neural networks are implicit and unable to use human language.They require a wide range of tested input and output data for successful modeling [30].Although ANN is a powerful technique for modeling various problems in the real world, it has its own weaknesses.If the input data is vague or the uncertainty is relatively high, a fuzzy system such as ANFIS may be a better option [31].
In order to evaluate the potential of ANFIS in predicting the energy efficiency indicators of the drive wheel, experiments were performed in the soil reservoir.Input parameters included tire load, wheel speed and slip, each at three different levels.ANFIS was used with a combined method of descending gradient and least squares to find the optimum learning parameters using different membership functions (MFS).The results of this study showed high accuracy (MSE = 0.0166 and R 2 = 0.98) in prediction [31].A fuzzy logic model was developed to describe soil fragmentation for seedbed preparation in a combination of primary and secondary tillage tools from subsoil, moldboard plow and disc.In this research, an intelligent model based on the Mardani fuzzy modeling principles approach was developed to predict soil fragmentation during tillage operations.Model inputs were: soil moisture, tractor speed and sampling depth.The fuzzy model consisted of 50 rules in which the three parameters of mean root mean square error (RMSE), relative error (e), and coefficient of determination (R 2 ) were used to evaluate fuzzy models.These parameters were calculated to be 0.17%, 3.95% and 0.99%, respectively.According to the results of this study, the fuzzy model can be used as a method to predict soil fragmentation during tillage operations [32].In the study, the ANFIS model was used to predict soil compaction under tractor wheels.Model inputs included four factors: manure rate, number of passes, moisture and depth.The network structure was selected from the network segmentation type.The number of membership functions for each input was considered to be three and the membership functions of bell, triangle, trapezoid and Gaussian were tested for best performance.The results of this study showed that the ANFIS model has a higher accuracy in predicting soil density than the regression model [33].
The main objectives of this research are as follows: 1.
Investigating the density created in different depths of the soil due to the passage of tractor wheels of different size; 2.
Investigating the effect of inflation pressure changes of drive wheel and machine traffic on the soil density created at different depths; 3.
Modelling the effect of tire size, tire pressure, and machine traffic on soil density at different depths using statistical and adaptive neuro-fuzzy inference system (ANFIS) methods

Materials and Methods
Field experiments were performed to measure the bulk density of the soil under two common and narrow radial tires.The first type of tire was a standard 8 layer tire with a size of 14-38.Three inflation pressures of 80, 120 and 160 kPa were used to evaluate the effect of tire pressure on soil compaction.The second type of tire was a 10 layer radial narrow tire with a size of 9.5-48.Inflation pressures of 240, 275 and 310 kPa were investigated.For each of these tire pressures, three traffics of 1, 3 and 5 were tested.In all tests, soil samples were obtained from four depths of 10, 20, 30 and 40 cm.The farm had been plowed already, and an ITM70 tractor equipped with a bedder-ridger on a three-point hitch was used for conducting tests.In this case, the weight measured on the rear axle was 2740 kg, and on the front axle 770 kg.In all experiments, the tractor speed was constant, equal to 3 km/h.The process of adjusting the desired speed was done in such a way that the tractor was moved to a distance of 30 m on the field to achieve a speed of 3 km/h with two gears and at engine speed of 1400 and 1600 rpm. the desired speed was obtained using a calibration diagram.
A hydrometric test (ASTM D-422) was used to determine the soil texture by determining the relative amount of sand, silt and clay.It was found that there were 43.92, 32.06 and 24.02% sand, soil and clay, respectively.According to the soil texture triangle, the soil texture was loamy.A standard proctor test was used to determine the soil critical moisture.According to Figure 1, a moisture level of 25% is the critical moisture content, and the highest value of compaction occurs at this humidity.
To measure the bulk density after the tire passes, standard cylinders were placed in groups of 3 at different depths of 10, 20, 30 and 40 cm measured from the center of the cylinders to the soil surface and in line with tire path direction in the center line of the tire track (Figure 2).In order that the cylinders were placed exactly at the desired height and in one direction inside the profile, a string was drawn from the middle of the profile and from the beginning to the end.Figure 3 shows how the cylinders were exerted inside the soil vertically in the direction of the string.Each group of cylinders was about 40 cm apart so that there were no problems when removing the cylinders.It should be noted that the soil surface appears loose in Figure 3; this is because the soil had been dug to place the cylinder at a predetermined depth.To implement the fuzzy neural adaptive model, MATLAB 2014a software toolbox under Windows was used.According to the performed experiments, soil density was defined as system output, while tire type, tire inflation pressure, traffic and soil depth were independent inlet parameters (Figure 4).The data is usually divided into three categories: training data, validation data, and test data.The training data were used in the training process to calculate the gradient and optimize the parameters.One of the points that should be considered during the training is aiming to prevent the network from becoming specialized and the occurrence of over-fitting phenomena.For this purpose, a series of data was considered for validation during the training of the network.These data are actually part of the training data.In this way, in the regular intervals of the optimization process, the data obtained from the network were checked using the validation data.In this case, network training continued until the optimization error related to the evaluation data began to increase, and as soon as this error increased to a certain value or a certain number of repetition, network training was stopped.Finally, the test data was used to test the performance of the network after training.The test data was not used during network training, but it was used to compare the error rate.Accordingly, in this research, 70% of the data were considered for training, 15% for evaluation, and 15% for testing the network.A network separation structure was considered to create the network.Four commonly used membership functions, namely triangular, trapezoidal, Gaussian, and bell, were considered as membership functions to represent the inputs (Figure 5), and four different models were constructed using these membership functions.
The number of membership functions for each variable was considered proportional to the levels tested for that variable.Table 1 shows the membership functions for each variable.
The output membership function in this network is linear.A hybrid optimization method was used for network training.The number of rules created by the network was 54.Given that the number of rules created was very high, only a few of them are presented in Table 2.  Finally, the performance of the models was evaluated based on the statistical criteria presented in Equations ( 1) and (2). where: MAPE: the absolute mean percentage of system error R 2 : is the coefficient of determination N: the number of samples Y predicted : the predicted values Y measured : the measured values Y measured : the average of measured values.

Results and Discussion
Table 3 shows the analysis of variance of soil bulk density for two types of ordinary and narrow tire.Data were analyzed in a 4 × 3 × 4 factorial design based on a randomized complete block design.The data coefficient of variation (CV) was 0.44%.Average initial soil bulk density was 1.017 g/cm 3 , and based on data of initial bulk density it was uniform in all plots.In comparison with initial bulk density, statistical analysis showed that the main effects of tire type, traffic and depth on soil density change were significant at the level of 1% probability.The interactions of depth with traffic, depth with tire type, and traffic with tire type on soil bulk density were found to be significant at the 1% probability level.However, the combined interaction of depth, tire type and traffic on soil bulk density was not found to be significant at the 1% probability level.

The Binary Effect of Tire Type and Depth on Soil Bulk Density
The binary effect of tire type and depth on soil bulk density was evaluated using a Duncan test at 1% probability level, and is presented in Figure 6.According to the results, the maximum value of bulk density was 1.217 g/cm 3 for narrow tire at a depth of 10 cm, and the lowest value was 1.061 g/cm 3 which occurred for common tires at a depth of 40 cm.For both tires, a significant decrease in the soil bulk density was seen with increasing depth, but the intensity of these changes in the narrow tire and especially in shallow depths was greater than for the common tire.The bulk density of soil for the narrow tire, even at a depth of 40 cm, was greater than the bulk density of soil for a normal tire at a depth of 10 cm.The effect of tire size on stress applied to different soil layers was examined [34].Their results showed that a narrow tire created much more stress in different layers of soil than a common tire, and as a result, it compacted the soil more.Stress for a narrow tire at a depth of 50 cm was higher than the stress of a common tire at depth of 30 cm.

The Binary Effect of Tire Type and Traffic on Soil Bulk Density
The results of comparing the mean interaction of tire type and depth on soil bulk density using a Duncan test at 1% probability level are presented in Figure 7.In both tires, the soil bulk density increased with increasing traffic.According to the results, the highest mean value of soil density was found for the narrow tire after five passes with a value of 1.185 g/cm 3 and the lowest for a normal tire after one pass with a value of 1.112 g/cm 3 .It is noteworthy that the bulk density of soil for a narrow tire, even at one pass, is more than the bulk density of soil for a normal tire in five passes.

Results of Modeling Soil Bulk Density Using ANFIS
Four different models were developed to predict soil bulk density.Table 4 presents the structural parameters of the models along with their statistical criteria to evaluate their performance.According to the results presented in Table 4, it is clear that all models have a high ability (R 2 ≥ 0.99 and ε (%) ≤ 2) to predict.The best model for predicting soil bulk density is a model that uses the bell membership function (Gbellmf) (R 2 = 0.9949 and ε (%) = 1.05).Table 5 presents the statistical characteristics of the stepwise regression model for predicting soil bulk density.The ANFIS models have poorer performance than the regression model due to a low coefficient determination of 0.96.However, regression models have valuable advantages, such a the fact that ANFIS models do not provide any information about the internal structure of the model and the relationships between independent and dependent variables, while the regression model directly deals with the impact of each factor and tries to provide a model in which the importance and impact of each factor is clear.According to Table 5 and the standard coefficients mentioned therein, type of tire, soil depth, inflation pressure, and finally the number of passes have the greatest impact on the bulk density.Figure 8 shows the relationship between measured and predicted values using ANFIS models under different operating conditions.In addition, the deviation between the measured values and the predicted values was calculated and plotted by ANFIS and multiple regression models.Figure 9 shows that the deviation distance of the predicted values of the ANFIS model (−0.011 to 0.0085) was much less than the deviation distance of the predicted values of the multiple regression model (−0.021 to 0.033).Figure 10 shows the three-dimensional variation of soil bulk density with the interaction of tire type and soil depth.It indicates that with increasing depth, the bulk density of the soil decreased, indicating that the compaction created by the tire traffic of the tractor in the topsoil was higher than in the sub layers.It was found that at all depths the density change by the passage of common tire was less than the density created by the narrow tire.Such a low pressure at the soil surface and soil behavior was investigated using a flat tire [3]. Figure 11 shows the binary effect of tire pressure and tire traffic on the soil bulk density.As the tire pressure decreases, the bulk density of the soil decreased.As the pressure decreased, the contact area between the tire and the soil increased, which reduces the pressure on the soil, resulting in decreased soil compaction and this was reported by other researchers [35][36][37].With increasing traffic, the bulk density of the soil increased [33,38,39].Figure 12 illustrates with a surface plot of the interaction of tire type and traffic on the bulk density of the soil.With increasing traffic, when using a narrow tire, the soil bulk density increased sharply, but using a common tire was less effective on soil bulk density increment.This indicates that narrow tires are very destructive on soil structure.Figure 13 shows the interactive effect of tire type and tire inflation pressure on soil compaction.The soil bulk density changed strongly with increasing pressure when using a narrow tire, but less changes occurred when using a wide tire with increasing pressure.The results showed that using a narrow tire increased the soil bulk density in different layers of soil much more in comparison with a common tire.A narrow tire compacted the soil more and increased the soil density, even at a depth of 50 cm, more than a common tire at a depth of 30 cm.Additionally, the soil bulk density of soil for a narrow tire at a depth of 40 cm was greater than that of a normal tire at a depth of 10 cm. 2.
In both tires, the soil bulk density of the soil increased with increasing tire wheeling.According to the results, the highest mean value of soil density for narrow tire was after five passes, with a value of 1.185 g/cm 3 .The lowest value was for a common tire after one pass, with a value of 1.112 g/cm 3 .It is noteworthy that the bulk density of soil for a narrow tire after one pass was more than the bulk density of soil for a common tire after five passes.

3.
Using a narrow tire with tire wheeling increased the soil bulk density very sharply, but when using a wide tire, with increasing passes, fewer changes in the bulk density of the soil occurred.This indicates that using narrow tires should be considered to be very destructive on the soil structure.Additionally, with increasing tire pressure, and when using a narrow tire, the soil bulk density changed much more than when the common tire's pressure changed.

Figure 1 .
Figure 1.Determining the critical soil moisture using the proctor test.

Figure 2 .
Figure 2. Passing the tractor with two types of tires on the field.

Figure 3 .
Figure 3. Placement of the cylinders inside the soil profile.

Figure 6 .
Figure 6.The main effect of tire type and depth on soil bulk density (different English alphabets shows significant difference between treatments).

Figure 7 .
Figure 7.The main effect of tire type and traffic on soil bulk density (different English alphabets shows significant difference between treatments).

Figure 8 .
Figure 8. Relationship between measured and predicted values of soil bulk density.

Figure 9 .
Figure 9. Deviation distance of predicted values of ANFIS and regression models from measured values.

Figure 10 .
Figure 10.Binary effect of tire type and soil depth on the soil bulk density.

Figure 11 .
Figure 11.Binary effect tire inflation pressure and traffic on the soil bulk density.

Figure 12 .
Figure 12.Tire type and traffic effect on the soil bulk density.

Figure 13 .
Figure 13.Three-dimensional surface curve of the main effect of tire type pressure and tire pressure on the soil bulk density of the soil.

Table 1 .
Membership functions of input variables.

Table 2 .
Some rules created in the ANFIS model.

Table 3 .
Analysis of variance of soil bulk density change due to tire type, traffic and depth.

Table 4 .
Specifications and evaluation results of adaptive fuzzy-neural inference system models (ANFIS).

Table 5 .
Statistical characteristics of stepwise regression model for predicting soil bulk density.