1. Introduction
During rapid economic and social development, large-scale natural forest harvesting and conversion of natural forests into plantations commonly occur, resulting in forest degradation, soil erosion, and other forest and environmental problems [
1]. To circumvent land degradation, the establishment of ecological welfare forestland was initiated in China, which is conducive to improving the ecological environment and the ecosystem services and functions of forests, which is of great significance to the restoration of fragile ecosystems and expansion of forest plantations [
2,
3]. Since China implemented classification-based forest management, a number of plantations with highly important ecological locations or extremely fragile ecological conditions has been established as ecological welfare forests, including water conservation forests, soil and water conservation forests, windbreak and sand dune fixation forests, and shore restoration forests [
4]. These ecological welfare forests play an important role in territorial ecological security, biodiversity conservation, and sustainable economic and social development [
5]. However, the quality of some ecological welfare forests in China is generally poor, and they have become low-yielding and inefficient with poor ecosystem stability and ecological functions, making it difficult to give full play to the ecological protection function of forests [
6].
The practice of artificial afforestation shows that continuous growing of monoculture forests on the same land results in decline of soil fertility, reduction in forest stand productivity and the reduction of forest ecological functions [
7]. Thus, mixed-species forest is opined to be a good model for afforestation and forest stand quality improvement in today's ecological welfare forest. Several studies have shown the advantages with mixed-species forests over monoculture forests. For instance, Maire et al. [
8] analyzed the light energy utilization efficiency of pure and mixed forests of
Eucalyptus grandis Hill ex Maiden and
Acacia mangium Willd, and found that the mixed forest canopy was more favorable to the growth of plants in the stands that use sunlight to provide nutrients for growth. Forrester et al. [
9] established 1:1 stand of
E. robusta and
Acacia mearnsii De Wilde and found that mixed stand improved stand water use efficiency, canopy photosynthetic capacity, and aboveground and underground carbon distribution. Liao et al. [
10] found that a mixed forest of
Eucalyptus grandis x
urophylla and
A. mangium was conducive to forest growth in a relatively barren red-soil mountainous area, and its ecological benefits, wind resistance, and disease resistance, were better than those of a pure forest was. Chen et al. [
11] found that the mixture of
E. urophylla and
A. mangium can effectively increase the soil nutrient content and that the nutrient content of litter in mixed forests exceed that of pure forests. The study by Zang et al. [
12] showed that the number and diversity of micro-organisms in the mixed forest of
Acacia crassicarpa Benth and
E. robusta were greater than those in the pure forest, and afforestation of mixed
A. crassicarpa and
E. robusta was beneficial to improve soil fertility and vegetation diversity.
Acacia cincinnata F. Muell is an evergreen tree of the legume family, native to the northeastern coastal region of Australia, introduced to China in the 1970s, and widely planted in Guangdong, Fujian, Hainan Provinces [
13].
A. cincinnata has the advantages of fast growth, strong nitrogen fixing ability, resistance to drought and barrenness, wide timber use, and short rotation of 6 - 8 years. Therefore, it has been widely used for montane afforestation, soil and water conservation afforestation, and coastal conservation afforestation in southern Fujian, China [
14]. Owing to its production of a large leaf litter,
A. cincinnata has a good effect on improving soil fertility. It can also be used in crop rotation planting and suitable for establishment of mixed forest with
E. robusta in southern Fujian to effectively improve the consumption of soil nutrients, fertilizer, and water on forestland due to the rapid growth of eucalyptus trees [
15]. At the same time,
A. cincinnata is also an excellent choice for establishment of mixed -forests with acacia species such as
A. mangium, to effectively improve the stability and diversity of the forest stand structure [
16]. Due to the long-term focus on the expansion of afforestation areas of
A. cincinnata and insufficient attention to the improvement of stand management practices, especially neglecting stand-tending operation, many
A. cincinnata forests in southern Fujian have become low-yielding, with poor ecosystem stability and ecological functions [
17]. It is, therefore, difficult to truly achieve the multi-functional benefits of ecological welfare forests in maintaining water conservation, soil fertility, biodiversity conservation and recreational landscape effects [
18].
Therefore, the establishment of A. cincinnata plantations mixed with E. robusta and A. mangium may be conducive to improving the soil fertility and water conservation capacity of barren land, which is of great significance for improving the ecological and economic benefits of forestland in southern China. However, the effects of these mixtures on growth of tree species, undergrowth vegetation diversity, soil physico-chemical and biological properties of A. cincinnata plantation are largely unexplored. Thus, the aims of this study were to examine whether mixture of A. cincinnata with A. mangium or E. robusta will have an impact on growth of tree species, if so which species, A. mangium or E. robusta has the largest impact; and to evaluate the effects of mixed species planting on understory vegetation, soil and bacterial community structure and diversity. The soil bacterial community was explored by sequencing of the 16S rRNA gene in soils from different forest stand types. Our research study endeavored to answer four main questions: (i) Does mixed-species planting improve growth of tree species? (ii) Does mixed-species planting increased the diversity and biomass of the understory vegetation as well as C, N and P concentrations in the soil, (iii) Is the bacterial community of pure and mixed-species stands different due to their higher nutrient inputs (litter fall and biomass) in the mixed stands? (iv) Is there any link between bacterial community and soil attributes (C, N, P, C:N, C:P and N:P ratios? The study will provide an insight into the scientific management of A. cincinnata mixed forest to meet a range of economic and ecological benefits in southern China.
2. Materials and Methods
2.1. Study site
The study was conducted in Daqiutian Wutai Mountain Forest Farm, Nan'an county, Quanzhou, Fujian, China (118°24′20″N, 25°16′45″E), which is located at an altitude of 80~1000 m. The area has a subtropical monsoon climate with an annual average temperature of 19.5 to 21.0 °C with a frost-free period of 350 d. The soil is red soil developed from sandstone shale with low soil fertility. The plantations were established in 2014, and the stand density was 1110 trees hm-2 of E. robusta (40%; 444 trees hm-2) and A. cincinnata (60%; 666 trees hm-2) forest, 1850 trees hm-2 of A. mangium (25%; 463 trees hm-2) and A. cincinnata (75%; 1387 trees hm-2) forest and 1230 trees hm-2 of A. cincinnata forest (100%), with a total area of about 6.13 hm2, 7.6 hm2, and 7.06 hm2, respectively.
2.2. Measurement of growth traits
In January 2021 (7 years after planting), three 20 × 20 m standard plots were set- up in each forest stand type. Thereafter, the diameter at breast height (DBH) and tree height of each tree were measured in the standard plot and the average volume and stand volume of each plot were calculated. The DBH of trees was divided into DBH class of 2 cm interval [
19,
20]. The number of tree in each diameter class was counted, and the average DBH of trees in each diameter class was calculated separately for each forest stand as follows:
Here, N is the total number of trees in the i-th diameter class, and di is the DBH of the i-th diameter class.
To improve the accuracy of the tree height measurements, the measured mean DBH was substituted into the growth model and then we fitted the tree height as follows:
where y and x is the fitted height and DBH, respectively
The volume per tree was calculated as follows:
For Acacia, V = 0.00005276×D1.882161×H1.009317
For E. robusta, V = 0.000109154×D(1.87892-0.00569186×(D+H))×H(0.652598+0.00784754×(D+H))
Where D is the average DBH, H is the fitted tree height.
To calculate the volume of each diameter class, the average volume per tree of each diameter class was multiplied by the number of trees in the diameter class. To calculate the stand volume, the volume of each diameter class was added up in the sample plot.
2.3. Biomass and nutrient concentrations of understory vegetation
The species richness of the understory vegetation was determined by conducting an inventory in three 1 m × 1 m plots in each forest stand type. All species were identified in situ using plant identification manual and its cover was determined. The aboveground biomass of understory vegetation (herbaceous plants and shrubs) was determined by harvesting all individuals in three 1m×1 m plots that were set-up in each forest stand. The underground (root) biomass was determined by digging each plant roots. The fresh biomass samples were weighed in the field and then the samples brought back to the laboratory where they were oven-dried at 105 °C to a constant mass for the determination of dry mass. After crushing the samples and passing through 0.15 mm sieve to remove impurities, the carbon (C) and nitrogen (N) concentrations were determined by Vario Max Elemental Carbon and Nitrogen Analyzer (Elementar, Germany), and the Molybdenum-antimony colorimetry method was used to determine the total phosphorus concentration in the biomass of understory vegetation.
2.4. Soil analysis
Soil samples were collected from tree plots in each stand at different soil depth: 0-10 cm, 10-20 cm and 20-40 cm. The soil physical properties, including soils water content, bulk density, water holding capacity and non-capillary and capillary porosity, were determined using the ring knife method, keeping the soil structure intact. The chemical properties were measured using air-dried soil (sieved at 0.15 mm). Carbon and nitrogen contents were measured with an elemental analyzer (Vario MAX CNS, Elementar, Hanau, Germany) while P was determined after digestion by HF, HClO4, and HNO3 using the molybdenum-blue method.
2.5. Extraction and Amplification of Soil Bacterial DNA for Sequencing
The DNA in soil samples was extracted using the CTAB method, and the quality of the extracted DNA was checked using electrophoresis in agarose gels (1% w/v in TAE buffer). After electrophoresed for 30 min, the samples were thawed on ice, centrifuged and thoroughly mixed. The quality of the samples was detected by Nanodrop, and 30 ng sample were taken for PCR amplification. PCR products were purified using the Agencourt AMPure XP Nucleic Acid Purification Kit. The V3–V4 hypervariable regions of the 16S rDNA of the bacteria were amplified. The primer sequences were ACTCCTACGGGAGGCAGCAG and GGACTACHVGGGTWTCTAAT. The amplification was conducted under the following reaction conditions: 5 min of initial denaturation at 94 °C, 30 s of denaturation at 94 °C, 30 s of annealing at 55 °C, and 28 cycles of 1 min of elongation at 72 °C. The amplified PCR products were sequenced on an Illumina MiSeq (PE300) sequencing platform. The original sequencing data were deposited in the NCBI SRA database under the accession number PRJNA954280. The high-quality Clean Tags sequences were obtained after quality control filtering, and the Clean Tags sequences were clustered to produce operational taxonomic units (OTUs) based on 97% sequence similarity. The UNITE taxonomic database and the RDP Classifier database were used for species annotation of OUT representative sequences and analysis of bacterial community composition at each taxonomic level to obtain information on bacterial species in soils of different forest stand types.
2.6. Statistical analysis
The obtained OUT clustering results and species information were analyzed using Mothur software, and the abundance index of Chao 1 and Observed species, phylogenetic diversity (PD_whole_tree) and Shannon’s diversity index were calculated. One-way analysis of variance (ANOVA) was conducted to examine significant differences in growth traits of A. cincinnata, biomass of understory vegetation, and bacterial diversity among the different forest stand types while two-way ANOVA was conducted to determine significant differences in soil physico-chemical properties among different soil depths and stand types. Means that exhibited significant differences were compared using Tukey’s honestly significant test using SPSS software (version 22, IBM, New York, USA). The relationship between soil C, N and P content and the relative abundance of the dominant bacterial phyla was analyzed by redundancy analysis using Canoco software (version 5). All data were expressed as the mean ± standard error (SE).