Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Biochar particle size, shape, and porosity act together to influence soil water properties

  • Zuolin Liu ,

    zl17@rice.edu

    Affiliation Department of Earth Science, Rice University, MS, Houston, Texas, United States of America

  • Brandon Dugan,

    Current address: Department of Geophysics, Colorado School of Mines, Golden, Colorado, United States of America

    Affiliation Department of Earth Science, Rice University, MS, Houston, Texas, United States of America

  • Caroline A. Masiello,

    Affiliations Department of Earth Science, Rice University, MS, Houston, Texas, United States of America, Departments of Chemistry and BioSciences, Rice University, MS, Houston, Texas, United States of America

  • Helge M. Gonnermann

    Affiliation Department of Earth Science, Rice University, MS, Houston, Texas, United States of America

Abstract

Many studies report that, under some circumstances, amending soil with biochar can improve field capacity and plant-available water. However, little is known about the mechanisms that control these improvements, making it challenging to predict when biochar will improve soil water properties. To develop a conceptual model explaining biochar’s effects on soil hydrologic processes, we conducted a series of well constrained laboratory experiments using a sand matrix to test the effects of biochar particle size and porosity on soil water retention curves. We showed that biochar particle size affects soil water storage through changing pore space between particles (interpores) and by adding pores that are part of the biochar (intrapores). We used these experimental results to better understand how biochar intrapores and biochar particle shape control the observed changes in water retention when capillary pressure is the main component of soil water potential. We propose that biochar’s intrapores increase water content of biochar-sand mixtures when soils are drier. When biochar-sand mixtures are wetter, biochar particles’ elongated shape disrupts the packing of grains in the sandy matrix, increasing the volume between grains (interpores) available for water storage. These results imply that biochars with a high intraporosity and irregular shapes will most effectively increase water storage in coarse soils.

Introduction

Biochar is charcoal made for the purpose of soil amendment [1]. Amending soil with biochar is an approach to mitigate climate change [2] and to improve crop productivity [1, 3]. Once mixed with soil, biochar can affect plant growth by altering soil hydrologic properties [47] and nutrient availability [8].

Water movement and storage in soils are crucial for nutrient delivery and plant productivity. Biochar has the potential to alter soil hydrology and to drive shifts in the amount of water stored in soils. To understand how biochar amendment may influence water delivery to plants, we must understand how biochar affects soil hydrologic processes. However, while many studies report effects of specific biochars on specific soil water properties [9, 10], there is a dearth of mechanistic information available, and mechanisms are needed to predict under what circumstances biochar will have beneficial effects on soils.

Understanding how the amount of water held at field capacity (θfc) and at permanent wilting point (θpwp), and the amount of plant available water (θpaw) of soil change with biochar amendment is an efficient way to quantify how biochar affects soil water conditions and plant growth. We use the water retention curve, which defines equilibrium water content (θ) at a given soil water potential (ψ), to extract the key parameters of water content at saturation (θs), θfc (water content at ψ = -33kPa), θpwp (water content at ψ = -1500kPa), and θpaw (= θfc - θpwp) [11, 12]. Water held at field capacity is also defined as the water present in the soil after gravity-driven drainage. Water held at and beyond permanent wilting point is assumed to be held at a pressure too high for plants to extract from soil [13].

Previous studies have shown that biochar increased water retention of soil [14]; however, the mechanisms controlling these observations remain elusive. Sandy soils are a particularly appealing target for biochar amendment because studies on sand and sandy loam often show an increase in θpaw after biochar amendment [10, 1518]. However, few studies focused on the mechanism of how biochar increase θpaw. Without understanding the mechanisms that control biochar-driven changes of water retention of soil, it is difficult to predict when and by how much biochar will improve soil water retention.

Biochar’s particle size, shape, and internal structure likely play important roles in controlling soil water storage because they alter pore characteristics. For instance, biochar has pores inside of particles (intrapores), which may provide additional space for water storage beyond the pore space between particles (interpores) [19]. Particle size may influence both intrapores and interpores through different processes because the size and connectivity of these particles likely differ. In addition, when applied in the field, biochar particles may have different sizes and shapes compared to soil particles. This addition of biochar grains with different shapes and sizes will change interpore characteristics (size, shape, connectivity, and volume) of soil and thus will affect water storage and mobility. For instance, fine biochar particles can fill pores between coarse soil particles, decreasing pore size and changing interpore shape. Conversely, high aspect ratio biochar particles may interfere with packing of low aspect ratio soil grains, leading to increased interpore sizes. Both of these cases can be expected to change soil water retention.

Here we develop a mechanistic understanding of how and when biochar application affects water retention in a sandy matrix. Sandy matrices are a particularly important system to constrain because biochar application has the potential to increase the resilience of agriculture in sandy systems. We determined θs, θfc, θpwp, and θpaw by measuring water retention curves of sand mixed with three particle-size ranges of biochar at 2 wt% (kg biochar/kg total dry mixture x100%). In addition, we conducted control experiments measuring water retention curves of sand plus fine sand (<0.251 mm, volume of fine sand was equal to volume of fine biochar at 2 wt% biochar rate) and sand plus coarse sand (0.853–2.00 mm, volume of coarse sand was equal to volume of coarse biochar at 2 wt% biochar rate). We constrained a suite of physical properties (skeletal density, envelope density, and biochar intraporosity) that influence water retention. Last, we qualitatively and quantitatively characterized the size and shape of biochar particles. Using these constraints, we develop a conceptual model of how these physical properties drive changes in water retention of biochar-sand mixtures.

Materials and methods

Sample preparation

We pyrolyzed mesquite feedstock (2.00–2.30 mm) at 400°C for 4 hours to form biochar (Table 1) as described in Kinney et al. [20] and Liu et al. [21]. Ash content, pH, electrical conductivity, and carbon, nitrogen, and hydrogen content of biochar were reported in Liu et al. [21]. We ground and sieved biochar into three sizes: fine (<0.251 mm, NO. 60 U.S. Std. mesh), medium (0.251–0.853 mm) and coarse (0.853–2.00 mm, NO. 20-NO. 10 U.S. Std. mesh). To obtain accurate mass fractions, all sand and biochar were oven dried at 60°C for 72 hours to remove any water absorbed during storage. We then mixed 2 wt% biochar into the sand (silica sand [Pavestone, US] sieved into 0.251–0.853 mm, NO. 60-NO. 20 U.S. Std. mesh) (Table 1). We created controls by mixing medium sand with fine sand (<0.251 mm, volume of fine sand was equal to volume of fine biochar at 2 wt% biochar rate) and coarse sand (0.853–2.00 mm, volume of coarse sand was equal to volume of coarse biochar at 2 wt% biochar rate).

thumbnail
Table 1. Descriptions of silica sand and biochar used in this study.

https://doi.org/10.1371/journal.pone.0179079.t001

To understand pore systems and water storage, we measured the skeletal density (density of the solids without intrapores, ρs) of biochar and sand by helium pycnometry in a 1cm3 sample chamber (AccuPyc II 1340, Micromeritics, Norcross, GA) and the envelope density of biochar (density including intrapores, ρe) (Geopyc 1360, Micromeritics, Norcross, GA) (Table 2). AccuPyc measures biochar’s skeletal volume by detecting the pressure change due to the change of helium volume that is displaced by biochar’s skeleton. It is assumed that helium molecules penetrate all biochar intrapores. Skeletal density was then obtained using biochar mass divided by its skeletal volume. Geopyc measures biochar’s envelope volume by subtracting volume of a consolidated quasi-fluid composed of small, rigid spheres (DryFlo) from the volume of the same consolidated DryFlo after biochar has been added. Enveloped density is the result of biochar mass divided by the envelope volume. We only measured ρe of biochar before grinding due to instrumental limitation on the minimum particle size measurable. However, grinding biochar into smaller size may result in reduction of intraporosity and thus may affect water storage. For details on measurement of ρs and ρe, see Brewer et al. [22]. We then used these two measurements to calculate porosity of biochars (ϕb = 1 - ρe/ρs). The use of ρs and ρe measurements to calculate porosity compares favorably with porosity determined by mercury (Hg) porosimetry and with N2-sorption based techniques (e.g. BET). Benefits to density-based porosity measurements include ease, speed, low cost, and no involvement of toxic materials. Like Hg porosimetry, density-based measurements detect the entire range of pore sizes in a sample, compared to N2-sorption based techniques, which can only detect very small pores and can miss >90% of the pore volume of biochars [22]. In addition, while Hg porosimetry measures total porosity (pores inside plus pores between particles), the combined total porosity measured through density analysis measures only the porosity inside of particles (intrapores) and does not detect the porosity between particles (interpores). Two disadvantages of density-based techniques include: (a) the inability to make measurements on particles smaller than 2 mm, requiring us to assume that the porosity of small biochar was approximately equal to that of the large biochar; (b) density-based porosity measurements provide only total porosity, unlike Hg porosimetry, which can provide information on the entire spectrum of pore characteristics.

thumbnail
Table 2. Particle size, skeletal density (ρs), and envelope density (ρe) of sand and biochar used in this study.

We report average and standard deviation of at least three measurements.

https://doi.org/10.1371/journal.pone.0179079.t002

Measurement of water retention curves

We measured water retention curves at room temperature (22. 3 ± 0.2°C) using a Hyprop for ψ of +2 to -440 kPa and a WP4C device for ψ of -100 to -300,000 kPa (both pieces of equipment were made by Decagon Devices Inc., Pullman, WA).

For Hyprop measurements, each sample was poured without intentional compaction into a stainless-steel cylinder (2.5 x 10−4 m3 volume, diameter 8cm, height 5cm) with a piece of fabric filter and a plastic cap on the bottom. We then put the cylinder with the sample into a beaker with de-aired, purified water (18.2 MΩ-cm, PURELAB® Ultra Laboratory Water Purification Systems, SIMENS, Germany). The use of purified water allowed us to exclude osmotic potential effects from the measurement. In the beaker, purified water rose from the bottom into the sample through a fabric filter and pushed air out of the sample through the top. We considered samples saturated by this technique after at least 24 hours [23]. We installed two tensiometers with heights of 0.5 and 3.5 cm and the cylinder with sample was clamped to the tensiometer assembly. We then removed the fabric filter and the plastic cap to allow water to evaporate from the sample. We monitored ψ during evaporation using the Hyprop tensiometers and sample mass by a mass balance (Kern EG 2200, Balingen, Germany) [24]. The water retention curve was defined by the average ψ measured by two tensiometers and the water content by volume [θ = (M-Md)V/ρw, m3/m3] calculated from sample mass with water (M, kg), dry sample mass (Md, kg), sample volume (V = 2.5 x 10−4 m3) and water density (ρw = 1000 kg/m3).

For biochar-sand mixtures and sand samples, we measured three water retention curves by the Hyprop (S1 Fig). We reported the average and the standard deviation of these replicates. From the Hyprop data, we extracted θs and θfc and reported the average and standard deviation of these replicated measurements (Table 3). In addition, we measured three replicated water retention curves for fine sand-sand mixtures and coarse sand-sand mixtures using the Hyprop and reported the average and the standard deviation of the replicates.

thumbnail
Table 3. Bulk density (ρb), total porosity (ϕT), and saturated water content (θs) of samples for measuring water retention curves by the Hyprop device.

We report average and standard deviation of at least three measurements.

https://doi.org/10.1371/journal.pone.0179079.t003

For WP4C measurements we added de-aired, purified water to sand and biochar-sand mixtures to make samples with water content from 0.000 m3/m3 to a value that is near the field capacity measured by the Hyprop for each sample. We then placed about 7.5 ml of each sample into a 15 ml stainless steel chamber. The chamber was covered with a plastic cap for 2–3 hours to allow moisture equilibration across the whole chamber. We then inserted the sample chamber into the WP4C, removed the plastic cap, and measured the ψ by a dew point hygrometer. Sample mass was measured to calculate θ [(M-Md)V/ρw, m3/m3]. Based on a suite of experiments, it is difficult to prepare several samples with same ψ. Therefore, we used the WP4C data to estimate one θpwp (without error bars) for our biochar-sand mixtures and sand.

With Hyprop measurements we can determine θfc and WP4C measurements we can determine θpwp, which then allows us to calculate plant available water (θpaw = θfc.- θpwp).

Bulk density and total porosity

We calculated the dry bulk density (ρb = Md/V) of each sample (sand or biochar-sand mixture) (Table 3) using measured, dry sample mass (Md) and total sample volume (V = 2.5 x 10−4 m3), which is the volume of the stainless-steel cylinder for the Hyprop.

We also calculated total porosity (ϕT, volume fraction of intrapores plus interpores) (Table 3) of the biochar-sand mixtures using Eq 1. (1) where Ms (kg) is mass of sand, Mb (kg) is mass of biochar, ρss (kg/m3) is skeletal density of sand, and ρsb (kg/m3) is skeletal density of biochar (Table 2).

Bimodal van Genuchten model

Soil water potential is composed of pressure potential, gravitational potential, osmotic potential, and perhaps by other potential terms [11], although these first three terms are understood to control most systems. Pressure potential (mainly capillary pressure) is a function of soil pore size distribution. Gravitational potential depends on the elevation reference. Osmotic potential is a function of solute concentration. Therefore, in a soil with an absence of solutes, soil water potential is controlled by soil pore size distribution.

To properly describe the water retention characteristics of our biochar-sand mixtures with interpores and intrapores, we used a bimodal van Genuchten model (VGbi) [25] to fit the measured water retention curve for each treatment.

(2)

In Eq 2, θr is the residual water (equal to 0 in our experiments) and θs is the saturated water content (Table 3), k is the number of “pore systems” (i.e. interpores and interpores) that form the total pore size distribution (total water retention curve). In our mixtures with two types of pores (interpores and intrapores), we assume that k is equal to 2, and wi is the weighting factors of each sub soil water retention curve for each pore system where 0 < wi < 1 and Ʃwi = 1. Where αi (>0) is the inverse of the air entry pressure and ni (>1) is a measure of the pore-size distribution the for each sub soil water retention curve, respectively. With these constraints, we simplify Eq 2 into Eq 3.

(3)

The fitted parameters and goodness of fit were determined using the MATLAB Curve Fitting toolbox (Table 4).

thumbnail
Table 4. Bimodal van Genuchten model parameters and goodness of fit include: the weighting factors of soil water retention curve for interpores (w1) and intrapores (w2), the inverse of the air entry pressure for interpores (α1) and intrapores (α2), the measure of the pore-size distribution for interpores (n1) and intrapores (n2), R-square (R2) and root-mean-square error (RMSE).

https://doi.org/10.1371/journal.pone.0179079.t004

Biochar and sand particle shape

To qualitatively examine size and shape of biochar and sand particles, we photographed biochar particles and sand particles under a microscope with a maximum zoom of 1:20 (Stereo Discovery.V20, Zeiss, Germany) (S2 Fig). To quantitatively characterize particle size and shape we used a Camsizer (Retsch Technology, Germany). We measured particle size distribution (S2 Fig) of fine, medium, and coarse biochars and sands and reported median diameter of particles’ shortest chord (D50) and aspect ratio (AR) which is D50 divided by median of the maximum distance between two parallel tangential lines of a particle projection (Table 5). We used these images and the Camsizer data to develop a conceptual model of how biochar particles break and mix with sand particles.

thumbnail
Table 5. Median diameter of particles’ shortest chord (Dmin50) and particles’ aspect ratio (AR defined as Dmin50 divided by Dmax50) of biochar and sand used in this study.

We made measurements through dynamic image analysis (Camsizer, Retsch Technology, Germany).

https://doi.org/10.1371/journal.pone.0179079.t005

Statistical analyses

We used Levene’s test to confirm equality of variances of θs, θfc, ρb, and ϕT between replicated measurements. We then performed statistical comparisons of θs, θfc, ρb, and ϕT between different biochar-sand mixtures and sand by the one-way analysis of variance (ANOVA), followed by Tukey-Kramer’s post hoc test if differences were deemed significant at a p-value less than 0.05. We also computed Pearson correlation coefficients (R) to assess the relationships between θs, θfc, θpwp, θpaw, and ρb as well as the relationships between θs, θfc, θpwp, θpaw, and ϕT.

Results

Biochar grain size played an important role in the water retention of sand-biochar mixtures. Field capacity, permanent wilting point, and plant available water of sand-biochar mixtures all increased with biochar particle size.

The water content at saturation (θs) of fine (0.39 ± 0.03 m3/m3) and coarse (0.37 ± 0.04 m3/m3) biochar-sand mixtures were statistically the same (p>0.05) as that for sand (0.34 ± 0.02 m3/m3); however, the water content at saturation of the medium biochar-sand mixtures (0.41 ± 0.01 m3/m3) was 21% higher (p<0.01) than that for sand (Table 3). While we do not have a definitive explanation for the higher water contents at saturation in the medium biochar-sand mixture, we hypothesize that this is due to changes in packing when particles are combined of similar size, but differing aspect ratio. The differences of θ between biochar-sand mixtures and sand became smaller when ψ became lower. When ψ was less than -5000 kPa, water retention curves of biochar-sand mixtures and sand merged (Fig 1).

thumbnail
Fig 1.

Comparisons of water retention curves (water content, θ, versus soil water potential, ψ) between sand and sand plus 2 wt% (a) fine (b) medium and (c) coarse biochar showed that biochar addition increased water content at given soil water potential. Data indicated with the dots were measured by the Hyprop and the WP4C and data indicated by the lines were fitted by bimodal van Genuchten model (VGbi). We report average and standard deviation of at least three measurements.

https://doi.org/10.1371/journal.pone.0179079.g001

Compared to the θfc of sand (0.025 ± 0.005 m3/m3), the field capacity of medium (0.042 ± 0.002 m3/m3) and coarse (0.050 ± 0.005 m3/m3) biochar-sand mixtures increased by 68% (p<0.01) and 100% (p<0.01), respectively (Fig 2). Field capacity of the fine biochar-sand mixture (0.028 ± 0.001 m3/m3) increased by 12% relative to sand but was not statistically significant (p = 0.25). Similarly, permanent wilting point increased from 0.005 m3/m3 for sand to 0.007 m3/m3 (40% increase), 0.010 m3/m3 (100% increase), and 0.010 m3/m3 (100% increase) for fine, medium, and coarse biochar-sand mixtures, respectively (no p value due to lack of replicates). These increases in θfc and θpwp resulted in increases in θpaw for medium and coarse biochar-sand mixtures. Sand had θpaw of 0.018 ± 0.005 m3/m3, whereas θpaw was 0.021 ± 0.001 m3/m3 for fine biochar-sand mixtures (17% increase, but not different within error). Plant available water was higher for medium and coarse biochar-sand mixtures: 0.032 ± 0.002 m3/m3 (78% increase) and 0.040 ± 0.005 m3/m3 (122% increase) (Fig 2).

thumbnail
Fig 2. Field capacity (θfc), permanent wilting point (θpwp), and plant available water (θpaw) of sand, sand plus 2 wt% fine, medium, and coarse biochar indicated that θfc, θpwp, and θpaw increased with biochar addition as well as biochar particle size.

Values and error bars for θfc were the average and standard deviation of at least three replicates conducted for each treatment. Values for θpwp and θpaw are only one replicate. Error bars of θpaw are the same as error bars of θfc.

https://doi.org/10.1371/journal.pone.0179079.g002

Compared with sand, there were no significant changes of bulk density for fine, medium and coarse biochar-sand mixtures (p = 0.55, 0.22, and 0.08, respectively). Total porosity increased 9.3% (p<0.05) for all biochar-sand mixtures (Table 3).

The addition of fine sand and coarse sand increased water content at higher soil water potential values (Fig 3). This may be the result of changes in particle packing when combining sand particles of differing diameters. As water potential values dropped, the water content of sand, fine sand + sand, and coarse sand + sand mixtures became closer. At saturation the water contents for the fine sand-sand mixture (0.36 ± 0.0 m3/m3) and coarse sand-sand mixture (0.37 ± 0.01 m3/m3) were the same to that of sand alone (0.34 ± 0.02 m3/m3) within error. The water retention curves of fine sand-sand mixtures and coarse sand-sand mixtures continued to overlap with that of the sand sample for all ψ <-1.8kPa (Fig 3).

thumbnail
Fig 3.

(a). Measured water retention curves (water content, θ, versus soil water potential, ψ, measured by the Hyprop) and (b) bimodal van Genuchten model (VGbi) of data from Fig 3A. Sand, fine sand plus sand (volume of fine sand is equal to volume of fine biochar at 2 wt% biochar rate), and coarse sand plus sand (volume of coarse sand is equal to volume of coarse biochar at 2 wt% biochar rate). These three curves overlapped with each other indicating that addition of small fraction of different sizes of sand did not cause significant change in soil water retention at such low rate.

https://doi.org/10.1371/journal.pone.0179079.g003

Compared with sand (bulk density = 1520 ± 20 kg/m3), the bulk density of fine sand-sand mixture (1600 ± 00 kg/m3) and coarse sand-sand mixture (1580 ± 30 kg/m3) were increased (p<0.05) by 5.2% and 4.0%, respectively. Correspondingly, the total porosity of the fine sand-sand mixture (0.40 ± 0.0 m3/m3) and coarse sand-sand mixture (0.43 ± 0.01 m3/m3) were 7.0% and 4.7% lower (p<0.05) than that of sand (0.41 ± 0.01 m3/m3).

Discussion

Our results suggest that the pores inside biochar (intrapores) and the pores created between biochar particles and soil particles (interpores) play fundamentally different roles in soil water retention when capillarity pressure is the main component of soil water potential. Intrapores control water retention at lower soil water potential values causing an increase in field capacity, permanent wilting point, and plant available water for medium and coarse biochar-sand mixtures. However, interpores control water retention at higher soil water potential values for fine biochar-sand mixtures.

Intrapores dominate water retention at lower soil water potential; interpores dominate water retention at higher soil water potential

To understand how pore type and size act to control soil water retention in biochar-sand mixtures, we used a simple calculation to estimate pore diameters that correlate with our observed increases in water retention. We assumed that capillary pressure (Pc) was the major component of soil water potential.

(4)

Here γ is surface tension for the water-air interface at 20°C (equal to 0.072 N/m), θc (°) is contact angle between the water-air interface and biochar/sand surface and d (m) is the pore diameter. Contact angle describes the hydrophilicity or hydrophobicity of a solid surface [26]. The contact angle of biochar surfaces varies with pyrolysis conditions and feedstock type due to the presence of C-H functional groups on the surface of biochar particles [20]. Meanwhile, the measurement of contact angle can be affected by factors like measurement method [27, 28], liquid type [29], particle size [27], and surface morphology [30]. This leads to the complexity of measuring θc of biochar resulting in uncertainties of biochar’s θc. While acknowledging these uncertainties, we assumed θc = 55° for biochar, which is the minimum reported contact angle of fresh biochar in existing studies using direct and indirect methods [7, 2729, 31]. This assumes biochar is hydrophilic, allowing water to penetrate its intrapores. If biochar is hydrophobic (θ >90°), then water entry pressure is positive [32]. In this case, an applied pressure exceeding the entry pressure is needed for water to enter the biochar intrapores. Lack of this external pressure will prevent saturation of biochar intrapores. However, the hydrophobicity of biochar could be reduced by exposure to water [20], as would happen in virtually all environmental conditions, decreasing the contact angle of biochar. Therefore, we assume biochar is hydrophilic in this study, with the understanding that our results are representative of biochar that has had at least some environmental exposure.

Most biochar intrapores have diameters (d) <0.01 mm [33, 34]. Based on these constraints, Eq 4 provides ψ < -16.5 kPa when d is less than 0.01 mm. This suggests that the pores smaller than 0.01 mm control water retention of our biochar-sand mixtures when ψ is less than -16.5 kPa. Given the small size of these pores, we assume this represents intrapores (pores inside biochar particles). We did not use θc = 0° for biochar because there is unlikely that this biochar is fully hydrophillic. However, if θc = 0° for biochar, then ψ is equal to -28.8 kPa which is close to -16.5 kPa considering ψ spans several orders of magnitude.

For a mono-dispersed sand (0.251–0.853 mm), the interpore size (d) is larger than 0.1 mm if we assume d is 40% of particle diameter (>0.251 mm) for packed spheres [35]. We used the contact angle of a hydrophilic sand (θc≈0°) because previous studies showed that biochar only causes a small degree of change in contact angle in soil [27, 28]. Based on this, we then estimated that ψ will be greater than -2.88 kPa when interpore size >0.1 mm. Therefore, interpores would be more likely control water retention curves at higher ψ.

Biochar intrapores increase field capacity, permanent wilting point and plant available water

The high intraporosity (ϕb) of the parent biochar (= 0.6 m3/m3, ϕb = 1-ρe/ρs) (Table 2) suggests that biochar intrapores have the capacity to increase soil water storage, and statistical analyses support this conclusion. As total porosity (ϕT) increased, the amount of water held at a number of pressures (θfc, θpwp, and θpaw) increased for sand amended with medium and coarse biochars. We found positive relationships between θi and ϕT, θfc and ϕT, θpwp and ϕT, θpaw and ϕT (R = 0.63, 0.78, 0.85 and 0.75) (Fig 4). Based on these observations and our calculation that intrapores control water retention at ψ less than -16.5 kPa, we conclude that the increase of θfc (at ψ = -33kPa), θpwp (at ψ = -1500kPa), and θpaw (= θfc- θpwp) (Fig 3) caused by coarser biochar addition is controlled by biochar intrapores (Fig 5).

thumbnail
Fig 4.

Negative Pearson correlation coefficients (R) between bulk density (ρb) and (a) initial water content (θi), (b) field capacity (θfc), (c) permanent wilting point (θpwp), (d) plant available water (θpaw) and positive R between total porosity (ϕT) and (e) initial water content, (f) filed capacity, (g) permanent wilting point, (h) plant available water showed that soil water retention decreased with ρb increase but increased with ϕT increase.

https://doi.org/10.1371/journal.pone.0179079.g004

thumbnail
Fig 5.

Schematic of (a) and (b) sand (dark gray); (c) and (d) sand plus medium biochar (black) on a plot of water retention curves for these two samples. Pores inside of biochar particles were filled with water (light gray) thus increased in water content at saturation as well as field capacity.

https://doi.org/10.1371/journal.pone.0179079.g005

We also found that θfc, θpwp, and θpaw decreased as biochar particle size decreased (Fig 2). We interpret this as the result of destruction of intrapores caused by grinding biochar into smaller particles. This would decrease biochar’s internal porosity, and should be associated with an increase biochar’s skeletal density which was indeed observed (Table 2). As a result, finer biochar-sand mixtures have lower water content at a given soil water potential (Figs 1 and 2).

How much water can biochar intrapores hold? Simple calculations to understand biochar internal water-holding capacity

We used the intraporosity of our parent biochar to calculate a realistic upper estimate of how much water biochar intrapores can hold. We then calculated the increase of water content by biochar intrapores in our experiments. By comparing these two calculations, we can better understand biochar intrapores’ actual role in water retention of soil.

The intraporosity of our parent biochar is 0.6 m3/m3, which means that the parent biochar intrapores can store up to 0.6 m3 water/m3 biochar.

The water content of biochars determined from water retention curves were lower than the water content that parent biochar can store. In section 4.1 we documented that intrapores control water retention of our biochar-sand mixtures when soil water potential (ψ) is less than -16.5 kPa. The water content held by medium and coarse biochars (m3 water/m3 biochar) at -16.5 kPa from measured water retention curves (θb) is showed in Eq 5.

(5)

Where θdiff (m3 water/m3 total mixture) is the difference in the water content at ψ = -16.5 kPa between the biochar-sand mixture and the sand, or water content increased by biochar at ψ = -16.5 kPa; and fvb (m3 biochar/ m3 total mixture) is volume fraction of biochar which can calculated from Eq 6. (6) where Mb (kg) is the biochar mass; ρeb (assumed to be 570 kg/m3 for all biochars, Table 2) is the envelope density of biochar and V (= 2.5 x 10−4 m3) is the sample volume of the mixture. Therefore, by combing Eqs 5 and 6, we obtained Eq 7.

(7)

For medium biochar, θdiff = 0.016 m3/m3 and Mb = 7.5 x 10−3 kg result in a θb of 0.30 m3 water/m3 biochar. Similarly, for coarse biochar, θdiff = 0.027 m3/m3 and Mb = 7.5 x 10−3 kg result in a θb of 0.52 m3 water/m3 biochar. Therefore, the amount of water (0.52 m3 water/m3 biochar) held by coarse biochar intrapores at ψ = -16.5 kPa is slightly less than the maximum water (0.6 m3 water/m3 biochar) that can be stored by parent biochar intrapores. However, water held by medium biochar intrapore at ψ = -16.5 kPa is half of the water content that can be stored by parent biochar intrapores. We interpret this to mean that the decrease in water stored in medium biochar is due to destruction of intrapores caused by grinding biochar into smaller particles.

Biochar particle shape affects interpore volume driving changes in soil water properties

We observed that fine biochar addition to sand increased ϕT (Table 3) as well as water content for ψ greater than -33 kPa (Fig 1). However, there was no significant change of water content when ψ is less than -33 kPa. Based on our interpretation that interpores control water content when ψ is greater than -16.5 kPa when capillarity pressure is the main component of soil water potential, we interpret these results to mean that adding fine biochar into sand increased interpore volume.

Both the size and the shape of biochar particles can impact interpore volume in biochar-sand mixtures. The water retention curves of fine sand-sand mixtures and coarse sand-sand mixtures overlapped with that of sand sample for ψ<-2.6kPa (Fig 3). This indicated that grain size did not play an important role in this soil water potential range. However, the water content of the fine biochar-sand mixture was higher than that of the fine sand-sand mixture (volume of fine sand is equal to volume of fine biochar at 2 wt% biochar rate). Fine biochar particles were more elongated than sand particles as documented microscopic images (S1 Fig) and as quantified by lower AR measured by the Camsizer (Table 5). Through numerical simulation, Deng and Davé [36] found that when elongated particles contacted each other perpendicularly or with angles other than aligned parallel, porosity increased in comparison to when all particles are spheres. Therefore, it is possible that perpendicular contacts between elongated biochar particles and sand particles created more space between particles, causing increased interporosity (Fig 6). As a result, the fine biochar-sand mixture had higher water content than the sand sample at ψ greater than -33 kPa. Because this soil water potential range is above field capacity, the increase of water retention is not likely to increase plant-available water under dry conditions. Instead, it would provide more storage of water on the landscape under wet conditions. For instance, the increase of θi by fine and medium biochar shows that biochar intrapores are likely to hold more water near the surface during a rain event, which may help reduce runoff. Doan et al. [37] reported that presence of biochar significantly reduced water runoff for three years of application in mesocosms. Depending on the scale of application, reduction of runoff could change local rivers’ hydrographs.

thumbnail
Fig 6.

Schematic of (a) and (b) sand (dark gray) plus fine sand and (c) and (d) sand plus fine biochar (black) on a plot of water retention curves for these two samples. Biochar particles are more elongated which creates more pore space when packing. This may increase the distance between particles resulting in increased of interporosity. Sand plus fine biochar had a higher water (light gray shade) content than that of sand plus fine sand at higher soil water potential, probably due to its higher interporosity. However, the two water retention curves merged at lower soil water potential values (less than -33kPa) indicating that the intrapores of the fine biochar does not contribute to soil water retention as discussed in section 4.4.

https://doi.org/10.1371/journal.pone.0179079.g006

Biochar production can be optimized to produce favorable soil water effects

Our results showed that biochar intrapores played an important role in increasing the water retention of sand-biochar mixtures water potentials less than -16.5 kPa. Water retention improvements are most useful when they impact the plant available water. We show here that intraporosity increases plant available water, suggesting that biochar with high intraporosity will be most useful. Feedstock type, pyrolysis temperature, and charring residence time influence biochar’s internal porosity [22]. For instance, biochars with low intraporosity such as wastewater sludge biochar and poultry litter biochar are less favorable for soil water storage at low water potentials (<-16.5 kPa) because their internal porosity is very low [38, 39]. Grass biochar should be better for water storage than wastewater sludge biochar because grass biochar has higher intraporosity than that of wastewater sludge biochar [22]. These interpretations are supported by existing studies. For example, Sun and Lu [40] observed an increase of plant available water by straw biochar and no effect on plant available water by wastewater sludge biochar. In their study, straw biochar increased soil pore volume of pores <10 μm but wastewater sludge biochar did not cause a significant change in soil pore volume in this size range. Higher pyrolysis temperatures produce more porous biochar [22]. Depending on feedstock type, characteristics of biochar intrapores also vary with charring residence time [41]. Therefore, an optimal charring temperature and residence time should be selected to produce biochar with high intraporosity.

The efficiency of biochar for improving soil water retention will be reduced if biochars are hydrophobic, but hydrophobicity can likely be managed by pretreatment. Hydrophobic biochar has positive water entry pressure [32], meaning that an applied force is required for water to enter intrapores. Lack of this external force would prevent water from entering intrapores thus preventing saturation of biochar intrapores and limiting water retention benefits of biochar. Jeffery et al. [31] reported that grass species biochar did not improve soil water retention; this is probably due to its high hydrophobicity (average contact angle of 118°), although it is notable that grass biochar has lower hydrophobicity compared to leaf or wood biochars [20]Biochar’s hydrophobicity varies with production temperature and feedstock [20, 29], but is usually eliminated by brief environmental exposure. Pretreating biochar either by initially wetting it, or by composting, is likely to significantly reduce problems associated with hydrophobicity [20].

Effect of biochar on soil water retention may change over time

Our experiments documented significant increases (up to 127%) in plant-available water after mixing coarser biochar with sand at a laboratory timescale. Over the timescale of field application, biochar particle size, intraporosity, and hydrophobicity might change, likely altering soil water retention. For instance, biochar particle size can be reduced by natural forces such as freeze and thaw cycles [42], plant root penetration [43], and bioturbation [44]. Biochar’s intraporosity can also be reduced by sorption of minerals [45, 46], adsorption of organic matter [47], and microorganism growth [48]. Using x-ray photoelectron spectroscopy, LeCroy et al. [49] found evidence of increased surface oxidation on biochar particles suggesting that the first stage of biochar patina development involves sorption of dissolved organic compounds in soil. In addition, microscopic and pycnometric data from recent field trials point to blockage of biochar intrapores by either organics, minerals, or a combination [46]. Biochar hydrophobicity can prevent water from penetrating into biochar intrapores, prohibiting an improvement of soil water retention [31]. However, Ojeda et al. [27] observed a 69.5% decrease of contact angle of biochar after one year of its addition to soil suggesting that initial biochar hydrophobicity disappeared within one year. This decrease in hydrophobicity will improve soil water retention.

Conclusion

In this study, we used a simple sand-biochar system to develop a mechanistic understanding of how biochar’s internal pores (intrapores) and the pores between biochar and sand particles (interpores) affect soil water retention. In our experiments the addition of biochar to sand increased initial water content and field capacity. Our controlled particle size and porosity conditions allowed the development of conceptual models that connect biochar properties to soil water benefits. We propose that the increase of water retention of sandy soils by biochar addition is caused by biochar intraporosity at lower ψ and by increasing interporosity due to elongated biochar particle shape increasing interpores space at higher ψ when capillarity pressure is the main component of soil water potential. This suggests that to increase plant-available water (θpaw) in sandy soils, biochar with high intraporosity and an irregular shape will be most effective. Various production factors (i.e. feedstock, pyrolysis temperature and charring residence time) may be useful to produce biochar with varying porosity. Further studies are needed to address how long biochar intraporosity, particle size, and particle shape will last after field application.

Supporting information

S1 Fig.

Photomicrograph of (a) fine biochar (<0.251 mm), (b) fine sand (<0.251 mm), (c) medium biochar (0.251–0.853 mm), (d) sand (0.251–0.853 mm), (e) coarse biochar (0.853–2.00 mm), (f) coarse sand (0.853–2.00 mm) and (g) parent biochar (2.00–2.30 mm).

https://doi.org/10.1371/journal.pone.0179079.s001

(TIF)

S2 Fig. Particle diameter (D) distribution (% by volume) of materials used in this study, as measured by a Retsch Camsizer.

Particle diameter (D) is the shortest chord of a particle projection (results close to screening/sieving).

https://doi.org/10.1371/journal.pone.0179079.s002

(TIF)

S1 Table. Nomenclature describes all symbols used in this study.

https://doi.org/10.1371/journal.pone.0179079.s003

(XLS)

Acknowledgments

We thank J. Nittrouer for providing access to the Camsizer and T. Dong for help measuring particle-size distributions. We also benefited from the collaborative support of the Rice University Biochar group.

Author Contributions

  1. Conceptualization: ZL BD CAM.
  2. Data curation: ZL.
  3. Formal analysis: ZL.
  4. Funding acquisition: BD CAM.
  5. Investigation: ZL BD CAM.
  6. Methodology: ZL BD CAM.
  7. Project administration: BD CAM.
  8. Resources: BD CAM HG.
  9. Software: ZL.
  10. Supervision: BD CAM.
  11. Validation: ZL BD CAM HG.
  12. Visualization: ZL.
  13. Writing – original draft: ZL.
  14. Writing – review & editing: ZL BD CAM HG.

References

  1. 1. Lehmann J, Joseph S. Biochar for environmental management: science and technology. UK: Earthscan; 2009.
  2. 2. Woolf D, Amonette JE, Street-Perrott FA, Lehmann J, Joseph S. Sustainable biochar to mitigate global climate change. Nature communications. 2010;1(56):1–9. pmid:20975722
  3. 3. Glaser B, Lehmann J, Zech W. Ameliorating physical and chemical properties of highly weathered soils in the tropics with charcoal-a review. Biol Fertil Soils. 2002;35:219–30.
  4. 4. Githinji L. Effect of biochar application rate on soil physical and hydraulic properties of a sandy loam. Archives of Agronomy and Soil Science. 2013:1–14.
  5. 5. Tryon EH. Effect of Charcoal on Certain Physical, Chemical, and Biological Properties of Forest Soils. Ecological Monographs. 1948;18(1):81–115.
  6. 6. Streubel JD, Collins HP, Garcia-Perez M, Tarara J, Granatstein D, Krguger CE. Influence of contrasting biochar types on five soils at increasing rates of application. Soil Biology & Biochemistry. 2011;75:1402–13.
  7. 7. Herath HMSK, 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(0):188–97.
  8. 8. Major J, Steiner C, Downie A, Lehmann J. Biochar Effects on Nutrient Leaching. Biochar for Environmental Management: Science and Technology. London, UK: Earthscan; 2009. p. 271–88.
  9. 9. Ma N, Zhang L, Zhang Y, Yang L, Yu C, Yin G, et al. Biochar Improves Soil Aggregate Stability and Water Availability in a Mollisol after Three Years of Field Application. PloS one. 2016;11(5):e0154091. pmid:27191160
  10. 10. Omondi MO, Xia X, Nahayo A, Liu X, Korai PK, Pan G. Quantification of biochar effects on soil hydrological properties using meta-analysis of literature data. Geoderma. 2016;274:28–34. https://doi.org/10.1016/j.geoderma.2016.03.029.
  11. 11. Hillel D. Environmental soil physics C.A.: Academic Press; 1998.
  12. 12. Marshall TJ, Holmes JW, Rose CW. Soil physics: Cambridge University Press, UK; 1996.
  13. 13. Brady NC, Weil RR. The nature and properties of soils. 14 ed: Prentice Hall; 2008. 975 p.
  14. 14. Scott HL, Ponsonby D, Atkinson CJ. Biochar: an improver of nutrient and soil water availability—what is the evidence? CAB Reviews. 2014;9(019):1–19.
  15. 15. Briggs C, Breiner JM, Graham RC. Physical and chemical properties of pinus ponderosa charcoal: Implications for soil modification. Soil Sci. 2012;177:263–8. Epub 268.
  16. 16. Brockhoff SR, Christians NE, Killorn RJ, Horton R, Davis DD. Physical and Mineral-Nutrition Properties of Sand-Based Turfgrass Root Zones Amended with Biochar. Agron J. 2010;102(6):1627–31.
  17. 17. Abel S, Peters A, Trinks S, Schonsky H, Facklam M, Wessolek G. Impact of biochar and hydrochar addition on water retention and water repellency of sandy soil. Geoderma. 2013;202:183–91. pmid:WOS:000320423300020.
  18. 18. Hansen V, Hauggaard-Nielsen H, Petersen CT, Mikkelsen TN, Müller-Stöver D. Effects of gasification biochar on plant-available water capacity and plant growth in two contrasting soil types. Soil and Tillage Research. 2016;161:1–9. https://doi.org/10.1016/j.still.2016.03.002.
  19. 19. Masiello C, Dugan B, Brewer C, Spokas K, Novak J, Liu Z. Biochar effects on soil hydrology2015.
  20. 20. Kinney TJ, Masiello CA, Dugan B, Hockaday WC, Dean MR, Zygourakis K, et al. Hydrologic properties of biochars produced at different temperatures. Biomass and Bioenergy. 2012;41:34–43.
  21. 21. Liu Z, Dugan B, Masiello CA, Barnes RT, Gallagher ME, Gonnermann H. Impacts of biochar concentration and particle size on hydraulic conductivity and DOC leaching of biochar–sand mixtures. Journal of Hydrology. 2016;533:461–72.
  22. 22. Brewer CE, Chuang VJ, Masiello CA, Gonnermann H, Gao X, Dugan B, et al. New approaches to measuring biochar density and porosity. Biomass and Bioenergy. 2014;(0).
  23. 23. ASTM International. Standard test methods for determination of the soil water characterisitic curve for desorption using hanging column, pressure exactor, chilled mirror hygrometer, or centrifuge (Standard D6836). West Conshohocken, PA: American Society for Testing and Materials; 2008.
  24. 24. Schindler U, Durner W, von Unold G, Mueller L, Wieland R. The evaporation method: Extending the measurement range of soil hydraulic properties using the air-entry pressure of the ceramic cup. Journal of Plant Nutrition and Soil Science. 2010;173(4):563–72.
  25. 25. Durner W. Hydraulic conductivity estimation for soils with heterogeneous pore structure. WATER RESOURCES RESEARCH. 1994;30(2):211–23.
  26. 26. Letey J, Carrillo MLK, Pang XP. Approaches to characterize the degree of water repellency. Journal of Hydrology. 2000;231–232:61–5.
  27. 27. Ojeda G, Mattana S, Avila A, Alcaniz JM, Volkmann M, Bachmann J. Are soil-water functions affected by biochar application? Geoderma. 2015;249:1–11. pmid:WOS:000352671000001.
  28. 28. Smetanova A, Dotterweich M, Diehl D, Ulrich U, Fohrer N. Influence of biochar and terra preta substrates on wettability and erodibility of soils. Zeitschrift Fur Geomorphologie. 2012;57:111–34. pmid:WOS:000317328200008.
  29. 29. Wang D, Zhang W, Hao X, Zhou D. Transport of biochar particles in saturated granular media: Effects of Pyrolysis temperature and particle size. Environ Sci Technol. 2013;47:821–8. pmid:23249307
  30. 30. Shang J, Flury M, Harsh JB, Zollars RL. Comparison of different methods to measure contact angles of soil colloids. Journal of Colloid and Interface Science. 2008;328(2):299–307. pmid:18930239
  31. 31. Jeffery S, Meinders MBJ, Stoof CR, Bezemer TM, van de Voorde TFJ, Mommer L, et al. Biochar application does not improve the soil hydrological function of a sandy soil. Geoderma. 2015;251–252:47–54.
  32. 32. Wang Z, Wu L, Wu QJ. Water-entry value as an alternative indicator of soil water-repellency and wettability. Journal of Hydrology. 2000;231:76–83. pmid:WOS:000087736400007.
  33. 33. Hardie M, Clothier B, Bound S, Oliver G, Close D. Does biochar influence soil physical properties and soil water availability? Plant and Soil. 2014;376(1–2):347–61.
  34. 34. Sun H, Hockaday WC, Masiello CA, Zygourakis K. Multiple controls on the chemical and physical structure of biochars. Ind Eng Chem Res. 2012;51:3587–97.
  35. 35. Villaume JF. Investigations at Sites Contaminated with Dense, Non-Aqueous Phase Liquids (NAPLs). Ground Water Monitoring & Remediation. 1985;5(2):60–74.
  36. 36. Deng XL, Dave RN. Dynamic simulation of particle packing influenced by size, aspect ratio and surface energy. Granular Matter. 2013;15(4):401–15. pmid:WOS:000323280500003.
  37. 37. Doan TT, Henry-des-Tureaux T, Rumpel C, Janeau J- L, Jouquet P. Impact of compost, vermicompost and biochar on soil fertility, maize yield and soil erosion in Northern Vietnam: A three year mesocosm experiment. Science of The Total Environment. 2015;514:147–54. pmid:25659313
  38. 38. Novak JM, Lima I, Xing B, Gaskin JW, Steiner C, Das KC, et al. Characterization of Designer Biochar Produced at Different Temperatures and Their Effects on a Loamy Sand. 2009. 2009;3. Epub 2014-08-20.
  39. 39. Waqas M, Li G, Khan S, Shamshad I, Reid BJ, Qamar Z, et al. Application of sewage sludge and sewage sludge biochar to reduce polycyclic aromatic hydrocarbons (PAH) and potentially toxic elements (PTE) accumulation in tomato. Environmental Science and Pollution Research. 2015;22(16):12114–23. pmid:25877899
  40. 40. Sun FF, Lu SG. Biochars improve aggregate stability, water retention, and pore- space properties of clayey soil. Journal of Plant Nutrition and Soil Science. 2014;177(1):26–33. pmid:WOS:000331950400003.
  41. 41. Wang Y, Hu Y, Zhao X, Wang S, Xing G. Comparisons of Biochar Properties from Wood Material and Crop Residues at Different Temperatures and Residence Times. Energy & Fuels. 2013;27(10):5890–9.
  42. 42. Carcaillet C. Are Holocene wood-charcoal fragments stratified in alpine and subalpine soils? Evidence from the Alps based on AMS 14C dates. The Holocene. 2001;11(2):231–42.
  43. 43. Lehmann J, Pereira da Silva J Jr., Steiner C, Nehls T, Zech W, Glaser B. Nutrient availability and leaching in an archaeological Anthrosol and a Ferralsol of the Central Amazon basin: fertilizer, manure and charcoal amendments. Plant and Soil. 2003;249(2):343–57.
  44. 44. Ponge JF, Topoliantz S, Ballof S, Rossi JP, Lavelle P, Betsch JM, et al. Ingestion of charcoal by the Amazonian earthworm Pontoscolex corethrurus: A potential for tropical soil fertility. Soil Biology & Biochemistry. 2006;38(7):2008–9. pmid:WOS:000239946100053.
  45. 45. Lehmann J. Bio-energy in the black. Frontiers in Ecology and the Environment. 2007;5(7):381–7.
  46. 46. Sorrenti G, Masiello CA, Dugan B, Toselli M. Biochar physico-chemical properties as affected by environmental exposure. Science of The Total Environment. 2016;563:237–46. pmid:27135586
  47. 47. Kaiser K, Guggenberger G. The role of DOM sorption to mineral surfaces in the preservation of organic matter in soils. Organic Geochemistry. 2000;31(7–8):711–25.
  48. 48. Thies JE, Rillig MC. Characteristics of biochar: biological properties. Biochar for environmental management: Science and technology. 2009:85–105.
  49. 49. LeCroy C, Masiello CA, Rudgers JA, Hockaday WC, Silberg JJ. Nitrogen, biochar, and mycorrhizae: Alteration of the symbiosis and oxidation of the char surface. Soil Biology & Biochemistry. 2013;58:248–54. pmid:WOS:000317158300030.