Can application of biochar improve the soil water characteristics of silty loam soil?

This study investigates if a biochar application in 2014 and its reapplication in 2018 had affected the selected physical and hydrophysical soil properties of silty loam Haplic Luvisol at Dolná Malanta experimental site (Nitra, Slovakia) during the studied period of 2018–2020. Biochar was produced from cereal husks and paper fiber sludge by pyrolysis at temperature of 550 °C for 30 min and was applied to the soil in doses of 0, 10, and 20 t ha−1. The biochar treatments were combined with three levels of N fertilization (N0, N1, and N2). Biochar significantly decreased the bulk density up to − 12% and increased the soil porosity up to + 22% even 4–6 years after the biochar application. We also observed an increase in saturated hydraulic conductivity, which was reflected in better precipitation retention and a higher average soil water content by 12–15% depending on the fertilization level. Depending on the level of N fertilization, the interval of plant available water enlarged by 8–51% (10 t ha−1 of biochar), 18–21% (20 t ha−1 of biochar) after single biochar application, and by 18–34% (10 t ha−1 of biochar) and 19–31% (20 t ha−1 of biochar) after reapplication of biochar, which has positive implications regarding providing soil water to plants. The results indicate that biochar amendments have the potential to substantially improve the soil water characteristics, and therefore, the quality status of silty loam agricultural soils.


Introduction
Soil, with its water retention ability, plays an important role in influencing the natural water cycle. This important property can be best observed during extreme droughts under the conditions of climate change. Agriculture soil drought occurs when soil moisture drops below the limiting value of the wilting point when the plants are unable to absorb water sufficiently from the soil profile. The soil water constants such as wilting point delineate the specific soil moisture conditions that are defined at certain water potential values (Scherer et al. 2022) and are used also at monitoring the long-term changes in the course of soil moisture to define the soil water regime.
Currently, the improvement of water retention in the landscape and within the soil profile is an often-discussed topic, especially in relation to climate change mitigation (Basso et al. 2013). Some studies suggest that biochar can improve the quality of degraded soils (Hseu et al. 2014;Ndor et al. 2015;Karim et al. 2020;Xing et al. 2021) and the soil water regime (Hardie et al. 2013;Yu et al. 2013;Mukherjee and Lal 2013;Arthur et al. 2015;Castellini et al. 2015;Ma et al. 2016;Salinas et al. 2018).
Biochar is a carbon-rich product made from biomass by pyrolysis process that takes place at high temperatures Responsible editor: Lu Zhang ranging from 300 up to 900 °C under low oxygen or no oxygen conditions (Spokas 2010;Ok et al. 2016;Lehmann and Stephen 2015;Antunes et al. 2017;Inyang et al. 2016). The properties of biochar, to a large extent depend on the characteristics of the input raw materials and on the conditions of pyrolysis (Antunes et al. 2017). The higher pyrolysis temperatures (> 500 °C) result with higher pH, porosity, specific surface area, carbon content and stability of the final product. The sorption capacity of the biochar increases at lower pyrolysis temperatures (≤ 400 °C) (Wu et al. 2012;Zhao et al. 2013). Initially, the contribution of biochar application was discussed in terms of carbon sequestration and the biochar potential to reduce greenhouse gas emissions into the atmosphere (Lehmann et al. 2011;Rizhiya et al. 2019;Kotuš and Horák 2021). It is within the interest of the scientific community to discuss the benefits of biochar application also in terms of its impact on other agro-environmental parameters of soils such as pH, organic carbon content, organic matter content, soil structure, and soil aggregate stability) (Peng et al. 2011;Horák 2015;Šimanský et al. 2016;Ajayi and Horn 2016;Juriga and Šimanský 2019;Horák et al. 2020aHorák et al. , 2020b. Biochar has the potential to increase the soil water retention (Karhu et al. 2011) and the amount of water available to plants ( Baronti et al. 2014;Blanco-Canqui 2020;Medyńska-Juraszek et al. 2021). This means that the soils amended with biochar could potentially retain more water from precipitation, which could in turn lead to an increase of crop production in non-irrigated areas (Jeffery et al. 2011) and higher irrigation efficiency in irrigated areas. According to Wang et al. (2019) high dosage (≥ 10 t ha −1 ) of high pore volume biochar with bulky particle size (≥ 1 mm) can improve water retention of coarse-textured soil with limited water storage capacity and may improve soil's resilience during hydrological extremes. Biochar has the potential to influence the physical and hydrophysical characteristics of soil (Castellini et al. 2015;Makó et al. 2020) and thereby alter water balance in the ecosystems (Spokas 2010;Jeffery et al. 2015;Lehmann and Stephen 2015;Shackley et al. 2016;Antunes et al. 2017;Alghamdi 2018;Yu et al. 2019). However, the effect of biochar application is evaluated for a wide range of soil textures from sandy (Wiersma et al. 2020) to clay soils (Wong et al. 2022). It is believed that in clay soils, biochar can reduce bulk density and increase soil porosity. Bulk density is one of the most important physical properties of the soil, which influences the infiltration of precipitation (Ueckert et al. 1978). A decrease in soil bulk density (Laird et al. 2010;Githinji 2014;Are et al. 2017;Karim et al. 2020) increases soil porosity and soil aeration (Verheijen et al. 2010;Sun and Lu 2014;Walters and White 2018) and may have a positive effect on root and microbial respiration (Rutigliano et al. 2014). In addition, the effect of biochar on soil moisture (Novák et al. 2012), soil hydraulic conductivity (Lei and Zhang 2013;Esmaeelnejad et al. 2017), soil infiltration capacity (Abrol et al. 2016), field capacity (Jones et al. 2010), available soil water capacity (Brockhoff et al. 2010), and the relationship between soil water potential and soil moisture has been studied in the literature.
Most published works studying the effect of biochar on the physical and hydrophysical properties of soil are only short-term, especially the experiments that are carried out in laboratory conditions. There are only a limited number of studies that were conducted in field conditions over 5 years after biochar application in the soil. The impact of biochar reapplication on soil properties is an emerging topic, and to our knowledge, studies in this area are very rare (Nguyen et al. 2018). Therefore, this paper aims to analyze the impact of biochar application and reapplication in doses of 10 and 20 t ha −1 without or with nitrogen fertilizer on selected physical soil properties (bulk density, soil porosity) and hydrophysical soil properties (saturated soil hydraulic conductivity, soil moisture, and water availability for plants) from the 4th up to 6th year after biochar application at the experimental site in Dolná Malanta (Slovakia). We hypothesized that biochar reapplication will increase soil porosity and soil water content and that biochar application will affect the evaluated soil hydrophysical properties even 6 years after its application.

Study area
The experimental site in Dolná Malanta is located approximately 5 km north-east from the city of Nitra (Slovakia) (latitude 48°19′23.41″, longitude 18°09′0.7″) with an altitude of 175 m above sea level. The soil reference group at the site was characterized as Haplic Luvisol (IUSS Working Group WRB 2015). The soil on average contained 360.4 g kg −1 of sand, 488.3 g kg −1 of silt, and 151.3 g kg −1 of clay and is classified as a silty loam texture (Soil Portal 2019). The climate is characterized as warm, slightly dry, with an average annual air temperature of 9.8 °C and an average annual precipitation of 540 mm (according to the 30-year climate normal 1960-1991) (Šiška et al. 2005). The average annual air temperature during the studied years 2018, 2019, and 2020 was 11.4 °C, 10.9 °C, and 8.8 °C, respectively. The total precipitation in 2018, 2019, and 2020 was 537.4 mm, 625.4 mm, and 669.4 mm, respectively.
Prior establishment of the field experiment, the site was used for conventional agricultural crop production. Before the initialization of the experiment with biochar, soil samples were taken at ten random locations from a depth of 0-10 cm in order to obtain the initial soil properties. The average soil organic carbon content, soil pH (KCl), and dry bulk density was 10.2 g kg −1 , 5.58 and 1.39 g cm −3 , respectively (Kondrlova et al. 2018). The field experiment was established in spring 2014 to study the impact of biochar on the greenhouse gas emissions (Horák 2017;Horák et al. 2017), soil hydrophysical properties and soil quality (Horák 2015;Igaz et al. 2018;Šimanský et al. 2018;Tarnik 2019;Vitková and Šurda 2019;Toková et al. 2020;Horák et al. 2020a,b) and finally crop yields (Kondrlova et al. 2018;Aydin et al. 2020). Since the establishment of the experiment, common field crops were grown at the site; namely spring barley (Hordeum vulgare L.), corn (Zea mays L.) and peas (Pisum sativum L.) in 2018, 2019 and 2020, respectively.

Biochar application
Biochar used in this field experiment was made from paper fiber sludge mixed grain husks at a 1:1 ratio (according to their weight). The feedstock was processed by pyrolysis at 550 °C for 30 min in a Pyreg reactor (Pyreg GmbH, Gӧrhe, Germany) and provided by Sonnenerde company (Austria) as a final commercial product. The size of biochar particles was under 5 mm (Šimanský et al. 2019a). Table 1 provides further detailed information on the biochar parameters. Biochar was applied manually to the soil surface at doses of 0, 10, and 20 t ha −1 in the spring 2014 (Table 2) and it was incorporated into the soil with rakes. The study was conducted on plots of 4 × 6 m with three repetitions per treatment. In the spring 2018, the biochar-amended plots were divided into two separate 4 × 3 m plots. While one half of the original plot had the initial amount of biochar applied (2014), the other half received additional biochar (biochar reapplication in 2018, abbreviated "reap") at the same dose as when the field experiment was established in 2014 (Table 2).

N fertilizer application
An overview of the specific doses of N fertilizer in the studied years 2018, 2019, and 2020 is given in Table 2. Nitrogen industrial fertilizer (calcium ammonium nitrate LAD 27) was applied at the N1 and N2 application levels each year. The specific doses of N fertilizer at the N1 application level were calculated based on the requirements of grown crop each year using the balance method. An application dose at the N2 fertilization level was 50% higher than at the N1 level. An exception was made for spring barley in 2018, when the application dose at N2 fertilization level was doubled compared to the N1 fertilization level (Table 1).

Determination of the soil physical and hydrophysical properties
The basic soil physical properties (bulk density, particle density, porosity) and soil hydrophysical properties (hydraulic  conductivity, plant available water) were determined from undisturbed soil samples (100 cm 3 ). In total, 4 soil sampling events were conducted in the autumn 2018, spring 2019, autumn 2019, and spring 2020. Each time, 3 undisturbed soil samples were taken from a depth of 0-10 cm from each individual plot, which means 9 representative soil samples per each treatment (a grand total of 135 soil samples per sampling event). For more accurate statistical evaluation of data, one value with the biggest deviation from the data subset was excluded from further analysis (n = 8). Disturbed soil samples for determination of gravimetric water content were taken from a depth of 0-10 cm on biweekly intervals. A total of 13, 13, and 12 sampling events were performed in 2018, 2019, and 2020, respectively. At each sampling event, one representative soil sample was taken from each plot; therefore, there were 3 samples per treatment and total of 39, 39, and 36 samples per treatment in 2018, 2019, and 2020, respectively. The dry bulk density (BD) was calculated according to Eq. (1): The particle density (PD) of the dried undisturbed soil samples was determined indirectly using the air pycnometer method. This method operates on the principle of utilizing the vacuum created by the mercury column after inserting the undisturbed soil sample into an air pycnometer (Eijkelkamp Soil and Water 2020). This procedure determines the volume of the solid phase of the soil. PD was calculated according to Eq. (2) : Total porosity (P) was calculated according to the Eq. (3): The saturated hydraulic conductivity (K sat ) of the soil was determined in the laboratory using an empty extension cylinder placed on the top of saturated undisturbed soil sample and filled with tap water creating a variable hydraulic slope. Using this setting, non-stationary flow was created in a saturated soil sample due to the changing hydraulic slope Mohsenipour and Shaid 2016). K sat was calculated using Eq. (4): (1) BD = mass weight of dry soil (g) total volume of undisturbed soil sample (cm 3 ) (g cm −3 ) (2) PD = mass weight of dry soil (g) total volume of solid phase of the soil (cm 3 ) (g cm −3 ) (3) P = particle density g cm −3 − bulk density (g cm −3 ) where L is the height of the metal cylinder containing the soil sample (cm), t is time of water drop from the height H 2 to height H 1 (h), H 2 is the initial level of water in the extension cylinder (cm) and H 1 is the level of water in the extension cylinder after the water drop for time t (cm). The pressure plate apparatus was used to determine water content at specific pressure potentials. Plant available water (θ PAW ) was then calculated as the difference between field capacity (θ FC ) measured at a pressure potential of − 20 kPa and permanent wilting point (θ PWP ) measured at a pressure potential of − 1500 kPa (Eq. 5).
Mass soil water content (SWC) was calculated gravimetrically from disturbed soil samples dried at 105 °C (Eq. 6):

Statistical analysis
The effect of biochar application and subsequent reapplication on the physical and hydrophysical properties of the soil was evaluated using a one-way analysis of variance (one-way ANOVA). Statistically significant effects at p < 0.05 were determined by the least significant difference (LSD) test. All analyses were performed in Statgraphics Centurion XV.I software (Statpoint Technologies, Inc., Warrenton, VA, USA).

Effect of applied and reapplied biochar on soil bulk density
When looking at non-fertilized treatments (N0 level of fertilization), biochar did not have any significant effect on BD in the autumn 2018 (Table 3). A significant (p < 0.05) reduction of BD by 12 and 11% was observed in the treatments with a higher dose of biochar (B20 + N0 and B20reap + N0, respectively) in the spring 2019 (Table 4). In the autumn 2019, biochar application did not affect BD values, except for the treatment B10reap + N0 where BD increased significantly by 6% (Table 5). In the spring 2020, BD decreased weight of fresh soil (g) − weight of oven dry soil (g) weight of oven dry soil (g) . 100 (% mass) significantly (p < 0.05) in all treatments with biochar by 4 up to 6% (Table 6).
In fertilized treatments at the first level of fertilization (N1), only insignificant changes were observed in autumn 2018 and BD decreased slightly in the treatment B10reap + N1 compared to the control B0 + N0 (Table 3). In the spring 2019, there was a trend towards BD reduction in all treatments with biochar application, but a significant (p < 0.05) decrease by 8.5 and 9% was observed only in the treatments with a lower dose of biochar (B10 + N1 and B10reap + N1 respectively) ( Table 4). In the autumn 2019, BD generally increased slightly except a decrease in the B20 + N1 treatment. However, significant change (+ 5%) was observed only in the treatment B10reap + N1 compared to the control B0 + N0 (Table 5). In the spring 2020, BD decreased significantly (p < 0.05) in all treatments with biochar application in the range from 6 to 9% when compared to the control treatment B0 + N0 (Table 6).
A general decreasing trend in BD was observed after biochar application at the second level of fertilization (N2) during all three monitored years. In the autumn 2018, a significant decrease of BD was observed in the treatments B10reap + N2, B20 + N2, and B20reap + N2 compared to the treatment without biochar with nitrogen fertilizer (B0 + N2) (Table 3). However, BD in this control treatment (B0 + N2) was slightly higher (but insignificant) compared to the unfertilized B0 + N0 control (Table 3). In the spring 2019, there was a trend towards BD reduction in all treatments with biochar, but a significant decrease (p < 0.05) was only observed in the treatments B10reap + N2, B20 + N2, and B20reap + N2 (Table 4). In the autumn 2019, lower BD was again observed in the treatments B10reap + N2, B20 + N2, and B20reap + N2, but this reduction was significant (p < 0.05) only in comparison to fertilized control treatment (B0 + N2) ( Table 5).
In the case of B0 + N2 treatment, BD increased slightly compared to the B0 + N0 control (Table 5). A significant reduction (p < 0.05) of BD ranging from 4 to 10% was observed in all treatments with biochar addition in spring 2020 compared to the control B0 + N0. Soil bulk density is an important indicator of soil compaction because it is directly affected by soil structure (Hillel 1998). A decrease in BD after biochar application could be caused by the lower BD of biochar (Table 1) compared to BD of silty loam soil at the experimental site. Biochar itself is a porous material (Adekiya et al. 2020), and when added to soil, it decreases bulk density and increases soil porosity (Nyambo et al. 2018). The change in porosity in soils treated with biochar is mainly due to the formation of macropores and rearrangement of soil particles (Hseu et al. 2014), aeration of the soil and reduced mechanical impedance of root growth in soils (Blanco-Canqui 2017). Biochar has the ability to form soil aggregates in combination with the soil particles. The improvement of structural properties of soil in this field experiment has already been observed in the past and it was published in other works (Šimanský 2016; Šimanský et al. 2019b). Some authors argue that the basic cations present in biochar can bond with clay and organic particles in soil via cation bridges to form soil aggregates (Bronick and Lal 2005). Gradual oxidation of biochar in soil also increases the number of functionality groups on its surface (Obia et al. 2016), which can absorb soil particles and minerals, which in turn contributes to the formation of soil aggregates and to reduction of BD. In our work, however, we also observed an increase in BD in some treatments (autumn sampling) amended only with nitrogen fertilizer (B0 + N1 and B0 + N2). This phenomenon occurred because mineral nitrogen in the soil can accelerate the mineralization of organic matter in the soil, which can lead to an increase in BD (Šimanský et al. 2016). In contrast, biochar combined with nitrogen fertilizer generally resulted with a reduction of BD. This is because nitrogen has a beneficial effect on the incorporation of biochar into the soil. In our study, the use of a higher dosage of biochar (20 t ha −1 ) in combination with a higher level of fertilization (N2) seemed to be more effective.

Effect of applied and reapplied biochar on soil porosity
In general, we observed a trend in a decrease of BD and an increase of PD during all three monitored years in the treatments with biochar that led to an increase of P (Tables 3, 4, 5, and 6).
At the first level of fertilization (N1), P increased insignificantly in the autumn 2018 in the treatments B10reap + N1, B20 + N1, and B20reap + N2 (Table 3). In the spring 2019, porosity values increased significantly (p < 0.05) in all treatments with biochar in the range from 12 to 17% compared to the control treatment at the fertilization level B0 + N1. At the same time, however, P in B0 + N1 was lower by 8.6% than in the non-amended control treatment B0 + N0 (Table 4). Autumn samples in 2019 showed an insignificant slightly increasing trend in P in the treatments B10reap + N1, B20 + N1, and B20reap + N1 compared to the control B0 + N1 (Table 5). Spring samples in 2020 showed a significant increase (p < 0.05) of P in all biochar treatments ranging from 12 to 18% compared to B0 + N0, except treatment B20reap + N1 (Table 6).
At the second level of fertilization (N2), we repeatedly observed a trend towards an increase of P in all monitored years. In the autumn 2018, P increased significantly (p < 0.05) by 16 up to 25% in the treatments B10reap + N2, B20 + N2, and B20reap + N2 compared to the fertilized control B0 + N2 (Table 3). At the same time, P in B0 + N2 control treatment was lower by 8% than in the non-amended control treatment B0 + N0 (Table 3). In the spring 2019, biochar amendment generally increased P in all treatments. The effect was insignificant when compared to B0 + N0 treatment; however, in comparison to B0 + N2, an increase of P was significant ranging from 10 to 12% (Table 4). Application of biochar showed positive effect also in the autumn 2019, as P increased in all treatments. In the case of B10reap + N2, B20 + N2, and B20reap + N2 treatments, the increase by 10, 6.7, and 10% was significant (p < 0.05) compared to the control B0 + N2 (Table 5). In the spring 2020, a significant increase (p < 0.05) of soil porosity ranging from 13 to 20% was observed in all biochar treatments compared to the treatment B0 + N0 (Table 6).
As mentioned above, biochar alone is a highly porous material with a low bulk density (Jones et al. 2010;Duarte et al. 2019;Adekiya et al. 2020), and therefore, its application to soil can increase soil porosity. In our work, we generally observed a decreasing trend in BD, especially in treatments with a higher dose of biochar and also increasing trend in P after application of biochar into the soil. Our observations are consistent with other studies (Chan et al. 2008;Githinji 2014;Mukherjee and Lal 2013;Nelissen et al. 2015). Since P was calculated from BD and PD, its increase is mainly related to a decrease of BD and increase of PD. It is assumed that P increase may also be the result of improvements in soil structure after biochar application (Šimanský 2016).

Effect of applied and reapplied biochar on saturated hydraulic conductivity
Application of biochar alone or in combination with fertilizer slightly increased the values of K sat almost in all treatments during the whole studied period, except for the spring 2020.
In the case of the unfertilized treatments (N0), an increase of K sat in B10reap + N0, B20 + N0, and B20reap + N0 treatments ranged from 5 to 560% in the autumn 2018. In the case of B20reap + N0, the difference was significant (p < 0.05) compared to B0 + N0 (Table 3). In the spring 2019, an insignificant trend of increased K sat was observed in all treatments with biochar in the range from 6 to 464% compared to the control (Table 4). In the autumn 2019, K sat increased in treatments B10 + N0, B20 + N0, and B20reap + N0 in the range Table 6 Effect of biochar application and reapplication on soil properties in spring 2020 (means ± standard deviations) (n = 8) Different letters (a, b, c, d) indicate that treatment means are significantly different at p < 0.05 according to least significant difference test (oneway ANOVA). For explanation of treatment acronym see from 127 to 284%. In the case of B20reap + N0, the difference was significant (p < 0.05) ( Table 5). In the spring 2020, K sat increased slightly only in the treatment B10 + N0 (Table 6); however, all observed changes were insignificant (p > 0.05). At the first level of fertilization (N1), K sat increased insignificantly in all treatments with biochar in the range from 55 to 145% in the autumn 2018 (Table 3) and in the range from 20 up to 350% in the spring 2019. Only in the case of B20reap + N1 treatment in spring 2019, the increase was significant (p < 0.05) ( Table 4). In the autumn 2019, K sat increased slightly only in the treatments B10 + N1 and B10reap + N1 compared to the control B0 + N0 (Table 5); however, all observed values for biochar treatments were significantly lower when compared to fertilized control B0 + N1. In the spring 2020, K sat increased slightly only in the treatments B10reap + N1 and B20 + N1, however, all observed changes were insignificant (Table 6).
At the second level of fertilization (N2), we observed a slight increase of K sat in the autumn 2018 in all treatments with biochar ranging from 23 to 333%. However, the difference was significant (p < 0.05) only in case of the treatment B10 + N2 compared to control treatment (Table 3). In the spring 2019, K sat generally increased in the range from 100 to 299% but a significant increase (p < 0.05) was observed only in the treatment B20 + N2 compared to B0 + N0 (Table 4). In the autumn 2019, K sat increased insignificantly in all treatments with biochar in the range from 63 to 369% compared to control B0 + N0 (Table 5). However, all observed values for biochar treatments were significantly lower when compared to fertilized control B0 + N2. In the spring 2020, the application of biochar did not generally affect K sat . The only exception was B10 + N2 treatment where we observed a significant increase of 121% compared to the control B0 + N0.
Increasing the velocity at which water enters the soil is important for capturing rainfall, soil water retention, and overall soil management (Blaco-Canqui 2017). From this point of view, K sat is an important soil property in solving the problems with soil water regime. In our study, we observed a trend of a K sat increase in loamy soil in almost all cases with applied biochar, either in combination with N fertilizer or without N fertilizer. Some studies suggest that biochar can increase K sat in clay soils and sandy-loam soils (Dan et al. 2015) by increasing the number of connected macropores. Another study suggests that biochar could effectively suppress water loss in coarse sandy soils (Uzoma et al. 2011) and thus reduce K sat due to the absorption capacity of biochar. The changes in K sat are associated with changes in soil porosity, aggregation, and water retention capacity (Nelissen et al. 2015). The size of the biochar particles itself and the particle size of the soil to which the biochar was applied also play an important role (Hardie et al. 2013;Jeffery et al. 2015;Lehmann and Stephen 2015;Esmaeelnejad et al. 2017). The fact that the average K sat values varied in such a wide range compared to the control is not surprising, since this value has a high variability. The standard deviation between the observed average values was rather large (Tables 3, 4, 5, and 6). K sat depends on soil structure, which can vary in both space and time. According to Chirico et al. (2007), the relationship between K sat and other soil properties, such as soil texture, soil bulk density, and organic matter, is not strong enough to permit accurate estimations of K sat . Any change in the surface soil layer caused by natural factors, such as rainfall, growth and decay of plant roots, shrinking and swelling, or human interventions, such as cultivation treatments and wheel-traffic compaction can lead to changes in soil porosity. These changes can seriously affect K sat and consequently soil water storage (Kargas et al. 2021). According to Alletto and Coquet (2009), K sat of agricultural soil changes following a seasonal pattern. A decreasing tendency is expected from tillage to crop harvest due to soil compaction and progressive pore clogging.
Regardless of the determination method of K sat in the field or laboratory, the measured value is representative for the specific place and time at which the measurement was carried out (Kargas et al. 2021). The measurements of K sat can be affected by several factors. The first factor is the quality of the undisturbed soil sample and the content of the preferred water flow pathways in the soil sample. In one treatment of the field experiment, 9 soil samples were taken (8 used for further analysis), each of which could have a different composition (regarding the amount, size, and space distribution of pores). Another factor influencing the measurement could be the number of repetitions of the measurement. In our case, we performed two measurements per sample.

Effect of applied and reapplied biochar on plant available water
In the non-fertilized treatments (N0), we observed a trend towards an increase of the θ PAW in the autumn 2018 in treatments with reapplied biochar compared to B0 + N0. In the case of treatment B20reap + N0, the difference was significant (p < 0.05) (+ 16%) compared to control B0 + N0. In both treatments with biochar reapplication the reason was a slight increase of θ FC and in the case of treatment B10reap + N0 also a decrease of the θ PWP (Table 3). In the spring 2019, θ PAW increased significantly (p < 0.05) by 76, 37, and 65% in the treatments B10reap + N0, B20 + N0 and B20reap + N0 compared to the control treatment B0 + N0 (Table 4). This change was caused by a significant increase of θ FC (+ 5.6% in the treatment B10reap + N0) and a significant decrease of θ PWP (− 6.4 and − 8% in the treatments B20 + N0 and B20reap + N0, respectively). Although we did not observe any significant effect of biochar amendment on θ FC and θ PWP in the autumn 2019 compared to B0 + N0, θ PAW increased significantly (p < 0.05) in the treatments B10 + N0, B10reap + N0, and B20reap + N0 in the range from 16.7 up to 27% compared to B0 + N0 (Table 5). In the spring 2020, a general increase of θ PAW was observed in all treatments with biochar compared to B0 + N0. Significant (p < 0.05) changes of 26% and 29% were observed in the treatments B10 + N0 and B10reap + N0 compared to control B0 + N0 (Table 6).
At the first level of fertilization (N1), θ PAW increased insignificantly in the treatments B10 + N1, B10reap + N1, and B20reap + N1 in the autumn 2018 (Table 3). In the spring 2019, an increase of θ PAW was observed in all treatments with biochar compared to B0 + N0. Moreover, in most of the treatments the changes were significant (p < 0.05) (except for the treatment B20 + N1) and ranged from 34 to 62% compared to control treatment B0 + N0 (Table 4). At the same time, we observed an insignificant increase of θ FC and general decrease of θ PWP in all treatments with biochar ( Table 4). The values of θ PAW in the autumn 2019 increased only insignificantly in all treatments with biochar ( Table 5). The same trend was observed also in the spring 2020, when θ PAW increased after biochar application. Moreover, in the treatment B10 + N1 the difference was significant (p < 0.05) compared to control B0 + N0 (Table 6). At the same time, we did not observe any effect of biochar on θ FC , while θ PWP decreased significantly in the treatments B10 + N1 and B10reap + N1 by 4.7 and 5.4%, respectively (Table 6).
At the second level of fertilization (N2), θ PAW increased in the autumn 2018 in the treatments B10 + N2, B10reap + N2, and B20reap + N2 in the range from 5 up to 17% (Table 3). In the spring 2019, an increase of θ PAW in the range from 18 to 94% was observed in all treatments with biochar; however, the values were significantly different (p < 0.05) compared to the control B0 + N2 only in the treatments with reapplied biochar B10reap + N2 and B20reap + N2 (Table 4). At the same time, θ FC increased and θ PWP decreased significantly only in the treatment B20reap + N2 (by + 6.6%) and B10 + N2 (by − 7.5%) respectively compared to control treatment B0 + N0, while the effect of biochar in other treatments was insignificant (Table 4). In the autumn 2019, θ PAW increased in all treatments with biochar, but a significant increase of θ PAW (p < 0.05) was observed only in the treatment B20 + N2 (+ 19%) (Table 5). In the spring 2020, θ PAW increased slightly in all treatments with biochar but the difference (+ 21%) was significant (p < 0.05) only in the treatment B20 + N2 when compared to the control B0 + N0 (Table 6).
Plant available water in the soil is determined according to the range between the values of the soil water constants θ FC and θ PWP . The larger this interval is the more water can plants draw from the soil through their roots. In our work, we observed a general trend of an increase in θ FC and a decrease of θ PWP . This fact indicates that the interval between these two water constants has enlarged, which in turn has led to an increase of θ PAW . In summary, it can be concluded that in our observations, biochar increased the interval of available water for plants by up to tens of percent. Depending on the amount of applied biochar and the amount of nitrogen fertilizer, the average increase in the interval of plant available water ranged from + 10 up to + 34%. A more pronounced increase occurred in the treatments with biochar reapplication at 10 t ha −1 (18-34%) and 20 t ha −1 (19-31%). Another positive finding is that biochar application in 2014 still had a favorable effect on θ PAW in 2020. An average increase in θ PAW depending on the amount of applied nitrogen fertilizer ranged from 8 up to 51% (at 10 t ha −1 of biochar) and from 18 up to 21% (at 20 t ha −1 of biochar) when compared to the control.

Effect of applied and reapplied biochar on soil water content
Soil water content (SWC) was determined at 13 sampling events in 2018 and 13 and 12 sampling events in 2019 and 2020, respectively. Table 7 presents the mean observed values of SWC during the whole studied period. Considering a trend that SWC took in the unfertilized treatments (N0) and treatments at the first level of fertilization (N1), biochar application at 10 t ha −1 less often increased SWC than at the application rate of 20 t ha −1 . The most significant improvement was observed in the reapplication of biochar in a dose of 20 t ha −1 at all fertilization levels (N0, N1, N2), while the average SWC for the entire "vegetation" period was higher by 12-15% in 2018, 3-6% in 2019, and 2-7% in 2020 when compared to the control. It seems that the positive impact of biochar on SWC decreases over time due to agrotechnical operations Table 8 Percentual change in soil water content (%) and number of events with higher and lower soil water content during 2018 when compared to the control treatment The percentual difference in soil water content between the specific treatment and the control is accompanied by the color gradient from dark red (the lowest value in comparison to control) to dark blue (the highest value in comparison to control). + and − stands for the number of events with higher, respectively lower soil water content in comparison to the control treatment. For explanation of treatment acronym see Table 2   Table 9 Percentual change in soil water content (%) and number of events with higher and lower soil water content during 2019 when compared to the control treatment The percentual difference in soil water content between the specific treatment and the control is accompanied by the color gradient from dark red (the lowest value in comparison to control) to dark blue (the highest value in comparison to control). + and − stands for the number of events with higher, respectively lower soil water content in comparison to the control treatment. For explanation of treatment acronym see and admixing of biochar into deeper layers. However, we assume that the positive effect of biochar on SWC increases in the deeper layers of the soil profile where biochar gradually reaches over time.
The positive effect of a higher dose of reapplied biochar on SWC is confirmed when analyzing individual terms of soil moisture measurement. Overall biochar amendment did not significantly affect SWC at any fertilization level, except for reapplication of biochar in a dose of 20 t ha −1 for all fertilization levels (N0, N1, N2).
Overall, SWC was not significantly affected by biochar application at any fertilization level, except for reapplication of biochar at a dose of 20 t ha −1 at all fertilization levels (N0, N1, N2). Biochar reapplication at that rate (B20reap) had the highest potential to generally increase SWC as it was observed during the individual sampling events (Tables 8, 9, and 10). In the treatment B20reap + N0, a trend to increase SWC was observed at 10, 10, and 7 sampling events in 2018, 2019, and 2020, respectively. At the first level of fertilization the same application dose of biochar (B20reap + N1) tended to increase SWC at 12, 10, and 9 sampling events in 2018, 2019, and 2020, respectively. In the treatment B20reap + N2, a trend to increase SWC was observed at 10, 10, and 10 sampling events in 2018, 2019, and 2020, respectively.
One of the main advantages of applying biochar into soil is its ability to retain water by its natural internal porosity (Yu et al. 2013;Rasa et al. 2018) or its ability to aggregate soil (Herath et al. 2013;Laird et al. 2010), which is indirectly related to soil water content (Rasa et al. 2018). Our findings on the positive effect of biochar on soil water content are consistent with the work of Leelamanie (2014) and Agegnehu et al. (2015).

Conclusions
Based on the results observed in our field experiment with biochar at the experimental site in Dolná Malanta (Nitra, Slovakia) we can conclude that the application of biochar in 2014 and its subsequent reapplication in 2018 resulted in positive trends in changes of studied soil hydrophysical characteristics and storage of available water for cultivated plants. The results showed that biochar can reduce bulk density in silty loam Haplic Luvisols and possibly in similar soil types with loamy texture.
The most significant trend towards a reduction of bulk density was observed in the treatments with biochar application and reapplication at dose of 20 t ha −1 in combination with all levels of fertilization. At the same time, depending on the amount, the application of biochar increased soil porosity in the range from + 10 up to + 20%. Subsequent reapplication of biochar after 4 years supported this positive trend of increasing porosity in loamy soil. Increasing soil porosity by the biochar application results in an improvement of the soil structure, the quality of the soil pores, and their interconnection. This has a positive effect on the observed increase in saturated hydraulic conductivity of the soil. The increase in hydraulic conductivity was reflected in better retention of precipitation, increased soil water content and in overall improved hydrological soil balance. The highest average increase in soil water content was observed in the treatments with reapplied biochar at a dose of 20 t ha −1 (12-15%) in combination with all levels of fertilization (N0, N1, and N2). We assume, that the repeated biochar application influenced soil pore size redistribution towards an Table 10 Percentual change in soil water content (%) and number of events with higher and lower soil water content during 2020 when compared to the control treatment The percentual difference in soil water content between the specific treatment and the control is accompanied by the color gradient from dark red (the lowest value in comparison to control) to dark blue (the highest value in comparison to control). + and − stands for the number of events with higher, respectively lower soil water content in comparison to the control treatment. For explanation of treatment acronym see increase of the proportion of capillary pores in the soil. Biochar directly led to an enlargement of the interval of plant available water by decreasing the value of permanent wilting point and increasing the field capacity and therefore had a favorable effect on providing sufficient amount of water to plants. Depending on the level of N fertilization, the interval of plant available water increased by 8-51% (biochar application dose of 10 t ha −1 ), 18-21% (biochar application dose of 20 t ha −1 ), 18-34% (biochar reapplication dose of 10 t ha −1 ) and 19-31% (biochar reapplication dose of 20 t ha −1 ).