Next Article in Journal
S-Like Ribonuclease T2 Genes Are Induced during Mobilisation of Nutrients in Cotyledons from Common Bean
Next Article in Special Issue
A Sustainable Approach for Improving Soil Properties and Reducing N2O Emissions Is Possible through Initial and Repeated Biochar Application
Previous Article in Journal
Effects of Mechanically Transplanting Methods and Planting Densities on Yield and Quality of Nanjing 2728 under Rice-Crayfish Continuous Production System
Previous Article in Special Issue
Quanti-Qualitative Response of Swiss Chard (Beta vulgaris L. var. cycla) to Soil Amendment with Biochar-Compost Mixtures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrological Properties of a Clay Loam Soil as Affected by Biochar Application in a Pot Experiment

Department of Science of Agriculture, Food, Natural Resources and Engineering, University of Foggia, 71122 Foggia, Italy
*
Author to whom correspondence should be addressed.
Submission received: 8 January 2021 / Revised: 21 February 2021 / Accepted: 2 March 2021 / Published: 5 March 2021

Abstract

:
Improving soil-water relations by amending soil with biochar might play a significant role in increasing water availability for agricultural crops as well as decreasing water loss through drainage or runoff. While the effects of biochar on the hydrological properties on coarse-textured soils are generally positive and well-documented in the literature, studies on biochar effects on fine-textured soils are rather scarce and even contradictory. Therefore, the aim of this paper was to investigate the impact of biochar on the bulk density, water retention curve (together with several water capacitive indicators) and water infiltration rate in a clay loam soil. A pot experiment was carried out under lab conditions in which biochar was mixed with soil at rates of 0 (B0 or control), 2, 4, 6, 8 and 10% dw (B2, B4, B6, B8 and B10, respectively). Water retention of soil–biochar mixtures at different matrix potentials was determined using a pressure plate apparatus. From these measurements, a series of capacitive indicators was derived and the fitting of the van Genuchten model was also performed. Water infiltration into soil–biochar mixtures was measured by means of a mini-disk infiltrometer and the obtained data were analyzed both directly and by fitting the Philip’s model. Biochar significantly affected the considered soil properties. As the biochar rate increased, the bulk density decreased and water retention increased (B6, B8 and B10 > B2, B4 and B0), while the infiltration rate decreased (B0 > B2, B4, B6, B8 and B10). Although the experiment was performed on sieved and repacked soil samples under lab conditions, the results confirmed that biochar has the potential to increase plant-available water, while possibly reducing drainage water in a clay loam soil by lowering the infiltration rate.

1. Introduction

Over the last decade, the use of biochar as a soil organic amendment has attracted increasing interest for its environmental and agronomic benefits, which is well-documented in several reviews and meta-analyses [1,2,3,4]. Biochar is a carbon-rich (60–80% C) material that can be obtained from different feedstocks (forest wood, manure, sewage sludge, agriculture residues) under a high-temperature treatment (from 300 to 1200 °C) and in an oxygen-limited environment (pyrolysis or gasification) [5,6,7]. With regard to environmental sustainability, biochar, due to its long-term stability in soil, has been reported to have great climate mitigation potential [8] for its ability to sequester atmospheric CO2 and reduce greenhouse gas (GHG) emission [9]. Concerning its agronomic implications, biochar has been shown to increase physical, chemical and microbiological soil fertility as well as reduce nutrient leaching and improve crop yield response [10,11,12,13]. In recent years, biochar application has also been suggested for controlling soil and water pollution by holding contaminants that are both inorganic (nitrogen, phosphorous, heavy metals) and organic (pesticides, herbicides, antibiotics) [14]. These biochar properties are mainly due to its highly porous structure, large inner surface area, greater negative surface charge and charge density [15,16].
Experimental evidence suggest that the physical properties of biochar (porosity and surface area) could improve not only soil physical properties (e.g., porosity, bulk density and aggregate stability) but also hydrological properties, including water retention and hydraulic conductivity [2,17]. Increases in soil water retention following biochar addition have been extensively reported in the literature, under both controlled (laboratory, lysimeters and pots) and open-field conditions. Several authors observed the rise in water retention on coarse- and medium-coarse textured soils amended with both lower [18,19] and higher biochar doses (>15 Mg ha1) [20,21]. Conversely, other authors reported no changes in water retention in fine-textured soils following biochar addition [22,23]. This clearly suggests that the biochar effect is largely dependent on soil type (texture and structure) and application rate. Moreover, it was reported that the effect of biochar on water retention is also linked to original feedstock and applied thermochemical conditions [23]. In this regard, an improvement of water retention was observed on coarse-textured soils after the application of two types of gasification biochar [19], while no effects were detected after applying biochar obtained from straw and pine wood at different pyrolysis temperatures [24]. A number of studies also report contrasting effects of biochar on saturated and unsaturated soil hydraulic conductivity (Ksat and Kunsat, respectively). Biochar sometimes increased, sometimes decreased Ksat, while limited results about the effects on Kunsat were reported. Similar to water retention, the biochar effect on hydraulic conductivity is dependent on soil texture [1]. Several studies [25,26] reported a Ksat decrease in sandy soil after biochar addition. Similar results were found by other authors [27] in coarse sand and fine sand soils, although they found a Ksat increase in clay loam soil. Likewise, Ksat increase was reported in poorly drained silty loam, silty clay and clay soil by several authors [16,28,29]. All these studies were conducted under controlled conditions (column, laboratory and greenhouse), but Ksat increase was also observed in field experiments on clay soils [30,31]. Soil hydraulic conductivity is also influenced by biochar type and application rate. Indeed, Kunsat increase was found after wood biochar addition to soil [25]; in laboratory experiments on clay soil, a significant Ksat increase was observed at higher biochar rates with no difference for Kunsat [32]. It was also reported that biochar hydrophobicity could induce water repellency and reduce hydraulic conductivity [33]. The hydrophobic nature of biochar decreases with increasing pyrolysis temperature [34], and biochar obtained at lower temperatures could be very water repellent due to alkyl groups on its surface [35].
If applying biochar to cultivated soil increases water retention, it could, in turn, increase plant-available water. This effect can be attributed to water retained in the pores within the biochar structure (intrapores) and those formed at the interface between biochar and soil particles or within small aggregates (interpores) [36]. Micropores can hold water strongly through capillary and adhesive forces, and mesopores contain water readily available to plants, while macropores are easily drained [1]. The increased micro- and meso-porosity due to biochar addition could explain why coarse-textured soils are more benefitted than fine-textured ones. Apart from soil porosity, also soil structure may play a relevant role [37]. In sandy loam, biochar promoted soil aggregation and aggregate stability due to the addition of organic carbon. Conversely, in clay soil, larger biochar doses strengthened repulsive forces between particles and monovalent cations, producing colloid dispersion and aggregate breakdown [37,38].
Moreover, biochar decreases hydraulic conductivity in sandy soils, thereby reducing the risk of nutrient leaching, whereas it increases hydraulic conductivity in fine-textured soils, thus enhancing water infiltration and limiting water runoff [1]. The mechanisms generally reported to explain these effects on conductivity in sandy soil are the filling or clogging of macropores by fine biochar particles and the increase in the number of micropores that bind water by strong capillary and adsorptive forces [27]; in fine-texture soil, the mechanisms are the improvement of aggregation and aggregate stability by biochar [39]. However, some exceptions to this general scheme were also detected [40,41,42,43].
In a previous paper from our research team [13], we reported the results of biochar application to a loam-texture soil in a 2-year field experiment made of four consecutive crop-growing cycles. The applied biochar played an important role in reducing water drainage loss, thereby limiting nutrient leaching. However, a “threshold effect” was observed: when a lower biochar amount was applied, the increase in water retention prevailed over the increase in water conductivity, causing diminished drainage; conversely, when a higher biochar amount was added, the reverse effect occurred, and drainage was greater. This must be considered our starting study hypothesis that must be verified. We did not perform any direct measurements of hydraulic soil properties in that previous field experiment, and our supposed explanation cannot be demonstrated until a new field experiment is performed. Nevertheless, prior to performing this further field trial, we considered working out a pot experiment in order to shed some light on this specific issue. Therefore, strictly controlled lab conditions were applied to a fine-textured soil similar to the one previously used. We are aware that, compared to in situ soils, the use of sieved and repacked soil in pots can be of concern when studying soil physical properties. Indeed, soil structure, pore architecture and pore size distribution are deeply altered when sieved and repacked soil is used, and the consequent soil hydraulic properties can change dramatically. However, under these strictly controlled experimental conditions, more extensive and accurate measurements can be performed and the effects of different biochar addition rates can be better compared and discriminated.
The aim of the present study was to evaluate whether biochar application could affect soil hydraulic properties in a clay loam soil in pots under lab conditions. The effect of different rates of biochar addition on bulk density, water retention curve, together with a set of water-capacitive indicators and water infiltration rates, was analyzed. No attempts were made to extrapolate the obtained results to “in-field” soil conditions, but deeper insights into how biochar operates in a soil mixture at different addition rates can surely promote more focused research activities in the open field.

2. Materials and Methods

Figure 1 shows the sequential phases applied when performing our experiment: from soil and biomass sampling, through biochar production and preparing soil–biochar mixtures, to measurements and data interpretation.

2.1. Soil and Biochar Preparation and Characterization

Soil samples were collected from the 0–0.20 m upper soil layer in Foggia, South Italy (latitude 41°27′ N, longitude 15°32′ E). The soil was oven dried at 105 °C until reaching a constant weight, after which it was crushed and passed through a 2 mm sieve. The particle size distribution was analyzed by the pipette-gravimetric method and the soil type was determined according to the USDA (United States Department of Agriculture) soil classification [44]. The pH was measured on a 1:2.5 (w/v) soil/water suspension with deionized water by a digital pH meter (GLP 22 pH-meter, Crison Instruments, Barcelona, Spain) and the electrical conductivity (EC) on saturated soil paste extract by a digital conductivity meter (GLP 31 EC-meter, Crison Instruments, Barcelona, Spain). Total nitrogen (total N) was determined by the Kjeldahl method [45] and the available phosphorus (P2O5) by extraction with sodium bicarbonate [46]. Organic carbon was determined by the dichromate oxidation method [47] and converted to organic matter (OM) by the conventional factor 1.724. The soil was of a clay loam texture, poor in nitrogen but rich in available phosphorus, with low organic-matter content, an alkaline pH and low electrical conductivity (Table 1).
Biochar applied to soil was produced at the STAR*Facility Centre, Foggia University, from residual vine biomasses (Vitis vinifera L.) collected from a local vineyard. The pruning residues (15% humidity) were chipped into particles of approx. 50 mm, mixed and then pyrolized at a temperature of 750 °C for 8 h, in a pilot scale with a fixed-bed tubular reactor (30 L capacity). The heating rate was 10 °C min−1. Once cooled, it was ground and passed through a 2 mm sieve.
The pH and electrical conductivity (EC) were determined on a 1:20 (w/v) soil/biochar suspension with deionized water, after shaking and waiting an equilibrium time of 90 min before measurement by a pH meter (GLP 22+ Crison Instruments, Barcelona) and an EC-meter (GLP 31+ Crison Instruments, Barcelona), respectively. The proximate analysis—the fixed carbon, volatile solids, ash and water content––was carried out by a TGA Analyzer (LECO-TGA701). The ultimate analysis of the carbon (C), hydrogen (H) and nitrogen (N) content was performed by dry combustion using a CHN Elemental Analyzer (CHN LECO628). In the case of Corg, dry combustion was preceded by inorganic C removal with HCl. Similarly, sulfur (S) was determined using a S module (S LECO628) combined with the CHN Elemental Analyzer. Oxygen (O) was calculated by difference. The carbon stability of biochar was evaluated indirectly by the molar ratios of hydrogen to organic carbon (H/Corg) and oxygen to organic carbon (O/Corg). The micro- and macro-elements were determined by digesting the dried biochar sample in concentrated HNO3, in a closed vessel microwave digester (CEM-Mars6). The metals in the solution were then analyzed by inductively coupling plasma optical-emission spectroscopy (ICP-OES Agilent 720). The physicochemical characteristics of the obtained biochar are listed in Table 2. It was characterized by the typical alkaline pH value [48] and resulted in a higher EC than the receiving soil. The prepared biochar fully complied with the standards established by both EBC (European Biochar Certificate) [49] and IBI (International Biochar Initiative) [50] because the C content was higher than 50%, the Corg higher than 60%, the H/Corg ratio lower than 0.7 and the O/Corg ratio lower than 0.4. Taking these properties into account, the available biochar was assigned to the Class 1 grade, while it also offered long-term carbon stability and persistence.

2.2. Experimental Treatments

Soil and biochar, both dried and sieved, were used to prepare soil–biochar mixtures. Biochar was added to the soil at 0, 2, 4, 6, 8 and 10% rates. The applied formula to set the correct amount of biochar in the mixture was the following:
B i = S   X i 100 X i
where B is the biochar dry weight, S the soil dry weight, X the percentage of biochar to be added, while subscript i refers to the soil–biochar specific treatment.
Six experimental treatments (soil–biochar mixtures) were arranged as follows: B0 (soil + 0% biochar or control), B2 (soil + 2% biochar), B4 (soil + 4% biochar), B6 (soil + 6% biochar), B8 (soil + 8% biochar) and B10 (soil + 10% biochar).
After soil–biochar mixture preparation, 5 L cylindrical plastic pots (internal Ø = 20 cm; h = 16 cm) were filled with the same mixtures. An equal volume of each soil–biochar mixture was placed in every pot without intentional compaction and each experimental treatment was done three times. Therefore, 18 pots made up the experimental setting to address the bulk-density measurement (Section 2.3) and the infiltration trial (Section 2.7). The other 18 pots (prepared according to the same procedure) were used to extract smaller soil samples to address the pressure plate apparatus (Section 2.4). Moreover, three identical pots were filled with an equal volume of biochar only. All these pots were sealed with a plastic sheet and kept at room temperature (22 ± 0.5 °C) for three months until the experimental tests were performed.

2.3. Bulk Density of Soil–Biochar Mixtures

The bulk density of both the soil–biochar mixtures and biochar alone was determined immediately before the start of the hydrological tests. After three months of storage, the soil–biochar mixtures in the pots did not show any kind of shrinkage. By gently trimming at the top it was possible to bring all the pots to full and equal volume without causing any disturbance to the underlying mixture. Pots were weighed and the bulk density (g cm−3) of each soil–biochar mixture (ρm) was determined as the ratio of the net soil–biochar mixture weight to the pot volume. Considering that treatment B0 represented the non-amended soil, its bulk density corresponded to the bulk density of the soil alone (ρs), that is, without biochar addition. The same procedure was also applied to the three pots filled with biochar only. In this way, it was possible to estimate the bulk density of biochar alone (ρb).

2.4. Water Content of Soil–Biochar Mixtures at Different Water Potentials

A pressure-plate apparatus equipped with a porous ceramic membrane (Soil Moisture Equipment Corp., Santa Barbara, CA, USA) was used to measure water content (ω, g g−1) of soil–biochar mixtures, according to the Richards’ method [51]. Samples of an equal volume of dry soil–biochar mixtures were excerpted from the pots using PVC sampling rings (Ø = 52 mm; h = 1 cm; volume = 21 cm3). On theses samples, the bulk density of the soil–biochar mixtures was further calculated in parallel with the procedure already reported in Section 2.3. Then, each sample was saturated with distilled water through partial soaking and by imbibition from below. As soon as saturation was reached, every sample was weighed and then placed on a porous ceramic plate before being placed into the pressure chamber. By using an air compressor, a constant pressure was established inside the chamber. The selected pressure corresponded to a specific value of the assigned matrix potential to be tested. Seven subsequent matric potentials (ψ, kPa), namely, saturation (ψ = 0 kPa) and ψ = −10, −30, −50, −100, −500 and −1500 kPa were selected.
The mixture ring samples were finally weighted and then dried at 105 °C for 24 h.
At every pressure step, the mixture water content was determined gravimetrically with respect to the dried weight. Therefore, the water content at each matric potential was expressed as percentage by weight of the soil–biochar mixtures (ω, g g−1).
The previously determined bulk density of each soil–biochar mixture (ρm) allowed converting the water content from a gravimetric (ω, g g−1) to a volumetric basis (θ, cm3 cm−3), according to the following equation:
θ = ω   ρ m ρ w
where ρw is the water density currently assumed equal to 0.998 g cm−3 at 20 °C.

2.5. Water Capacitive Indicators of Soil–Biochar Mixtures

A set of capacity-based soil water indicators [52] was applied to each of the six soil–biochar mixtures to compare their hydraulic properties. These indicators derive directly from the volumetric water contents (θ, cm3 cm−3) previously measured at characteristic matrix potentials, specifically: saturation (θSATψ = 0 kPa), field capacity (θFC: → ψ = −10 kPa), matrix capacity (θMRX: → ψ = −100 kPa) and wilting point (θWP: → ψ = −1500 kPa). It is generally remarked that the threshold value ψ = −100 kPa approximately discriminates structural from matric soil pores [53].
Therefore, the following capacitive indicators (Figure 2) were applied:
Total available water: TAW (%) = (θSATθWP) * 100;
Drainage water: DRW (%) = (θSATθFC) * 100;
Plant-available water: PAW (%) = (θFCθWP) * 100;
Readily available water: RAW (%) = (θFCθMRX) * 100;
Limited available water: LAW (%) = (θMRXθWP) * 100.
It can be easily observed that:
TAW = DRW + PAW;
PAW = RAW + LAW;
Therefore, the following further relative indicators were applied:
Relative available water (%): (PAW/TAW) * 100;
Relative readily available water (%): (RAW/PAW) * 100.

2.6. Soil Water Desorption Curve of Soil–Biochar Mixtures

Water content data on volumetric basis were used to describe the water desorption curve of the six soil–biochar mixtures. For this purpose, the van Genuchten model [54] was fitted to the experimental data by applying the following equation:
θ ( ψ )   =   θ r   +   θ s     θ r [ 1   + ( α   ψ ) n ] m  
where: θ (ψ) is the actual soil water content (cm3 cm−3) at the matric potential ψ (kPa); θr is the residual soil water content (cm3 cm−3) at limit values of matric potential (ψ → ∞); θs is the saturated water content (cm3 cm−3), assuming equivalence with total porosity; α is a dimensionless parameter related to the inverse value of the air entry potential (kPa−1); n and m are dimensionless curve shape parameters. In particular, m characterizes the asymmetry of the curve and the Mualem constraint m = 1 − 1/n was applied [55]. Therefore, Equation (10) only contains the following four unknown parameters: θs, θr, α, and n to be determined by a statistical fitting procedure (Section 2.5). Once specific estimates are assigned to the model parameters, further indicators useful in the soil–biochar mixture comparisons are the following:
  • matric potential of the inflection point of the water desorption curve (ψi, kPa):
ψ i =   1 a   ( n 1 n ) 1 n
  • water content at the inflection point, (θi, cm3 cm−3):
θ i = θ r + ( θ s θ r ) ( 1 + 1 m ) m
  • slope of the water retention curve at the inflection point (Si):
S i = n   ( θ s θ r ) ( 2 n 1 n 1 ) ( 1 n 2 )
It should be specified that the soil water release or desorption curve is expressed as volumetric water content (θ, cm3 cm−3) versus the absolute value of pore water tension (ψ, kPa) is expressed as natural logarithm: ln|ψ|.
The pore volume distribution function, P(ψ), may be defined as the slope (first derivative dθ/dψ) of the water release curve expressed as volumetric water content θ (cm3 cm−3), versus ln(ψ), and plotted against the equivalent pore diameter, δe (μm), on a logarithmic scale [52]. The equivalent pore diameter (δe) may be determined according to the Young–Laplace equation, and assuming that the pores are perfectly cylindrical, uniform, and equally drained, according to the following formula [41,52]:
δ e =   292 ψ
with δe (μm) and ψ (kPa).
Therefore, the van Genuchten curves were converted into pore volume distribution functions, and this can be considered the ultimate way to compare the six soil–biochar mixtures with respect to their hydraulic properties.

2.7. Water Infiltration of Soil–Biochar Mixtures

Water infiltration in the pots filled with the six soil–biochar mixtures was measured using a mini-disk infiltrometer (Meter Group Inc., Pullman, WA, USA). A constant pressure head of −2 cm was applied [56]. The amount of water (cm) that infiltrated into the soil–biochar mixtures was recorded manually at intervals of 10 s until all water from the reservoir chamber was infiltrated. The infiltration test was applied to every available pot and; therefore, three infiltration measurements were performed for each soil–biochar mixture.
Cumulative infiltration (IC, cm) was measured over time and the infiltration rate (IR, cm s−1) was determined as the difference between two cumulative infiltration (IC) readings over the corresponding time interval (t):
I R = ( I C ) 2 ( I C ) 1 t 2 t 1
In particular, three values of average infiltration rate were calculated. The initial infiltration rate (IRstart) was determined over the first 60 s of the infiltration run, the final infiltration rate (IRend) determined over the last 60 s of the infiltration run and, finally, the total infiltration rate (IRtotal), determined as the ratio of the total infiltrated water volume over the total time of the experimental infiltration process.
Cumulative infiltration (IC) data were analyzed by fitting the Philip’s infiltration model [56,57,58], as follows:
I C = C 1 t   +   C 2 t
where: IC (cm) is the cumulative infiltration; C1 (cm/s1/2) is the soil sorptivity; C2 (cm/s) is the parameter related to soil hydraulic conductivity; t is the time (s). A mathematical fitting procedure was applied to estimate the two coefficients, C1 and C2, and their assigned values were used as further indicators to compare the hydraulic properties of the soil–biochar mixtures. The instantaneous infiltration rate was determined by differentiating, with respect to time, the cumulative infiltration as reported in the following equation:
I R = 0.5   C 1 / t   +   C 2

2.8. Statistical Analysis

The soil–biochar bulk-density data were processed according to an ANCOVA model considering the biochar addition rate (from zero to 10%) as the regressing variable. Both “expected” and “observed” data were the response variables. The former being the data theoretically compatible with a simple “dilution” effect of biochar in the mixtures, while the latter the actually measured data. In other words, a “lighter” component (the lower density biochar) when added to a “heavier” medium (the higher density soil) produced a decrease in the mixture bulk density proportional to the amount (weight) of the added component, according to the following equation:
ρm(i) = ρs (100 − Xi)/100 + ρb X/100
where, as already stated, X is the percentage of the added biochar, while the subscript i refers to the soil–biochar specific treatment. According to Equation (18), the expected values of bulk density can be calculated. Therefore, the ANCOVA model allowed testing if the decreasing effect of biochar on bulk density (i.e., the slope of the regression line) was statistically different in the two compared data sets (expected vs. observed).
A one-way ANOVA test was applied to statistical process all the water-capacitive indicators as well as the average infiltration rates with respect to the experimental treatments (the six soil–biochar mixtures). A mean comparison was performed via a Tukey honestly significant difference (HSD) test, and a 0.05 probability level was applied. A multivariate analysis was also performed on the water-capacitive indicators, which were processed as a whole through a cluster analysis. The “average linkage method” was applied, and the resulting dendrogram was represented considering the detected clusters.
The parameter estimates of the two models, van Genuchten and Philip’s, were obtained through an iterative, nonlinear, least-square fitting procedure. In this respect, a “stepwise” approach was applied according to the following procedure:
  • The model was fitted independently to the data from every experimental Bi treatment (from B0 to B10) to obtain specific parameter estimates. The overall experimental error (SSe—Sum of Square error) was obtained as the sum of SSe related to each treatment curve (ΣSSi).
  • The model was then fitted to the complete dataset, regardless of experimental treatments, to obtain a single, general estimate of each model parameter.
  • The general model (point 2) was gradually relaxed, one parameter at a time, keeping the others fixed. The statistical significance of an independent estimates of each Bi treatment i was evaluated, one parameter at a time, by means of an F Fisher test against the error term (SSe).
  • We proceeded with point 3 in order of relevance, from the most to the least influential parameter. The additional SS, accounted for by the model when fitting Bi-independent estimates of each parameter, represented its relevance.
  • A significant F-test allowed the conversion of the model from one single estimate of the considered parameter to different estimates, according to a multiple range comparison among the Bi treatments. A non-significant F-test allowed the preservation of a unique estimate of the considered parameter.
  • The same procedure was repeated each time considering the remaining model parameters; every time, a new fitting procedure assigned new estimates to the remaining parameters.
All data processing and statistical tests were performed using the JMP statistical software package, version 11 (SAS Institute Inc, Cary, NC, USA).

3. Results and Discussion

3.1. Bulk Density of Soil–Biochar Mixtures

Two independent and parallel bulk-density estimates of the six soil–biochar mixtures were obtained. The first estimates derived from weighing the entire pots (Section 2.3), while the second from weighing the smaller ring samples placed in the pressure-plate apparatus (Section 2.4). No statistical differences were detected between the two datasets (data not showed). Since the first estimates came from a greater volume of mixtures (the pots) than the second ones (ring samples), the former were considered more reliable and taken as a reference.
The bulk density of the soil–biochar mixtures decreased linearly with the biochar fraction added to the soil (Figure 3).
Our results were consistent with other experimental studies [1] as well as meta-analysis [59]. Most studies showed a reduction in bulk density as a consequence of biochar addition although few cases also reported increase in bulk density [17]. A biochar addition roughly higher than 2% by weight should be enough to create a significant decrease in bulk density [60]. The rate of biochar application as well as the density and porosity of the original soil are critical for predicting the effects of a biochar addition to soil. In our experiment, the biochar bulk density was 0.321 ± 0.074 g cm−3, while the soil bulk density was 1.191 ± 0.079. Therefore, considering an addition of biochar (10% by weight), a bulk density decrease of 12.6% was actually detected compared to an expected decrease of 7.3% if only a “dilution” effect were supposed. This decrease, therefore, does not follow the proportions resulting from a simple “dilution” effect (red line in Figure 3) since the reduction observed (blue line in Figure 3) was significantly stronger than expected. This was proved by the result of the ANCOVA test reported in Table 3.
As far as we know, only one other paper [53] compared the expected “dilution” data with the actual observed data of bulk density, thereby showing that the effect of biochar in reducing the soil bulk density was more than proportional to the simple contribution of its weight. Consequently, it is appropriate to call into question the effects of the spatial arrangement of biochar particles in relation to soil particles or small aggregates. This implies an increase in total porosity not only due to the intrinsic biochar microporosity, but also due to a soil–biochar interporosity. Apart from reasons of simple particle spatial arrangement, the effect of soil–biochar interaction may also lead to the development of structural aggregates. Considering that biochar increases the mixture cation exchange capacity (CEC), the bridging of clays and organic matter may be promoted; on the other hand, the increase in soil reaction (pH) may also increase surface nonacidic cations bridging soil colloids [53]. A better particle aggregation may promote soil structuring and a decrease in bulk density. Considering the limited time available for soil structuring in this experiment, we believe the former factor should be more effective than the latter.

3.2. Water Capacitive Indicators of Soil–Biochar Mixtures

The one way ANOVA test performed on each capacitive indicator (Table 4) gave quite clear results concerning the different hydraulic properties of the six soil–biochar mixtures. TAW, PAW, RAW, and LAW increased progressively with increasing biochar application fractions; conversely, DRW decreased. In general, biochar effectively contributed to improved soil quality with respect to water plant–soil relationships. When comparing the maximum range between treatment pairs, water holding capacity (TAW) and water availability to plants (PAW, RAW and LAW) increased significantly. Mostly importantly, the highest increase was observed in RAW (+73%), followed by PAW (+61%) and LAW (+54%), while DRW was reduced by approximately 46%. Not only the amount of available water to plants was significantly expanded, but also the amount of drainage water (lost water or only temporary available to plants) was significantly reduced. Concerning the relative indicators, PAW/TAW grew significantly larger with the biochar addition (+32%, the maximum range between treatment pairs), while RAW/PAW did not significantly change.
Several experimental evaluations [17], both in the field and in the lab, demonstrated that soil amendment with biochar is associated with a significant increase in PAW. The magnitude of this increase is related to the soil’s being higher in coarse-textured and lower in fine-textured soils, as well as roughly proportional to the biochar application rate. Similar to what was previously observed with bulk density, this effect is mostly interpreted in reference to a remarkable increase in the total soil porosity and specific surface area [59]. According to the same meta-analysis [59] the effect of biochar on PAW increases as its porosity increases. Biochar has been seen to increase not only total soil porosity but also pore connectivity and the number of pores [61,62], thus increasing soil water retention. As remarked already, biochar intrapores (pores inside the biochar particles) together with soil–biochar interpores jointly contribute to increasing water retention [63]. According to our results, the increase in PAW due to the biochar addition was mostly determined by a significant increase in water content at FC (+27.9%, the maximum range between treatment pairs).
The treatment with the highest biochar fraction did not always show the most noticeable results; indeed, with respect to PAW and RAW, it was B8, not B10, that offered the best soil conditions although it was not statistically different. A sort of “threshold effect” in biochar application was detected in this lab experiment. Interestingly, the same effect was also observed in our previous field-conducted experiment [13]. Indeed, B8 and B10 showed a general trend (not statistically confirmed) to swap places. This was particularly clear when the result of the cluster analysis jointly applied to all the data reported in Table 4 are shown (Figure 4). A clear split into two clusters (B0–B2–B4 vs. B6–B8–B10) was shown. Moreover, a closer “friendship” between B6 and B10 rather than B8 and B10 was observed. It is difficult to provide an explanation of the observed swapping between B10 and B8, but it seemed that no further porosity was added to the mixture; that is to say, the “useful” porosity to the plant water supply was partially lost in B10 with respect to B8.

3.3. Water Retention Curves of Soil–Biochar Mixtures

The fitting procedure of the van Genuchten equation, performed according to a “stepwise” approach, provided a model consisting of only two curves significantly different from each other. The first curve included the following treatments: B0, B2 and B4, while the second one the other three treatments: B6, B8 and B10. The pair of curves differed significantly in the value of θs and n (Table 5), while neither α nor θr had statistically different results among the Bi treatments (Table 5). In particular, θr never differed significantly from zero in any Bi treatment, and for this reason its value was zeroed (Table 5).
The most influential effects of biochar application was a significant increase in the total volume of water at saturation (θs). As a consequence, the soil water content also increased at lower water potentials, thus improving plant water availability. The increase in the value of the n parameter further operated in this direction. Indeed, the increase in the n “shape” parameter accentuated the steepness of the curve, simultaneously increasing the water volume at FC and, by contrast, decreasing it at WP. Consequently, a remarkable increase in PAW was obtained and this effect represented another relevant result of biochar application: a significant enhanced capacity of the amended soil to hold water not only at soil saturation but also at unsaturated soil conditions.
Beyond the rigorous statistical separation performed with the “stepwise” approach, both θs and n showed a gradual, progressive increase in their values from B0 to B10 (Table 6). In contrast, no kind of trend was detected with respect to the α parameter (Table 6).
The effect of the biochar application on both θs and n allowed us to deduce that its main consequence is a substantial, additional contribution to soil porosity that can increase water sorption capacity. While the increasing trend on θs seemed quite regular, the values related to the n parameter reached a maximum at B8 and did not increase at B10, thus limiting any further increase in soil water-holding capacity (Table 6). Again, as previously noticed, a peculiar effect of biochar at a higher application rate (although not statistically confirmed) was observed.
According to the statistical results reported in Table 5, Figure 5A shows the two desorption curves with respect to the first (B0, B2 and B4) and the second (B6, B8, B10) group of treatments. The two curves differed significantly, visually confirming what was previously observed with regard to the capacitive indicators and the van Genuchten parameters. They showed similar values around and beyond water potential (WP) (N.B.: ψ in absolute terms); conversely, they showed a divergent trend at increasingly lower water potentials (N.B.: ψ in absolute terms), in the range of PAW as well of DRW.
Figure 5B represents, even more clearly, the conditions occurring in the soil as a consequence of biochar addition from lower (B0, B2 and B4) to higher quantities (B6, B8, B10). What could be observed in the pore-size distribution of the soil–biochar mixtures was a significant increase in that fraction of porosity straddling the size of the pores that marks the transition from “matrix” porosity (intra-aggregate voids) to “structural” porosity (interaggregate voids). Indeed, among aggregates, pores are formed that hold water primarily by capillary forces associated with the clumping of aggregates; these can be termed “structural” pores. In contrast, “matrix” pores are intra-aggregate, formed among elementary textural particles where interparticle adhesive forces principally hold water. In this respect, it is usually assumed that the former fraction of soil water is more easily available for the plant than the latter.
The inflection points of the two curves in Figure 5A corresponded to the maximum value of the two curves in Figure 5B. They represented the slope of the water-retention curves at their respective inflection points (Si). The Si value was significantly higher in the blue line (B6, B8, B10) than in the red one (B0, B2 and B4), although they were both achieved at exactly the same water potential (170 kPa), corresponding to an equivalent pore size equal to 1.72 μm. As expected, the soil water content at the inflection point (θi) was higher in the blue line (B6, B8, B10) than in the red one (B0, B2 and B4), further stressing the improved soil quality offered by biochar application with respect to soil-water relationships.
The symmetry observed in the two pore size-distribution curves suggests that biochar addition at higher rates did not produce a curve shifting and significantly increased the amount of the existing pore classes, both intra- and inter-aggregates (i.e., at both sides of MRX, see Figure 5B). This general increase in porosity can explain both the reduction in bulk density of the mixtures, which was more than proportional to the bulk density of the added biochar, and the increase in the mixture water-retention capacity. Indeed, the first effect can be predominantly accounted for by the increase in the structural (interaggregate) porosity, while the second effect can be predominantly associated with interporosity with respect to RAW and intraporosity with respect to LAW. Adding RAW and LAW gives PAW.
If we consider that the peak of the two curves in Figure 5B roughly corresponds to the edge that separates the micro- from the meso-porosity, it is possible to understand the beneficial effect exerted by the addition of biochar in favouring water retention, precisely at the most convenient potentials for plant-water uptake (i.e., those in the range WP–FC).

3.4. Water Infiltration of Soil–Biochar Mixtures

The average infiltration rate of water into the six different soil–biochar mixtures showed decreasing values as the biochar application rates from B0 to B10 (Table 7) increased. This trend was observed with respect to the initial (first 60 s of the process), the final (last 60 s) and the total average value of infiltration rates (the entire process duration). Considering the B8 treatment, a lower value for the final and total rate was detected compared to that of B10. From a statistical perspective, the ANOVA results determined only two distinct treatment groups, namely, the first, represented by only the B0 treatment, and the second, which comprises all other treatments (Table 7).
The ANCOVA model, applied to the infiltration data in order to estimate the coefficients of the Philip’s model, confirmed the same statistical results shown in Table 7; the two model coefficients, C1 and C2, were significantly affected by the biochar presence or absence but not significantly by the biochar amount (Table 8).
Both the “sorptivity” coefficient (C1) and the “conductivity” coefficient (C2), showed higher values in soils without biochar application and lower values when biochar was added, regardless of the biochar amount applied. These values fully confirmed the results of the ANOVA test performed on the initial, final and total average infiltration rates, as reported in Table 6.
The application of the Philip’s model considering both the cumulative infiltration and the instantaneous infiltration rate is represented in Figure 6.
Considering Equation (16), the first-term of the model (controlled by the “sorptivity” coefficient C1) strongly dominated the process at the initial stage of infiltration. At this stage the soils were almost dried and the capillary and matrix forces at the soil pore surfaces largely prevailed over the gravitational forces because the larger pores were still empty. Therefore, vertical water movement proceeded at almost the same rate as horizontal water movement. The “gravity” component, controlled by the second term of the model (the “conductivity” coefficient C2), exerted a negligible effect at the very beginning of the infiltration process, but became progressively more important as infiltration proceeded and the soils became wetter. Finally, the second term of the Philip’s equation strongly dominated the process when the soils were almost under saturated condition and the gravitational forces exerted a prevailing effect, while capillarity and matrix forces were, this time, negligible.
Considering the overall infiltration results (Table 7 and Table 8 and Figure 6), it appeared that the predominant effect was the presence or absence of biochar, rather than the quantity added to soil. Biochar slowed infiltration into the soil compared to soil without biochar applied. This effect was significantly greater than the effect exerted by the different biochar concentrations in the soil. These results are not in agreement with those generally observed in fine- and fine-loam soils. According to a review paper [1], biochar application to a fine-textured soils at 5 and 10% by weight increased unsaturated hydraulic conductivity at all matric potentials and generally increased water infiltration. According to the same review [1], the application of biochar to fine-textured soils can improve water flow (both infiltration and hydraulic conductivity) because of an increase in pore volume and an improved soil aggregation. It is generally recognized that the addition of biochar significantly enhanced the formation of macro-aggregates that slightly increased the saturated hydraulic conductivity of the amended soils [64]. However, it should be remembered that in our experiment a sieved and repacked soil was used, and these particular conditions may have profoundly affected the results, which would be very different if the soil had been in situ. Indeed, it was also observed [59] that, as a general trend, biochar addition reduced soil-saturated hydraulic conductivity in laboratory experiments compared to those in field studies but without considering specific soil textures. It should be expected that under our lab experiment there was neither the time nor the proper conditions to promote the development of a good soil structure with macro-aggregates because of the contribution of biochar as a soil amendment.
Concerning the effects displayed by the biochar application to soil, on the one hand, the augmented intraporosity conferred by biochar (with small pore size) significantly reduced soil water conductivity when soil was wet and almost saturated [27]. On the other hand, when the soil was still almost dry, we can suppose that biochar showed its hydrophobic property, thus exerting a water-repellent action that significantly reduced the matrix and capillary effect of water attraction by pores and surfaces. This could be the reason that the “sorptivity” coefficient was lower in biochar rather than in the non-amended soil. Biochar hydrophobicity has been largely discussed in the literature [1,33,35]. It was observed, indeed, that biochar’s microstructure is often highly hydrophobic, suggesting that when added to soil it can induce water repellency. This higher water repellency, in turn, can be the possible explanation for reduced infiltration and hydraulic conductivity. It was also observed that, over time, the biochar in the soil in situ undergoes oxidation involving organic compounds on the surface of biochar particles. These oxidations can progressively mitigate biochar’s hydrophobic property [1], thus favoring the restoration of higher water conductivity. Results from the literature indicate that “aged” biochar is more hydrophilic then “fresh” biochar and that after additional wetting/drying treatments biochar become less hydrophobic [65,66,67].

4. Conclusions

The results obtained from a lab experiment with soil in pots confirmed that biochar has the potential to improve the hydrological properties of clay loam soil. We observed, on the one hand, a substantial increase in water retention capacity; on the other, a reduced water infiltration rate. The first effect agrees with the majority of experimental evidence retrieved from the literature, while the second one is somehow not aligned with the prevailing literature. Our interpretation is based on the relevant increase in porosity (intrapore and interpore) that biochar addition brings to soil, permitting higher water adsorption. The altered soil structural condition, the limited time elapsed between biochar application and the experimental measurements, together with the storage of the soil in pots under a sealed covering, prevented the development of macrostructural soil aggregates able to generate larger pores and facilitate internal water fluxes. Furthermore, water infiltration and its conductivity were hampered by the hydrophobic action exerted by “fresh” biochar (i.e., obtained only recently by pyrolysis and applied to the soil).
The positive effect on water-holding capacity showed that biochar has the potential to effectively improve the water available for plants and contribute to reducing the frequency and volume of irrigation in fine-textured soils. These aspects are of great agronomic relevance and can be confirmed in “open field” conditions, too.
The observed effects of biochar with respect to slower water infiltration and limited water conductivity, if confirmed under “open filed” conditions, may have a dual consequence. When water infiltration is hindered, the risk of water runoff and soil erosion strongly increases, as does possible waterlogging; this is a negative factor. Conversely, a reduced water movement due to limited hydraulic conductivity significantly slows deep water drainage and consequent nutrient loss through leaching, the latter being a positive factor. Which of the two should be considered as the prevailing factor and which kind of tradeoff should be established largely depends on the environmental farming conditions and operations. Therefore, further experiments under field conditions are needed to clarify these aspects.
Particularly in agricultural regions characterized by climatic water shortage and frequently droughts, such as the Mediterranean, greater water availability and reduced water loss through drainage or runoff due to biochar application could have an agronomic benefit that is perfectly in line with the sustainable management of water resources for irrigation.

Author Contributions

Conceptualization, A.L. and M.M.; methodology, A.L. and M.M.; formal analysis, A.L., M.F. and M.M.; investigation, A.L. and M.M.; data curation, A.L., M.F. and M.M.; writing—original draft preparation, A.L. and M.M.; writing—review and editing, A.L. and M.M.; visualization, A.L. and M.M.; supervision, A.L. and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors are grateful to Teresa Gentile and Valentina Del Vento for the technical assistance in the laboratory experiment and analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Blanco-Canqui, H. Biochar and soil physical properties. Soil Sci. Soc. Am. J. 2017, 81, 687–711. [Google Scholar] [CrossRef] [Green Version]
  2. Ding, Y.; Liu, Y.; Liu, S.; Huang, X.; Li, Z.; Tan, X.; Zeng, G.; Zhou, L. Potential benefits of biochar in agricultural soils: A Review. Pedosphere 2017, 27, 645–661. [Google Scholar] [CrossRef]
  3. Kavitha, B.; Reddy, P.V.L.; Kim, B.; Lee, S.S.; Pandey, S.K.; Kim, K.H. Benefits and limitations of biochar amendment in agricultural soils: A review. J. Environ. Manag. 2018, 227, 146–154. [Google Scholar] [CrossRef]
  4. Semida, W.M.; Beheiry, H.R.; Setamou, M.; Simpson, C.R.; Abd El-Mageed, T.A.; Rady, M.M.; Nelson, S.D. Biochar implications for sustainable agriculture and environment: A review. S. Afr. J. Bot. 2019, 127, 333–347. [Google Scholar] [CrossRef]
  5. Zabaniotou, A.; Rovas, D.; Libutti, A.; Monteleone, M. Boosting circular economy and closing the loop in agriculture: Case study of a small-scale pyrolysis-biochar based system integrated in an olive farm in symbiosis with an olive mill. Environ. Dev. 2015, 14, 22–36. [Google Scholar] [CrossRef]
  6. Monlau, F.; Francavilla, M.; Sambusiti, C.; Antoniou, N.; Solhy, A.; Libutti, A.; Zabaniotou, A.; Barakat, A.; Monteleone, M. Toward a functional integration of anaerobic digestion and pyrolysis for a sustainable resource management. Comparison between solid-digestate and its derived pyrochar as soil amendment. Appl. Energy 2016, 169, 652–662. [Google Scholar] [CrossRef]
  7. Zabaniotou, A.; Rovas, D.; Delivand, M.K.; Francavilla, M.; Libutti, A.; Cammerino, A.R.B.; Monteleone, M. Conceptual vision of bioenergy sector development in Mediterranean regions based on decentralized thermochemical systems. Sustain. Energy Technol. Assess. 2017, 23, 33–47. [Google Scholar] [CrossRef]
  8. Woolf, D.; Amonette, J.E.; Street-Perrott, F.A.; Lehmann, J.; Joseph, S. Sustainable biochar to mitigate global climate change. Nat. Commun. 2010, 1, 56. [Google Scholar] [CrossRef] [Green Version]
  9. Wang, J.; Xiong, Z.; Kuzyakov, K. Biochar stability in soil: Meta-analysis of decomposition and priming effects. GCB Bioenergy 2015, 8, 512–523. [Google Scholar] [CrossRef] [Green Version]
  10. Agegnehu, G.; Nelson, N.; Bird, M.I. Crop yield, plant nutrient uptake and soil physicochemical properties under organic soil amendments and nitrogen fertilization on nitisols. Soil Tillage Res. 2016, 160, 1–13. [Google Scholar] [CrossRef]
  11. Demiraj, E.; Libutti, A.; Malltezi, J.; Rroço, E.; Brahushi, F.; Monteleone, M.; Sulçe, S. Effect of organic amendments on nitrate leaching mitigation in a sandy loam soil of Shkodra district, Albania. Ital. J. Agron. 2018, 13, 1136. [Google Scholar] [CrossRef] [Green Version]
  12. Libutti, A.; Mucci, M.; Francavilla, M.; Monteleone, M. Effect of biochar amendment on nitrate retention in a silty clay loam soil. Ital. J. Agron. 2016, 11, 273–276. [Google Scholar] [CrossRef] [Green Version]
  13. Libutti, A.; Cammerino, A.R.B.; Francavilla, M.; Monteleone, M. Soil amendment with biochar affects water drainage and nutrient losses by leaching: Experimental evidence under field-grown conditions. Agronomy 2019, 9, 758. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, C.; Zeng, G.; Huang, D.; Lai, C.; Chen, M.; Cheng, M.; Tang, W.; Tang, L.; Dong, H.; Huang, B.; et al. Biochar for environmental management: Mitigating greenhouse gas emissions, contaminant treatment, and potential negative impacts. Chem. Eng. J. 2019, 373, 902–922. [Google Scholar] [CrossRef]
  15. Wang, S.; Gao, B.; Zimmerman, A.R.; Li, Y.; Ma, L.; Harris, W.G.; Migliaccio, K.W. Physicochemical and sorptive properties of biochars derived from woody and herbaceous biomass. Chemosphere 2015, 134, 257–262. [Google Scholar] [CrossRef]
  16. Herath, H.M.S.K.; Camps-Arbestain, M.; Hedley, M. Effect of biochar on soil physical properties in two contrasting soils: An Alfisol and an Andisol. Geoderma 2013, 209–210, 188–197. [Google Scholar] [CrossRef]
  17. Omondi, M.O.; Xia, X.; Nahayo, A.; Liu, X.; Korai, P.K.; Pan, G. Quantification of biochar effects on soil hydrological properties using meta-analysis of literature data. Geoderma 2016, 274, 28–34. [Google Scholar] [CrossRef]
  18. Głąb, T.; Palmowska, J.; Zaleski, T.; Gondek, K. Effect of biochar application on soil hydrological properties and physical quality of sandy soil. Geoderma 2016, 281, 11–20. [Google Scholar] [CrossRef]
  19. Hansen, V.; Nielsen, H.H.; Petersen, C.T.; Mikkelsen, T.N.; Stöver, D.M. Effects of gasification biochar on plant-available water capacity and plant growth in two contrasting soil types. Soil Tillage Res. 2016, 161, 1–9. [Google Scholar] [CrossRef] [Green Version]
  20. Xiao, Q.; Zhu, L.; Zhang, H.; Li, X.; Shen, Y.; Li, S. Soil amendment with biochar increases maize yields in a semi-arid region by improving soil quality and root growth. Crop Pasture Sci. 2016, 67, 495–507. [Google Scholar] [CrossRef]
  21. Yang, C.D.; Lu, S.G. Effects of five different biochars on aggregation, water retention and mechanical properties of paddy soil: A field experiment of three-season crops. Soil Till. Res. 2021, 205, 104798. [Google Scholar] [CrossRef]
  22. Carvalho, M.; Madari, B.; Bastiaans, L.; Oort, P.V.; Leal, W.; Heinemann, A.; da Silva, M.A.S.; Maia, A.H.N.; Parson, D.; Meinke, H. Properties of a clay soil from 1.5 to 3.5 years after biochar application and the impact on rice yield. Geoderma 2016, 276, 7–18. [Google Scholar] [CrossRef]
  23. Wang, D.; Li, C.; Parikh, S.J.; Scow, K.M. Impact of biochar on water retention of two agricultural soils—A multi-scale analysis. Geoderma 2019, 340, 185–191. [Google Scholar] [CrossRef]
  24. Burrell, L.D.; Zehetner, F.; Rampazzo, N.; Wimmer, B.; Soja, G. Long-term effects of biochar on soil physical properties. Geoderma 2016, 282, 96–102. [Google Scholar] [CrossRef]
  25. Uzoma, K.C.; Inoue, M.; Andry, H.; Zahoor, A.; Nishihara, E. Influence of biochar application on sandy soil hydraulic properties and nutrient retention. J. Food Agric. Environ. 2011, 9, 1137–1143. [Google Scholar]
  26. Zhang, J.; Chen, Q.; You, C.F. Biochar effct on water evaporation and hydraulic conductivity in sandy soil. Pedosphere 2016, 26, 265–272. [Google Scholar] [CrossRef]
  27. Lim, T.J.; Spokas, K.A.; Feyereisen, G.; Novak, J.M. Predicting the impact of biochar additions on soil hydraulic properties. Chemosphere 2016, 142, 136–144. [Google Scholar] [CrossRef]
  28. Ajayi, A.; Horn, R. Modification of chemical and hydrophysical properties of two texturally diffrentiated soils due to varying magnitudes of added biochar. Soil Tillage Res. 2016, 164, 34–44. [Google Scholar] [CrossRef]
  29. Barnes, R.T.; Gallaghe, M.E.; Masiello, C.A.; Liu, Z.; Dugan, B. Biochar induced changes in soil hydraulic conductivity and dissolved nutrient flxes constrained by laboratory experiments. PLoS ONE 2014, 9, e108340. [Google Scholar] [CrossRef] [Green Version]
  30. Asai, H.; Samson, B.K.; Stephan, H.M.; Songyikhangsuthor, K.; Homma, K.; Kiyono, Y.; Inoue, Y.; Shiraiwa, T.; Horie, T. Biochar amendment techniques for upland rice production in Northern Laos. Field Crops Res. 2009, 111, 81–84. [Google Scholar] [CrossRef]
  31. Major, J.; Lehmann, J.; Rondon, M.; Goodale, C. Fate of soil-applied black carbon: Downward migration, leaching and soil respiration. Glob. Chang. Biol. 2010, 16, 1366–1379. [Google Scholar] [CrossRef]
  32. Kameyama, K.; Miyamoto, T.; Shiono, T.; Shinogi, Y. Influence of sugarcane bagasse-derived biochar application on nitrate leaching in calcaric dark red soil. J. Environ. Qual. 2012, 41, 1131–1137. [Google Scholar] [CrossRef] [PubMed]
  33. Githinji, L. Effect of biochar application rate on soil physical and hydraulic properties of a sandy loam. Arch. Agron. Soil Sci. 2014, 60, 457–470. [Google Scholar] [CrossRef]
  34. Gray, M.; Johnson, M.G.; Dragila, M.I. Water uptake in biochars: The roles of porosity and hydrophobicity. Biomass Bioenergy 2014, 61, 196–205. [Google Scholar] [CrossRef]
  35. Kinney, T.J.; Masiello, C.A.; Dugan, B.; Hockaday, W.C.; Dean, M.R.; Zygourakis, K.; Barnes, R.T. Hydrologic properties of biochars produced at different temperatures. Biomass Bioenergy 2012, 41, 34–43. [Google Scholar] [CrossRef]
  36. Liu, Z.; Dugan, B.; Masiello, C.A.; Gonnermann, H.M. Biochar particle size, shape, and porosity act together to influence soil water properties. PLoS ONE 2017, 12, e0179079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Pituello, C.; Dal Ferro, N.; Francioso, O.; Simonetti, G.; Berti, A.; Piccoli, I.; Pisi, A.; Morari, F. Effects of biochar on the dynamics of aggregate stability in clay and sandy loam soils. Eur. J. Soil Sci. 2018, 69, 827–842. [Google Scholar] [CrossRef]
  38. Kumari, K.G.I.D.; Moldrup, P.; Paradelo, M.; Elsgaard, L.; de Jonge, L.W. Effect of Biochar on Dispersibility of Colloids in Agricultural Soils. J. Environ. Qual. 2017, 46, 143–152. [Google Scholar] [CrossRef]
  39. Ghorbani, M.; Asadi, H.; Abrishamkesh, S. Effects of rice husk biochar on selected soil properties and nitrate leaching in loamy sand and clay soil. Int. Soil Water Conserv. Res. 2019, 7, 258–265. [Google Scholar] [CrossRef]
  40. Major, J.; Rondon, M.; Molina, D.; Riha, S.J.; Lehmann, J. Nutrient leaching in a Colombian savanna oxisol amended with biochar. J. Environ. Qual. 2012, 41, 1076–1086. [Google Scholar] [CrossRef]
  41. Hardie, M.; Clothier, B.; Bound, S.; Oliver, G.; Close, D. Does biochar influence soil physical properties and soil water availability? Plant Soil 2014, 376, 347–361. [Google Scholar] [CrossRef]
  42. Jeffery, S.; Meinders, M.B.J.; Stoof, C.R.; Bezemer, T.M.; van de Voorde, T.F.J.; Mommer, L.; van Groenigen, J.W. Biochar application does not improve the soil hydrological function of a sandy soil. Geoderma 2015, 251–252, 47–54. [Google Scholar] [CrossRef]
  43. Castellini, M.; Giglio, L.; Niedda, M.; Palumbo, A.D.; Ventrella, D. Impact of biochar addition on the physical and hydraulic properties of a clay soil. Soil Tillage Res. 2015, 154, 1–13. [Google Scholar] [CrossRef]
  44. Soil Survey Staff. Soil Survey Staff. Soil taxonomy: A basic system of soil classification for making and interpreting soil surveys. In Natural Resources Conservation Service. Handbook 436, 2nd ed.; United States Department of Agriculture: Washington, DC, USA, 1999. [Google Scholar]
  45. Bremner, J.M. Nitrogen-total. In Methods of Soil Analysis Part 3—Chemical Methods; SSSA Book Ser. No. 5; Sparks, D.L., Page, A.L., Johnston, C.T., Summ, M.E., Eds.; SSSA: Madison, WI, USA, 1996; pp. 1058–1121. [Google Scholar]
  46. Olsen, S.R.; Cole, C.V.; Watanabe, F.S.; Dean, L.A. Estimation of Available Phosphorus in Soil by Extraction with Sodium Bicarbonate; USDA Circular 939; USDA: Washington, DC, USA, 1954; pp. 1–19.
  47. Walkley, A.; Black, I.A. An examination of the Degtjareff method fordetermining soil organic matter and a proposed modification of the chromicacid titration method. Soil Sci. 1934, 37, 29–38. [Google Scholar] [CrossRef]
  48. Chan, K.Y.; Xu, K. Biochar: Nutrient properties and their enhancement. In Biochar for Environmental Management. Science and Technology; Lehmann, J., Joseph, S., Eds.; Earthscan: London, UK, 2009; pp. 67–84. [Google Scholar]
  49. EBC. European Biochar Certificate Guidelines for a Sustainable Production of Biochar; European Biochar Certificate (EBC): Arbaz, Switzerland, 2012. [Google Scholar]
  50. IBI. International Bichar Inititative, Standardized Product Definition and Product Testing Guidelines for Biochar That Is Used in Soil; IBI-STD-2.0; IBI: Toronto, ON, Canada, 2014. [Google Scholar]
  51. Klute, A.; Dirksen, C. Hydraulic conductivity and diffusivity: Laboratory methods. In Methods of Soil Analysis-Part 1. Physical and Mineralogical Methods; Klute, A., Ed.; Wiley: Hoboken, NJ, USA, 1986; pp. 687–734. [Google Scholar]
  52. Reynolds, W.D.; Drury, C.F.; Tan, C.S.; Fox, C.A.; Yang, X.M. Use of indicators and pore volume-function characteristics to quantify soil physical quality. Geoderma 2009, 152, 252–263. [Google Scholar] [CrossRef]
  53. Walters, R.D.; White, J.G. Biochar in Situ Decreased Bulk Density and Improved Soil-Water Relations and Indicators in Southeastern US Coastal Plain Ultisols. Soil Sci. 2018, 183, 99–111. [Google Scholar] [CrossRef]
  54. van Genuchten, M.T. A closed form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  55. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 1976, 12, 513–522. [Google Scholar] [CrossRef] [Green Version]
  56. METER Group. Mini Disk Infiltrometer-User’s Manual; Decagon Devices, Inc.: Pullman, WA, USA, 2014; pp. 10–14. [Google Scholar]
  57. Zhang, R. Infiltration models for the disk infiltrometer. Soil Sci. Soc. Am. J. 1997, 61, 1597–1603. [Google Scholar] [CrossRef]
  58. Philip, J.R. The theory of infiltration: 4. Sorptivity and algebraic infiltration equations. Soil Sci. 1957, 84, 257–264. [Google Scholar] [CrossRef]
  59. Edeh, I.G.; Mašek, O.; Buss, W. A meta-analysis on biochar’s effects on soil water properties–New insights and future research challenges. Sci. Total Environ. 2020, 714, 136857. [Google Scholar] [CrossRef]
  60. Kayode, S.A. 2019 Biochar and Soil Physical Health. In Biochar–An Imperative Amendment for Soil and the Environment; Vikas, A., Peeyush, S., Eds.; IntechOpen: London, UK, 2019; Available online: https://www.intechopen.com/books/biochar-an-imperative-amendment-for-soil-and-the-environment/biochar-and-soil-physical-health (accessed on 20 December 2020). [CrossRef] [Green Version]
  61. Quin, P.R.R.; Cowie, A.L.; Flavel, R.J.; Keen, B.P.; Macdonald, L.M.; Morris, S.G.; Singh, B.P.; Young, I.M.; Zwieten, L. Oil mallee biochar improves soil structural properties—A study with x-ray micro-CT. Agric. Ecosyst. Environ. 2014, 191, 142–149. [Google Scholar] [CrossRef]
  62. Zhao, S.X.; Ta, N.; Wang, X.D. Effect of temperature on the structural and physico-chemical properties of biochar with apple tree branches as feedstock material. Energies 2017, 10, 1293. [Google Scholar] [CrossRef] [Green Version]
  63. Hyväluoma, J.; Kulju, S.; Hannula, M.; Wikberg, H.; Källi, A.; Rasa, K. Quantitative characterization of pore structure of several biochars with3D imaging. Environ. Sci. Pollut. Res. 2018, 25, 25648–25658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Ouyang, L.; Wang, F.; Tang, J.; Yu, L.; Zhang, R. Effects of biochar amendment on soil aggregates and hydraulic properties. J. Soil Sci. Plant Nutr. 2013, 13, 991–1002. [Google Scholar] [CrossRef] [Green Version]
  65. Aller, D.; Rathke, S.; Laird, D.; Crusea, R.; Hatfield, J. Impacts of fresh and aged biochars on plant available water and water use efficiency. Geoderma 2017, 307, 114–121. [Google Scholar] [CrossRef]
  66. Mia, S.; Dijkstra, F.A.; Singh, B. Long-term aging of biochar: A molecular understanding with agricultural and environmental implications. Adv. Agron. 2017, 141, 1–51. [Google Scholar]
  67. Günala, E.; Erdem, H.; Çelik, İ. Effects of three different biochars amendment on water retention of silty loam and loamy soils. Agric. Water Manag. 2018, 208, 232–244. [Google Scholar] [CrossRef]
Figure 1. Flow chart summarizing the phases in which the experimental procedure adopted in this work is divided.
Figure 1. Flow chart summarizing the phases in which the experimental procedure adopted in this work is divided.
Agronomy 11 00489 g001
Figure 2. Graphical representation of the soil water capacitive indicators (TAW, PAW, DRW, LAW and RAW) as determined by the volumetric water content (%) corresponding to the characteristic water potentials (SAT, FC, MRX and WP) along the soil adsorption curve. (Read the text for explanations of acronyms and initialisms.)
Figure 2. Graphical representation of the soil water capacitive indicators (TAW, PAW, DRW, LAW and RAW) as determined by the volumetric water content (%) corresponding to the characteristic water potentials (SAT, FC, MRX and WP) along the soil adsorption curve. (Read the text for explanations of acronyms and initialisms.)
Agronomy 11 00489 g002
Figure 3. Linear regression of bulk density (g cm−3) against increasing fractions (%) of biochar in the soil to obtain the experimental soil–biochar mixtures. Red line represents the expected trend based on a simple “dilution” effect, while the blue line regressed the observed data as experimentally observed. The slopes of the two lines are significantly different according to the ANCOVA results. Vertical bars are standard errors.
Figure 3. Linear regression of bulk density (g cm−3) against increasing fractions (%) of biochar in the soil to obtain the experimental soil–biochar mixtures. Red line represents the expected trend based on a simple “dilution” effect, while the blue line regressed the observed data as experimentally observed. The slopes of the two lines are significantly different according to the ANCOVA results. Vertical bars are standard errors.
Agronomy 11 00489 g003
Figure 4. Cluster analysis (average linkage method) applied on data reported in Table 4. A clear split into two clusters (B0–B2–B4 vs. B6–B8–B10) is shown. Moreover, a closer “friendship” between B6 and B10 rather than B8 and B10 can be observed. This particular behavior suggests a possible threshold effect of biochar application to soil. (Read the text for an in-depth explanation).
Figure 4. Cluster analysis (average linkage method) applied on data reported in Table 4. A clear split into two clusters (B0–B2–B4 vs. B6–B8–B10) is shown. Moreover, a closer “friendship” between B6 and B10 rather than B8 and B10 can be observed. This particular behavior suggests a possible threshold effect of biochar application to soil. (Read the text for an in-depth explanation).
Agronomy 11 00489 g004
Figure 5. (A) Desorption curves of the two statistically different treatment clusters (B0, B2 and B4) vs. (B6, B8 and B10) according to the van Genuchten model. (B) pore size distribution functions of the same former curves, obtained by calculating their slope (first derivative dθ/dψ) and plotting the data on an x-axis where soil water potential values were converted into pore equivalent diameter.
Figure 5. (A) Desorption curves of the two statistically different treatment clusters (B0, B2 and B4) vs. (B6, B8 and B10) according to the van Genuchten model. (B) pore size distribution functions of the same former curves, obtained by calculating their slope (first derivative dθ/dψ) and plotting the data on an x-axis where soil water potential values were converted into pore equivalent diameter.
Agronomy 11 00489 g005
Figure 6. Cumulative infiltration (dashed lines) and instantaneous infiltration rate (solid line) derived from the application of the Philip’s model. Biochar treatments (B2, B4, B6, B8 and B10) were grouped because they were not statistically different from teach other. B0 remained statistically separated from the other treatments.
Figure 6. Cumulative infiltration (dashed lines) and instantaneous infiltration rate (solid line) derived from the application of the Philip’s model. Biochar treatments (B2, B4, B6, B8 and B10) were grouped because they were not statistically different from teach other. B0 remained statistically separated from the other treatments.
Agronomy 11 00489 g006
Table 1. Main physicochemical characteristics of soil samples.
Table 1. Main physicochemical characteristics of soil samples.
CharacteristicValue 1
Clay (%)37.09 ± 0.91
Silt (%)36.23 ± 1.08
Sand (%)26.68 ± 0.28
pH (-)8.26 ± 0.06
2 EC (dS m−1)0.47 ± 0.01
P2O5 (mg kg−1)61.40 ± 2.70
Corg (g kg−1)1.44 ± 0.03
Total N (‰)0.60 ± 0.10
C/N (-)2,40 ± 0.48
3 OM (%)2.47 ± 0.05
1 Values are means (n = 3) ± standard errors. 2 EC: Electrical Conductivity. 3 OM: Organic Matter.
Table 2. Main physicochemical characteristics of biochar added to the soil.
Table 2. Main physicochemical characteristics of biochar added to the soil.
CharacteristicValue 1CharacteristicValue 1
pH10.61±0.06Ca (mg kg−1)101,986.40±1791.84
EC (mS m−1)249.00±4.84K (mg kg−1)34,358.80±1622.78
Fixed carbon (% dw)74.79±0.33Mg (mg kg−1)12,789.57±32.75
Volatile solids (% dw)15.31±0.31P (mg kg−1)5479.17±169.99
Ash (% dw)9.90±0.04Fe (mg kg−1)4014.62±365.69
Water content (% a.r.)3.35±0.03Si (mg kg−1)2868.42±586.87
C (% dw)67.69±1.08S (mg kg−1)1914.16±49.53
H (% dw)2.06±0.05Na (mg kg−1)1738.70±181.87
N (% dw)1.01±0.02Cu (mg kg−1)636.74±92.02
S (% dw)0.18±0.01Zn (mg kg−1)240.16±3.29
O (% dw)17.94±1.46Mn (mg kg−1)195.21±12.32
Corg67.01±1.07H/Corg ratio (-)0.37±0.01
C/N66.19±0.28O/Corg ratio (-)0.20±0.02
1 Values are means (n = 3) ± standard errors. dw: dry weight; a.r.: as received.
Table 3. Results of the ANCOVA test performed on the expected and observed bulk density values of the six soil–biochar mixtures.
Table 3. Results of the ANCOVA test performed on the expected and observed bulk density values of the six soil–biochar mixtures.
Source of VariationEstimateStd ErrProb
Intercept1.1910.001<0.0001
Average slope−1.1890.023<0.0001
Expected vs. Observed−0.3190.014<0.0001
Table 4. Capacitive indicators (%) for the six soil–biochar mixtures (treatments).
Table 4. Capacitive indicators (%) for the six soil–biochar mixtures (treatments).
TreatmentTAWDRWPAWRAWLAWPAW/TAWRAW/PAW
B041.95c14.74a27.21b17.48b9.73b65.06b64.19a
B243.82c14.46a29.36b18.72b10.64b67.02b63.76a
B445.01bc15.02a29.99b19.23b10.76b66.70b64.12a
B649.81ab8.93ab40.88a26.99a13.88a82.05a65.97a
B851.00a7.18b43.82a30.29a13.53a85.93a69.10a
B1051.62a7.97b43.65a28.71a14.94a84.73a65.73a
Std Error1.09 1.33 0.83 0.89 0.40 2.30 1.26
CV (%)2.31 11.68 2.32 3.78 3.26 3.06 1.92
Within columns, means followed by the same letter are not significantly different via Tukey HSD at α = 0.05.
Table 5. Results of the “stepwise” procedure applied in performing the parameter fitting of the van Genuchten model to the six soil–biochar mixtures (Bi treatments). Standard error of the estimates, together with their confidence intervals at 95% of probability are also reported.
Table 5. Results of the “stepwise” procedure applied in performing the parameter fitting of the van Genuchten model to the six soil–biochar mixtures (Bi treatments). Standard error of the estimates, together with their confidence intervals at 95% of probability are also reported.
ParameterTreatmentEstimateStd ErrLower Conf.Upper Conf.
θs (%)B0, B2, B469.180.4468.3070.07b
B6, B8, B1075.080.5673.9776.19a
n (-)B0, B2, B41.380.011.351.41b
B6, B8, B101.470.021.441.51a
α (kPa−1)B6, B8, B100.00660.00040.00580.0075n.s.
B0, B2, B40.00570.00040.00500.0064
θr (%)All Bi0.000- - -
In the last column, different letters indicate estimate values significantly different via Tukey HSD at α = 0.05; n.s, not significant.
Table 6. Estimates of the van Genuchten model parameters to the six soil–biochar mixtures (treatments) without the clustering operations applied in the stepwise approach. (Compare Table 6 with Table 5).
Table 6. Estimates of the van Genuchten model parameters to the six soil–biochar mixtures (treatments) without the clustering operations applied in the stepwise approach. (Compare Table 6 with Table 5).
TreatmentθsStd ErrnStd ErrαStd Err
B067.480.801.340.020.00650.0007
B269.560.721.390.020.00550.0006
B470.510.701.410.030.00510.0005
B673.510.761.450.030.00760.0008
B875.171.001.490.030.00670.0007
B1076.560.961.490.030.00560.0006
Table 7. Average infiltration rates (IR cm sec−1) in the soil–biochar mixtures (treatments) at the beginning (first 60 s), at the end (last 60 s) and for the entire duration of the process.
Table 7. Average infiltration rates (IR cm sec−1) in the soil–biochar mixtures (treatments) at the beginning (first 60 s), at the end (last 60 s) and for the entire duration of the process.
TreatmentIR InitialIR FinalIR Total
B00.0381a0.0220a0.0281a
B20.0248b0.0134b0.0158b
B40.0243b0.0124b0.0150b
B60.0225b0.0108b0.0136b
B80.0215b0.0098b0.0125b
B100.0210b0.0105b0.0130b
Std Err0.0013 0.0015 0.0012
CV (%)5.13 11.32 7.35
Within columns, means followed by the same letter are not significantly different via Tukey HSD at α = 0.05.
Table 8. Results of the stepwise procedure in application of the Philip’s cumulative infiltration model to the six soil–biochar mixtures. Biochar treatments (B2, B4, B6, B8 and B10) were grouped because they were not statistically different from each other. B0 remained statistically separated from the other treatments. The standard error of the estimates, together with their confidence intervals at 95% of probability are also reported.
Table 8. Results of the stepwise procedure in application of the Philip’s cumulative infiltration model to the six soil–biochar mixtures. Biochar treatments (B2, B4, B6, B8 and B10) were grouped because they were not statistically different from each other. B0 remained statistically separated from the other treatments. The standard error of the estimates, together with their confidence intervals at 95% of probability are also reported.
Sourse of VariationEstimateStd ErrLower ConfUpper Conf
C1: [X = sqrt (time)]0.15310.00900.13550.1706**
Treat [No Biochar]+0.02490.00900.00740.0425a
Treat [Biochar]−0.02490.0090−0.0425−0.0074b
C2: [X2 = time]0.01160.00080.01010.0131**
Treat [No Biochar]+0.00330.00080.00180.0048a
Treat [Biochar]−0.00330.0008−0.0048−0.0018b
Within columns, means followed by the same letter are not significantly different via Tukey HSD at α = 0.05. **, F test significant at p ≤ 0.01.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Libutti, A.; Francavilla, M.; Monteleone, M. Hydrological Properties of a Clay Loam Soil as Affected by Biochar Application in a Pot Experiment. Agronomy 2021, 11, 489. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11030489

AMA Style

Libutti A, Francavilla M, Monteleone M. Hydrological Properties of a Clay Loam Soil as Affected by Biochar Application in a Pot Experiment. Agronomy. 2021; 11(3):489. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11030489

Chicago/Turabian Style

Libutti, Angela, Matteo Francavilla, and Massimo Monteleone. 2021. "Hydrological Properties of a Clay Loam Soil as Affected by Biochar Application in a Pot Experiment" Agronomy 11, no. 3: 489. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11030489

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop