Seasonal Dynamics of Soil Moisture in an Integrated-Crop-Livestock-Forestry System in Central-West Brazil

: Integrated-crop-livestock-forestry (ICLF) systems are currently promoted as a measure for sustainable intensiﬁcation of agricultural production. However, due to complex interactions among ICLF components, we are still lacking evidence about the system’s resilience regarding water availability, especially for regions characterized by pronounced wet and dry seasons and frequent droughts. For a mature ICLF system in the Cerrado biome of central-west Brazil comprising rows of eucalyptus trees ( Eucalyptus grandis x Eucalyptus urophylla , H13 clone) at a spacing of 22 m in combination with Brachiaria brizantha cv. BRS Piat ã pasture we continuously measured soil moisture (SM) until 1 m depth and supported this data with measurements of photosynthetically active radiation (PAR) and aboveground green grass biomass (AGBM) across transects between the tree rows for almost two years. Across the seasons a distinct gradient was observed with SM being lower close to the tree rows than in the space between them. During winter SM decreased to critical values near the tree lines in the topsoil. During spring and summer, incident PAR was 72% and 86% lower close to the trees than at the center point. For autumn and winter PAR was more evenly distributed between the tree rows due to inclination with notably up to four times more radiation input near the tree lines compared to spring and summer. AGBM showed a clear distribution with maximum values in the center and about half of the biomass close to the tree rows. Our data suggest that, restrictions in AGBM accumulation shifted among seasons between water limitations in winter and light limitations during summer. Interestingly, SM changes during wetting and drying events were most pronounced in subsoils near the tree rows, while the topsoil showed much less ﬂuctuations. The subsoil in central position showed the lowest SM dynamics in response to drought maintaining a relative high and constant SM content, therefore functioning as important water reservoirs likely improving the resilience of the system to drought stress. Results of this study could help to improve management and the design of ICLF systems in view of sustainability and resistance to (water) crises but should be further supported by in depth analysis of soil water dynamics as affected by climate gradients, soil types and different management practices.


Introduction
The Brazilian Savannah, also known as the Cerrado, covers an area of 204.7 million ha (about 22% of the Brazilian territory). In the last few decades, about 50% of the natural Cerrado vegetation has been replaced by agricultural crops and pastures [1][2][3]. According to Macedo [4], the planted pasture area in the Cerrado amounts to 60 million ha, of which According to the World Reference Base for Soil Resources [19], the soil of the study area is classified as a Ferralsol. Soil samples were taken in 10,20,30,40,60 and 100 cm depth with coring rings (10 cm −3 ). The soil profile was created using an excavator in February 2017. The 1.2 m-deep and 1 m-wide hole was dug in Rep 2 ( Figure 2). For each depth, soil samples were taken with a coring ring in three of the four soil profile walls, resulting in three replicates per depth. The soil had a sandy clay texture (55% sand, 6% silt, 39% clay), a bulk density of 1.15 g cm −3 , a usable field capacity between soil moisture values of 15.5 Vol% (pF 4.2, permanent wilting point) and 42.5 Vol% (pF 1.8, field capacity), a hydraulic conductivity of 59.2 cm day −1 and a pH value (measured in water) of 4.5.

ICLF System
The ICLF system was established in 2008. The area was prepared using light to heavy tillage and limestone/gypsum application (3 and 1 t ha −1 , respectively) to increase the pH and to improve P availability for the crops. In early November 2008, NPK fertilizer (5:25:15) was applied at a rate of 400 kg ha −1 , followed by levelling with a disk harrow. The amount of fertilizer was determined through chemical soil analysis and applied according to the experience and guidelines of EMBRAPA Beef cattle. In late November 2008, the crop component (Glycine max cv. BRS 255RR), and in January 2009, the tree component (Eucalyptus grandis x Eucalyptus urophylla, H13 clone), were planted. The eucalyptus seedlings were transplanted in east-west oriented rows. Each replicate had an area of 1.4 ha, and consisted of two single tree rows at a distance of 22 m and 2 m distance between the trees within each row, resulting in a tree density of 227 trees ha −1 (Figure 2). After the soybean harvest in April 2009, the pasture component (Brachiaria brizantha cv. BRS Piatã) was established and beef cattle (Nellore heifers) were introduced to the system as the trees reached a diameter at breast height, at least 60 mm. The cattle grazed at varying stocking rates depending on biomass production in order to keep the height of the pasture at approximately 35 cm (put and take method). The management of the ICLF system involved According to the World Reference Base for Soil Resources [19], the soil of the study area is classified as a Ferralsol. Soil samples were taken in 10,20,30,40,60 and 100 cm depth with coring rings (10 cm −3 ). The soil profile was created using an excavator in February 2017. The 1.2 m-deep and 1 m-wide hole was dug in Rep 2 ( Figure 2). For each depth, soil samples were taken with a coring ring in three of the four soil profile walls, resulting in three replicates per depth. The soil had a sandy clay texture (55% sand, 6% silt, 39% clay), a bulk density of 1.15 g cm −3 , a usable field capacity between soil moisture values of 15.5 Vol% (pF 4.2, permanent wilting point) and 42.5 Vol% (pF 1.8, field capacity), a hydraulic conductivity of 59.2 cm day −1 and a pH value (measured in water) of 4.5.

ICLF System
The ICLF system was established in 2008. The area was prepared using light to heavy tillage and limestone/gypsum application (3 and 1 t ha −1 , respectively) to increase the pH and to improve P availability for the crops. In early November 2008, NPK fertilizer (5:25:15) was applied at a rate of 400 kg ha −1 , followed by levelling with a disk harrow. The amount of fertilizer was determined through chemical soil analysis and applied according to the experience and guidelines of EMBRAPA Beef cattle. In late November 2008, the crop component (Glycine max cv. BRS 255RR), and in January 2009, the tree component (Eucalyptus grandis x Eucalyptus urophylla, H13 clone), were planted. The eucalyptus seedlings were transplanted in east-west oriented rows. Each replicate had an area of 1.4 ha, and consisted of two single tree rows at a distance of 22 m and 2 m distance between the trees within each row, resulting in a tree density of 227 trees ha −1 (Figure 2). After the soybean harvest in April 2009, the pasture component (Brachiaria brizantha cv. BRS Piatã) was established and beef cattle (Nellore heifers) were introduced to the system as the trees reached a diameter at breast height, at least 60 mm. The cattle grazed at varying stocking rates depending on biomass production in order to keep the height of the pasture at approximately 35 cm (put and take method). The management of the ICLF system involved a crop rotation strategy every four years (i.e., cultivation of soybean as a crop for five months) followed by three and a half years of Brachiaria brizantha cv. BRS Piatã pasture. For this study, data was collected in the third year after pasture establishment during the second pasture cycle of the ICLF system from November 2015 until August 2017 at tree height of about 25 to 27 m. In January 2016, 50 kg N ha −1 (in the form of urea) and 300 kg NPK ha −1 (0:20:20) were applied to the pasture.

PAR, SM and AGBM
Photosynthetically active radiation (PAR), soil moisture (SM) and aboveground green grass biomass (AGBM) were measured along a transect of five sampling points between tree rows to assess the shading and soil moisture gradient created by the trees. These five sampling points were: P1S: 1 m and P6S: 6 m from the southern row; P1N: 1 m and P6N: 6 m from the northern row; P11: 11 m from the southern and the northern row respectively ( Figure 2). PAR (µmol m −2 s −1 ) was measured monthly with a Decagon AccuPAR LP-80 at clear sky, in the morning and in the afternoon directly above the grass canopy at each sample point between the tree rows, and on an open field (without any trees) beside the ICLF system. PAR values shown in this study represent means of measurements performed in the morning and afternoon.
At each sampling point in each replication, a fiberglass access tube (DeltaT Type: ALT1) of 1 m length was vertically installed for volumetric SM (Vol%) measurements with a portable DeltaT FDR (frequency domain reflectometry) PR2/6 profile probe. FDR measurements were taken weekly in depths of 10, 20, 30, 40, 60 and 100 cm. The FDR soil moisture measurements were validated against gravimetrical soil moisture measurements. Samples for gravimetrical soil moisture measurements were taken with an auger right next to every FDR access tube at the same depths where the FDR probe measurements were taken.
AGBM (g DW m −2 ) was quantified monthly. At each sampling point in each replication, 1 m 2 of grass biomass was harvested 5 cm from the ground. All harvested biomass samples were separated into green biomass and dead biomass, dried for 48 h at 70 • C and weighed to determine the dry weight biomass per area in g m −2 .

Data Analysis
Data was analyzed in a randomized complete block (three replicates) design using the generalized least squares model of R [20][21][22], with repeated measurements over time, and the seasons (spring, summer, autumn, winter) as a repeated factor. For the analysis of SM and AGBM sample points, season and block were considered as fixed effects. The most appropriate correlation structure for repeated measurements was chosen according to the Akaike Information Criterion (AIC). For AGBM the first order autoregressive (ar(1)) correlation structure was chosen. Fisher's LSD (least significant difference) test and significance at p < 0.05 were used for the comparison of mean values. The Pearson correlation coefficient (r) between PAR and AGBM and SM and AGBM was calculated also using R [20]. To normalize data for the correlation analysis, relative values of PAR, SM and AGBM along the gradient between the tree rows were calculated.

PAR between the Tree Rows during Different Seasons
PAR received by the grass canopy varied widely between the seasons, but also between the tree rows, especially in spring and summer ( Figure 3). In comparison to open field measurements, P11 received 9% less PAR in spring and 21% less PAR in summer, respectively. Furthermore, in both seasons, the sample points closest to the tree rows (P1S and P1N) received with about 200 µmol m −2 s −1 even less PAR than all other sample points in winter. In winter, PAR was more evenly distributed between the sample points, with P1S receiving 43% and P6S 74% less PAR compared to the open field. In autumn, P1S and P6S received 33% less PAR compared to the open field, P11 and P1N about 50% less, and P6N 78% less. In contrast to spring and summer, in autumn and winter, the sample points P6S and P6N received even slightly more PAR than the center sample points.

SM in ICLF between the Tree Rows during Different Seasons
Throughout the measurement period, monthly mean SM in the top 1 m soil layer was highest in the rainy season (November-April), and up to 10 Vol% lower in the dry season (May-October) for all sample points ( Figure 4). On average over the sample points, SM was lowest in August 2016, with 21.1 Vol%, and highest in May 2016, with 30.5 Vol%. Between the sample points, SM was lowest at P1N, with 23.2 Vol%, followed by P1S, with 25.0 Vol%, and highest at P11, with 30.1 Vol%, on average over the measurement period. Little difference was found between P6S and P6N, which had 27.2 Vol% and 27.9 Vol% SM on average over the measurement period.

SM in ICLF between the Tree Rows during Different Seasons
Throughout the measurement period, monthly mean SM in the top 1 m soil layer was highest in the rainy season (November-April), and up to 10 Vol% lower in the dry season (May-October) for all sample points ( Figure 4). On average over the sample points, SM was lowest in August 2016, with 21.1 Vol%, and highest in May 2016, with 30.5 Vol%. Between the sample points, SM was lowest at P1N, with 23.2 Vol%, followed by P1S, with 25.0 Vol%, and highest at P11, with 30.1 Vol%, on average over the measurement period. Little difference was found between P6S and P6N, which had 27.2 Vol% and 27.9 Vol% SM on average over the measurement period.
Analysis of variance revealed highly significant differences (p ≤ 0.001) between the months and between the sample points, but no significant difference between both years (Table 1).  Analysis of variance revealed highly significant differences (p ≤ 0.001) between the months and between the sample points, but no significant difference between both years (Table 1). During the first year and the second year, and at all SP, SM increased with depth, except at P6S, P11, and P6N in the first year and spring of the second year ( Figure 5). At these sample points, SM decreased at depths between 40 m and 60 cm, but increased at greater depths. Throughout the profile, SM was lower during the second year compared to the first year. During both years and for all seasons, SM next to the trees (P1S and P1N) was lower compared to the center point P11, as well as P6S and P6N. Further, SM at P1S was always slightly lower than at P1N. Throughout the profile, SM was lowest during  During the first year and the second year, and at all SP, SM increased with depth, except at P6S, P11, and P6N in the first year and spring of the second year ( Figure 5). At these sample points, SM decreased at depths between 40 m and 60 cm, but increased at greater depths. Throughout the profile, SM was lower during the second year compared to the first year. During both years and for all seasons, SM next to the trees (P1S and P1N) was lower compared to the center point P11, as well as P6S and P6N. Further, SM at P1S was always slightly lower than at P1N. Throughout the profile, SM was lowest during winter followed by spring, summer, and autumn during the first year, whereas during the second year, SM was lowest in spring followed by winter, summer and autumn. However, in the topsoil (10 cm to 30 cm depth), SM was lowest during winter in both years. the second year, SM was lowest in spring followed by winter, summer and autumn. However, in the topsoil (10 cm to 30 cm depth), SM was lowest during winter in both years.  Analysis of variance showed significant differences in SM between years (except for 60 cm depth), seasons and sample points at all measured depths (Table 2). Whereas in the topsoil at 10 cm depth, season had the largest effect on SM, the influence of the sample point became larger with increasing soil depth. Analysis of variance showed significant differences in SM between years (except for 60 cm depth), seasons and sample points at all measured depths (Table 2). Whereas in the topsoil at 10 cm depth, season had the largest effect on SM, the influence of the sample point became larger with increasing soil depth.

AGBM between the Tree Rows during Different Seasons
Grazing intensity differed largely between both years, with grazing activities from December until July, and the highest stocking rate (3.2 AU ha −1 ) in February during the first year and grazing in December only during the second year ( Figure 6). Consequently, in the first year, AGBM was with 114 g DW (dry weight) m −2 on average, substantially lower than in the second year, with 297 g DM m −2 .

AGBM between the Tree Rows during Different Seasons
Grazing intensity differed largely between both years, with grazing activities from December until July, and the highest stocking rate (3.2 AU ha −1 ) in February during the first year and grazing in December only during the second year ( Figure 6). Consequently, in the first year, AGBM was with 114 g DW (dry weight) m −2 on average, substantially lower than in the second year, with 297 g DM m −2 .  The ANOVA showed highly significant differences (p ≤ 0.001) between both years, months and sample points (Table 3). In both years, the highest AGBM was measured in January, but in the first year, AGBM drastically decreased from February to March, whereas in the second year, AGBM decreased more slowly. The largest difference between both years was found between March and August. In both years, between most of October The ANOVA showed highly significant differences (p ≤ 0.001) between both years, months and sample points (Table 3). In both years, the highest AGBM was measured in January, but in the first year, AGBM drastically decreased from February to March, whereas in the second year, AGBM decreased more slowly. The largest difference between both years was found between March and August. In both years, between most of October and February, AGBM was highest at P11, and lowest at P1S and P1N. However, between March and May in the second year, the highest AGBM was found at P1S and P6S. AGBM was highest during summer at all sampling points in both years (Table 4). In the first year, AGBM was lowest in autumn and winter, while in the second year, it was lowest in spring and winter. In the first year, AGBM at P11 was significantly higher than at P1S and P1N during summer, autumn and winter. During spring, no significant differences between the sample points were found. In year 2, AGBM was only significantly higher at P11 than at both points close to the trees during summer, while during spring and winter, significant differences were only found between P11 and one of both measurement points close to the trees.

Correlation between PAR and AGBM and SM 0-100 cm and AGBM between the Tree Rows
In spring and summer, AGBM was closely correlated with both, PAR and SM, with higher correlation coefficients with PAR (Table 5). In autumn and winter, the correlation between PAR and AGBM was negative and not significant, whereas the correlation between SM and AGBM was positive and in winter, also highly significant. Table 5. Pearson correlation coefficient for correlations between PAR (µmol m −1 s −1 ) and AGBM (g DW m −2 ) and SM 0-1000 mm (Vol%) and AG BM (g DW m −2 ) in spring, summer, autumn, winter. ***, **: significant at p-value ≤ 0.001, ≤ 0.01, respectively.

PAR, SM and AGBM Distribution throughout the Year
The trees in the ICLF system affected the amount of PAR received by the pasture, with a relative PAR reduction depending on the season and on the distance to the trees. The fact that in silvopastoral systems the pasture is shaded by the trees, is a well-known phenomenon with negative effects on both, light quantity and quality as red/far red ratio, received at the grass canopy [17,23]. The relative PAR reduction depends mostly on tree height, which was about 26 m at the beginning of our measurements and thus, led to a relatively large mean PAR reduction of 53% compared to above canopy PAR values (Table 6). Further, the magnitude of the effects depends on tree spacing, with spacing of more than 20 m between tree rows, as in our case, being beneficial for the relative light intensity in the median between tree rows [3,17]. In addition to tree height and spacing, time of day and time of year influence PAR reception at grass canopy level, related to the inclination of the sun. Light transmission at the center point between the tree rows was observed to be lower close to sunrise and sunset and during autumn and winter months [14]. In our experimental plots, PAR was reduced by more than 50% at the center point P11 during autumn and winter. Since PAR values represent means of measurements performed in the morning and afternoon, the sun did not rise above the tree canopy yet at the time of measurement during autumn and winter. However, a reduction in PAR near the trees was higher during spring and summer, whereas during autumn and winter, PAR was more evenly distributed between the tree rows. In absolute values, PAR received by the grass canopy close to the trees was up to 4 times higher in autumn and winter, than during spring and summer. Feldhake [23] observed peak levels of solar radiation under the trees at 25% cloud cover because of increased diffuse radiation. In our experiment, PAR measurements were performed only under clear sky conditions. However, we hypothesized that a high ratio of diffuse radiation between the tree rows during autumn and winter led to both, a low relative PAR reduction close to the trees in comparison to the center point P11, and higher absolute PAR value close to the trees during autumn and winter compared to spring and summer. As seasons were defined meteorologically, and not by days of solstice and equinox, spring (September-November) and summer (December-February) resulted in relatively similar PAR values. Equally, PAR distribution during autumn (March-May) and winter (June-August) were very similar. SM over a soil depth of 0-100 cm showed a distinct seasonal dynamic closely related to rainfall pattern (Figure 4). SM was significantly lower in the second year, since precipitation was substantially lower than in the first year (1020 mm vs. 1553 mm). Further, SM in the topsoil (0-30 cm) was generally lower and more dynamic, which has also been found in a similar experimental setup by Bosi et al. [15], who related this effect to water uptake by the pasture and evaporation from the soil, especially at the beginning of a drying period. Furthermore, in the topsoil, permanent wilting point was almost reached at some measurement points. For Brachiaria brizantha, more than 80% of the roots were found in the upper 30 cm of the soil, both under full irrigation and under water deficit [24], and plant biomass was significantly decreased by water deficit, when soil water levels reached 25% of field capacity [25]. At our study site, a soil water content of 25% of field capacity equals a volumetric soil moisture of 19.4%. During winter, soil moisture at 0-30 cm was often below this threshold and, thus, negatively affected plant biomass.
Intra-and inter-annual AGBM values varied largely to intra-and inter-annual climate variability and development related senescence processes of the grass. Depending on temperature, daylength, solar radiation, and rainfall [26], canopy senescence is faster than canopy renewal during the cool and dry winter period. With regard to development, in Brazil, the cultivar Piatã starts flowering in January and February [27], influenced by changes in photoperiod that stimulates floral induction [28]. Flowering and seed development generally serve as a cue for tiller senescence in grasses [29].
Despite the lower precipitation, more AGBM was observed in the second year, due to a significantly lower grazing intensity. In the first year, the pasture was grazed starting in December, which lead to a strong decline in AGBM until March, without recovery during autumn and winter. In the second year, grazing took place only during summer and a very low stocking rate was maintained. Under the conditions of the second year AGBM still peaked during January, but declined very slowly until August. The much higher AGBM in the second year could explain the lower top-soil values. Greater grass leaf area lead to increased transpiration, combined with less rainfall.

PAR, SM and AGBM Gradients between the Tree Rows
During all seasons and in both years, the highest AGBM was found in the center point (P11), with a gradient of decreasing values towards the tree rows. According to Wilson and Wild [30], most tropical grasses decrease AGBM under shading almost proportionally to the amount of shade, provided water and nutrients are not limiting. Lower grass biomass as a result of lower light intensity under and next to the trees in silvopastoral systems has been reported by several authors [3,29,30]. In our study, a PAR gradient between the tree rows was found during spring and summer, but with large PAR reductions under the trees relative to the center point P11 of 72% and 86% in spring and summer, respectively, whereas AGBM was reduced under the trees by 48% in spring and 42% in summer. Combining measurements from different seasons and at different tree densities, Santos et al. [3] found in a similar ICLF system under similar climatic conditions that for every 1% reduction of PAR there was a 1.35% decrease of forage dry mass. Since grass growth is subject to strong seasonal variations related to both weather and development, we suggest that the effect of PAR reduction on biomass should be analyzed for each season individually. In our case, every 1% reduction of PAR resulted in a 0.67% reduction in AGBM during spring, but only a 0.49% reduction during summer. However, during autumn and winter the pattern of PAR distribution between the tree rows was by far less pronounced, and the observed differences in AGBM could not be explained by PAR reductions under the trees. Pasture growth responses to shading also depend on N management. At high N doses, Brachiaria species showed a stronger response to light intensity in terms of tiller number [31]. In this study, 50 kg/ha N was applied at the beginning of the measurements. Therefore, we assume there was a medium response of biomass to light intensity.
Although potential water loss via transpiration from the grass canopy was highest at the center point, due to higher PAR (especially during spring and summer) and greater AGBM, SM was highest at P11. Soil water loss from the upper 1 m during a dry spell and soil water recharge during a rainy period was higher closer to the trees at points P1S and P1N (Figure 7). Soil drying and wetting cycles lead to the greatest SM variability occurring in the subsoil near the tree lines, rather than the topsoil, the interface to the atmosphere and the main rooting zone for grasses. Surprisingly, the center subsoil areas at 60-100 cm depth were the most conservative water holding areas of the ICLF. Topsoil SM variability showed an intermediate response to a drying and wetting cycle. Although Eucalyptus species are considered an intensive consumer of water [32], our results show that highest fluctuations of SM are initiated in subsoil horizons, which are often considered as buffers to fluctuations created by drying and wetting cycles [33,34]. This buffer function was still visible in the center subsoil area between tree lines, where the lowest fluctuations of SM were observed throughout drying and wetting cycle (Figure 7).
Shaded forage grasses can increase specific leaf area, leaf elongation rate, and leaf length to increase light interception [31,35,36], leading to a larger leaf area, despite a lower AGBM. However, it is unlikely that the higher soil water loss close to the trees is caused by higher transpiration of the grass next to the trees. Probably, the Eucalyptus trees depleted the soil water at P1S and P1N due to the large share of their fine roots found in the topsoil [37]. Eucalyptus tree roots can grow deeper than 3 m within one year, but fine root density decreases sharply with depth [38]. Eucalyptus is a C3-species and thus, has a lower water use efficiency than the C4 species Brachiaria [39]. Higher water uptake near the trees in silvopastoral systems was found by several authors for several species e.g., in ICLF systems with Brachiaria brizantha and Eucalyptus urograndis [15], or with Brachiaria decumbens and Brazilian native trees [14], in an intercropping system with switchgrass and loblolly pine [40] and in a Panicum maximum pasture with Eucalyptus argophloia [41].
However, not only water loss, but also soil water recharge during rainy periods was found to be higher closer to the trees. It has been argued that this effect is a result of large water interception by the trees and its subsequent deposition into the soil [15]. However, measurement of intercepted rainfall between the tree rows could not confirm this hypothesis as cumulative precipitation was roughly equal across the points of measurement between the tree rows ( Figure 8).   lolly pine [40] and in a Panicum maximum pasture with Eucalyptus argophloia [41]. However, not only water loss, but also soil water recharge during rainy periods was found to be higher closer to the trees. It has been argued that this effect is a result of large water interception by the trees and its subsequent deposition into the soil [15]. However, measurement of intercepted rainfall between the tree rows could not confirm this hypothesis as cumulative precipitation was roughly equal across the points of measurement between the tree rows ( Figure 8). Since soil water content and thus, soil water potential, was generally lower closer to the tree rows, we hypothesize that lateral subsurface water movements led to differences in soil water recharge. However, hydraulic conductivity of the soil at the experimental site was estimated to be 59.2 cm day −1 , which is probably too low to fully explain our observations. The experimental site has a slope of 2% in north-south direction, which makes it perpendicular to the tree rows, and means that lateral subsurface water movements could be facilitated by gravity, which could also explain the slightly higher soil moisture at P1S than at P1N.
During summer and spring, AGBM was strongly correlated with PAR, whereas the correlation between AGBM and SM was weaker, which indicates a stronger limitation of AGBM by radiation than by SM. Considering the relatively low SM during spring, a mainly radiation-driven growth of Brachiaria brizantha illustrates the good drought tolerance of this species [42]. A positive relationship between grass biomass and radiation in Since soil water content and thus, soil water potential, was generally lower closer to the tree rows, we hypothesize that lateral subsurface water movements led to differences in soil water recharge. However, hydraulic conductivity of the soil at the experimental site was estimated to be 59.2 cm day −1 , which is probably too low to fully explain our observations. The experimental site has a slope of 2% in north-south direction, which makes it perpendicular to the tree rows, and means that lateral subsurface water movements could be facilitated by gravity, which could also explain the slightly higher soil moisture at P1S than at P1N.
During summer and spring, AGBM was strongly correlated with PAR, whereas the correlation between AGBM and SM was weaker, which indicates a stronger limitation of AGBM by radiation than by SM. Considering the relatively low SM during spring, a mainly radiation-driven growth of Brachiaria brizantha illustrates the good drought tolerance of this species [42]. A positive relationship between grass biomass and radiation in silvopastoral systems was also found by Santos et al. [3] and Silva-Pando et al. [43]. However, we found that, during the winter, AGBM was solely correlated with SM, but not with PAR, indicating a water limitation. It has been suggested that in silvopastoral systems with a comparable climate to our site, a PAR limitation is the major constraint to production of dry matter when there are no limitations caused by soil water deficit [18]. As during winter, SM next to the trees was below 25% of the field capacity, which has been shown to be critical for Brachiaria growth [25]. The positive correlation between AGBM and SM between the tree rows supports the finding of a water limitation of the pasture next to the trees during the dry season of the year.
Trees provide shade for animals, which means that ICLF systems can also improve animal thermal comfort indices by reducing full sun exposure [44]. Assuming animals are seeking comfort in the shade, they likely spend more time close to the trees during spring and summer. Therefore, we suggest that apart from lower light intensity and lower soil moisture, grass biomass next to the trees is also impaired by trampling and probable higher grazing intensity. However, a considerable accumulation of nutrients from dung and urine would also accumulate next to tree lines. Additionally, biomass inputs from foliage, including allelopathic substances depending on the species, could affect understory plant growth. During the second year, grazing intensity in summer was much lower than in the first year, while during autumn and winter no grazing took place. To support the idea of higher grazing and trampling intensity under the trees, in the second year, smaller differences in AGBM between the sampling points P11 and P1S/P1N were found from summer to winter. While during the first year, compared to P11, AGBM next to the trees was 47%, 55% and 52% lower during summer, autumn and winter, respectively, it was only 37%, 0% and 24% lower during summer, autumn and winter of the second year. As we can only compare two years of different grazing intensities, no valid conclusions can be drawn as to the impact of animal presence on the grass biomass gradient between the tree rows, including possible interactions with soil water dynamics.
ICLF systems offer a promising land use strategy to address the central challenges to future agricultural production systems [8]. The available management options regarding combinations of different components and practices for ICLF systems are numerous, and indicators are highly requested to guide farmers towards a sustainable land use.
While several studies show the beneficial effects and synergies created by integrated systems with attributes of reduced soil degradation [10], greenhouse gas emission (GHG) mitigation [45,46], and increased soil carbon sequestration [47], the effects on the water balance have yet to be explored. In view of a changing climate, water and its management will be seen as a critical aspect for integrated agricultural production systems. However, a sustainable management approach that takes into account water-related challenges requires a detailed knowledge base, including a comprehensive overview of how land use practices and environments are altering water pathways within ICLF systems.
Our results from a mature ICLF system in central-west Brazil suggest that the system's seasonality and resource cycles, such as those observed in episodic wetting and drying events, strongly interact with a management impact. Consequently, for plant water availability, temporary water stress may develop in deeper soil horizons near tree rows, where soil water content fluctuations were surprisingly more pronounced compared to topsoils across the analyzed gradient. These results stand in contrast to the generally held view that topsoils are usually subject to greater shifts in SM compared to the subsoil [15,34]. Deeper soil horizons in the center position between tree rows, however, showed a conservative response during a drying and wetting cycle, thus likely functioning as buffer zone for soil water resources. For the management of soil water reserves, improving the system's resilience to seasonal drought or periods of water stress, tree row distances and height, age (i.e., shading) as well as root architecture should be carefully considered and adapted to local conditions. Management of water resources will certainly at the same time translate into processes triggering key-parameters affecting sustainability such as greenhouse gas emissions and C-sequestration of ICLF systems. The identification of tipping points, where stress and/or land use practices push systems towards a nonsustainable response is essential in the development of best practices for a responsive or regenerative management of ICLF systems. This would require monitoring of key parameters in high spatial and temporal resolution. Therefore, we advocate increased research on water-related functions and processes of ICLF systems to support and fully develop the promising land use potentials offered by these systems to improve sustainability of agricultural production.

Conclusions
Spatial and temporal water availability in mature ICLF systems is dynamic, affected by climate seasonality, management practices, and system design. Throughout the year soil moisture was lower closer to the tree rows, but critical values reaching the wilting point were only observed during winter (dry season) and in the topsoil. Seasonal soil moisture pattern and radiation input suggest that grass productivity is shifting from light limitation during summer to a water limitation during winter, especially near tree lines. Notably, soil drying and wetting cycles lead to highest SM variability in the subsoils near the tree lines, and not, as expected, in topsoils, which are considered the interface for evaporation to the atmosphere. Relatively wet subsoils in the center position between tree lines were most conservative regarding SM fluctuations, suggesting these areas function as important water reservoirs throughout the year, and during drought events. The system's resilience to