10-Years Studies of the Soil Physical Condition after One-Time Biochar Application

: The ten-year experiment on the soil physical properties of biochar-amended Podzol was studied. Biochar was applied to the soil in the following rates: treatment BC10—10 Mg × ha − 1 , treatment BC20—20 Mg × ha − 1 , treatment BC30—30 Mg × ha − 1 and treatment BC0—Control (soil without the addition of biochar). Biochar was mixed the soil arable layer (0–20 cm). Soil samples were collected ten times, once a year—after harvest rye. They were taken from layers: 0–10 cm and 10–20 cm, in six replicates, using 100 cm 3 metal cylinders. The soil physical properties were determined: particle size distribution, particle density, bulk density, total porosity, air capacity and permeability (at − 15.5 kPa), water content at sampling, ﬁeld water capacity (at − 15.5 kPa), available and unavailable water content, and the ratio of ﬁeld water capacity and total porosity was calculated. It was found that biochar application causes changes in the soil physical condition. The soil density decreased, while the porosity, aeration and water retention increased; the ratio of ﬁeld water capacity and total porosity was favorable. These changes cannot be considered as permanent. Most of the analyzed properties showed a durability of no more than 3–4 years. We found that biochar incorporation into soil is a good method for environmental management of waste biomass.


Introduction
The data on global soil resources show the continuous increase in the area of degraded soils [1][2][3][4]. Therefore, an important element in management of soils is to prevent their degradation. Taking preventive actions is considered to be a better approach than the subsequent implementation of soil remediation measures [5][6][7][8][9]. Soil is an irreplaceable element of the environment and it is subjected to the pressure of various degrading factors. Excessively compacted soils occur commonly. The studies by some authors, e.g., Jones [10], McQueen and Shepherd [11], Drewry et al. [12], Reynolds et al. [13,14], and Pranagal [15], indicate that high soil compaction may limit the growth of crop plants.
Soil is a complex medium that undergoes dynamic changes both in time and on its surface [16][17][18][19][20]. The properties of soil can be modified by incorporating diverse mineral and/or organic materials into it [6,7,15,[21][22][23][24]. Waste materials originating from different industries, including agriculture, are used very frequently for this purpose. According to many authors, application of organic wastes to soil inhibits soil degradation, improves soil quality, and has a beneficial impact in the context of safe food. However, these authors also admit that most organic wastes decompose quickly and their next applications are required [8,9,21,22,25]. Management of waste by its application to soil allows two objectives to be achieved: (i) to dispose of waste, and (ii) to improve soil properties [6,8,9,15,21,22,[25][26][27]. very low oxygen level (1-2%) and at a maximum combustion temperature of 650 • C. The characteristics of the biochar used in the research are presented in the paper Pranagal et al. [8].
The field experiment was set up in the second 10 days of August 2010 in a field where winter rye (Secale cereale L., cv. "Dańkowskie Diament") was grown in monoculture. Monoculture cropping of winter rye continued for another 10 years (2010-2019). The maintenance of monoculture was intended to facilitate interpretation of results of determinations of soil physical properties. In this way, additional factors were avoided, which would have affected the soil physical condition, e.g., different agronomic practices and/or the effects of the root systems of crop plants. The experiment was set up in a randomized block design in three replicates. The area of each single plot was 18 m 2 . The spacing between plots fertilized with the different rates of biochar was 2 m. The biochar was applied according to the following experimental design: treatment BC10-10 Mg × ha −1 ; treatment BC20-20 Mg × ha −1 ; and treatment BC30-30 Mg × ha −1 . No biochar was applied to the soil in one of the treatments and this was the control treatment (BC0). The biochar was initially spread over the soil surface, and subsequently ploughing was done at a depth of about 20 cm in order to mix it with the soil. Then, the soil surface was leveled with a power harrow. In the third 10 days of September, winter rye was sown using a cultivating seed drill. Winter rye was sown every year in the third 10-day period of September. Mineral fertilizers were applied every year of the experiment at the following rates: 70 kg ha −1 N (ammonium nitrate), 26 kg ha −1 P (triple superphosphate), 66 kg ha −1 K (muriate of potash, KCl). Phosphorus and potassium fertilizers as well as 20 kg ha −1 N were applied before sowing. In spring, the remaining portion of the nitrogen (N) rate was applied before plant growth began (30 kg ha −1 ) and at the stem elongation stage (20 kg ha −1 ).
Undisturbed soil samples were collected once during every growing season-after harvest of rye in the first 10 days of August (sampling dates: I, II, III, . . . IX). In order to determine the initial condition of the soil, samples were taken before biochar application-sampling date 0 (2010). Soil samples were collected from the layers of 0-10 cm and 10-20 cm to 100 cm 3 metal cylinders in six replicates.

Analyses
Physical soil properties, such as particle size distribution (PSD), particle density (PD), bulk density (BD), total porosity (TP), air capacity (at −15.5 kPa) (FAC), air permeability (at −15.5 kPa) (FAP), water content at sampling (SM), field water capacity (FC), available water content (AWC), unavailable water content-permanent wilting point (UWC) and FC/TP ratio were studied. Soil water retention curve was determined with the use of a pressure plate apparatus (Soil-moisture Equipment Corp., Santa Barbara, California, USA). The level of field soil saturation with water was calculated for soil moisture level at a potential value of −15.5 kPa and permanent wilting point (UWC) of −1550.0 kPa. The content of organic carbon (TOC) was measured with the use of a Shimadzu TOC-VCSH analyzer with an SSM-5000A adapter for solid sample combustion.
The soil physical properties were determined according to the following procedures: -soil texture-particle size distributions (PSD)-with Casagrande method modified by Prószyński [53], -particle density (PD)-with the pycnometric method [54] (Mg × m −3 ), -bulk density (BD)-with the gravimetric method, from the ratio of the mass of soil dried at 105 • C to the initial soil volume of 100 cm 3 [55] (Mg × m −3 ), -total porosity (TP) was calculated from the results of particle density (PD) and bulk density (BD), air capacity at the potential of −15.5 kPa (FAC) was derived from the results of total porosity (TP) and field water capacity (FC) (−15.5 kPa), FAC = TP−FC [56] (m 3 × m −3 ), -air permeability at the potential of −15.5 kPa (FAP) was measured using an apparatus for the measurement of permeability of molding sand, LPiR-2e [7][8][9]. The measurements were conducted at vertical (upward) airflow through the soil sample. The pressure head in the measurement was not taken into account in the measurement results. The apparatus was produced by MULTISERW-Morek (Poland), -water content at sampling (SM) was calculated from the ratio of the mass of water contained in the soil during the sampling to the dry matter of soil dried at 105 • C [57] (kg × kg −1 ), -field water capacity (FC) was calculated from the ratio of the volume of water contained in the soil at the potential of −15.5 kPa to the soil volume [58,59] (m 3 × m −3 ), -available water content (AWC) was obtained from FC (−15.5 kPa) and unavailable water content (−1550.0 kPa)-permanent wilting point (UWC), AWC = FC−UWC [59,60] (m 3 × m −3 ), -unavailable water content (UWC) was calculated from the ratio of the volume of water contained in the soil at the potential of −1550.0 kPa to the soil volume [58,59] (m 3 × m −3 ).
Air-water relations of the soil were analyzed by determining the values of FC/TP ratio [14,60,61].

Statistical Analysis
The results were statistically evaluated with analysis of variance (ANOVA). All pairs of means between treatments were compared with the Tukey's test and the lowest significant difference test (LSD). The analysis (two-way ANOVA-biochar application rate (BC0, BC10, BC20, BC30) × soil layer and one-way ANOVA-soil from treatments BC0, BC10, BC20, BC30 and sampling dates 0-IX- Table 1) was performed for the results concerning the soil from the layer 0-10 cm, 10-20 cm ( Table 2) and 0-20 cm ( Table 3). In that manner mean values for each soil property under analysis were compared. The statistical evaluations (ANOVA-LSD) were conducted assuming the significance level of α = 0.05. Statistica 11 by Statsoft and ARSTAT by University of Life Sciences in Lublin were used to statistical analyses.

Soil Texture (PSD) and Density (PD and BD)
As a basic characteristic, soil texture (PSD) is considered to be very stable. Significant changes in soil particle size distribution occur very rarely. The studied soil was also characterized by stable PSD. The biochar incorporated into it at different rates (BC10, BC20, and BC30) did not cause any changes in soil texture. In the case of application of the higher biochar rates (BC20 and BC30), only minor changes, of about ±0.5%, in the content of the silt and clay fraction were recorded. Throughout the entire measurement period (2010-2019), the studied soil was classified as loamy sand (SL). As shown in the studies by Blott and Pye [53], Mikheeva [62] as well as Carter and Bentley [63], ±0.5% differences in PSD do not cause significant changes in soil functioning. Such changes can only marginally affect, e.g., soil compaction or/and air-water properties. The authors of some papers [5,7,8] have stressed that high rates of different soil-applied materials, both mineral and organic ones, result in a significant change in particle size distribution. The soil texture is then restructured towards the one that is determined by the type of the material applied to the soil [15,24,42].
Similarly to particle size distribution (PSD), particle density (PD) is a trait characterized by substantial stability [24,25,43,54]. In the present study, PD varied slightly from 2.56 Mg × m −3 (BC30; layer 0-10 cm; sampling date I) to 2.64 Mg × m −3 (BC0; layer 10-20 cm; sampling date 0) ( Figure 1). Particle density of soil without biochar ranged from 2.62 Mg × m −3 to 2.64 Mg × m −3 . PD in the 0-10 cm layer was higher than in the deeper layer of 10-20 cm almost in every case. Biochar application caused a distinct decrease in PD, which was proportional to the rate applied. The differences in PD results that were found immediately after biochar incorporation into the soil (sampling date I) between treatment BC0 and treatments BC10, BC20, and BC30 remained at a similar level until the end of the experiment (sampling date IX). It can therefore be concluded that biochar application caused persistent changes in PD ( Figure 1, Table 1). These changes were found to be persistent in the case of both soil layers Agronomy 2020, 10, 1589 5 of 22 (0-10 cm and 10-20 cm) and each biochar rate. According to the arithmetic means, the biochar rates of 20 Mg × ha −1 (BC20) and 30 Mg × ha −1 (BC30) resulted in a significant decrease in PD in comparison with the control soil-BC0 (ANOVA-LSD) (Tables 2 and 3). Similar changes in soil particle density (PD) have also been obtained by other authors, e.g., Pranagal [8], Meena et al. [23], and Githinji [64]. In their studies, these authors used very high biochar rates-25-100%, v/v. As a result, they found decrease to 40% PD. Nonetheless, these authors stressed that such a clear effect was short-lived (2-3 years). The relationships of PD with other soil properties, e.g., PSD, SOM, BD, FC, etc., are well known and have been frequently described by, among others, Jones [10], Reynolds et al. [24], Herath et al. [43], Blake and Hartge [54,55], Carter and Bentley [63], and Ball et al. [65].  (Tables 2 and 3). Similar changes in soil particle density (PD) have also been obtained by other authors, e.g., Pranagal [8], Meena et al. [23], and Githinji [64]. In their studies, these authors used very high biochar rates-25-100%, v/v. As a result, they found decrease to 40% PD. Nonetheless, these authors stressed that such a clear effect was short-lived (2-3 years). The relationships of PD with other soil properties, e.g., PSD, SOM, BD, FC, etc., are well known and have been frequently described by, among others, Jones [10], Reynolds et al. [24], Herath et al. [43], Blake and Hartge [54,55], Carter and Bentley [63], and Ball et al. [65].
; the FC/TP ratio. Each letter (a, b, c) means a significant difference (treatment × layer) according to Tukey's the lowest significant difference (LSD 0.05 ) and see Table 1.
Explanations as for Tables 1 and 2.
Bulk density (BD) is an important physical parameter considered to be one of the soil compaction measures [10][11][12][13][14][15]. The BD results obtained in this research were within a very wide range from 1.19 Mg × m −3 (BC30; layer 0-10 cm; sampling date I) to 1.73 Mg × m −3 (BC0; layer 10-20 cm; sampling date 0) ( Figure 2). As in the case of PD measurements, the topsoil layer (0-10 cm) was characterized by lower compaction than the 10-20 cm layer. The incorporated biochar definitely contributed to a decrease in BD. The higher the biochar rate, the more the BD decreased. The changes in bulk density related to both layers analyzed (0-10 cm and 10-20 cm). The differences between the soil in treatment BC0 and the soil in treatments BC10, BC20, and BC30, respectively, persisted only over three growing seasons (sampling dates I, II, and III). At sampling date IV, these differences disappeared completely (Table 1). In the next years (sampling dates V-IX), only differences that resulted rather from seasonal variation were observed ( Figure 2). An analysis of variance (ANOVA-LSD) showed that only the differences in mean BD values between treatment BC0 and treatment BC30 could be considered significant (Tables 2 and 3). The obtained BD results were similar to those found by other authors, e.g., Pranagal et al. [8,15], Drewry et al. [12], Reynolds et al. [13,14], Arshad and Martin [66], Logson and Karlen [67], and McQueen and Shepherd [11]. They indicated that the bulk density of soils with the particle size distribution of loamy sands most frequently ranges from 1.50 Mg × m −3 to 1.70 Mg × m −3 . These researchers also emphasized that favorable conditions for the growth of crop plants occur when BD ≤ 1.60 Mg × m −3 . Based on the classification by Paluszek [68], the investigated soil was most often classified in the following BD classes: "medium" and "high". Therefore, one could expect excessive soil compaction and disturbances in soil-plant-atmosphere relations according to the specialists involved in research on the soil physical condition [8][9][10][11][12][13][14][15]. Such a situation occurred before the establishment of the experiment (sampling date 0) and also at successive sampling dates from the fifth year of the experiment (sampling date IV). Such soil condition was predominantly found in the deeper layer (10-20 cm) ( Figure 2). The authors of similar studies [7][8][9][21][22][23][24]69] have also reported a deterioration in air-water conditions in the soil caused by a high BD. They underlined that the persistence of changes in the physical condition of the soil after its amendment with different additional materials was usually transient.
The TP results exhibited the same trends in changes as in the case of bulk density (BD). The addition of biochar, particularly at its highest rate (treatment BC30), had a significant effect on increasing total porosity (TP). The changes in TP applied to the entire arable soil layer (0-20 cm) to a similar degree. The differences due to biochar application persisted until the fifth year of the experiment (sampling date IV). At sampling date IV, the TP values equalized, while at subsequent dates (sampling dates V-IX) only small fluctuations in them were recorded (Table 1). According to the mean values, the TP of the soil amended with biochar at a rate of 30 Mg × ha −1 (treatment BC30) was significantly highest (ANOVA-LSD) (Tables 2 and 3). Numerous papers on soil porosity [9,[13][14][15]63,68] have demonstrated that an ideal situation occurs in the soil when TP ≥ 0.500 m 3 × m −3 . According to the above cited authors, such TP should provide proper soil-plant-atmosphere relations. The studied soil was usually characterized by TP that was lower than its recommended optimum value (Figure 3). Other authors [42][43][44][45][46][68][69][70][71] have claimed that TP ≥ 0.500 m 3 × m −3 is not sufficient to maintain optimal air-water conditions. They have also stressed that a favorable distribution of the network of soil pores and their patency are equally important as the total porosity (TP) of the soil.
The TP results exhibited the same trends in changes as in the case of bulk density (BD). The addition of biochar, particularly at its highest rate (treatment BC30), had a significant effect on increasing total porosity (TP). The changes in TP applied to the entire arable soil layer (0-20 cm) to a similar degree. The differences due to biochar application persisted until the fifth year of the experiment (sampling date IV). At sampling date IV, the TP values equalized, while at subsequent dates (sampling dates V-IX) only small fluctuations in them were recorded (Table 1). According to the mean values, the TP of the soil amended with biochar at a rate of 30 Mg × ha −1 (treatment BC30) was significantly highest (ANOVA-LSD) (Tables 2 and 3). Numerous papers on soil porosity [9,[13][14][15]63,68] have demonstrated that an ideal situation occurs in the soil when TP ≥ 0.500 m 3 × m −3 . According to the above cited authors, such TP should provide proper soil-plant-atmosphere relations. The studied soil was usually characterized by TP that was lower than its recommended optimum value ( Figure 3). Other authors [42][43][44][45][46][68][69][70][71] have claimed that TP ≥ 0.500 m 3 × m −3 is not sufficient to maintain optimal air-water conditions. They have also stressed that a favorable distribution of the network of soil pores and their patency are equally important as the total porosity (TP) of the soil.
Field air capacity (FAC) was calculated at a potential of −15.5 kPa. Such energy state of the soil (−15.5 kPa) was accepted as field water saturation of the soil. Then all pores with dimensions >20 µm-macropores are filled with air. Pores of this size perform an important role since they are responsible for gas exchange between the soil and the atmosphere. For this reason, they are often called aeration pores [56,[68][69][70][71]. Field air capacity (FAC) was calculated at a potential of −15.5 kPa. Such energy state of the soil (−15.5 kPa) was accepted as field water saturation of the soil. Then all pores with dimensions >20 μm-macropores are filled with air. Pores of this size perform an important role since they are responsible for gas exchange between the soil and the atmosphere. For this reason, they are often called aeration pores [56,[68][69][70][71].
In our study, the FAC value was in the range from 0.092 m 3 × m −3 (BCO; layer 10-20 cm; sampling date 0) to 0.245 m 3 × m −3 (BC30; layer 0-10 cm; sampling date II) ( Figure 4). The effect of biochar application was visible only in the two first years-sampling dates I and II. Subsequently, the above-mentioned differences continued to decrease (Table 1). In the control treatment (BC0), FAC was sometimes higher than in the biochar-amended soil (treatments BC10, BC20, and BC30). Such a situation occurred, for example, in: (i) layer 0-10 cm; sampling dates V and VI; and (ii) layer 10-20 cm; sampling dates III and V ( Figure 4). It can be presumed that air capacity was more affected by the seasonal dynamics than biochar application. The analysis of variance (ANOVA-LSD) of the means revealed that significantly the lowest FAC was found in the soil in treatment BC20 (Tables 2  and 3). Many authors [56,61,72,73] assume that FAC > 0.100 m 3 × m −3 in order to provide proper soil aeration. On the other hand, some authors think that the threshold value of FAC should be higher, e.g., >0.110 m 3 × m −3 [14,74,75]; >0.120 m 3 × m −3 [70], and >0.140 m 3 × m −3 [12,76,77]. The above indicated threshold values of FAC apply to full water saturation of the soil at a potential of −15.5 kPa (FC) [59,60]. Under natural conditions, such a state occurs only several times during a growing season and it is a quickly receding state [8,68,71,78]. In most measurements, the studied soil met the criterion of FAC >0.100 m 3   In our study, the FAC value was in the range from 0.092 m 3 × m −3 (BCO; layer 10-20 cm; sampling date 0) to 0.245 m 3 × m −3 (BC30; layer 0-10 cm; sampling date II) (Figure 4). The effect of biochar application was visible only in the two first years-sampling dates I and II. Subsequently, the above-mentioned differences continued to decrease (Table 1). In the control treatment (BC0), FAC was sometimes higher than in the biochar-amended soil (treatments BC10, BC20, and BC30). Such a situation occurred, for example, in: (i) layer 0-10 cm; sampling dates V and VI; and (ii) layer 10-20 cm; sampling dates III and V (Figure 4). It can be presumed that air capacity was more affected by the seasonal dynamics than biochar application. The analysis of variance (ANOVA-LSD) of the means revealed that significantly the lowest FAC was found in the soil in treatment BC20 (Tables 2 and 3). Many authors [56,61,72,73] assume that FAC > 0.100 m 3 × m −3 in order to provide proper soil aeration. On the other hand, some authors think that the threshold value of FAC should be higher, e.g., >0.110 m 3 × m −3 [14,74,75]; >0.120 m 3 × m −3 [70], and >0.140 m 3 × m −3 [12,76,77]. The above indicated threshold values of FAC apply to full water saturation of the soil at a potential of −15.5 kPa (FC) [59,60]. Under natural conditions, such a state occurs only several times during a growing season and it is a quickly receding state [8,68,71,78]. In most measurements, the studied soil met the criterion of FAC >0.100 m 3 × m −3 [56,61,72,73] as well as the minimum criterion proposed by Reynolds et al. [14,74] and Castellini et al. [75] (>0.110 m 3 × m −3 ), Walczak et al. [70] (>0.120 m 3 × m −3 ), Drewry et al. [12], Drewry [76], and Mueller et al. [77] (>0.140 m 3 × m −3 ). When analyzing FAC results, one needs to remember that in the case of sandy soils a greater problem can be an excess of air, not its lack. After adding the biochar to the soil, the porosity (TP) increased (Figure 3), and its substantial part was occupied by macropores (>20 µm). This situation was unfavorable because this caused excessive aeration of the soil. As a consequence, the soil could have become dry and soil organic matter mineralization could have accelerated [1,4,8,63,66]. It should be emphasized that in the case of the investigated soil the unfavorable changes only applied to sampling dates I and II. At the other sampling dates, i.e., dates III-IX, biochar was not found to have any negative impact on field air capacity (FAC).
Agronomy 2020, 10, x FOR PEER REVIEW 11 of 22 results, one needs to remember that in the case of sandy soils a greater problem can be an excess of air, not its lack. After adding the biochar to the soil, the porosity (TP) increased (Figure 3), and its substantial part was occupied by macropores (>20 μm). This situation was unfavorable because this caused excessive aeration of the soil. As a consequence, the soil could have become dry and soil organic matter mineralization could have accelerated [1,4,8,63,66]. It should be emphasized that in the case of the investigated soil the unfavorable changes only applied to sampling dates I and II. At the other sampling dates, i.e., dates III-IX, biochar was not found to have any negative impact on field air capacity (FAC). Field air permeability (FAP) at full water saturation (15.5 kPa) is an important parameter describing soil air properties. It characterizes the actual vertical gas flow in the soil and gas exchange between the soil and the atmosphere. FAP is a physical indicator of pore openness [79][80][81]. Measurements of FAP showed this parameter to vary significantly. Its value was within a wide range from 4.5 × 10 −8 × m 2 × Pa −1 × s −1 (BC0; layer 10-20 cm; sampling date 0) to 78.9 × 10 −8 × m 2 × Pa −1 × s −1 (BC30; layer 0-10 cm; sampling date V). As in the case of FAC, higher values were recorded in the topsoil layer (0-10 cm), while smaller ones in the lower layer of 10-20 cm. The biochar's effect increasing the air permeability remained visible for eight years of the experiment (sampling dates I-VIII). In the last (tenth) year of the experiment (sampling date IX), the observed differences in FAP disappeared completely ( Figure 5). Field air permeability (FAP) at full water saturation (15.5 kPa) is an important parameter describing soil air properties. It characterizes the actual vertical gas flow in the soil and gas exchange between the soil and the atmosphere. FAP is a physical indicator of pore openness [79][80][81]. Measurements of FAP showed this parameter to vary significantly. Its value was within a wide range from 4.5 × 10 −8 × m 2 × Pa −1 × s −1 (BC0; layer 10-20 cm; sampling date 0) to 78.9 × 10 −8 × m 2 × Pa −1 × s −1 (BC30; layer 0-10 cm; sampling date V). As in the case of FAC, higher values were recorded in the topsoil layer (0-10 cm), while smaller ones in the lower layer of 10-20 cm. The biochar's effect increasing the air permeability remained visible for eight years of the experiment (sampling dates I-VIII). In the last (tenth) year of the experiment (sampling date IX), the observed differences in FAP disappeared completely ( Figure 5). The greatest differences in FAP values between the control treatment BC0 and the biochar addition treatments (BC10, BC20 and BC30) were found in the 0-10 cm layer at sampling dates V and VI ( Figure 5). Given the above, a question arises why the largest range of FAP occurred not at sampling date I but at sampling date V. A probable reason for such FAP results could have been the interaction of a number of soil physical characteristics, e.g., PD, BD, FAC, SM, FC, etc. The comparison of the mean FAP values in the analysis of variance and of significant differences (ANOVA-LSD) revealed that such differences occurred both between the layers (0-10 cm and 10-20 cm) and treatments (BC0, BC10, BC20 and BC30) (Tables 1-3). FAP is a sensitive parameter for changes in the physical condition of the soil environment. It quickly responds to, e.g.: (i) modifications in the composition of individual soil components; (ii) fluctuations in soil moisture; and (iii) soil compaction. FAP also exhibits significant variations both within a single growing season and between seasons [82][83][84]. The studies by Paluszek [68] and Pranagal [71] demonstrated that proper gas flow in soil occurs when the vertical air permeability is at least 35.0 × 10 −8 × m 2 × Pa −1 × s −1 . This criterion was met by only 22/80 cases of means from the FAP measurements in the studied soil. Based on the division of soils into quality classes of permeability proposed by the above-mentioned authors [68,71], the studied soil was most frequently (53/80 cases) classified in the "medium" class (15-75 × 10 −8 × m 2 × Pa −1 × s −1 ), while only in 2/80 cases this was the "high" class (75-150 × 10 −8 × m 2 × Pa −1 × s −1 ). In the other cases, FAP was "low" (2.5-15.0 × 10 −8 × m 2 × Pa −1 × s −1 ). An increase in FAP resulting from biochar application proved to be quite persistent, unlike the parameters describing soil compaction (BD, TP and FAC), for which the changes showed lower persistence. Most probably, this was primarily attributable to the greater patency of soil pores than in the control treatment (BC0) and the presence of free interaggregate spaces and large intergrain pores. The authors of studies on soil air permeability [68,71,[79][80][81][82][83][84] have often indicated such a finding. The above cited papers also stressed some risks arising from a substantial increase in air permeability. Conditions for excessive aeration of the soil, and thus its drying, can then be created. The greatest differences in FAP values between the control treatment BC0 and the biochar addition treatments (BC10, BC20 and BC30) were found in the 0-10 cm layer at sampling dates V and VI ( Figure 5). Given the above, a question arises why the largest range of FAP occurred not at sampling date I but at sampling date V. A probable reason for such FAP results could have been the interaction of a number of soil physical characteristics, e.g., PD, BD, FAC, SM, FC, etc. The comparison of the mean FAP values in the analysis of variance and of significant differences (ANOVA-LSD) revealed that such differences occurred both between the layers (0-10 cm and 10-20 cm) and treatments (BC0, BC10, BC20 and BC30) (Tables 1-3). FAP is a sensitive parameter for changes in the physical condition of the soil environment. It quickly responds to, e.g.: (i) modifications in the composition of individual soil components; (ii) fluctuations in soil moisture; and (iii) soil compaction. FAP also exhibits significant variations both within a single growing season and between seasons [82][83][84]. The studies by Paluszek [68] and Pranagal [71] demonstrated that proper gas flow in soil occurs when the vertical air permeability is at least 35.0 × 10 −8 × m 2 × Pa −1 × s −1 . This criterion was met by only 22/80 cases of means from the FAP measurements in the studied soil. Based on the division of soils into quality classes of permeability proposed by the above-mentioned authors [68,71], the studied soil was most frequently (53/80 cases) classified in the "medium" class (15-75 × 10 −8 × m 2 × Pa −1 × s −1 ), while only in 2/80 cases this was the "high" class (75-150 × 10 −8 × m 2 × Pa −1 × s −1 ). In the other cases, FAP was "low" (2.5-15.0 × 10 −8 × m 2 × Pa −1 × s −1 ). An increase in FAP resulting from biochar application proved to be quite persistent, unlike the parameters describing soil compaction (BD, TP and FAC), for which the changes showed lower persistence. Most probably, this was primarily attributable to the greater patency of soil pores than in the control treatment (BC0) and the presence of free interaggregate spaces and large intergrain pores. The authors of studies on soil air permeability [68,71,[79][80][81][82][83][84] have often indicated such a finding. The above cited papers also stressed some risks arising from a substantial increase in air permeability. Conditions for excessive aeration of the soil, and thus its drying, can then be created.

Soil Water Properties (SM, FC, AWC and UWC) and FC/TP Ratio
The actual soil moisture (SM) is a very dynamic property. It predominantly depends on the amount of precipitation, evapotranspiration, and agronomic practices. During drought periods, it also tends to be dependent on: (i) soil texture; (ii) structure; and (iii) compaction [68][69][70][71]85]. The soil moisture (SM) during the experiment (2010-2019) was characterized by very high interseasonal variability. The SM values ranged from 0.059 kg kg −1 (BC0; layer 0-10 cm; sampling date IV) to 0.149 kg kg −1 (BC20; layer 10-20 cm; sampling date V). The effect of applied biochar on SM was particularly visible at sampling dates II, V, and VIII in both layers analyzed ( Figure 6).

Soil Water Properties (SM, FC, AWC and UWC) and FC/TP Ratio
The actual soil moisture (SM) is a very dynamic property. It predominantly depends on the amount of precipitation, evapotranspiration, and agronomic practices. During drought periods, it also tends to be dependent on: (i) soil texture; (ii) structure; and (iii) compaction [68][69][70][71]85]. The soil moisture (SM) during the experiment (2010-2019) was characterized by very high interseasonal variability. The SM values ranged from 0.059 kg kg −1 (BC0; layer 0-10 cm; sampling date IV) to 0.149 kg kg −1 (BC20; layer 10-20 cm; sampling date V). The effect of applied biochar on SM was particularly visible at sampling dates II, V, and VIII in both layers analyzed ( Figure 6). A state can therefore be drawn that the effect of biochar application was only periodically positive ( Table 1). The analysis of variance (ANOVA) and comparison of pairs of mean SM values (LSD) showed that significantly the least water was found in treatments BC0 and BC10 at sampling (Tables 1-3). Pandian et al. [86] informed that the effect of biochar incorporation into the soil was a 2.5% increase in its moisture. Recorded differences in SM can also be due to spatial variations that this characteristic is subject to. This applies to even small areas. It has been shown in papers by, among others, White [85], Kutílek [87], Leśny [88], Petrosyants et. al. [89] as well as Usowicz and Usowicz [90]. It should be underlined that this experiment was conducted in an area with dominant rainfall-based water management. Hence, SM was primarily dependent on the amount and distribution of precipitation and air temperature [85,88,89]. These factors cause recurrent soil wetting and drying cycles, while during the winter period freeze and thaw cycles. The cyclicity of these processes promotes the formation of soil aggregates [68,71,[91][92][93].
In the case of the determined field water capacity-FC (−15.5 kPa), the effect of adding the biomass processed through pyrolysis to the soil proved to be consistent with the soil compaction results (BD and TP) and the aeration properties (FAC and FAP).  A state can therefore be drawn that the effect of biochar application was only periodically positive ( Table 1). The analysis of variance (ANOVA) and comparison of pairs of mean SM values (LSD) showed that significantly the least water was found in treatments BC0 and BC10 at sampling (Tables 1-3). Pandian et al. [86] informed that the effect of biochar incorporation into the soil was a 2.5% increase in its moisture. Recorded differences in SM can also be due to spatial variations that this characteristic is subject to. This applies to even small areas. It has been shown in papers by, among others, White [85], Kutílek [87], Leśny [88], Petrosyants et al. [89] as well as Usowicz and Usowicz [90]. It should be underlined that this experiment was conducted in an area with dominant rainfall-based water management. Hence, SM was primarily dependent on the amount and distribution of precipitation and air temperature [85,88,89]. These factors cause recurrent soil wetting and drying cycles, while during the winter period freeze and thaw cycles. The cyclicity of these processes promotes the formation of soil aggregates [68,71,[91][92][93].
In the case of the determined field water capacity-FC (−15.5 kPa), the effect of adding the biomass processed through pyrolysis to the soil proved to be consistent with the soil compaction results (BD and TP) and the aeration properties (FAC and FAP). Small but visible variations in FC values were recorded from sampling date I to sampling date VII. However, they were usually greater in the deeper layer of 10-20 cm. This was also confirmed by the analysis of variance (ANOVA-LSD) which revealed significantly the highest FC values in the 10-20 cm layer. On the other hand, in the comparison of pairs of means in the treatments (BC0, BC10, BC20, and BC30), significantly the highest FC was found for the soil in treatment BC20 (Tables 2 and  3). The differences in FC results due to biochar application persisted for seven years (sampling dates I-VII). Subsequently, they began to disappear definitely from sampling date VII (Figure 7, Table 1). When assessing the field water capacity (FC) of the studied soil, it can be classified in the "medium" class according to the classification by Walczak et al. [70]. Nonetheless, when we take into account the recommendations given by Reynolds et al. [14], the obtained results were close to the determined optimum FC (0.300-0.350 m 3 × m −3 ) only in several situations. The obtained results confirmed the strong dependence of FC on soil compaction (Figures 2, 3 and 7) [63]. This dependence has also been described in experiments carried out in soils of other types [8,13,21,22,42,43].
In the context of ensuring stable yields, a very important parameter of the soil physical condition is available water content (AWC) [14,68,70,94,95]. Over the 10-year measurements (2010-2019), the AWC values were in a rather wide range from 0.181 m 3 × m −3 (BC30; layer 10-20 cm; sampling date VII) to 0.290 m 3 × m −3 (BC20; layer 0-10 cm; sampling date II). The nature of the changes was the same as in the case of FC (Figures 7 and 8, Table 1). Small but visible variations in FC values were recorded from sampling date I to sampling date VII. However, they were usually greater in the deeper layer of 10-20 cm. This was also confirmed by the analysis of variance (ANOVA-LSD) which revealed significantly the highest FC values in the 10-20 cm layer. On the other hand, in the comparison of pairs of means in the treatments (BC0, BC10, BC20, and BC30), significantly the highest FC was found for the soil in treatment BC20 (Tables 2 and 3). The differences in FC results due to biochar application persisted for seven years (sampling dates I-VII). Subsequently, they began to disappear definitely from sampling date VII (Figure 7, Table 1). When assessing the field water capacity (FC) of the studied soil, it can be classified in the "medium" class according to the classification by Walczak et al. [70]. Nonetheless, when we take into account the recommendations given by Reynolds et al. [14], the obtained results were close to the determined optimum FC (0.300-0.350 m 3 × m −3 ) only in several situations. The obtained results confirmed the strong dependence of FC on soil compaction (Figures 2, 3 and 7) [63]. This dependence has also been described in experiments carried out in soils of other types [8,13,21,22,42,43].
In the context of ensuring stable yields, a very important parameter of the soil physical condition is available water content (AWC) [14,68,70,94,95]. Over the 10-year measurements (2010-2019), the AWC values were in a rather wide range from 0.181 m 3 × m −3 (BC30; layer 10-20 cm; sampling date VII) to 0.290 m 3 × m −3 (BC20; layer 0-10 cm; sampling date II). The nature of the changes was the same as in the case of FC (Figures 7 and 8, Table 1).  The analysis of variance and comparison of pairs of AWC means (ANOVA-LSD) showed that the 10-20 cm soil layer in treatments BC20 and BC30 had significantly the best retention properties, while as regards the entire arable layer (0-20 cm), the soil in treatment BC20 retained most water (Tables 2 and 3). Visible differences in AWC for the analyzed treatments were also observed during measurements at sampling date VII (Figure 8). In the subsequent years of the experiment (sampling dates VIII and IX), the previously observed differences disappeared. Similar observations regarding an increase in stocks of water and its availability to plants as a direct result of biochar incorporation into the soil have also been reported by Abel et al. [42]-sandy soil, Pranagal et al. [8]-loamy sand soil, Githinji [64]-sandy loam soil, Herath et al. [43]-loamy soils (Alfisol and Andisol), Lu et al.
[96]-clayey soil (Vertisol) as well as by Omondi et al. in a meta-analysis [46]. Evaluation of the obtained AWC results allow us to state that the studied soil was characterized by good retention properties. According to Walczak et al. [70], in terms of AWC values the investigated soil could be classified in the "high" class (≥0.210 m 3 × m −3 -in 43/80 cases) and also in the "medium" class (0.120-0.210 m 3 × m −3 -37/80). In accordance with the proposals of Reynolds et al. [14] and Castellini et al. [75], on the other hand, the analyzed soil was assessed as "ideal" (≥0.200 m 3 × m −3 -in 49/80 cases), while in many cases also as "good" (0.150-0.200 m 3 × m −3 -31/80). The changes in AWC caused by biochar application were found over a period of seven years. They were most frequently beneficial, though small.
The unavailable water content (UWC) was typical of soils with the particle size distribution of loamy sand (LS) and with a low organic carbon content [8,63,[68][69][70][71]. The biochar incorporated into the soil resulted in a visible and persistent increase in UWC in micropores (<0.2 μm). During the measurement period, UWC ranged from 0.037 m 3 × m −3 (BC0; layer 0-10 cm; sampling date 0) to 0.049 m 3 × m −3 (BC30; layer 0-10 cm; sampling date I). The range of UWC results of similar magnitude (0.06-0.010 m 3 × m −3 ) persisted over the nine years of the experiment (sampling dates I-IX). The changes in the unavailable water content applied to both soil layers analyzed and were proportional to the amount of biochar applied (10, 20, and 30 Mg × ha −1 ) (Figure 9). This was confirmed by the analysis of variance (ANOVA-LSD) which revealed that significantly the highest UWC was in treatments BC20 and BC30 (Tables 1-3). The studied soil was characterized by a small The analysis of variance and comparison of pairs of AWC means (ANOVA-LSD) showed that the 10-20 cm soil layer in treatments BC20 and BC30 had significantly the best retention properties, while as regards the entire arable layer (0-20 cm), the soil in treatment BC20 retained most water (Tables 2 and 3). Visible differences in AWC for the analyzed treatments were also observed during measurements at sampling date VII (Figure 8). In the subsequent years of the experiment (sampling dates VIII and IX), the previously observed differences disappeared. Similar observations regarding an increase in stocks of water and its availability to plants as a direct result of biochar incorporation into the soil have also been reported by Abel et al.  [46]. Evaluation of the obtained AWC results allow us to state that the studied soil was characterized by good retention properties. According to Walczak et al. [70], in terms of AWC values the investigated soil could be classified in the "high" class (≥0.210 m 3 × m −3 -in 43/80 cases) and also in the "medium" class (0.120-0.210 m 3 × m −3 -37/80). In accordance with the proposals of Reynolds et al. [14] and Castellini et al. [75], on the other hand, the analyzed soil was assessed as "ideal" (≥0.200 m 3 × m −3 -in 49/80 cases), while in many cases also as "good" (0.150-0.200 m 3 × m −3 -31/80). The changes in AWC caused by biochar application were found over a period of seven years. They were most frequently beneficial, though small.
The unavailable water content (UWC) was typical of soils with the particle size distribution of loamy sand (LS) and with a low organic carbon content [8,63,[68][69][70][71]. The biochar incorporated into the soil resulted in a visible and persistent increase in UWC in micropores (<0.2 µm). During the measurement period, UWC ranged from 0.037 m 3 × m −3 (BC0; layer 0-10 cm; sampling date 0) to 0.049 m 3 × m −3 (BC30; layer 0-10 cm; sampling date I). The range of UWC results of similar magnitude (0.06-0.010 m 3 × m −3 ) persisted over the nine years of the experiment (sampling dates I-IX). The changes in the unavailable water content applied to both soil layers analyzed and were proportional to the amount of biochar applied (10, 20, and 30 Mg × ha −1 ) (Figure 9). This was confirmed by the analysis of variance (ANOVA-LSD) which revealed that significantly the highest UWC was in treatments BC20 and BC30 (Tables 1-3). The studied soil was characterized by a small percentage of UWC in FC. On average, this percentage did not exceed 20%. The comparison of the FC, AWC, and UWC results allowed finding that the biochar added to the soil had a positive effect. The increase in FC was primarily due to the increase in AWC arising from an increased content of mesopores (0.2-20.0 µm). It is worth stressing that the volume of micropores (<0.2 µm), in comparison with the volume of macropores (>20.0 µm) and mesopores (0.2-20.0 µm), showed much smaller variations. Similar observations have also been described by the authors of other experiments [5,8,22,24,25,46,78,88]. Nevertheless, the above-mentioned authors emphasized that the effect of biochar application was short-lived and most frequently limited to the topsoil layer (0-10 cm).
percentage of UWC in FC. On average, this percentage did not exceed 20%. The comparison of the FC, AWC, and UWC results allowed finding that the biochar added to the soil had a positive effect. The increase in FC was primarily due to the increase in AWC arising from an increased content of mesopores (0.2-20.0 μm). It is worth stressing that the volume of micropores (<0.2 μm), in comparison with the volume of macropores (>20.0 μm) and mesopores (0.2-20.0 μm), showed much smaller variations. Similar observations have also been described by the authors of other experiments [5,8,22,24,25,46,78,88]. Nevertheless, the above-mentioned authors emphasized that the effect of biochar application was short-lived and most frequently limited to the topsoil layer (0-10 cm). During the ten-year experiment (2010-2019), air-water conditions in the soil were also evaluated. To this end, the FC/TP ratio was determined, which underwent quite numerous changes. They applied to all treatments (BC0, BC10, BC20 and BC30). These changes resulted not only from the addition of biochar to the soil, but also from seasonal variation. The FC/TP value ranged from 0.53 (BC10; layer 0-10 cm; sampling date I) to 0.74 (BC20; layer 10-20 cm; sampling date III) ( Figure  10). In the control treatment (BC0), air-water conditions were most stable-the FC/TP ratio value ranged from 0.56 to 0.73. The calculated mean FC/TP values for all treatments were within the optimum range, i.e., 0.60-0.70 [14,60,61,63,74,75]. The analysis of variance (ANOVA-LSD) showed that significantly the lowest FC/TP values were found for the soil in treatment BC10 (Tables 1-3). It should however be stressed that periods when the FC/TP ratio diverged from its optimum value (0.60-0.70) were found with regard to the soil in all experimental treatments. Both the appearance of aeration impairment conditions (FC/TP > 0.70) and the condition of the soil signifying its excessive aeration, and thus a periodic water deficit (FC/TP < 0.60), could be found. During the ten-year experiment (2010-2019), air-water conditions in the soil were also evaluated. To this end, the FC/TP ratio was determined, which underwent quite numerous changes. They applied to all treatments (BC0, BC10, BC20 and BC30). These changes resulted not only from the addition of biochar to the soil, but also from seasonal variation. The FC/TP value ranged from 0.53 (BC10; layer 0-10 cm; sampling date I) to 0.74 (BC20; layer 10-20 cm; sampling date III) ( Figure 10). In the control treatment (BC0), air-water conditions were most stable-the FC/TP ratio value ranged from 0.56 to 0.73. The calculated mean FC/TP values for all treatments were within the optimum range, i.e., 0.60-0.70 [14,60,61,63,74,75]. The analysis of variance (ANOVA-LSD) showed that significantly the lowest FC/TP values were found for the soil in treatment BC10 (Tables 1-3). It should however be stressed that periods when the FC/TP ratio diverged from its optimum value (0.60-0.70) were found with regard to the soil in all experimental treatments. Both the appearance of aeration impairment conditions (FC/TP > 0.70) and the condition of the soil signifying its excessive aeration, and thus a periodic water deficit (FC/TP < 0.60), could be found. Agronomy 2020, 10, x FOR PEER REVIEW 17 of 22

Conclusions
The analysis of the obtained results allowed us to verify the hypothesis that one-time biochar amendment contributes to an improvement in soil physical properties and that resultant changes in soil properties are persistent. To this end, the range and persistence of changes in soil properties caused by agricultural use of biochar produced from waste winter wheat straw were evaluated in this study. A conclusion can be drawn that biochar application of BC10, BC20, and BC30 did improve soil physical properties, that is, the particle density (PD) and bulk density (BD) decreased, while the total porosity (TP) increased. The positive changes in PD, BD and TP also had an effect on increasing the following parameters: water content at sampling (SM), air capacity at −15.5 kPa (FAC), air permeability at −15.5 kPa (FAP), field water capacity at −15.5 kPa (FC), available water content (AWC), unavailable water content (UWC), and FC/TP ratio. It should however be remembered that an increase in soil air capacity and air permeability is not always positive. A large increase in the value of these traits can lead to unfavorable changes in air-water relations, i.e., to excessive aeration of the soil or even to its drying. Such a situation may apply in particular to poor quality soils with the sandy texture. The changes in the physical properties of soil described in this study, although visible, cannot be considered as permanent. Statistical analysis confirmed that only two soil properties has changed permanently: particle density (PD) and unavailable water content (UWC). Most of the analyzed properties showed a durability of no more than 3-4 years. We also found that biochar incorporation into soil is a good way for environmental management of waste biomass. The results presented by us and the observations collected during the experiment demonstrate that environmental management of biochar requires very great caution. When incorporating this type of waste into soil, one needs to remember not to disturb the balance in soil functions, i.e., the "container" and "filter" functions. To determine the timing for repeated agricultural, but safe use of waste, we propose to continue soil monitoring in subsequent years.

Conclusions
The analysis of the obtained results allowed us to verify the hypothesis that one-time biochar amendment contributes to an improvement in soil physical properties and that resultant changes in soil properties are persistent. To this end, the range and persistence of changes in soil properties caused by agricultural use of biochar produced from waste winter wheat straw were evaluated in this study. A conclusion can be drawn that biochar application of BC10, BC20, and BC30 did improve soil physical properties, that is, the particle density (PD) and bulk density (BD) decreased, while the total porosity (TP) increased. The positive changes in PD, BD and TP also had an effect on increasing the following parameters: water content at sampling (SM), air capacity at −15.5 kPa (FAC), air permeability at −15.5 kPa (FAP), field water capacity at −15.5 kPa (FC), available water content (AWC), unavailable water content (UWC), and FC/TP ratio. It should however be remembered that an increase in soil air capacity and air permeability is not always positive. A large increase in the value of these traits can lead to unfavorable changes in air-water relations, i.e., to excessive aeration of the soil or even to its drying. Such a situation may apply in particular to poor quality soils with the sandy texture. The changes in the physical properties of soil described in this study, although visible, cannot be considered as permanent. Statistical analysis confirmed that only two soil properties has changed permanently: particle density (PD) and unavailable water content (UWC). Most of the analyzed properties showed a durability of no more than 3-4 years. We also found that biochar incorporation into soil is a good way for environmental management of waste biomass. The results presented by us and the observations collected during the experiment demonstrate that environmental management of biochar requires very great caution. When incorporating this type of waste into soil, one needs to remember not to disturb the balance in soil functions, i.e., the "container" and "filter" functions. To determine the timing for repeated agricultural, but safe use of waste, we propose to continue soil monitoring in subsequent years.
Author Contributions: J.P., P.K., conceptualized the research; P.K., J.P., designed the research; J.P., P.K., performed the experiments; J.P. analyzed the data; J.P., P.K., wrote the paper. All authors have read and agreed to the published version of the manuscript.