Understanding the Role of Humic Acids on Crop Performance and Soil Health

Humic acids (HA) are organic molecules that play essential roles in improving soil properties, plant growth, and agronomic parameters. The sources of HA include coal, lignite, soils, and organic materials. Humic acid-based products have been used in crop production in recent years to ensure the sustainability of agriculture production. Reviewed literature shows that HA can positively affect soil physical, chemical, and biological characteristics, including texture, structure, water holding capacity, cation exchange capacity, pH, soil carbon, enzymes, nitrogen cycling, and nutrient availability. This review highlights the relevance of HA on crop growth, plant hormone production, nutrient uptake and assimilation, yield, and protein synthesis. The effect of HA on soil properties and crops is influenced by the HA type, HA application rate, HA application mode, soil type, solubility, molecular size, and functional group. This review also identifies some knowledge gaps in HA studies. HA and its application rate have not been tested in field experiments under different crops in rotation, nitrogen fertilizer forms, sites and climatic conditions. Furthermore, HA chemical and molecular structures, their water and alkaline soluble fractions have not been tested under field experiments to evaluate their effects on crop yield, quality, and soil health. The relationship between soil-plant nutrient availability and plant nutrient uptake following HA application should also be further studied.


INTRODUCTION
Humic substances (HS) are remains of decomposed plant and animal materials such as lignin, tannins, cellulose, and cutins (Tan et al., 2000;Billingham, 2012;Hayes and Swift, 2020). High quantities of HS are present in the soil after incorporating harvested residues (Wiesler et al., 2016). Increased animal and biogas production have reduced the amount of harvested residues on most arable land, resulting in decreased HS in the soil. Over the past decades, researchers have attempted to replenish the decreased HS with external applications (Rose et al., 2014;Gerke, 2018). The external sources of HS are mostly commercially produced from soils, coal, lignite, and organic materials (Gollenbeek and Van Der Weide, 2020;Yang et al., 2021).
HS are classified as humic acids (HA), fulvic acids (FA), and humins based on their solubility in water, acidic or alkaline solutions (De Melo et al., 2016). Due to the non-degrading nature of the humin fraction in HS, researchers have focused on the HA and FA fractions because they are capable of improving soil fertility and health within short time frames. The HA and FA fractions of HS are chemically reactive and able to resist microbial reactions, thereby performing beneficial roles in soils and plants (Billingham, 2012). The ability of HA to withstand degradation for long periods and their amphiphilic properties enable them to form complex cations (Wood, 1996). HA fraction contains about 60% organic carbon (C), which plays an important role in the growth of soil microorganisms (Sible et al., 2021). In addition to C, HA also contain nitrogen (N), oxygen (O), hydrogen (H), and sulfur (S).
Humic acids play several important roles such as: increase soil physical and biochemical activities by improving structure, texture, water holding capacity (WHC), and microbial population (Nardi et al., 2017(Nardi et al., , 2021Fuentes et al., 2018;Shah et al., 2018); increase soil nutrients availability, especially micronutrients by chelating and co-transporting micronutrients to plants (Yang et al., 2021); reduce the transportation of toxic heavy metals by precipitating them, thus reducing toxic heavy metals intake by plants (Wu et al., 2017). Humic acids also increase crop growth by increasing plant growth promoting hormones such as auxin and cytokinin, which aid in stress resistance, nutrients metabolism, and photosynthesis (Billingham, 2012;Rose et al., 2014;Canellas et al., 2020;Laskosky et al., 2020;Nardi et al., 2021;van Tol de Castro et al., 2021). Some studies have also reported no effects on crop growth and soil health following HA application (Albiach et al., 2001;Bybordi and Ebrahimian, 2013;El-Bassiouny et al., 2014;Mukherjee et al., 2014;Kelapa and Banyuasin, 2016). Although high HA doses are associated with enhanced soil physical characteristics (Gollenbeek and Van Der Weide, 2020), their effects on soil chemical characteristics and crops are still uncertain (Rose et al., 2014). In the review by Rose et al. (2014), among the factors analyzed in mostly greenhouse experiments, HA source had significant effects on both root and shoot growth while application rate only significantly affected shoot growth. A review by De Melo et al. (2016) highlighted carboxylic (COOH) and phenolic (OH) groups as predominant HA features that are largely responsible for their functions in the soil. A recent review by Nardi et al. (2021) showed that HA chemical and molecular structures, sources, and application rates are critical for determining their effects on crops and soil. Importantly, HA application can have inconsistent results on yield, possibly due to the different HA biological origins (Sible et al., 2021).
In view of inconsistent results of HA application on crop agronomic performance due to differences in HA sources and experimental conditions, and the paucity of literature on field experiments compared to laboratory trials, this review sought to understand HA application in agricultural production. The objectives of this review are to (1) identify the effects of HA on crop agronomic performance and soil health parameters in both laboratory and field experiments; (2) identify the factors that affect the efficiency of HA; (3) identify knowledge gaps in HA application on crop performance and soil health.

RELATIONSHIP BETWEEN HUMIC ACIDS STRUCTURE AND FUNCTIONS
The functions of HA are associated with their structures, which are source dependent (Rupiasih, 2005;Garciá et al., 2016;García et al., 2019;Nardi et al., 2021;van Tol de Castro et al., 2021).
Although HA structure contains many functional groups, the most predominant are phenolic (OH), and carboxylic (COOH) groups (Figure 1) (Nardi et al., 2021). The COOH and OH functional groups are mainly responsible for HA functions such as improving soil physical and chemical properties as well as plant growth (Figure 1) (De Melo et al., 2016;Nardi et al., 2021). Dissociation of these functional groups creates polar and nonpolar ends, which are the hydrophilic and hydrophobic parts, respectively (Mirza et al., 2011); both ends play roles in the mechanisms that confer useful HA functions (Figure 2A). The hydrophilic end is primarily involved in chelating functions, while the hydrophobic end is connected with repelling purposes (Billingham, 2012). Once the OH and COOH groups dissociate, the polar end of the anionic part forms complexes with cationic metals through electrostatic bonding in the soil, thus retaining these metals in the soil ( Figure 2B). The hydrophilic part, which is also water-loving, forms micelle that increases soil WHC. On the other hand, the non-polar end repels water molecules reducing water infiltration and improving clay aggregate stability (Billingham, 2012). A recent study by van Tol de Castro et al. (2021) reports that the aromatic and aliphatic functional groups of HA were responsible for increasing N uptake and soluble sugars, which resulted in a corresponding yield increase in rice (Figure 1); meanwhile an earlier finding by Garciá et al. (2016) showed that HS aliphatic and aromatic functional groups stimulated root growth in rice seedlings.
Humic acids with low molecular weight (LMW) contain more phenolic and carboxylic functional groups than HA with high molecular weight (HMW) (De Melo et al., 2016). The chelating ability of HA has also been attributed to LMW, which is efficacious in altering the biochemical characteristics of the FIGURE 1 | Graphical representation of humic acid, its chemical and molecular constitutes and functions. soil, while HMW is efficient in improving the soil physical conditions (Yang and Antonietti, 2020). HA with HMW have also been found to stimulate plasma membrane H + ATPase, allowing LMW HA to co-transport nutrients and perform other biological activities in plants (Figure 1) (Nardi et al., 2021). Enhancement in root growth was also found when both HMW and LMW fractions of HA from vermicompost were applied to both Arabidopsis and maize seedlings (Canellas et al., 2010). A recent review by Nardi et al. (2021) identified that HA molecular size could alter organic acid release in the rhizosphere. A summary of HA chemical and molecular components and functions are presented in Figure 1. More research needs to be conducted to elucidate the interactions of HA structure and organic acids, especially in the rhizosphere, where HA have different interactions with different crops. Furthermore, studies need to be carried out on how the molecular size of HA affects the release of root exudates considering that HA stimulating effect on roots' plasma membrane has been reported (Nardi et al., 2021).

EFFECT OF HUMIC ACIDS ON SOIL AND PLANTS
Soil Texture, Structure, and Water Holding Capacity The continuous practice of tillage and planting on the same land every year has a negative effect on soil texture and structure. The application of HA has been reported to have positive effects on the texture and structure of degraded soils (Billingham, 2012;Yang et al., 2021). The effect of HA on soil properties is summarized in Table 1. Soil structural stability has been attributed to increased adsorption of HA onto clay surfaces (Chen et al., 2017). The addition of HA leads to the formation of chelates with cationic metals (Yamaguchi et al., 2004;Billingham, 2012). These metals act as a bridge between HA and clay surfaces, thereby forming complexes ( Figure 2B). For example, the application of bentonite-humic acid increased the macro-aggregates of degraded sandy soil in a 7-year continuous cropping with maize (Zhou et al., 2019). The application of potassium humate increased aggregate stability in both bulk loamy acidic and sodic soils in a controlled study (Imbufe et al., 2005). In a similar experiment conducted by Piccolo et al. (1997), coal-derived humic acid increased soil structure stability irrespective of soil type in controlled wetting/drying experiments for 24 h.
Lack of positive responses following HA application have also been reported. Albiach et al. (2001) found that continuous application of commercially-produced HA had no significant effect on soil aggregate stability after 5 years of continuous cropping on sandy-silty-loam soil. Similarly, the application of coal-produced HA did not improve soil aggregate stability after two growing cycles of corn (Mukherjee et al., 2014). In both aforementioned experiments, the amount of HA applied was not sufficient to cause a change in soil texture and structure, as suggested by the authors. Furthermore, the soil used in the aforementioned studies were neutral to alkaline soils, which could negatively affect bridging humic molecules in HA, as demonstrated in the review by Gerke (2018). The HA source influences soil texture and structure (Rose et al., 2014). Therefore, it will be imperative for researchers to test the efficiency of a single HA source without extrapolating the results to other sources.
Humic acids have also been reported to increase soil WHC (Billingham, 2012;Yang et al., 2021). The water-attracting hydrophilic part of HA and improved soil structure increase soil WHC. It has also been reported that the combined application of HA and FA has a high chance of forming colloids or humicclay complexes, resulting in increased WHC (Billingham, 2012). In a 7-year field experiment to test the effect of bentonite HA on WHC, the application of 30 Mg ha −1 significantly increased soil WHC and the increase was more pronounced after the fourth year of the experiment (Zhou et al., 2019). Application of HA has also been shown to increase compatible solutes such as proline and glycine betaine in plants, which is an

Soil Cation Exchange Capacity
The ability of soil to hold nutrients depends on how much cations it can retain. HA have been shown to increase the cation exchange capacity (CEC) of the soil ( Figure 2D) (Billingham, 2012). Yang et al. (2021) summarized the contribution of HA in increasing CEC as follows: (1) increasing adsorption of exchangeable cations by providing large surface area for inorganic colloids; (2) dissociation of COOH and OH groups to produce polar ends that form complex with cations; and (3) increase dissolution of soil minerals, which create large surface area for chemical reaction. In an incubation experiment in bulk soil to study the effect of 26 different HA produced from peat and coal, all the samples tested increased CEC of amended soil from 1 to 58% (Giannouli et al., 2009). Surprisingly, there was no linear relationship between the initial CEC of tested samples and the corresponding percentage increase in CEC of amended soils. This may indicate that the CEC of the source material does not directly translate into increased CEC of the amended soil, rather HA quality has an important role on soil CEC. In another study by Laskosky et al. (2020), the effect of humalite, peat, and biochar on CEC of Orthic Gray Luvisol degraded soil was tested in a pot experiment growing barley plants; the residual soil CEC was higher in humalite compared to biochar amended soil. There is no evidence in literature evaluating the effect of HA application on CEC in bulk soil vs. rhizosphere in a single experiment. Most of the research on HA effects on CEC have been conducted under controlled environments and in short-term studies. Therefore, research needs to be conducted to fill in the knowledge gap on HA effect on soil CEC under long term field experiments involving several crops.

Soil pH
Soil pH influences nutrient availability. The ability of HA to affect soil pH changes depends on the amount of carboxylic and phenolic functional groups they contains (Rupiasih, 2005). Few studies have reported changes in pH with respect to HA application in post-harvest soil analysis. With the limited number of studies, inconsistent results on the effects of HA on soil pH are reported. In a barley-grown pot experiment with different HA application rates (0-26.2 g/kg), post-harvest soil pH analysis showed a decreasing trend in pH with increased humalite rates (Laskosky et al., 2020). In a 3-year continuous field cropping experiment involving peanut, HA application did not significantly affect soil pH . In a hydroponics study involving wheat to study HA effect on pH buffering capacity, the application of HA failed to increase the pH buffering capacity of the nutrient solution with an initial pH of 5.3 (Mackowiak et al., 2001); the authors attributed this failure to low HA concentration applied. Pertusatti and Prado (2007) in a lab assessment showed that HA buffer pH changes at pH between 5.5 and 8. Overall, the effect of HA on pH is dependent on the experimental conditions, plant grown and HA source. Studies needs to be conducted to elucidate the ideal condition specific for every HA source to identify its effect on soil pH changes.

Soil Carbon and Enzymatic Activities
Soil C content can be a direct measure of soil health. Carbon emission resulting in climate change and its adverse effects on the environment has resulted in research focused on soil C sequestration. HS have been reported to be an important pool of short-and long-lived soil C (Rupiasih, 2005). Although HA and FA fractions are decomposable in nature, their break down is slow, resulting in a constant supply of C to the soil. In the review by Sible et al. (2021), HA contain a high amount of C compared to FA. This indicates that the application of HA may provide additional C to soil microorganisms, essential for soil biological activities. The supply of C to the soil after HA application depends on the decomposition rate, turnover, and residence time in the soil (Fontaine et al., 2007). HA are supposed to undergo further transformational changes into smaller molecular sizes after soil application (Grinhut et al., 2007). It has been reported that HA turnover is contingent upon how synergistically plants and microorganisms combine in offering the needed priming effect (Dungait et al., 2012). In addition to the interaction between plants and microorganisms, environmental factors such as pH, moisture, oxygen, and HS properties affect HA decomposition rate (Dungait et al., 2012).
There have been numerous studies on the effect of HA application on soil C storage. In a 3-month soil incubation study with 13 C labeling technique to study the effect of HA addition on C sequestration, up to 58% C added was retained (Spaccini et al., 2002). The increase in C sequestration was a function of the HA chemical properties; the higher the hydrophobicity of the HA material, the higher the increase in C sequestration. Using incubation studies, Gümüs and Seker (2015) also found that the application of HA increased organic C content in clayey soil, where improvement in soil C depends on the HA application rate. In another study involving brachiaria/livestock and pearl millet/no livestock rotational system on clayey soil in Brazil, it was found that management systems that had a higher HA fraction had higher C sequestration (Loss et al., 2013). The application of HA-rich vermicompost significantly increased C sequestration in sandy-loam soil compared with NPK and soil alone treatments, through the increase in microbial biomass C and population, after 60 days of pea grown indoors (Maji et al., 2017). Most studies report on HA effects on total carbon pool but not labile and recalcitrant C. To our knowledge, only one study has reported the effect of HA from organic waste on labile and recalcitrant C in humus-rich soil under field conditions (Hu et al., 2019). They showed that soil supplements high in HA negatively correlated with labile C but had a significant and positive correlation with recalcitrant C in the soil. The paucity of literature on HA effects in long term field studies involving crops in rotations and post-soil analysis of C pool fractions warrants further research (Olk et al., 2018).
Humic acid application has been reported to increase microbial population and activities (Maji et al., 2017;Li et al., 2019). In a 3-year continuous cropping with pea to investigate HA effects on soil enzymatic activity, application of 1,000 Kg ha −1 HA significantly increased urease, phosphatase and sucrase activity after 140 days of plant growth under greenhouse conditions . In the same vein, the application of 9,000 Kg ha −1 HA-rich vermicompost in a pot experiment with pea significantly increased urease activity after 12 days of plant growth (Maji et al., 2017). In both experiments, the microbial population and C/N ratio increased, resulting in increased urease activity. In contrast, Shen et al. (2020b) found in a soil incubation experiment that HA application from weathered coal inhibited urease activity. Similarly, urease activity was inhibited when HA produced from Leonardite was investigated in soil incubation experiments (AL-Kanani et al., 1990). The microbial biomass was not determined in the previous two experiments by AL-Kanani et al. (1990) and Shen et al. (2020b), but less biological activity in bulk soil may have reduced microbial biomass and subsequent reduction in urease activity (Elmajdoub et al., 2014). Tomar and MacKenzie (1984) also suggested that HMW and enzyme binding by HA carboxylic and phenolic functional groups may inhibit urease activity. More research needs to be conducted to better understand HA effects on soil enzymatic activities involving different crops under field conditions in different soil types.

Urea Hydrolysis, Ammonification, and Nitrification
Urea is the most common N fertilizer used by crop producers due to its cost-effectiveness. For urea to be available for plant uptake, it has to be hydrolyzed into ammonium. Subsequent nitrification and de-nitrification convert ammonium into nitrate and nitrogen gases, respectively (Shen et al., 2020b). Unlike nitrate that is prone to leaching in the soil due to its hydrophobic nature, ammonium is stable in the soil (Oelmann et al., 2007). Plants can take up N in the form of ammonium and nitrate (Ups et al., 1990). The rapid conversion of ammonium to nitrate has been a major concern for agronomists as the latter has been shown to decrease nitrogen use efficiency (NUE) (AL-Kanani et al., 1990;Barth et al., 2020). The hydrolysis of urea is catalyzed by urease enzymes, which are produced by microorganisms in the soil (Tomar and MacKenzie, 1984). The rate of urea hydrolysis has also been linked to high concentrations of nickel (Tan et al., 2000). Due to the chelating properties of HA, it is able to form complex nickel, thus slowing down urea hydrolysis (Sible et al., 2021).
However, various studies have shown contradictory results on HA effects on urea hydrolysis (AL-Kanani et al., 1990;Maji et al., 2017;Li et al., 2019). In a 90-day soil incubation experiment to study the effect of coal-produced HA on urea transformation, HA application reduced urea hydrolysis by twofolds compared to the control (Shen et al., 2020b). In the same study, soil ammonium concentration was stable until the third day, after which it reduced drastically; the reduction in soil ammonium concentration reflected a corresponding increase in soil nitrate concentration. In a similar experiment by AL- Kanani et al. (1990), the application of up to 3.4% HA produced from Leonardite resulted in a 5-22% increase in urea hydrolysis in Typic Hapludoll and Typic Cryochrept soils. In all the aforementioned experiments, there was a decrease in pH in HAtreated soils. It has also been shown that HA application reduces soil nitrifying microorganisms, thereby increasing ammonium concentration and reducing soil nitrification rate (Dong et al., 2009). In a pot experiment grown with barley to study the effect of biochar, peat and humalite on urea hydrolysis and nitrification in the rhizosphere, humalite containing higher HA increased hydrolysis of urea into ammonium, while reducing nitrification of ammonium to nitrate (Laskosky et al., 2020). The high ammonium retention in the humalite treatment further resulted in increased CEC and decreased pH. It is important to note that HA used in different studies had different chemical and structural composition and application rates. Therefore, it will be interesting to investigate the effect of HA sources and application rate on urea hydrolysis and subsequent nitrification and denitrification in different crops under defined growth conditions.

Soil Nutrients Availability and Uptake by Plants
HA ability to increase soil nutrients availability and uptake by plants have been summarized in Figures 2A-F. Many studies have shown the ability of HA to stabilize ammonium, which improves soil N availability (Ahmed et al., 2006;Dong et al., 2009;Rose et al., 2014;Zhang et al., 2019;Laskosky et al., 2020;Shen et al., 2020b). Nitrogen is also present in HA molecules which becomes available to plants after soil application (Billingham, 2012). Plants mostly take up N in the form of inorganic ammonium and nitrate, but also N in the form of amino acids (Nardi et al., 2002). HA N uptake in higher plants depends on their molecular size and carboxylic group (Nardi et al., 2000). Piccolo et al. (1992) showed that N uptake was positively correlated with LMW and the amount of carboxylic functional group of HA. Tavares et al. (2019) studied the effect of vermicompost HA (80 mg L −1 ) on N uptake after exposing rice plants to ammonium and nitrate for 48 h; vermicompost significantly increased ammonium and nitrate uptake.
Apart from N, phosphorus (P) is also an essential nutrient, which synergistically increases crop growth and yield. The application of HA and FA has been shown to increase phosphatase activity by soil microorganisms, resulting in increased soil P solubilization (Sharma et al., 2013). HA also reduce sorption and increase desorption of soil phosphate ions, thereby increasing P in soil solution (Zhu et al., 2018). HA chelates soil micronutrients and co-transport them into plants ( Figure 2C) (Sible et al., 2021). However, HA ability to chelate micronutrients and co-transport them into plants depends on HA molecular weight (Zanin et al., 2019). HA also increase the permeability of the plasma membrane, which serves as the point of nutrient absorption ( Figure 2C) (Nardi et al., 2002). The higher the rate of HA applied, the higher CEC of the amended soil (Laskosky et al., 2020). While high doses of HA may improve the physical properties of the soil, its high binding capacities can also render some nutrients, especially micronutrients, unavailable for plant uptake. HA carboxylic and phenolic groups act as metalchelating agents that form metal-humic complexes in the soil, which may increase the presence of soil micronutrients but decrease their availability for plant uptake (Yang et al., 2021). For example, it is suggested by Shen et al. (2020a) that high doses of HA bind strongly to heavy metals such as iron, zinc and manganese in the soil, thereby limiting them for plant uptake. The application of potassium humate (up to 20 kg ha −1 ) in a rice field experiment increased soil micronutrients content, but these micronutrients were not available for plant uptake (Nandakumar et al., 2004). In another study, the application of 10 ml/L HA produced from Leonardite plus 30 ml Hoagland solution, did not significantly increase micronutrients and calcium uptake in broad bean (Bulut and Akinci, 2010). HA effects on soil micronutrients availability and plant uptake cannot be assumed to be a positive linear relationship due to differences in HA chemical and physical composition. Therefore, research needs to be conducted to elucidate the relationship between HA rates, nutrient availability, and plant uptake.

Nitrogen Assimilation and Protein Content in Plants
Nitrogen is very important nutrient for plants growth, yield and yield quality. Plants' ability to take up N is the first critical step of N assimilation processes in plants. N metabolism in plants is dependent on enzymes, which reduce N forms taken up by plants into final assimilatory products. In a 2-week hydroponics study, the application of different HA concentrations (0, 1, 5 mg L −1 ) significantly increased enzymes involved in the reduction and assimilation processes of N in maize in a dose-dependent manner (Vaccaro et al., 2015). Surprisingly, the dose-dependent increase in these enzymes did not reflect in the leaf protein content of maize; the low HA level had significantly higher protein content. The genes encoding these enzymes are mediated by nitrate supply (Vaccaro et al., 2015). In a similar study involving an 8-week pot experiment with lettuce by Haghighi et al. (2012), the application of up to 1,000 mg L −1 HA increased nitrate concentration and nitrate reductase in a dose-dependent manner. Unlike the study by Vaccaro et al. (2015), Haghighi et al. (2012) showed increase in protein content that was positively correlated with HA application rates. Inconsistent results have been reported on the application of HA on grain protein content. The application of up to 400 g liquid HA without base fertilizer in a field experiment significantly increased millet protein concentration relative to untreated control (Saruhan et al., 2011); when the mode of application was compared, soil and leaf application had the highest protein concentration than seed application. However, foliar application of up to 400 mg L −1 HA did not have significant effect on millet protein content in field experiments (Shen et al., 2020a). In another 2-year field experiment in Iran to evaluate the effect of foliar application of HA and different urea levels, the treatments did not significantly increase wheat protein and gluten content (Nasiroleslami et al., 2021). Soil application of 1,000 kg ha −1 HA significantly increased peanut protein content in the first growing year compared with only urea treatment but did not have a significant effect on protein content in the last 2 years of field experiments . Literature suggests that HA positively affect crop N uptake irrespective of the amount and form of N applied, but the assimilation process depends on the N form. Reviewed literature also suggests that HA have different effects on protein content depending on the amount of HA applied, application mode, and crop type. The lack of a clear pattern on N assimilation and crop protein content may be a function of HA type on enzymes that are involved in the metabolic process but more research is warranted in order to get a clear picture.

Plant Agronomic Parameters
The effect of HA on crop agronomic parameters has been summarized in Table 2. Various studies have been conducted to evaluate HA effects on plant growth and agronomic parameters such as root and shoot growth, leaf chlorophyll content, and yield. HA have been shown to stimulate root and shoot growth by enhancing the production of plant growth-promoting hormones such as auxin and cytokinin, and metabolic enzymes (Rose et al., 2014;Olaetxea et al., 2020). The improved uptake of macro and micronutrients following HA application increase the leaf chlorophyll concentration, which positively affect shoot growth (Chen et al., 2004;Fan et al., 2014;Sible et al., 2021). The production of plant hormones and enzymes as well as the increase in root and shoot weight, chlorophyll content, and photosynthetic rate following HA application, have been shown to improve yields (Delfine et al., 2005;Bybordi and Ebrahimian, 2013). Atiyeh et al. (2002) found that, the application of HA-rich vermicompost at rates between 0.15 and 0.4 g kg −1 significantly increased the height and leaf area of cucumber seedlings; higher concentrations failed to increase the root and shoot dry weight. Also, the application of HA-rich vermicompost significantly increased the length/height and weight of roots and shoots of pea compared to NPK treatment alone (Maji et al., 2017). In a related study, HA application up to 40 ml L −1 did not significantly affect the "chlorophyll a" concentration but increased "chlorophyll b" concentration, fruit weight, and yield of pepper (Karakurt et al., 2009). Application of Leonarditederived HA in greenhouse experiments increased the root and shoot biomass of Gerbera and canola, respectively (Akinremi et al., 2000;Nikbakht et al., 2008). The application of HA derived from organic waste significantly increased agronomic parameters of chrysanthemum including the chlorophyll content, leaf area, and root and shoot dry weight (Fan et al., 2014). Arjumend et al. (2015) reported a significant effect of different HA rates and recommended NPK mineral fertilizer on shoot and root weight, chlorophyll content, thousand grain weight, and grain yield in wheat. In a field experiment, Khan et al. (2010) reported that different HA rates and NPK significantly increased wheat plant height, spike length, thousand grain weight, and total yield; a major increase was found at higher HA application rates and moderate NPK rates. Based on field experiments conducted by El-Bassiouny et al. (2014) and Mohammed et al. (2019), HA applied with NPK fertilizer significantly improved the shoot dry weight and chlorophyll content of wheat and Stevia rebaudiana; there was no significant difference in wheat plant height, spike length and yield in this trial. In a 2-year field experiment conducted in Iran to evaluate the effect of foliar application of HA and different urea levels, HA treatment significantly increased the spike number, biomass, and yield of wheat at a urea application rate of 150 Kg ha −1 compared to 225 Kg ha −1 (Nasiroleslami et al., 2021). In another 2-year field experiments carried out in Iran with canola, the interaction of different urea and HA-rich zeolite rates did not have an effect on the agronomic performance of canola, but urea and zeolite alone significantly increased agronomic parameters (Bybordi and Ebrahimian, 2013).
The inconsistent results observed in the aforementioned studies show that recommendation for HA use in improving crop agronomic performances can only be reliable after being tested under specific conditions. Literature suggests that combined application of HA and mineral fertilizers form complexes that slowly release nutrients and subsequent crop uptake but the interaction effect depends mostly on the HA source, application rate, and crop type (Rose et al., 2014). Therefore, it is imperative to elucidate for a particular HA source, the optimum application and N fertilizer rates on crop agronomic parameters in multiple crops under defined growth conditions.

FACTORS AFFECTING HUMIC ACIDS EFFICIENCY
Humic Acid Source HA effects on soils and crops depends on the HS source (Rose et al., 2014;Gollenbeek and Van Der Weide, 2020). The source of HA applied depends on various factors such as nutritional composition, mode of production, functional group composition, and intended purpose. Among five different HA sources that were analyzed for their effectiveness on crop agronomic parameters, they followed the decreasing order; compost from manure < compost from green waste < soil < brown coal < peat (reviewed by Rose et al., 2014). HA extracted from different organic materials have different bioactivity potentials (Martinez-balmori et al., 2014). It was found that commercially-produced HA were less effective than HA produced from waste materials (Arancon et al., 2006). Jindo et al. (2020) also found that HA produced from composted materials were efficient in increasing plant agronomic and physiological activities. However, Khan et al. (2018) found no significant differences on wheat yield between treatments with HA produced from plants and coal. Different sources of HA contain different nutrient compositions and chemical structures, which can influence their performance in the soil. García et al. (2019) found different quantities of functional groups in HA obtained from Elliot soil, peat, leonardite, Su Wanee river, and Hill soil. In a pot experiment conducted by Hamad and Tantawy (2018) with three different sources of HA extracted from clayey soil, podrite and compost, sorghum root, and shoot growth correlated with amount of aromatic, aliphatic and carboxyl functional groups present in the various HA sources. Furthermore, N uptake was positively correlated with the high amounts of HA carboxylic functional group (Piccolo et al., 1992;Nardi et al., 2000). However, when several HA sources (lignite, soil, compost, coal, and peat) were tested for their fungicidal function, HA sources with high aromatic functional group related negatively to fungicidal function (Wei et al., 2018). This is an indication that HA source selection should be geared at a particular objective. Laskosky et al. (2020) conducted a pot experiment using three different HA-containing sources with different chemical properties; humalite, and peat which had high initial concentration of N and P compared to biochar, resulted in barley plants with significantly increased N and P concentrations. There are limited studies that evaluate and compare the effect of different HA sources on crop agronomic parameters under laboratory and field conditions, warranting further research.

Application Rate
It has been suggested that HA application rates are most effective under severe stress conditions (Rose et al., 2014). The effectiveness of HA application rate is also contingent on the source and crop type being grown (Olk et al., 2018). Under waterstressed conditions, plants triggered water deficit responses such as proline production. Under water-stressed field conditions, the effects of different HA application rates had a significant effect on millet yield, but the increase was not dose-dependent (Shen et al., 2020a). Similarly, HA elicited catalase and proline activity of maize seedlings under water-stressed conditions, but the increase was not dependent on the rate applied . However, an earlier study by Lotfi et al. (2015) reported a HA dose-dependent increase in proline and catalase activity in rapeseed under water-stressed conditions. Under saltstressed condition, HA increased agronomic parameters (e.g., plant height, leaf area, stem diameter, chlorophyll content and yield), and proline content of bean plants (Taha and Osman, 2018). Under similar salt-stressed conditions, Yousif et al. (2020) found that, increase in agronomic parameters and proline content in maize correlated with increase in HA application rates. Mohammed et al. (2019) found an increase in the growth and agronomic parameters of Stevia rebaudiana in a field experiment under optimal soil water and salt (Na + , Ca 2+ , and Mg 2+ ) conditions; the increase in growth and agronomic parameters correlated with increased HA application rates.

Factors Conditions References
HA sources Type of material extracted from (soil, peat, coal, lignite, organic residues etc.), nutrients composition, chemical structure and molecular weight. Rose et al., 2014;Hamad and Tantawy, 2018;García et al., 2019;Laskosky et al., 2020 Application rate Stress conditions (saline, drought, acidic, heavy metal concentration), nutrients composition of HA, source of HA, type of crop, type of soil, and environmental conditions. Bybordi and Ebrahimian, 2013;Rose et al., 2014;Lotfi et al., 2015;Olk et al., 2018;Mohammed et al., 2019;Baía et al., 2020;Shen et al., 2020a;Yousif et al., 2020;Sible et al., 2021 Soil type Type of clay, sand, and adsorption capacity of HA Feng et al., 2005;Zhang et al., 2013;Chen et al., 2017;Singh et al., 2017;Nardi et al., 2021;Sarlaki et al., 2021 Solubility Water, alkaline and acid extractable fraction of HA Pinton et al., 1999;Schmidt and Santi, 2007;Liu et al., 2008;De Melo et al., 2016;Savy et al., 2017;Sible et al., 2021 Frontiers in Agronomy | www.frontiersin.org Bybordi and Ebrahimian (2013) did not find any significant difference in canola agronomic parameters following different HA application rates. In nutrient-rich HA, the rate applied will have an effect on soils and plants since an additional source of nutrients from HA will be added to the soil (Sible et al., 2021). Karakurt et al. (2009) found an increase in pepper yield after moderate HA application, but similar yields were observed between the untreated control and higher HA application rates. HA application rates is dependent on environmental and soil conditions (Table 3), source and composition, as well as crop type, making it difficult to predict its effect on different crops.

Soil Type
Soil type plays an important role in HA adsorption and decomposition. HA are efficient when retained in the soil after application without leaching (Chen et al., 2017). Sandy soils have large textures and poor structure and therefore have poor retention of external application of nutrients and other soil amendments (Sarlaki et al., 2021). Soil clay fraction, which plays a vital role in retaining HA, differs among soil types . The differences in the binding capacities of different clay minerals in soils affect the rate at which HA are adsorbed onto soil surfaces. Kaolinite, which is 1:1 clay interacts effectively with HA, thereby retaining them onto its surface (Al-Essa, 2019). Feng et al. (2005) and Chen et al. (2017) found increased adsorption of HA onto kaolinite surfaces compared to montmorillonite. Chen et al. (2017) also found that kaolinite adsorbed HA more than montmorillonite due to their contrasting physical and chemical characteristics. Another study by Zhang et al. (2013) evaluated the effects of HA on three clay minerals (kaolinite, montmorillonite, and illite); they found that the specific surface area of montmorillonite decreased, thereby decreasing CEC compared to other clay minerals. However, hydrogen bonding in montmorillonite increased, thereby increasing adsorption capacity for NH + 4 . The efficacy of HA depends on how they are able to adsorb onto clay surfaces; hence differences in clay fractions at different sites will affect HA function and subsequent soil properties and crop performance. For example, Tahir et al. (2011) found that HA application increased wheat agronomic parameters in non-calcareous relative to calcareous soils. Khan et al. (2018) found higher wheat spike weight and grain yield on clayey loam compared to sandy loam soil. The recent review by Nardi et al. (2021) showed that different soils have different effects on maize nitrate and ammonium uptake as well as N metabolism after HA application. Although Rose et al. (2014) had indicated that the soil type had little impact on HA performance, their sample size was not adequate to extrapolate the results to represent all soil types. Literature on the effects of soil types on HA efficiency have been summarized in Table 3.

Solubility
Humic acid solubility depends on the pH of the medium (MacCarthy et al., 1990). HA are partially soluble in water and alkaline medium but precipitate under very low pH (De Melo et al., 2016). Researchers claim that the extraction of HA with alkaline change their structure thereby making alkaline solutions inappropriate for studies (Kleber and Johnson, 2010), but Olk et al. (2019, and references therein) claim the opposite. HA can form complexes with soil cationic nutrients, and the solubility of these complexes can influence the release of these nutrients to crops (Sible et al., 2021). Therefore, the proportion of HA that is soluble in acid, alkaline, and water will affect the growth of crops. Pinton et al. (1999) found that watersoluble HA increased nitrate uptake by activating root plasma membrane H + ATPase. Treatment of maize root with nitrate and water-extractable HA increased nitrate uptake, resulting in upregulation of nitrate assimilation enzymes (Zanin et al., 2018). Savy et al. (2017) found that water-soluble HA extracted from giant reed increased gibberellin activity in watercress (Lepidium sativum L.) seedlings due to the presence of phenolic moieties. The application of water-soluble HA increased root surface area of Arabidopsis due to an increase in root hairs, cortical cells, and tangential walls in the endodermal cell layer (Schmidt and Santi, 2007).
Furthermore, the acid extractable fraction of HA can form a stable complex with cations in the soil (Liu et al., 2008;De Melo et al., 2016), thereby increasing soil nutrients availability and improving soil physical and chemical characteristics (Billingham, 2012). There are limited studies on the effect of extractable fractions of HA in laboratory experiments and no extensive research on the effect of alkali and water-extractable fraction of HA on the growth and agronomic performance of crops under field conditions. Therefore, more research is needed to address this knowledge gap.

RECENT ADVANCES IN HUMIC ACID RESEARCH Extraction Methods
Efficient extraction methods of HA are still a challenge to many scientists. HA produced from peat, soil, and lakes limit larger field applications due to their lower quantities (Yang and Antonietti, 2020). Most commercially-produced HA are from lignite and coal due to their more extensive deposits, especially in oil producing regions. Humic acids from these sources tend to have more carbon but less oxygen and nitrogen content as well as low carboxylic and phenolic functional groups, which mostly contribute to HA physical, chemical and biological activities (Fatima et al., 2021). Recent development in HA research has focused on efficient extraction methods that produce high HA yield with many crop growthpromoting functional groups (phenolic and carboxylic groups). Pre-treatment of lignite and coal with acid breaks down the complex compounds, thus reducing them to weaker acidic functional groups (Barhoumi et al., 2019). Nitric acid (HNO 3 ) has been used to increase the yield of HA produced from lignite and bituminous coal following alkaline extraction (Zara et al., 2017;Boral et al., 2021;Fatima et al., 2021). In addition to increased HA yield, they also found an increase in nitrogen and oxygen contents and functional groups such as carboxylic, phenolic and amine groups. Sabar et al. (2020) have also used HNO 3 , H 2 O 2 , and fungal strains to increase the HA yield, molecular weight, and aromatic functional group in coal-derived HA. The development of efficient extraction methods is keen on increasing HA functional groups, thereby increasing HA bioavailability.

Molecular Mechanisms of Humic Acids
Humic acids are known over the past decades for their immense contribution to plant growth, agronomic parameters, stress tolerance, and soil health. The molecular mechanisms behind these positive effects on plants were not known until recently, when researchers have begun elucidating how plants respond to HA at the molecular level. Shah et al. (2018) detailed how HA contribute to plants' molecular response to growth and stress conditions in their review. In the current review, we highlight a few of the most recent developments regarding plants' molecular responses to HA. Cha et al. (2020) showed that Arabidopsis plant grown in HA up-regulated heattolerant gene and Heat-Shock Proteins (HSP); when HSP was knockout, Arabidopsis failed to withstand heat stress. Recently, Zandonadi et al. (2019) identified plant growth promoting hormone, Alkamides in HA, which is responsible for overexpression of protein genes inducing cell division and cytokinin production in maize. The application of HA caused down-regulation of GRF gene, a drought susceptible gene in wheat, indicating that HA can serve as a signaling molecule and trigger wheat growth under drought conditions (Arslan et al., 2020). Treating maize roots with HA up-regulated putative VHS/GAT and 2-cys peroxiredoxin BAS1 genes that are involved in antioxidant functions and N assimilatory pathway, respectively (Nunes et al., 2019). Identifying and understanding these molecular responses triggered by HA will assist plant breeders to target these genes in breeding involving HA.

Knowledge Gaps and Future Research Needs
The global effort to reduce the amount of N fertilizers in food production systems requires the optimization of N fertilizer FIGURE 3 | Model showing known and unknown research that have been conducted under pot and field experiments on crop yield and quality, soil quality and research gaps that need to be elucidated. COOH, carboxyl group; OH, phenolic group; HMW, high molecular weight; LMW, low molecular weight; WHC, water holding capacity; CEC, cation exchange capacity. application rates in different crops, soil types, and under unpredictable climatic conditions. HS are a promising tool to further optimize fertilizer application and nitrogen use efficiency in crops. We have identified some knowledge gaps (Figure 3) that warrants further research. There are no studies that have been conducted to test a specific HA sources, their application rate and mode of application either in a series of pot or field experiments in different crops. Neither has the interaction between HA, N forms, sites, and climatic conditions on crop yield and quality as well as soil health and quality been fully understood. We also identified research gaps in the relationship between HA-soil-plant nutrient availability and plant uptake under different experimental conditions. Such knowledge will help agronomists and crop producers understand how HA interacts with different crops, forming the basis for planning sustainable cropping systems. This review identified that the water and alkaline soluble fractions as well as the chemical and molecular structures of HA have not been assessed under both pot and field experiments for their effects on crop yield/ quality and soil health. Furthermore, there is a lack of information on how HA chemical and molecular structures (carboxyl, phenolic, aliphatic, aromatic functional groups, HMW, and LMW compounds) affect yield and crop quality, soil quality parameters, soil nutrient availability, plant uptake, and rhizosphere root exudates. Research on these topics will benefit industries involved in HA production to focus on HA fraction and functional groups capable of conferring significant benefits to crops. Most of the laboratory and field trials on HA effects on protein concentration have been conducted by analyzing the grain protein concentration but data are still limited, preventing solid conclusions. Also, data on the mechanisms on how HA contribute to N assimilation and protein production are still lacking considering that most of the cereal and pulse crops cultivated are required to attain high protein contents. Therefore, understanding how HA increases protein concentration in crops will be crucial for HA industry players and crop producers.

CONCLUSIONS
This review has revealed that HA application has potential significant effects on crop agronomic performance and soil quality parameters. This review identified several factors that affect HA performance in crops and soils; the most influential of them is the HA source. The HA chemical and molecular structure, solubility, and other factors such as application rate, soil, and crop type also affect HA effects on crop performance. This review via the evaluation of both laboratory and field experiments, identified the effects of HA on crop agronomic performance and soil health parameters. Knowledge gaps in HA studies have been identified in this review. More research is needed to optimize the combined effect of different HA application rates and mineral fertilizers on crop performance and soil quality parameters under defined field conditions; more importantly, long-term studies involving different soil types, crops and weather patterns in warranted to truly exploit the benefits of HS.