Abstract
The Qinghai-Xizang Plateau has abundant sunlight resources. Making rational use of this advantage can effectively increase the forage yield and nutrient accumulation of mixed grassland sown with gramineous and leguminous plants. External fertilization is a key measure to further enhance the productivity and quality of the mixed grassland. This study focused on the Avena sativa+Pisum sativum mixed sowing system and set up four treatment groups: single application of nitrogen fertilizer (A1), single application of phosphorus fertilizer (A2), nitrogen and phosphorus combined application (A3), and no fertilization control (CK). The study evaluated the effects of different fertilization treatments on the spectral indices and leaf morphology of the oat-pea mixed sowing leaves and their impact on forage quality. The results showed that under the nitrogen and phosphorus combined application treatment, the total chlorophyll content of oat plants was 15.14% higher than that of the CK; the soil-plant analysis development value was 26.44% higher than that of the CK, 17.05% higher than A1, and 26.85% higher than A2. In terms of leaf morphology, the nitrogen and phosphorus combined application significantly increased leaf area, while single application of nitrogen fertilizer made leaf length, leaf width, and perimeter higher than that of the CK. In terms of nutritional quality, the nitrogen and phosphorus combined application significantly increased total phosphorus, total nitrogen, crude protein, soluble sugar, and total carbon content, and simultaneously reduced neutral detergent fiber and acid detergent fiber by 7.78% and 25.25%, respectively, compared to the CK. In summary, the combined application of nitrogen and phosphorus significantly affected the spectral indices, morphology, and plant nutrient indicators of the leaves. This fertilization strategy and mixed sowing combination can be applied in the experimental area to further enhance the yield and quality of artificial grassland.
Similar content being viewed by others
Plant community restoration time is a key factor affecting plant nutrient and species richness
Productivity, soil fertility and enzyme activity of mixed forage grasslands improved by alfalfa and nitrogen addition in Horqin Sandy Land, China
Effects of different afforestation years on soil moisture and nutrient content in Maxian Mountain of the Loess Plateau
Introduction
The Qinghai-Xizang Plateau is located in the central part of Asia, with an average altitude of approximately 4,000 m. It is one of the sources of the Yangtze River, the Yellow River, and the Lancang River1. Its complex geographical features have not only shaped a unique ecological environment for the plateau grasslands but also played an important role in regulating the global climate. For a long time, grassland animal husbandry has been the basic industry in the Sanjiangyuan region and the fundamental means of survival for the herdsmen2. However, with the intensification of climate change and population growth, traditional grassland animal husbandry is facing a severe challenge of unbalanced forage supply. The long winter with a grass shortage period of up to seven or eight months leads to a significant decline in the nutritional value of the forage. The available land area in the Sanjiangyuan region is small, and how to utilize the limited land area to balance ecological protection and economic development has become one of the urgent problems to be solved3,4. Therefore, the establishment of high-quality and high-yield artificial grasslands proposed in this context becomes the key to solving this problem.
Exogenous nutrient addition is a key management measure for enhancing the productivity of artificial grasslands. The rational application of nitrogen fertilizer and phosphorus fertilizer has a significant impact on the growth, quality and yield of forage grass5,6. Studies have shown that when using nitrogen fertilizer, the yield of alfalfa will increase by 8% to 25%, while when using phosphorus fertilizer, it will increase by 122% to 172%7[,8. There is a large amount of phosphorus in the soil, but it usually exists in the form of insoluble mineral phosphates or is combined with organic matter, which makes it difficult for plants to directly absorb or utilize. This becomes a major limiting factor for plant growth and yield9. In agricultural and livestock production, phosphate fertilizers are usually used to meet the phosphorus requirements of plants. Due to the tendency of chemical phosphate fertilizers to highly combine with ions such as calcium, iron and aluminum, they will deposit in the form of inorganic phosphates in the soil, so the effective utilization efficiency of plants is lower than 25%10,11. The increase in nitrogen deposition in the atmosphere will lead to soil acidification, which will affect soil respiration and ultimately affect the photosynthesis of plant leaves12,13. Nitrogen and phosphorus in the leaves are powerful factors for photosynthesis. When the nitrogen intake increases and chlorophyll increases, the concentration and activity of enzymes involved in photosynthesis will increase14,15, thereby improving the photosynthetic capacity of plant leaves16,17.
The competition among different species and the proportion of leguminous plants in the forage have significant differences in their ability to fix nitrogen18,19, so the plant nutrients and community yield may be affected by various leguminous forage species and their proportions in the legume-cereal mixed planting combinations. The net photosynthetic rate, water use efficiency (WUE), and leaf nitrogen content of the legume-mixed planting combinations are all higher than those of monoculture20,21,22. Yellow clover (Medicago sativa ssp.) was added to the forage to increase the content of calcium (Ca), potassium (K), and magnesium (Mg), which increased the yield of the forage and also improved its nutritional components23,24. However, due to geographical limitations, there are low nutrient content soil conditions in the Sanjiangyuan region of the Qinghai-Tibet Plateau. Therefore, some exogenous fertilizers must be added to provide nutritional support for the growth and development of the grass-planting legume-cereal combinations. The accumulation of plant nutrients may benefit from the photosynthesis of plants, as the duration of sunlight in the mountain environment is relatively long. However, no existing research has demonstrated that applying nitrogen and phosphorus fertilisers in high-altitude regions promotes nutrient accumulation through photosynthesis in plants. Therefore, this study investigates the relationship between plant photosynthesis and fertiliser-derived nutrient accumulation by cultivating an annual grass-legume mixed cropping system in high-altitude areas. The authors propose the following points: 1: Fertilisers improve leaf physiological and biochemical indicators via photosynthesis; 2: Fertilisers enhance plant quality by accumulating nitrogen and phosphorus nutrients; 3: Identifying the most effective fertiliser types for high-altitude regions.
Materials and methods
Details about the test site in general
The field scientific observation and research station of the Ministry of Education (located at 33°24’30” north latitude, 97°18’00” east longitude, and at an altitude of 4270 m) is situated within the Sanjiangyuan grassland ecosystem in Chengduo County, Yushu Prefecture, Qinghai Province (Fig. 1). By collecting the climate data of this area and consulting relevant literature, it is determined that this region has a typical plateau continental climate. The winter is cold with low temperatures. The annual average maximum temperature is 24.1 degrees Celsius, the annual average minimum temperature is −26.6 degrees Celsius, and the annual average temperature is −1.4 degrees Celsius. The annual average precipitation is 406.00 millimeters, with 83.79% of the total precipitation concentrated from June to October. The number of cold weather days is approximately 260, and the average daily sunshine duration is about 2650.5 hours25. The annual average precipitation is 406.00 millimeters, with 83.79% of the total precipitation concentrated from June to October. The pH value of the soil is 6.92, the organic matter content is 2.36%, the total nitrogen content is 9.50 g per kilogram, the available nitrogen content is 14.0 milligrams per kilogram, the total phosphorus content is 8.20 g per kilogram, the available phosphorus content is 7.0 milligrams per kilogram, the total potassium content is 13.50 g per kilogram, and the available potassium content is 76.5 milligrams per kilogram.
Experimental sites in Chengduo County, Yushu Prefecture, Qinghai Province, China.(A) Avena sativa L planting trial test site.(B) Pisum sativum L planting trial test site.
Test materials
The Grassland Institute of the Qinghai Provincial Academy of Animal Husbandry and Veterinary Medicine supplied the oat and forage pea varieties for testing in this experiment, QingTian No1 and Qing Jian No1, respectively. The nitrogen fertilizer used in this experiment was urea (containing N 46%) and the phosphate fertilizer was calcium superphosphate (containing P2O5 12%).
Design of the experiment
This experiment was conducted from May 12th to September 20th, 2024. Based on the test results from the existing regional studies and the screening of the previous pre-experiment, the 70:30 mixed planting ratio of oats (Avena sativa L.) and forage peas (Pisum sativum L.) achieved the highest yield and quality26. Based on the above experimental results, this experiment added four fertilization treatments (Table 1): applying nitrogen fertilizer alone (A1), applying phosphorus fertilizer alone (A2), applying a nitrogen-phosphorus compound fertilizer (A3), and the no-fertilization control (CK) of this study, with the fertilizer amount referring to the results of Feng et al. A randomized complete block design was adopted, and each fertilization treatment had six replicates. No additional fertilization or irrigation was carried out during the entire growth period.
The experimental field measured 3 m×5 m (15 m² total area). One day prior to sowing, the field was plowed and leveled. Seeds were sown in parallel rows at a depth of 3–4 cm, with 10 manually dug rows per plot spaced 30 cm apart. Fertilizers were broadcast at the sowing sites. Harvesting was completed on September 25th. Both plant species were sown simultaneously, with the seeding density of grasses and legumes in each row determined by the monoculture seeding rates of individual planting treatments.
Methods and indicators
Description of leaf morphology characteristics
On September 20, 2024, scans were conducted on the leaves of herbaceous and leguminous forage plants in each experimental plot. A total of six measurements were taken for each plot. The measurement work was carried out using the Yaixin-1241 Leaf Area Measuring Instrument, which measured the leaf area, leaf perimeter, leaf length, and leaf width data of both grass and legume plants, and the results were recorded27,28.
Analysis of leaf spectral index characteristics
On September 22, 2024, using the CI-710 S plant spectrometer, the physiological and biochemical indicators of grass and legume plants in each treatment group were measured29. Each experimental field underwent 6 repetitions of experiments. The following parameters were measured: green index (G), normalized vegetation index (NDVI), chlorophyll A (CPHLA), chlorophyll B (CPHLB), total chlorophyll (CPHLT), soil-plant analysis development index (SPAD), chlorophyll B (CPHLB), plant senescence reflection index (PSRI), and water band index (WBI), and the data of grass and legume experiments were recorded separately30.
The formulae for the indicators are as follows
- (1)
: (NDVI = frac{{R800 – R680}}{{R800 + R680}}). The numerical range is from − 1 to 1, and r represents the reflectance value in the red light band.
- (2)
: (PSRI = frac{{R680 – R500}}{{R750}}) .The numerical range is from − 1 to 1, and r represents the reflectance value in the red light band.
- (3)
: (WBI = frac{{R900}}{{R970}}) .The numerical range is from 0.8 to 1.2, and Rr represents the reflectance value in the red light band.
- (4)
: (CPHLA = 12.72A663 – 2.59A645) .A663 and A645 are the absorbance at 663 nm and 645 nm.
- (5)
: (CPHLB = 22.88A645 – 4.67A663) ,A663 and A645 are the absorbance at 663 nm and 645 nm.
- (6)
:(CPHLT = CPHLA + CPHLB)
- (7)
: (G = frac{{R554}}{{R677}}), R is the reflection value in the red light band.
- (8)
: (SPAD = ( – 42.9 + (42.1*Math.pow(T931/T653,0.215))))
Feed quality inspection
On September 25, 2024, three 1-meter-long sample sections were randomly selected from each treatment area on both sides, and they were trimmed to be level with the ground. Subsequently, the grass samples of the mixed species of Gramineae and Fabaceae were placed in a drying oven and dried at 110 °C for 30 min, then dried to constant weight at 80 °C. The neutral detergent fiber (NDF) and acid detergent fiber (acidic detergent fiber, ADF) were analyzed using the Van Soest method; total phosphorus (TP) was determined using the sulfuric acid-peroxide digestion technique; crude protein (CP) was determined using the sulfuric acid-peroxide digestion method; soluble sugar (SS) was determined using the anhydroxyl method; total carbon (TC) and total nitrogen (TN) were determined using an elemental analyzer.
Statistical analysis
The data were summarized using Microsoft Excel 2016, and a one-way ANOVA was conducted using IBM SPSS 23.0. The data satisfied normality (Shapiro-Wilk, P > 0.05) and homogeneity of variance (Levene test, P = 0.12). The groups with significant differences were identified through Duncan’s multiple comparison test (P < 0.05). The relevant icons were drawn using the software named Origin 2021. The geographical map of the experimental area location was drawn using ArcMap 10.8.
Results
Changes in leaf morphology under different fertilizer treatments
Each of the three fertilization techniques improved the physiological morphological characteristics of the oat leaf blades. The results in Fig. 2 show that the fertilization treatments (A1, A2, and A3) improved the morphological characteristics of the oat tees more generally. A3 leaf area rose by 31.81% (P < 0.05) (Fig. 2a), A1 leaf perimeter increased by 25.61% (P < 0.05) (Fig. 2b), A1 leaf length increased by 17.71% (P < 0.05) (Fig. 2c), and A1’s leaf width expanded by 62.50% (P < 0.05) (Fig. 1d) in comparison to the no-fertilization (CK) treatment. Out of the three fertilizer treatments, the oat leaf morphology in the A1 treatment exhibited a superior trend, and there was no discernible difference between the A2 and A3 treatments overall.
The effects of different fertilization treatments on the leaf morphology of oat. The X-axis: The capital letters of the four treatment groups represent the separate application of nitrogen fertilizer (A1), separate application of phosphorus fertilizer (A2), nitrogen-phosphorus compound fertilizer (A3), and the control group without fertilizer application (CK). (a) The differences in leaf area caused by different fertilization treatments on the Y-axis. (b) The differences in leaf circumference under different fertilization treatments on the Y-axis. (c) The differences in leaf length under different fertilization treatments on the Y-axis. (d) The differences in leaf width under different fertilization treatments on the Y-axis. The “*” in figures a, b, c, and d indicate the values where there are significant differences among the treatments at the P < 0.05 level (Duncan test).;“ns” indicates no significant difference (P > 0.05).
The improvement of the forage pea’s leaf area, girth, and length was more strongly influenced by different fertilization treatments (Fig. 3). The forage pea leaf area increased by 27.08%, 28.50%, and 34.69% (P < 0.05) in A1, A2, and A3 when compared to the no-fertilizer (CK) treatment; the forage pea leaf perimeter increased by 76.43% (P < 0.05) in A1 treatment; and the forage pea leaf length increased by 9.90%, 12.58%, and 13.40% (P < 0.05) in A1, A2, and A3 treatments. However, there was no discernible difference in the forage pea leaf width. Additionally, there were no notable variations among the three fertilizer applications.
The effects of different fertilization methods on the leaf morphology of forage peas. The X-axis: The capital letters in the 4 treatment groups represent the separate application of nitrogen fertilizer (A1), separate application of phosphorus fertilizer (A2), nitrogen-phosphorus compound fertilizer (A3), and the control group without fertilizer application (CK). (a) The Y-axis represents the differences in leaf area under different fertilization treatments. (b) The Y-axis represents the differences in leaf perimeter under different fertilization treatments. (c) The Y-axis represents the differences in leaf length under different fertilization treatments. (d) The Y-axis represents the differences in leaf width under different fertilization treatments. The “*” in figures a, b, c, and d indicate the values where there are significant differences among the treatments at the P < 0.05 level (Duncan test).;“ns” indicates no significant difference (P > 0.05).
Variations in the spectral index properties of leaves treated with various fertilizers
The different fertilizer treatments had a substantial impact on the soil-vegetation analysis development, but not on the oat spectral index. The results evaluated in Fig. 4 showed that the total chlorophyll amount in A3 was 17.07% (P < 0.05) higher than in A1 treatment and 15.14% (P < 0.05) higher than in CK treatment. In comparison to the A1 and A2 treatments, the A3 treatment improved the soil-vegetation analysis development indexes by 17.05% and 26.85% (P < 0.05), respectively, and by 26.44% (P < 0.05) when compared to the CK treatment. The water band index index increased by 4.81% (P < 0.05) in A2 compared to A3.Together, these findings showed that fertilizer treatment did not significantly affect the overall oat spectral index; however, the soil-vegetation analysis development indices varied significantly.
The influence of fertilization on the spectral index characteristics of oat leaves. X-axis: The capital letters in the four treatment groups respectively represent the separate application of nitrogen fertilizer (A1), the separate application of phosphorus fertilizer (A2), the application of nitrogen-phosphorus compound fertilizer (A3), and the control group without fertilizer application (CK) (a) Y-axis represents the change in chlorophyll A content. (b) Y-axis represents the change in chlorophyll B content. (c) Y-axis represents the change in total chlorophyll content. (d) Y-axis represents the change in green index content. (e) Y-axis represents the change in NDVI (Normalized Difference Vegetation Index) content. (f) Y-axis represents the change in plant senescence reflection index content. (g) Y-axis represents the change in soil-vegetation analysis development index content. (h) Y-axis represents the change in water band index content. In figures a, b, c, d, e, f, g, h, “*”, indicates significant differences among the treatments at the P < 0.05 level (Duncan’s test); “ns” indicates no significant difference (P > 0.05).
The results in Fig. 5 showed that there was no significant difference in the spectral index features of fodder pea leaves across the different fertilization treatments. While the other indexes were not substantially differ from one another, the CK treatment showed a superior trend in the total chlorophyll index, which was 17.97% and 20.25% greater than that of the A1 and A3 treatments (P < 0.05). A detailed analysis of the previously given data showed that fertilizer treatments had no appreciable effect on the overall spectral index of forage peas, but the CK treatment demonstrated a superior trend in total chlorophyll.
The influence of different fertilization treatments on the spectral index characteristics of forage pea leaves. The X-axis: The capital letters in the four treatment groups represent the separate application of nitrogen fertilizer (A1), separate application of phosphorus fertilizer (A2), application of nitrogen-phosphorus compound fertilizer (A3), and no fertilizer application as the control group (CK). (a) The Y-axis represents the change in chlorophyll A content. (b) The Y-axis represents the change in chlorophyll B content. (c) The Y-axis represents the change in total chlorophyll content. (d) The Y-axis represents the change in green index content. (e) The Y-axis represents the change in NDVI (Normalized Difference Vegetation Index) content. (f) The Y-axis represents the change in plant senescence reflection index content. (g) The Y-axis represents the change in soil-vegetation analysis development index content. (h) The Y-axis represents the change in water band index content. In figures a, b, c, d, e, f, g, and h, “*” indicates significant differences among the treatments at the P < 0.05 level (Duncan’s test); “ns” indicates no significant difference (P > 0.05).
Variations in the quality of the feed with various fertilizer treatments
Fertilization procedures significantly improved the fodder’s quality. The results of the analysis showed that the three fertilization treatments significantly reduced the contents of acid detergent fiber and neutral detergent fiber compared to the CK treatments by 8.67% and 19.28% (P < 0.05) in the A1 treatment, 11.25% and 28.92% (P < 0.05) in the A2 treatment, and 7.78% and 25.25% (P < 0.05) in the A2 treatment (Fig. 6). Forage soluble sugar, crude protein, total phosphorus, and total nitrogen all showed a discernible improvement as compared to CK therapy. When the three fertilization treatments were compared, the A3 treatment had lower amounts of neutral and acid detergent fiber and higher levels of soluble sugar, crude protein, total phosphorus, total nitrogen, and total carbon.
Changes in plant quality under different fertilizer treatments. X-axis: The capital letters in the four treatment groups represent the separate application of nitrogen fertilizer (A1), separate application of phosphorus fertilizer (A2), application of nitrogen-phosphorus compound fertilizer (A3), and no fertilizer application as the control group (CK). (a) The Y-axis represents the change in neutral detergent fiber content. (b) The Y-axis represents the change in acidic detergent fiber content. (c) The Y-axis represents the change in crude protein content. (d) The Y-axis represents the change in soluble sugar content. (e) The Y-axis represents the change in total nitrogen content. (f) The Y-axis represents the change in total phosphorus content. (g) The Y-axis represents the change in total carbon content. * indicates significant differences, ns indicates no significant differences. In figures a, b, c, d, e, f, and g, “*” indicates significant differences among the treatments at the P<0.05 level (Duncan’s test);“ns” indicates no significant differences (P>0.05).
Changes in plant quality under different fertilizer treatments. X-axis: The capital letters in the four treatment groups represent the separate application of nitrogen fertilizer (A1), separate application of phosphorus fertilizer (A2), application of nitrogen-phosphorus compound fertilizer (A3), and no fertilizer application as the control group (CK). (a) The Y-axis represents the change in neutral detergent fiber content. (b) The Y-axis represents the change in acidic detergent fiber content. (c) The Y-axis represents the change in crude protein content. (d) The Y-axis represents the change in soluble sugar content. (e) The Y-axis represents the change in total nitrogen content. (f) The Y-axis represents the change in total phosphorus content. (g) The Y-axis represents the change in total carbon content. * indicates significant differences, ns indicates no significant differences. In figures a, b, c, d, e, f, and g, “*” indicates significant differences among the treatments at the P < 0.05 level (Duncan’s test); “ns” indicates no significant differences (P > 0.05).
Analysis of correlations between fodder quality, morphological traits, and leaf spectral index characteristics
Figure 6’s results show that while fodder quality and oat forage leaf spectral features and overall leaf shape are positively connected, NDF and ADF indexes are substantially inversely connected. Furthermore, there is a strong negative correlation between NDF and ADF indicators and leaf shape, although there is a significant positive correlation with forage soluble sugars, crude proteins, total phosphorus, total nitrogen, and total carbon (Fig. 7). Forage soluble sugars, crude protein, total phosphorus, total nitrogen, and total carbon all showed a strong negative connection with the NDF and ADF indexes. Furthermore, the morphology of the unique blades as a whole and the spectral properties of forage pea forage blades were positively correlated, as was the forage quality.
Correlation analysis of spectral index characteristics, morphological characteristics of leaf area and forage quality. These indices are leaf area, perimeter, length and width. The green index (G), normalized vegetation index (NDVI), chlorophyll A (CPHLA), chlorophyll B (CPHLB), total chlorophyll (CPHLT), soil-plant analysis development index (SPAD), plant senescence reflection index (PSRI), water band index (WBI), neutral detergent fiber (NDF), acid detergent fiber (ADF), crude protein (CP), soluble sugar (SS), total nitrogen (TN), total carbon (TC), total phosphorus (TP). Red indicates positive correlation, blue indicates negative correlation, the color shade represents the intensity (-1 to +1). The flatter the ellipse in the figure, the stronger the correlation. The direction of inclination indicates positive correlation (right upper – left lower inclination) / negative correlation (left upper – right lower inclination).
Construction of structural equation model for forage quality, morphological characteristics and leaf spectral index features
Based on the spectral index characteristics of the leaves, the morphological characteristics of the leaves, and the quality indicators of the forage grass, a random forest model was constructed (Fig.8A) to identify the main driving factors affecting the change in crude protein. The results showed that variables such as GLA, LLC, TP, SPAD, CPHLT, fertilizer application amount (F), NDVI, PSRI, WBI, G, etc., in the grass family and leguminous plants were the main variables affecting the crude protein index. According to the analysis results of the structural equation model (SEM)(Fig.8B), the overall fit of the model was good: Fisher’s C=12.693, degrees of freedom df = 8, P= 0.123>0.05, indicating that there was no significant difference between the model and the data, and it was acceptable. The explanatory power (R²) of each endogenous variable showed that variables such as WBI-G (R² = 0.70), WBI-L (R² = 0.62), F (R² = 0.61), CP (R² = 0.55), and GLA (R² = 0.42) had a relatively high proportion of variance explanation, while the explanatory power of CPHLB-L (R²= 0.19), NDVI-L (R²=0.23*), CPHLT-G (R²=0.25*), and LLC (R²=0.26) was relatively low but partially reached a significant level (p<0.05,***p<0.001). F, WBI-G, CPHLT-L, TP, WBI-L had direct effects on CP and were significantly positively correlated (P<0.05), CPHLB-L indirectly affected CP through F, and LLC indirectly affected CP through WBI-L.
Structural equation model of spectral index characteristics, leaf morphology characteristics and forage quality. (A) is the result analysis diagram of the random forest model, (B) is the structural equation model constructed based on the results of the random forest. These indicators include leaf area, perimeter, length and width of herbaceous and leguminous plant leaves. Green index (G), normalized vegetation index (NDVI), chlorophyll A (CPHLA), chlorophyll B (CPHLB), total chlorophyll (CPHLT), soil-plant analysis development index (SPAD), plant senescence reflection index (PSRI), water band index (WBI), neutral detergent fiber (NDF), acid detergent fiber (ADF), crude protein (CP), soluble sugar (SS), total nitrogen (TN), total carbon (TC), total phosphorus (TP). Red indicates negative correlation, green and blue indicate positive correlation, and the depth of color represents the strength of correlation.
Discussion
Effect of fertilizer treatments on morphological characteristics of leaves
Leaves are one of the important organs of forage. The overall morphology of leaves has a large surface area, which is conducive to the input of water and nutrients as well as the output of photosynthetic products31. The results of the leaf morphology indicators selected in this study indicate that fertilization treatment has a significant impact on the morphology of oat leaves. Among them, the nitrogen and phosphorus mixed application treatment results in larger leaf area for oat, while the single application of nitrogen fertilizer treatment has higher leaf circumference, length, and width, and there are significant differences between this treatment and the non-fertilization control treatment. These differences are mainly attributed to the specificity of the physiological functions of different elements32. Nitrogen is the core component of chlorophyll and photosynthetic proteins, mainly driving cell division and elongation, so single application of nitrogen fertilizer significantly increases the circumference, length, and width of the leaves; while phosphorus is involved in the construction of energy carriers (ATP) and the transportation of photosynthetic products, and is more focused on influencing the differentiation of mesophyll cells and the expansion of single leaf area. When nitrogen and phosphorus are mixed applied, they form a synergistic effect, providing a protein framework and ensuring energy supply, enabling oat to overcome nutrient limitation and form the largest leaf area33. The differences between species can also be analyzed. Oats, as a grass plant, does not have nitrogen fixation ability and is highly sensitive to soil nitrogen. Exogenous nitrogen application directly stimulates the comprehensive increase of their leaf morphology; while for forage peas, as a leguminous plant, it can fix nitrogen through symbiotic rhizobia, and exogenous nitrogen fertilizer may trigger the “nitrogen repression” effect, causing it to prioritize resources allocation to oat plants for stem elongation to compete for light34. This explains why fertilization has no significant effect on the width of pea leaves35. In addition, fertilization promotes the rapid growth and height increase of oats, causing forage peas to be in a lower spatial position in the community and being shaded36,37. The shading stress induces peas to increase leaf area and length to enhance light capture efficiency, which is consistent with the leaf response trend after adding phosphorus. This further amplifies the phenotypic differences between grass and leguminous plants among different fertilization treatments38. Based on these findings, we conclude that the combination of nitrogen and phosphorus has the most significant impact on leaf morphology and is the best fertilization strategy for high-altitude regions.
Effect offertilization treatments on spectral index characteristics of leaves
The strength of photosynthetic capacity of leaves is a key parameter determining crop yield, and premature leaf senescence is an important factor limiting high yield. Its main characteristic is the functional loss of the photosynthetic apparatus in the leaf, leading to a decline in photosynthetic function. Reasonable fertilization can reduce the nutrient output of leaves after flowering, delay leaf senescence, and maintain a higher photosynthetic rate39. The research results show that in oat plants under different fertilization treatments, the CPHLT and SPAD indices of the nitrogen-phosphorus combined application treatment are at a relatively high level, which is significantly improved compared to the single application of nitrogen fertilizer and no fertilization treatment. The WBI index of leaves under the single application of phosphorus fertilizer is relatively high, significantly superior to the nitrogen-phosphorus combined application treatment. Other leaf spectral characteristic indicators do not have significant advantages. The fertilization treatments have no significant impact on the spectral indices of forage peas except for CPHLT. This is mainly because oat is a C₃ grass plant, and its photosynthetic apparatus is highly sensitive to nutrient supply40. The addition of exogenous nitrogen and phosphorus fertilizers can directly promote chlorophyll synthesis and delay the degradation of chloroplasts, so the total chlorophyll content of oat under the nitrogen-phosphorus combined application treatment significantly increases, the light energy utilization efficiency is enhanced, the leaf senescence process is effectively delayed, and this reflects the positive regulatory effect of fertilization on grass crops41. For forage peas as a leguminous plant, their roots can symbiotically fix nitrogen with rhizobia and have a lower dependence on soil exogenous nitrogen. When exogenous fertilizers are applied, peas may preferentially utilize the nitrogen fixed by themselves rather than the increased mineral nitrogen in the soil, resulting in their sensitivity to fertilization responses being much lower than that of oat42. This explains why the spectral characteristic indicators of pea leaves do not show significant improvement under fertilization treatments. Moreover, in the high-altitude environment background, oat plant height is already limited by natural conditions such as low temperature and strong wind. In the oat and pea mixed-sowing system, the nitrogen-phosphorus combined application promotes oat plants to obtain more sufficient nutrient supply, resulting in a significant increase in plant height and canopy density43. Although forage peas are climbing plants, their upward growth height is limited by the height of oat plants. In the high-altitude conditions, the absolute value of oat plant height is relatively low, and the climbing space of peas is physically compressed, causing them to be forced to be in the lower part of the canopy44. Therefore, under fertilization treatment, oat leaves are more lush, the canopy density is significantly increased, and they have a stronger shading effect on the lower peas, resulting in a significant reduction in the photosynthetic effective radiation received by pea leaves and a lower light energy acquisition ability than the no-fertilization treatment. This shading effect not only reduces the photosynthetic rate of pea leaves but may also induce a compensatory increase in the total chlorophyll content of peas, resulting in a phenomenon where the total chlorophyll content of peas under fertilization treatment is higher than that under no-fertilization treatment. This is essentially an adaptation mechanism in a low-light environment, rather than a superiority in nutrient supply45. At the same time, the effects of fertilization treatment on the distribution of canopy water and metabolic activity also show species-specific differences. Under the single application of phosphorus fertilizer treatment, the water band index of leaves is higher, indicating that phosphorus fertilizer promotes the increase of water content in the plant canopy and the enhancement of metabolic rate. However, in the oat and pea mixed-sowing system, it may be affected by interspecific competition. Oat, with stronger nutrient absorption ability and plant height advantage, occupies a dominant position in water and nutrient acquisition, while peas, due to their inferior position in the canopy, have their leaf water content, metabolic activity, and photosynthetic efficiency inhibited. In addition, this study only assesses the accumulation and conversion efficiency of nutrients based on leaf spectral index characteristics46. Although spectral parameters such as total chlorophyll content, water band index, and soil-vegetation analysis development index can well reflect the canopy structure, chlorophyll content, and water status, they cannot deeply determine the content of various substances in leaves and the activity of key enzymes47. Therefore, spectral characteristics are difficult to fully reveal the internal mechanisms of the interaction between oat and pea under fertilization treatment in terms of photosynthetic products and canopy competition dynamics. Future research needs to combine physiological and biochemical measurements at the leaf scale to systematically clarify the underlying mechanisms of the species-specific differences in the fertilization effect in the mixed-sowing system.
Effect offertilization treatments on forage quality
The main purpose of adding exogenous fertilizers is to achieve higher target yields, corresponding forage quality, and to enhance soil fertility48. The results of this study indicate that there are significant differences in the regulatory effects of different fertilization treatments on forage fiber components and nutritional quality. All three fertilization treatments significantly reduced the neutral detergent fiber and acid detergent fiber contents of forage, with the single application of phosphate fertilizer having the most prominent effect, reducing by 11.25% and 28.92% respectively. This is mainly because phosphorus directly participates in the regulation of plant cell wall metabolism, and phosphorus is an important component of ATP and nucleic acids, which can promote energy metabolism and protein synthesis in plants, accelerate cell division and elongation49. Plants tend to allocate more photosynthetic products to protoplasm synthesis, resulting in a decrease in the proportion of fiber components in dry matter. Moreover, the application of phosphate fertilizer can promote root development and nutrient absorption efficiency, enhance the overall metabolic activity of the plant, accelerate the transport and transformation of carbohydrates, and increase the accumulation of soluble sugars and other non-structural carbohydrates, thereby diluting the relative content of fiber components50. At the same time, the fertilization treatments significantly increased the contents of soluble sugars, crude protein, total phosphorus, and total nitrogen, which reflects the comprehensive promoting effect of exogenous nutrient addition on the productivity of grassland ecosystems. Among them, the nitrogen-phosphorus combined application treatment had a more significant effect compared to the no-fertilization treatment, because grass plants lack biological nitrogen fixation ability, and their growth and development are highly dependent on soil mineral nitrogen. The addition of exogenous nitrogen directly alleviates nitrogen limitation and promotes the tillering, leaf expansion, and photosynthetic carbon assimilation of grass plants51. Moreover, the addition of nitrogen activates the nitrogen fixation activity of leguminous plants through the “nitrogen stimulation effect”, and the leguminous plant root nodule nitrogen fixation system gradually plays a supplementary role in the early growth stage, forming a dual supply pattern of exogenous nitrogen and biological nitrogen fixation52. Based on the above analysis, nitrogen-phosphorus combined application can simultaneously meet the large demand for nitrogen by grass plants and the sensitive dependence of leguminous plants on phosphorus, avoiding the possible nutrient imbalance caused by single fertilization. Moreover, the mixed planting system itself has the nitrogen-fixing function of leguminous plants and the exogenous nitrogen-phosphorus supply forms a synergistic effect. The addition of nitrogen not only directly promotes the growth of grass plants but also indirectly improves the utilization efficiency of phosphorus through optimizing the canopy structure and extending the photosynthetic functional period, forming a positive feedback loop53. The formation of the differences in fertilization effects in this study is essentially the result of the spatial heterogeneity distribution of phosphorus, the differentiation of nutrient acquisition strategies of grass and leguminous plants, and the dynamic competition and complementarity of resources among species in the mixed planting system. The enrichment characteristic of phosphorus in the surface soil determines its direct and efficient regulation of fiber components, while nitrogen-phosphorus combined application achieves multi-target synergy in improving forage yield, quality, and soil fertility through coupling the rapid response of grass plants to mineral nitrogen and the ecological service function of leguminous plants’ nitrogen fixation. This forms a positive feedback loop54. This study shows that the mixed application of nitrogen and phosphorus can simultaneously meet the large demand for nitrogen by grass plants and the sensitive dependence of leguminous plants on phosphorus, avoiding the possible nutrient imbalance caused by single fertilization. Moreover, the mixed planting system itself has the nitrogen-fixing function of leguminous plants and the exogenous nitrogen-phosphorus supply forms a synergistic effect. The addition of nitrogen not only directly promotes the growth of grass plants but also indirectly improves the utilization efficiency of phosphorus through optimizing the canopy structure and extending the photosynthetic functional period, forming a positive feedback loop. This study provides theoretical basis and practical feasibility for its promotion and establishment in high-altitude regions.
Conclusions
The results of the different fertilization treatments conducted at high altitudes on the feed quality produced by the mixed planting of leguminous plants and forage grasses show that fertilization can improve the physiological and biochemical indicators of leaves through photosynthesis and enhance the feed quality. Compared with the untreated treatment, the effect is better. In all treatment groups, the group that applied a mixture of nitrogen and phosphorus showed better performance in leaf morphology, photosynthetic capacity, and nutritional components than the untreated group, and was able to improve plant quality by accumulating nitrogen and phosphorus nutrients. Therefore, for high-altitude regions, the most effective fertilizer type is the treatment of combined application of nitrogen and phosphorus. However, it should be noted that this conclusion only applies to the mixed planting combination of oats and forage pea and their respective ratios. Additionally, this research trial was a one-year short-term experiment and lacked long-term longitudinal studies. Long-term related research will be conducted in the future.
Data availability
Data sets are available from corresponding authors on reasonable request.
References
Gou, W. L. et al. Effects of mixing ratio and mowing stubble on root morphology and soil nutrients in annual grass-bean mixed-grassland. Grassland Turf. 40, 1–7 (2020).
Google Scholar
Xiang, X. M. et al. Response of vegetation and soil to exogenous nitrogen inputs in alpine meadow grassland in the Sanjiangyuan area. J. Grasslands. 29, 2010–2016 (2021).
Google Scholar
Coban, O. et al. Soil microbiota as game-changers in restoration of degraded lands. Science 375, 0725 (2022).
Google Scholar
Li, W. L. et al. Impact of restoration measures on plant and soil characteristics in the degraded alpine grasslands of the Qinghai Tibetan Plateau: A meta-analysis. Agric. Ecosyst. Environ. 2, 108394 (2023).
Google Scholar
John, P. et al. Organic agriculture in the twenty-first century. Nat. Plants 2, 15221 (2016).
Google Scholar
Zhao, N. et al. Cropping practices manipulate soil bacterial structure and functions on the Qinghai–Tibet Plateau. J. Plant Physiol. 271, 153666 (2022).
Google Scholar
Xie, Y. et al. fertilizer effects of nitrogen, phosphorus and potassium on alfalfa in the dam area and recommended fertilizer application rates. Chin. J. Grassland. 34, 52–57 (2012).
Ian, J. et al. The worldwide leaf economics spectrum. Nature 428, 821–827 (2004).
Google Scholar
Mahmoud, G. et al. Pre-cecal phosphorus digestibility for corn, wheat, soybean meal,and corn gluten meal in growing Japanese quails from 28 to 32 d of age. Anim. Nutr. 5, 148–151 (2019).
Google Scholar
Terry, J. et al. Phytotoxicity of soilborne glyphosate residues is influenced by the method of phosphorus fertilizer application. Plant Soil 442, 445–465 (2017).
Bicharanloo, B. et al. Nitrogen and phosphorus availability have stronger effects on gross and net nitrogen mineralisation than wheat rhizodeposition. Geoderma 405, 26–35 (2022).
Google Scholar
Evans, J. R. Photosynthesis and nitrogen relationships in leaves of C3 plants. Oecologia 78, 9–19 (1989).
Google Scholar
X, R. et al. Grass-legume mixture and nitrogen application improve yield, quality, and water and nitrogen utilization efficiency of grazed pastures in the loess plateau. Front. Plant Sci. https://doi.org/10.3389/fpls.2023.1088849 (2023).
Google Scholar
Arif, I. et al. Plant microbiome engineering: Expected benefits for improved crop growth and resilience. Trends Biotechnol. 38, 1385–1396 (2020).
Google Scholar
Wang, C. G. et al. Effects of nitrogen, phosphorus and potassium fertilizer dosage on rice yield and quality. Jilin Agricultural Sci. 35, 28–33 (2010).
Bedison, J. E. et al. Is the growth of temperate forest trees enhanced along an ambient nitrogen deposition gradient?. Ecology 90, 1736–1742 (2009).
Google Scholar
Walters, M. B. et al. Intraspecific growth and functional leaf trait responses to natural soil resource gradients for conifer species with contrasting leaf habit. Tree Physiol. 33, 297–310 (2013).
Google Scholar
Zhang, Y. L. et al. Effects of fertilization and mixed sowing ratio on the production perfor-mance of grass-alfalfa mixed forage. Chin. J. Grassland. 42 (6), 115–124 (2020).
Gou, W. L. et al. Progress of researches on yield increase and efficiency of grass-bean forage mixture. Herbology 3, 16–23 (2020).
Fan, Y. et al. Optimizing cropping systems of cultivated pastures in the mountain–basin systems in Northwest China. Appl. Sci. (Basel) 10, 6949 (2020).
Google Scholar
Dong, K. W. et al. More fertilizer and impoverished roots required for improving wheat yields and profits under climate change. Field Crops Res. 249, 107–116 (2020).
Erin, M. et al. Structural equation modeling reveals complex relationships in mixed forage swards. Crop Prot. 78, 106–113 (2015).
Google Scholar
Letícia, D. F. et al. Application of Azospirillum on seeds and leaves, associated with Rhizobium inoculation, increases growth and yield of common bean. Arch. Microbiol. 203, 1033–1038 (2021).
Google Scholar
Mortenson, M. C. et al. Carbon sequestration in rangelands interseeded with yellow-flowering alfalfa (Medicago sativa ssp. falcata). Environ. Manage. 33, S475–S481 (2004).
Google Scholar
Feng, T. X. et al. Response of Leaf Morphology and Yield of Annual Grass-Bean Mixed Forage to Fertiliser Application in the Sanjiangyuan Area. PLoS ONE. 32, 3951–3961 (2025).
Feng, T. X. et al. Nitrogen and phosphorus fertilisers optimise root morphology and soil nutrients in mixed annual grass and bean sown grassland in alpine regions. J. Grassland Sci. 3, 15–20 (2025).
Korva, J. T. et al. A simple and low-cost method for leaf area measurement of detached leaves. Exp. Agric. 07, 10–15 (2002).
Yang, Z. S. Low-cost 3D phenotype measurement method of leafy vegetables using video recordings from smartphones. Sensors 10, 15–19 (2020).
Google Scholar
Kupčinskienė, A. et al. Vegetation indices for early grey mould detection in lettuce grown under different lighting conditions. Plants 11, 18–26 (2023).
Zhang, X. et al. Hyperspectral estimation of canopy chlorophyll of winter wheat by using the optimized vegetation indices. Comput. Electron. Agric. 1, 1–12 (2022).
Osnas, J. L. D. et al. Global leaf trait relationships: mass, area, and the leaf economic spectrum. Science 340, 741–744 (2014).
Google Scholar
Rivas-Ubach, A. et al. Strong relationship between elemental stoichiometry and metabolome in plants. Proc. Natl. Acad. Sci. U.S.A. 109, 4181–4186 (2012).
Google Scholar
Daniel, E. et al. Effect of soil phosphorus availability and residue quality on phosphorus transferfrom crop residues to the following wheat. Plant. Soil. 416, 1–2 (2017).
Google Scholar
Khuyen, T. K. H. P pools after seven-year P fertilizer application are influeced by wheat straw addition and wheat growth. J. Soil Sci. Plant Nutr. 19, 603–610 (2019).
Google Scholar
Aleric, K. M. et al. Growth and photosynthetic re-sponses of the federally endangered shrub, Lindera melissi-folia(Lauraceae), to varied light environments. Am. J. Bot. 92, 682–689 (2005).
Google Scholar
Xie, H. T. et al. Leaf non-structural carbohy-drate allocation and C:N:P stoichiometry in response to light acclimation in seedlings of two subtropical shade-toler-ant tree species. Plant Physiol. Biochem. 124, 146–154 (2018).
Google Scholar
Myers, J. A. et al. Carbohydyate storage enhances seedling shade and stress tolerance in a neotropical forest. Joural Ecol. 95, 383–395 (2007).
Google Scholar
Rezvan, K. et al. Soil carbon dynamics in wheat plots establishedon grassland in 1911 as influenced by nitrogen and phosphorus fertilizers. Can. J. Soil Sci. 8, 158–160 (2018).
Shi, H. T. et al. Modulation of auxin content in Arabidopsis confers improved drought stress resistance. Plant Physiol. Biochem. 82, 209–217 (2014).
Google Scholar
Zhang, K. et al. Effects of low light at seedling stage on photosynthetic characteristics of peanut. Chin. Agricultural Sci. 43, 65–71 (2010).
Google Scholar
Ye, Z. P. Progress in modelling photosynthesis response to light and CO2. J. Plant. Ecol. 34, 727–740 (2010).
Google Scholar
Gabriele, B. et al. Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 1, 1–13 (2010).
Deyn, G. B. D. et al. Plant species identity and diversity effects on different trophic levels of nematodes in the soil food web. Oikos 106, 576–586 (2004).
Google Scholar
Ye, Y. et al. Responses of soil nematode community to monoculture or mixed culture of a grass and a legume forage species in China. Pedosphere 30, 791–800 (2020).
Google Scholar
Nyfeler, D. et al. Grass-legume mixtures can yield more nitrogen than legume pure stands due to mutual stimulation of nitrogen uptake from symbiotic and nonsymbiotic sources. Agric. Ecosyst. Environ. 140, 155–163 (2011).
Google Scholar
Lin, W. S. et al. Meta-analysis of the effects of nitrogen and phosphorus additions on carbon, nitrogen and phosphorus stoichiometric characteristics of soils in alpine meadows of Qinghai-Tibet. J. Grassland Sci. 30, 3345–3354 (2022).
Google Scholar
Renzen, W. M. et al. Analysis of the effect of nitrogen and phosphorus addition on the production and ecological function of Pennisetum lanceolatum artificial grassland. Plateau Sci. Res. 4, 55–61 (2020).
Halim, R. A. et al. Productivity, relative yield and plant growth of forage corn intercropped with soybean under different crop combination ratio. Legume Res. 39, 558–564 (2016).
Wang, L. J. et al. Effects of ALA on photosynthesis and chlorophyll fluorescence characteristics of radish leaves at different leaf positions. Northwest. J. Bot. 3, 12 (2005).
Zhang, X. Z. et al. Effects of nitrogen, phosphorus and potassiumfertilization rates on yield, quality and economic benefits of mixed-seeded grassland. Grassland Turf. 34, 56–60 (2014).
Song, J. C. et al. Effects of nitrogen and phosphorusfertilization on forage yield and nutritive quality of Pisum sativum in alpine zone. J. Grassland Sci. 30, 732–742 (2022).
LI, J. K. et al. Effects of different mixing ratios of oats and arrow end peas on forage yield and quality in Tibet. J. Grassland Sci. 19, 830–833 (2011).
Google Scholar
Luiz, G. D. O. et al. Mixed grass–legume pastures in integrated crop-livestock systems: A strategy to improve soil health and soybean yield in the Brazilian Cerrado. Plant Soil 2, 1–12 (2025).
Feng, T. X. et al. Effects of different mix combinations and proportions of oats and peas on forage yield and quality in alpine regions. J. Grassland Sci. 30, 487–494 (2022).
Google Scholar
Acknowledgements
Besides the financial support, the authors, would also like to acknowledge the support of all the stakeholders and staff involved in the research task, without whom this project would not have been possible!
Funding
Comprehensive Demonstration of Feeding Grass Supply and Winter Livestock Utilization Quality Improvement and Efficiency Enhancement in Chengduo Area (2024-NK-137); Research and Demonstration of Nutritional Regulation Technology of Artificial Grasslands in Chengduo Area (2025-NK-P27); Funded by the Key Laboratory of Three River Source Ecosystem of Qinghai University.
Author information
Authors and Affiliations
Contributions
T.X.F(Ting-xu Feng) Conducting experiments, collating and analysing data, producing tables and images, writing original drafts of papers, writing- reviewing and editing. F.L(Fei Li) Help with data collection and analysis.W.S.L(Wei-shan Lin) Help with experiments and data collection. L.Z(Lin Zhang) Conduct article supervision and essay writing checks. X.M.X(Xue-mei Xiang) Provided ideas for article writing and helped with experiments. X.J.W(Xi-jie Wei) Perform experimental material acquisition and test setup. K.J.D(Ke-jia De) Performed project management, provided funding, supervised the review of draft documents and approved final drafts.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Consent for publication
All the authors involved in this article have agreed that “The quality of artificial grasslands in alpine regions is controlled by crop photosynthesis under different fertilization conditions” will be published in the Scientific Reports.
Commitment
The author promises that all data appearing in the text have been obtained by actual testing. That all claims expressed in the text are those of the authors only and do not necessarily represent the claims of their affiliated organisations, or those of the publisher, editors and reviewers. Any product that may be evaluated in this text, or claims that may be made by its manufacturer, are not warranted or endorsed by the publisher.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Reprints and permissions
About this article
Cite this article
Feng, Tx., Li, F., Lin, Ws. et al. The quality of artificial grasslands in alpine regions is controlled by crop photosynthesis under different fertilization conditions.
Sci Rep 16, 19540 (2026). https://doi.org/10.1038/s41598-026-53159-2
Received:
Accepted:
Published:
Version of record:
DOI: https://doi.org/10.1038/s41598-026-53159-2
Keywords
- Alpine region
- Fertilizer containing nitrogen and phosphate
- Grass-bean mix
- Characteristics of the leaf spectral index
- Leaf shape
- Forage quality
Source: Ecology - nature.com
