Oil-palm management alters the spatial distribution of amorphous silica and mobile silicon in topsoils

. Effects of oil-palm ( Elaeis guineensis Jacq.) management on silicon (Si) cycling under smallholder oil-palm plantations have hardly been investigated. As oil palms are Si accumulators, we hypothesized that management practices and topsoil erosion may cause Si losses and changes in spatial Si concentration patterns in topsoils under oil-palm cultivation. To test this hypothesis, we took topsoil samples under mature oil-palm plantations in well-drained and riparian areas of Jambi Province, Indonesia. The samples were taken from four different management zones within each oil-palm plot: palm circles, oil-palm rows, interrows, and below frond piles. We quantiﬁed mobile Si (Si M ) and Si in amorphous silica (Si Am ) by the extraction of CaCl 2 and NaCO 3 , respectively. Both fractions are important Si pools in soils and are essential for plant–


Introduction
The transformation of lowland rainforests into cash-crop plantation systems (e.g.timber, rubber, and oil palm) involves vast expansion of oil-palm monocultures in Jambi Province, Indonesia (Drescher et al., 2016;Tsujino et al., 2016).To date, smallholder farmers manage ∼ 40 % of oilpalm plantations in Jambi Province (Euler et al., 2016), and palm oil remains a tropical cash crop with high demand on the global market (Schaller et al., 2020).Oil-palm cultivation has improved the livelihoods of many smallholder farmers yet at the expense of the natural environment (Clough et al., 2016;Grass et al., 2020;Qaim et al., 2020), leading to a decrease in biodiversity (Drescher et al., 2016;Meijaard et al., 2020) and ecosystem services (Dislich et al., 2017).Due to these "trade-offs" (Grass et al., 2020) and a global interest in reducing deforestation (Tsujino et al., 2016), much research focuses on identifying ways to increase land-use sustainability while keeping current oil-palm plantations profitable (Darras et al., 2019;Luke et al., 2019).
Under humid tropical climate conditions, intense silicate weathering and element leaching from soils take place, including leaching of silicon (Si), i.e. desilication (Haynes, 2014).Farmers commonly apply nitrogen-phosphorouspotassium (NPK) fertilizers and lime to maintain an adequate plant nutrition and soil pH (Darras et al., 2019).However, Si also plays an important role in terrestrial biogeochemical cycling (Struyf and Conley, 2012) and enhances crop production in several ways (Epstein, 2009;Guntzer et al., 2012).In soils, silicic acid can mobilize phosphate by occupying anion adsorption sites.Si also mitigates plant toxicity by binding toxic cations such as aluminium (Al), cadmium (Cd), and arsenic (As) that become mobile at low soil pH (Street-Perrott and Barker, 2008;Schaller et al., 2020).Furthermore, Si can increase drought resistance of plants (Schaller et al., 2020).Silica precipitates in cell walls, cell lumen, and intercellular spaces of leaves and can reduce transpiration (Epstein, 2009).In Si-depleted soils, some crops, including oil palms, can thus benefit from Si fertilization (Klotzbücher et al., 2018).
In terrestrial ecosystems, Si cycling is mostly driven by two Si pools: mobile Si in soil solution (Si M ) and Si present in amorphous silica (Si Am ) (Struyf et al., 2010;de Tombeur et al., 2020).Si M is the Si fraction that is readily available to plants and usually present as monomeric silicic acid (H 4 SiO 4 ) in terrestrial environments (Georgiadis et al., 2013).Si Am is the largest non-mineral Si pool in soils (Barão et al., 2014;Unzué-Belmonte et al., 2017).Its solubility exceeds that of silicate minerals by several orders of magnitude (Iler, 1979;Fraysse et al., 2009).Si Am in soils can be subdivided into Si Am of biogenic origin and of pedogenic origin.The first mainly includes Si in phytoliths, i.e. small bio-opal bodies precipitated in plant tissues that are released during plant-litter decomposition (Barão et al., 2014;Clymans et al., 2015;Schaller et al., 2021).Soil microorganisms (testate amoebae, sponges, diatoms) contribute to a lesser extent (Schaller et al., 2021).Si Am of pedogenic origin, i.e. silica precipitated from soil solution, mainly occurs as soil-particle coatings and void infillings (Schaller et al., 2021).Si Am in topsoils is predominantly of biogenic origin (Clymans et al., 2015;Schaller et al., 2021), whereas Si Am in subsoils is mostly of pedogenic origin (Schaller et al., 2021).
Ecosystem Si cycling can be altered by human impact such as deforestation (Conley et al., 2008), land-use and landcover (LULC) change (Struyf et al., 2010;Barão et al., 2020), and fire (von der Lühe et al., 2020;Schaller and Puppe, 2021).After LULC transformation from forest to arable land, Si can be lost from the system through harvest, topsoil erosion, and increased Si leaching in soil.Si leaching in soil is triggered by reduced interception which results in increased percolation (Keller et al., 2012;Vandevenne et al., 2012;Kraushaar et al., 2021).Si-accumulating plants such as rice, wheat, barley, maize, and oil palm (Ma and Takahashi, 2002;Munevar and Romero, 2015) are characterized by Si accumulation of >1 wt % in dry leaf tissue and a Si / Ca ratio >1 (Ma and Takahashi, 2002).Such Si accumulators may accelerate Si turnover at the soil-plant interface by taking up high amounts of Si from soil solution and returning Si-rich litter to soils (Struyf andConley, 2009, 2012).In oil-palm plantations, we therefore expected Si losses by harvest and topsoil erosion (Vandevenne et al., 2012;Munevar and Romero, 2015).In addition, we expected that the spatial arrangement of oil-palm rows and interrows -with frond piles (frond pile) or without ("empty" interrow) -results in a corresponding spatial Si concentration pattern in topsoils.
Oil palms are planted in rows (Kotowska et al., 2015) (Fig. 1a).A distance is kept between the rows to ensure sufficient light exposure (Corley and Tinker, 2016).The space between two oil-palm rows is referred to as an interrow.They serve either as harvesting paths or as deposition sites for cutoff palm fronds that are stacked up in long, flat piles (Corley and Tinker, 2016).Fertilizers are only applied within a circle of ∼ 1.5-2 m around the palm stem (palm circle) (Munevar and Romero, 2015;Formaglio et al., 2020).In addition, nutrients are released from decaying plant litter.Thus, we hypothesized that Si is mainly released and returned to soils in the form of biogenic Si Am under frond piles, leading to higher topsoil Si Am concentrations, while other management zones (including palm circles, oil-palm rows, and empty interrows) might be at risk of Si depletion.
Furthermore, we hypothesized that in oil-palm plantations established on sloping terrain, Si is removed by topsoil erosion in scarcely vegetated interrows.We assumed that phytoliths might be even more prone to erosion than mineral soil particles because of their lower density, leading to a disproportionately high Si Am loss through topsoil erosion.Such additional Si Am loss from interrows would be unfavourable, as interrows may serve as new planting sites in a subsequent plantation cycle after ∼ 25 years (Corley and Tinker, 2016).Thus, our study aimed at assessing the impact of management practices in smallholder oil-palm plantations on Si cy-cling.In addition, we considered it important to account for potential differences in the intensity of natural desilication in different landscape positions.Therefore, we carried out the same study in two different landscape positions associated with differing water regimes: (i) in well-drained areas with presumably high desilication rates and (ii) in riparian areas where we assumed that regular flooding might involve an input of Si dissolved in stream water into the system, partially compensating for desilication.

Study area and sites
The study was associated with the interdisciplinary Collaborative Research Centre (CRC) 990, funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), addressing environmental and socioeconomic impacts of rainforest conversion into plantation systems in Sumatra, Indonesia (Drescher et al., 2016;Dislich et al., 2017).Thus, it was conducted on CRC 990 plots in smallholder oil-palm plantations in the Harapan landscape of Jambi Province, Sumatra, Indonesia (1 • 55 0 S, 103 • 15 0 E; 50 ± 5 m a.s.l.).Geologically, this lowland landscape is characterized by pre-Paleogene metamorphic and igneous bedrock that is overlain by lacustrine and fluvial sediments (De Coster, 1974) in which predominantly loamy mineral soils have formed (Allen et al., 2016).Preliminary results showed that quartz, kaolinite, and Fe-Al oxides are the most abundant minerals in these highly weathered soils.In our study area, Acrisols are present in well-drained areas, found at higher elevation and on sloping terrain.Stagnosols and stagnic Acrisols dominated in seasonally flooded riparian plots, i.e. in floodplains (Hennings et al., 2021).The Harapan region is characterized by a humid tropical climate (Af in the Köppen-Geiger classification) with a mean annual temperature of 26.7 • C and a mean annual precipitation of 2230 mm (Drescher et al., 2016).The rainy season has two precipitation maxima: one in December and another one in March.A dry period lasts from July to August (Drescher et al., 2016).The natural vegetation is a mixed dipterocarp lowland rainforest (Laumonier, 1997) which is nearly only preserved in the Harapan rainforest, an ecosystem restoration area in the south of Jambi Province (Harrison and Swinfield, 2015), and in the Barisan mountains in the west of Jambi Province (Drescher et al., 2016).In addition to oil-palm plantations, other important land-use systems in Jambi Province include rubber plantations and agroforestry systems (Dislich et al., 2017).

Topsoil samples
From April to August 2018, topsoil sampling was conducted in four well-drained (HO1-HO4) and four riparian plots (HOr1-HOr4).Oil palms were planted between 1997 and 2001 in well-drained plots and between 1998 and 2008 in riparian areas (Hennings et al., 2021) following a triangular planting scheme with ∼ 9 m distance between the stems (Fig. 1a).Interrows were used either as harvesting paths or to stack cut-off palm fronds (frond pile) (Kotowska et al., 2015).In plot HO1, every interrow contained frond piles.Thus, topsoil samples of interrows were obtained only from three well-drained plots.The understorey vegetation of all well-drained plots was occasionally weeded.Two riparian plots (HOr1 and HOr2) had a well-maintained grass cover between the oil palms.
In each of the eight plots, topsoil samples were taken with steel cylinders (height = 4 cm, volume = 100 cm 3 ) at five locations along the slope.At each location, topsoils were sampled from four different management zones, i.e. (1) palm circle, (2) oil-palm row, (3) interrow, and (4) frond pile, to assess spatial patterns of Si Am and Si M concentrations in topsoils within the oil-palm plantations (Appendix Tables A1  and A2).Interrow topsoil samples were taken at a maximum distance between oil palms.The samples were dried (40 • C, 24 h) and sieved (≤ 2 mm) prior to Si analyses.An aliquot of each sample was dried at 105 • C to determine the water content of the samples dried at 40 • C.

Sediment traps
Sediment traps were installed in sets of two in interrows of the well-drained plots HO1-HO4 on 8-12 • sloping land (Sinukaban et al., 2000).Each trap consisted of a rectangular aluminium frame (2 × 1 m, 2 m 2 ).Its downslope-facing short side was funnel-shaped, directing surface runoff and eroded soil material into a bucket (Fig. 1b and Appendix Table A3).A second bucket was connected to the first bucket by a 2 cm thick tube to catch potential overflow.The traps were checked and maintained weekly from the beginning of September 2018 to the end of August 2019.The understorey vegetation in the sediment traps was kept in place to ensure that the understorey vegetation was representative of the oilpalm plantations.Both sediment traps in HO1 were manually weeded after 6 months because inside the traps, vegetation covered nearly 100 % of the soil surface, impeding topsoil erosion.Eroded soil material was collected whenever present, dried (40 • C, 48 h), sieved (≤ 2 mm), and weighed prior to Si analyses.Samples of eroded soil material from plot HO2 were excluded from further analysis because both traps were contaminated by crude oil.Losses of Si Am were calculated for each sediment trap by multiplying the concentration of Si Am of each sediment sample by the amount of eroded soil material collected by each trap (Appendix Table B3).Erosion estimates were determined for each trap by summing up the amount of eroded soil material for the 12month period from the beginning of September 2018 until the end of August 2019.Precipitation data of the two closest meteorological stations were used for correlating the observed soil erosion with precipitation.Distances between meteorohttps://doi.org/10.5194/soil-9-169-2023SOIL, 9, 169-188, 2023  et al., 2000).The downslope funnel-shaped part of the aluminium frame (F) directs surface runoff, together with eroded soil material, into a bucket (B1) that is connected to a second bucket (B2) by a 2 cm wide tube (T).Photos of the sediment traps are shown in the Appendix Table A3.
logical stations and plots comprised ∼ 2 km for HO1, ∼ 3 km for HO2, ∼ 8 km for HO3, and ∼ 6 km for HO4.At each meteorological station, precipitation was measured by two automated precipitation transmitters (Thies Clima, Göttingen, Germany), at a height of 1.5 m and a horizontal distance of about 6 m.
2.3 Determination of silicon pools in topsoils 2.3.1 Silicon in amorphous silica (Si Am ) Silicon in amorphous silica (Si Am ) was extracted from topsoil samples and eroded soil material by 1 % Na 2 CO 3 solution (Meunier et al., 2014).At 85 • C, amorphous silica dissolves within 2-3 h in 1 % Na 2 CO 3 solution, thereby rapidly raising the Si concentration in solution.Once amorphous silica is completely dissolved, the release of Si to solution is only sustained by the slower dissolution of silicate minerals, which follows a linear trend.The Si concentration was measured four times during the linear dissolution phase.A linear equation was fitted to the data.The Si Am concentration was inferred from the y intercept of the linear regression.
In detail, 40 mL of 1 % Na 2 CO 3 solution was added to approximately 30 mg of soil material.The samples were then placed into a shaking water bath at 85 • C. To ensure steady Si release from topsoils, the samples were manually shaken at time intervals of 45 min.Aliquots were taken after 3, 3.75, 4.5, and 5.25 h.For this purpose, the samples were taken out of the water bath, cooled in a cold-water basin (10 min), and centrifuged (5 min, 3000 rpm).A 0.25 mL aliquot was taken from the supernatant of each sample and neutralized with 2.25 mL 0.021 M HCl.The Si concentrations in the aliquots were analysed by the molybdenum blue method (Grasshoff et al., 2009) using a UV-VIS spectrophotometer (Lambda 40, Perkin Elmer, Germany) at 810 nm.We chose 1 % Na 2 CO 3 as an extractant and used the extraction method by Meunier et al. (2014) instead of the stronger extractant of 0.1 M NaOH used by Barão et al. (2015) because we assumed that most Si in topsoils is of biogenic origin and dissolved well by Na 2 CO 3 (Meunier et al., 2014).

Mobile silicon (Si M )
Mobile silicon (Si M ) was extracted by a CaCl 2 solution, which provides electrolytes resembling natural soil solutions (Sauer et al., 2006;Georgiadis et al., 2013).From each sample, 1 g of soil material was mixed with 5 mL of 0.01 M CaCl 2 and left for 24 h, shaking for 1 min h −1 on an overhead shaker.Samples were centrifuged (5 min, 3000 rpm) and the supernatant was filtered through ash-free paper filters (1-2 µm).The Si concentrations were analysed in filtrates by the molybdenum blue method.We transformed the measured Si concentration (µg g −1 ) into the amount of Si M per gram 105 • C dried soil.

Statistical analyses
Statistical analyses were conducted on the grand means of topsoil Si concentrations in each water regime and management zone.The two latter were grouped into (i) palm circles in well-drained/riparian areas (each, n = 4), (ii) oil-palm rows in well-drained/riparian areas (each, n = 4), (iii) interrows in well-drained (n = 3)/riparian areas (n = 4), and (iv) frond piles in well-drained/riparian areas (each, n = 4).The four management zones were tested for significant differences in topsoil Si concentrations, within both the welldrained and riparian areas.In addition, we tested the welldrained and riparian areas for significant differences in topsoil Si concentrations by comparing the same management zone under two different water regimes.The data were log transformed to assert normal distribution (Shapiro-Wilk test) and homogeneity of variances (Levene test).Both criteria were met for all groups except for Si M in topsoils of oilpalm rows in well-drained areas (Appendix Table B4).We conducted a one-way analysis of variance (ANOVA) to detect if Si Am and Si M concentrations in topsoils of different management zones differed significantly within well-drained and riparian areas, as well as between well-drained and riparian areas.Then we used the Tukey-Kramer post hoc test to identify which management zones differed significantly.The level of significance was set at p ≤ 0.05.We used the open-source software R version 3.6.2and R CRAN packages ggplot2 (Wickham, 2016), ggpubr (Kassambara, 2022), car (Fox and Weisberg, 2019), and psych (Revelle, 2022) to perform these statistical analyses.
In well-drained plots, mean topsoil Si M concentrations were about twice as high under frond piles (13.68 ± 6.54 µg g −1 ) and in palm circles (11.17 ± 5.42 µg g −1 ) compared to oil-palm rows (6.38 ± 2.85 µg g −1 ) and interrows (5.62 ± 0.10 µg g −1 ) (Fig. 2c).Only plot HO1 showed exceptionally high topsoil Si M concentrations in oil-palm rows (outlier), which could be attributed to the dense vegetation throughout that smallholder plantation.In riparian plots, mean topsoil Si M concentrations were twice as high under frond piles (19.56 ± 6.13 µg g −1 ) compared to mean topsoil Si M concentrations in palm circles, oil-palm rows, and interrows, which all range around 11 µg g −1 (Fig. 2d).Mean topsoil Si M concentrations did not differ significantly (p ≤ 0.05) between the other management zones within the same water regime (well-drained/riparian), nor did mean topsoil Si M concentrations (in the same management zone) differ between water regimes.

Topsoil erosion and associated losses of Si Am
In plots HO3 and HO4, median Si Am concentrations in topsoils of interrows (1.53-1.57mg g −1 ) were roughly twice as high as in eroded soil material (0.66-0.88 mg g −1 ) (Table 1).In plot HO1, the median Si Am concentration in eroded soil material (1.61 mg g −1 ) was twice as high as in eroded soil material of plots HO3 and HO4 (0.66-0.88 mg g −1 ).Over the entire sampling period of 12 months, the four sediment traps in plots HO1 and HO4 indicated erosion rates of ∼ 4-5 Mg ha −1 yr −1 (Table 2).In plot HO3, a similar erosion rate was obtained from trap 1 (∼ 6 Mg ha −1 yr −1 ), whereas the erosion rate observed in trap 2 of plot HO3 was twice as high (∼ 12 Mg ha −1 yr −1 ).The Si Am losses through topsoil erosion amounted to 6-9 kg ha −1 yr −1 in the four sediment traps of HO1 and HO3, and 5-7 kg ha −1 yr −1 in both sedi- In plot HO1, a dense cover of mosses, grasses, and 20-50 cm high understorey vegetation prevented soil loss from September 2018 until end of January 2019 (Table 1, Fig. 3a).After manually weeding plot HO1 at the end of January 2019, the vegetation coverage was kept minhttps://doi.org/10.5194/soil-9-169-2023SOIL, 9, 169-188, 2023 imal (around 5 %).Noticeable losses of soil and corresponding losses of Si Am occurred between February (13-21 g m −2 of sediment, 16-53 mg m −2 of Si Am ) and the end of May 2019 (16-100 g m −2 of sediment/38-192 mg m −2 of Si Am ) (Fig. 3a and b).In plot HO3, scarce understorey vegetation of herbaceous plants (no grasses and mosses) covered about a third of the sediment traps (Table 1).Soil and corresponding Si Am losses were recorded continuously from September 2018 to May 2019 (Fig. 3).Each week, losses of topsoil material amounted to 4-62 g m −2 (corresponding to 1-90 mg m −2 Si Am ) (Fig. 3a, b).At three sampling dates, one in December 2018 and two in February 2019, peak soil losses ≥ 150 g m −2 occurred.The corresponding Si Am losses of these sampling dates were ≥ 90 mg m −2 , hence also representing among the highest Si Am losses throughout the sampling period.In plot HO4, vegetation coverage in the traps increased from 40 % in September 2018 to 60 % in May 2019 (Table 1).Soil loss occurred from mid-September 2018 to the end of May 2019 (Fig. 3).Losses of eroded soil material barely exceeded 50 g m −2 of sediment.However, an event with approximately 20 g m −2 of soil loss had corresponding Si Am losses ranging from 5-160 mg m −2 , thus showing a large variability.

Spatial topsoil Si Am concentration patterns
In oil-palm plantations, cut-off palm fronds stacked in every second interrow represent the main source of phytoliths (Albert et al., 2015;Clymans et al., 2015;Huisman et al., 2018), although these frond pile areas may occupy less than 15 % of the plantation area (Tarigan et al., 2020).Once phytoliths are released into topsoils, they can replenish the topsoil Si Am pool.Therefore, we had hypothesized that Si Am is mainly returned to soils under frond piles (see Sect. 1), leading to a spatial topsoil Si Am pattern with higher Si Am concentrations under frond piles.This hypothesis was corroborated for welldrained plots as topsoil Si Am concentrations were indeed significantly higher (2-fold) under frond piles (∼ 4 mg g −1 ) than in all other management zones (∼ 2 mg g −1 ) (Fig. 2a, Appendix Table B2a).Lower Si Am levels in palm circles, oilpalm rows, and interrows may reflect the pedogenic Si Am pool with only minor contributions of biogenic Si Am , e.g. from grass phytoliths.A possible reason for this is that decaying palm fronds are not returned to these management zones.
In the riparian plots, topsoil Si Am concentrations were equally high under frond piles and in interrows (∼ 3 mg g −1 ).This can only be explained by an additional important source of topsoil Si Am in interrows that was present in the riparian plots.The only potential Si Am source includes litter of grasses (Poaceae) and sedges (Cyperaceae) which also release considerable amounts of phytoliths upon their decomposition.Grasses and sedges are considered effective Si ac-cumulators as well (Blecker et al., 2006;Quigley et al., 2017).This explanation is further supported by the significantly lower (p ≤ 0.05) topsoil Si Am concentrations in palm circles (∼ 1.4 mg g −1 ) (Fig. 2b, Appendix Table B2b).Palm circles are weeded and treated with herbicides regularly.Thus, this management zone also lacks litter return and with that, a principal source of Si.The significant difference in topsoil Si Am concentrations between interrows and palm circles can only be explained by the presence or absence of grasses as phytolith sources.This observation highlights the importance of grasses and sedges in oil-palm plantations as they can also maintain soil-plant-Si cycling in the system.Thus, our original hypothesis that Si Am is mainly returned to the soils under frond piles, while topsoils in other management zones tend to be depleted in Si Am , is valid only in oil-palm plantations with a negligible grass cover.
The absence of any significant differences in topsoil Si Am concentrations between the two differing water regimes suggests that there was no noticeable Si supply by stream water to topsoils in riparian areas.In fact, release rates of biogenic Si Am from decaying oil-palm and other litter must be similar in both water regimes, so too the rate at which oil palms take up Si from soil solution and form phytoliths.This result contrasts with Vander Linden and Delvaux (2019) and Georgiadis et al. (2017), who found that soil type and soil properties affect Si cycling when comparing well-drained and floodplain soils: Vander Linden and Delvaux (2019) observed that flooding can temporarily increase the soil pH of acidic soils to ∼ 6.5-7.2.In this pH range, phytoliths dissolve faster (Vander Linden and Delvaux, 2019).Georgiadis et al. (2017) pointed to the effect of alternating redox conditions in soils as caused, e.g. by perched water, a fluctuating groundwater table, and flooding.When pedogenic Fe, Al, or Mn oxides and hydroxides are exposed to reducing conditions during flooding, they dissolve and release occluded Si into soil solution.After flooding, when the soil is exposed to oxidizing conditions again, Si can be occluded in and adsorbed to the surfaces of newly formed pedogenic oxides and hydroxides (Georgiadis et al., 2017).Georgiadis et al. (2017) found that these redox-induced dynamics affected mainly the following Si fractions: Si in soil solution (mobile Si), Si adsorbed to Fe oxides and hydroxides, and Si occluded in Fe oxides and hydroxides.Many other researchers found soil pH, soil texture, and soil chemistry to govern biogenic and pedogenic Si pools down to at least 1 m soil depth (Alexandre et al., 1997;Struyf and Conley, 2009;Li et al., 2020).The latter soil characteristics were kept constant in our study as we focused specifically on the effect of flooding.The reason that we did not detect any effect of flooding might be the advanced weathering and desilication status of the soils in our study area, which may have led to overall low Si levels of the investigated soils.2015), who compared δ 13 C values in soil profiles on the same well-drained oil-palm plantations of our study region.High erosion rates are to be expected, as oil-palm plantations have a rather open canopy compared to rainforests, permitting raindrops to directly hit the ground (Oliveira et al., 2013;Corley and Tinker, 2016).
During heavy rainfalls, raindrops release kinetic energy that breaks up soil aggregates, especially when hitting bare soil.Mobilized fine and broken aggregates can fill soil pores, thereby reducing infiltration, promoting surface runoff (Oliveira et al., 2013;Tarigan et al., 2020), and hence promoting erosion.Besides, soil compaction may be substantial in oil-palm interrows which are frequently used as harvesting paths and are therefore kept vegetation-free, making them particularly prone to surface runoff and erosion (Comte et al., 2012;Guillaume et al., 2016).This explanation is further supported by our sediment-trap data: traps with a low vegetation cover (e.g.HO3 whole year and HO1 as of February 2019) exposed to daily rainfalls exceeding 25 mm d −1 showed higher losses of soil (∼ 50-100 g m −2 , Table 1 and Fig. 3a) than traps (e.g.HO4 whole year and HO1 prior to February 2019) that had less than 50 % of bare soil at similar rainfall intensities.This again highlights the importance of cover crop in oil-palm plantations countering soil erosion (Guillaume et al., 2016(Guillaume et al., , 2015;;Corley and Tinker, 2016;Luke et al., 2019).Furthermore, stacked palm fronds, especially if aligned perpendicularly to the slope, may reduce soil erosion on oil-palm plantations (Corley and Tinker, 2016).Some of the questions to be answered in this study included the extent to which soil erosion reduces the topsoil Si Am pool in oil-palm plantations and whether the lower denhttps://doi.org/10.5194/soil-9-169-2023SOIL, 9, 169-188, 2023 sity of phytoliths compared to mineral soil particles caused proportionally greater losses of Si Am through soil erosion.
To our knowledge, only a few studies exist in which the effect of soil erosion on the topsoil Si Am pool has been addressed.Almost all of them focused on arable soils (Clymans et al., 2015;Unzué-Belmonte et al., 2017;Kraushaar et al., 2021).Clymans et al. (2015) determined mean topsoil Si Am concentrations of 1.76 mg −1 g −1 in arable fields in sloping terrain and temperate climate.This compares well to topsoil Si Am concentrations in interrows from our study (Table 1).
In contrast, Si Am concentrations in eroded soil material differed by a factor of 2 between plots HO1 (1.61 mg g −1 ) towards HO3 (0.88 mg g −1 ) and HO4 (0.66 mg g −1 ) (Table 1 and Appendix Table B3).A possible explanation could be the differing maintenance of the cover crop.High median Si Am concentrations in eroded soil material were determined in previously vegetated traps (e.g.HO1 until the end of January 2019), whereas lower median Si Am concentrations in eroded soil material were measured in traps with less vegetation (HO3 and HO4, whole year, Table 1).We may infer from this observation that the cover crop in plot HO1 maintained higher Si levels in the topsoil through continuous phytolith release from litter.After weeding and keeping the vegetation cover at around 5 %, this phytolith-enriched topsoil was eroded, leading to higher Si Am concentrations in the eroded soil material.In contrast, plots HO3 and HO4 had lower and more dispersed Si Am concentrations in eroded soil material as they lacked an additional Si source.Further, low Si Am concentrations suggest that topsoil with originally high amounts of biogenic Si Am has already been eroded over time, leaving mainly pedogenic Si Am .A greater variability in Si Am concentrations in eroded soil material in plots HO3 and HO4 was probably caused by a slight increase in vegetation cover during the year and, secondarily, by varying daily rainfalls (Fig. 3).Thus, these observations could provide a basis to state that phytoliths are preferably eroded from topsoils.This in turn would assert our hypothesis.Nevertheless, further field experiments and observations are required to confirm this statement.

Spatial topsoil Si M concentration patterns
Topsoil Si M concentrations in well-drained plots were highest under frond piles (∼ 14 µg g −1 ), followed by palm circles (∼ 11 µg g −1 ), and lowest in oil-palm rows and interrows (∼ 6 µg g −1 ).However, these differences were not statistically significant (Fig. 2c, d).Higher topsoil Si M concentrations under frond piles can be explained by the high solubility of Si Am that is released from decaying palm fronds in the form of phytoliths.In addition to Si Am (Barão et al., 2014;Unzué-Belmonte et al., 2017), Si associated with soil organic matter (SOM) also represents a readily soluble Si fraction in topsoils (Alexandre et al., 1997;Georgiadis et al., 2013;von der Lühe et al., 2020).Such readily soluble Si fractions usually contribute most Si to the soil solution (Struyf et al., 2010;de Tombeur et al., 2020).
The plots in the riparian areas showed high topsoil Si M concentrations under the frond piles (∼ 20 µg g −1 ).All other management zones had lower topsoil Si M concentrations in the range of ∼ 11-12 µg g −1 (Fig. 2b).However, this difference was also not statistically significant.In riparian plots, flooding may lead to a redistribution of Si M across the oilpalm plantation, hence explaining similar topsoil Si M concentrations in palm circles, oil-palm rows, and interrows.In riparian areas that are flooded during the rainy season, dissolved Si in stream water (Cornelis et al., 2011;Dürr et al., 2011) may be another source of topsoil Si M alongside Si Am .Therefore, we had hypothesized that Si input from stream water may lead to higher topsoil Si M levels in riparian areas compared to well-drained areas.Indeed, topsoil Si M concentrations under frond piles in riparian plots (∼ 20 µg g −1 ) tended to be higher compared to well-drained plots (∼ 14 µg g −1 ).Likewise, topsoil Si M concentrations in oil-palm rows and interrows in riparian plots (∼ 11-12 µg g −1 ) also tended to be higher compared to well-drained plots (∼ 6 µg g −1 ).However, these differences were not statistically significant, so our hypothesis cannot be fully asserted.

Conclusions and recommended measures
Based on the differing topsoil Si Am concentrations observed in the different management zones, we conclude that current oil-palm management practices cause a distinct spatial topsoil Si Am concentration pattern.Especially the stacking of cut-off palm fronds in long piles and subsequent decomposition promotes Si Am return to soils.Thus, the highest topsoil Si Am concentrations occur below frond piles.Similarly, high concentrations may be found in interrows if additional sources of biogenic Si Am such as Si-accumulating plants (grasses, sedges) are present.Lower topsoil Si Am concentrations in oil-palm rows and unvegetated interrows reflect a lack of Si Am return to soils through plant litter in these management zones.Moreover, pronounced topsoil erosion in unvegetated interrows involves Si Am losses and may therefore cause additional Si Am depletion in this management zone.A dense cover of grasses and mosses in interrows may efficiently reduce erosion and associated Si Am losses.
Topsoil Si M concentrations in the different management zones showed that biogenic Si Am was an important readily available source of Si M .Thus, analogous to topsoil Si Am concentrations, the highest topsoil Si M concentrations also occurred under frond piles.Our hypothesis that regular flooding involves an input of Si dissolved in stream water into the system in riparian areas, partially replenishing the Si M pool, could not be statistically proven in this study.Although topsoil Si M concentrations tended to be higher in riparian areas, the differences between well-drained and riparian plots were not statistically significant.
In conclusion, our findings suggest that erosion could be reduced efficiently, and Si cycling could be maintained within the system if smallholders followed some suggested measures such as (i) maintaining a grass cover in oil-palm rows and interrows, (ii) incorporating oil-palm litter into plantation management, and (iii) preventing soil compaction and surface-crust formation.It would be advisable to raise awareness on topsoil erosion and its potential causes.Furthermore, any logistical efforts and costs involved for implementing these measures (e.g. for ameliorating soil compaction) would have to be feasible.

Figure 1 .
Figure 1.(a) Triangular planting arrangement on smallholder oilpalm plantations in the study area (sketch by Britta Greenshields).Topsoil samples were taken in four distinct management zones: (A) palm circles, (B) oil-palm rows, (C) interrows, and (D) below frond piles.(b) This is a 2 m 2 sediment trap with scarce understorey vegetation installed in pairs in four well-drained plots (sketch by Britta Greenshields followingSinukaban et al., 2000).The downslope funnel-shaped part of the aluminium frame (F) directs surface runoff, together with eroded soil material, into a bucket (B1) that is connected to a second bucket (B2) by a 2 cm wide tube (T).Photos of the sediment traps are shown in the Appendix TableA3.

Figure 2 .
Figure 2. Concentrations of mobile Si (Si M ) and Si in amorphous silica (Si Am ) in topsoils of four different management zones: palm circles (n = 4), oil-palm rows (n = 4), interrows (n = 3), and below frond piles (n = 4) on smallholder oil-palm plantations in two different landscape positions with differing water regimes (well-drained and riparian).Boxes indicate interquartile ranges and whiskers extend 1.5 times the interquartile range below or above the box.If lower-case letters (a, b) differ from one another, this indicates a significant difference between management zones within a water regime (p ≤ 0.05).The Si concentrations were calculated with oneway analysis of variance (ANOVA) and Tukey-Kramer post hoc test.

Figure 3 .
Figure 3. (a) Weekly losses of topsoil and (b) Si in amorphous silica (Si Am ) in eroded topsoil, collected from sediment traps (n = 6) of oil-palm plantations in well-drained areas.

Table 1 .
Topsoil Si Am concentrations in interrow and sediment-trap samples.
a Si Am concentrations for plot HO1 as of February 2019 (after manual weeding).b Si Am concentrations for plot HO3 and HO4 for the whole sampling duration.c Replicates for interrow topsoil samples/replicates for eroded soil samples.d Every interrow on plot HO1 contained stacked frond piles with no sampling possible.NA: not available.

Table 2 .
Annual losses of soil and Si Am through erosion.

Table B1 .
Mean topsoil Si concentrations in different management zones of oil-palm plantations.

Table B2 .
(a) Topsoil Si concentrations in different management zones of well-drained oil-palm plantations.(b) Topsoil Si concentrations in different management zones of riparian oil-palm plantations.

Table B3 .
Weekly loss of eroded soil material and Si Am in eroded soil material from sediment traps under oil-palm plantations (sloping terrain).
NA: not available.

Table B4 .
Mean Si concentrations and statistical analyses on log-transformed data.± standard deviation.Statistics were conducted by one-way ANOVA and Tukey-Kramer post hoc test.Normally distributed data, whereby the homogeneity of variances was asserted.
Mean * Italics suggest that the homogeneity of variances was not asserted.