Description of study area
Field investigations were conducted along a 1,300-km transect in the southeastern Tibetan Plateau during July and August 2020. The study area spanned latitudes of 29°20'N to 29°53'N and longitudes of 91°00'E to 95°46'E, with elevations ranging from 2,647 to 4,898 meters above sea level. Climatic data were obtained from the Xizang Meteorological Stations via the China Meteorological Network, covering the period from 1990 to 2020. Based on these data, we calculated the mean annual temperature (MAT), mean annual precipitation (MAP), and mean annual sunshine hours for each sampling site. Across all sites, the average MAT was 5.03°C, the average MAP was 647.78 mm, and the mean annual sunshine duration totaled 2,556.52 hours.
The Qinghai-Tibet Plateau is dominated by alpine meadows and alpine steppes, with temperate grasslands occurring at lower elevations. Due to the harsh climatic conditions, the plant growing season is notably short. Dominant species in the study area include Kobresia spp., Potentilla spp., and Poa spp.
Field transect survey and sampling methodology
We established 34 sampling sites along the transect, with geographic coordinates (latitude, longitude) and elevation recorded at each site using a GPS device (eTrex Venture, Garmin, USA). To ensure spatial independence, sites were spaced more than 50 km apart. Within each site, we systematically placed five quadrats (1m×1m) with a 100m interval between different quadrats. We conducted comprehensive vegetation surveys in each quadrat, measuring community and species-level parameters. We recorded the plant density, mean height, coverage, and aboveground biomass of each plant species. Based on their widespread distribution across sampling sites, we selected four plant species for stoichiometric analysis: Potentilla saundersiana (PS), Carex moorcroftii (CM), Kobresia humilis (KH), and Leontopodium leontopodioides (LL).
When the four target species did not co-occur at a site, we collected only the available species, and stoichiometric traits were analyzed for each species across sites. Overall, we gathered data from 14 sites for PS, 17 sites for CM, 13 sites for KH, and four sites for LL (Fig. 1). At each study site, we collected 30 intact specimens per target species for stoichiometric characterization. This sampling design provided three biological replicates per species per site, with each replicate consisting of ten individual plants. Following plant collection, we extracted soil samples (0–15 cm depth) from the rhizosphere of each species. Soil from 10 individuals was composited per replicate. In total, we analyzed three elements across 138 individuals, yielding 433 samples (site locations and climatic data are detailed in Table S1).
Laboratory analysis
All plant samples were oven-dried at 65°C for 48 hr until reaching constant weight, then finely ground using a ball mill (Retsch MM 400, Germany). Soil samples were air-dried, homogenized, and sieved (2.0 mm mesh) to remove coarse debris. The N concentrations in leaves, stems, roots, and soil were quantified via the Kjeldahl digestion-distillation method. For total P concentrations, plant and soil samples underwent persulfate-mediated oxidative digestion, followed by ammonium molybdate spectrophotometric detection. The organic carbon contents were determined using the Springer-Klee wet oxidation method, involving dichromate (K₂Cr₂O₇)-sulfuric acid (H₂SO₄) digestion under controlled heating.
Data calculation and statistical analysis
Data normality was assessed using the Kolmogorov-Smirnov test before conducting statistical analyses. When required, log-transformation was applied to normalize the data distribution. We initially conducted Pearson's correlation analysis to examine the influence of elevation gradients on the stoichiometric traits of four plant species. Subsequently, a two-way ANOVA was performed to assess the effects of species and organs on these stoichiometric traits. Additionally, multiple comparisons of means were performed using Tukey’s HSD test for roots and soil stoichiometric traits, while paired t-tests were applied to assess stoichiometric traits of shoots for KH and LL, and stem and leaves for CM and PS. Significance levels were set at P<0.05, P<0.01, and P<0.001. We performed the principal component analysis (PCA) to quantify the biogeochemical niche of four plant species based on three elemental traits. The first three principal components (PC1, PC2, and PC3) were extracted using the stats package in R. Biogeochemical niche differentiation was further assessed by calculating hypervolumes from the PCA scores using the hypervolume package. To examine how climatic factors (MAT and MAP), soil properties (C, N, and P concentrations), and biological factors (species and organ type) influenced stoichiometric variation along the elevation gradient, we employed both Generalized Linear Mixed Models (GLMMs) and Linear Mixed Effects Models (LMEMs), with elevation included as a random effect. Because of missing soil phosphorus data for LL, the GLMM and LMEM analyses were restricted to PS, CM, and KH. All statistical analyses were performed in R version 4.1.1 (R Core Team, 2015), using the lmer package for LMEMs and glmm.hp for GLMMs. All the graphs were generated using OriginPro 2024 (Education Edition) or ArcGIS (version 10.8, ESRI, Redlands, CA).