Quantifying Mechanistic Detachment Parameters Due to Humic Acids in Biological Soil Crusts

: Humic acid (HA) is a material that could be used to decrease erosion and improve soil structure. It is also known that biological soil crusts (biocrusts) have a major role in soil stabilization, but the mechanism is not well understood in the presence of HA, especially with mechanistic soil detachment rate parameters ( b 0 and b 1 ) of the Wilson model, where b 0 is the dimensional soil detachment parameter and b 1 is the dimensional soil threshold parameter. Therefore, this study intends to (1) investigate the effect of different humic acid (HA) concentrations (0%, 4%, and 8%) on mechanistic soil detachment rate parameters ( b 0 and b 1 ,) in the Tigris Riversides of the Gheraiat region, Baghdad City, Iraq, of the crusted versus uncrusted soils using a small scale model of the jet erosion tests (mini-JET) at different curing periods (1 week, 2 weeks, and 3 weeks), and (2) examine the impact of HA on b 0 and b 1 parameters versus some soil characteristics, such as electrical conductivity, cation exchange capacity, and soil organic matter for uncrusted and crusted soils. Thirty-six undisturbed soil specimens (18 for crusted soils and 18 for uncrusted soils) were acquired from the Tigris Riverbank. On these specimens, the mini-JET was used to determine the mechanistic cohesive soil erodibility parameters b 0 and b 1 . The results showed that the value of b 0 decreased up to 60% with an increase in curing times for crusted soils until they reached their optimum values at 2 weeks. There was no consistent pattern for b 1 at different curing times. As the concentration of HA increased, the value of b 0 decreased up to 86% and 99% for crusted and uncrusted soils, respectively. HA signiﬁcantly improved electrical conductivity, exchange capacity, and organic matter in the soil and reduced soil erodibility. This study provides the beneﬁts of adding HA to the soils as a soil stabilizer using a low-cost technique, which is the JET instrument. varied curing periods. In the presence of humic acid, the reduction in parameter b 0 reached 86 and 99% for crusted and uncrusted soils, respectively, in comparison with those without humic acid. These results were conﬁrmed by statistical analysis using paired t -test ANOVA. The EC and CEC were improved by humic acid and even more in crusted soils in comparison to uncrusted soils. This research found that the SOM signiﬁcantly increased in soils containing humic acid. The soil erodibility decreased as EC, CEC, and SOM increased. This work shows how employing a JET instrument can reduce testing times, save energy, and provide a low-cost method for observing crusted soil stabilization in the presence of humic acids.


Introduction
Humic acid (HA) is a major constituent of humic substances, which are the key organic components of soil, peat, coal, upland streams, dystrophic lakes, and ocean water. The addition of humic acid to the soil leads to an increase in the absorption of nutrients by the plant, as HA acts as a medium for transferring nutrients from the soil to the plant, especially in the case of soil drought [1,2]. By improving the soil properties and increasing the soil's ability to hold water, HA increases the number of beneficial microbes in the soil, dismantles heavy soil granules, and enhances its physical, chemical, and biological qualities [3].
A key factor in stabilizing soil aggregates, providing nutrients, and maintaining the microbiological balance of the earth is soil organic matter (SOM) [4]. The dynamics of SOM are motivated by the addition of organic compounds that are transformed by a combination of physical, chemical, and biological factors as time progresses [5,6]. When SOM decomposes, its composition, structure, and functional groupings change. Humic acids are made up of multiple fractions of SOM [7]. Fulvic and humic acids, for example, play a role in the transport of cations between the organic−mineral complexes exchanged [7]. Fulvic and humic acids also have an impact on the soil's organic cation exchange capacity (CEC) [8][9][10][11]. HAs have a dark color, are hydrophilic, and stable, and their components can vary as follows: 40 to 60% carbon; 30 to 50% oxygen; 0 to 4% nitrogen; 2 to 4% hydrogen; 0 to 2% sulfur [12]. The HA compounds interact with ions to produce complexes with various stabilities and structural properties [13]. All of these characteristics can be used to identify the soil management stage in agricultural areas.
Humic acids have been discovered to be connected with aggregate particle size distribution and soil carbon preservation in order to improve particle cohesion and aggregate stability [4]. By studying the spectrum properties of different fractions of SOM, it can provide significant information on how to improve soil quality [14,15]. The use of species diversity in agricultural management substantially increases the quality of SOM, thus effectively influencing the relationship between soil aggregates [16]. As a result, the chemical composition of the soil changes, and as part of soil quality, HA increases [17]. One technique to decrease erosion is to improve the soil structure by modifying soil qualities such as porosity and structural stability, as well as infiltration and other factors that contribute to increase water absorption in the soil [18]. This method improves soil water absorption. As a result, this reduces surface runoff and prevents soil particles from being eroded [1,2]. HA provides the same benefits as solid organic aggregates and can be used in a variety of irrigation methods [19].
Furthermore, the activity of soil microorganisms is an indicator of soil health because microorganisms contribute to soil fertility in many ecosystems, produce various biologically active substances, and have a higher efficiency in the biosorption of heavy metals [14,17]. Biological soil crusts (biocrusts) are a type of microorganism that may be found all over the world. Biocrusts have a significant influence on the biotic content of soils. Biocrusts are located in the top layer of the soil and have the potential to bind the soil particles together physically and chemically [20]. Biocrusts include bacteria, microalgae, cyanobacteria, fungus, protozoa, viruses, and certain nematodes. Biocrusts serve a variety of ecosystem activities and can also serve as a bioindicator of soil quality. Green algae (Chlorophyte), blue-green algae (Cyanobacteria), yellow-green algae (Xanthophyte), and diatoms are all part of biocrusts [21].
Soil erodibility is a natural process in which the natural erosive agents such as water, wind, gravity, or human-induced artificial processes remove and transport soil. Water erosion in many parts of the world is a serious and ongoing environmental concern that causes decreases in soil production and other soil functions and damages the environment [22,23]. Soil erodibility, in particular, threatens regional and even global ecological safety patterns. The influence of humans on soil erodibility is mainly due to industrial mining, processing, construction, smelting, road construction, and farm reconstruction activities. All of these destroy vegetation cover, increase runoff, and cause erosion. Prevention of soil erodibility is therefore fundamental to natural resource management and conservation [24].
There are a number of methods for determining soil erodibility in the presence of biocrusts caused by surface runoff [20,[25][26][27][28]. Using an electronic micro-penetrometer, Drahorad and Henningsen [25] investigated the biocrust stability of two research sites. They came to the conclusion that the micro-penetrometer device can detect areas of altering crust stability in connection to biocrust properties. Xiao et al. [26] measured soil erodibility in manufactured crusts versus bare soils in the field using a V-shaped runoff collection area over an eight-year period. They discovered that artificial crusts grew entirely after 4 years and reduced surface runoff by twenty-seven percent in the first 3 years. In addition, Gao et al. [27] investigated the erodibility of two crusted soils in China using nozzles to simulate rainfall from a height of 16 m onto experimental boxes (containing biocrusts). They found that biocrusts enhanced SOM content and reduced soil detachment and soil loss. Gao et al. [28] used a flume to investigate soil resistance owing to runoff scouring on Land 2021, 10, 1180 3 of 22 eight levels of crusted versus uncrusted soils. They discovered that soils with more than 36% moss cover are completely protected from runoff erodibility.
All of the preceding methods were expensive and investigated soil stability or erodibility due to surface runoff in the presence of biocrusts. Hashim et al. [20] used a miniature version of the jet erosion test (mini-JET) device to assess soil erodibility parameters in a convenient manner. The JET device uses a jet of water in submerged conditions. This jet of water is generated by constant water pressure, and it impacts on the soil surface. The jet creates a scour hole on the soil surface by applying a certain shear force [29][30][31][32][33][34][35]. The mini-JET is one of the approaches for determining soil erosion parameters that is gaining popularity [29][30][31][32][33][34][35]. The United States Department of Agriculture's Research Service (USDA-ARS) in Stillwater, Oklahoma, developed the mini-JET instrument [29]. The mini-JET device is simple to set up and use, it uses less water, and weighs less compared to the other techniques. Hashim et al. [20] utilized mini-JET on the soil specimens of two crusted versus uncrusted soils, collected from the Gheraiat region of Iraq. The specimens were air dried at different water contents to derive the erodibility parameters of the linear model. They found a reduction in erodibility coefficient for both biocrusts, while an increase in erodibility coefficient was observed for uncrusted soil. Hashim et al. [20] found that JET is a useful instrument for simulating rainfall-induced surface runoff. In comparison to prior methods, the JET technique involves less effort, uses less testing time, and is a more economical technique. However, till now, no study has used the JET device to explore the influence of HA on soil erosion parameters of crusted versus uncrusted soils.
There are two approaches for calculating cohesive soil erosion rates: a linear model and a mechanistic soil detachment rate model (nonlinear model). Wilson [36,37] developed the nonlinear mechanistic detachment rate model in 1993, and the dimensional detachment parameter (b 0 ) and the dimensional threshold parameter (b 1 ) are soil erodibility parameters in his model. Any additional forces that effect soil detachment can be incorporated, or considered, in the Wilson model parameters (b 0 and b 1 ). The Wilson model was widely utilized by several researchers to determine the rate of soil erosion using the JET device [29][30][31][32][33][34][35]. The nonlinear model parameters (b 0 and b 1 ) have never been determined using a JET device in the presence of humic acids for crusted versus uncrusted soils at different curing durations.
The current research exclusively explored the influence of HA on mechanistic soil detachment parameters b 0 and b 1 in biocrust and uncrusted soils at the Gheraiat region of the Tigris riverbanks of northern Baghdad City, Iraq. The farmers in this area employ HA to improve soil nutrition, resulting in high-quality plant outputs. There has been no previous research in this region or/and in the world to investigate the effect of HA on biocrust and to assess mechanistic soil erodibility parameters (b 0 and b 1 ). The behavior of crusted and uncrusted soils in the presence of HA was examined in the current study using varied curing durations. As a result, the major purpose of this work was to investigate how varying humic acid concentrations (0, 4%, and 8%) affected the mechanistic soil detachment parameters b 0 and b 1 of crusted versus uncrusted soils using the mini-JET at various curing durations (1 week, 2 weeks, and 3 weeks). This research also investigated the influence of HA on b 0 and b 1 parameters versus some soil characteristics such as electrical conductivity (EC), CEC, and soil organic matter (SOM). Figure 1 shows the research area, which is located in the region of Gheraiat, Tigris Riversides, Baghdad City, Iraq. The region of Gheraiat is a semiarid region with a large number of biocrusts that spread down the riverside. The high vegetation cover, moderate temperatures (especially in the spring), and wet underlying soils are the key causes for the increased amount of algae developing in this location [20]. Land 2021, 10, x FOR PEER REVIEW 4 temperatures (especially in the spring), and wet underlying soils are the key cause the increased amount of algae developing in this location [20]. The cover of biocrusts increased as a result of this protective plant system, con erably improving the soil environment. Different microenvironments distinguish t crusts. Cyanobacteria algae (blue-green) and Chlorophyte (green) biocrusts are the main forms of biocrusts found throughout the different areas of the soil surface in research region. The study region is primarily made up of filamentous Chlorophyt The cover of biocrusts increased as a result of this protective plant system, considerably improving the soil environment. Different microenvironments distinguish these crusts. Cyanobacteria algae (blue-green) and Chlorophyte (green) biocrusts are the two main forms of biocrusts found throughout the different areas of the soil surface in the research region. The study region is primarily made up of filamentous Chlorophyte biocrusts, which served as the study's focus. Eighty five percent of the Chlorophyte biocrust is made up of Zygogonium ericetorum algae, 10% of Spirogyra algae, and 5% of Oedogonium algae. These algae were measured using an optical microscope. The green (Chlorophyte) biocrust was compared to an uncrusted control soil from the same region, which had no impact from vegetation or organic matter. The two types of soil specimens (crusted soil and uncrusted soil) were collected using polyvinyl chloride molds (PVC), as shown in Figure 2.

Study Area
ocrusts, which served as the study's focus. Eighty five percent of the Chlorophyte biocrust is made up of Zygogonium ericetorum algae, 10% of Spirogyra algae, and 5% of Oedogonium algae. These algae were measured using an optical microscope. The green (Chlorophyte) biocrust was compared to an uncrusted control soil from the same region, which had no impact from vegetation or organic matter. The two types of soil specimens (crusted soil and uncrusted soil) were collected using polyvinyl chloride molds (PVC), as shown in Figure 2.

Humic Acid and Sampling Collection
DELTA K-Humate, made in USA, a brand of humic acid, was used with different concentrations of 0.0%, 4.0%, and 8.0% to be added to the collected soil samples. Organic HA granules are soil conditioners that may enhance nutrient uptake. Agro humic organic is thoroughly mixed with water to produce liquid products. The resulting organic solution is compatible with many liquid fertilizers. Humic acid is 99% sub-bituminous coal and 1% other ingredients. HA is composed of carbon, nitrogen, oxygen, and a small amount of phosphorous and sulfur. The concentration of 0% was just distilled water without humic acid. The concentration of 4% was a mix of 40 g of humic granules dissolved in 1 L of distilled water. The concentration of 8% was a mix of 80 g of humic granules dissolved in 1 L of distilled water. For crusted and uncrusted soil samples, 100 mL of HA concentration was added and distributed on the soil surface for each treatment.
To evaluate the effect of humic acid on soils containing biological crusts and bare soils, specimens were acquired at a moderate temperature during the spring season, when biological crust formation is at its peak. Thirty-six undistributed soil specimens were collected, with eighteen samples containing biological crusts and eighteen samples containing no biological crusts (uncrusted soils). Before extraction, the undisturbed soil specimens were slightly moistened to ensure that the biocrusts remained intact and were less likely to break up throughout the sample gathering process (Figure 2a,c). PVC sample molds were used to collect the undisturbed soil specimens, which had a diameter of 10.1 cm and a height of 11.8 cm. Then, the undisturbed soil specimens were transferred within one hour to the laboratory in a nylon bag inside a closed case ( Figure 2). The samples were kept in a glass house to ensure the survival of living organisms. However, by the third week of curing time, the living organisms had died and the soil surface was completely dry.

Analysis of Soil Physical-Chemical Characteristics
Physical and chemical testing of collected soil specimens was conducted according to ASTM standards [38], as stated in Table 1. The soil was classified as loam soil by the United States Department of Agriculture (USDA). Sieve analysis and hydrometer testing were carried out in accordance with ASTM Standard D422. The ASTM D854 standard was utilized to compute the basic soil gravity. Moisture content and SOM testing were performed in line with ASTM standard D2974. A pH meter and a portable conductivity meter were used to determine the pH and EC, respectively, in a 1:1 soil:water ratio slurry. The volume of CO 2 produced by the interaction of hydrochloric acid with soil carbonate was used to calculate the amount of CaCO 3 . The CaSO 4 determination was applied according to [39]. Exchangeable cations (CEC) were determined by the BaCl 2 method [40]. In addition, X-ray fluorescence (XRF) analyses were conducted to discover the elements of tested specimens. The Spector X-LAB 2000 XRF device was utilized by preparing tablets (pellets) from samples of the study soil by taking 4 g from the soil with diameters of less than 75 um and adding 0.9 g of a filler of wax for the purpose of the consistency of the model. The sample was then mixed well. The sample was placed in a tablet press machine with the use of air pressure of 15 bar for one minute, after which it was placed inside the XRF machine according to the method used in the laboratory. In this research, the uncrusted and crusted soil samples were analyzed using XRF analysis for a case without humic acid after one week of curing time and a case with humic acid after three weeks of curing time.

Testing Procedure
The thirty-six undistributed soil specimens (18 for crusted and 18 for uncrusted soils) were arranged and tested as follows. At the time of acquisition, humic acids with concentration levels of 0%, 4%, and 8% were added to each of the two samples of each group (group of six crusted samples and group of six uncrusted soils), as shown in Figure 2b,c. The zero concentration presented the samples with no humic acid. Then, these samples were tested using a mini-JET device (Figure 2d) in the same manner and after one week of humic acid addition, in order to investigate the influence of humic acids on soil samples after one week of curing time. The same procedure was repeated for curing times of two weeks and three weeks. Soil samples were tested at different curing times to investigate the effect of humic acids on mechanistic soil erosion parameters (b 0 and b 1 ) for both types of soils (crusted and uncrusted soils) at different soil moisture levels. Three levels of humic acid concentration were used to find the best concentration that could be used to stabilize the soil.

Deriving of Wilson Model Parameters
Wilson [36,37] created a comprehensive framework for estimating the rate of erosion based on dislodging and stabilizing forces, as well as corresponding moment lengths for particle detachment. Particle separation, according to Wilson [36,37], occurs when the driving force exceeds the resisting moment. Wilson [36] proposed that the time averaged velocity may be used to estimate incipient motion. Wilson [36] used a probability framework for turbulent forces created by Einstein [41] and Partheniades [42] to find the particle detachment. Because there was insufficient information on some parameters of the Wilson model, Wilson [36,37] employed a calibration approach to obtain some of these unknown parameters. Wilson [36,37] proposed a single method for calculating the overall detachment rate by include the effect of particle sizes in the calibration methods indirectly. As a result, the dimensional parameters b 0 and b 1 are used to express the entire erosion or detachment rate (ε r ) as follows: where, b 0 represents the dimensional soil detachment parameter in g/m.s.N 0.5 ; b 1 represents the dimensional soil threshold parameter in Pa; τ is the boundary shear stress in the JET environment in Pa; R s is a particle density; R w is the density of water; G r = 2 3 is the spherical particle's geometry ratio; P e denotes the lower particle parameter's exposure; P n is a composite of particle and fluid properties; P dd is the detachment distance parameter equal to 2; e v is the coefficient of variation equal to 0.35; P ls is a dimensionless parameter that is affected by particle size, bed orientation, and slope; f c is the dimensionless parameter based on soil cohesion; P o is known as a velocity flow parameter in the Wilson model; and d is an equivalent particle diameter. Al-Madhhachi et al. [29] created a spreadsheet to calculate the parameters b 0 and b 1 using curve-fitting approaches to reduce the error in these functions when compared to JET erosion data. Al-Madhhachi et al. [29] presented the operation, methodology, parts, description, and function of a mini-JET apparatus. The JET erosion data are the scour depth readings versus different time intervals, with a total test length of 120 min for each test.

Statistical Analysis
The outcomes of computed statistical differences for mechanistic erosion parameters (b 0 and b 1 ) of the two groups (crusted and uncrusted soils) utilizing normality tests (Shapiro−Wilk) of paired t-tests for analysis of variance (ANOVA) were presented. The mean and standard deviation of parameters b 0 and b 1 at different humic acid concentrations were also reported using a paired t-test ANOVA. The parameters b 0 and b 1 of crusted and uncrusted soils were also tested using the Tukey test method. When compared to the ANOVA with a significance level of 0.05, this revealed a significant discrepancy. The normalized objective function (NOF) was also used to assess how the predicted Wilson model parameters matched the observed parameters and was expressed as follows: where, V i is the observed data, D i is the predicted data, M is the observation number, j is an iteration, and M av is the observed mean data [34].

Mechanistic Detachment Parameters in Presence of Humic Acid
The influence of curing times of 1, 2 and 3 weeks on the dimensional soil detachment parameter (b 0 ) is shown in Figure 3, where the soil samples were measured with different doses of humic acid. For a case without humic acid (Figure 3a) for uncrusted soil, the value b 0 increased as curing time increased, whereas for biocrusts, the parameter b 0 decreased until it reached its low peak after 2 weeks. The filaments of Chlorophyte, which were able to generate surface aggregates and improve soil stability by creating extracellular polysaccharides, were responsible for the decrease in b 0 values of biocrusts.
A reduction in values b 0 for crusted and uncrusted soils was observed in the presence of HA at 4% and 8% concentrations, as shown in Figure 3b,c, respectively. At varying curing times and in the presence of HA, there was no consistent behavior of b 0 . The findings in Figure 3 were confirmed by statistical analysis using paired t-test ANOVA, which compared the mechanistic detachment parameters (b 0 and b 1 ) of crusted and uncrusted soils at different humic acid concentrations (Table 2). In the absence of humic acid, there was a statistically significant difference in b 0 between uncrusted and crusted soils (p = 0.031). For the cases of humic acid (4% and 8%), the results indicated that there was no statistically significant difference of b 0 for uncrusted versus crusted soils (p > 0.050). Table 2. Statistical analysis of the influence of curing times (1 week, 2 weeks, and 3 weeks) using paired t-test ANOVA on mechanistic detachment parameters (b 0 and b 1 ) for undisturbed soil specimens (n = 36) of crusted and uncrusted soils at different humic acid concentrations (0%, 4%, and 8%).  Figure 4 shows the effect of curing durations of 1, 2, and 3 weeks on the dimensional soil threshold parameter (b 1 ), which was tested with varied concentrations of humic acid. In the absence of humic acid, there was no consistent behavior of b 1 at different curing times for crusted and uncrusted soils (Figure 4a). Hashim et al. [20] found similar results for critical shear stress in a linear model. An occasion with humic acid produced similar outcomes for parameter b 1 (Figure 4b,c). Table 2 confirms the findings of b 1 values in Figure 4. There was no statistically significant difference between crusted and uncrusted soils for humic acid, and no humic acid cases with a p-value of more than 0.05 (Table 2).  A reduction in values b0 for crusted and uncrusted soils was observed in the presence of HA at 4% and 8% concentrations, as shown in Figure 3b,c, respectively. At varying curing times and in the presence of HA, there was no consistent behavior of b0. The findings in Figure 3 were confirmed by statistical analysis using paired t-test ANOVA, which compared the mechanistic detachment parameters (b0 and b1) of crusted and uncrusted soils at different humic acid concentrations ( Table 2). In the absence of humic acid, there was a statistically significant difference in b0 between uncrusted and crusted soils (p = 0.031). For the cases of humic acid (4% and 8%), the results indicated that there     Figure 5 illustrates the behavior of b 0 for a case with no HA, a case with 4% HA, and a case with 8% HA, where the soil samples were measured at different curing times, in order to explore the influence of HA concentrations on parameter b 0 . For both crusted and uncrusted soils, and for all curing durations, Figure 5 indicates that the value b 0 declined as HA concentration increased. Statistical analysis using paired t-test ANOVA in Table 2 confirms the findings in Figure 5. There was no statistically significant difference in b 0 between uncrusted and crusted soils, except for the case with no humic acid.   The influence of HA concentrations on parameter b 1 , which were tested with different curing durations, are shown in Figure 6. There was no consistent behavior of b 1 for crusted and uncrusted soils at varied curing times in all cases ( Figure 6). As mentioned above, Table 2 confirms the findings of b 1 values in Figure 6.

Influnce of Humic Acids on Soil Charecteristics
In order to investigate the HA influence on soil characteristics, X-ray fluorescence (XRF) analyses were performed to determine the components of tested samples. Table 3 shows that the CaO was found in concentrations ranging from 16.24 to 14.76%, whereas SiO 2 was found in concentrations ranging from 39.27 to 34.36% as determined by XRF analysis. This is owing to the calcareous structure of the crust. The average content of Al 2 O 3 was roughly 8.62%, according to an estimate from XRF analysis. Fe 2 O 3 accounted for 5.25% of the total elements. The average MgO content was 4.37%. All of the samples analyzed had simple but substantial changes in their chemical composition (oxides). The CaO and SiO 2 concentrations were higher in samples with humic acid than in samples without humic acid, especially with crusted soils. In order to examine the impact of HA on erosion parameters b 0 and b 1 versus some soil characteristics, the influence of HA was first investigated on several soil characteristics in order to determine which soil characteristics would be chosen. With crusted soils, the results showed that there was a slight decrease in soil pH of 7.37, 7.10, and 6.98 for a 4% concentration of HA, and 7.09, 6.78, and 6.32 for an 8% concentration of HA, as compared with soil pH with no humic acids of 7.51, 7.49, and 7.45 for one week, two weeks, and three weeks of curing times, respectively (Figure 7b). Different results were observed for uncrusted soils (Figure 7a). Various functional groups in humic acid, which neutralize the soil pH, caused a slight increase in pH. Similar results were reported by SÖNMEZ anḑ Sevket [43] and Çelik et al. [44] for uncrusted soils.
With crusted soil, as shown in Figure 7d, EC(electrical conductivity) values increased significantly in comparison to a high HA application rate. According to a 4% dose of HA, the EC was 1.8, 2.32, and 2.98 d.s/m and for an 8% dose of HA, the EC was 2.5, 3.1, and 3.42 d.s/m as compared with the control sample of 1.5, 1.9, and 2.4 d.s/m for 1 week, 2 weeks, and 3 weeks, respectively. The same results were observed for uncrusted soils, revealing that soil EC linearly increased in response to increments in HA dose. Even higher values of EC for crusted soils in comparison to uncrusted soils (Figure 7c,d). As a result of the present and previous studies, the increase in HA dose is associated with an increase in soil EC level. Experiments with varied dosages of HA may have higher EC values because the HA materials added have some counter ions in solution (such as K + or Na + ) along with negatively charged HA molecules that increase EC organic particles [45].
There is a significant calcium carbonate content (CaCO 3 ) in the soil. The results of the crusted soil in Figure 8b demonstrate a slight decrease in CaCO 3 content as the percentage of humic acids increases. After one week, two weeks, and three weeks at 4% HA, the lime content was 229.12, 227.94, and 225.38 g/kg of soil, respectively, and 228.3.1, 226.09, and 222.17 g/kg of soil, respectively, at 8% HA. For soils with no humic acid, the CaCO 3 content was 231.23, 241.24, and 248.6 g/kg of soil at one week, two weeks, and three weeks, respectively. This reduction can be explained by the humic acid's release of carbon dioxide (CO 2 ) from calcium carbonate. As a result, carbon dioxide in the soil solution reacts with water to generate carbonic acid, which subsequently separates into its ions, with the hydrogen ions reducing the soil pH as proposed by SÖNMEZ andŞevket [43]. SÖNMEZ andŞevket [43] gave the same results. For uncrusted soil, there was no influence either of humic acid or of curing times on CaCO 3 (Figure 8a). Another key chemical component in soils is calcium sulfate (CaSO 4 ). Both crusted and uncrusted soils showed a decrease in CaSO 4 as the humic acid dose increased (Figure 8c,d). In comparison to uncrusted soils, crusted soils had an even higher CaSO 4 content. the results showed that there was a slight decrease in soil pH of 7.37, 7.10, and 6.98 for a 4% concentration of HA, and 7.09, 6.78, and 6.32 for an 8% concentration of HA, as compared with soil pH with no humic acids of 7.51, 7.49, and 7.45 for one week, two weeks, and three weeks of curing times, respectively (Figure 7b). Different results were observed for uncrusted soils (Figure 7a). Various functional groups in humic acid, which neutralize the soil pH, caused a slight increase in pH. Similar results were reported by SÖNMEZ and Şevket [43] and Çelik et al. [44] for uncrusted soils.   Figure 9a,b indicate that soil CEC increased as humic acid increased, reaching a peak after three weeks. At three weeks of curing time, the CEC of the control soil (soil with no humic acid) was 3.48 cmol/kg in crusted soil (Figure 9b). The CEC for a 4% HA treatment was 4.46 cmol/kg, and for an 8% HA concentration, it was 5.64 cmol/kg. In comparison to the control soil, HA treatments considerably enhanced CEC. The same results were observed with crusted samples (Figure 9a).
The different concentration treatments of HA significantly affected the soil organic matter content (SOM) based on error bars of standard deviation (Figure 9c,d). The results of the crusted soils (Figure 9d) showed that SOM increased with the concentration of HA. The values of SOM were 7.2%, 8.4%, and 9.6% for the 4% dose of HA, 9.7%, 10%, and 10.76% for the 8% dose of HA, and 2.56%, 2.87%, and 3.12% for no humic acid soil at 1 week, 2 weeks, and 3 weeks, respectively. The same results were obtained for samples without a crust (Figure 9c). It is interesting to note that the level of SOM was the same for uncrusted soil in the case with no HA in comparison to crusted soil. This was due to the absence of biological shell and humic acids. The different concentration treatments of HA significantly affected the soil organic matter content (SOM) based on error bars of standard deviation (Figure 9c,d). The results of the crusted soils (Figure 9d) showed that SOM increased with the concentration of HA. The values of SOM were 7.2%, 8.4%, and 9.6% for the 4% dose of HA, 9.7%, 10%, and 10.76% for the 8% dose of HA, and 2.56%, 2.87%, and 3.12% for no humic acid soil at 1 week, 2 weeks, and 3 weeks, respectively. The same results were obtained for samples without a crust (Figure 9c). It is interesting to note that the level of SOM was the same for The results from Figures 7-9 indicated that HA had an influence on EC, CEC, and SOM for crusted and uncrusted soils. Accordingly, Figure 10 presents the relationships between mechanistic soil erodibility parameters (b 0 and b 1 ) and EC, CEC, and SOM, for crusted and uncrusted soil samples with and without humic acids. Parameter b 0 decreased as EC increased, while there was no general trend for parameter b 1 (Figure 10a). The EC was 0.9 d.S/m initially at one week of curing time for uncrusted soil (without HA) and reached 4.9 d.S/m at the third week of curing time for crusted soil (with 8% of HA). Similarly, parameter b 0 decreased as CEC and SOM increased, while there was no general trend for parameter b 1 (Figure 10b,c). The CEC and SOM were 2.41 cmol/kg and 0.7% initially at one week of curing time for uncrusted soil (without HA), and reached 5.64 cmol/kg and 10.5% at the third week of curing time for crusted soil (with 8% of HA), respectively. uncrusted soil in the case with no HA in comparison to crusted soil. This was due to the absence of biological shell and humic acids.
The results from Figures 7-9 indicated that HA had an influence on EC, CEC, and SOM for crusted and uncrusted soils. Accordingly, Figure 10 presents the relationships between mechanistic soil erodibility parameters (b0 and b1) and EC, CEC, and SOM, for crusted and uncrusted soil samples with and without humic acids. Parameter b0 decreased as EC increased, while there was no general trend for parameter b1 (Figure 10a). The EC was 0.9 d.S/m initially at one week of curing time for uncrusted soil (without HA) and reached 4.9 d.S/m at the third week of curing time for crusted soil (with 8% of HA). Similarly, parameter b0 decreased as CEC and SOM increased, while there was no general trend for parameter b1 (Figure 10b,c). The CEC and SOM were 2.41 cmol/kg and 0.7% initially at one week of curing time for uncrusted soil (without HA), and reached 5.64 cmol/kg and 10.5% at the third week of curing time for crusted soil (with 8% of HA), respectively.

Discussion
For crusted soils with no humic acid, the erodibility coefficient of the linear model from the Hashim et al. [20] study decreased as curing time increased. In the present research, the parameter b 0 of the nonlinear model (Wilson model) behaved similarly to the erodibility coefficient of the linear model until it reached its low peak at 2 weeks, after which the b 0 increased (Figure 3a). As a result, the b 0 of the nonlinear model behaved similarly to the erodibility coefficient of the linear model at different curing times. Hashim et al. [20] Land 2021, 10, 1180 18 of 22 conducted their experiments for a maximum of 12 days of the curing period. The current study, on the other hand, was carried out with a maximum curing time of three weeks. The increase in the parameter b 0 after two weeks is due to the breaking down of dead algae, which resulted in an increase in scour depth readings (for a single test by mini-JET) as shown in Figure 11.

Discussion
For crusted soils with no humic acid, the erodibility coefficient of the linear from the Hashim et al. [20] study decreased as curing time increased. In the pre search, the parameter b0 of the nonlinear model (Wilson model) behaved similarl erodibility coefficient of the linear model until it reached its low peak at 2 week which the b0 increased (Figure 3a). As a result, the b0 of the nonlinear model b similarly to the erodibility coefficient of the linear model at different curing Hashim et al. [20] conducted their experiments for a maximum of 12 days of the period. The current study, on the other hand, was carried out with a maximum time of three weeks. The increase in the parameter b0 after two weeks is due breaking down of dead algae, which resulted in an increase in scour depth reading single test by mini-JET) as shown in Figure 11.  When HA is applied to the soil, it has a considerable impact on soil physical ical, and biological features [46]. An increase in SOM content due to humic acid re increased soil aggregation, water-holding capacity, infiltration rate, and water re The effect of HA concentrations on b0 for crusted and uncrusted soils (in the study) revealed a reduction in b0 with increasing HA dose, with the best stability f Figure 11. Comparison of scour depth readings using the mini-JET device for curing times of 2 weeks and 3 weeks for crusted soils with no humic acid and how the Wilson model fits the observed data. Figure 11 shows a comparison of scour depth readings and how the Wilson model fits the observed data at 2 and 3 weeks of curing times. The NOF values were 0.242 for 3 weeks and 0.083 for 2 weeks.
When HA is applied to the soil, it has a considerable impact on soil physical, chemical, and biological features [46]. An increase in SOM content due to humic acid results in increased soil aggregation, water-holding capacity, infiltration rate, and water retention. The effect of HA concentrations on b 0 for crusted and uncrusted soils (in the current study) revealed a reduction in b 0 with increasing HA dose, with the best stability found at an HA concentration of 8% ( Figure 5). These results can be explained by humic acid's ability to accumulate soil aggregates, resulting in decreased soil erosion. These findings were consistent with those of Khaled and Fawy [47], who discovered that soil stability increased as the HA dose increased. Khaled and Fawy [47] claim that HA reduces soil erosion by enhancing the cohesive forces between soil particles. HA improves the exchange capacity and refining features of soil, as well as its structure and physical properties, to increase the absorption of several elements and make them available to plants. The humic acid improved aggregate stability by creating clay-metal-humic connections, which were facilitated by the O-containing, hydrophilic groups of HA, according to Piccolo and Mbagwu [48,49].
Gümüş andŞeker [50] investigated the stability of soil aggregates with varying HA doses (0, 0.5, 1, 2, and 4%) over three curing periods (21,42, and 62 days). They found that aggregate stability increased slowly with the HA curing period. The same results were obtained in our study for uncrusted soils (Figures 3 and 5). Incubation periods for different organic amendments were examined by Piccolo and Mbagwu [49] for their effect on aggregate stability. Adding HA to soil after 7 and 15 days resulted in the highest aggregate stability values [43]. Aggregate stability decreased by 79% for HA after 40 days of incubation compared to the 7-day values [49]. Soil efficacy improved after 7 and 15 days of HA treatment, but after 40 days there was no significant difference. As a result, the efficacy of HA alone peaked at 15 days.
In comparison to a high HA dose, Gümüş andŞeker [50] discovered that EC values increased dramatically. Our research yielded the same results, as shown in Figure 7c. Hashim et al. [20] found that the erodibility coefficient of a linear model decreased as salt concentration increased. Similar behavior for parameter b 0 was found in the present study ( Figure 10a). Drahorad and Henningsen [51] found that the contact points in the biocrusts were cemented between the surface and soil particles as salt concentration increased. According to Hansen et al. [52] and Yu et al. [53], soil CEC increases as humic acid increases, peaking after three weeks in soil samples without crusts. Figure 9a of the current investigation confirmed these findings. Interactions between SOM and clay minerals resulted in higher soil CEC values, according to Yu et al. [53]. The results of Figure 10b found that parameter b 0 decreased as CEC increased for crusted and uncrusted soils. Figure 9c,d showed that differing HA concentrations had a substantial impact on SOM for crusted and uncrusted soils. Humic acid treatments resulted in a considerable increase in SOM. These findings matched those of Sootahar et al. [54]. Gümüş andŞeker [50] discovered similar results, with SOM values increasing dramatically as the amount of HA applied increased. The SOM increased significantly after a 21-day curing period with 0.5% of HA [50]. At 42 and 62 days of curing, the soil SOM increased linearly when the HA dose was increased, with the highest effect at 2 and 4% concentrations of HA. Hashim et al. [20] found that SOM tended to increase in soils containing biocrusts. It was clearly indicated that this increment was due to the application of an algal crust. Hashim et al. [20] reported that the erodibility coefficient of a linear model decreased as SOM increased. Our research reported similar results for parameter b 0 related to the values of SOM (Figure 10c). There was a gap in the values of SOM from 2 to 5% related to the influence of HA on parameter b 0 and SOM. The soil samples of SOM below 2% are the crusted and uncrusted soil samples without humic acid, while the others above 5% of SOM are the crusted and uncrusted soil samples with humic acid (Figure 10c).

Conclusions
A mini-JET apparatus was employed on the undisturbed soil specimens in the presence of humic acids (concentrations of 0 to 8%) of crusted and uncrusted soils acquired from the Gheraiat region of Iraq, and air dried at different curing times (one to three weeks) to determine the mechanistic soil erosion parameters (b 0 and b 1 ) of the Wilson model. At two weeks of curing time, the parameter b 0 for crusted soils with no humic acid was reduced by 60%, whereas at three weeks, the parameter b 0 increased by 90% for uncrusted soil. The breaking down of dead algae caused a slight increase in b 0 after two weeks on the crusted soils. For crusted and uncrusted soils, there was no consistent behavior of b 1 at varied curing periods. In the presence of humic acid, the reduction in parameter b 0 reached 86 and 99% for crusted and uncrusted soils, respectively, in comparison with those without humic acid. These results were confirmed by statistical analysis using paired t-test ANOVA. The EC and CEC were improved by humic acid and even more in crusted soils in comparison to uncrusted soils. This research found that the SOM significantly increased in soils containing humic acid. The soil erodibility decreased as EC, CEC, and SOM increased.