Changes of Technosol properties and vegetation structure along a chronosequence of dredged sediment deposition in areas with alluvial gold mining in Colombia

Alluvial gold mining leaves a vast amount of dredged sediment that covers the natural soil, destroys riparian ecosystems, and impacts riverbeds and valleys. In Colombia, more than 80,000 ha are covered by dredged sediment. Technosols developed from such deposits undergo strong transformations during the early stages of formation due to the diverse nature of the constituents. The aims of this study are (i) to investigate changes in soil morphological and physicochemical properties at early stages of pedogenesis, (ii) to estimate the variability of size structure distribution of the new plant cover as affected by the time period since establishment and deposit type, and (iii) to determine which soil parameters correlate best with the size structure distribution of vegetation. Soil characterization through morphological and physicochemical analysis was conducted on Technosols developed from gravel and sandy deposits of different ages (0–12 years). A vegetation survey was conducted to measure tree dimensions and identify the most common species. Data were analyzed using non-parametric tests to avoid the loss of information due to data smoothing. The studied Technosols are subject to processes similar to those occurring in natural parent materials such as mineral transformations, changes in redox conditions, organic matter accumulation in the topsoil, and a rapid differentiation of horizons. However, the Technosols in our study sites show an unusually fast development not often observed in natural soils, with drastic changes observed at very early stages of formation. Gravel deposits offer more favorable conditions for plant growth, and marked changes in pedogenic processes can be observed compared with sandy deposits, which are reflected in changes of chemical properties. Factors such as the technology used for mining and the continuous deposition of dredged sediments in a cumulative way during the exploitation period, as well as the disturbances after deposition, result in a high diversity of constituents and high heterogeneity of dredged sediment deposits. Gravel and sand deposits, originated from the same sections of the river and after similar separation processes for ore extraction, undergo divergent pedogenic processes at different rates, which could be explained by a remarkable effect of particle size distribution. For vegetation establishment, the selection of tree species should be based on their capacity for nutrient pumping through deep rooting to sequester carbon and to adapt morphologically to heterogeneity in nutrient availability by growing roots in nutrient-rich zones.

in 2016, an area of 78,939 ha had been affected by alluvial gold mining by 2014 (UNODC 2016). In addition, 24,450 ha of natural forest and secondary vegetation were lost in 2014 due to alluvial gold mining activities. An analysis of the temporal dynamics shows that between 2001 and 2014, an area of 44,746 ha was degraded by alluvial gold mining over the whole country (UNODC 2016). Once the exploitation of alluvial deposits using suction and bucket dredges has taken place, the dredges are transferred to a new working front leaving behind three elements: (i) gravel deposits that consist of rock fragments coming from material extracted by bucket dredges, (ii) sandy deposits that consist of sand from material extracted by suction dredges, and (iii) artificial wetlands that remain after the exploitation of the deep alluvium of the river and surrounding swamps. Deposits created by alluvial gold mining usually have a low macronutrient content and an acidic pH that tend to disrupt soil-forming processes and plant growth (Cooke and Johnson 2002). There is also a discontinuity between the upper and lower parts of the soil profiles due to the superimposition of the dredged sediment over the natural soil (Wahsha and Al-Rshaidat 2014).
The dredged sediment deposits formed after alluvial exploitation are the parent material of soils which can be classified as Technosols given the human influence exerted on them and the fact that their properties and pedogenesis are dominated by their technical origin (IUSS Working Group WRB 2015). Soils that are the result of a combination of materials, such as the materials deposited by dredges, undergo strong transformations during the early stages of formation. Even if the composition of Technosols is in principle different to the composition of natural soils, their pedogenesis is controlled by the effect of aggregation or decarbonisation factors, i.e., soil fauna, microorganisms, plant roots, binding agents, and environmental conditions which may be similar to those in natural soils (Jangorzo et al. 2013). However, technogenic materials often display highly diverse constituents, high spatial variability, and frequently reported more intense and rapid pedogenic processes in the first stages of weathering than in natural soils (Huot et al. 2015), possibly due to imbalances between soil artifacts (IUSS Working Group WRB 2015) and environmental conditions (Séré et al. 2010).
Technogenic materials changed by soil-forming processes may create (i) a high diversity of constituents and the potential for simultaneous occurrence of processes that would have taken place in contrasting climatic conditions, (ii) high spatial variability of the materials and changes in topography resulting from localized processes, anthropogenic influence, or distinct rates of soil development over short distances, and (iii) successive and simultaneous processes occurring at different time scales (Huot et al. 2015). Therefore, it is important to analyze the degree of advanced soil-forming processes in Technosols over periods of time based on the analysis of soil morphological properties and conventional measurements of soil physicochemical properties that are the result of successional changes in the vegetation cover.
One approach to characterize transformations of Technosols typically involves the identification of a chronosequence consisting of sampling and measurements in recently deposited materials and contrasting the results with those of Technosols developed a long time ago. Yet, Technosols are very heterogeneous, and the properties of each area of the deposits depend on the technology used for gold mining or the intensity of weathering processes. Therefore, transformation processes have to be interpreted with care (Uzarowicz et al. 2017). Furthermore, comparing plant community resemblance between different sites or the progressive change over time can provide a sensitive measure of ecologically relevant changes in the environment (Philippi et al. 1998). One approach to achieve this is to determine plant community resemblance based on size structure. According to De Cáceres et al. (2013), in plant communities, the most natural structural variables for wood vegetation are plant height and trunk diameter. Abundance profiles of the plant community for a given area can be calculated with the values of the structural variables for each individual.
The main objective of this study was to analyze the early pedogenic trends and patterns of vegetation establishment in Technosols developed in dredged sediment deposits through the analysis of changes in morphological and physicochemical properties after different time periods since the deposition of dredged sediments, at different depths and for the two types of deposits found in the study site, i.e., gravel and sandy deposits. The specific aims of this study were (i) to investigate changes in soil morphological and physicochemical properties at early stages of pedogenesis, (ii) to estimate the variability of size structure of the new plant cover resulting from the time since dredged sediment disposal and deposit type, and (iii) to determine which soil parameters correlate best with the size structure distribution of the vegetation.

Site description
Both gravel and sandy deposits were established by hydraulic filling. The sediment deposits consisted of bottom sediments dredged from river sections in the neighborhood of the landfill site. The first stage in alluvial gold mining is to take the dredged riverbed material and separate the finer fractions where the gold is found from the larger mineral fraction. In a typical alluvial mining process, physical separation methods such as screen filters and gravity separation are used to separate the gold from the mineral fraction. After separation, the larger mineral fraction is pumped over the deposition sites. Gravel material is deposited by bucket dredges in strip-like patterns that are used to contain the sandy material pumped by suction dredges. At the time of the sampling, the ages of the sites varied between 1 and 12 years since deposition.
To mitigate the environmental impact of alluvial gold mining, frequent attempts have been made to establish vegetation on gravel and sandy deposits left by dredges. The Colombian company Mineros S.A. established farm forestry systems on gravel and sandy deposits created after alluvial gold exploitation with suction and bucket dredges. Farm forestry systems are defined as the growth of woody perennials on the same management unit as used for agricultural crops and animal husbandry, with both ecological and economical interactions between the woody and non-woody components of the system (Nair 1993). Vegetation establishment begins by flattening the topography with bulldozers in order to reduce the slope to below 30% for the subsequent establishment of a vegetation cover. Multipurpose trees are often planted on cropland, around homesteads, as shelterbelts, or for vegetation establishment on less fertile areas. Farmers commercialize the timber, use it as firewood, and often do tree beekeeping. The most common tree species found in the area are Cecropia peltata (Urticaceae), Acacia mangium (Fabaceae), Guazuma u l m i f o l i a ( M a l v a c e a e ) , M i c o n i a m i n u t i f l o r a (Melastomataceae), Sapium glandulosum (Euphorbiaceae), Tectona grandis (Lamiaceae), as well as shrubs such as Eupatorium vitalbae (Asteraceae), Vernonanthura patens (Asteraceae), and Rubus ulmifolius (Rosaceae). The most abundant herbaceous species in the area are Calathea lutea (Marantaceae) and Cyperus ferax (Cyperaceae). The aim of the company is to support agricultural land use by establishing crops in more fertile areas of the deposits while creating productivity on less fertile areas by planting trees. So far, under the Mineros S.A. initiative, more than 700 ha have been restored in this way and more than 40 farmers and their families have benefitted by the entitlement of 14-ha plots.
Based on a detailed survey of the Techosols derived from dredged sediments, four areas were chosen for this study based on the time since deposition of dredged sediments, as well as on the presence of both types of deposits. The study sites are located in northwest Colombia in the gold mining area of El Bagre, Antioquia ( Fig. 1), in the humid tropical forest zone with an average temperature of 28°C and an annual precipitation of 2000-4000 mm. The region has a dry season from November to March and a rainy season from April to October. The topography is mostly flat with the characteristics of a low-lying forest that remains flooded most of the year, thereby preventing productive activities in agriculture, animal husbandry, or forestry. Such areas are the most affected by gold mining.

Soil sampling and analysis
Areas where dredged sediments had been deposited since 2002, 2006, 2010, and 2014 were selected for this study. Six sampling sites were chosen within each of the four areas giving a total of 24 sampling sites with 12 located in gravel deposits and 12 in sandy deposits. Soil samples were collected from September 2015 to February 2016 at depths of 0-20, 20-40, and 40-60 cm. For soil sampling, an area of 10 × 10 m was selected at each sampling site to describe the morphology of the soil profile, to perform field measurements of hydraulic conductivity, and to take samples for physicochemical soil analyses. For soil sampling, an auger with 10-cm diameter was used for sandy deposits and a shovel for gravel deposits. Gravel is defined as any particle with 2-60 mm diameter.
Samples for bulk density were taken as undisturbed soil cores (10 cm diameter, 7 cm length) from each sampled surface layer. Bulk density was estimated with oven-dry cores (105°C for 48 h). To estimate the water retention curve (field capacity, saturated water content, permanent wilting point (%, v/v)), undisturbed cores (4 cm diameter, 6 cm length) were taken from each sampled surface layer. Saturated water content, field capacity (applied pressure 0.33 bar), and permanent wilting point (applied pressure 15 bar) were determined according to Richards and Weaver (1944). The structural stability index was computed following Yoder (1936). Soil texture was analyzed using the Bouyoucos hydrometer method. Hydraulic conductivity was measured in the field using a tension infiltrometer (Eijkelkamp, Netherlands) and a Guelph constant head permeameter (Eijkelkamp, Netherlands). Soil pH was determined by the potentiometric method 1:1 water, total organic carbon (TOC) by Walkey-Black, and total nitrogen (TN) by the Kjeldahl method. Available phosphorus (P) was assessed by a modified Bray-II method. Contents of potassium (K), calcium (Ca), magnesium (Mg), iron (Fe), manganese (Mn), zinc (Zn), copper (Cu), and aluminum (Al) were quantified by atomic absorption spectroscopy after extraction with HCl (37%, m/v). Altitude, longitude, and latitude were recorded using a portable GPS (Garmin eTrex 10, Germany).

Soil morphology
Soil morphology (Figs. 2 and 3) was studied following the guidelines developed by the US Soil Science Division Staff (2017). Soil pits were dug and the different horizons of the soil profiles identified. The morphology of each soil horizon was described according to the field-description handbook of Loaiza et al. (2015) and the FAO Guidelines for soil description (Jahn et al. 2006) for the following properties: presence of roots, root size, particle shape, particle size, color, structure, consistence, and stoniness. For the soil color description, Munsell color charts were used (Munsell 2013), and for the mineralogy analysis in the field, × 10-20 lenses were used.

Vegetation survey
The measurement of the size structure of vegetation was conducted on the 24 sampling sites according to González et al.  (2006). Trees and shrubs were georeferenced within the boundaries of the 10 × 10 m plots, and diameter at breast height (DBH), as well as total height of each living plant was measured. Tree heights were measured with a Vertex IV hypsometer (Haglöf, Sweden).

Data analysis
In this study, data were strongly non-normal and resistant to transformation. Non-parametric tests rely on the assumption that the samples are independent and come from the same distribution and are usually less powerful than the corresponding parametric tests when the normality assumption holds. Nonparametric tests often require modifying the hypothesis, as the tests about the population center are about the median instead of the mean, which does not answer the same question as the corresponding parametric procedure (McArdle and Anderson 2001;Anderson 2001). In spite of this drawback, the use of non-parametric procedures was chosen as the best option for this study, as the sample size was in some cases not large enough to perform parametric analyses with strongly non-normal data. A three-way permutational multivariate analysis of variance (PERMANOVA) was conducted in which soil physicochemical properties were examined as a function of time period since dredged sediment deposition, sampling depth, and deposit type (gravel or sandy deposit). The P values were calculated on 10,000 permutations of the data using the adonis function in the R vegan package (Oksanen et al. 2017). Kruskal-Wallis tests were used to independently compare each soil parameter associated with each deposit type. To visualize the results of the PERMANOVA, multivariate homogeneity of group variances was assessed by means of principal coordinate analysis. Each soil sampling point is depicted with its mean distance to the centroid of each grouping factor, which in this case is deposit type, time period since establishment of vegetation, and sampling depth. These figures were elaborated using the betadisp function in the R vegan package.
A linear discriminant analysis (LDA) was performed (Huberty and Morris 1989) to identify the soil properties that contributed most to differences between gravel and sandy deposits. For this purpose, the values of soil properties for both types of deposits were standardized and weights were assigned to each soil property in terms of their discriminating power between gravel and sand. The soil properties with large weights were considered as those that contributed the most to differences between deposit types. This analysis was performed using the lda function of the R package MASS (Venables and Ripley 2002). A graphical representation of the two-dimensional partition done by LDA analysis was elaborated using the partimat function of the R package klaR (Weihs et al. 2005). Abundance profiles of living plants were calculated for each plot to analyze the distribution of plant sizes. Based on the method developed by Clarke and Ainsworth (1993), an analysis was performed to identify the combination of soil properties that best explained the patterns observed in the size distribution of the plant community. For this, the bioenv function of the R package vegan was used to identify the most suitable subset of soil-physicochemical properties that best correlated with plant community similarity matrices constructed based on size distribution measurements (Oksanen et al. 2017). The most suitable subset of soil properties found after the bioenv analysis was further subjected to a permutation test using the adonis function in R to determine significance.

Soil morphology
The morphological description of the Technosols for both gravel and sandy deposits (Tables 1 and 2) shows a superimposition of sub-horizontal layers, in some cases with wavy boundaries within layers. In the gravel deposits, fine roots were rare in the C horizons. However, frequent fine and medium roots could be observed in O horizons (0-2 cm). Formation of O horizons, dominated by organic material, was observed for gravel deposits > 4 years after dredged sediment deposition, and for sandy deposits with > 8 years after deposition. Different kinds of material were recorded in both gravel and sandy deposits > 4 years after deposition, with changes in particle-size distribution and color. The upper layers of the gravel deposits exhibited a light brown color (10 YR4/3), and rounded and sub-rounded gravel, with a few fine roots in the younger gravel deposits, and abundant fine and medium roots in the gravel deposits with > 8 years after deposition. Deeper in the profile, layers with various shades of yellow could be identified (5Y6/4, 5Y7/6). Coarse Non-revegetated 4 years 8 years 12 years Fig. 3 Soil profile of sandy deposits in non-revegetated areas and 4, 8, and 12 years since dredged sediment deposition roots were noted in the gravel deposits 12 years after deposition. The upper layers of the sandy deposits were dark yellow and brown (2.5Y 5/4, 10 YR 4/3), with fine and medium roots in areas with > 4 years after deposition. Abundant fragments of angular and sub-rounded quartz, as well as abundant coarse pores, occurred in the upper layers of the sandy deposits. Redoximorphic features (Lindbo et al. 2010) were observed in the profile of sandy deposits > 4 years after dredged sediment deposition, where the water table was found at a depth of 58 cm. Accumulation of iron oxide in the root channels was found in sandy deposits with > 4 years after deposition, indicating temporary changes in the redox potential. In addition, the presence of sapric material was found in the sandy deposit 8 years after deposition. Moderately decomposed organic matter occurred in gravel deposits with > 4 years after deposition, and biological activity of mesofauna was found in the profile of sandy deposits  Other features: sapric material; slightly decomposed organic matter; structure and tissue of litter are conserved > 4 years after deposition. As expected, in both sand and gravel deposits, the abundance and size of roots increased from younger to older areas.

Physicochemical properties of Technosols
Descriptive statistics of the physical properties of the deposits are shown in Table 3. The gravel fraction is usually considered as having no effect on the chemical and biological functioning of the soil and was therefore removed before the laboratory analysis. The gravel content represents 8.4 ± 3.7% (SD) of the total weight of the samples taken in gravel deposits. The percentage of sand observed in the fraction with < 2 mm particle size of gravel deposits shows a moderate to high content of sand in non-revegetated areas and very high content in areas > 4 years after dredged sediment deposition. A moderate silt content was found in non-revegetated gravel deposits and low to moderate content in gravel deposits > 4 years after deposition. The clay content was highest in non-revegetated gravel deposits and lowest in non-revegetated sandy deposits. A very high sand content was observed in non-revegetated sandy deposits but was lower in sandy deposits > 4 years after deposition. Low silt contents were observed in nonrevegetated sandy deposits and low to moderate silt contents in sandy deposits 4 years after deposition, but these were lower in areas > 8 years after deposition.
Stability index, bulk density, and hydraulic properties were measured for undisturbed samples without removing the gravel. A poor structural stability was observed in both gravel and sandy deposits. The highest stability indices were observed in non-revegetated sandy deposits. In general, moderate to low bulk density values were observed in both deposit types. Gravel and sandy deposits exhibited low water retention capacity and water contents at wilting point, characteristic of soils with a coarse texture. The water retention properties of the sand and gravel deposits showed similar values for both deposit types and all stages of vegetation establishment. The saturated hydraulic conductivity was higher in gravel deposits than in sandy deposits. A high variability of saturated hydraulic conductivity was observed for both gravel and sandy deposits for all time periods after dredged sediment deposition, ranging from very low to moderate values in both deposit types.
Soil acidity (pH) and selected soil chemical properties are shown in Table 4 for both gravel and sandy deposits. Following up on the significant differences in chemical properties of the gravel deposits after different time periods since dredged sediment deposition, Kruskal-Wallis tests and Kruskal-Wallis multiple comparisons were conducted for each individual chemical property of the gravel deposits. Results show significant changes in pH, P, K, Ca, Cu, and Na. Gravel deposits > 4 years after deposition have significantly lower pH and lower P, K, Ca, Cu, and Na values compared with gravel deposits that were still without vegetation cover. The pH values for gravel deposits range from nearly neutral in non-revegetated areas to strongly acidic in areas 12 years after deposition. The TN content of gravel deposits is 20% higher in areas 12 years after deposition compared with non-revegetated areas. The content of TOC is 23% higher for gravel deposits 12 years after deposition compared with newly established areas. The content of exchangeable cations is considerably lower in gravel deposits revegetated for > 4 years compared with non-revegetated areas. Fe and Zn content is higher in gravel deposits with > 4 years after deposition. The Cu content is lower in gravel deposits > 4 years after deposition.
The results of the PERMANOVA did not show significant differences in the chemical properties of sandy deposits after different time periods since dredged sediment deposition. Descriptive statistics (Table 3) of individual properties show that pH values are strongly acidic in most cases for the sandy deposits in the four studied areas and do not exhibit high variability in the three studied sampling depths. No accumulation of TOC and TN can yet be observed > 12 years after deposition. Exchangeable cations and Fe content are higher in sandy deposits > 4 years after deposition. Similar Zn values are found for the sandy deposits of the nonrevegetated and revegetated areas. The Cu contents are similar in sandy deposits at all times after sediment deposition.
Linear discriminant analysis (LDA) shows that Ca and P content contributes the most to differentiating gravel and sandy deposits, and that Zn, Al, Cu, and Fe are also of considerable importance (Fig. 5). Differences in Ca content explain 26.7% of the variability between gravel and sandy deposits, being the parameter that has the highest differentiation capacity between both deposit types. Phosphorus and Zn content each explain 12% of the variability between gravel and sandy deposits, followed by contents of Al (11%), Cu (10.7%), and Fe (10.4%). Two-dimensional partitions (Fig. 5) of the soil properties with highest discriminant capacity (Ca, P, Zn, Al) show that the lowest error rates for classification between gravel and sandy deposits based on soil properties are obtained when Ca and Zn are used for partition, followed by Ca and Al. This indicates that the most relevant chemical properties that differentiate gravel and sandy deposits are the Ca, Zn, and Al content.

Development of the vegetation cover
Height distribution of trees and shrubs was determined in different years of vegetation establishment for both deposit types (Fig. 6). Profiles show that smaller trees < 5 m in height predominated on the plots > 4 years after dredged sediment deposition. A few taller trees (> 10 m height) were found < 4 years after deposition and were more abundant in areas > 8 years after deposition. The height profiles also reveal that the number of living trees in each height range was higher for gravel than for sandy deposits.
For shrubs, the height distribution reveals similar numbers of individuals among the different height ranges for both deposit types > 4 years after dredged sediment deposition. However, gravel deposits > 8 years had higher numbers in all height ranges compared with younger deposits. Analysis of similarity based on the Bray-Curtis index (Clarke 1993) shows that the time period since deposition had a greater effect (R = 0.27, p = 0.001) on changes in size distribution (height) of trees than the deposit type (R = 0.04, p = 0.001). However, for height distribution of shrubs, it was observed that time period since establishment had a similar effect (R = 0.19, p = 0.001) to that of deposit type (R = 0.19, p = 0.001).
Analysis of the density of trees and shrubs in terms of living individuals per square meter shows significant changes in the density of trees, with significantly higher values (p = 0.01) for areas > 4 years after dredged DF, degrees of freedom *Significant differences sediment deposition. No significant differences are observed for shrubs among the different time periods since deposition or between gravel and sandy deposits. Linking the size structure of the plant community to the variations of the chemical properties of the Technosol shows that pH (p = 0.034), P (p = 0.023), and Mg (p = 0.039) had the strongest influence (Spearman correlation = 0. 4365) on the height distribution of the plant community compared with the other chemical properties, as their variability follows similar patterns to those of height distribution of plant community.

Discussion
The Technosols developed from dredged sediment deposits are originally highly diverse and heterogeneous because of the fact that the materials are deposited during cumulative events across time and the disturbances after deposition such as mixing and flattening. This heterogeneity can increase with time, as the nature of the parent materials, their particle size, and soil porosity affect pedogenesis, possibly at distinct rates at smaller scales especially during early stages of formation (Huot et al. 2015). Furthermore, local heterogeneities may change the hydrodynamic of the soil and constrain the development of organisms, which can result in constraints for plant colonization and uneven vegetation coverage (Huot et al. 2015). Despite their young age, Technosols developed from dredged sediment deposits produced by alluvial gold mining are subject to pedogenic processes under the influence of climatic conditions and the development of a vegetation cover. The Technosols undergo processes that are similar to those occurring in natural materials such as mineral transformations, changes in redox conditions, organic matter accumulation at the surface, and a relatively rapid differentiation of horizons. The studied Technosols show an unusually high rate of development not often observed in natural soils, as drastic changes can be observed at very early stages of formation. Just like in natural soils, the colonization by organisms is the main driver of development of soils formed from technogenic materials (Jangorzo et al. 2015). Based on the analysis of size distribution of tree and shrub communities, gravel deposits seem to offer more favorable conditions for the growth of woody plant species than sandy deposits. Therefore, gravel deposits are likely to have a higher rate of pedogenic processes compared with sandy deposits. The evolution of gravel deposits could be dominated by the disequilibrium between the dredged sediments that form the deposits and the environment, which could generate intense changes during initial stages of pedogenesis. In addition, the colonization by organisms in gravel deposits could start a new phase of development characterized by processes driven by biological activity (Huot et al. 2015). Sandy deposits are likely to have a polycyclic evolution at slower rates, dominated by an acidification process unless soil amendments (i.e., lime, ash) are applied, and will prevent vegetation development and slow down the biological pedogenic processes even further.
In the studied Technosols, there is clear evidence of the rapid development of new horizons, visually distinguished according to color and root density in both types of deposits. The formation of O horizons in soils > 4 years after dredged sediment deposition evidences the development of an increasingly abundant plant cover that consequently increases the amounts of TOC in topsoil with increasing age (Uzarowicz et al. 2017). The presence of sapric material, as well as the biological activity > 4 years after deposition, suggests that porosity and aggregation processes Fig. 4 Principal coordinate analysis of soil properties with grouping factors: a years since dredged sediment deposition, b deposit type, and c sampling depth. Dashed lines represent mean distance to centroid of each set of soil properties representing a sampling point to each level of the grouping factor. First and second axes explain 84.4 and 12.7% of the total variability, respectively increased with time under the influence of biological factors. Therefore, the influence of biological factors on pedogenic development of the studied Technosols can be similar to that of natural soils (Jangorzo et al. 2015). However, the lack of significant differences in soil properties at different sampling depths can be explained by the mixing pattern of the material, which initially creates a homogeneous profile. Nevertheless, it is possible that a deeper sampling can reveal differences, as the subsurface is expected to be less influenced by soil forming processes (Huot et al. 2015).
No significant differences were found for changes in soil physical properties of both types of deposits over the different time periods since dredged sediment deposition. The extreme values observed for certain soil chemical properties may limit the power of the chronosequence approach used for this study to elucidate significant trends on Technosol developments. In addition, the lack of significant changes in soil-physical properties could be explained by the relatively short time of establishment of the plots, which could indicate that the time period was too short for environmental conditions and climate to have significantly influenced the physical properties.
The analysis of bulk density in the context of vegetation establishment is crucial given that high bulk densities might limit rooting depth. Severely compacted soils (bulk density > 1.7 g/ cm 3 ) cannot hold enough plant-available water to sustain plants during dry seasons (Maiti and Ghose 2005). In addition, it is expected that while reclaiming mining land, the soil becomes compacted with the subsequent overall decrease in groundwater recharge through precipitation (Ahirwal and Maiti 2016). For the study site, the bulk density values cannot be considered as those of severely compacted soils, although the gravel deposits had higher bulk densities than the sandy deposits (Table 3). Furthermore, the field capacity for both deposits was below 20%. On the same lines, very rapid infiltration was observed for both types of deposits, as well as high saturated hydraulic conductivities. Even though gravel is considered an inert Characters in red represent sets of soil samples that correspond to G wrongly classified as sand deposits and sets of soil samples that corresponds to S wrongly classified as gravel component of a soil, it has a diluting effect on the amount of stored water and nutrients, affects infiltration and runoff rates, and introduces heterogeneity into soils (Rytter 2012).
A slight decrease in pH of the topsoil is a frequently observed process in Technosols forming from dredged sediment deposits (Huot et al. 2015). A moderately low pH was observed in both gravel and sandy deposits, while extremely low pH values that might cause metal toxicity to plants were not registered. The pH values of the gravel deposits were found to vary from neutral for samples corresponding to the non-revegetated area to low values for samples > 12 years after deposition. The change in the pH can be explained by oxidation of the minerals and root growth. Roots typically change the pH of the rhizosphere by imbalances in uptake of cations and anions and by release of organic acids and chelating agents that alter the availability of some nutrients (Darrah et al. 2006). The predominance of an acid reaction has also been documented for technogenic surface formations of dredged dumping grounds in areas of placer gold mining in Russia (Shlyakhov and Osipov 2004).
During mining and reclamation activities, TOC is lost in various ways, i.e., through erosion of soil by water and wind, reduction of inputs due to lack of vegetation, heavy washing of soil by dredges. However, TOC plays a crucial role in vegetation growth and ecosystem development through an improved water-holding capacity, nutrient availability, and soil structure (Huot et al. 2015). For our study site, the levels of TOC were remarkably high. In addition, both types of deposits exhibited higher values of TOC 12 years after dredged sediment deposition. The high values might be attributable to the fact that the deposits were former bottoms of sedimentation basins (Shlyakhov and Osipov 2004). In addition, the higher values found for older deposits might be attributable to vegetation development in reclaimed sites, higher biomass production, and litter accumulation with time (Tetteh et al. 2015). Fig. 6 Abundance profiles of a trees and b shrubs calculated for different years since dredged sediment deposition (0, 4, 8, and 12 years) for gravel and sandy deposits. Total tree/shrub height was used as a structural variable to calculate the profiles Phosphorus levels in gravel deposits > 4 years after dredged sediment deposition are lower than in gravel deposits with no vegetation cover. Previous studies conducted by Castellanos and Leon (2010) in the same research area show that P was the nutrient that limited the productivity of A. mangium used to reforest sand and gravel deposits (Castellanos and Leon 2010). Therefore, as a consequence of vegetation reestablishment, the P concentration is reduced.
In dredged sediment deposits, metals and nutrients can enter the O and A horizons through dry and wet atmospheric depositions, stem flow, litterfall, and direct upward migration from the C horizon (Vandecasteele et al. 2009). Lindsay and Norvell (1978) mention values rated as highly sufficient for ecologically sustainable reclamation as being 4.5 mg/kg Fe, 1.0 mg/kg Mn, 1 mg/kg Zn, and 0.4 mg/kg Cu. In 97% of the samples from the gravel deposits, Zn exceeded the concentration recommended by the above authors, while in contrast only 6% of the sandy deposit samples had Zn concentrations above the threshold. All samples taken from the gravel and sandy deposits showed levels of Fe, Mn, and Cu that exceeded the threshold recommended by the authors. Critical levels of Mn that can begin to affect the yield of sensitive species are > 65 mg/kg Mn (Hazelton and Murphy 2007). For the gravel deposits, 52% of the samples had Mn contents above the critical level and 46% of the sandy deposits also showed Mn values above the critical level.
The analysis of vegetation development through the calculation of size distribution profiles of the trees and shrubs in four different areas with different time periods since revegetation provides evidence of a progressive vegetation growth. In the specific case of Technosols, a progressive plant community succession is only achievable through a favorable combination of soil type, soil physical regime, and type of initial plant community (Ciarkowska et al. 2016). The fact that trees and shrubs were taller in the gravel deposits suggests that the plant roots managed to proliferate into the nutrient-rich zones of these deposits. This can be explained by the high nutrient content of the < 2-mm fraction of the gravel deposits. Plant growth, in spite of the observed uneven distribution of the < 2mm nutrient-rich fraction within the deposit layers, suggests that plants can respond morphologically to heterogeneity in nutrient availability of these deposits by preferentially growing roots in nutrient-rich zones or by shifting growth allocation to or from roots (Robinson 2005). Therefore, changes in root architecture can be mediating adaptation of the plants to these deposits (López-Bucio et al. 2003).
The lack of soil physicochemical properties required for plant growth in some areas of the Technosols can be compensated for through the application of amendments such as green manure or compost (Vergnes et al. 2017). In acidic areas with low concentrations of exchangeable cations, liming or especially wood ash application would be a suitable strategy due to their neutralizing capacity and availability in most households in the area (Neina et al. 2017). However, it is difficult to overcome limiting factors such as soil texture, stoniness or poor structure. Limiting factors for vegetation establishment in gravel and sandy deposits are the lower P contents and exchangeable cations observed in the older deposits, and the high Fe, Mn, Zn, and Al contents, low water holding capacity and poor structural stability. Therefore, a careful selection of the plant species to be planted is required. For example, the use of legumes can increase soil-N content, increase TOC, and also improve soil porosity and structure (Li et al. 2016). In addition, grasses characterized by fast growth can also help to stabilize the deposits and cover the land, and therefore can be used as a suitable initial cover for the barren deposits (Fullen et al. 2006). Vegetation succession will change the dominant plant cover of recently formed deposits from grasses to trees and shrubs, which will likely improve the structure of the dredged sediment deposits through root binding (Santini and Fey 2016). Thus, exclusively planting fastgrowing trees may not be as beneficial as mixed stands (Neina et al. 2017).

Conclusions
The Technosols developed from dredged sediment deposits undergo processes similar to those occurring in natural materials such as mineral transformations, changes in redox conditions, organic matter accumulation at the surface, and a relatively rapid differentiation of horizons. In addition, the studied Technosols show an unusually high rate of pedogenic processes with strong changes at very early stages of formation. These changes can be explained by the technology used for mining, the deposition of material during cumulative events across time, and disturbances after deposition that lead to very heterogeneous technogenic materials. Therefore, the extreme values observed for some soil chemical properties may limit the power of the chronosequence approach used for this study to elucidate significant trends on Technosol developments.
The particle size distribution has a significant effect on the pedogenesis of the studied Technosols, as gravel and sand deposits originated from the same sections of the river and subject to similar processes of removal of finer fractions could undergo pedogenic processes at different rates and with distinct characteristics. The pedogenic processes in gravel deposits occur faster than in sandy deposits. Based on the analysis of size distribution of tree and shrub communities, gravel deposits offer more favorable conditions for plant growth than sandy deposits. Sandy deposits are likely to have a polycyclic evolution at slower rates and may acidify, preventing vegetation development and slowing down pedogenic processes resulting from biological factors. Biological activity of mesofauna observed in sandy deposits suggests that the colonization by organisms could drive a new phase of development.
The farmers living in the study area face the challenge of identifying suitable areas and management strategies for successful cropping. Limiting factors for cropping in gravel and sandy deposits are the low P contents and exchangeable cations observed in the older deposits, and the high Fe, Mn, Zn, and Al contents, low water holding capacity, and poor structural stability. Therefore, a careful selection of the species to be planted is required. For crop production, gravel as well as sandy deposit areas can be used, but organic amendments, such as organic matter or lime would be required for satisfactory establishment and growth of cultivated plants. Their establishment and growth can be supported through an initial coverage with natural topsoil. Such an intervention would pay off most probably only with high-value cash crops. The fact that sand and gravel deposits have remarkably different properties implies that for crop establishment they should be managed differently. Consequently, surveys are required to characterize the chemical and physical properties of the deposits to develop the best suited management strategies for agricultural use and reforestation. Furthermore, before the cultivation of food crops the levels of toxic elements should be assessed to ensure that crops growing in deposit areas do not pose a threat to human or livestock health.