Next Article in Journal
Effective Pollen-Fertility Restoration Is the Basis of Hybrid Rye Production and Ergot Mitigation
Previous Article in Journal
An Overview of Cryo-Scanning Electron Microscopy Techniques for Plant Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forage Yield, Canopy Characteristics, and Radiation Interception of Ten Alfalfa Varieties in an Arid Environment

1
State Key Laboratory of Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China
2
College of Grassland Agriculture, Northwest A&F University, Yangling 712100, China
*
Author to whom correspondence should be addressed.
Plants 2022, 11(9), 1112; https://doi.org/10.3390/plants11091112
Submission received: 3 March 2022 / Revised: 15 April 2022 / Accepted: 16 April 2022 / Published: 20 April 2022
(This article belongs to the Section Crop Physiology and Crop Production)

Abstract

:
An increasing demand for new and improved livestock forage products is driving the development of forage systems in arid regions. Our study evaluated the productivity of 10 alfalfa (Medicago sativa L.) varieties and the relationship between forage yield and canopy structure traits, such as plant height, branch number, stem/leaf ratio, and leaf area index in the arid Hexi Corridor, north-west China. Here, plant height, primary branch number per plant, and stem/leaf ratio were positively correlated with forage yield. In terms of a two-year total yield, Gannong No. 5 produced the highest value (13,923 kg ha−1), followed by the WL342HQ (12,409 kg ha−1), Phabulous (11,928 kg ha−1), and Xinjiang Daye (11,416 kg ha−1) varieties. Therefore, these four alfalfa varieties are suitable for large-scale cultivation in the Hexi Corridor region and other arid areas where the effect of precipitation is even larger than that of temperature. These results provide valuable information for the selection and cultivation of alfalfa varieties, which could improve forage yield and the production of livestock in arid regions.

1. Introduction

Arid and semi-arid environments cover about 30% of the world’s total land area and 52% of that in China [1]. Additionally, the area of semi-arid land increased by 25% from 1986 to 2010 [2]. Livestock production is a major economic contributor in arid and semi-arid regions [3]. However, water is often considered to be the most important determinant that restricts forage yield in these environments [4]. Forage crops usually suffer from drought stress in these regions, even under irrigation, owing to high annual evaporation. Therefore, it is imperative to screen suitable forage crops and varieties which are more adaptable to water-limited conditions to improve livestock production sustainability in arid and semi-arid regions.
Alfalfa (Medicago sativa L.) is a widely planted, perennial forage legume in arid and semi-arid areas because of its high productivity, nutritional value, nitrogen fixation capacity, and adaptability to various environments [5,6]. The canopy structure has a significant influence on the forage yield of alfalfa owing to the optimization of light energy utilization and accumulation. Some traits that influence canopy development and productivity are plant height (PH), the stem to leaf ratio, branch number, leaf area index (LAI), and light interception (LI) [7,8]. Previously, a study estimated the PH of different alfalfa crop canopies, and established that PH is positively correlated with forage yield [9]. The branch number is a vital yield component of forage crops [10,11]. It exerts a strong positive impact on forage yield [12,13], and has been included in plant growth models to improve predictive ability [14]. Generally, dry matter accumulation was found to be reduced in proportion to the level of the water deficit to which the plants were exposed, because such conditions reduced shoot density and leaf area [15]. Another study shows that the number of shoots per plant varies between genotypes differing in fall dormancy [16]. Elsewhere, research has indicated that there is a positive genetic correlation between the stem/leaf ratio and forage yield [17]. A decrease in the stem/leaf ratio was found to be accompanied by reduced forage yield [18]. Thus, the stem/leaf ratio is regarded as a critical indicator of inferred forage yield. However, Annicchiarico [19] considered forage yield and leaf/stem ratio to be genetically independent, and this is consistent with genetic correlation results for alfalfa cultivars [20]. Thus, there is limited and somewhat inconsistent information on the relationship between the stem/leaf ratio and forage yield of alfalfa. The LAI is a useful indicator that affects forage yield; it has been used indirectly as an input variable for primary production models, crop growth, and yield forecasting models [21,22,23,24,25]. The LI reflects the photosynthetically active radiation intercepted and accumulated during the crop growth cycle which is converted into biomass [26,27]. Therefore, LAI and LI should be considered when investigating forage productivity.
Recently, several studies have assessed alfalfa varieties for their morphological characteristics and forage yield under drought conditions. For instance, among Mediterranean cultivars, Ameristand801S and ABT805, for example, were found to be more drought tolerant than Magali, Melissa, and Coussouls, owing to higher persistence values and autumn dormancy [28]. Other cultivars such as Anchor and WL 318 are considered drought-resistant in the U.S. owing to their greater root length and water retention [29]. In France, the TC cultivar has been planted in arid regions owing to its high water use efficiency, leaf area ratio, specific leaf area, and dry matter content [30]. In China, cultivars such as Xinmu No. 1, Longdong, and Aohan showed a high degree of tolerance to drought stress [31,32,33]. Nevertheless, there is a need to find other varieties that have better performance in order to sustain optimal productivity of foraging agricultural systems in these regions.
The Hexi Corridor is an arid region located in the Gansu Province (north-western China) with abundant light and heat resources for crop production [34,35]. However, the forage growing season is generally short, resulting in a severe seasonal imbalance in forage yield, and challenging livestock production. This is further aggravated by the degradation of natural grassland caused by high livestock density grazing and the lack of research-based planning during grazing. These factors result in low and unstable forage production that cannot meet the yield and nutritional demands of livestock. Therefore, it is necessary to select forage species and varieties with high yield and quality. In 2013, the local government identified alfalfa as a suitable forage option; it was subsequently grown on about 650,000 ha in Gansu province [36]. Currently, the Xinjiang Daye variety is widely grown in this region [37]. However, information on other varieties suitable for planting in the region is insufficient. Thus, selecting more suitable alfalfa varieties is essential for the development of forage and improvement of livestock production.
To understand the conditions of alfalfa growth, we hypothesized that canopy structure traits were critical indicators for forage yield. This study explored the morphological traits and productivity of 10 alfalfa varieties in the Hexi Corridor. The study objectives were to: (1) evaluate the effects of canopy structure traits on forage yield; and (2) identify alfalfa varieties which are most suitable for arid regions. The results will guide alfalfa production and variety development in dry areas.

2. Results

2.1. Climatic Conditions

Precipitation varied in the three study years. In 2017, the annual precipitation was 127 mm, i.e., higher than the long-term average (117 mm; 2005–2020), while in 2016 and 2018 (69 and 99 mm, respectively) it was lower. From 2016 to 2018, precipitation within the growing season (April to October) was 61, 126, and 93 mm, respectively. However, 75–79% of precipitation occurred during four months (June to September) within the experimental years (Figure 1). In terms of pan evaporation, the values were 2116 mm in 2016, 2055 mm in 2017, and 2187 mm in 2018. Overall, the annual pan evaporation was 20 times greater than the annual precipitation between 2016 and 2018.
In 2016, 2017, and 2018, the average air temperatures were 18.4 °C, 18.0 °C, and 17.8 °C, respectively. The maximum average air temperature occurred in July 2016 (33.4 °C), 2017 (31.8 °C), and 2018 (32.8 °C), and the warmest months were June, July, and August. During the growing seasons, the accumulated radiation was 6160 MJ m−2, 6131 MJ m−2, and 6021 MJ m−2 in 2016, 2017, and 2018, respectively. As demonstrated in Figure 1, the maximum monthly cumulative radiation value was recorded in June for 2016 and in May for 2017 and 2018.

2.2. Forage Yield

The interaction between year and variety significantly influenced annual yield and the yield of each cutting (p < 0.001) (Table 1). From 2017 to 2018, a significant increase in performance was observed in WL168HQ (p < 0.001), WL343HQ (p < 0.001), Xinjiang Daye (p < 0.001), Phabulous (p < 0.001), and Xinmu No. 2 (p < 0.001) (Figure 2, Table S1); the rates of increase were 146%, 146%, 76%, 63%, and 20%, respectively (Figure 2). In contrast, we detected an annual yield reduction rate of 41% and 20% in the Gannong No. 3 (p = 0.004) and Wudi (p < 0.001) varieties, respectively. Lastly, there was no significant difference in the annual yields of the Derby, Gannong No. 5, and Sanditi varieties (p = 0.063, 0.197, and 0.353, respectively) between 2017 and 2018.
In 2017, the annual yield of Gannong No. 5 (6924 kg ha−1) was significantly higher than those of all the other alfalfa varieties, followed by Gannong No. 3 (5242 kg ha−1) (Figure 2, p < 0.001). In addition, we observed the lowest annual yield in WL168HQ (2680 kg ha−1) (p < 0.001). In 2018, the annual yield varied from 3057 kg ha−1 in Wudi to 8820 kg ha−1 in WL343HQ (p < 0.001). In terms of the two-year total yield, Gannong No. 5 (13,923 kg ha−1) had the greatest significant value among the 10 varieties, followed by WL343HQ (12,409 kg ha−1), Phabulous (11,928 kg ha−1), Xinjiang Daye (11,416 kg ha−1), and Derby (9762 kg ha−1) (p < 0.001), whereas the lowest value was observed with the Wudi variety (6876 kg ha−1) (p < 0.001).
The forage yield in each cutting decreased with harvest time (Figure 3). In 2017, the proportion of first cutting yield to annual yield of 10 varieties ranged from 42 to 67%, whereas in 2018, it ranged from 35 to 60%. In addition, the first cutting yield of Derby (3163 kg ha−1) was significantly the greatest among the 10 varieties, followed by Gannong No. 5 (2883 kg ha−1), Sanditi (2748 kg ha−1), and Phabulous (2716 kg ha−1) (p < 0.001) in 2017. The second cutting yield varied from 615 kg ha−1 in WL168HQ to 1633 kg ha−1 in Gannong No. 5 (p < 0.001) in 2017. On average, the third and fourth cutting yields were equivalent to 27% (619 kg ha−1) and 20% (457 kg ha−1) of the first cutting yield (2334 kg ha−1) of the 10 alfalfa varieties in 2017. In 2018, the top two yields of the first cutting were WL343HQ (5251 kg ha−1) and Xinjiang Daye (3648 kg ha−1) varieties, whereas the lowest value was observed with Gannong No. 3 (1694 kg ha−1) and Wudi (1657 kg ha−1) varieties (p < 0.001). Here, we observed that the second cutting yield of Xinjiang Daye (2318 kg ha−1) was higher than those of the other nine varieties, while the lowest values occurred with Sanditi (713 kg ha−1) (p < 0.001).

2.3. Plant Height

Plant height at spring (PHS) and at flowering (PHF) were measured in this study. These parameters were found to be significantly affected by the varieties and the interaction between years and varieties (p < 0.001) (Table 1). In 2017, the PHSs of Gannong No. 5 (26.6 cm) and Derby (25.3 cm) were significantly greater than those of other alfalfa varieties, whereas the lowest values were observed in WL168HQ (17.6 cm) and Wudi (17.8 cm). In 2018, the greatest and lowest PHSs occurred in Gannong No. 5 (29.2 cm) and Sanditi (16.8 cm), respectively. Overall, Phabulous and WL343HQ presented a significant increase in PHS from 2017 to 2018 (Table S1, p = 0.047, and =0.012, respectively). In contrast, Sanditi showed a significant decrease (p = 0.02).
In terms of PHF, the values of Gannong No. 5, Phabulous, WL168HQ, WL343HQ, and Wudi were significantly higher in 2018 than in 2017 (Table S1, p = 0.008, 0.039, 0.007, 0.014, and 0.045, respectively). However, the PHF of Sanditi was significantly lower (p = 0.013). There were no significant differences for Derby, Gannong No. 3, Xinjiang Daye, and Xinmu No. 2 between years (Figure 4, Table S1, p = 0.408, 0.871, 0.621, and 0.399, respectively). In 2017, the PHF of Wudi (67.2 cm) was considerably lower than those of the other nine alfalfa varieties (p < 0.001). In 2018, the lowest and greatest PHFs occurred in Sanditi (67.9 cm) and Gannong No. 5 (98.5 cm), respectively (p < 0.001). Additionally, the PHF of Wudi (75.5 cm) was lower than those of Gannong No. 3 (83.3 cm), Phabulous (88.1 cm), WL168HQ (83.5 cm), and WL343HQ (87.8 cm) (p < 0.001).

2.4. Branch Number

As illustrated in Figure 5, the number of branches varied among the 10 varieties in 2017. The primary branch number per plant (PB) of Sanditi (18.0) was higher than those of WL168HQ (12.7), Wudi (13.0), and WL343HQ (12.0) (p = 0.079). The secondary branch number per stem (SB) of Derby (19.0) and Wudi (19.0) were significantly greater than that of WL168HQ (16.3) (p = 0.131). Lastly, the total branch number per plant (TB) of Derby (311), Gannong No. 3 (306), and Sanditi (324) were significantly higher than those of WL168HQ (205) and WL343HQ (205) (p = 0.026).

2.5. Leaf Area Index and Light Interception

The LAI varied significantly among the 10 alfalfa varieties (p < 0.001) (Figure 6). The LAI of Xinjiang Daye, Gannong No. 5, and Xinmu No. 2 (2.62, 2.43, and 2.32, respectively) were significantly higher than those of Derby (1.78), Gannong No. 3 (1.56), Sanditi (1.42), and Wudi (1.65) (p < 0.001). Alternatively, the LI of the Xinjiang Daye varieties was significantly greater than that of Xinmu No. 2 (p < 0.001). However, the LI of Gannong No. 3 and Sanditi (0.56 and 0.57, respectively) were significantly lower than those of Gannong No. 5 (0.76), Phabulous (0.71), WL343HQ (0.68), Xinmu No. 2 (0.72), and Xinjiang Daye (0.82) (p < 0.001).

2.6. Stem to Leaf Ratio

As shown in Figure 7, the stem to leaf ratio varied among the 10 varieties in 2017. The stem/leaf ratio of Gannong No. 3 (3.69) was significantly greater than those of the other eight alfalfa varieties, except Phabulous (p < 0.001). The lowest ratio was observed in WL168HQ (1.71) (p < 0.001).

2.7. The Central Leaflet Length, Width, and Area

The leaf size attributes (length, width, and area) differed significantly among the 10 varieties (p = 0.035, 0.013, and <0.001, respectively) (Figure 7). The central leaflet length varied from 24.0 mm (WL343HQ) to 28.3 mm (Derby) (p = 0.035). Its width varied from 8.7 mm (Wudi) to 10.8 mm (Sanditi) (p = 0.013). The central leaflet area of Xinmu No. 2 had the greatest value (208 mm2) among 10 alfalfa varieties (p < 0.001). Moreover, the area of central leaflet of WL343HQ (150 mm2) was significantly lower than those of other alfalfa varieties (p < 0.001).

2.8. Correlation Analysis for First Cutting Yield

Considering that forage yield was influenced by the morphological traits of the plants, we analyzed correlations between forage yield and the above-mentioned traits using data (sample size of 30 for each trait) in the same growing year (2017), which are summarized in the correlation coefficient matrix (Table 2) and PCA biplot based on the correlation matrix (Figure 8). Before starting PCA, the applicability of PCA was tested using the Kaiser–Meyer–Olkin (KMO) and Barlett tests. The calculated results were KMO = 0.66 (>0.5) and Barlett test value = 0 (<0.05), indicating that the dataset was appropriate for a factor analysis.
The first two principal components explained 70.8% of the first yield variation (Figure 8). Thus, we focused on these two components to represent the PCA ordination biplot loadings. The angles between the vectors represent the correlations, where small angles (close to 0°) indicate high positive correlations; wide angles (close to 180°) indicate high negative correlations; and angles close to 90° indicate a lack of correlation [38]. As illustrated in Figure 8 and Table 2, PCA and correlation analysis showed that the yield of first cutting had a significant positive association with PHS, PHF, and TB (p = 0.006, 0.002, and <0.001, respectively) in 2017. The total yield of 2017 was strongly related to PHS, PHF, and stem/leaf ratio (p < 0.001, =0.005, and <0.001, respectively). The two-year total yield was significant positive affected by PHF, LAI, and LI (p < 0.001, =0.002, and <0.001, respectively). Moreover, there was a highly significant positive relationship between LAI and LI (r = 0.8745, p < 0.001) (Table 2). The stem/leaf ratio showed a significant positive association with PHF and PB (p < 0.001 and =0.027, respectively).

3. Discussion

Stem traits, including plant height (PH), branch number, and stem to leaf ratio, are essential indicators related to the growth and development of forage, and are associated with the distribution of solar radiation in the canopy of forage crops [39]. In this study, PH (including PHS and PHF) was shown to significantly affect forage yield, which is consistent with previous studies [9,20,40]. In addition, PH demonstrated a positive correlation to the central leaflet area and stem to leaf ratio, but no effects on LAI and LI.
Branch number is an important morphological characteristic that could enhance the forage yield in alfalfa grown in arid environments [16]. In this study, the results indicated that PB was significantly corelated with forage yield, whereas SB had no effect on yield. Additionally, SB expressed less internal variation among the 10 alfalfa varieties. These results demonstrated that PB plays a greater role in plant architecture construction, determining forage yield. These findings are consistent with those of previous studies [41]. Therefore, we suggest that PB should be viewed as a crucial indicator of high productivity varieties for plant breeding research in the future.
The PCA and correlation analysis results revealed that a higher stem/leaf ratio indicates a better yield performance. This could be attributed to a higher stem/leaf ratio accompanied by higher PH and PB. These results are consistent with those reported by other researchers [17,18,42].
The leaf is considered to be an important functional unit of a plant, contributing to its yield formation. Traditionally, LAI and LI demonstrated a positive correlation with forage yield [25]. In this study, LAI and LI had a significant positive correlation with two-year total yield. However, they did not affect the yield of first cutting and total yield from four cuttings in 2017. The first reason could be that the effects of leaflet density and number to the LAI and LI were more significant than the central leaflet size in an arid environment [43]. Thus, in arid environments, the positive yield effect of canopy structure traits should not be attributed simply to the leaf traits, such as LAI and LI; rather, it could also be due to its effects of stem traits, like PH and branch number, via an increase in the stem dry weight ratio, as proposed by Hakl (2021) [7] and Ta (2020) [8]. Secondly, forage yield is related to light interception that accumulates throughout the entire growth cycle. Therefore, only one LAI and LI measurement before harvest does not truly reflect the efficiency of light energy utilization. It is suggested that researchers consider the significance of dynamic and continuous measurements of LAI and LI in future experiments.
Choosing an alfalfa variety is arguably one of the most important decisions to make in an alfalfa operation. Alfalfa yield reveals the suitability of a given variety to arid regions. On average, the forage annual yield in the second harvest year increased by 30% (1326 kg ha−1) compared with that of the first harvest year. Similarly, a study showed that the forage yield of alfalfa reached peak performance in the second and third harvest years [44]. Therefore, the advantage of accumulation of perennial alfalfa biomass appeared in the second harvest year. From our perspective, we should not only consider the yield of first harvest year, but also pay close attention to long-term productivity.
As one of the most widely planted varieties in the Hexi Corridor, Xinjiang Daye showed good forage yield in 2017 and 2018. In addition, each yield of its four cuttings significantly increased between the two growing years. In addition, it showed outstanding performance regarding the accumulation of light energy owing to its high leaf area index and excellent light interception.
Three varieties, i.e., Gannong No. 5, WL343HQ, and Phabulous, performed better in terms of their two-year total yield than Xinjiang Daye. Gannong No. 5 gave the highest value, i.e., 22% (2507 kg ha−1) higher than that of Xinjiang Daye. This could be due to its excellent PH, LAI, LI, central leaflet area, and stem to leaf ratio. WL343HQ presented the most significant increase in forage yield from 2017 to 2018 among the 10 alfalfa varieties, and reached the highest yield in 2018. The two-year total yield of WL343HQ was more than 9% (993 kg ha−1) higher than that of Xinjiang Daye. WL343HQ had a low PH, low PB, low stem/leaf ratio and low central leaflet area, implying that this cultivar has substantial nutritional potential. A higher leaf/stem ratio has been largely used as a positive indicator of alfalfa forage quality in previous studies [19], owing to its strict positive correlation with digestibility and intake of forage [18] that arises from the higher digestibility of leaves relative to stems [45,46]. The two-year yield of Phabulous was 4% (512 kg ha−1) more than that of Xinjiang Daye. Phabulous possessed a high potential for forage yield because of its high stem/leaf ratio and PH features. However, the two-year total yield of Derby was slightly below 14% (1653 kg ha−1) compared with that of Xinjiang Daye. Considering its high branch number, PH, and the central leaflet area, Derby could be used as a supplementary variety to reduce forage shortages for animals.
In 2017 and 2018, Wudi displayed poor performances in its annual yield. Gannong No. 3 showed poor productivity, significantly reducing each cutting yield from 2017 to 2018. Although the yield of Xinmu No. 2 and WL168HQ increased from 2017 to 2018, they presented poor performances in their two-year total yields. Moreover, Wudi, Sanditi, Gannong No. 3, Xinmu No. 2, and WL168HQ showed poor over-wintering performances (data not shown) in this study. When winter hardiness is poor for alfalfa, the productivity is reduced in an arid environment [47]. Generally, the yield should be the greatest determining factor for variety selection. In this study, Wudi, Gannong No. 3, Xinmu No. 2, Sanditi, and WL168HQ were found to be the least suitable varieties and are not recommended for planting in the Hexi Corridor region.
The study evaluated the correlation of forage yield with canopy structure traits and analyzed the productivity of 10 alfalfa varieties in an arid environment. Our results suggest that we should focus more intensively on stem traits, such as PH, PB, and stem/leaf ratio, which contribute significantly to canopy structure, and consequently, produce more forage yield in an arid environment. Moreover, Gannong No. 5, WL343HQ, Phabulous, and Xinjiang Daye achieved the highest forage yields. However, more field experiments with multiple locations are necessary to test these varieties over a broader area. In addition, a closer investigation of forage quality and related indicators of the selected varieties is necessary in follow-up research. These results provide more insights into alfalfa characteristics and could be used to improve forage yield and livestock production sustainability in the world’s ever-expanding arid regions.

4. Materials and Methods

4.1. Site Description

The field experiment was conducted from 2016 to 2018 at the Linze experimental station of Lanzhou University in the Hexi Corridor, Gansu Province, China (100°02′ E, 39°15′ N, elevation 1390 m a.s.l.), a typical temperate desert climate area. We collected the climate datasets from 2005 to 2020. The annual mean air temperature was 7.6 °C, absolute maximum and minimum were 39.1 °C and −27 °C, and average annual precipitation was 117 mm. Moreover, the mean annual pan evaporation was 2390 mm, i.e., 20 times greater than the annual average precipitation. The frost-free season was 165 days, on average. The annual accumulated temperature of ≥10 °C was 3088 °C, on average [48].
The experimental site has a meadow soil mapped as an Entisol according to the U.S. soil classification system [49], a pH of 7.6, 0.806g kg−1 total nitrogen, and 0.703 g kg−1 total phosphorus at sowing in 2016. Previously, the experimental site had been utilized to grow yellow sweet clover [Melilotus officinalis (L.) Pall.] and maize (Zea mays L.).

4.2. Experimental Design and Agronomic Management

As outlined in Table 3, a total of 10 alfalfa cultivars were used in this experiment (Table 3). Alfalfa seeds were sown in nursery substrate (containing more than 40% organic matter and 15% humic acid) and acclimatized in a greenhouse for 45 days. Each plant was then transplanted into the experimental site on 20 June 2016. All experimental plots were arranged in a randomized complete block design with three replicates. The plot size was 3.6 m × 2.4 m, and the seedlings were planted at a spacing of 60 cm within and between rows.
In 2016, the experimental field received rotary tillage to a depth of 0.2 m before sowing. In addition, no fertilizer was applied at sowing or during the growing seasons of this experiment. The field received 150 mm of irrigation on 24 April, 25 June, 25 August, and 25 October during each growing season based on conventional local practices. However, irrigation could not compensate for the water loss caused by high pan evaporation (over 2390 mm per year). Weeds were controlled manually as necessary.

4.3. Measurements

4.3.1. Meteorological Data

Daily precipitation, incoming solar radiation, and maximum and minimum air temperature data were obtained from the Linze Farming Ecosystem experimental weather station of the Chinese Academy of Sciences [50].

4.3.2. Plant Height

Plant Height was measured one month after the spring emergence and at flowering (first cutting). Three plants were randomly selected from each plot, and their height was measured as the distance from ground to the top end using a ruler (accuracy 1 mm).

4.3.3. Branch Number

During the initial growth phase, a primary axis develops from the seed, and secondary shoots subsequently develop from the axillary buds of cotyledons and primary axes leaves [51]. The primary branch number per plant (PB) and secondary branch number per stem (SB) were counted on five randomly selected plants from each plot at the flowering stage before the first cutting in 2017. Lastly, the total branch number per plant (TB) was calculated as the product of PB and SB.

4.3.4. Leaf Area Index and Light Interception

The leaf area index (LAI) and light interception (LI) were measured at the flowering stage before the first cutting in 2017 during clear sunshine from 11:00 to 13:00 using a Ceptometer (AccuPAR LP-80, Decagon Devices Inc., Pullman, WA, USA). One above and four below the canopy (close to the ground) randomly chosen readings on points were collected from each plot. Notably, the average below canopy reading was used. The LI was calculated [52] as:
L I = 1 P A R b l P A R a b   ,
where PARbl and PARab are photosynthetically active radiation below and above the canopy, respectively.

4.3.5. Forage Yield and Stem to Leaf Ratio

Samples were harvested manually using scythes on the plant base and measured four times annually in 2017 and 2018 at an interval of 4–5 weeks, starting at the initial flowering stage, specifically, 2 June, 8 July, 9 August, and 14 September in 2017, and 12 June, 14 July, 17 August, and 18 September in 2018. In the first cutting of 2017, the stem/leaf ratio was recorded on a single plant basis (after separation and oven drying of the stems and leaves). Subsequently, 10 plants were randomly chosen from each plot, cut at ground level and immediately weighed, before being oven-dried at 64 °C for 48 h to a constant weight in order to determine the dry weight.

4.3.6. The Central Leaflet Length, Width, and Area

The central leaflet was a mature leaf taken from the main stem. To measure the length, width, and area, 20 randomly selected central leaflets from each plot (each species a total of 60 leaflets) were scanned, and the image was analyzed with an Epson GT-15000 and imaging software (WinSEEDLETM, Quebec City, QC, Canada).

4.4. Statistical Analysis

All the statistical analyses described below were performed using the R software (version 4.0.3). Plots were produced using the ggplot2 [53] and Factoextra [54] packages. A one-way analysis of variance (ANOVA) was conducted to assess differences of all measured traits among the 10 alfalfa varieties during the experimental years. The main effects and interactions of year and cultivar were obtained using two-way ANOVA. All tests were considered significantly different at p < 0.05, and the means were compared using Duncan’s multiple comparison test. In addition, an independent samples t-test was conducted to assess the differences of traits between 2017 and 2018. Finally, we evaluated the relationships between canopy traits and forage yield using principal component analysis (PCA) based on a correlation matrix and Pearson’s correlation analysis. Before starting PCA, we checked the sampling adequacy to detect whether or not the data would factor well. Sampling adequacy was measured by the Kaiser–Meyer–Olkin criterion (KMO). KMO varied from 0 to 1 and the overall KMO had to be not less than 0.5 to proceed to the principal component analysis. The Bartlett Test is another criterion in this sampling procedure which is likely also significant [55].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11091112/s1, Table S1: p-values for independent sample t-tests of 10 alfalfa varieties between 2017 and 2018.

Author Contributions

Conceptualization, X.J.; methodology, X.J.; formal analysis, X.J.; investigation, X.J.; resources, X.J.; data curation, X.J.; writing—original draft preparation, X.J. and Z.Z.; writing—review and editing, X.J. and Z.Z.; visualization, X.J.; supervision, Y.W.; project administration, Y.W.; funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Strategic Priority Science and Technology Project of Chinese Academy of Sciences (XDA26030103) and the National Regional Test of Grass Varieties of China (21301060001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data supporting this study are included in the article.

Acknowledgments

The authors thank Yafeng Huang for providing the AccuPAR LP-80 plant canopy analyzer, and Jeffrey A. Coulter, Xiang Lu, Xueli Niu, Rui Dong, and Jianfeng Wang for providing helpful feedback on a previous draft of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Salinger, J.; Sivakumar, M.V.; Motha, R.P. Increasing Climate Variability and Change: Reducing the Vulnerability of Agriculture and Forestry; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2005; Volume 70, pp. 1–3. [Google Scholar]
  2. Yin, Y.; Ma, D.; Wu, S. Enlargement of the semi-arid region in China from 1961 to 2010. Clim. Dyn. 2019, 52, 509–521. [Google Scholar] [CrossRef]
  3. Dixon, J.A.; James, D.E.; Sherman, P.B. Economics of Dryland Management; Routledge: London, UK, 2013; Volume 2, p. 27. [Google Scholar]
  4. Blum, A. Effective use of water (EUW) and not water-use efficiency (WUE) is the target of crop yield improvement under drought stress. Field Crops Res. 2009, 112, 119–123. [Google Scholar] [CrossRef]
  5. Arshad, M.; Feyissa, B.A.; Amyot, L.; Aung, B.; Hannoufa, A. MicroRNA156 improves drought stress tolerance in alfalfa (Medicago sativa) by silencing SPL13. Plant Sci. 2017, 258, 122–136. [Google Scholar] [CrossRef] [PubMed]
  6. Avci, M.A.; Ozkose, A.; Tamkoc, A. Determination of yield and quality characteristics of alfalfa (Medicago sativa L.) varieties grown in different locations. J. Anim. Vet. Adv. 2013, 12, 487–490. [Google Scholar] [CrossRef]
  7. Hakl, J.; Kunzová, E.; Tocauerová, Š.; Menšík, L.; Mrázková, M.; Pozdíšek, J. Impact of long-term manure and mineral fertilization on yield and nutritive value of lucerne (Medicago sativa) in relation to changes in canopy structure. Eur. J. Agron. 2021, 123, 126219. [Google Scholar] [CrossRef]
  8. Ta, H.T.; Teixeira, E.I.; Brown, H.E.; Moot, D.J. Yield and quality changes in lucerne of different fall dormancy ratings under three defoliation regimes. Eur. J. Agron. 2020, 115, 126012. [Google Scholar] [CrossRef]
  9. Payero, J.; Neale, C.; Wright, J. Comparison of eleven vegetation indices for estimating plant height of alfalfa and grass. Appl. Eng. Agric. 2004, 20, 385. [Google Scholar] [CrossRef] [Green Version]
  10. Chen, B.; Xu, K.; Li, J.; Li, F.; Qiao, J.; Li, H.; Gao, G.; Yan, G.; Wu, X. Evaluation of yield and agronomic traits and their genetic variation in 488 global collections of Brassica napus L. Genet. Resour. Crop Evol. 2014, 61, 979–999. [Google Scholar] [CrossRef]
  11. Li, F.; Chen, B.; Xu, K.; Gao, G.; Yan, G.; Qiao, J.; Li, J.; Li, H.; Li, L.; Xiao, X.; et al. A genome-wide association study of plant height and primary branch number in rapeseed (Brassica napus). Plant Sci. 2016, 242, 169–177. [Google Scholar] [CrossRef]
  12. Agudamu; Yoshihira, T.; Shiraiwa, T. Branch development responses to planting density and yield stability in soybean cultivars. Plant Prod. Sci. 2016, 19, 331–339. [Google Scholar] [CrossRef] [Green Version]
  13. Ramesh, R.; Ramprasad, E. Effect of plant growth regulators on morphological, physiological and biochemical parameters of soybean (Glycine max L. Merrill). In Biotechnology and Bioforensics; Springer: New York, NY, USA, 2015; Volume 7, pp. 61–71. [Google Scholar]
  14. Wang, X.; Li, L.; Yang, Z.; Zheng, X.; Yu, S.; Xu, C.; Hu, Z. Predicting rice hybrid performance using univariate and multivariate GBLUP models based on North Carolina mating design II. Heredity 2016, 118, 302–310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Pembleton, K.G.; Rawnsley, R.P.; Donaghy, D.J.; Volenec, J.J. Water deficit alters canopy structure but not photosynthesis during the regrowth of alfalfa. Crop Sci. 2009, 49, 722–731. [Google Scholar] [CrossRef]
  16. Volenec, J.J.; Cherney, J.H.; Johnson, K.D. Yield components, plant morphology, and forage quality of alfalfa as influenced by plant population. Crop Sci. 1987, 27, 321–326. [Google Scholar] [CrossRef]
  17. Ray, I.M.; Townsend, M.S.; Muncy, C.M.; Henning, J.A. Heritabilities of water-use efficiency traits and correlations with agronomic traits in water-stressed alfalfa. Crop Sci. 1999, 39, 494–498. [Google Scholar] [CrossRef]
  18. Kephart, K.D.; Buxton, D.R.; Hill, R.R. Digestibility and cell-wall components of alfalfa following selection for divergent herbage lignin concentration. Crop Sci. 1990, 30, 207–212. [Google Scholar] [CrossRef]
  19. Annicchiarico, P. Alfalfa forage yield and leaf/stem ratio: Narrow-sense heritability, genetic correlation, and parent selection procedures. Euphytica 2015, 205, 409–420. [Google Scholar] [CrossRef]
  20. Julier, B.; Huyghe, C.; Ecalle, C. Within- and among-cultivar genetic variation in alfalfa: Forage quality, morphology, and yield. Crop Sci. 2000, 40, 365–369. [Google Scholar] [CrossRef] [Green Version]
  21. Brisson, N.; Gary, C.; Justes, E.; Roche, R.; Mary, B.; Ripoche, D.; Zimmer, D.; Sierra, J.; Bertuzzi, P.; Burger, P. An overview of the crop model STICS. Eur. J. Agron. 2003, 18, 309–332. [Google Scholar] [CrossRef]
  22. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. AquaCrop—The FAO crop model to simulate yield response to water: I. Concepts and underlying principles. Agron. J. 2009, 101, 426–437. [Google Scholar] [CrossRef] [Green Version]
  23. Fang, H.; Liang, S.; Hoogenboom, G. Integration of MODIS LAI and vegetation index products with the CSM–CERES–Maize model for corn yield estimation. Int. J. Remote Sens. 2011, 32, 1039–1065. [Google Scholar] [CrossRef]
  24. Bolton, D.K.; Friedl, M.A. Forecasting crop yield using remotely sensed vegetation indices and crop phenology metrics. Agric. For. Meteorol. 2013, 173, 74–84. [Google Scholar] [CrossRef]
  25. Kross, A.; McNairn, H.; Lapen, D.; Sunohara, M.; Champagne, C. Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops. Int. J. Appl. Earth Obs. Geoinf. 2015, 34, 235–248. [Google Scholar] [CrossRef] [Green Version]
  26. Mattera, J.; Romero, L.A.; Cuatrín, A.L.; Cornaglia, P.S.; Grimoldi, A.A. Yield components, light interception and radiation use efficiency of lucerne (Medicago sativa L.) in response to row spacing. Eur. J. Agron. 2013, 45, 87–95. [Google Scholar] [CrossRef]
  27. Monteith, J.L. Climate and the efficiency of crop production in Britain. Philos. Trans. R. Soc. Lond. B 1977, 281, 277–294. [Google Scholar] [CrossRef]
  28. Bouizgaren, A.; Farissi, M.; Ghoulam, C.; Kallida, R.; Faghire, M.; Barakate, M.; Al Feddy, M.N. Assessment of summer drought tolerance variability in Mediterranean alfalfa (Medicago sativa L.) cultivars under Moroccan fields conditions. Arch. Agron. Soil Sci. 2013, 59, 147–160. [Google Scholar] [CrossRef]
  29. Carter, P.R.; Sheaffer, C.C.; Voorhees, W.B. Root growth, herbage yield, and plant water status of alfalfa cultivars. Crop Sci. 1982, 22, 425–427. [Google Scholar] [CrossRef]
  30. Erice, G.; Louahlia, S.; Irigoyen, J.J.; Sanchez Diaz, M.; Avice, J.C. Biomass partitioning, morphology and water status of four alfalfa genotypes submitted to progressive drought and subsequent recovery. J. Plant Physiol. 2010, 167, 114–120. [Google Scholar] [CrossRef] [PubMed]
  31. Ni, Y.; Guo, Y.; Han, L.; Tang, H.; Conyers, M. Leaf cuticular waxes and physiological parameters in alfalfa leaves as influenced by drought. Photosynthetica 2012, 50, 458–466. [Google Scholar] [CrossRef]
  32. Quan, W.; Liu, X.; Wang, H.; Chan, Z. Comparative physiological and transcriptional analyses of two contrasting drought tolerant alfalfa varieties. Front. Plant Sci. 2016, 6, 1256. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, W.; Kim, Y.; Lee, H.; Kim, K.; Deng, X.; Kwak, S. Analysis of antioxidant enzyme activity during germination of alfalfa under salt and drought stresses. Plant Physiol. Biochem. 2009, 47, 570–577. [Google Scholar] [CrossRef]
  34. Li, X.; Tong, L.; Niu, J.; Kang, S.; Du, T.; Li, S.; Ding, R. Spatio-temporal distribution of irrigation water productivity and its driving factors for cereal crops in Hexi Corridor, Northwest China. Agric. Water Manag. 2017, 179, 55–63. [Google Scholar] [CrossRef]
  35. Bao, C.; Fang, C. Water resources constraint force on urbanization in water deficient regions: A case study of the Hexi Corridor, arid area of NW China. Ecol. Econ. 2007, 62, 508–517. [Google Scholar] [CrossRef]
  36. Wen, Z.; Duan, T.; Christensen, M.J.; Nan, Z. Microdochium tabacinum, confirmed as a pathogen of alfalfa in Gansu province, China. Plant Dis. 2015, 99, 87–92. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Luo, Y.; Li, G. The effect of water stress on growth and biomass of Medicago sativa cv. Xinjiangdaye. Acta Pratacult. Sin. 2014, 23, 213–219. [Google Scholar] [CrossRef]
  38. Avola, G.; Riggi, E.; Gresta, F.; Sortino, O.; Onofri, A. Random effects models, BLUPs and redundancy analyses for grain legume crops in semi-arid environments. Eur. J. Agron. 2018, 93, 18–26. [Google Scholar] [CrossRef]
  39. Keating, B.A.; Carberry, P.S. Resource capture and use in intercropping: Solar radiation. Field Crops Res. 1993, 34, 273–301. [Google Scholar] [CrossRef]
  40. Monirifar, H. Path analysis of yield and quality traits in alfalfa. Not. Bot. Horti Agrobot. Cluj-Na. 2011, 39, 190–195. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, F.; Kang, J.; Long, R.; Yu, L.-X.; Wang, Z.; Zhao, Z.; Zhang, T.; Yang, Q. High-density linkage map construction and mapping QTL for yield and yield components in autotetraploid alfalfa using RAD-seq. BMC Plant Biol. 2019, 19, 165. [Google Scholar] [CrossRef] [Green Version]
  42. Avci, M.; Hatipoglu, R.; Çinar, S.; Kiliçalp, N. Assessment of yield and quality characteristics of alfalfa (Medicago sativa L.) cultivars with different fall dormancy rating. Legume Res. 2017, 41, 369–373. [Google Scholar] [CrossRef] [Green Version]
  43. Baldissera, T.C.; Frak, E.; Carvalho, P.C.; Louarn, G. Plant development controls leaf area expansion in alfalfa plants competing for light. Ann. Bot. 2014, 113, 145–157. [Google Scholar] [CrossRef] [Green Version]
  44. Theron, J.F.; Snyman, H.A. Productivity evaluation of Medicago sativa cultivars under irrigation in a semi-arid climate. Afr. J. Range Forage Sci. 2015, 32, 161–171. [Google Scholar] [CrossRef]
  45. Albrecht, K.A.; Wedin, W.F.; Buxton, D.R. Cell-wall composition and digestibility of alfalfa stems and leaves. Crop Sci. 1987, 27, 735–741. [Google Scholar] [CrossRef]
  46. Marten, G.; Buxton, D.; Barnes, R. Feeding value (forage quality). In Alfalfa and Alfalfa Improvement; Amer Society of Agronomy: Madison, WI, USA, 1988; Volume 14, pp. 463–491. [Google Scholar]
  47. Castonguay, Y.; Laberge, S.; Brummer, E.C.; Volenec, J.J. Alfalfa Winter Hardiness: A Research Retrospective and Integrated Perspective. In Advances in Agronomy; Academic Press: Cambridge, MA, USA, 2006; Volume 90, pp. 203–265. [Google Scholar]
  48. Su, P.; Liu, X.; Zhang, L.; Zhao, A.; Li, W.; Chen, H. Comparison of δ 13C values and gas exchange of assimilating shoots of desert plants Haloxylon ammodendron and Calligonum mongolicum with other plants. Isr. J. Plant Sci. 2004, 52, 87–97. [Google Scholar] [CrossRef]
  49. Su, Y.; Wang, J.; Yang, R.; Yang, X.; Fan, G. Soil texture controls vegetation biomass and organic carbon storage in arid desert grassland in the middle of Hexi Corridor region in Northwest China. Soil Res. 2015, 53, 366–376. [Google Scholar] [CrossRef]
  50. Linze Farming Ecosystem Experimental Station. Linze Inland River Basin Research Station, Gansu Province. National Science & Technology Infrastructure. Available online: http://lzd.cern.ac.cn/meta/metaData (accessed on 1 February 2021).
  51. Hanson, A.A.; Barnes, D.K.; Hill, R.R.; Fick, G.W.; Holt, D.A.; Lugg, D.G. Environmental physiology and crop growth. In Alfalfa and Alfalfa Improvement; Amer Society of Agronomy: Madison, WI, USA, 1988; Volume 5, pp. 163–194. [Google Scholar]
  52. Purcell, L.C. Soybean canopy coverage and light interception measurements using digital imagery. Crop Sci. 2000, 40, 834–837. [Google Scholar] [CrossRef]
  53. Wickham, H. ggplot2: Create Elegant Data Visualisations Using the Grammar of Graphics. Version 3.3.5. 25 June. 2021. Available online: https://cran.r-project.org/web/packages/ggplot2/index.html (accessed on 17 February 2022).
  54. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. Version 1.0.7. 1 April. 2020. Available online: https://cran.r-project.org/web/packages/factoextra/index.html (accessed on 17 February 2022).
  55. Wania, A.; Kühn, I.; Klotz, S. Plant richness patterns in agricultural and urban landscapes in Central Germany—Spatial gradients of species richness. Landsc. Urban Plan. 2006, 75, 97–110. [Google Scholar] [CrossRef]
Figure 1. Monthly maximum (Tmax), and minimum (Tmin) average air temperature (°C), precipitation (mm), pan evaporation (mm), and cumulative radiation (MJ m−2) at the experimental station in 2016–2018.
Figure 1. Monthly maximum (Tmax), and minimum (Tmin) average air temperature (°C), precipitation (mm), pan evaporation (mm), and cumulative radiation (MJ m−2) at the experimental station in 2016–2018.
Plants 11 01112 g001
Figure 2. The annual yield of 10 alfalfa varieties in 2017 and 2018, and the two-year total yield. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Figure 2. The annual yield of 10 alfalfa varieties in 2017 and 2018, and the two-year total yield. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Plants 11 01112 g002
Figure 3. Forage yield for each cutting (1st to 4th) of 10 alfalfa varieties in 2017 and 2018. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Figure 3. Forage yield for each cutting (1st to 4th) of 10 alfalfa varieties in 2017 and 2018. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Plants 11 01112 g003
Figure 4. The plant heights of 10 alfalfa varieties at spring and at flowering of first cutting in 2017 and 2018. Different letters indicate significantly different means at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Figure 4. The plant heights of 10 alfalfa varieties at spring and at flowering of first cutting in 2017 and 2018. Different letters indicate significantly different means at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Plants 11 01112 g004
Figure 5. The Primary Branch number per plant (PB), Secondary Branch number per stem (SB), and Total Branch number per plant (TB) of 10 alfalfa varieties at the flowering stage before the first cutting in 2017. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Figure 5. The Primary Branch number per plant (PB), Secondary Branch number per stem (SB), and Total Branch number per plant (TB) of 10 alfalfa varieties at the flowering stage before the first cutting in 2017. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Plants 11 01112 g005
Figure 6. Leaf area index (LAI) and light interception (LI) of 10 alfalfa varieties at the flowering (first cutting) in 2017. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Figure 6. Leaf area index (LAI) and light interception (LI) of 10 alfalfa varieties at the flowering (first cutting) in 2017. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Plants 11 01112 g006
Figure 7. Stem/Leaf Ratio and central Leaflet Length, Width, and Area at flowering (first cutting) in 2017. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Figure 7. Stem/Leaf Ratio and central Leaflet Length, Width, and Area at flowering (first cutting) in 2017. Different letters indicate means that are significantly different at p < 0.05. Error bars indicate the standard deviation of the means. DB: Derby, G3: Gannong No. 3, G5: Gannong No. 5, PA: Phabulous, SA: Sanditi, W1: WL168HQ, W3: WL343HQ, WD: Wudi, X2: Xinmu No. 2, XD: Xinjiang Daye.
Plants 11 01112 g007
Figure 8. Principal component analysis of alfalfa growth traits and forage yield. PHS: plant height at spring, PHF: plant height at flowering of first cutting, TB: total branch number per plant, LAI: leaf area index, LI: light interception, S/L: stem to leaf ratio, Y1C1: first cutting yield in 2017, Y1T: total yield from four cuttings in 2017, YT: two-year total yield.
Figure 8. Principal component analysis of alfalfa growth traits and forage yield. PHS: plant height at spring, PHF: plant height at flowering of first cutting, TB: total branch number per plant, LAI: leaf area index, LI: light interception, S/L: stem to leaf ratio, Y1C1: first cutting yield in 2017, Y1T: total yield from four cuttings in 2017, YT: two-year total yield.
Plants 11 01112 g008
Table 1. The two-way ANOVA p-values for main effects and interactions of year and cultivar on alfalfa yield and plant height.
Table 1. The two-way ANOVA p-values for main effects and interactions of year and cultivar on alfalfa yield and plant height.
DFPHSPHF1st Yield2nd Yield3rd Yield4th YieldAnnual Yield
Year10.2240.125<0.001<0.001<0.001<0.001<0.001
Variety9<0.001<0.001<0.001<0.001<0.001<0.001<0.001
Year × variety9<0.001<0.001<0.001<0.001<0.001<0.001<0.001
Note: DF: degrees of freedom, PHS: plant height at spring, PHF: plant height at flowering of first cutting, 1st Yield: first cutting yield, 2nd Yield: second cutting yield, 3rd Yield: third cutting yield, 4th Yield: fourth cutting yield, Annual Yield: total yield from four cuttings.
Table 2. Correlation coefficient matrix among morphological traits of 10 alfalfa varieties.
Table 2. Correlation coefficient matrix among morphological traits of 10 alfalfa varieties.
PHSPHFPBSBLAILIS/LY1C1Y1TYT
PHS1
PH0.577 ***1
PB0.0180.3001
SB0.125−0.0110.1001
LAI0.3860.242−0.213−0.3431
LI0.4250.206−0.119−0.1570.875 ***1
S/L0.2950.574 ***0.403 *0.129−0.023−0.0391
Y1C10.493 **0.548 **0.55 **0.273−0.070.0180.4621
Y1T0.648 ***0.495 **0.3320.2510.0670.1370.652 ***0.702 ***1
YT0.587 ***0.45−0.092−0.1920.552 **0.576 ***0.1110.320.4551
Note: *, ** and *** denote significance difference at p < 0.05, 0.01 and 0.001, respectively. PHS: plant height at spring, PHF: plant height at flowering of first cutting, PB: primary branch number per plant, SB: secondary branch number per stem, LAI: leaf area index, LI: light interception, S/L: stem leaf ratio, Y1C1: first cutting yield in 2017, Y1T: total yield from four cuttings in 2017, YT: two-year total yield.
Table 3. Details of 10 alfalfa varieties used in this study.
Table 3. Details of 10 alfalfa varieties used in this study.
NameTypeCountry of OriginInstitution-Provided Seed
DerbyCultivarUnited KingdomLanzhou University
Gannong No. 3CultivarChinaGansu Agriculture University
Gannong No. 5CultivarChinaGansu Agriculture University
PhabulousCultivarUnited StatesUSDA-NPGS
SanditiCultivarFranceLanzhou University
WL168HQCultivarUnited StatesUSDA-NPGS
WL343HQCultivarUnited StatesUSDA-NPGS
WudiLand raceChinaLanzhou University
Xinjiang DayeLand raceChinaLanzhou University
Xinmu No. 2CultivarChinaLanzhou University
Note: USDA-NPGS: United States Department of Agriculture, National Plant Germplasm System.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jia, X.; Zhang, Z.; Wang, Y. Forage Yield, Canopy Characteristics, and Radiation Interception of Ten Alfalfa Varieties in an Arid Environment. Plants 2022, 11, 1112. https://doi.org/10.3390/plants11091112

AMA Style

Jia X, Zhang Z, Wang Y. Forage Yield, Canopy Characteristics, and Radiation Interception of Ten Alfalfa Varieties in an Arid Environment. Plants. 2022; 11(9):1112. https://doi.org/10.3390/plants11091112

Chicago/Turabian Style

Jia, Xitao, Zhixin Zhang, and Yanrong Wang. 2022. "Forage Yield, Canopy Characteristics, and Radiation Interception of Ten Alfalfa Varieties in an Arid Environment" Plants 11, no. 9: 1112. https://doi.org/10.3390/plants11091112

APA Style

Jia, X., Zhang, Z., & Wang, Y. (2022). Forage Yield, Canopy Characteristics, and Radiation Interception of Ten Alfalfa Varieties in an Arid Environment. Plants, 11(9), 1112. https://doi.org/10.3390/plants11091112

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop