Identifying Mechanisms for Successful Ecological Restoration with Salvaged Topsoil in Coastal Sage Scrub Communities

Although aboveground metrics remain the standard, restoring functional ecosystems should promote both aboveground and belowground biotic communities. Restoration using salvaged soil—removal and translocation of topsoil from areas planned for development, with subsequent deposition at degraded sites—is an alternative to traditional methods. Salvaged soil contains both seed and spore banks, which may holistically augment restoration. Salvaged soil methods may reduce non-native germination by burying non-native seeds, increase native diversity by adding native seeds, or transfer soil microbiomes, including arbuscular mycorrhizal fungi (AMF), to recipient sites. We transferred soil to three degraded recipient sites and monitored soil microbes, using flow cytometry and molecular analyses, and characterized the plant community composition. Our findings suggest that salvaged soil at depths ≥5 cm reduced non-native grass cover and increased native plant density and species richness. Bacterial abundance at recipient sites were statistically equivalent to donor sites in abundance. Overall, topsoil additions affected AMF alpha diversity and community composition and increased rhizophilic AMF richness. Because salvaged soil restoration combines multiple soil components, including native plant and microbial propagules, it may promote both aboveground and belowground qualities of the donor site, when applying this method for restoring invaded and degraded ecosystems.


Introduction
Ecological restoration aims to repair damage caused to degraded ecosystems, which typically are threatened by reduced biological diversity [1]. Restoration goals, therefore, include promoting diverse suites of native organisms, as well as these species' interactions with biotic communities and subsequently deposited at one or more degraded recipient sites. Transporting salvaged topsoil has the advantage of adding seed, potentially reintroducing a diverse native seed bank (which is crucial to the success of a restoration project [26]) to a degraded area, authentically replicating the collective plant community of a cohesive native habitat. The addition of a layer of salvaged topsoil to a degraded site may also bury the existing non-native seeds, reducing competition for natives that germinate from the salvaged seed bank, as well as potentially reducing maintenance needs. This method may also transfer the beneficial effects of the soil microbial community, including mutualistic AMF and soil bacteria, facilitating the development of a native plant and viable microbial community. Land managers have adapted this method from the mining industry, where it has proven to be more effective in mitigating ecological damage than restoration conducted on sites where healthy topsoil is not added [27][28][29]. Although the use of this technique is widespread in the mining industry, its application for non-mining restoration work (although successful [7,30,31]) is understudied. It is still unclear exactly what factors may contribute to the success of salvaged soil restoration.
This study aimed to identify the contributing mechanism(s) to successful restoration with salvaged soil. Specifically, we manipulated the depth at which salvaged topsoil was deposited at multiple degraded sites, to evaluate how salvaged topsoil helps restore systems. We asked the following questions: (1) Does salvaged topsoil reduce competition with non-natives via burial of seeds of non-native, invasive plants? (2) Does salvaged topsoil increase native diversity with the addition of a native seed bank from the donor site? and (3) Does salvaged topsoil result in successful transfer of a cohesive soil microbiome from the donor site? All three of these mechanisms have been cited as possible reasons for the success of restoration with salvaged topsoil, and practitioners recommend different amounts of topsoil to be deposited based on their understanding of the site-specific context, as the relative importance of these factors [32][33][34][35]. If burial of seeds of non-native, invasive species is important, then more soil would need to be deposited; although, if deposition of native seeds is more important, then it might require less soil, depending on the density of native seeds at the donor site. Finally, if propagules from the donor microbiome drive the success of salvaged topsoil, then perhaps only a light dusting of soil is required for the successful transfer of the donor soil microbiome. It is also possible that success may vary depending on conditions at the recipient sites or that some combination of non-native seed burial, native seeds deposition, and microbiome transfer may contribute to the success of restoration with salvaged topsoil.

Study System
This study was conducted within the Nature Reserves of Orange County, in Southern California. The donor site consisted of a diverse coastal sage scrub community that was slated for removal for the development of housing. Coastal sage scrub is a diverse plant community dominated by drought-deciduous shrubs found in Coastal and Baja California [36]. Home to at least 60 plants and 30 animal taxa that are considered rare, this community has seen continuous decline due to urbanization, frequent fires, drought, and atmospheric nitrogen deposition, and is a common target of restoration projects in Southern California [37,38]. Degraded sites once occupied by these coastal sage scrub plants (such as the three recipient sites utilized in this study) are typically dominated by non-native invasive annual species, such as grasses from Bromus or Avena genera, as well as mustards such as Brassica nigra and Hirschfeldia incana.
Characterized by a Mediterranean climate, these sites experience hot, dry summers and mild winters. Rainfall tends to be variable; below average (16.

Site Selection
The Orchard Hills donor site was mapped in August of 2015, prior to soil salvage operations, such that areas with high native shrub diversity were visually selected as potentially high-quality donor sections, while areas with high non-native cover were mapped as sections to be avoided. Three and a half hectares out of approximately 16 hectares were deemed high quality based on these criteria, and they were therefore selected as appropriate to serve as donor material. Soil was salvaged in December of 2015. After existing vegetation was removed, bulldozers were used to scrape soil from approximately the top 15 cm of the site. Scraping depth was visually approximated, and therefore likely spans a range of 10-20 cm of depth. Soil was first piled up at the base of the donor site and subsequently moved into dump trucks, to allow for its transfer to the recipient sites.
Three recipient sites were selected based on their high level of degradation, as well as accessibility. Recipient sites are in relative proximity to the donor site, with a maximum distance of 6.21 km (Figure 1, Table S1). Sites were prepared for deposition of salvaged soil by hand-weeding existing non-native vegetation, as well as herbicide treatment, using low-dose glyphosate targeted at large non-native clusters.

Experimental design
At each of the three recipient sites, three blocks of five experimental plots (3 m × 3 m) were established ( Figure S1). One block represented all deposition-depth treatment combinations, which were randomly assigned to each plot: (a) control (no salvaged soil); (b) light dusting of <0.5 cm thick salvaged soil; (c) 5 cm salvaged soil deposition; (d) 10 cm salvaged deposition; and (e) 15 cm salvaged deposition. A dusting of cover was intended to provide an inoculation of native microorganisms rather than act as a seed-dispersal method ( Figure 2).
Soil deposition occurred at recipient sites in December 2015, in proportion to their total area (Table S1). Soil was initially deposited by dump trucks at the recipient site, in numerous piles, which were evenly distributed throughout the site. These piles were then spread at a uniform 10 cm thickness across the entire recipient site, during which the control and dusted experimental plots were kept clear by using tarps. The even spread of soil was visually monitored. Tarps were subsequently removed to allow for dusting treatment plots to be prepared. The 5 and 15 cm plots were achieved through the transfer of the top 5 cm of soil from the 5 cm plot to the 15 cm plot, and depths were confirmed via measurement by ruler.
Excluding the experimental plots, the recipient sites were weeded regularly via a combination of hand weeding and herbicide. Experimental plots received no explicit maintenance treatment, so as to allow for data collection on non-native species.

Donor Site Seed-Bank Characterization
To allow for analysis of the depth of the donor seed bank, stratified soil samples were collected from six locations corresponding with soil salvage. Samples were taken at the following depths: (a) 0-5 cm, (b) 5-10 cm, (c) 10-15 cm, and (d) 15-20 cm, using a 6.3 cm in diameter bulb-planting tool. Soil samples were separated into organic and inorganic layers through a water floatation method, in which floating organic material is continuously poured off into a vacuum filtration system, until no more is present in the sample [40]. Organic layers were subsequently separated by size, using 2000-, 500-, and 250-micron soil sieves. Sieved layers were analyzed under a dissecting microscope, where seeds were isolated from remnant organic material and counted.

Figure 1.
Map of study sites, including the donor site, Orchard Hills, from which soil was collected, and the three recipient sites, West Loma, Hick's Haul, and Portola, where salvaged soil was deposited.

Experimental design
At each of the three recipient sites, three blocks of five experimental plots (3 m × 3 m) were established ( Figure S1). One block represented all deposition-depth treatment combinations, which were randomly assigned to each plot: (a) control (no salvaged soil); (b) light dusting of <0.5 cm thick salvaged soil; (c) 5 cm salvaged soil deposition; (d) 10 cm salvaged deposition; and (e) 15 cm salvaged deposition. A dusting of cover was intended to provide an inoculation of native microorganisms rather than act as a seed-dispersal method ( Figure 2).

Plant Communities Data Collection
The percent cover of all plants in each study plot was estimated by using a point intercept method in February of 2016. The center of the plot was marked with a flag, and a 1 × 1 m quad pre-marked with 12 intersect points was laid in one quarter of the experimental plot (with the corner of the quad anchored at the pin flag). A pin flag dowel was dropped at each intersect within the frame, and each species touching the dowel was recorded. This procedure was repeated for each of the four quarters of the experimental plot, for a total of 48 total points ( Figure S2). Percent cover was estimated per species Soil deposition occurred at recipient sites in December 2015, in proportion to their total area (Table S1). Soil was initially deposited by dump trucks at the recipient site, in numerous piles, which were evenly distributed throughout the site. These piles were then spread at a uniform 10 cm thickness across the entire recipient site, during which the control and dusted experimental plots were kept clear by using tarps. The even spread of soil was visually monitored. Tarps were subsequently removed to allow for dusting treatment plots to be prepared. The 5 and 15 cm plots were achieved through the transfer of the top 5 cm of soil from the 5 cm plot to the 15 cm plot, and depths were confirmed via measurement by ruler.
Excluding the experimental plots, the recipient sites were weeded regularly via a combination of hand weeding and herbicide. Experimental plots received no explicit maintenance treatment, so as to allow for data collection on non-native species.

Donor Site Seed-Bank Characterization
To allow for analysis of the depth of the donor seed bank, stratified soil samples were collected from six locations corresponding with soil salvage. Samples were taken at the following depths: (a) 0-5 cm, (b) 5-10 cm, (c) 10-15 cm, and (d) 15-20 cm, using a 6.3 cm in diameter bulb-planting tool. Soil samples were separated into organic and inorganic layers through a water floatation method, in which floating organic material is continuously poured off into a vacuum filtration system, until no more is present in the sample [40]. Organic layers were subsequently separated by size, using 2000-, 500-, and 250-micron soil sieves. Sieved layers were analyzed under a dissecting microscope, where seeds were isolated from remnant organic material and counted.

Plant Communities Data Collection
The percent cover of all plants in each study plot was estimated by using a point intercept method in February of 2016. The center of the plot was marked with a flag, and a 1 × 1 m quad premarked with 12 intersect points was laid in one quarter of the experimental plot (with the corner of the quad anchored at the pin flag). A pin flag dowel was dropped at each intersect within the frame, and each species touching the dowel was recorded. This procedure was repeated for each of the four quarters of the experimental plot, for a total of 48 total points ( Figure S2). Percent cover was estimated per species by dividing the number of hits by the total number of intersect points. This sampling method focused on recording data in the center 4 m 2 of each plot, in order to reduce edge effects.
Due to the germination of natives following the initial data collection effort, as well as to determine the effect of treatment over time, additional plant community data were collected in March 2016, and then again in March 2017. Count data were taken in 1 × 1 m quadrats placed in the center of each study plot, avoiding areas with disturbances. For each native species, all individual plants were counted and recorded. The abundance of non-native species was so high that it was not possible Due to the germination of natives following the initial data collection effort, as well as to determine the effect of treatment over time, additional plant community data were collected in March 2016, and then again in March 2017. Count data were taken in 1 × 1 m quadrats placed in the center of each study plot, avoiding areas with disturbances. For each native species, all individual plants were counted and recorded. The abundance of non-native species was so high that it was not possible to count individual plants. Thus, for non-native species, percent cover was estimated within each quadrat.

Bacterial Abundance
Soil samples used to determine bacterial abundance were collected (using a 2.5 cm diameter, 10 cm deep soil core) at six different areas of the donor site in January 2016, corresponding to locations from which soil was salvaged. Soil samples were also collected in triplicate from each treatment plot in January 2016, following initial soil deposition, and again in April 2017, near the end of the growing season, for the second year of the study. Samples were homogenized, intermittently stored at 4 • C and subsequently analyzed via an optimized method described by Khalili et al. (2019) in which soil samples are combined with a detergent, shaken, and layered onto a Nycodenz solution. Flow cytometry was performed on the resulting cell fraction [41].

Arbuscular Mycorrhizal Fungi
To characterize the structure and functional attributes of the AMF community, initial soils were collected during December 2015, in triplicate, from each of six locations, using a 2.5 cm diameter soil core sampler, at a depth of 10 cm, within the selected high-quality hectarage at the donor site, as well as from three locations at each of the three recipient sites. These samples were collected prior to soil scraping at the Orchard Hills donor site but after existing vegetation was removed from both donor and recipient sites. Soils were transported in sterile plastic bags, within a cooler containing dry ice, to the University of California Irvine; these soil samples were subsequently frozen at −20 • C, within 3 h of collection. Prior to AMF analyses, the triplicate samples were composited and homogenized per sampling location, and separated into three aliquots. For measuring mass, one aliquot was dried; for fungal hyphal analyses, we stored another aliquot at 18 • C. Our final aliquot was used for DNA extraction and AMF community analyses; this aliquot was stored at −20 • C. We collected additional soil samples in January 2016, at six different areas of the donor site, and in January 2016 and March 2017, from experimental plots at each of the three recipient sites, to characterize bacterial and fungal abundance and examine how AMF communities changed, following the salvage soil deposition through time, and as a result of the experimental manipulations. These soils were collected in triplicate, as per the initial soil collection, and transported, stored, and processed, as above, for subsequent analyses.
As a metric for soil fungal abundance, or fungal biomass, we extracted fungal hyphae from soil. We measured fungal hyphal length, as per procedures modified from both Sylvia (1992) and Brundrett et al. (1994) [42,43]. We extracted fungal hyphae from each initial soil sample collected in December 2015, and then from post-treatment soils collected in April 2017, close to the end of the growing season, to quantify AMF hyphal length as an index of AMF abundance. Specifically, we extracted hyphae from 4 g of soil (fresh weight). Soils were placed in a beaker containing 1.5 M sodium hexametaphosphate solution. The solution was stirred vigorously on a stir plate for 5 min, as a vortex emerged; soil slurries were transferred to another beaker, which remained stirring for an additional 2 min. Next, we passed the solution through a 0.2 µm nylon filter, vacuumed through fritted glass, which retained fungal hyphae on the filter as it allowed the solution to pass through. Hyphae were stained by immersing the filter in 0.01% (w/v) acid fuchsin dye; we air-dried each filter for at least 15 min, at room temperature, on a glass microscope slide, after which cover slips were affixed, using polyvinyl alcohol/lactic acid/glycerol medium [44]. Slides were dried at 60 • C for a minimum of 24 h. Then we examined each slide at 40 × magnification, using a phase contrast compound microscope (Nikon, Melville, NY, USA). Arbuscular mycorrhizal hyphal lengths were quantified by examining each filter at 100 random locations, for the presence or absence of AMF hyphae. Arbuscular mycorrhizal hyphae were distinguished from other fungal hyphae based on morphology. Generally, AMF hyphae feature angular branching, lack of septa, and irregular walls [45].

Fungal Library Construction and Sequencing
Samples were kept frozen in a −20 • C freezer, at the University of California Irvine, and then transported on dry ice, to the University of California Riverside, for DNA extraction. From each sample, we weighed and subsequently homogenized 0.25 g of soil on a PowerLyzer 24, following the PowerSoil PowerLyzer DNA Extraction Kit protocol (Mo Bio Laboratories, Carlsbad, CA, USA). We quantified DNA extracts with a NanoDrop 2000/2000c UV-Vis spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). Then, we purified our extracts with magnetic beads, using 1:1 sample to beads, with a PEG-bead protocol, as per Rohland and Reich (2012) [46]. Next, we quantified DNA concentrations in our extracts via nanodrop and standardized the concentrations of our DNA templates to 10 ng/µL. We amplified our samples in duplicate by polymerase chain reaction (PCR), for the 18S region, using the Glomeromycotina-specific AML2 and the universal eukaryote WANDA primer set [47,48]. We constructed our libraries by using a dual-indexing approach, as per Berry et al. (2011) [49]. Our first-step PCR reactions (i.e., PCR1) used primers containing universal tails, which were synthesized 5' to the locus specific sequences [50]. Our PCR reactions contained our template DNA with 1 × Phusion HF Buffer (Thermo Fisher Scientific Our thermal cycler began with a heated lid, and then our program was initiated for 2 min at 95 • C, followed by 35 cycles of 30 s at 95 • C, 30 seconds at 55 • C, and 4 min at 60 • C; then, it entered a cooling phase, which was held at 10 • C. We visualized our PCR1 products on a 1% agarose gel, and then subsequently purified the PCR1 products with a magnetic PEG-bead cleanup; purified PCR1 products were eluted in 20 µL Tris-HCl (pH 8.0). To dilute our locus-specific amplicons, we combined 9 µL of Tris-HCl (pH 8.0) with 1.0 µL of our magnetic-bead-purified sample, to generate diluted purified PCR1 product. We used this dilute-purified product as template in our second, indexing PCR reaction (i.e., PCR2), using indexing primers with sequences matching our universal tails at the 3' end; these matched the Illumina MiSeq flow-cell sequences at the 5' end. Our thermal cycler conditions for tailing (PCR2) reactions were similar to the first-round PCR1 reaction, except that our master mix contained 100 nM of each indexing primer, and only 15 total cycles were performed in the thermal cycler for PCR2. Although PCR1 reactions were set up by adding DNA template from each sample into duplicate wells, for PCR2, our reactions were conducted as single reactions, with only one locus-specific amplicon (PCR1 product) per sample. We used agarose gel electrophoresis to check our indexed PCR2 products, and then we purified the products with the magnetic PEG-bead based cleanup. We quantified our purified indexed amplicons in duplicate, using the Quant-iT PicoGreen dsDNA assay kit (Life Technologies, Grand Island, New York, NY, USA), and pooled in equimolar concentrations, to include every sample into our final pooled library. We concentrated our resulting pooled library with the magnetic-PEG-bead-based protocol and conducted qPCR with P5/P7 adapted standards (New England Biolabs Inc., Ipswich, NY, USA) for quantifying the Illumina adapted fraction of our library. We estimated average fragment sizes by using a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA) prior to sequencing, at the UC Riverside Genomics Core. This sequencing was performed on a MiSeq Desktop Sequencing instrument (Illumina Inc, San Diego, CA, USA) in 2 × 300 paired end mode.

Bioinformatics
The fungal sequences were demultiplexed by the MiSeq, using Illumina protocols. We used the forward read for the SSU locus and checked quality with FastQC [51]. We filtered demultiplexed files, using multiple_split_libraries_fastq.py command in QIIME 1.9.1 [52], with Phred scores of 20 (q = 19) as a quality control parameter; we did not allow any low-quality base calls (r = 0), and we only retained reads possessing 95% of initial sequence length following quality truncation (p = 0.95). We used VSEARCH [53] in uchime_denovo mode to screen for chimeras. We used Swarm [54] with a resolution of d4 for picking operational taxonomic units (i.e., OTUs). This d4 resolution for local clustering threshold collapses sequences with fewer than or equal to 4 differences into a single representative OTU, as long as it passed our stringent quality filtering threshold of q20. We used BLAST to assign taxonomy, with an e-value of 0.001 [55], via the MaarjAM database [56]. We truncated the reference databases, prior to analysis, to only include the region of interest, as well as to limit any spurious results. For all analyses, with the exception of alpha diversity, we used the metagenomeSeq package of Bioconductor [57] in R (R Core Team 2017), to normalize our OTUs, using cumulative sum scaling (CSS-normalization), prior to further analyses; CSS normalization attempts to avoid biases associated with marker gene surveys which could result from uneven sequencing depth. Read counts were rescaled against a quantile, which was determined by assessing the count deviations for each sample, as compared to the distribution of counts across all of our other samples [57]. For alpha diversity analyses, instead of using CSS normalization, we used a rarefied OTU table, with sequences rarefied to 7600 reads per sample.
To examine responses of the AMF community (18S) within a fungal functional group framework, we assigned families of Glomeromycotina to AMF functional groups: rhizophilic, edaphophilic, and ancestral, using AMF resource allocation patterns defined in previous studies [20,58]. We did not include sequences assigned to the taxon Geosiphon pyriformis, as reads reportedly identified as G. pyriformis are not considered as rhizophilic, ancestral, or edaphophilic AMF. This functional group framework is based on studies demonstrating that AMF are host-specific, and exhibit diverse resource allocation patterns [59][60][61]. Arbuscular mycorrhizal fungal families with high allocation to extraradical hyphae, including Gigasporaceae and Diversisporaceae, are members of the guild "edaphophilic"; those with high allocation to root colonization (i.e., Glomeraceae, Claroideoglomeraceae, and Paraglomeraceae) are characterized as belonging to the "rhizophilic" guild; families which allocate lower AM fungal hyphae to either root colonization or soil foraging than the aforementioned guilds are known as "ancestral" AMF families, including Ambisporaceae, Archaeosporaceae, and Acaulosporaceae [20]. Previous studies show that fungal taxa from the rhizophilic AMF functional group are important for protecting their host plant from pathogen colonization, while edaphophilic AM fungi putatively improve plant nutrient uptake [62].

Statistical Analyses
All statistical analyses for plant community data, seed bank characterization, and bacterial abundance were performed, using IBM SPSS Statistics Software v.22 (International Business Machines Corporation, Armonk, NY, USA) and SAS 9.4 (SAS Institute, Cary, NC, USA). The Kruskal-Wallis H test was used to test whether the abundance of seeds in the seed bank varied depending on the depth at which the soil was collected from the donor site (1-5 cm, 5-10 cm, 10-15 cm, 15-20 cm). ANOVA was used to determine the effects of salvaged soil treatment (control, dusting, 5, 10, or 15 cm salvaged soil deposition) on the percent cover of non-native forbs in February. The site was used as a random factor. Tukey's post hoc tests were used to detect pairwise differences. Separate repeated measures mixed-model ANOVAs were used to determine whether percent cover, abundance, and species richness of native and non-native plant species in March varied significantly depending on treatment group, recipient site (Hick's Haul, Portola, or West Loma), or the interaction between treatment, site, and year. Block was included as a random factor in the models; Tukey-Kramer tests were conducted, along with the mixed-model ANOVAs. To determine whether bacterial abundance was affected by treatment, sampling time, or an interaction between treatment and sampling time, ANOVA was used; the site was included as a random factor. Tukey's post hoc tests for pairwise differences were applied, as appropriate.
We used the Bray-Curtis index to compute distance matrices and visualized AMF community turnover and beta-diversity, using non-metric multidimensional scaling (NMDS). Bray-Curtis distances were generated from CSS-normalized data before we filtered our sequence-based dataset for functional group assignment. We used the ggplot2 package [63] and the "stat ellipse" function, with 95% confidence intervals, to visualize these NMDS plots in R (R version 3.2.1; R Core Team 2017). We assessed the fit of the data via calculating stress values associated with the NMDS, and we assigned stress values of less than 0.2, as acceptable. Using our rarefied OTU table, we calculated alpha diversity by using the core_diversity_analyses.py function in QIIME 1.9.1 [52] and calculated effective species number, to test for differences in AMF (18S) alpha diversity between sites, by sampling time, site, and topsoil treatment. We used the Shannon Diversity Index and performed linear models, in order to compare diversity across sites and treatment groups (α = 0.05). Furthermore, we examined AMF community composition across time, site, and topsoil thicknesses, by performing multivariate permutational ANOVA analyses (PERMANOVA), using the "adonis" function within the "vegan" package in R (999 permutations) [64], and evaluated differences in the OTU richness of AM fungal functional groups for the 18S locus, using linear models. We performed all statistical analyses in R version 3.4.4 (R version 3.4.4; R Core Team 2018).

Donor Site Seed Bank
A total of 434 total seeds were found in all soil samples collected at the donor site. Seeds were significantly more abundant (p = 0.025) per mL of soil in the topmost layer (0-5 cm) than in all lower layers (5-10 cm, 10-15 cm, and 15-20 cm). The top 5 cm of soil was found to contain 67% of total seeds, with the top 15 cm of soil contributing 94% of total seeds found. Since the top 15 cm of soil was initially collected from the donor site, we can assume that the majority of seeds from the donor site were transferred to the recipient site ( Figure 3).

Vegetation Cover and Species Richness at the Recipient Sites
Point-intercept data collected in February 2016 were dominated by non-native grasses and forbs, as non-native species typically germinate earlier than natives in coastal sage scrub communities [65]. Native forbs comprised, on average, only 6.19% of the total vegetation cover per plot, and they were not affected by salvaged soil treatment (F 4,30 = 1.637, p = 0.181). The percent cover of non-native grasses was significantly greater in control and dusted plots than in treatment plots in which 5, 10, and 15 cm of salvaged topsoil was added (F 4,30 = 24.598, p < 0.001), indicating that the addition of at least 5 cm of salvaged soil significantly reduced the germination of early season non-native grasses. There were no detectable effects of soil treatments on percent cover of non-native forbs (F 4,30 = 2.589, p = 0.118; Figure 4).

Donor Site Seed Bank
A total of 434 total seeds were found in all soil samples collected at the donor site. Seeds were significantly more abundant (p = 0.025) per mL of soil in the topmost layer (0-5 cm) than in all lower layers (5-10 cm, 10-15 cm, and 15-20 cm). The top 5 cm of soil was found to contain 67% of total seeds, with the top 15 cm of soil contributing 94% of total seeds found. Since the top 15 cm of soil was initially collected from the donor site, we can assume that the majority of seeds from the donor site were transferred to the recipient site ( Figure 3).

Vegetation Cover and Species Richness at the Recipient Sites
Point-intercept data collected in February 2016 were dominated by non-native grasses and forbs, as non-native species typically germinate earlier than natives in coastal sage scrub communities [65]. Native forbs comprised, on average, only 6.19% of the total vegetation cover per plot, and they were not affected by salvaged soil treatment (F4,30 = 1.637, p = 0.181). The percent cover of non-native grasses was significantly greater in control and dusted plots than in treatment plots in which 5, 10, and 15 cm of salvaged topsoil was added (F4,30 = 24.598, p < 0.001), indicating that the addition of at least 5 cm of salvaged soil significantly reduced the germination of early season non-native grasses. There were no detectable effects of soil treatments on percent cover of non-native forbs (F4,30 = 2.589, p = 0.118; Figure 4).  March data from 2016 and 2017 yielded various native germinants. The density of native plants was significantly affected by soil salvage treatments, with higher native density observed in plots treated with at least 5 cm of salvaged soil (F4,59 = 21.42, p < 0.001). Effects of site, as well as interactions among both year-by-site, and site-by-treatment were also observed, driven by higher success at West Loma, as compared to other recipient sites. Tukey-Kramer post hoc tests revealed that, overall, mean native density was highest in treatments with 5, 10, or 15 cm of salvaged soil than in control plots or in plots that received a dusting of soil (Figure 5a; Table S2). Although post hoc tests reveal 15 cm of March data from 2016 and 2017 yielded various native germinants. The density of native plants was significantly affected by soil salvage treatments, with higher native density observed in plots treated with at least 5 cm of salvaged soil (F 4,59 = 21.42, p < 0.001). Effects of site, as well as interactions among both year-by-site, and site-by-treatment were also observed, driven by higher success at West Loma, as compared to other recipient sites. Tukey-Kramer post hoc tests revealed that, overall, mean native density was highest in treatments with 5, 10, or 15 cm of salvaged soil than in control plots or in plots that received a dusting of soil (Figure 5a; Table S2). Although post hoc tests reveal 15 cm of salvaged soil to yield higher native density than plots treated with only 5 cm of soil (with 10 cm of soil resulting in intermediate density), results obtained with 5, 10, and 15 cm of salvaged soil are not statistically unique. Mean native species richness was also significantly affected by treatment, increasing with the amount of salvaged soil added (F4,59 = 20.55, p < 0.0001). Effects of year-by-treatment, site, year-bysite, and year-by-treatment were also found. Similar to mean native density, significant site effects were caused by the greater number of native plant species at West Loma, when compared to other sites. Native species richness decreased at all sites, from the first to the second year of the study, (as did non-native species richness; see Figure 6b) According to post hoc tests, highest native species richness values were recorded in treatment plots with 5, 10, and 15 cm of salvaged soil, when compared to those that received no soil or a dusting of soil (Figure 5b; Table S3). Mean native species richness was also significantly affected by treatment, increasing with the amount of salvaged soil added (F 4,59 = 20.55, p < 0.0001). Effects of year-by-treatment, site, year-by-site, and year-by-treatment were also found. Similar to mean native density, significant site effects were caused by the greater number of native plant species at West Loma, when compared to other sites. Native species richness decreased at all sites, from the first to the second year of the study, (as did non-native species richness; see Figure 6b) According to post hoc tests, highest native species richness values were recorded in treatment plots with 5, 10, and 15 cm of salvaged soil, when compared to those that received no soil or a dusting of soil (Figure 5b; Table S3). Plant community data collected in March 2016 and 2017 revealed that non-native cover tended to be lower in plots with more salvaged topsoil, but this result was not significant (F4,59 = 2.19 p = 0.0808). There were, however, significant differences in non-native cover, depending on year, site, and the year-by-site interaction, potentially due to varying rainfall between the two years of the study. Average non-native cover values were lower in all treatments receiving at least 5 cm of salvaged soil than in control plots in 2016. This trend did not continue into the following year of the study (Figure 6a; Table S4).
Non-native species richness was significantly higher in plots treated with 10 cm of salvaged soil than in other treatments (F4,59 = 4.78, p = 0.0021). Tukey-Kramer post hoc tests only differentiate the Plant community data collected in March 2016 and 2017 revealed that non-native cover tended to be lower in plots with more salvaged topsoil, but this result was not significant (F 4,59 = 2.19, p = 0.0808). There were, however, significant differences in non-native cover, depending on year, site, and the year-by-site interaction, potentially due to varying rainfall between the two years of the study. Average non-native cover values were lower in all treatments receiving at least 5 cm of salvaged soil than in control plots in 2016. This trend did not continue into the following year of the study (Figure 6a; Table S4).
Non-native species richness was significantly higher in plots treated with 10 cm of salvaged soil than in other treatments (F 4,59 = 4.78, p = 0.0021). Tukey-Kramer post hoc tests only differentiate the 10 cm salvage soil treatment as significantly different from the control and dusted treatments, with the 10 cm salvaged soil treated plots containing a higher number of unique non-native species. Plots treated with 5 or 15 cm of salvaged soil had an intermediate non-native species richness (Figure 6b; Table S5).

Soil Microbial Content
Although there was no effect on fungal biomass as AM hyphal length (M/g dry wt soil; p > 0.05), mean bacterial abundance of treatment plots varied significantly based on treatment (F 5,76 = 19.37, p < 0.0001), as well as sampling occurrence. Plots that received 5, 10, or 15 cm of salvaged soil were statistically indistinguishable from the donor site in terms of bacterial abundance when the site was first sampled in January of 2016, with the donor site, as well as plots that received at least 5 cm of salvaged soil, having the lowest bacterial abundance. The same results were observed in April 2017, although bacterial abundance was lower in all plots. During the first sampling occurrence in 2016, the plot treated with a dusting of salvaged soil, meant to serve as an inoculum, had significantly higher bacterial abundance than the donor site or any other treatment plot (Figure 7). 10 cm salvage soil treatment as significantly different from the control and dusted treatments, with the 10 cm salvaged soil treated plots containing a higher number of unique non-native species. Plots treated with 5 or 15 cm of salvaged soil had an intermediate non-native species richness (Figure 6b; Table S5).

Soil Microbial Content
Although there was no effect on fungal biomass as AM hyphal length (M/g dry wt soil; p > 0.05), mean bacterial abundance of treatment plots varied significantly based on treatment (F5,76 = 19.37, p < 0.0001), as well as sampling occurrence. Plots that received 5, 10, or 15 cm of salvaged soil were statistically indistinguishable from the donor site in terms of bacterial abundance when the site was first sampled in January of 2016, with the donor site, as well as plots that received at least 5 cm of salvaged soil, having the lowest bacterial abundance. The same results were observed in April 2017, although bacterial abundance was lower in all plots. During the first sampling occurrence in 2016, the plot treated with a dusting of salvaged soil, meant to serve as an inoculum, had significantly higher bacterial abundance than the donor site or any other treatment plot (Figure 7).

Arbuscular Mycorrhizal Fungal Communities
Although the arbuscular mycorrhizal fungal community composition was altered by topsoil treatment (PERMANOVA, p = 0.002), these effects differed by site (p = 0.03) and treatment (p < 0.001; Figure 8). In 2017, we detected a shift in AMF communities at Portola in plots with at least 5 cm of salvaged soil, as compared to both control and dusted treatments. However, AMF communities at both Hick's Haul and West Loma were equivalent across treatments.

Arbuscular Mycorrhizal Fungal Communities
Although the arbuscular mycorrhizal fungal community composition was altered by topsoil treatment (PERMANOVA, p = 0.002), these effects differed by site (p = 0.03) and treatment (p < 0.001; Figure 8). In 2017, we detected a shift in AMF communities at Portola in plots with at least 5 cm of salvaged soil, as compared to both control and dusted treatments. However, AMF communities at both Hick's Haul and West Loma were equivalent across treatments.
Overall, we detected lower AMF alpha-diversity prior to topsoil delivery (2015) than was detected during our 2017 sampling event (p < 0.03). However, topsoil thickness had variable effects on alpha diversity (p = 0.50). For instance, at Hick's Haul, AMF alpha-diversity was greatest in plots containing the thickest layer of topsoil (15 cm; p = 0.04); dusted treatments at Hick's Haul also contained greater alpha-diversity than was detected in control plots (p = 0.05). In contrast, at West Loma, we detected lower AMF alpha-diversity in plots containing 5 cm of topsoil (p = 0.02) than was found in control plots within the same site. Overall, we detected lower AMF alpha-diversity prior to topsoil delivery (2015) than was detected during our 2017 sampling event (p < 0.03). However, topsoil thickness had variable effects on alpha diversity (p = 0.50). For instance, at Hick's Haul, AMF alpha-diversity was greatest in plots containing the thickest layer of topsoil (15 cm; p = 0.04); dusted treatments at Hick's Haul also contained greater alpha-diversity than was detected in control plots (p = 0.05). In contrast, at West Loma, we detected lower AMF alpha-diversity in plots containing 5 cm of topsoil (p = 0.02) than was found in control plots within the same site.
During this initial pre-treatment sampling event in 2015, recipient plots had statistically equivalent rhizophilic AMF richness to the donor site. However, in 2017, after topsoil treatments, we detected an additive effect of salvaged topsoil additions on taxa from the rhizophilic AMF functional group, such that taxa richness of rhizophilic AMF was significantly higher in recipient plots treated with topsoil than these same recipient sites prior to topsoil delivery (p = 0.001; Figure 9). It is worth noting that no effects of topsoil treatments were detected for either edaphophilic or ancestral AMF functional groups. Figure 9. Richness of rhizophilic AMF taxa. "Donor site" refers to OTU (operational taxonomic unit, equivalent to molecular species concept) richness found at the donor site, while "Recipient sites" refers to all other sampling locations at the recipient sites. Rhizophilic AMF richness was measured During this initial pre-treatment sampling event in 2015, recipient plots had statistically equivalent rhizophilic AMF richness to the donor site. However, in 2017, after topsoil treatments, we detected an additive effect of salvaged topsoil additions on taxa from the rhizophilic AMF functional group, such that taxa richness of rhizophilic AMF was significantly higher in recipient plots treated with topsoil than these same recipient sites prior to topsoil delivery (p = 0.001; Figure 9). It is worth noting that no effects of topsoil treatments were detected for either edaphophilic or ancestral AMF functional groups. Overall, we detected lower AMF alpha-diversity prior to topsoil delivery (2015) than was detected during our 2017 sampling event (p < 0.03). However, topsoil thickness had variable effects on alpha diversity (p = 0.50). For instance, at Hick's Haul, AMF alpha-diversity was greatest in plots containing the thickest layer of topsoil (15 cm; p = 0.04); dusted treatments at Hick's Haul also contained greater alpha-diversity than was detected in control plots (p = 0.05). In contrast, at West Loma, we detected lower AMF alpha-diversity in plots containing 5 cm of topsoil (p = 0.02) than was found in control plots within the same site.
During this initial pre-treatment sampling event in 2015, recipient plots had statistically equivalent rhizophilic AMF richness to the donor site. However, in 2017, after topsoil treatments, we detected an additive effect of salvaged topsoil additions on taxa from the rhizophilic AMF functional group, such that taxa richness of rhizophilic AMF was significantly higher in recipient plots treated with topsoil than these same recipient sites prior to topsoil delivery (p = 0.001; Figure 9). It is worth noting that no effects of topsoil treatments were detected for either edaphophilic or ancestral AMF functional groups. Figure 9. Richness of rhizophilic AMF taxa. "Donor site" refers to OTU (operational taxonomic unit, equivalent to molecular species concept) richness found at the donor site, while "Recipient sites" refers to all other sampling locations at the recipient sites. Rhizophilic AMF richness was measured Figure 9. Richness of rhizophilic AMF taxa. "Donor site" refers to OTU (operational taxonomic unit, equivalent to molecular species concept) richness found at the donor site, while "Recipient sites" refers to all other sampling locations at the recipient sites. Rhizophilic AMF richness was measured in the donor site and the recipient site, "pre-treatment", in December 2015 and recipient sites, "post-treatment" were measured in April 2017. Rhizophilic AMF richness was significantly higher in recipient sites, post-treatment, in 2017, than in the recipient sites, pre-treatment, in 2015 (p < 0.001).

Discussion
We found evidence supporting our proposed mechanisms by which salvaged topsoil can result in successful ecological restoration: (1) reduction of competition with non-natives via burial of invasive seeds, (2) the addition of native plants through establishment of a native seed bank, and (3) successful transfer of a cohesive soil microbiome. While other studies have focused on how salvaged topsoil can be a useful technique for mining reclamation, or investigated whether the salvaged soil may be stored prior to use [29,[66][67][68], this is the first study that we are aware of in which the depth of deposition is compared in terms of its influence on both native and non-native plants, as well as the belowground microbial community. Ecologists are increasingly finding evidence that aboveground and belowground interactions are important for structuring community composition [8,69,70]. Here, we demonstrate that such interactions are also important for understanding ecological restoration.

Reduction of Competition with Non-Natives via Burial of Invasive Seeds
The initial decrease in non-native grasses with deposition of salvaged soil, compared to treatments containing only a dusting or no salvaged soil, supports the premise that salvaged soil is important in blocking the germination of the existing non-native, invasive seed bank. There was a similar pattern for non-native forbs (although no significant effects were observed by treatment). This could be attributed not to the failure of the topsoil to block germination of the existing degraded seed bank of the recipient site, but rather to the presence of non-native forb seeds in the salvaged soil itself. Although seeds found in soil collected from the donor site were not able to be identified by species, the plant communities' data indicate the presence of non-native seeds in the salvaged soil, as evidenced by the higher non-native species richness present in plots containing 5, 10, and 15 cm of salvaged soil in March (especially during the first year of the study). Non-native percent cover data collected in March during the same time were not significantly affected by treatment, but the presence of a significantly higher amount of unique non-native species in plots treated with 10 cm of salvaged soil indicates that non-native seed originated from soil from the donor site. A closer look at vegetation data from March 2016 allows us to see that 44% of the total non-native species found in the experimental plots are present only in at least one of the three salvage treatments (5, 10, and 15 cm) and not in the control or dusted plots (including species like Melilotus indicus, Lysimachia arvensis, and Erodium cicutarium); this indicates that these non-native species were transferred to the recipient site, within the salvaged soil (Table S5b).
In contrast, three non-native species were only found in the control and/or dusted treatment plots (Hordeum murinum, Malva parviflora, and Sisymbrium irio) and other specific non-native species' cover was greatly reduced by the addition of at least 5 cm of salvaged soil (e.g., Brassica nigra), which indicates that restoration with salvaged soil does indeed successfully reduce germination of the existing seed bank by burying seeds at the degraded site. Although more than half of the non-native species were present in only one of the two soil sources (in either the recipient site soil (the control and dusted plots) or the donor site soil (5, 10, and 15 cm of salvaged soil addition), some non-native species were found in all treatment plots. (Table S5b) Rather than indicating unsuccessful reduction of germination of seeds existing at the recipient site, these species were most likely originally present in the seed bank of both the recipient site and the donor site and are therefore present in all plots.
We found that a salvaged soil application of at least 5 cm is successfully able to hinder germination of seed existing at the recipient site, but may also result in our restored sites emulating the characteristics of the donor site's reference ecosystem, including its seed composition. As was the case in our study, this may result in the unintentional transplantation of non-native species to the recipient sites.

Addition of Native Plants via Establishment of a Native Seed Bank
Coastal sage scrub communities are known to form soil seed banks [71]. An intact seed bank was found at the donor site, with 94% of total seeds present in the top 15 cm of salvaged soil, indicating that a majority of seeds were salvaged along with the soil and thereby deposited at recipient sites. Our results, that native density and native species richness were greater for plots treated with salvaged soil than those that received only a dusting or no salvaged soil, provide support for the idea that salvaged topsoil works by transferring native seeds to the recipient sites. In March 2016, during the first year of the study, a total of only three native species were present in control and/or dusted plots, while 21 native species were present in plots receiving either 5, 10, or 15 cm of salvaged soil ( Table S5a). The increase in native density and native species richness were stronger during the first year (2016) of the study than during the second year (2017), and effects varied by site-with West Loma experiencing the most success, according to these metrics. The success of West Loma compared to the other two recipient sites may be due to varying abiotic factors (slope, sun explore, etc.) at this site that benefit plant growth. Nevertheless, potential site differences or site-specific characteristics must be taken into consideration when using this method, and further study identifying criteria associated with sites most receptive to this method of restoration would be a particularly valuable avenue for future research.
The overall weakening of the effect of the initial salvage treatment over time, in terms of native species success, could be due to a variety of factors. First, the decrease of native species richness from 2016 to 2017 could be due to unique species present in the salvaged seed bank, which may have germinated but not reproduced successfully, due to the lack of rainfall during the first year of the study. Alternatively, a lack of maintenance inside experimental plots during the first year may have allowed for the non-natives to compete with and hinder native species [65].
In terms of rebuilding the aboveground plant community, restoration with salvaged soil proved successful in our trial, as the degraded recipient seed bank was covered up by donor site soil, whose attributes were subsequently expressed in the aboveground vegetation. This method, however, remains only as successful as the quality of the donor site. The aboveground vegetation is not necessarily reflected in the seed bank, with non-native seed having the potential to dominate the belowground landscape [71]. In a different restoration experiment involving the deposition of salvaged soil, high percentages of non-native cover were unexpectedly recorded [31]. This method is therefore best suited to salvaging from areas with very high native cover, as well as an intact seed bank. Careful aboveground mapping of vegetation, but also sampling of the seed bank, would likely aid in adequately selecting areas of high enough quality to serve as donor sites. Donor sites with some initial native aboveground or belowground presence of non-natives may still potentially serve as adequate donor sites if this restoration technique is coupled with regular non-native removal.

Transference of the Soil Microbiome
We found that benefits of salvaged soil can influence soil biota. Mycorrhizal fungal communities can shift, and soil bacteria can become less abundant above a particular threshold of salvaged soil thickness employed by practitioners. Although salvaged soil influenced bacterial abundance, total and mycorrhizal fungal abundance either remained unchanged or was highly variable across salvaged soil treatments at recipient sites, likely increasing fungal:bacterial abundance in plots treated with at least 5 cm of salvaged soil. Although there are inherent challenges to comparing across microbial biomass metrics and methodologies, fungal:bacterial dominance putatively has implications for biotic communities' response to environmental change and has also been correlated with ecosystem functions such as carbon sequestration and litter decomposition [72]. However, measuring fungal decomposition or AMF nutrient cycling in salvaged soils was outside the scope of this study.
Deploying AMF inoculum in restoration studies often includes adding root fragments, AMF spores, or material containing hyphal fragments to growth media in the greenhouse or seedlings planted in pots within the nursery, prior to out-planting in the field [73] or directly to the field in restoration experiments [16,32]. Salvaged topsoil from donor sites contains AMF spores, along with plant root and AMF hyphal fragments, and may therefore directly manipulate AM community structure via translocating these propagules to recipient sites. Our study evaluated whether transferring components of the soil microbiome would alter the diversity, composition, and functional capacity of biotic communities at recipient sites. The addition of AMF propagules may improve the encounter rate of germinating seedlings and mycorrhizal spores or hyphae, which potentially increases the degree to which these fungi penetrate [74], and exchange resources with, their host plant's roots; this may subsequently promote their host plant's growth [75].
Previous studies show that arbuscular mycorrhizal inoculation improves seedling survival and plant growth within arid lands [76]. Our study, in a semi-arid ecosystem, transferred AMF propagules, along with the intact seed bank, from a donor site to three recipient sites, and we not only detected increased native plant density and richness, but also found greater AMF alpha-diversity at a majority of plots, following topsoil delivery. Previous studies have shown that AMF diversity declines in invaded systems [77] and that edaphic factors could also play a role in hindering the diversity or determining the structure of AMF communities [78]. Although some plants are generalists (i.e., can associate with a variety of AMF taxa), other plant species require more specific associations of only one, or few, AMF taxa [79]. Restoration projects characterized by diverse assemblages of native and endemic AMF, as may be found when salvaging topsoil from local intact reference ecosystems, may ultimately be more effective at promoting the performance of rare native taxa, and therefore at restoring both the structure and the functioning of intact ecosystems. Our study showed an increase in AMF diversity after topsoil delivery, which corresponded to increased propagule pressure of plant seeds and AMF propagules, as well as changes in vegetation communities in our experimental plots. Since AMF are obligate plant symbionts, changes in vegetative cover and plant community structure could feed back to promote the proliferation of a variety AMF taxa via changes in host plant dynamics.
A site's microbiome may change with severe or monodominant invasions [80], as was characteristic of our recipient sites prior to salvaged topsoil deployment, via mechanisms such as increasing pathogen levels, plant-soil feedbacks, or by altering root-associated microbial symbioses [81]. However, ecological restoration techniques can improve the likelihood of successful restoration outcomes by supporting native plant establishment and helping to establish the AMF community in certain ecosystems [82][83][84]. In our study, we found that salvaged topsoil delivery elicited shifts in AMF communities in plots receiving at least 5 cm of salvaged soil. However, at one recipient site (West Loma), 5 cm of salvaged soil corresponded with low AMF diversity. Although light dusting treatments ostensibly added biotic propagules, AMF alpha-diversity was only higher at one site (Hick's Haul) in these dusted plots, than in untreated controls. Given that bacterial abundance was highest in these same dusted treatments, it is likely that fungal-bacterial interactions, such as microbial antagonism, coupled with limited substrate availability, as was likely characteristic of this light dusting treatment, may have resulted in conditions that promoted bacterial growth over plant root mutualists, such as AMF. However, there were more rhizophilic AMF taxa detected overall in recipient sites receiving topsoil than in recipient locations prior to salvaged topsoil delivery. These diverse assemblages of rhizophilic AMF may have contributed to the initial success of native plants in our treatment plots, as rhizophilic AMF are known to colonize plant roots and assist in establishment via protecting their plant host from damage caused by root pathogens [82].
Although plant invasions and land-use practices, such as tillage and herbicide applications, have consequences for soil microbial biomass and activity, some studies have found that soil microbial communities are largely resilient to ecological restoration techniques. In our study, bacterial abundance was significantly indistinguishable at the donor site versus treatment plots during the two study years. Control and dusting had significantly higher abundance in 2016 and in 2017. In 2016, dusting had significantly higher abundance than all other treatment plots, suggesting that the inoculation via dusting may have worked, and these treatment plots may contain microbes from the donor site, as well as the recipient site, contributing to this high abundance.
Although soil profiles are rarely homogenous, soil thickness or depth in the soil profile has been an important gradient structuring root distribution and bacterial community composition and functioning [85]. A recent study by Brewer et al. (2019) illustrated that bacterial attributes associated with microbial survival in oligotrophic conditions may promote their survival and capacity to thrive in subsurface soil horizons. Furthermore, microbial communities at depth were more dissimilar to surface microbial communities, along with greater distance from the surface within a soil profile [86]. Our study examined how microbial (AMF) communities and bacterial abundance shifts with topsoil thickness in recipient sites receiving topsoil from a reference coastal sage scrub ecosystem. Treatments with the smallest amount of inoculum (light dusting) yielded the greatest effects on increasing bacterial abundance after one year. Furthermore, treatments with the greatest topsoil depth (≤15 cm) corresponded to the lowest values of bacterial biomass. This suggests that inoculum size may play a role in affecting microbial activity, as well as potentially influencing fungal:bacterial abundance, along with associated ecosystem functioning. Previous studies examining bacteria and soil depth have shown that microbial functions, including microbial biomass normalized potential extracellular enzyme activity [87], and the proportion of Gram + bacteria [82] in the microbial community increased with depth. Although our study extended the soil profile to include heightened topsoil thickness, rather than extending deeper into soil horizons, our findings suggest that adding topsoil could have consequences for functions driven by microbial processes, with associated implications on both aboveground and belowground ecosystems.

Methodological Implications
Salvaged soil restoration resulted in successful ecological restoration, through a combination of three mechanisms: the burial of existing non-native seeds, thereby hindering their germination; the addition of a native seed bank as contained within the donor soil, resulting in higher native diversity; and the transfer of soil bacteria, as well as AMF diversity and rhizophilic AMF taxa richness.
Methodologically, there seems to be a critical threshold of 5 cm, above which the application of thicker layers of salvaged soil, whether 10 or 15 cm, provides minimal additional benefits. Both native density and native species richness were significantly different from dusted and control treatments at this critical 5 cm threshold of added salvaged soil; the same trend was observed in the ability of salvaged soil to reduce the cover of non-native grasses. Although a dusting of salvaged soil was sufficient for affecting bacterial abundance, AMF diversity was more variable and may require the addition of more topsoil. The most diverse AMF community was found within plots containing greater volumes of salvaged topsoil (15 cm treatments); however, these effects were only evident at one recipient site. Our study demonstrated variability in AMF communities and diversity that was likely context-dependent and driven by site-specific characteristics. However, across all sites, we detected shifts in AMF communities at ≥5 cm of added salvaged soil.
Given the increasing invisible costs associated with transporting and deploying greater masses of soil from donor sites, determining a threshold is critical for maximizing our restoration efforts for restoring the greatest area of degraded land with the least volume of outside inputs. Therefore, this study provides valuable information about how ≥5 cm thick layers of salvaged soil were generally sufficient for obtaining beneficial effects of this technique. Through the combined effect of the three contributing mechanisms, the successful application of this singular method allows for the collective and holistic restoration of diverse aboveground and belowground ecosystems.
Supplementary Materials: The following are available online at http://www.mdpi.com/1424-2818/12/4/150/s1. Table S1: Soil deposition volume. Figure S1: Experimental block layout at recipient sites. Figure S2: Representation of point-intercept data collection method. Table S2: ANOVA results for plant community data. Table S3: Tukey-Kramer post hoc results for plant community data. Table S4: Tukey post hoc results for bacterial abundance. Table S5 Funding: This research was not supported by funding agencies.