with the collaboration of Iranian Society of Mechanical Engineers (ISME)

Using the Response Surface Methodology to Predict the Effect of Different Moisture Levels on the Bulk Density and Penetration Resistance of Soil Under Different Operating Conditions

Document Type : Research Article- En

Authors

1 Department of Agriculture Machines and Equipment, College of Agriculture, University of Basrah, Iraq

2 Department of Soil and water Science, College of Agriculture, University of Basrah, Iraq

Abstract
Soil properties play a fundamental role in the success of agricultural operations through their impact on crop growth and quality, as they determine their ability to retain water and absorb nutrients, and affect soil aeration and the root system. The aim of this study is to predict bulk density and resistance to soil penetration under different moisture levels during tillage operations. It includes four moisture levels: 7, 14, 22, and 28%, and three types of plows: the moldboard plow, chisel plow, and disc plow. Moreover, soil samples were collected at two depths: 15 cm and 30 cm. The change in the physical properties of the studied soil is also measured during the growth periods of wheat crop (after tillage, beginning of the season and end of the season). The study is conducted in Al-Qurna district, north of Basra Governorate, Iraq, in clay loam soil. The results are analyzed and mathematical equations are obtained to predict the studied properties using the response surface methodology. The obtained results indicate that soil moisture during plowing, plow type, soil depth, and crop growth periods have a significant effect on soil bulk density and penetration resistance. The 14% moisture treatment is superior, recording the lowest bulk density and lowest penetration resistance of 1.12 Mg m-3 and 1133 kN m-2, respectively. While the 28% moisture treatment provided the highest bulk density and highest penetration resistance of 1.22 Mg m-3 and 1379 kN m-2, respectively. The results also show that increasing the soil depth from 15 to 30 cm increases the bulk density and soil penetration resistance, by 12 and 45.70%, respectively. Plowing with a disc plow improves soil properties, giving the lowest bulk density and penetration resistance of 1.12 Mg m-3 and 1074 kN m-2, respectively. While using the chisel plow leads to recording the highest bulk density and penetration resistance, which reached 1.22 Mg m-3 and 1442 kN m-2, respectively. As for the moldboard plow, the bulk density and soil penetration resistance reached 1.18 Mg m-3 and 1282 kN m-2, respectively. The growth periods have a significant effect on the studied soil properties where the beginning of the growing season provided the lowest bulk density. The bulk density reached 1.17, 1.13, and 1.23 Mg m-3 for the periods after plowing, at the beginning of the season and its end, respectively. While the penetration resistance after plowing is superior with the lowest resistance compared to the beginning of the season and its end, as it reached 897, 1327, and 1573 kN m-2, respectively. The results of data analysis show that the obtained mathematical models accurately and efficiently predict bulk density and soil resistance to penetration under the experimental conditions, with a high coefficient of determination (R2) of 0.6460 and 0.8114 for the bulk density and penetration resistance, respectively.

Keywords

Subjects

Introduction

Tillage greatly affects many physical properties of soil, such as bulk density and porosity. It breaks up the soil to reduce its compaction and increase its porosity, thereby aiding in weed control and increasing crop production. However, improper tillage may lead to soil hardening and deterioration of its physical properties, which negatively affects aeration, root growth, and microorganism activity, thus reducing production. Therefore, choosing the appropriate type of tillage is essential to achieve the best productivity (Boydas and Turgut, 2007; Shabanpour, Fekri, Bagheri, Payman, and Rahimi-Ajdadi Shabanpour, Fekri, Bagheri, Payman, and Rahimi-Ajdadi 2022). Agricultural work greatly affects the physical properties of soil and the moisture level during tillage. Tillage under conditions of high or low moisture can lead to the formation of large soil clods and deterioration of the physical properties of soil (Shittu, Oyedele, and Babatunde, 2017). On the other hand, tillage contributes to improving the physical properties of soil, such as reducing bulk density, increasing porosity, and improving soil resistance to penetration.

Bulk density is an important physical property of soil, and it is greatly affected by tillage and moisture level. According to a study by Nassir (2018), the optimum soil moisture content of 16.47% achieved the best results for bulk density (1.16 Mg m-3) and soil penetration resistance (678.57 kN m-2), compared to moisture levels of 10.23% and 24.68%, which resulted in bulk density of 1.36 and 1.20 Mg m-3 and penetration resistance of 788.16 and 835.86 kN m-2, respectively. Ahmadi and Ghaur (2015) showed that soil bulk density increases with soil moisture at 12, 15, 17, 19, and 21%. Soil compaction is influenced by various factors, such as tractor movement across the field, the number of passes made, the type of tillage employed, the inherent properties of the soil, and its moisture content during tillage. Soil compaction is usually expressed in terms of bulk density, porosity, or soil resistance to penetration (Javadi and Spoor, 2006; Rashidi, Tabatabaeefar, Keyhani, and Attarnejad, 2007). According to 2Ahmadi and Mollazade (2009), tillage at 13-15% soil moisture reduced soil resistance to penetration by 40%, whereas at 15-18% moisture, the reduction was only 4.9%. Soil resistance to penetration depends on soil type, water content, clay content, bulk density, soil depth, and tillage system. Tillage equipment has a significant impact on soil physical properties, such as bulk density and penetration resistance (Naderi-Boldaji, Azimi-Nejadian, and Bahrami, 2024; Tahmasebi, Gohari, Sharifi Malvajerdi, and Hedayatipour, 2023). A study by 14Kostić, Rakić, Savin, Dedović, and Simikić (2016) showed that the type of tillage affects the bulk density of soil, with density being 1.50, 1.47, and 1.45 Mg m-3 for the moldboard plow, chisel plow, and disc plow, respectively. Bulk density increases with increasing soil depth due to higher soil strength, with bulk density ranging from 1.33 to 1.38 Mg m-3 when the depth increases from 15 to 50 cm (Salim, Almaliki, and Nedawi, 2022).

Soil penetration resistance is an indicator of soil hardness, as soil with high resistance can hinder root spread, lead to waterlogging, and decrease aeration, which negatively affects crop growth. Therefore, tillage operations are carried out to break up the soil and reduce penetration resistance, which promotes root spread and improves soil physical properties (Kuroyanagi, Kaneko, Watanabe, Fujita, and Odahara, 1997). Several studies have shown that tillage reduces soil penetration resistance compared to no tillage (Hajabbasi, 2010; Kahlon, Lal, and Varughese, 2013), and that the plow contributes to increased penetration resistance compared to other conventional tillage methods. In addition, increasing tillage depth increases soil penetration resistance (Biberdzicet al., 2020; Dekematiet al., 2019; Kuhwaldet al., 2016).

Neural networks have been used in several studies on agricultural tillage equipment to predict energy requirements and evaluate the performance of tillage equipment based on variables such as moisture, tillage depth, and plow type (Almaliki, Himoud, and Al-Khafajie, 2019), showing high agreement with field experimental data. This method is fast, accurate, and low-cost compared to conventional methods. Therefore, these techniques can be used to predict soil properties under different conditions. Neural networks have also been used to predict soil disintegration during tillage and its effects on water movement, bulk volume, water drainage, moisture content, and soil bulk density (Taghavifar and Mardani, 2014). The tillage process is influenced by both the type of plow used and the soil moisture content at the time of tilling. Given that assessing soil properties after tillage and throughout the growing season can be both labor-intensive and costly, this research seeks to predict two critical soil characteristics—bulk density and penetration resistance. These factors are essential indicators of tillage quality and favorable growth conditions. This study will investigate the impact of varying moisture levels on soil conditions, utilizing three types of plows including moldboard, chisel, and disc and examining two soil depths of 15 cm and 30 cm. Measurements will be taken at three key intervals: immediately after tillage, at the start of the growing season, and at its conclusion.

Materials and Methods

Field experiments

The field experiment was conducted in Al-Qurna district in Basra governorate in Iraq on clay loam soil. The work began with determining the moisture content of the soil at plowing by experimenting with enclosing a certain area of the soil and flooding it with water, then samples are taken every two days to measure the change in soil moisture. Based on the data obtained, the required moisture levels for the experiment are determined. The field is divided into four sectors, each with an area of 1600 m2, and each sector is irrigated at different intervals according to the specified moisture levels, which are 7%, 14%, 22%, and 28% (depending on the limits of plasticity). Three types of plows are used for each sector: a three-furrow moldboard plow with a working width of 1 m, a three-furrow disc plow with a working width of 1.0 m, and a chisel plow with 11 shanks arranged in three rows with a working width of 2.2 m. Plowing speed of 3.06 km h-1. Soil samples are taken to measure the apparent density and penetration resistance after plowing at two depths of 15 and 30 cm. After preparing the field for cultivation, it is divided into 36 experimental units. Each unit area is 12 m2 (6 × 2 m), suitable for using four moisture levels, three types of plows, and three replicates for each treatment. The field is planted with wheat (Triticum aestivum L.) of the research variety 22. Soil samples are collected after plowing, at the beginning of the growing season, and at the end of the season before harvest to evaluate the changes in the studied physical properties during the season.

Studied characteristics

FBulk density

Bulk density is measured by taking undisturbed soil samples using a core sampler, following the method described by Black, Evans, White, Ensminger, and Clark (1965). The soil samples are weighed before drying, then dried in an oven at a temperature of 105°C until a constant weight is reached. Bulk density (ρb) is calculated using Equation (1):

ρb=MSV(1)

where: ρb = Bulk density of the soil (Mg m-3); MS = Mass of the solid particles (Mg); V = Total volume of the soil, which is the volume of the cylinder (m³).

Soil penetration resistance

To assess soil penetration resistance, we utilize a Dutch-made field cone penetrometer from Eijkelkamp Agrisearch Equipment. This device applies variable pressure vertically onto the soil surface, and each treatment is tested using three replicates. The cone index (CI) is calculated mentioned in ASABE Standards (2009) as:

Cone Index(CI)=Penetration forceCone base area(2)

where Cone Index (kN m-2); Penetration force (kN); Cone base area (m2).

Mathematical model

The response surface methodology is used to develop mathematical models and analyze data to predict the bulk density and soil penetration resistance. In this study, 216 experiments are conducted, including the use of three types of tillage machines (moldboard plow, chisel plow, and disc plow), four moisture levels (7%, 14%, 22%, and 28%), and three crop growth stages (after tillage, beginning of the season, and end of the season), and measurements are made at two different soil depths. The study aims to develop accurate models for the bulk density and soil penetration resistance to evaluate the effect of these factors on soil properties during the growing season.

Results and Discussion

Bulk density

The results of the statistical analysis are shown in Table 2, demonstrating a significant effect of soil moisture on the bulk density of the soil. Figure 1 shows that the bulk density of the soil increases with increasing moisture content from 7% to 28%. Soil with 14% moisture exhibited the lowest bulk density, measuring 1.12 Mg m-3, while soil with 7% moisture had a slightly higher density of 1.17 Mg m-3. While there is no significant difference between the moisture at 22% and 28%, as the bulk density reached 1.20 and 1.22 Mg m-3, respectively. The superior bulk density achieved by the soil at 14% moisture is due to the improvement of the mechanical properties of the soil, such as reduced cohesion and adhesion, which facilitated the disintegration of the soil during plowing, thus reducing its apparent density. As for the moisture content of 22% and 28%, the cohesion and adhesion of the soil increased, which led to soil compaction and an increase in its bulk density. This is in line with the results of the study by Nassir (2018), which indicated that higher moisture levels lead to increased soil cohesion.

Characteristics Bulk density (Mg m-3) Penetration resistance (kN m-2) Electrical conductivity (ds m-1) Cohesion (kN m-2) Adhesion (kN m-2)
15 cm 30 cm 15 cm 30 cm 15 cm 30 cm 15 cm 30 cm
1 1.42 1.44 1700 1800 14.98 14.96 10.7 10.78 0.0867
2 1.11 1.25 1200 1333 13.3 7.55 5.34 6.99 0.1263
3 1.24 1.28 1066 1133 10.71 7.36 7.71 8.75 0.1362
4 1.35 1.45 820 850 3.66 9.03 9.58 9.92 0.304
Table 1. Primary soil properties
Source Sum of squares df F-Value p-value (Prob > F)
Model 2.85 19 18.82 < 0.0001
A-Moisture content 0.15 1 19.45 < 0.0001
B-Depth 1.19 1 150.89 < 0.0001
C-Growing season 0.40 2 25.38 < 0.0001
D-Plow type 0.38 2 24.19 < 0.0001
AB 0.025 1 3.12 0.789
AC 1.647E-003 2 0.10 0.9011
AD 0.018 2 1.15 0.3175
BC 0.55 2 34.75 < 0.0001
BD 0.087 2 5.51 0.0047
CD 0.016 4 0.50 0.73771
Table 2. Analysis of variance for the effect of study factors on the bulk density of soil

Fig. 1. Effect of soil moisture on the bulk density of soil (Mg m-3)

Figure 2 and the variance analysis table (Table 2) show that soil depth has a significant effect on bulk density. At a depth of 15 cm, the lowest bulk density was observed, measuring 1.10 Mg m-3, in contrast to the 30 cm depth where the density increased to 1.25 Mg m-3. This difference is due to the effects of tillage, crop growth, and root spread at a depth of 15 cm, which contributes to soil loosening and helps reduce bulk density.

Fig. 2. Effect of soil depth on soil bulk density (Mg m-3)

In contrast, the depth of 30 cm is relatively far from the root zone, and smoothing equipment did not reach it, which led to an increase in soil density at this depth. These results are consistent with the findings of Salimet al. (2022), where they found that the bulk density of soil increases with increasing depth from 15 to 50 cm, ranging between 1.33 and 1.38 Mg m-3.

The crop growth periods clearly affect the bulk density of the soil. As shown in Table 2 and Figure 3, the growth period has a significant effect on the change in bulk density. The soil recorded the lowest bulk density at the beginning of the growing season, reaching 1.13 Mg m-3, while this density increased to 1.23 Mg m-3 at the end of the season. After the plowing process, the density reached 1.17 Mg m-3. The decrease in density at the beginning of the season is due to the effect of smoothing and leveling processes carried out after plowing, in addition to the spread of crop roots, which contributed to reducing the bulk density. On the other hand, the bulk density increased at the end of the growing season as a result of repeated irrigation processes, which led to the movement of soil particles and their settlement in the pores, in addition to the stability of the soil over time. The bulk density following plowing is higher than at the start of the season, because the soil surface remains uneven from the plowing process. These results are consistent with the findings of Shabanpouret al. (2022), where an increase in the bulk density of the soil is observed after harvest compared to the beginning of the growing season.

Fig. 3. Effect of growth periods on soil bulk density (Mg m-3)

The results of the statistical analysis in Table 2 show a significant effect of the type of plow on the bulk density of the soil. As shown in Figure 4, plowing with a disc plow recorded the lowest bulk density of 1.12 Mg m-3, which is attributed to the nature of the disc plow's work, which is characterized by its ability to work in different field conditions. As it works to split and loosen the soil by rotating the discs, which leads to raising, turning, and loosening the soil. In contrast, the moldboard plow recorded a higher density of 1.18 Mg m-3, due to its method of operation that depends on turning the soil using the plow, which leads to an increase in the weight applied to the soil and the formation of more cohesive blocks compared to the disc plow. As for the chisel plow, it recorded the highest bulk density of 1.23 Mg m-3, due to its work on splitting the soil without turning it, which leads to loosening the soil locally and increasing its density compared to the reversible plows. These results are consistent with those of AbdulSada and Almaliki (2023).

Fig. 4. Effect of plow type on soil bulk density (Mg m-3)

The analysis presented in Table 2 indicates that there are no significant effects arising from the interactions between soil moisture and soil depth, soil moisture and growth periods, or soil moisture and plow type. Additionally, there is no significant interaction between growth periods and plow type with respect to bulk density. However, the results in Table 2 and Figure 5 indicated that there is a significant effect on the interaction between soil depth and growth periods on the bulk density. The depth of 15 cm at the beginning of the growing season recorded the lowest bulk density of 0.99 Mg m-3, while the depth of 30 cm at the beginning of the growing season recorded the highest bulk density. There is no significant difference between the depth of 15 cm and the depth of 30 cm at the end of the growing season, where the density reached 1.28 and 1.26 Mg m-3, respectively.

Fig. 5. Effect of the interaction between growth periods and soil depth on the bulk density of soil (Mg m-3)

The results of the statistical analysis show a significant effect on the interaction between the type of plow and soil depth. As shown in Figure 6, plowing with a disc plow at a depth of 15 cm recorded the lowest bulk density of 1.07 Mg m-3, without a significant difference compared to plowing with a moldboard plow at the same depth (1.08 Mg m-3). On the other hand, plowing with a chisel plow at a depth of 30 cm recorded the highest bulk density, with measurements of 1.29 Mg m-3. This value was not significantly different from that obtained with a moldboard plow at the same depth, which reached 1.28 Mg m-3.

Fig. 6. Effect of interaction between plow type and soil depth on soil bulk density (Mg m-3)

Table 3 shows the mathematical models for each plow during the crop growth periods to predict the bulk density of the soil under different field conditions. Through these equations, the bulk density of the soil can be predicted by entering the variables of soil moisture and soil depth.

Measurement time Plow type Bulk density equation
After plowing Moldboard 0.94716 − 1.30326E-003 × Soil moisture + 7.44398E-003 × Depth + 1.79299E-004 × Soil moisture × Depth
Chisel 1.08638 + 4.07646E-004 × Soil moisture + 2.63842E-003 × Depth + 1.79299E-004 × Soil moisture × Depth 
Disc 1.05377 − 2.40310E-003 × Soil moisture + 1.17546E-003 × Depth + 1.79299E-004 × Soil moisture × Depth
Start of the growing season Moldboard 0.63644 − 8.87002E-004 × Soil moisture + 0.019731 × Depth + 1.79299E-004 × Soil moisture × Depth
Chisel 0.76336 + 8.23900E-004 × Soil moisture + 0.014925 × Depth + 1.79299E-004 × Soil moisture × Depth
Disc 0.73159 − 1.98685E-003 × Soil moisture + 0.013462 × Depth + 1.79299E-004 × Soil moisture × Depth
End of the growing season Moldboard 1.09289 − 4.52339E-004 × Soil moisture + 4.08287E-003 × Depth + 1.79299E-004 × Soil moisture × Depth
Chisel 1.18898 + 1.25856E-003 × Soil moisture − 7.22686E-004 × Depth + 1.79299E-004 × Soil moisture × Depth
Disc 1.19387 − 1.55219E-003 × Soil moisture − 2.18565E-003 × Depth + 1.79299E-004 × Soil moisture × Depth
Table 3. Equations for predicting bulk soil density depending on the type of plow and growing season periods

Soil penetration resistance

The results of the analysis of variance given in Table 4 display a significant effect of soil moisture on soil resistance to penetration. As shown in Figure 7, the soil recorded the lowest resistance to penetration at 14% moisture, reaching 1133 kN m-2. The resistance increased at soil moistures of 7%, 22%, and 28%, reaching 1257, 1294, and 1379 kN m-2, respectively. The decrease in resistance at 14% moisture is due to the decrease in soil strength and resistance as a result of reducing molecular cohesion and cohesion of water films in the brittle state of the soil at this moisture, which makes the cohesion between soil particles weak and easy to disintegrate and penetrate. In contrast, resistance increases at 7% moisture due to the increase in molecular cohesion, which enhances the strength and resistance of the soil to penetration. As for moistures of 22% and 28%, the increase in resistance is due to the increase in cohesion resulting from water films and soil pressure resulting from the overlap of its particles and the blockage of pores, which increases the soil resistance to penetration. These results are consistent with those of Ahmadi and Mollazade (2009), who found that soil moisture between 13% and 15% reduced soil resistance to penetration by 40%.

Source Sum of squares df F-Value p-value (Prob > F)
Model 6.177E+007 19 44.37 < 0.0001
A-Moisture content 7.501E+005 1 10.24 0.0016
B- Depth 3.032E+007 1 413.79 < 0.0001
C-Growing season 1.686E+007 2 115.04 < 0.0001
D-Plow type 4.895E+006 2 33.41 < 0.0001
AB 12585.54 1 0.17 0.6790
AC 3.117E+006 2 21.27 < 0.0001
AD 1.875E+005 2 1.28 0.2804
BC 4.614E+006 2 31.49 < 0.0001
BD 3.113E+005 2 2.12 0.1223
CD 7.023E+005 4 2.40 0.0517
Table 4. Analysis of variance for the effect of study factors on soil resistance to penetration

Fig. 7. Effect of soil moisture on soil resistance to penetration (kN m-2)

The analysis results given in Table 4 illustrate a significant effect of soil depth on soil penetration resistance. As shown in Figure 8, soil penetration resistance increases with increasing soil depth from 15 to 30 cm, where the resistance reached 891 and 1641 kN m-2, respectively. This is attributed to the increase in soil strength and cohesion with depth, in addition to the effect of smoothing and root spread processes at a depth of 15 cm, which reduces soil density and thus reduces its penetration resistance. These results are consistent with the findings of Aminet al. (2014), who found that soil penetration resistance increases with increasing soil depth.

Fig. 8. Effect of soil depth on soil resistance to penetration (kN m-2)

The effect of crop growth period on soil penetration resistance is significant, as given in Table 4 and Figure 9. The results show that the lowest penetration resistance is recorded after the tillage process, reaching 897 kN m-2. As the growth period progressed, the resistance increased at the beginning and end of the season, reaching 1327 and 1573 kN m-2, respectively.

Fig. 9. Effect of growth periods on soil resistance to penetration (kN m-2)

The decrease in resistance after tillage is attributed to soil disintegration, increased porosity, and decreased density, which reduces its resistance to penetration. However, after planting and irrigation, wetting and drying increased soil density, soil aggregates were broken, and pores were clogged, resulting in increased soil penetration resistance during the growing season. These results are consistent with the findings of Martinset al. (2021), who observed an increase in soil penetration resistance at the end of the growing season compared to the beginning.

The results shown in Figure 10 and Table 4 display a significant effect of the type of plow on soil penetration resistance. It is found that plowing with a disc plow under field conditions recorded the lowest penetration resistance, reaching 1074 kN m-2. It is followed by plowing with a moldboard plow, which recorded a penetration resistance of 1282 kN m-2, while plowing with a chisel plow recorded the highest penetration resistance, reaching 1442 kN m-2. This is attributed to the fact that the disc plow contributed to reducing the bulk soil density due to its efficiency in working under field conditions compared to the moldboard plow and chisel plow. These results are consistent with what was indicated by Dekematiet al. (2019) and Boydas and Turgut (2007).

Fig. 10. Effect of plow type on soil penetration resistance (kN m-2)

Table 4 shows that the interaction between soil moisture and soil depth, the interaction between soil moisture and plow type, the interaction between soil depth and plow type, and the interaction between growth periods and plow type, do not have a significant effect on soil penetration resistance. However, the table shows a significant effect to the interaction between soil moisture and growth periods. As shown in Figure 11, the lowest penetration resistance is recorded at soil moisture 28% after plowing, reaching 728 kN m-2, which is attributed to the high moisture content after plowing, as soil penetration resistance is inversely affected by moisture at the time of work. In contrast, the highest penetration resistance is recorded at soil moisture 28% at the end of the growing season, reaching 1871 kN m-2.

Fig. 11. Effect of interaction between soil moisture and growth periods on soil resistance to penetration (kN m-2)

The results of the statistical analysis in Table 4 also show that there is a clear effect to the interaction between soil depth and growth periods. It is noted from Figure 12 that the 15 cm depth treatment after plowing recorded the lowest soil penetration resistance, reaching 454 kN m-2, while the 30 cm depth at the end of the growing season gave the highest penetration resistance, reaching 2083 kN m-2.

Fig. 12. Effect of interaction between soil moisture and growth periods on soil resistance to penetration (kN m-2)

Table 5 shows the mathematical models for each plow during the crop growth periods to predict the soil resistance to penetration under different field conditions. Through these equations, it is possible to predict the soil resistance to penetration by entering the variables of soil moisture and soil depth.

Measurement time Plow type Soil penetration equation
After plowing Moldboard − 170.41239 − 19.80071 × Soil moisture + 59.81069 × Depth +0.12804 × Soil moisture × Depth
Chisel 1.41900 − 18.80060 × Soil moisture + 60.92180 × Depth + 0.12804 × Soil moisture × Depth
Disc − 237.80372 − 11.48524 × Soil moisture + 49.67180 × Depth + 0.12804 × Soil moisture × Depth
Start of the growing season Moldboard 585.87333 + 9.89211 × Soil moisture + 23.60699 × Depth + 0.12804 × Soil moisture × Depth
Chisel 624.37138 + 10.89222 × Soil moisture + 24.71810 × Depth + 0.12804 × Soil moisture × Depth
Disc 535.14867 + 18.20759 × Soil moisture + 13.46810 × Depth + 0.12804 × Soil moisture × Depth
End of the growing season Moldboard − 220.63425 + 14.19066 × Soil moisture + 68.65329 × Depth + 0.12804 × Soil moisture × Depth
Chisel − 80.05287 + 15.19077 × Soil moisture + 69.76440 × Depth + 0.12804 × Soil moisture × Depth
Disc − 484.90058 + 22.50614 × Soil moisture + 58.51440 × Depth + 0.12804 × Soil moisture × Depth
Table 5. Equations for predicting soil resistance to penetration depending on the type of plow and the growing season periods

Conclusion

The study concludes that the use of smart computing programs such as Design Expert shows a high ability to predict the bulk density and penetration resistance of soil with great accuracy, as the coefficient of determination (R2) reached 0.8460 for the bulk density and 0.8114 for the penetration resistance, indicating the efficiency of mathematical models in predicting soil properties compared to field results. The results show that soil moisture at 14% recorded the lowest bulk density and penetration resistance, reaching 1.12 Mg m-3 and 1133 kN m-2, respectively, followed by soil moisture at 7%, then 22% and 28%. The disc plow also outperformed in reducing the bulk density and penetration resistance, recording 1.12 Mg m-3 and 1074 kN m-2, followed by the moldboard and then the chisel. The results indicate that increasing the soil depth leads to an increase in the bulk density and penetration resistance by 12% and 45.70% when moving from a depth of 15 cm to 30 cm. It also shows that the beginning of the growing season is associated with the lowest bulk density of 1.13 Mg m-3, followed by after tillage and end of season. While the lowest penetration resistance is recorded after tillage, reaching 897 kN m-2, followed by the beginning of the season and end of season.

It is recommended that further studies be conducted on soils of different textures, under different climatic conditions, and for other crops to predict changes in soil properties during the growing season.

Acknowledgments

We extend our sincere thanks and gratitude to the College of Agriculture at the University of Basra, the Department of Agricultural Machines and Equipment, and the Department of Soil Sciences and Water Resources for providing support and research requirements.

Conflict of Interest: The authors declare no competing interests.

Author Contributions

M. Almoosa: Conceptualization, Data acquisition, Data pre and post-processing, Validation, Text mining, Review and editing services.

S. Al-Atab: Supervision, Methodology, Technical advice.

S. Almaliki: Supervision, Statistical analysis, Numerical/computer simulation, Software services, Visualization.

References

  1. AbdulSada, A. J., and Almaliki, S. (2023). Prediction of Soil Compaction using Conventional Tillage Systems under Different Operating Conditions. In IOP Conference Series: Earth and Environmental Science (Vol. 1259, No. 1, p. 012127). IOP Publishing.DOI
  2. Ahmadi, H., and Mollazade, K. (2009). Effect of plowing depth and soil moisture content on reduced secondary tillage. Agricultural Engineering International: The CIGR EJournal, 11, 1-9.
  3. Ahmadi, I., and Ghaur, H. (2015). Effects of soil moisture content and tractor wheeling intensity on traffic-induced soil compaction. Journal of Central European Agriculture, 16(4), 489-502.DOI
  4. Almaliki, S., Himoud, M., and Al-Khafajie, A. (2019). Artificial neural network and stepwise approach for predicting tractive efficiency of the tractor (CASE JX75T). The Iraqi Journal of Agricultural Science, 50, 1008-1017.DOI
  5. Amin, M., Khan, M. J., Jan, M. T., Rehman, M. U., Tariq, J. A., Hanif, M., and Shah, Z. (2014). Effect of different tillage practices on soil physical properties under wheat in semi-arid environment. Soil and Environment, 33(1), 33-37. https://www.cabidigitallibrary.org/doi/full/10.5555/20143226515
  6. American Society of Agricultural and Biological Engineers. (2009). ASAE D497.6: Agricultural machinery management data. ASABE.
  7. Biberdzic, M., Barac, S., Lalevic, D., Djikic, A., Prodanovic, D., and Rajicic, V. (2020). Influence of soil tillage system on soil compaction and winter wheat yield. Chilean Journal of Agricultural Research, 80(1), 80-89.DOI
  8. Black, C. A., Evans, D. D., White, J. L., Ensminger, L. E., and Clark, F. E. (Eds.). (1965). Methods of soil analysis (Agronomy Monograph No. 9, Parts I and II). American Society of Agronomy.
  9. Boydas, M. G., and Turgut, N. )2007(. Effect of tillage implements and operating speeds on soil physical properties and wheat emergence. Turkish Journal of Agriculture and Forestry, 31, 399-412. https://journals.tubitak.gov.tr/agriculture/vol31/iss6/6/
  10. Dekemati, I., Bogunovic, I., Kisic, I., Radics, Z., Szemők, A., and Birkás, M. (2019). The effects of tillage-induced soil disturbance on soil quality. Polish Journal of Environmental Studies, 28(5), 3665-3673.DOI
  11. Hajabbasi, M. A. (2010). Tillage effects on soil compactness and wheat root morphology. Journal of Agricultural Science and Technology, 3, 67-77. https://dor.isc.ac/dor/20.1001.1.16807073.2001.3.1.4.7
  12. Javadi, A., and Spoor, G. (2006). The effect of spacing in dual wheel arrangements on surface load support and soil compaction Journal of Agricultural Science and Technology, 8, 119-131. https://dor.isc.ac/dor/20.1001.1.16807073.2006.8.2.3.3
  13. Kahlon, M., Lal, R., and Varughese, M. (2013). Twenty-Two Years of Tillage and mulching impacts on soil physical characteristics and carbon sequestration in central Ohio. Soil and Tillage Research, 126, 151-158.DOI
  14. Kostić, M. M., Rakić, D. Z., Savin, L. Đ., Dedović, N. M., and Simikić, M. Đ. (2016). Application of an original soil tillage resistance sensor in spatial prediction of selected soil properties. Computers and Electronics in Agriculture, 127, 615-624.DOI
  15. Kuhwald, M., Blaschek, M., Minkler, R., Nazemtseva, Y., Schwanebeck, M., Winter, J., and Duttmann, R. (2016). Spatial analysis of long‐term effects of different tillage practices based on penetration resistance. Soil Use and Management, 32(2), 240-249.DOI
  16. Kuroyanagi, N., Kaneko, A., Watanabe, T., Fujita, A., and Odahara, K. (1997). Effect of long- term application of organic matters on upland field. (2) yield of upland crop and physical properties of soil. ( Fukuoka Agricultural Research Center, Chikushino, Fukuoka 818 Japan) Bull. Fukuoka Agriculture Research Center, 16, 63-66. https://cir.nii.ac.jp/crid/1571417124296650112
  17. Martins, R. N., Portes, M. F., e Moraes, H. M. F., Junior, M. R. F., Rosas, J. T. F., and Junior, W. D. A. O. (2021). Influence of tillage systems on soil physical properties, spectral response and yield of the bean crop. Remote Sensing Applications: Society and Environment, 22, 100517.DOI
  18. Naderi-Boldaji, M., Azimi-Nejadian, H., and Bahrami, M. (2024). A Finite Element Model of Soil-Stress Probe Interaction under a Moving Rigid Wheel. Journal of Agricultural Machinery, 14(1).DOI
  19. Nassir, A. J. (2018). Effect of moldboard plow types on soil physical properties under different soil moisture content and tractor speed. Basrah Journal of Agricultural Sciences, 31(1), 48-58.DOI
  20. Rashidi, M., Tabatabaeefar, A., Keyhani, A., and Attarnejad, R. (2007). Non-linear amodeling of pressure-sinkage behaviour in soils using the finite Element method. Journal of Agricultural Science and Technology, 9, 1-13. https://dor.isc.ac/dor/20.1001.1.16807073.2007.9.1.7.2
  21. Salim, A. E. A., Almaliki, S. A., and Nedawi, D. R. (2022). Smart Computing Techniques for Predicting Soil Compaction Criteria under Realistic Field Conditions. Basrah Journal of Agricultural Sciences, 35(1), 188-211.DOI
  22. Shabanpour, M., Fekri, S., Bagheri, I., Payman, S. H., and Rahimi-Ajdadi, F. (2022). Effects of tillage method and drainage management on some soil physical properties. Journal of Agricultural Sciences, 24-24.DOI
  23. Shittu, K., Oyedele, D., and Babatunde, K. (2017). The effects of moisture content at tillage on soil strength in maize production. Egyptian Journal of Basic and Applied Sciences, 4(2), 139-142.DOI
  24. Taghavifar, H., and Mardani, A. (2014). Applying a supervised ANN (artificial neural network) approach to the prognostication of driven wheel energy efficiency indices. Energy, 68, 651-657.DOI
  25. Tahmasebi, M., Gohari, M., Sharifi Malvajerdi, A., and Hedayatipour, A. (2023). Development and field evaluation of a variable-depth tillage tool based on a horizontal pneumatic sensor measurement. Journal of Agricultural Machinery, 13(1), 85-99.DOI

Authors retain the copyright. This is an open access article distributed under Creative Commons Attribution 4.0 International License (CC BY 4.0)

Send comment about this article
Enter Name.
Enter a valid email address.
Enter a vaid affiliation.
Enter comments (At leaset 10 words)
CAPTCHA Image
Enter Security Code Correctly.

  • Receive Date 29 September 2024
  • Revise Date 19 November 2024
  • Accept Date 03 December 2024
  • First Publish Date 07 July 2025