ORIGINAL RESEARCH article

Front. Microbiol., 15 December 2022

Sec. Terrestrial Microbiology

Volume 13 - 2022 | https://doi.org/10.3389/fmicb.2022.1065899

Photovoltaic panels have altered grassland plant biodiversity and soil microbial diversity

  • 1. College of Animal Science and Technology, Northeast Agricultural University, Harbin, China

  • 2. College of Animal Science and Veterinary Medicine, Heilongjiang Bayi Agricultural University, Daqing, China

Abstract

Introduction:

Human concerns about fossil fuel depletion, energy security and environmental degradation have driven the rapid development of solar photovoltaic (PV) power generation. Most of the photovoltaic power generation plants are concentrated in desert, grassland and arable land, which means the change of land use type. However, there is still a gap in the research of the PV panel layout on grassland plant species diversity and ecological function.

Methods:

In this study, Illumina high-throughput sequencing technology was used to investigate the effects of PV panel arrangement on grassland plant species diversity and soil microbial diversity. In view of the differences in the microclimate at different sites of the PV panels, quadrates were arranged in front edge (FE), beneath the center of each panel (BP), back edge (BE), the uncovered interspace adjacent to each panel (IS) and the undisturbed grassland around the PV panels (Control), respectively.

Results:

PV panels (especially FE) significantly increased the total aboveground productivity (total AGB) and plant species diversity in grasslands. FE increased precipitation accumulation and plant species diversity directly and indirectly changed the diversity of soil bacterial and fungal communities. PV panels decreased the relative abundance of Actinobacteriota, while increased the relative abundance of Proteobacteria, Acidobacteriota, and Methylomirabilota. EC, Margalef’ s richness and total AGB were the main factors affecting the composition of bacterial communities, while alkaline hydrolysis nitrogen (AN) and available phosphorus (AP) were the main factors affecting the composition of fungal communities.

Discussion:

In conclusion, the arrangement of PV panels increased the plant species diversity and soil microorganisms in grassland. This study provides important information for further understanding the impact of PV panels on grassland ecosystem function and is of great significance for maintaining grassland ecosystem function.

Introduction

Human concerns over fossil fuel depletion, energy security and environmental degradation have led to an increasing demand for clean renewable energy (Ding et al., 2016). The two outstanding characteristics of zero pollution and zero emissions make solar photovoltaic power (PV) a better energy source and an ideal alternative to traditional fossil fuels (Joshi and Arora, 2017). By the end of 2019, the global total installed capacity of PV had reached 505 GW (REN21, 2019), and it is still growing rapidly. PV generation has become an important pillar of national energy development strategies (Wu et al., 2021). However, PV is spatially intensive, large-scale and non-integrated, and PV deployment is estimated to require more land area than coal nuclear or natural gas technologies (Fthenakis and Kim, 2009). At present, the construction of PV plants is mostly concentrated in desert, grassland and arable land, which means the change of land use types (Mac Kay, 2013). Grassland is the world’s largest terrestrial ecosystem and an important part of the global natural ecosystem. They play an important role in the development of animal husbandry, windbreak and sand control, soil and water conservation, biodiversity protection and ecological balance (Li et al., 2006; Wen et al., 2013). Although PV generation has many social functions, such as helping to achieve low emission rate, alleviating energy crisis, improving the efficiency of resource utilization, there is still a gap in the research of PV panel layout on grassland diversity and ecological function.

The physical presence of PV panels will affect solar radiation flux (temperature), wind speed and turbulence (potential evaporation) and precipitation distribution under PV panels (Armstrong et al., 2014). A large number of studies have shown that PV panels reduce the amount of solar radiation received by the local surface and the atmospheric temperature in the growing season (from April to August) by converting part of the solar energy into electric energy (Armstrong et al., 2014; Liu et al., 2019; Yue et al., 2021). The shielding effect of PV panels leads to uneven precipitation distribution (Elamri et al., 2018; Li Y. et al., 2018), the presence of PV panels can concentrate water at its lower edge, which increases the local heterogeneity of soil water distribution and creates more permanent water storage under PV panels (Adeh et al., 2018; Yue et al., 2021). In addition, the lower potential evaporation caused by the reduction of solar radiation directly below the PV panels will lead to a reduction in drought, an increase in soil water availability, and help reduce water loss in arid climates (Marrou et al., 2013; Adeh et al., 2018; Choi et al., 2020).

The linkage changes of species interactions under climate change affect the biodiversity and function of terrestrial ecosystems (Walther et al., 2002; Gottfried et al., 2012; Langley and Hungate, 2014; Rudgers et al., 2014; Crowther et al., 2016). Given that light, temperature and water are the three key factors for plant growth and ecosystem function (Martens et al., 2000; Dupraz et al., 2011; Li W. H. et al., 2018), individual plant species differ in their ability to feedback climate change or survive under changing conditions through adaptation or plasticity (Fernandez-Going et al., 2013; Shi et al., 2015). The interaction between local microclimate changes caused by PV panels may affect plant community structure (Cleland et al., 2004; Adler et al., 2006; Yang et al., 2011), and directly/indirectly affect microbial community structure through changes in the plant community composition (Ren et al., 2018; Toju et al., 2018; Shi et al., 2020). Therefore, it is necessary to better understand the human interference (i.e., the laid of pv panels), climate change and bio-abiotic interactions affect how grassland ecosystem structure and function. This will provide a theoretical reference for the selection and construction of scientific, reasonable and efficient PV power plants, and provide scientific basis for vegetation and soil restoration after PV power plants construction.

Grassland ecosystems are considered to be one of the most sensitive ecosystems to climate change (Piao et al., 2012; Chen et al., 2013). Microclimate change caused by human disturbance will have a profound impact on grassland ecosystem function. Therefore, understanding the impact of PV panels on grassland ecosystem is of great significance for maintaining grassland ecosystem function. In this study, the PV power plant is located in Datong District, Daqing City. In the past, many large areas of grassland were constructed as PV power stations. The objectives of this study were: (1) to investigate the effects of PV panels laying on grassland productivity, plant community composition and species diversity, soil properties, microbial community composition and diversity; (2) to clarify the intrinsic relationship between plant species diversity, soil properties and soil microbial changes.

Materials and methods

Study site

The study area is located in the western Songnen Plain, which belongs to Datong District of Daqing City (46°10 ‘11’ N, 124°53 ‘56’ E; 135 m asl), which has a mid-temperate continental monsoon climate. Light is sufficient, precipitation is less, wind and drought in the same period, rain and heat in the same season. The mean annual temperature is 4.2°C, the mean temperature of the coldest month is-18.5°C, and the extreme minimum temperature is-39.2°C. The average temperature of the hottest month is 23.3°C, and the extreme maximum temperature is 39.8°C. The average annual precipitation is 400-500 mm, and the precipitation period is mostly concentrated in June to September. The average annual frost-free period is 143 days. The average annual wind speed is 3.8 m/s. The annual evaporation is 1,635 mm, the annual dryness is 1.2, the annual sunshine duration is 2,726 h, and the annual total solar radiation is 491.4 kJ/cm2. The main soil types were meadow and marsh soil and the grassland type in the study area is warm meadow grassland, and the main plant types are Puccinellia tenuiflora (Griseb.) Scribn. et Merr., Phragmites australis (Cav.) Trin. ex Steud. and Suaeda glauca (Bunge) Bunge.

Experiment design

The PV power station was randomly divided into three areas, and three PV panels were randomly selected in each area, making a total of nine PV panels as the research object. In addition, an undisturbed grassland around each area was selected as control, and a total of three control grasslands were selected, with a spacing of 20 m between the control grassland and the study areas. In view of the differences in the microclimate at different sites of the PV panel, quadrates were arranged in front edge (FE), beneath the center of each panel (BP), back edge (BE), the uncovered interspace adjacent to each panel (IS) and the undisturbed grassland around the PV panels (Control), respectively (Figure 1). The size of each quadrat is 1 × 1 m, four quadrats are arranged under each PV panel, and one quadrat is arranged in each control area, a total of 39 quadrats.

Figure 1

Plant sampling and analysis

During the study, all plant data were measured at the end of July, when the plants were in full bloom. The plant community was divided into four functional groups:PG (perennial grasses); PS (perennial sedges); PF (perennial forbs); ABH (annual and biennial herbs). After litter removal, we measured total plant community coverage, height, coverage and density of each species, and density of each plant functional group within each quadrat. The height of each plant functional group was obtained by calculating the mean of plants within the group. Each plant was then trimmed to floor level and placed in marked paper bags by species. Then it was dried at 65°C for 48 h and weighed to obtain the aboveground biomass (AGB) of each species, and AGB of each plant functional group and total AGB were calculated. Plant Margalef’s index, Shannon-Wiener index, Simpson index were used to illustrate biodiversity in this study. Margalef’s richness index was calculated as follows:

where N is the total number of individuals of all species, and S is the number of species. Shnnon-Wiener index was calculating as follows:

where Pi is the proportion of individual species i representing the relative density of plant species (species density/total density for all species × 100). Simpson diversity index was calculated as follows:

Soil Sampling and Analysis

After measuring the AGB parameters, soil was sampled at five points with a soil coring machine. Soil samples were collected at a depth of 0–20 cm. Five samples from the same quadrat at the same depth were thoroughly mixed to remove root and gravels, quickly packed into sterile sealable bags, placed in ice boxes for refrigeration, and promptly shipped to the laboratory. The sample was divided into two parts. Part of the soil was air-dried, ground and sieved to determine the physical and chemical properties of the soil. The remaining parts were stored at-80°C for soil microbial analysis.

Soil nutrient determination

Soil pH and electrical conductivity (EC; an important indicator of salinity) were determined with soil (air-dried)/water (1:5, w/v) suspensions (Bao, 2008), and soil organic matter (SOM) was assessed by K2CrO7 volumetric method (External heating method) (Sparks et al., 1996). Total nitrogen (TN) and alkali-hydrolyzed nitrogen (AN) were measured by Kjeldahl method and alkali-hydrolyzed diffusion method (Lu, 1999). Phosphorus (P) was extracted by HClO4-H2SO4 method, then added the molybdenum-antimony anti-reagent, colorimetric was performed at 880 nm or 700 nm wavelength, and calculated the total phosphorus (TP) content (Bao, 2008). Available phosphorus (AP) was extracted with 0.5 M NaHCO3 (pH 8.5), and was determined using molybdenum blue method (UV-752 Shanghai, China; Bao, 2008). Determination of total potassium and available potassium by flame photometry (Bao, 2008).

Soil DNA extraction, illumina sequencing, and data analysis

Total DNA was extracted from 0.25 g of Soil using the FastDNA® Spin Kit for Soil (MP Biomedicals, USA) isolation kit according to the manufacturer’s requirements. After DNA extraction, the extracted genomic DNA was detected by 1% agarose gel electrophoresis. For MiSeq library construction, bacterial V3 and V4 hypervariable regions were amplified using primers 338F (5′ -ACTCCTACGGGAGGCAGCA-3′) and 806R (5′ -GGACTACHVGGGTWTCTAAT-3′). Primer ITS1F was used (5’-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5’-GCTGCGTTCTTCATCGATGC-3′) were used to amplify the hypervariable region of fungal ITS1. Specific primers with barcode were synthesized according to the specified sequencing region. The amplification process was pre-denaturation at 95°C for 5 min, 25 cycles of 95°C for 30 s, 50°C annealing for 30 s, 72°C for 40 s, and 72°C for 7 min of extension. PCR products of the same sample were mixed and detected by 2% agarose gel electrophoresis. PCR products were recovered by gel cutting using AxyPrepDNA gel recovery kit (AXYGEN Company) and eluted with Tris_HCl. Electrophoresis was performed on 2% agarose. PCR products were quantified with the QuantiFluor™ -ST Blue fluorescence Quantification System (Promega). The PCR products were purified with TruSeqTM DNA Sample Prep Kit to construct the Miseq library. Illumina Miseq platform was used for on-machine sequencing of the constructed library.

The reads of each sample were spliced using FLASH (version 1.2.11) software (Magoc and Salzberg, 2011), primers were removed, and sequences were trimmed to remove low-quality sequences. Sequences were clustered at 97% similarity level (OTU) and screened using USEARCH (version 10.0) with a threshold of 0.005% for all sequences (Edgar, 2013). The RDP software (Version 2.2) was then used for comparison with the Silva and UNITE databases, with confidence intervals of 80% for bacterial and fungal classification (Wang et al., 2007; Quast et al., 2013). Mothur software was used to calculate Chao1, Simpson diversity, Shannon’s evenness, and Good’s coverage to estimate bacterial and fungal diversity and richness.

Statistical analysis

All statistical analyses were performed using R (4.1.2) software. One-way analysis of variance (ANOVA) and TukeyHSD test in “Agricolae” package were used to analyze the significant differences between groups in plant diversity, soil properties, soil bacterial α diversity and relative abundance of dominant phyla, and ANOVA was used to test the effects of PV panels on plants, soil, bacteria and fungi. Spearman rank correlation coefficient was calculated to analyze the correlation between vegetation, soil properties and soil microorganisms. Based on Bray-Curtis distance matrix, the similarity analysis (ANOSIM) of 999 permutations was performed, and the PCoA visualization of plant community composition at different sites was performed. Based on Bray-Curtis distance matrix, NMDS was used to visualize the composition of bacterial and fungal communities at different sites. Linear discriminant analysis (LDA) effect size (LEfSe) method was used to identify significantly different characteristics of bacteria and fungi at different sites by biomarkers, and the impact size of each feature was assessed with a threshold of 3 and a significant α of 0.5. The “vegdist” package was used to conduct RDA analysis of plant diversity, soil properties and bacterial and fungal communities, and 999 permutation similarity analysis (ANOSIM) was performed to infer potential relationships between microbial community composition and measured plant diversity and soil properties. SEM models were constructed using Amos 26 (AMOS IBM, USA) to assess the direct and indirect effects of vegetation and soil properties on bacterial and fungal community diversity.

Results

Effects of PV panels on plant community and soil properties

PV panels had significant effects including the Margalef’s richness index, Shnnon-Wiener index and Simpson diversity index of the plants. From FE to IS, BP, BE and Control, most diversity indices decreased (p < 0.05). For different sites under the PV panel, the diversity of FE plants increased the most significantly, while the diversity of BE plants increased the least. From the perspective of homogeneity, compared with IS, most diversity indices increased in FE and decreased in BP and BP (Table 1; Supplementary Figure S1). PV panels had significant effects on total AGB and community total coverage (p < 0.05) (Table 2), and total AGB decreased from FE to BE, IS, BP, and Control. Compared with Control, total AGB under FE, BE, IS, and BP increased by 80.02, 52.01, 43.04, and 17.10%, respectively (Table 2; Supplementary Figure S2). PV panels significantly increased the total coverage of plant communities, which decreased from FE to BP, IS, BE and Control. Compared with Control, the total plant community coverage under FE increased by 80.02% (Table 2). Compared with IS, total AGB and total coverage increased significantly under FE and decreased significantly under BP, but there was no significant difference under BE. Compared with Control, the total coverage of plant communities increased most significantly under FE, increasing by 80.02%. However, the total coverage of plant communities increased the least under BP, with an increase of 17.10% (Table 2).

Table 1

SiteFP
ControlISFEBPBE
Plant α diversityMargalef’s richness0.43 ± 0.08 b0.95 ± 0.21 ab1.04 ± 0.35 a0.72 ± 0.33 ab0.71 ± 0.34 ab3.3430.021
Shnnon-Wiener0.56 ± 0.34 c1.32 ± 0.27 ab1.46 ± 0.53 a0.96 ± 0.35 abc0.87 ± 0.38 bc5.2070.002
Simpson diversity0.31 ± 0.23 c0.66 ± 0.07 ab0.70 ± 0.18 a0.49 ± 0.17 abc0.43 ± 0.24 bc4.9250.003
Bacteria α diversityChao1 richness2866.99 ± 150.05 c4023.72 ± 280.48 ab4680.33 ± 369.26 a3595.27 ± 714.17 bc3776.64 ± 670.66 bc8.781<0.001
Simpson diversity0.03 ± 0.00 a0.01 ± 0.01 b0.01 ± 0.00 b0.01 ± 0.01 b0.01 ± 0.01 b6.368<0.001
Shannon’s evenness4.87 ± 0.08 c6.05 ± 0.43 ab6.43 ± 0.43 a5.75 ± 0.51 b5.98 ± 0.47 ab7.474<0.001
Good’s coverage0.99 ± 0.00 a0.99 ± 0.00 ab0.98 ± 0.00 b0.99 ± 0.01 ab0.98 ± 0.00 ab2.7300.045
Fungi α diversityChao1 richness170.43 ± 36.09 b290.54 ± 87.39 ab447.78 ± 176.34 a292.31 ± 163.00 ab318.92 ± 115.74 ab3.0650.029
Simpson diversity3.09 ± 0.52 a3.15 ± 0.67 a3.35 ± 0.62 a3.22 ± 0.61 a2.93 ± 0.41 a0.4250.789
Shannon’s evenness0.14 ± 0.12 a0.13 ± 0.09 a0.13 ± 0.08 a0.11 ± 0.06 a0.17 ± 0.11 a0.5910.672
Good’s coverage1.00 ± 0.00 a1.00 ± 0.00 a1.00 ± 0.00 a1.00 ± 0.00 a1.00 ± 0.00 a2.2280.087

One-way ANOVA analysis of the effects of PV panels on plant species, soil bacterial, and fungal diversity (FE: front edge of each panel, BP: beneath the center of each panel; BE: back edge of each panel; IS: the uncovered interspace adjacent to each panel; Control: the undisturbed grassland around the PV panels).

Table 2

ControlISFEBPBEFP
Total AGB312.32 ± 53.99b446.74 ± 168.52ab565.24 ± 102.21a365.73 ± 130.26b474.77 ± 171.52ab3.0390.030
Total coverage50.00 ± 5.00b66.11 ± 9.28ab78.11 ± 7.04a68.33 ± 11.73ab65.00 ± 10.61b5.437<0.001
PG coverage46.33 ± 5.51a42.61 ± 16.96a47.62 ± 25.17a41.00 ± 22.06a38.11 ± 20.51a0.2740.893
PF coverage0.00 ± 0.00a3.00 ± 6.78a11.62 ± 13.13a5.89 ± 14.30a0.06 ± 0.17a1.8340.145
PS coverage0.00 ± 0.00a13.39 ± 13.93a12.01 ± 14.56a1.00 ± 1.80a7.67 ± 11.30a2.0570.108
ABH coverage5.33 ± 3.06a19.39 ± 17.49a34.14 ± 16.01a24.89 ± 24.84a26.96 ± 23.66a1.3600.268
PG height219.92 ± 194.42a70.50 ± 15.37b65.91 ± 17.74b50.31 ± 18.67b65.39 ± 11.92b7.114<0.001
PF height0.00 ± 0.00b6.67 ± 15.21ab21.00 ± 20.83a3.07 ± 6.13ab2.22 ± 6.67b3.1870.025
PS height0.00 ± 0.00b44.00 ± 23.69a29.39 ± 24.80ab9.44 ± 17.40b39.78 ± 23.77a4.6870.004
ABH height37.00 ± 24.64a28.88 ± 16.96a35.26 ± 8.74a25.51 ± 12.44a43.47 ± 17.84a1.8550.141
PG abundance275.33 ± 200.14a178.89 ± 198.77a180.89 ± 119.44a149.00 ± 149.12a118.22 ± 132.35a0.6430.636
PF abundance0.00 ± 0.00a10.11 ± 29.96a30.33 ± 38.80a13.00 ± 26.23a0.33 ± 1.00a1.6210.192
PSabundance0.00 ± 0.00a30.11 ± 29.54a28.11 ± 28.79a4.89 ± 10.34a28.22 ± 29.45a2.1170.100
ABH abundance24.00 ± 23.30a124.01 ± 147.08a239.11 ± 210.53a249.22 ± 261.02a195.89 ± 212.45a1.0390.402
PG frequency1.31 ± 0.00a1.60 ± 0.24a1.54 ± 0.35a1.27 ± 0.38a1.28 ± 0.46a1.6350.188
PF frequency0.00 ± 0.00a0.03 ± 0.06a0.20 ± 0.20a0.09 ± 0.18a0.02 ± 0.06a2.8790.037
PS frequency0.00 ± 0.00a0.55 ± 0.22a0.43 ± 0.33ab0.20 ± 0.30ab0.41 ± 0.28ab3.4440.018
ABH frequency0.41 ± 0.27a1.18 ± 0.57a1.06 ± 0.61a0.94 ± 0.61a1.01 ± 0.37a1.2110.324
PG AGB288.50 ± 54.22a207.13 ± 114.39a241.59 ± 100.44a165.19 ± 81.09a183.98 ± 107.06a1.2960.291
PF AGB0.00 ± 0.00a14.06 ± 41.69a24.45 ± 29.10a8.85 ± 20.67a0.32 ± 0.95a1.1090.368
PS AGB0.00 ± 0.00b77.87 ± 81.59a39.75 ± 56.85ab1.55 ± 2.81b31.99 ± 56.86ab2.4740.062
ABH AGB23.82 ± 19.87a137.42 ± 153.47a222.04 ± 180.98a190.14 ± 201.34a258.48 ± 276.91a0.9730.435

One-way ANOVA analysis of the effects of PV panels on plant communities and functional groups (PG: perennial grasses; PS: perennial sedges; PF: perennial forbs; ABH: annual and biennial herbs).

The results of PCoA based on Bray-Crutis showed that PCo1 and PCo2 explained 68.89 and 20.85% of the variation in plant community composition, respectively. The plant community composition was significantly separated between Control and PV panels, indicating that PV panels changed the plant community composition, and the plant composition at different sites under PV panels was significantly different (F = 26.235; p < 0.001). PV panels promoted the growth of PF, PS and ABH, while inhibited the growth of PG (R2 = 0.755, p = 0.001) (Figure 2; Table 2). PV panels had significant effects on the height and frequency of plant functional groups (p < 0.05). However, there were significant differences in different sites under the PV panels. The PF height under FE increased the most significantly, but there was no significant difference between BE and Control. The PS height increased most significantly under IS and BE, but there was no significant difference between Control under BP. There was no significant difference in PG height between different sites of PV panels. There was no significant difference in the coverage and richness of each functional group at different sites under the PV panels compared with the Control. FE, BP and BE decreased the frequency of PS, but there was no significant difference between IS and Control. The PV panel significantly increased PS AGB, and IS was the most significant, while BP was the least (Table 2).

Figure 2

Except for TP, PV panel had no significant effect on soil properties (p > 0.05). Although there was no statistical difference in soil properties between different sites under PV panels and Control. However, compared with the Control, the pH, EC, TK and AK contents at different sites under the PV panel were slightly decreased, while AN, TP and AP were slightly increased. For different sites under the PV panel, most soil properties were highest under BP (Table 3).

Table 3

SitepHEC (mS/cm)SOC (g/kg)TN (g/kg)AN (mg/kg)TP (g/kg)AP (mg/kg)TK (g/kg)AK (mg/kg)
Control8.90 ± 0.04a3.2 ± 1.24a10.06 ± 0.46a1.23 ± 0.47a150.07 ± 28.12a0.95 ± 0.38a175.52 ± 29.66a36.75 ± 1.24a144.26 ± 1.44a
IS8.59 ± 0.41a1.84 ± 0.65a9.98 ± 2.12a1.09 ± 0.36a159.72 ± 80.16a1.00 ± 0.23a197.89 ± 54.65a36.66 ± 1.49a142.97 ± 3.41a
FE8.54 ± 0.41a1.94 ± 0.95a11.10 ± 1.51a1.56 ± 0.56a156.29 ± 49.67a1.29 ± 0.46a189.22 ± 42.50a36.10 ± 1.86a143.21 ± 2.70a
BP8.59 ± 0.25a1.95 ± 1.13a11.15 ± 1.74a1.22 ± 0.43a193.34 ± 94.41a1.48 ± 0.41a185.51 ± 68.83a35.41 ± 2.194a141.62 ± 3.33a
BE8.62 ± 0.49a1.69 ± 0.65a10.89 ± 1.36a1.52 ± 0.72a183.17 ± 69.24a1.39 ± 0.23a190.99 ± 72.62a36.60 ± 1.15a140.31 ± 2.64a
F0.5191.6960.8591.2790.4772.9730.0970.8581.769
P0.7220.1740.4980.2980.7520.033*0.9830.4990.158

One-way ANOVA analysis of the effects of PV panels on soil properties.

Effects of PV panels on soil microbial community diversity and composition

A total of 2,751,094 and 4,879,890 effective bacterial and fungal sequences were obtained from 39 samples after Illumina high-throughput sequencing and filtering to remove low-quality sequences. The effective sequence length of bacteria was mainly distributed in the range of 400-440 bp, and that of fungi was mainly distributed in the range of 220-240 bp. There were 23,143 bacterial OTUs and 3,064 fungal OTUs in all samples. The bacterial OTUs of each sample ranged from 43,933 to 119,026, with an average of 70,541, and the fungal OTUs ranged from 104,110 to 148,006, with an average of 125,125, and the homology threshold was 97%. The Good’s coverage of bacteria and fungi of all samples were higher than 98.61 ± 0.48% and 99.97 ± 0.01%, respectively, indicating that the sequencing quality was good.

PV panels had significant effects including bacterial Chao1 richness, Simpson diversity, Shannon’s evenness and Good’s coverage, and fungi Chao1 richness (p < 0.05) (Table 1). As for bacteria, the Chao1 richness and Shannon’s evenness of bacteria increased by PV panels, including the most significant increase in FE, but the least increase in BP, and no significant difference between IS and BE. The Simpson diversity of bacteria was decreased by PV panels, and there was no significant difference among different sites of PV panels (p > 0.05). For fungi, Chao1 richness was significantly improved by PV panels, including FE which was the most significant, but there was no significant difference among IS, BP, and BE. The fungi Simpson diversity and Shannon’s evenness were not significantly affected by PV panels (Table 1; Supplementary Figure S3). Spearman rank correlation analysis showed that bacteria Chao1 richness was positively correlated with Margalef’s richness, Shnnon-Wiener and Simpson diversity. Shannon’s evenness was positively correlated with Margalef’s richness, Shnnon-Wiener and total AGB, and negatively correlated with EC. Simpson diversity was negatively correlated with Margalef’s richness, Shnnon-Wiener and total AGB, and positively correlated with EC. Fungi Chao1 richness was positively correlated with AP, Shnnon-Wiener and Simpson diversity, but negatively correlated with EC (Table 4).

Table 4

TNTPANAPPHECSOCTKAKMargalefShnnonSimpsonTotal AGB
BacteriaChao10.2720.037−0.0730.124−0.266−0.2350.053−0.061−0.0460.5450.6140.6460.183
Shannon0.0760.0560.0830.105−0.189−0.4160.102−0.026−0.0600.3670.4810.5460.346
Simpson−0.0690.016−0.029−0.0120.1860.390−0.0450.0320.024−0.379−0.525−0.579−0.356
Good’s Coverage−0.126−0.118−0.124−0.1340.2210.291−0.1510.0180.148−0.241−0.384−0.384−0.251
FungiChao1−0.0450.2490.1780.372−0.154−0.3380.2260.031−0.0990.2430.3200.3170.120
Shannon−0.119−0.1400.0760.195−0.0610.1130.0620.043−0.1550.2120.1600.036−0.192
Simpson0.1510.2670.003−0.0740.013−0.1090.068−0.1240.142−0.149−0.165−0.0760.193
Good’s Coverage0.076−0.103−0.145−0.423−0.0180.124−0.046−0.308−0.043−0.063−0.187−0.196−0.021

Spearman rank correlation analysis of soil bacteria and fungi with plant species diversity and soil properties.

Based on Bray-Curtis distance matrix, we performed NMDS visualization for bacteria and fungi in different sites. Bacterial and fungal Stress were 0.0906 and 0.1673, respectively. Permanova analysis results showed that the bacterial community composition was significantly different among different sites (R2 = 1.564, p = 0.030), while the fungal community composition was similar among different sites (R2 = 1.020, p = 0.409). However, PV panels did not significantly affect the bacterial and fungal community composition, and bacterial and fungal community composition was similar between Control and different PV panel sites (Supplementary Figure S4).

At the phylum level, the relative abundance of the top 10 dominant bacteria in all sites was higher than 91% (29.77% Actinobacteriota, 18.98% Proteobacteria, 14.88% Chloroflexi and 14.15% Acidobacteriota, 5.01% Firmicutes, 3.27% Gemmatimonadota, 2.86% Myxococcota, 2.15% Methylomirabilota, 1.56% Bacteroidota and 1.23% Desulfobacterota). PV panels had significant effects on the relative abundance of Actinobacteriota, Proteobacteria, Acidobacteriota and Methylomirabilota (p < 0.05), decreased the relative abundance of Actinobacteriota, while increased the relative abundance of Proteobacteria, Acidobacteriota, and Methylomirabilota. There was no significant difference in the relative abundance of Actinobacteriota among different sites of PV panels. The relative abundance of Proteobacteria in BE was the highest, Actinobacteriota in FE was the highest, Methylomirabilota in IS and BP were significantly higher than that of FE and BE (Figure 3; Supplementary Table S1). In addition, LEfSe analysis (threshold of 3) showed that there were differences in bacterial community composition between Control and different sites of PV panels. Specifically, Proteobacteria, Actinobacteriota, Gemmatimonadota, Bacteroidota and Patescibacteria were the most enriched in the Control. Chloroflexi, Gemmatimonadota, Cyanobacteria, Actinobacteriota and Methylomirabilota in IS. Unclassified_k__norank_d__Bacteria, Verrucomicrobiota, Bacteroidota, Actinobacteriota, Planctomycetota, Proteobacteria, Myxococcota and Proteobacteria in FE. Methylomirabilota in BP. Proteobacteria and Actinobacteriota in BE (Supplementary Figure S5).

Figure 3

At the phylum level, the top 10 dominant fungi with relative abundance in all sites accounted for more than 99% (76.16% Ascomycota, 7.80% Basidiomycota, 6.75% Unclassified_k__Fungi, 3.88% Mortierellomycota, 3.06% Chytridiomycota, 1.15% Glomeromycota, 0.56% Rozellomycota, 0.27% Monoblepharomycota, 0.12% Aphelidiomycota and 0.07% Zoopagomycota). PV panels had no significant effect on the relative abundance of fungal dominant phyla, and there was no significant difference between the relative abundance of dominant phyla at different sites of PV panels and Control (Figure 3; Supplementary Table S2). LEfSe analysis showed that Basidiomycota, Ascomycota, Mucoromycota and Glomeromycota were the most enriched in Control. Ascomycota in IS. Basidiomycota, Zoopagomycota, Ascomycota and Aphelidiomycota in FE. Ascomycota and Glomeromycota in BP, BE is Ascomycota (Supplementary Figure S6).

Relationships between biotic and abiotic properties and microbial communities

RDA analysis was performed to illustrate the relationship between bacterial and fungal community composition and soil properties. The first two RDA1 and RDA1 axes explained 61.62 and 12.82% of the variation in overall bacterial community composition, respectively. RDA results showed that EC, Margale and total AGB affected the composition of bacterial community (p = 0.004; p = 0.035; p = 0.025) (Figure 4; Supplementary Table S3). The first two RDA1 and RDA1 axes explained 29.11 and 14.59% of the variation in fungal community composition, respectively. RDA results showed that AN and AP significantly affected the composition of fungal communities (p = 0.034; p = 0.018) (Figure 4; Supplementary Table S3). SEM model was used to describe the effects of PV panels on plant diversity and soil microorganisms, and the interactions among plant diversity, soil properties and soil microorganisms. SEM model was used to describe the interaction between plant communities, soil properties and soil microorganisms (p = 0.003, χ2/df = 1.036, GFI =0.809, RMSEA = 0.031). SEM results showed that plant diversity directly affected soil bacterial diversity or indirectly affected soil bacterial diversity by affecting soil pH (p < 0.05) (Figure 5).

Figure 4

Figure 5

Spearman rank correlation analysis showed that Actinobacteriota was negatively correlated with TP and total AGB, but positively correlated with EC. Proteobacteria was negatively correlated with AK. Chloroflexi was positively correlated with total AGB. Acidobacteriota was negatively correlated with pH and EC, and positively correlated with TP, AN and Shnnon-Wiener. Firmicutes was positively correlated with Simpson diversity. Gemmatimonadota was positively correlated with pH. Bacteroidota was positively correlated with pH. Desulfobacterota is negatively correlated with pH and the positive correlation total AGB, and negatively correlated with Shnnon-Wiener and Simpson diversity (Figure 6).

Figure 6

In the dominant phylum of fungi, Ascomycota was negatively correlated with Shnnon-Wiener. Monoblepharomycota was positively correlated with AP and Shnnon-Wiener. Aphelidiomycota was positively correlated with Shnnon-Wiener and Simpson diversity. Zoopagomycota was negatively correlated with pH and positively correlated with TP (Figure 6).

Discussion

Effects of PV panels on vegetation communities and soil properties

Temperature and precipitation were the main moderators of grassland community composition and richness (Yang et al., 2011; Jonas et al., 2015). PV panels significantly increased the diversity of plant communities for the following reasons: on the one hand, grasses have shallow and fibrous roots, usually distributed in the soil surface (Mackie et al., 2019), while forbs and sedges have deep roots and are resource conservative (Yang et al., 2011). The PV panels reduced water collection and lower potential evaporation, which promoted soil water retention and improved the effective water content of deep soil. By affecting interspecific competition, it inhibited the growth of dominant PG and created ecological niches, and to a certain extent promoted the growth of subdominant PF, PS and ABH. On the other hand, numerous studies have demonstrated that temperature changes have different effects on different plant groups of PF, PS and PG (Peng et al., 2017; Ganjurjav et al., 2018). The absorption of solar radiation by PV panels effectively reduced the soil temperature, and the change of temperature effectively regulated the interspecific relationship and promoted the growth of different species. In addition, the adaptation or plasticity to shading varies greatly among plants (Semchenko et al., 2011). Faster growing and larger hybrid grasses have a disproportionate advantage in increasing their light capture capacity during shading through an increase in stem and leaf elongation and specific leaf area (SLA) (Ballare et al., 1994; Schwinning and Weiner, 1998) and resulted in competitive exclusion of smaller or slower-growing grass individuals (Hautier et al., 2009). In conclusion, PV panels effectively inhibited the growth of dominant PG by changing abiotic factors, promoted the growth of subdominant species, and then increased the diversity of grassland plant communities, and changed the composition and structure of plant communities.

Water availability in semi-arid regions limits plant growth and ecosystem productivity (Niu et al., 2008). Precipitation from January to July was reported to be the main cause of community biomass fluctuation in the same area (Bai et al., 2004). PV panels increased soil available water content, which not only directly promoted community photosynthesis and plant growth (Bai et al., 2010; Li et al., 2016), and indirectly increased aboveground biomass by increasing plant community diversity (Isbell et al., 2013). Numerous studies have demonstrated a positive correlation between biodiversity and biomass production (Isbell et al., 2013; Hautier et al., 2014; Hautier and Van der Plas, 2022). However, there were significant differences in total AGB among different sites of PV panels, and the higher soil water content at FE resulted in the most significant increase in total AGB. Total AGB under BP was lower than that at other sites of the PV panels, which may be attributed to, on the one hand, the shading of the PV panels reduced precipitation; on the other hand, the low light level and thick litter layer under BP inhibited plant growth (Seitz et al., 2016). It is worth noting that from the perspective of homogeneity, IS was least affected by PV panels in different sites under PV panels, compared with IS, the plant species diversity and total AGB of FE were significantly improved, and BP were significantly reduced, which may be that the PV panels were oblique arrangement, so that the soil moisture content of FE was significantly higher than other sites under the PV panels, and BP to precipitation and light inhibited the growth of plants.

Plants and soil are inseparable (Pugnaire et al., 2019), plant species replacement caused by climate change will change material input and decomposition by changing the form of plant litter production and root deposition, thus affecting soil properties (Bardgett et al., 2013: Davidson and Janssens, 2006; Kulmatiski et al., 2008; Glassman et al., 2018). However, our results show that changes in plant community composition did not significantly alter soil properties. The study of Lambert et al. (2021) also showed that the construction of PV panels did not reduce the C and N contents, nor did it reduce the soil chemical properties. The reasons may be that the PV panel laying time in the study site is short, and the soil properties have no obvious difference with the Control due to the hysteresis of the response to plant community changes. In addition, low light levels increased litter production (Seitz et al., 2016), and increased soil material input, resulting in slightly higher soil nutrient content at BP than at other sites under PV panels.

Effects of PV panels on soil microbial diversity and structure

Soil moisture strongly influenced microbial diversity (Orwin et al., 2015; Yang et al., 2019). Different sites under PV panels (especially FE) significantly increased soil available water content, and increased soil water availability promoted microbial growth. PV panels significantly increased the richness and diversity of soil bacterial communities, changed the composition of bacterial communities, and increased the richness of fungal communities. FE bacterial diversity increased significantly compared to IS, while BP decreased significantly, possibly due to the redistribution of precipitation by PV panels. NMDS analysis showed that the composition of bacterial and fungal communities of PV panels was similar to that of controls, which may be due to the delayed response of bacteria and fungi due to the short PV placement time.

Actinobacteriota and Proteobacteria have also been identified as the major bacterial communities in soils (Banerjee et al., 2018). In this study, Actinobacteriota and Proteobacteria were the most dominant phyla in all sites. Consistent with previous studies, Actinobacteriota plays a key role in desert grassland soils, followed by Proteobacteria, a more important phylum in typical grassland soils (Wang et al., 2017). It has been reported that Actinobacteriota is a group of eutrophic bacteria whose relative abundance increases with the size of soil C and N pools (Fierer et al., 2007; Li et al., 2014). Acidobacteria is an oligotrophic bacterium with poor tolerance and low temperature tolerance (Fierer et al., 2007; Yao et al., 2017). The relative abundance of Actinobacteriota at different sites under PV panels decreased significantly, while Acidobacteria increased significantly, possibly because the laying of PV panels reduced soil nutrition. The relative abundance of Proteobacteria increased significantly under the artificial interference of PV panels. Xiao et al. (2021) and Kolton et al. (2011) also reported that the soil disturbed by mining was significantly enriched in Proteobacteria. Proteobacteria is a fast-growing bacterium that can utilize a variety of carbon sources (i.e., syntrophic bacteria) (Yao et al., 2017). This may be related to the extensive biodegradative and metabolic properties of Proteobacteria bacteria and their ability to inhabit a variety of habitats (Sinkko et al., 2013). Proteobacteria may grow rapidly when available substrates are unstable (e.g., soil disturbance) (Zhang et al., 2016). LefSe analysis showed that there were significant differences in dominant phyla at different sites under PV panels, and the significantly different species in FE were significantly higher than those in other sites, which indicated that PV panels, especially FE, significantly changed the diversity and structure of soil bacteria.

The major fungal phyla in soils worldwide are Ascomycota and Basidiomycota (Tedersoo et al., 2014). This is consistent with our results that Ascomycota was the dominant phylum in different fungal communities, followed by Basidiomycota. Ascomycotion is an oligotrophic fungus that can tolerate harsh conditions (Chen et al., 2017), Basidiomycota is a eutrophic bacterium (Sterkenburg et al., 2015), The lack of resources may be the reason why the relative abundance of Ascomycota and Basidiomycota did not change under different sites of PV panels.

Biotic and abiotic characteristics in relation to microbial communities

Changes in precipitation and temperature affect soil microbial communities directly by altering their growth and activity (Bardgett et al., 2008) and also by altering their function on individual plants (e.g., root biomass production) (Orwin et al., 2010), or plant community composition indirectly affects soil microbial populations (Kardol et al., 2007; Luo et al., 2017; Ren et al., 2018). RDA, Spearman rank correlation and SEM showed that plant diversity was significantly correlated with the diversity and composition of bacterial and fungal communities. Plant diversity directly affected soil bacterial diversity or indirectly affected soil bacterial diversity by affecting soil pH. PV panels alter plant community composition, which changes the quality and/or quantity of root exudates and litter entering the soil, which in turn affects microbial community composition and activity (Bardgett and Wardle, 2010). Studies have shown that low-quality, refractory, slow-growing species determine fungus-dominated microbial community composition, while nutrient-rich, fast-growing species (such as miscellaneous grasses) determine bacteria-dominated microbial community composition (Orwin et al., 2010; de Vries et al., 2012; Pugnaire et al., 2019). PV significantly increased fast-growing species such as PF, which provided a good growth substrate for bacteria and drove the growth of bacterial communities. In addition, each plant has unique root exudate characteristics, and a richer plant community may have created a niche for the growth of different microorganisms (Zak et al., 2003; Bartelt-Ryser et al., 2005; Loranger-Merciris et al., 2006). Meanwhile, bacterial diversity and dominant phyla were significantly correlated with total AGB, which was consistent with the conclusion of Bardgett and McAlister (1999), suggesting that productivity may be a mechanism linking plant species richness and soil community structure.

Salinity (EC) was a major factor affecting bacterial community composition (Lozupone and Knight, 2007), which was consistent with our RDA results. Spearman rank correlation analysis showed that Actinobacteriota and Acidobacteriota were significantly correlated with EC. Previous studies have shown that Firmicutes, Bacteroidota and Proteobacteria have high salt tolerance, and their relative abundance is related to EC (Antunes et al., 2011; Geyer et al., 2014). With the increase of soil salinity, only a limited number of bacterial groups can withstand the enormous pressure of high salinity, resulting in changes in the relative abundance of microorganisms and community diversity (Rath et al., 2018). In addition, salinity is often correlated with OM inputs and soil pH (Wong et al., 2008; Rengasamy, 2010). SEM results showed that bacterial diversity was affected by soil pH. The correlation analysis shows that Desulfobacterota is positively correlated with pH. Since most bacterial taxa have relatively narrow growth tolerances, changes in soil pH may directly affect changes in bacterial communities (Lauber et al., 2008; Rousk et al., 2010). Proteobacteria was negatively correlated with AK. K is an essential element that regulates cellular pH in microorganisms and drives uptake of soil nutrients by regulating cellular osmotic pressure (Pan et al., 2020). Soil bacteria have been reported to influence K solubility and availability, which in turn influences the selection of specific bacteria associated with K concentrations (Pereira et al., 2014; Emmert et al., 2021). There is no general consensus on the effect of P on microbial community composition, as P has negative, neutral, or positive effects on soil microorganisms in terrestrial ecosystems (DeForest et al., 2012; Liu et al., 2012; Huang et al., 2016). RDA showed that the fungal community composition was affected by AP content. Correlation analysis showed that the relative abundance of Actinobacteriota, Acidobacteriota, Monoblepharomycota, and Zoopagomycota were significantly correlated with P content. This may be due to the decrease in soil pH by P application, which drives the change in the relative abundance of different bacteria (Jones et al., 2009). In addition, AN also significantly affected fungal community composition. Studies have shown that soil N content can directly affect soil microorganisms, or indirectly affect soil microorganisms by changing soil C availability, C/N ratio and soil pH (Treseder, 2008; Serna-Chavez et al., 2013). Taken together, our study provides further evidence that PV panels drive soil microbial changes by influencing abiotic and biotic conditions, and that soil (microenvironment) heterogeneity is a factor that maintains soil microbial community diversity.

Conclusion

In conclusion, our study found that PV panels significantly increased grassland plant community diversity by driving microclimate change. FE increased precipitation accumulation and plant diversity directly and indirectly changed the diversity of soil bacterial and fungal communities. PV panels decreased the relative abundance of Actinobacteriota, while increased that of Proteobacteria, Acidobacteriota, and Methylomirabilota. EC, Margalef’s richness and total AGB were the main factors affecting the composition of bacterial community, while AN and AP were the main factors affecting the composition of fungal community. Different microorganisms were affected by different soil properties, and the increase of plant diversity driven by abiotic factors was the main factor maintaining the diversity of soil microbial community and ecosystem function. The results provided important information for further understanding the effect of photovoltaic panel laying on grassland ecosystem function, and were of great significance for maintaining grassland ecosystem function.

Funding

This project was supported by the new energy field in Daqing of China “The open competition mechanism to select the best candidates” scientific and technological project (2021BD05) and the Special Fund Project for Returnees Studying Abroad of Heilongjiang Province (2021294).

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: Sequence Read Archive (SRA) of the NCBI database under the BioProjectID PRJNA895197.

Author contributions

ZyB, ZjB, and MW conceived and designed this study. ZyB drafted the original manuscript. MW provided very constructive suggestions for revisions. ZyB, MW, ZjB, SQ, MZ, LK, RS, and MW contributed to the sampling and data analysis. All authors read and approved the final manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1065899/full#supplementary-material

References

  • 1

    AdehE. H.SelkerJ. S.HigginsC. W. (2018). Remarkable agrivoltaic influence on soil moisture, micrometeorology and water-use efficiency. PLoS One13:e0203256. doi: 10.1371/journal.pone.0203256

  • 2

    AdlerP. B.Hille Ris LambersJ.KyriakidisP. C.GuanQ. F.LevineJ. M. (2006). Climate variability has a stabilizing effect on the coexistence of prairie grasses. Proc. Natl. Acad. Sci. U.S.A.103, 1279312798. doi: 10.1073/pnas.0600599103

  • 3

    AntunesA.AlamI.BajicV. B.StinglU. (2011). Genome sequence of Salinisphaera shabanensis, a gammaproteobacterium from the harsh, variable environment of the brine-seawater interface of the Shaban deep in the Red Sea. J. Bacteriol.193, 45554556. doi: 10.1128/JB.05459-11

  • 4

    ArmstrongA.WaldronS.WhitakerJ.OstleN. J. (2014). Wind farm and solar park effects on plant-soil carbon cycling: uncertain impacts of changes in ground-level microclimate. Glob. Change Biol.20, 16991706. doi: 10.1111/gcb.12437

  • 5

    BaiY. F.HanX. G.WuJ. G.ChenZ. Z.LiL. H. (2004). Ecosystem stability and compensatory effects in the Inner Mongolia grassland. Nature431, 181184. doi: 10.1038/nature02850

  • 6

    BaiY. F.WuJ. G.ClarkC. M.NaeemS.PanQ. M.HuangJ. H.et al. (2010). Tradeoffs and thresholds in the effects of nitrogen addition on biodiversity and ecosystem functioning: evidence from inner Mongolia grasslands. Glob. Chang. Biol.16, 358372. doi: 10.1111/j.1365-2486.2009.01950.x

  • 7

    BallareC. L.ScopelA. L.JordanE. T.VierstraR. D. (1994). Signaling among neighboring plants and the development of size inequalities in plant populations. Proc. Natl. Acad. Sci. U. S. A.91, 1009410098. doi: 10.1073/pnas.91.21.10094

  • 8

    BanerjeeS.SchlaeppiK.van der HeijdenM. G. (2018). Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol.16, 567576. doi: 10.1038/s41579-018-0024-1

  • 9

    BaoS. D. (2008). Soil and Agricultural Chemistry Analysis. (Beijing: Chinese Agriculture Press).

  • 10

    BardgettR. D.FreemanC.OstleN. J. (2008). Microbial contributions to climate change through carbon cycle feedbacks. ISME J.2, 805814. doi: 10.1038/ismej.2008.58

  • 11

    BardgettR. D.ManningP.MorrienE.De VriesF. T. (2013). Hierarchical responses of plant-soil interactions to climate change: consequences for the global carbon cycle. J. Ecol.101, 334343. doi: 10.1111/1365-2745.12043

  • 12

    BardgettR. D.McAlisterE. (1999). The measurement of soil fungal: bacterial biomass ratios as an indicator of ecosystem self-regulation in temperate meadow grasslands. Biol. Fertil.29, 282290. doi: 10.1007/s003740050554

  • 13

    BardgettR. D.WardleD. A. (2010). Aboveground Belowground Linkages: Biotic Interactions, Ecosystem Processes, and Global Change. (Oxford, UK: Oxford University Press).

  • 14

    Bartelt-RyserJ.JoshiJ.SchmidB.BrandlH.BalserT. (2005). Soil feedbacks of plant diversity on soil microbial communities and subsequent plant growth. Perspect. Plant Ecol. Evol. Syst.7, 2749. doi: 10.1016/j.ppees.2004.11.002

  • 15

    ChenY. L.XuT. L.VeresoglouS. D.HuH. W.HaoZ. P.HuY. J.et al. (2017). Plant diversity represents the prevalent determinant of soil fungal community structure across temperate grasslands in northern China. Soil Biol. Biochem.110, 1221. doi: 10.1016/j.soilbio.2017.02.015

  • 16

    ChenH.ZhuQ. A.PengC. H.WuN.WangY. F.FangX. Q.et al. (2013). The impacts of climate change and human activities on biogeochemical cycles on the Qinghai-Tibetan plateau. Glob. Chang. Biol.19, 29402955. doi: 10.1111/gcb.12277

  • 17

    ChoiC. S.CagleA. E.MacknickJ.BloomD. E.CaplanJ. S.RaviS. (2020). Effects of revegetation on soil physical and chemical properties in solar photovoltaic infrastructure. Front. Environ. Sci.8:140. doi: 10.3389/fenvs.2020.00140

  • 18

    ClelandE. E.SmithM. D.AndelmanS. J.BowlesC.CarneyK. M.Horner-DevineM. C.et al. (2004). Invasion in space and time: non-native species richness and relative abundance respond to interannual variation in productivity and diversity. Ecol. Lett.7, 947957. doi: 10.1111/j.1461-0248.2004.00655.x

  • 19

    CrowtherT. W.Todd-BrownK. E. O.RoweC. W.WiederW. R.CareyJ. C.MachmullerM. B.et al. (2016). Quantifying global soil carbon losses in response to warming. Nature540, 104108. doi: 10.1038/nature20150

  • 20

    DavidsonE. A.JanssensI. A. (2006). Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature440, 165173. doi: 10.1038/nature04514

  • 21

    de VriesF. T.ManningP.TallowinJ. R. B.MortimerS. R.PilgrimE. S.HarrisonK. A.et al. (2012). Abiotic drivers and plant traits explain landscape-scale patterns in soil microbial communities. Ecol. Lett.15, 12301239. doi: 10.1111/j.1461-0248.2012.01844.x

  • 22

    DeForestJ. L.SmemoK. A.BurkeD. J.ElliottH. L.BeckerJ. C. (2012). Soil microbial responses to elevated phosphorus and pH in acidic temperate deciduous forests. Biogeochemistry109, 189202. doi: 10.1007/s10533-011-9619-6

  • 23

    DingM.XuZ. C.WangW. S.WangX. L.SongY. T.ChenD. Z. (2016). A review on China′s large-scale PV integration: Progress, challenges and recommendations. Renew. Sust. Energ. Rev.53, 639652. doi: 10.1016/j.rser.2015.09.009

  • 24

    DuprazC.MarrouH.TalbotG.DufourL.NogierA.FerardY. (2011). Combining solar photovoltaic panels and food crops for optimising land use: towards new agrivoltaic schemes. Renew. Energy36, 27252732. doi: 10.1016/j.renene.2011.03.005

  • 25

    EdgarR. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods10, 996998. doi: 10.1038/nmeth.2604

  • 26

    ElamriY.ChevironB.MangeA.DejeanC.LironF.BelaudG. (2018). Rain concentration and sheltering effect of solar panels on cultivated plots. Hydrol. Earth Syst. Sci.22, 12851298. doi: 10.5194/hess-22-1285-2018

  • 27

    EmmertE. A. B.GeletaS. B.RoseC. M.Seho-AhiableG. E.HawkinsA. E.BakerK. T.et al. (2021). Effect of land use changes on soil microbial enzymatic activity and soil microbial community composition on Maryland’s eastern shore. Appl. Soil Ecol.161:103824. doi: 10.1016/j.apsoil.2020.103824

  • 28

    Fernandez-GoingB. M.HarrisonS. P.AnackerB. L.SaffordH. D. (2013). Climate interacts with soil to produce beta diversity in Californian plant communities. Ecology94, 20072018. doi: 10.1890/12-2011.1

  • 29

    FiererN.BradfordM. A.JacksonR. B. (2007). Toward an ecological classification of soil bacteria. Ecology88, 13541364. doi: 10.1890/05-1839

  • 30

    FthenakisV.KimH. C. (2009). Land use and electricity generation: a lifecycle analysis. Renew. Sust. Energ. Rev.13, 14651474. doi: 10.1016/j.rser.2008.09.017

  • 31

    GanjurjavH.GornishE. S.HuG.WanY.LiY.DanjiuL.et al. (2018). Temperature leads to annual changes of plant community composition in alpine grasslands on the Qinghai-Tibetan plateau. Environ. Monit. Assess.190:585. doi: 10.1007/s10661-018-6964-0

  • 32

    GeyerK. M.AltrichterA. E.Takacs-VesbachC. D.Van HornD. J.GooseffM. N.BarrettJ. E. (2014). Bacterial community composition of divergent soil habitats in a polar desert. FEMS Microbiol. Ecol.89, 490494. doi: 10.1111/1574-6941.12306

  • 33

    GlassmanS. I.WeiheC.LiJ. H.AlbrightM. B. N.LoobyC. I.MartinyA. C.et al. (2018). Decomposition responses to climate depend on microbial community composition. Proc. Natl. Acad. Sci. U. S. A.115, 1199411999. doi: 10.1073/pnas.1811269115

  • 34

    GottfriedM.PauliH.FutschikA.AkhalkatsiM.BarancokP.AlonsoJ. L. B.et al. (2012). Continent-wide response of mountain vegetation to climate change. Nat. Clim. Chang.2, 111115. doi: 10.1038/NCLIMATE1329

  • 35

    HautierY.NiklausP. A.HectorA. (2009). Competition for light causes plant biodiversity loss after eutrophication. Science324, 636638. doi: 10.1126/science.1169640

  • 36

    HautierY.SeabloomE. W.BorerE. T.AdlerP. B.HarpoleW. S.HillebrandH.et al. (2014). Eutrophication weakens stabilizing effects of diversity in natural grasslands. Nature508, 521525. doi: 10.1038/nature13014

  • 37

    HautierY.Van der PlasF. (2022). The Ecological and Societal Consequences of Biodiversity Loss. (London: Wiley-ISTE).

  • 38

    HuangJ. S.HuB.QiK. B.ChenW. J.PangX. Y.BaoW. K.et al. (2016). Effects of phosphorus addition on soil microbial biomass and community composition in a subalpine spruce plantation. Eur. J. Soil Biol.72, 3541. doi: 10.1016/j.ejsobi.2015.12.007

  • 39

    IsbellF.ReichP. B.TilmanD.HobbieS. E.PolaskyS.BinderS. (2013). Nutrient enrichment, biodiversity loss, and consequent declines in ecosystem productivity. Proc. Natl. Acad. Sci. U. S. A.110, 1191111916. doi: 10.1073/pnas.1310880110

  • 40

    JonasJ. L.BuhlD. A.SymstadA. J. (2015). Impacts of weather on long-term patterns of plant richness and diversity vary with location and management. Ecology96, 24172432. doi: 10.1890/14-1989.1

  • 41

    JonesR. T.RobesonM. S.LauberC. L.HamadyM.KnightR.FiererN. (2009). A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses. ISME J.3, 442453. doi: 10.1038/ismej.2008.127

  • 42

    JoshiP.AroraS. (2017). Maximum power point tracking methodologies for solar PV systems—a review. Renew. Sust. Energ. Rev.70, 11541177. doi: 10.1016/j.rser.2016.12.019

  • 43

    KardolP.CornipsN. J.MoniqueM. L.van KempenJ. M.Bakx-SchotmanT.van der PuttenW. H. (2007). Microbemediated plant-soil feedback causes historical contingency effects in plant community assembly. Ecol. Monogr.77, 147162. doi: 10.1890/06-0502

  • 44

    KoltonM.HarelY. M.PasternakZ.GraberE. R.EladY.CytrynE. (2011). Impact of biochar application to soil on the root-associated bacterial community structure of fully developed greenhouse pepper plants. Appl. Environ. Microbiol.77, 49244930. doi: 10.1128/AEM.00148-11

  • 45

    KulmatiskiA.BeardK. H.StevensJ. R.CobboldS. M. (2008). Plant-soil feedbacks: a meta-analytical review. Ecol. Lett.11, 980992. doi: 10.1111/j.1461-0248.2008.01209.x

  • 46

    LambertQ.BischoffA.CueffS.CluchierA.GrosR. (2021). Effects of solar park construction and solar panels on soil quality, microclimate, CO2 effluxes and vegetation under Mediterranean climate. Land Degrad. Dev.32, 51905202. doi: 10.1002/ldr.4101

  • 47

    LangleyJ. A.HungateB. A. (2014). Plant community feedbacks and long-term ecosystem responses to multi-factored global change. AoB Plants6:12. doi: 10.1093/aobpla/plu035

  • 48

    LauberC. L.StricklandM. S.BradfordM. A.FiererN. (2008). The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol. Biochem.40, 24072415. doi: 10.1016/j.soilbio.2008.05.021

  • 49

    LiX. R.JiaX. H.DongG. R. (2006). Influence of desertification on vegetation pattern variations in the cold semi-arid steppes of Qinghai-Tibet plateau. North-West China. J. Arid Environ.64, 505522. doi: 10.1016/j.jaridenv.2005.06.011

  • 50

    LiY.KalnayE.MotesharreiS.RivasJ.KucharskiF.Kirk-DavidoffD.et al. (2018). Climate model shows large-scale wind and solar farms in the Sahara increase rain and vegetation. Science361, 10191022. doi: 10.1126/science.aar5629

  • 51

    LiW. H.LiX.ZhaoY. J.ZhengS. X.BaiY. F. (2018). Ecosystem structure, functioning and stability under climate change and grazing in grasslands: current status and future prospects. Curr. Opin. Environ. Sust.33, 124135. doi: 10.1016/j.cosust.2018.05.008

  • 52

    LiY.NiuS. L.YuG. R. (2016). Aggravated phosphorus limitation on biomass production under increasing nitrogen loading: a meta-analysis. Glob. Chang. Biol.22, 934943. doi: 10.1111/gcb.13125

  • 53

    LiC. H.YanK.TangL. S.JiaZ. J.LiY. (2014). Change in deep soil microbial communities due to long-term fertilization. Soil Biol. Biochem.75, 264272. doi: 10.1016/j.soilbio.2014.04.023

  • 54

    LiuL.GundersenP.ZhangT.MoJ. M. (2012). Effects of phosphorus addition on soil microbial biomass and community composition in three forest types in tropical China. Soil Biol. Biochem.44, 3138. doi: 10.1016/j.soilbio.2011.08.017

  • 55

    LiuY.ZhangR. Q.HuangZ.ChengZ.Lopez-VicenteM.MaX.et al. (2019). Solar photovoltaic panels significantly promote vegetation recovery by modifying the soil surface microhabitats in arid sandy ecosystem. Land Degrad. Dev.30, 21772186. doi: 10.1002/ldr.3408

  • 56

    Loranger-MercirisG.BarthesL.GastineA.LeadleyP. (2006). Rapid effects of plant species diversity and identity on soil microbial communities in experimental grassland ecosystems. Soil Biol. Biochem.38, 23362343. doi: 10.1016/j.soilbio.2006.02.009

  • 57

    LozuponeC. A.KnightR. (2007). Global patterns in bacterial diversity. Proc. Natl. Acad. Sci. U. S. A.104, 1143611440. doi: 10.1073/pnas.0611525104

  • 58

    LuR. (1999). Soil Agricultural Chemical Analysis method. (Beijing: Chinese Agriculture and Sciences Press).

  • 59

    LuoY. Q.JiangL. F.NiuS. L.ZhouX. H. (2017). Nonlinear responses of land ecosystems to variation in precipitation. New Phytol.214, 57. doi: 10.1111/nph.14476

  • 60

    Mac KayD. J. C. (2013). Solar energy in the context of energy use, energy transportation and energy storage. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci.371:20110431. doi: 10.1098/rsta.2011.0431

  • 61

    MackieK. A.ZeiterM.BloorJ. M. G.StampfliA. (2019). Plant functional groups mediate drought resistance and recovery in a multisite grassland experiment. J. Ecol.107, 937949. doi: 10.1111/1365-2745.13102

  • 62

    MagocT.SalzbergS. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics27, 29572963. doi: 10.1093/bioinformatics/btr507

  • 63

    MarrouH.DufourL.WeryJ. (2013). How does a shelter of solar panels influence water flows in a soil-crop system?Eur. J. Agron.50, 3851. doi: 10.1016/j.eja.2013.05.004

  • 64

    MartensS. N.BreshearsD. D.MeyerC. W. (2000). Spatial distributions of understory light along the grassland/forest continuum: effects of cover, height, and spatial pattern of tree canopies. Ecol. Model.126, 7993. doi: 10.1016/S0304-3800(99)00188-X

  • 65

    NiuS. L.WuM. Y.HanY.XiaJ. Y.LiL. H.WanS. Q. (2008). Water-mediated responses of ecosystem carbon fluxes to climatic change in a temperate steppe. New Phytol.177, 209219. doi: 10.1111/j.1469-8137.2007.02237.x

  • 66

    OrwinK. H.BucklandS. M.JohnsonD.TurnerB. L.SmartS.OakleyS.et al. (2010). Linkages of plant traits to soil properties and the functioning of temperate grassland. J. Ecol.98, 10741083. doi: 10.1111/j.1365-2745.2010.01679.x

  • 67

    OrwinK. H.DickieI. A.WoodJ. R.BonnerK. I.HoldawayR. J. (2015). Soil microbial community structure explains the resistance of respiration to a dry–rewet cycle, but not soil functioning under static conditions. Funct. Ecol.30, 14301439. doi: 10.1111/1365-2435.12610

  • 68

    PanX. M.ZhangS. R.ZhongQ. M.GongG. S.WangG. Y.GuoX.et al. (2020). Effects of soil chemical properties and fractions of Pb, cd, and Zn on bacterial and fungal communities. Sci. Total Environ.715:136904. doi: 10.1016/j.scitotenv.2020.136904

  • 69

    PengF.XueX.XuM. H.YouQ. G.JianG.MaS. X. (2017). Warming-induced shift towards forbs and grasses and its relation to the carbon sequestration in an alpine meadow. Environ. Res. Lett.12:044010. doi: 10.1088/1748-9326/aa6508

  • 70

    PereiraL. B.VicentiniR.OttoboniL. M. M.MoustafaA. (2014). Changes in the bacterial community of soil from a neutral mine drainage channel. PLoS One9:e96605. doi: 10.1371/journal.pone.0096605

  • 71

    PiaoS. L.TanK.NanH. J.CiaisP.FangJ. Y.WangT.et al. (2012). Impacts of climate and CO2 changes on the vegetation growth and carbon balance of Qinghai-Tibetan grasslands over the past five decades. Glob Planet Change.98-99, 7380. doi: 10.1016/j.gloplacha.2012.08.009

  • 72

    PugnaireF. I.MorilloJ. A.PenuelasJ.ReichP. B.BardgettR. D.GaxiolaA.et al. (2019). Climate change effects on plant-soil feedbacks and consequences for biodiversity and functioning of terrestrial ecosystems. Sci. Adv.5:eaaz1834. doi: 10.1126/sciadv.aaz1834

  • 73

    QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res.41, D590D596. doi: 10.1093/nar/gks1219

  • 74

    RathK. M.FiererN.MurphyD. V.RouskJ. (2018). Linking bacterial community composition to soil salinity along environmental gradients. ISME J.13, 836846. doi: 10.1038/s41396-018-0313-8

  • 75

    RenC. J.ChenJ.LuX. J.DoughtyR.ZhaoF. Z.ZhongZ. K.et al. (2018). Responses of soil total microbial biomass and community comsites to rainfall declines. Soil Biol. Biochem.116, 410. doi: 10.1016/j.soilbio.2017.09.028

  • 76

    REN21. (2019). Renewables 2019 Global Status Report. REN21 Secretariat, Paris.

  • 77

    RengasamyP. (2010). Soil processes affecting crop production in saltaffected soils. Funct. Plant Biol.37, 613620. doi: 10.1071/FP09249

  • 78

    RouskJ.BaathE.BrookesP. C.LauberC. L.LozuponeC. A.CaporasoJ. G.et al. (2010). Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J.4, 13401351. doi: 10.1038/ismej.2010.58

  • 79

    RudgersJ. A.KivlinS. N.WhitneyK. D.PriceM. V.WaserN. M.HarteJ. (2014). Responses of high-altitude graminoids and soil fungi to 20 years of experimental warming. Ecology95, 19181928. doi: 10.1890/13-1454.1

  • 80

    SchwinningS.WeinerJ. (1998). Mechanisms determining the degree of size asymmetry in competition among plants. Oecologia113, 447455. doi: 10.1007/s004420050397

  • 81

    SeitzS.GoebesP.SongZ. S.BruelheideH.HardtleW.KuhnP.et al. (2016). Tree species and functional traits but not species richness affect interrill erosion processes in young subtropical forests. Soil2, 4961. doi: 10.5194/soil-2-49-2016

  • 82

    SemchenkoM.LepikM.GotzenbergerL.ZobelK. (2011). Positive effect of shade on plant growth: amelioration of stress or active regulation of growth rate?J. Ecol.100, 459466. doi: 10.1111/j.1365-2745.2011.01936.x

  • 83

    Serna-ChavezH. M.FiererN.BodegomP. M. (2013). Global drivers and patterns of microbial abundance in soil. Glob. Ecol. Biogeogr.22, 11621172. doi: 10.1111/geb.12070

  • 84

    ShiZ.SherryR.XuX.HararukO.SouzaL.JiangL. F.et al. (2015). Evidence for long-term shift in plant community composition under decadal experimental warming. J. Ecol.103, 11311140. doi: 10.1111/1365-2745.12449

  • 85

    ShiY.ZhangK. P.LiQ.LiuX.HeJ. S.ChuH. Y. (2020). Interannual climate variability and altered precipitation influence the soil microbial community structure in a Tibetan plateau grassland. Sci. Total Environ.714:136794. doi: 10.1016/j.scitotenv.2020.136794

  • 86

    SinkkoH.LukkariK.SihvonenL. M.SivonenK.LeivuoriM.RantanenM.et al. (2013). Bacteria contribute to sediment nutrient release and reflect progressed eutrophication-driven hypoxia in an organic-rich continental sea. PLoS One8:e67061. doi: 10.1371/journal.pone.0067061

  • 87

    SparksD. L.PageA.HelmkeP.LoeppertR.SoltanpourP.TabatabaiM.et al. (1996). Methods of Soil Analysis. Part 3-Chemical Methods. Madison: Soil Science Society of America Inc..

  • 88

    SterkenburgE.BahrA.DurlingM. B.ClemmensenK. E.LindahlB. D. (2015). Changes in fungal communities along a boreal forest soil fertility gradient. New Phytol.207, 11451158. doi: 10.1111/nph.13426

  • 89

    TedersooL.BahramM.PolmeS.KoljalgU.YorouN. S.WijesunderaR.et al. (2014). Global diversity and geography of soil fungi. Science346:1256688. doi: 10.1126/science.1256688

  • 90

    TojuH.PeayK. G.YamamichiM.NarisawaK.HirumaK.NaitoK.et al. (2018). Core microbiomes for sustainable agroecosystems. Nature4, 247257. doi: 10.1038/s41477-018-0139-4

  • 91

    TresederK. K. (2008). Nitrogen additions and microbial biomass: a meta-analysis of ecosystem studies. Ecol. Lett.11, 11111120. doi: 10.1111/j.1461-0248.2008.01230.x

  • 92

    WaltherG. R.PostE.ConveyP.MenzelA.ParmesanC.BeebeeT. J. C.et al. (2002). Ecological responses to recent climate change. Nature416, 389395. doi: 10.1038/416389a

  • 93

    WangQ.GarrityG. M.TiedjeJ.ColeJ. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol.73, 52615267. doi: 10.1128/AEM.00062-07

  • 94

    WangX. B.LuX. T.YaoJ.WangZ. W.DengY.ChengW. X.et al. (2017). Habitat-specific patterns and drivers of bacterial Beta-diversity in China’s drylands. ISME J.11, 13451358. doi: 10.1038/ismej.2017.11

  • 95

    WenL.DongS. K.LiY. Y.WangX. X.LiX. Y.ShiJ. J.et al. (2013). The impact of land degradation on the C pools in alpine grasslands of the Qinghai-Tibet plateau. Plant Soil368, 329340. doi: 10.1007/s11104-012-1500-4

  • 96

    WongV. N. L.DalalR. C.GreeneR. S. B. (2008). Salinity and sodicity effects on respiration and microbial biomass of soil. Biol. Fert. Soils44, 943953. doi: 10.1007/s00374-008-0279-1

  • 97

    WuZ. Q.NjokeM. L.TianG. N.FengJ. (2021). Challenges of investment and financing for developing photovoltaic power generation in Cameroon, and the countermeasures. J. Clean. Prod.299:126910. doi: 10.1016/j.jclepro.2021.126910

  • 98

    XiaoE. Z.WangY. Q.XiaoT. F.SunW. M.DengJ. M.JiangS. M.et al. (2021). Microbial community responses to land-use types and its ecological roles in mining area. Sci. Total Environ.775:145753. doi: 10.1016/j.scitotenv.2021.145753

  • 99

    YangH. J.LiY.WuM. Y.ZhangZ.LiL. H.WanS. Q. (2011). Plant community responses to nitrogen addition and increased precipitation: the importance of water availability and species traits. Glob. Chang. Biol.17, 29362944. doi: 10.1111/j.1365-2486.2011.02423.x

  • 100

    YangF.ZhangD. D.WuJ. J.ChenQ.LongC. Y.LiY. H.et al. (2019). Anti-seasonal submergence dominates the structure and composition of prokaryotic communities in the riparian zone of the three gorges reservoir. China. Sci Total Environ.663, 662672. doi: 10.1016/j.scitotenv.2019.01.357

  • 101

    YaoF.YangS.WangZ. R.WangX.YeJ.WangX. G.et al. (2017). Microbial taxa distribution is associated with ecological trophic cascades along an elevation gradient. Front. Microbiol.8:2071. doi: 10.3389/fmicb.2017.02071

  • 102

    YueS. J.GuoM. J.ZouP. H.WuW.ZhouX. D. (2021). Effects of photovoltaic panels on soil temperature and moisture in desert areas. Environ. Sci. Pollut. R.28, 1750617518. doi: 10.1007/s11356-020-11742-8

  • 103

    ZakD. R.HolmesW. E.WhiteD. C.PeacockA. D.TilmanD. (2003). Plant diversity, soil microbial communities, and ecosystem function: are there any links?Ecology84, 20422050. doi: 10.1890/02-0433

  • 104

    ZhangC.LiuG. B.XueS.WangG. L. (2016). Soil bacterial community dynamics reflect changes in plant community and soil properties during the secondary succession of abandoned farmland in the loess plateau. Soil Biol. Biochem.97, 4049. doi: 10.1016/j.soilbio.2016.02.013

Summary

Keywords

photovoltaic panel, grassland, aboveground biomass, plant community composition, soil microbial diversity

Citation

Bai Z, Jia A, Bai Z, Qu S, Zhang M, Kong L, Sun R and Wang M (2022) Photovoltaic panels have altered grassland plant biodiversity and soil microbial diversity. Front. Microbiol. 13:1065899. doi: 10.3389/fmicb.2022.1065899

Received

10 October 2022

Accepted

21 November 2022

Published

15 December 2022

Volume

13 - 2022

Edited by

Bin Huang, Tobacco Research Institute (CAAS), China

Reviewed by

Tao Zhang, Northeast Normal University, China; Gang Fu, Institute of Geographic Sciences and Natural Resources Research (CAS), China

Updates

Copyright

*Correspondence: Mingjun Wang,

This article was submitted to Terrestrial Microbiology, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics