Physicochemical Changes in Loam Soils Amended with Bamboo Biochar and Their Inﬂuence in Tomato Production Yield

: Soil degradation and water stress in Costa Rica challenge the production of highly sensitive crops. This work is aimed at evaluating the physical and chemical changes in sandy loam (SL) and a silt loam (SiL) soil when amended with bamboo biochar while estimating the enhancement of tomato productivity. Biochar, obtained from Guadua Angustifolia bamboo feedstock, was mixed into sieved bulk soil substrate from the topsoil, from Andosol and Umbrisol groups, at application rates of 1, 2.5, and 5% (dry mass). Physicochemical and morphological properties of biochar such as pH, hydrophobicity, scanning electron microscopy images, helium picnometry, speciﬁc surface area by the Brunauer–Emmett–Teller (BET) method, CHNS, and ash content were determined. Soil hydrophobicity, acidity, electrical conductivity, cation exchange capacity and water retention, available water content, and air capacity were analyzed for the amended soils. Tomato yield was quantiﬁed after a harvest period of two months. The admixture of biochar did not signiﬁcantly increase soil cation exchange capacity but increased water retention in the range of available water content. Class A (>200 g) tomato yield increased 350% in the SL and 151% in the SiL. Class B (100–200 g) tomato yields increased 27% in the SL but decreased about 30% in the SiL. Tomato yield response seems attributable to variation of water retention capacity, available water content, and air capacity. These results support the use of adapted water management strategies for tomato production based on soil physical changes of biochar.


Introduction
Mankind is facing the challenge of feeding an increasing population by producing more food with optimal use of natural resources. The preservation of the quality of soils used for agricultural production is fundamental for this purpose, but climate change, water scarcity, and soil degradation have emerged as limiting factors for crop productivity [1]. Soil degradation in agricultural land is highly related to the depletion of its physical, chemical, and hydrological properties.
Tomato is one of the crops with the highest greenhouse production in Costa Rica. It is cultivated throughout the year; often in soils of the group Andosols [2] of volcanic origin,

Materials and Methods
The research was conducted in a greenhouse from September 2018 to February 2019 at the School of Agricultural Engineering at the Costa Rican Institute of Technology (ITCR) in Cartago, Costa Rica (9 • 51 26.09 N, 83 • 54 44.33 W). The study area was located 1470 m above sea level, with a humid tropical climate. During the year, temperature oscillates between 16 and 25 • C. The average annual precipitation is 1561 mm and the average temperature inside the greenhouse was 20 • C, with maximum recording temperatures of 30 • C during the sunny hours of the day and minimum temperatures of 12 • C at night.

Soil and Biochar Production
Two soil materials were collected from the topsoil (0-30 cm) in localities of the Central Valley in Costa Rica. We used soils from the groups Umbrisol and Andosol [2] with silt loam (SiL) and sandy loam (SL) textures, respectively. They were air dried, passed through a 2-mm sieve, and thoroughly mixed, to achieve homogeneity.
Biochar was obtained from bamboo stems (culms) of Guadua Angustifolia Kunth. The bamboo was collected from a plantation on the north-western side of Lake Arenal, Costa Rica. The bamboo culms were split by half ( Figure 1a) and air-dried in a greenhouse for over 30 days. The bamboo was pyrolyzed in absence of air, in sealed metallic kilns, utilizing the pyrolytic gases as fuel (Figure 1b) by conducting them via steel piping to a burner beneath the biomass compartment. These kilns were developed at the Costa Rica Institute of Technology (ITCR), and are described by Pérez Martínez [28]. The yield of biochar was generally 28% to 35%, relative to weight of air-dried material; if the weight of the firewood is considered, the yield was about 18% to 25%. The temperature inside the pyrolysis chamber of the kilns was measured using K-type thermocouples. During pyrolysis, the system reached temperatures over 400 • C and remained constant for 15 to 30 min, after which the kilns and their contents were left to cool for at least 3 h before opening. The cool Agronomy 2021, 11,2052 3 of 17 biochar was stored in airtight containers for several weeks to prevent spontaneous ignition, and later crushed and passed through a 3.5-mm sieve.
Agronomy 2021, 11, x FOR PEER REVIEW 3 of 17 reached temperatures over 400°C and remained constant for 15 to 30 min, after which the kilns and their contents were left to cool for at least 3 h before opening. The cool biochar was stored in airtight containers for several weeks to prevent spontaneous ignition, and later crushed and passed through a 3.5-mm sieve.
(a) (b) Figure 1. Carbonization kilns designed at the Costa Rica Institute of Technology for biochar production: (a) pyrolysis chambers of three kilns (behind) and halved bamboo stems (at left); (b) reactor during production.

Soil-Biochar Mixture Preparation
Soil-biochar mixtures were prepared by adding biochar at three application rates of 1, 2.5, and 5% (mass percent) to the two dry sieved soil materials, corresponding to equivalents of 10, 26, and 52 ton/ha, respectively, based on assumptions and application rates similar to those made in other studies [5,27,29,30]. Application rates were calculated considering a depth of 15 cm. These mixing rates were chosen to obtain significant results at the lowest possible cost for a farm; higher biochar contents tend to be dismissed by farmers as too costly. The treatments were set as: SL0 (control-unamended sandy loam), SL1 (sandy loam + 1% biochar), SL2.5 (sandy loam + 2.5% biochar), SL5 (sandy loam + 5%

Soil-Biochar Mixture Preparation
Soil-biochar mixtures were prepared by adding biochar at three application rates of 1, 2.5, and 5% (mass percent) to the two dry sieved soil materials, corresponding to equivalents of 10, 26, and 52 ton/ha, respectively, based on assumptions and application rates similar to those made in other studies [5,27,29,30]. Application rates were calculated considering a depth of 15 cm. These mixing rates were chosen to obtain significant results at the lowest possible cost for a farm; higher biochar contents tend to be dismissed by farmers as too costly. The treatments were set as: SL0 (control-unamended sandy loam), SL1 (sandy loam + 1% biochar), SL2.5 (sandy loam + 2.5% biochar), SL5 (sandy loam + 5% biochar), SiL0 (control-unamended silt loam), SiL1 (silt loam + 1% biochar), SiL2.5 (silt loam + 2.5% biochar), and SiL5 (silt loam + 5% biochar). Each treatment combination was replicated seven times, packed in truncated cone pots with 27.6 cm diameter at the bottom, 32.0 cm diameter at the top, with 15 cm height, and allocated in a completely random design in a greenhouse as shown in Figure 2. Each pot was filled with 8.4 kg of an air-dried soil mixture and packed to an initial soil density of 0.7 g/cm 3 . All pots were watered by drip irrigation for one week; then, one tomato seedling of the variety Mountain Fresh Plus was transplanted in each pot. Equal drip irrigation (2 L/h) was applied to all treatments three times a week. On irrigation days, each pot was watered for 10 min at 6 am, 3 pm, and 10 pm, in order to minimize evaporation and filtration losses.

Physicochemical Analysis of the Materials
Soil texture was determined using a modified Bouyoucos hydrometer method, by a combined sieve (>0.063 mm) [32] and the automated particle size analyzer (PARIO) from METER Group (Pullman, WA, USA) [33]. Soil organic matter was removed previous to texture analysis [34]. Soil classification was based on the FAO guidelines [2]. Particle size distribution of the bamboo biochar was determined from a representative sample of 55 g by wet sieve analysis using 4, 2, 1, 0.500, 0.250, 0.125, and 0.063 mm mesh size. The pH was measured in water using a 1:2 ratio (soil: water). Cation exchange capacity (CEC) was determined by the unbuffered salt extraction method [35]. Effective cation exchange capacity (ECEC) was assessed by the sum of bases, (SB = Ca 2+ + Mg 2+ + K + ) and acidity using ammonium acetate pH 7.0 [35]. Sodium (Na + ) analysis was not included in the ECEC, since it is negligible in Costa Rican soils due to high precipitation rates. Electrical conductivity (EC) was measured with a water conductivity meter (Hanna Instruments HI98312, R.I, USA) by adding 50 mL of deionized water to 10 g of soil sample, followed by thorough agitation for 5 min and rest for 30 min. Then, the mixture was filtered, and EC was measured in the extracted solution. Acidity was determined as the subtraction of ECEC-SB. Initially, the sandy loam (SL) was acidic (pH 4.3). To avoid the interference of this factor on tomato productivity, liming was performed by using quick lime (CaO) to level the pH to 5.0. SL was showing 6.6 cmol(+)/L acidity, from which 3.1 cmol(+)/L corresponded to exchangeable aluminum. After pH correction of SL with quick lime (CaO 105.5% calcium carbonate equivalent), the acidity was 0.74 cmol(+)/L, which is marginally above 0.5 cmol(+)/L, the threshold value where aluminum toxicity effects in plants are negligible [31]. Contrarily, SiL was not corrected with CaO, since the exchangeable aluminum was <0.16 cmol(+)/L, and the acidity reported was between 0.13 to 0.16 cmol(+)/L.
The same nutrient schedule was applied to all pots during a period of 16 weeks. Based on the tomato absorption curve and 56 plants, we applied two water soluble compound fertilizers NPK 12:42:12 and 10:5:25 with total doses of 0.08 and 7.88 kg, respectively. Furthermore, we added water soluble magnesium sulfate, potassium sulfate, and calcium nitrate with total doses of 2.70, 1.48, and 6.13 kg, respectively. Fertilizer doses were adjusted every 4 weeks depending on the plant requirements (Table 1).
Equal drip irrigation (2 L/h) was applied to all treatments three times a week. On irrigation days, each pot was watered for 10 min at 6 am, 3 pm, and 10 pm, in order to minimize evaporation and filtration losses. The four-week dose is calculated by multiplying the weekly dose by 4.

Physicochemical Analysis of the Materials
Soil texture was determined using a modified Bouyoucos hydrometer method, by a combined sieve (>0.063 mm) [32] and the automated particle size analyzer (PARIO) from METER Group (Pullman, WA, USA) [33]. Soil organic matter was removed previous to texture analysis [34]. Soil classification was based on the FAO guidelines [2]. Particle size distribution of the bamboo biochar was determined from a representative sample of 55 g by wet sieve analysis using 4, 2, 1, 0.500, 0.250, 0.125, and 0.063 mm mesh size. The pH was measured in water using a 1:2 ratio (soil: water). Cation exchange capacity (CEC) was determined by the unbuffered salt extraction method [35]. Effective cation exchange capacity (ECEC) was assessed by the sum of bases, (SB = Ca 2+ + Mg 2+ + K + ) and acidity using ammonium acetate pH 7.0 [35]. Sodium (Na + ) analysis was not included in the ECEC, since it is negligible in Costa Rican soils due to high precipitation rates. Electrical conductivity (EC) was measured with a water conductivity meter (Hanna Instruments HI98312, Smithfield, RI, USA) by adding 50 mL of deionized water to 10 g of soil sample, followed by thorough agitation for 5 min and rest for 30 min. Then, the mixture was filtered, and EC was measured in the extracted solution. Acidity was determined as the subtraction of ECEC-SB. Acid saturation percentage (AS%) was determined from the ratio of acidity to ECEC. Soil organic matter was quantified with the Dumas dry combustion method [36].
Elemental analysis of C, N, H, and S was carried out in biochar with a Vario MACRO cube analyzer, Elementar Americas Inc (Langenselbold, Germany). Oxygen content was determined by the mass balance of the previous elements, including ash content. Highly volatile, medium volatile matter, combustible material, and ash content of the biochar were determined according to the ASTM E1131-08 standard, by thermogravimetric analysis (TGA) with a TA SDT Q600, TA Instruments (New Castle, DE, USA).
Biochar pore structure was characterized by the scanning electron microscopy (SEM), electron microscope, TM 3000, Hitachi Tabletop Microscope (Tokyo, Japan). Specific surface area of the biochar was analyzed by the Brunauer-Emmett-Teller (BET) method, that uses a low-temperature N 2 adsorption-desorption analysis ASAP 2000 equipment, Micromeritics Inc. (Norcross, GA, USA) [37]. The skeletal density was determined by gas pycnometry using helium (Quantachrome Inc., Boynton Beach, FL, USA) at room temperature (25 • C) and 1.01 bar. Skeletal density values were obtained from five replicates and expressed as mean ± standard deviation. Biochar bulk density was determined from seven replicates, using a 100 mL container. Water repellency of the biochar was determined by the sessile drop method described by Masís-Meléndez et al. [37], and the molarity of ethanol droplet (MED) method was used for the soil and biochar mixtures.

Soil Water Holding Capacity Measurements
Core samples of ca. 30 cm 3 were taken from each pot, saturated from the bottom by capillary rise for 24 h and thereafter drained to matric potentials of −0.3, −1.0, −3.0, and −6.0 kPa by the hanging water column method, while samples for −130 kPa were desiccated in a Richards pressure plate apparatus. The weight of each sample was recorded after each matric potential to obtain the gravimetric water content. After the last matric potential, samples were oven-dried at 105 • C for 16 h. Volumetric water content was determined by multiplying the gravimetric water content by the bulk density of each sample. The gravimetric water content at the matric potential of −1500 kPa was determined in duplicates of bulk soil using a soil water potential meter (WP4C) from METER Group (Pullman, WA, USA), following the ASTM D6836 standard. Water retention curves were fitted with the closed-form Equation (1) [38]: where ϑ s and ϑ r are the saturated and residual water content (cm 3 /cm 3 ), respectively; ψ is the matric potential (kPa); α is the inverse of the air-entry value (1/cm); n and m are independent empirical parameters. Available water capacity (AWC) was calculated as a difference of the water content at a matric potential of −6 kPa (i.e., field capacity, FC, cm 3 /cm 3 ) minus the water content at a matric potential of −1500 kPa (i.e., permanent wilting point, PWP, cm 3 /cm 3 ) [39]. Air capacity (AC) was determined as the difference between the water content at saturation and the water content at a matric potential of −6 kPa. Pore size was defined as macropores (>−6 kPa) based on the cylindrical capillaries with an equivalent diameter greater than 50 µm, capillary pores comprising narrow coarse pores (−6 kPa to −30 kPa) with diameters between 50 and 10 µm, medium pores (−30 kPa to −1500 kPa) with diameters between 10 and 0.2 µm, and micropores (<−1500 kPa) with diameters <0.2 µm [40].
Hydraulic conductivity was determined with a minidisk infiltrometer model S from METER (Pullman, USA) by measuring infiltration with a suction of 1 cm H 2 O = 1 hPa.

Soil Water Storage and Air Capacity Indexes
To assess soil physical quality, we used the soil water storage capacity index (SWSCI) and air capacity index (ACI), proposed by [41], as shown in Equations (2) and (3): where TP corresponded to the total porosity, equivalent to the saturated water content obtained from the fitted soil water retention curve. The study from [41] suggested that optimal soil aeration occurs when 66% of the pore space is filled with water. Thus, the soil physical quality criteria to improve root growth was set as SWSCI = 0.66 and ACI = 0.34.

Fruit Yield and Classification
During two months (from 24 December 2018 to 22 February 2019), fruits were weekly harvested to determine fresh fruit weight, the number of fruits per plant, and average fruit diameter and height. This harvesting period was considered the best production time for this variety of tomato. Fruit shape was measured with a Vernier caliper. Fruits were graded into three classes: Class A tomatoes weighing more than 200 g; Class B tomatoes weighing between 100-200 g; Class C tomatoes weighing less than 100 g. These classes approximate to the market classes used in Costa Rica for tomatoes. Yield components were reported as average fruit weight per plant per class and average number of fruits per plant per class, from the tomatoes harvested in two months.

Statistical Analysis
For all parameters, mean and standard deviation were calculated. The one-factor analysis of variance (ANOVA) with p < 0.05 was used to assess differences in soil water storage and air capacity indexes, soil chemical properties, and tomato productivity per treatment. The differences between means were assessed using the HSD Tukey s test (p < 0.05). The Shapiro-Wilk test (p < 0.05) was used to determine the normality of the data. All the statistical analyses were performed with the freeware R version 3.6.1 using the Rcmdr package [42].

Characterization of the Biochar and Soil Amendments
Elemental analysis, specific surface area, skeletal density, and pH of the biochar used, with and without sieving, are listed in Table 2. It was found that total carbon (C) and nitrogen (N) content were mostly homogeneous in most particle sizes, except for 1.00 and 0.125 mm for C and 0.125 and 0.063 mm for N. Nevertheless, most of the particle sizes for C, H, O, N, and S did not show any significant difference. As expected, ash content was significantly higher at 0.063 mm biochar particle size (BPS). C:N ratio was mostly higher at bigger BPS, ranging from 0.5 to >2 mm. Likewise, O:C and H:C ratios remained homogeneous along with all BPS. Skeletal density (SD) values in the BPS range were homogeneous, suggesting that any porosity changes were due to the biochar internal porosity. Table 2. Elemental analysis, specific surface area (SSA), and skeletal density (SD) of the biochar after sieving and granulometric classification with biochar particle size (BPS) distribution ranging 3.5 to 0.063 mm. The BPS distribution in Figure 3 shows that 85% (>D 15 ) of the biochar particles were greater than those of the soil materials, thus increasing the fraction of larger particles in the treatments with biochar. On average, the BPS was seven times greater than that of SL (>D 15 ) and ten times larger than that of SiL. Below D 15 (i.e., 15% particles smaller than 0.0035 mm), corresponding to the range of silt and sand fractions, the BPS was smaller than that of the soil materials. The average diameter (D 50 ) in biochar was 0.41 mm, whereas in SL and SiL the average diameter was 0.07 and 0.03 mm, respectively. Figure 4 shows SEM images of the particle fractions of 0.063, 0.125, 0.250, 0.500, 1.0, and >2 mm (up to 3.5 mm). All particle sizes show a similar pore structure, with mainly tubular, parallel pores. The array of pores seen in Figure 4d-f is associated to fiber sclerenchyma material surrounding bundles of xylem and phloem vessels [43], with pores of about 10 µm and smaller (medium and fine pores) loosely distributed throughout an otherwise solid material surrounding the largest pores (see cross-sections in Figure 4b,c). The largest pores (>50 µm) are associated with the xylem, while smaller pores may be associated with phloem sieve vessels (Figure 4c,e), and parenchymal tissue, as seen in Figure 4a [44,45]. Xylem pores are continuous; phloem pores may be continuous or semi-continuous, and those in parenchymal tissue are discontinuous (Figure 4a,d). All size fractions contain nearly identical amounts of parallelly oriented pores, mostly about 30 µm to about 1 µm wide.
The biochar originating from parenchymal tissue (Figure 4a,d) clearly shows thinner walls between its discontinuous pores, and thus could have lower mechanical resistance and be the main source of very fine particles. The specific surface areas (SSA) of the biochar fractions were low (Table 2), indicating low microporosity (pore sizes below 50 nm). As seen in Figure 3, about 85% of this biochar was in a range over 3 µm; thus, its predominant contribution to porosity was between 1 µm and 30 µm diameters. Biochar maximum particle size of 3.5 mm was preferred to smaller sizes to minimize channel blockage by small biochar particles, promote water infiltration and solute transport, and preserve the effects of added macroporosity. Overall, the skeletal density values obtained were in the 1.4 to 1.6 g/cm 3 range, which is characteristic of biochar produced at low temperature. Likewise, bulk density of the biochar was between 0.20 and 0.34 g/cm 3 , which is typical of softwood biomass.   Figure 4 shows SEM images of the particle fractions of 0.063, 0.125, 0.250, 0.500, 1.0, and >2 mm (up to 3.5 mm). All particle sizes show a similar pore structure, with mainly tubular, parallel pores. The array of pores seen in Figure 4d-f is associated to fiber sclerenchyma material surrounding bundles of xylem and phloem vessels [43], with pores of about 10 µm and smaller (medium and fine pores) loosely distributed throughout an otherwise solid material surrounding the largest pores (see cross-sections in Figure 4b,c). The largest pores (>50 µm) are associated with the xylem, while smaller pores may be associated with phloem sieve vessels (Figure 4c,e), and parenchymal tissue, as seen in Figure 4a [44,45]. Xylem pores are continuous; phloem pores may be continuous or semi-continuous, and those in parenchymal tissue are discontinuous (Figure 4a,d). All size fractions contain nearly identical amounts of parallelly oriented pores, mostly about 30 µm to about 1 µm wide. The biochar originating from parenchymal tissue (Figure 4a,d) clearly shows thinner walls between its discontinuous pores, and thus could have lower mechanical resistance and be the main source of very fine particles. The specific surface areas (SSA) of the biochar fractions were low (Table 2), indicating low microporosity (pore sizes below 50 nm). As seen in Figure 3, about 85% of this biochar was in a range over 3 µm; thus, its predominant contribution to porosity was between 1 µm and 30 µm diameters. Biochar maximum particle size of 3.5 mm was preferred to smaller sizes to minimize channel blockage by small biochar particles, promote water infiltration and solute transport, and preserve the effects of added macroporosity. Overall, the skeletal density values obtained were in the 1.4 to 1.6 g/cm 3 range, which is characteristic of biochar produced at low temperature. Likewise, bulk density of the biochar was between 0.20 and 0.34 g/cm 3 , which is typical of softwood biomass. Chemical and physical properties of the silt loam (SiL) and the sandy loam (SL) treatmets, shown in Table 3, were analyzed after harvesting. It was found that the addition of biochar produced maximum increases of 0.2 pH units in the silt loam (SiL) and 0.6 pH units in the sandy loam (SL) treatments. SL mixtures showed a medium level of acidity. Acid saturation percentage (AS%), which is referred to as the availability of the exchangeable extractable aluminum concerning the exchangeable base cations, suggests a higher presence of soluble aluminum in SL mixtures than in SiL, which decreased with rising contents of biochar, although all values were within the optimum range below 10%.
In the amended sandy loam, the addition of biochar reduced salinity, although all EC values were below 2 mS/cm. We also found no significant changes (p > 0.05) in CEC as biochar content increased.
AWC decreased significantly, as biochar application rate increased, in SiL2.5 (25%) and SiL5 (21%) concerning the control (SiL0), whereas in the SL mixtures AWC increased but not significantly as the biochar dose increased (9% in SL2.5 and 14% in SL5 with respect to the SL0 control). AC did not register any significant change in the SiL mixtures, whereas biochar produced a significant increase in SL mixtures. Soil water retention at FC (−6 kPa) increased 11% in SiL5 and 9% in SL5 with respect to their corresponding control treatments. Although, initial bulk density was set to 0.7 g/cm 3 , final bulk density increased 14% for all treatments, independently from the biochar application rate. The impact on bulk density was negligible, since the differences between treatments were not statistically significant.  Chemical and physical properties of the silt loam (SiL) and the sandy loam (SL) treatmets, shown in Table 3, were analyzed after harvesting. It was found that the addition of biochar produced maximum increases of 0.2 pH units in the silt loam (SiL) and 0.6 pH units in the sandy loam (SL) treatments. SL mixtures showed a medium level of acidity. Acid saturation percentage (AS%), which is referred to as the availability of the exchangeable extractable aluminum concerning the exchangeable base cations, suggests a higher

Effect of Biochar on the Hydraulic Properties
The sessile droplet method and the MED test showed no evidence of hydrophobicity in the bamboo biochar, in the soils without biochar, nor in their mixtures. Therefore, water repellency was excluded from the analysis as a possible mechanism influencing the wettability of the amended soils. The addition of biochar produced opposite effects on water retention, AWC, AC, and pore size distribution, in the two soil types. Soil water retention increased with increasing biochar contents in both soils; however, this did not necessarily lead to an increase in AWC ( Figure 5). In the SiL treatments, a biochar application rate of 5% produced a significant increase in soil water retention, whereas AC decreased. Additionally, AWC decreased with higher biochar fractions, i.e., 25% and 21% for SiL2.5 and SiL5, respectively; compared to the unamended treatment. Contrarily, in the SL treatments, AC increased notoriously, probably due to the addition of biochar with a larger particle size than the soil. The increase of capillary pores led to higher water retention values at matric potentials between −6 kPa and −10 kPa (9% and 26% in SL2.5 and SL5, respectively). In addition, no significant water retention changes were observed at a matric potential close to −1500 kPa, consequently producing AWC increments of 15% in SL2.5 and 17% in SL5.
The parameters used to fit the soil water retention curves of the soil materials with the van Genuchten-Mualem model are shown in Table 4. The saturated hydraulic conductivity (K s ) increased in most of the soil mixtures, with the addition of biochar. However, we observed a higher increase of K s in the sandy loam, which coincides with larger pore formation at less negative matric potentials ( Figure 5).

Soil Physical Quality Indicators
The Soil Water Storage Index (SWSCI) and Air Capacity Index (ACI) based on [41] are shown in Figure 6. The SWSCI index defines the water-filled pore space compared to the total porosity while the ACI index indicates the aeration capacity of the soil. It was found that the unamended silt loam had a SWCI index over the recommended optimal value of 0.66, and none of the biochar application rates enhanced the SWCI or ACI indexes. Contrarily, the unamended sandy loam was close to the optimum SWCI index and biochar fractions of 1% and 5% increased the ACI index, corresponding to the formation of macropores.
Agronomy 2021, 11, x FOR PEER REVIEW 11 matric potential close to −1500 kPa, consequently producing AWC increments of 15% SL2.5 and 17% in SL5. The parameters used to fit the soil water retention curves of the soil materials w the van Genuchten-Mualem model are shown in Table 4. The saturated hydraulic ductivity (Ks) increased in most of the soil mixtures, with the addition of biochar. H ever, we observed a higher increase of Ks in the sandy loam, which coincides with la pore formation at less negative matric potentials ( Figure 5). Table 4. Soil hydraulic fitted parameters based on the van Genuchten-Mualem model (m = 1shown in Figure 5.   The Soil Water Storage Index (SWSCI) and Air Capacity Index (ACI) based on [41] are shown in Figure 6. The SWSCI index defines the water-filled pore space compared to the total porosity while the ACI index indicates the aeration capacity of the soil. It was found that the unamended silt loam had a SWCI index over the recommended optimal value of 0.66, and none of the biochar application rates enhanced the SWCI or ACI indexes. Contrarily, the unamended sandy loam was close to the optimum SWCI index and biochar fractions of 1% and 5% increased the ACI index, corresponding to the formation of macropores.

Variation of the Crop Productivity with the Biochar Fraction
Tomatoes were harvested for two months; a common practice among local producers of this specific variety. Figure 7 shows the average tomato yield by weight per plant, classified as class A, B, or C (bars), and the total fruit yield per plant is represented as a black continuous line. Additionally, the average number of fruits per plant by class, during the harvesting period, can be seen below the x-axis. This study was focused on classes A and B, due to their economic relevance to producers. In the silt loam treatments (SiL), class A tomato yields increased significantly with the 5% biochar addition (150% compared to the control), but class B yield decreased in average 25% at 2.5 and 5% biochar application rates. Overall, class B yield decreases were compensated by class A yield increments in the silt loam treatments. However, the total fruit yields of the SiL treatments were below the control. The sandy loam treatment with 5% biochar produced the highest fruit yield in class A, with less fruits per plant (352% compared to control). The 2.5% biochar fraction yielded the highest class; B tomatoes (27% compared to control). Overall, the sandy loam treatments produced a higher total fruit yield (3.8 kg/plant) than the silt loam treatments (3.0 kg/plant)

Variation of the Crop Productivity with the Biochar Fraction
Tomatoes were harvested for two months; a common practice among local producers of this specific variety. Figure 7 shows the average tomato yield by weight per plant, classified as class A, B, or C (bars), and the total fruit yield per plant is represented as a black continuous line. Additionally, the average number of fruits per plant by class, during the harvesting period, can be seen below the x-axis. This study was focused on classes A and B, due to their economic relevance to producers. In the silt loam treatments (SiL), class A tomato yields increased significantly with the 5% biochar addition (150% compared to the control), but class B yield decreased in average 25% at 2.5 and 5% biochar application rates. Overall, class B yield decreases were compensated by class A yield increments in the silt loam treatments. However, the total fruit yields of the SiL treatments were below the control. The sandy loam treatment with 5% biochar produced the highest fruit yield in class A, with less fruits per plant (352% compared to control). The 2.5% biochar fraction yielded the highest class; B tomatoes (27% compared to control). Overall, the sandy loam treatments produced a higher total fruit yield (3.8 kg/plant) than the silt loam treatments (3.0 kg/plant).

Discussion
This investigation revealed that the addition of biochar enhanced crop production by improving soil physical and chemical properties [47], with positive and negative effects depending on soil texture, which agrees with Sohi et al. [15]. Other studies, carried out in tomato crops, have also reported a relationship of crop yield to water and nutrient retention capabilities [5,6], soil aeration [8], and pH enhancement [48].
In our study, the addition of the bamboo biochar increased the pH less than one unit in the highest mixing rate (5%). Such small increments could be due to high levels of alkalis as Ca 2+ [49,50], carbonates, inorganic alkalis as silicates [51], and free hydroxyl (OH − ) ions due to the biochar ash content (>10%). The high concentration of Ca 2+ may have produced a cation imbalance, which is considered an indicator of low fertility [31,52].

Discussion
This investigation revealed that the addition of biochar enhanced crop production by improving soil physical and chemical properties [47], with positive and negative effects depending on soil texture, which agrees with Sohi et al. [15]. Other studies, carried out in tomato crops, have also reported a relationship of crop yield to water and nutrient retention capabilities [5,6], soil aeration [8], and pH enhancement [48].
In our study, the addition of the bamboo biochar increased the pH less than one unit in the highest mixing rate (5%). Such small increments could be due to high levels of alkalis as Ca 2+ [49,50], carbonates, inorganic alkalis as silicates [51], and free hydroxyl (OH − ) ions due to the biochar ash content (>10%). The high concentration of Ca 2+ may have produced a cation imbalance, which is considered an indicator of low fertility [31,52]. According to the results, we think that lime potentialized the effect of the liming capacity of the biochar. Lime CaO physically occupies the sites in the soil matrix normally occupied by ex-changeable aluminum while diminishing the solubility of aluminum (Al 3+ ) in the soil solution. Biochar alkalis, which according to Fidel et al. [51] include carbonates, inorganic alkalis, soluble alkalis, and structural alkalinity, may be more likely to neutralize the acidity of the soil solution. Our results did not reveal a significant CEC increase, hence there was no enhancement of nutrient retention capacity [15]. Similarly, the Ca/Mg/K ratios were not in their optimal ranges, therefore there was no improvement in the fertility of the soil [31]. Consequently, these properties did not provide any strong insight into our crop yield results. Nonetheless, these results contrast with other studies that reported a positive effect of biochar on CEC values of different soil amendments [53][54][55]. Low CEC values in our bamboo biochar were mainly due to the low SSA, where the low fractions of carboxylate groups did not contribute sufficiently to increase the negative surface charge on biochar particles [56]. Based on the low response of our crop yield to the chemical indicators, we could not set any strong driving mechanism to understand this yield output.
Contrarily, the effect of biochar on crop productivity was strongly associated with changes in soil water retention capacity, available water content (AWC), and air capacity (AC) conditions. Changes in soil water retention depend on the particle sizes of biochar and the soil texture, consistent with the work of Masiello et al. [57]. The larger particle size of the bamboo biochar compared to the soil materials led to a change in the grain size distribution in the mixtures. This, including the hydrophilicity of the biochar, strongly supports that the added biochar could have contributed to changes in AWC according to the particle size distribution of the different soil types. In the silt loam treatments, we found an opposite behavior between soil water retention and AWC. About 85% of the soil grains (smaller than the biochar particles) were able to fill the pores inside the biochar, forming new fine pores. This new pore formation in the intraparticle pore space increased soil water retention up to 10% and 8% of soil water retention in the SiL and SL treatments, respectively, though not necessarily AWC, especially in the SiL treatments. The new fine pores (<0.2 µm) were able to retain water inaccessible for plants at high matric potentials [57]. Additionally, the hydrophilic surfaces and the water-binding capacity of the biochar might have also played a role by increasing soil water retention at the permanent wilting point (PWP). Therefore, AWC, mainly in the silt loam treatments, did not significantly increase due to the formation of new fine pores able to hold soil water at high matric potentials. The increase of AWC in soils with high sand content agrees with previous studies [17,19,58] and differs with [24,26] that reported non-significant changes of biochar on water retention. Furthermore, it is contrary to available water content reports in silt loam soils by [59]. These contradictory results might be due to the different biochar particle size used in the studies [60].
Biochar kept high levels of the soil water storage index (SWSCI) in the silt loam treatments, and produced a close optimum equilibrium in the sandy loam mixtures, concerning the air capacity index (ACI). These results are consistent with Castellini et al. [61] who argue that coarse textures benefit more, in terms of their soil physical quality, than fine textures, from biochar additions, due to effects on its porosity and air capacity. Such a SWSCI/ACI relationship was directly correlated to fruit yield. Reynolds et al. [62] suggest SWSCI optimal values range between 0.6 and 0.7 and ACI optimal values surpass 0.14, in order to obtain optimal root-zone development. The formation of macropores (equivalent average pore diameter >50 µm) and capillary pores (pore diameter between 50 and 0.2 µm) in the sandy loam treatments produced optimal conditions for an ideal solute exchange, resulting in a higher total tomato yield, in contrast to the silt loam treatments, where AC conditions produced severe aeration deficits, non beneficial for plant growth [41,63]. These relations suggest that SWSCI/ACI imbalances may affect tomato production and is highly dependent on the soil type, biochar application rate, and its particle size [64].
Towards more efficient irrigation and management practices, we believe that differences in SWSCI/ACI indexes for the different soil textures may help promote adapted water management strategies based on soil physical changes provoked by biochar soil amendments. Adaptations in frequency and intensity of irrigation must be taken into consideration when adding biochar, to reduce water losses, from downward vertical processes, such as percolation, and to avoid excess water that might harm crop production.

Conclusions
Results demonstrate that the addition of a hydrophilic macroporous biochar to a sandy loam and a silt loam soil had dramatic effects on the soil's physical properties, leading to strong and contrasting impacts on the yield of tomatoes in a greenhouse experiment. The yields of A and B class tomatoes significantly increased in the sandy loam with the 2.5% biochar fraction (45% higher than control), whereas the contrary occurred in the silt loam treatments. For efficient water use and facing threats such as water scarcity and climate change, a water irrigation plan strategy for loam soils amended with hydrophilic macroporous biochars is recommended to control the water storage to air capacity ratio SWSCI/ACI.