The Responses of Stem and Leaf Functional Traits of Medicago sativa and Bromus inermis to Different Mixed Planting Patterns

: This study investigated the differences in stem and leaf growth characteristics of Medicago sativa and Bromus inermis in the Jiaozhou region of China during 2019–2020 under three different planting modes of the two forages: monoculture, mixed species sowing in the same rows, and mixed species sowing in alternating rows. No special management of the experimental plots was carried out in this study to simulate as much as possible the growth of forages in their natural state. The stem and leaf characteristics inﬂuencing the dry matter weight were calculated using grey correlation. These characteristics included leaf length, leaf width, leaf thickness, leaf area, leaf fresh weight, stem length, stem diameter, stem fresh weight, stem–leaf ratio, fresh matter yield, dry matter yield, and protein yield of M. sativa and B. inermis under different sowing methods in different years. The results showed that the weight pattern of the characteristics affecting the yield of M. sativa and B. inermis production was leaf area > stem diameter > leaf length > stem length > leaf width > leaf thickness, leaf area > leaf length > stem length > leaf width > leaf thickness > stem diameter. Considering all the growth factors, the production capacity was ranked as mixed sowing in alternating rows > mixed sowing in same rows > monoculture. Thus, the suitable mode for M. sativa – B. inermis sowing was mixed sowing in alternating rows.


Introduction
Selecting suitable legumes and grass forages to establish a high-yielding mixed grass system is an important way of solving the current grass-livestock imbalance since mixing perennial grass with legume forage improves the grass production and protein content per unit area [1][2][3][4].It also increases the nitrogen nutrient and organic matter content in the soil [5][6][7]; improves soil fertility [8]; and reduces the application of industrial nitrogen fertilizer, production costs, and environmental pollution [9,10].However, legumes are also known to contribute to the emission of N 2 O, especially when their above-ground residues are incorporated into the soil [11,12].This is why it is extremely important to find a suitable grass-legume combination that minimizes environmental impacts without compromising yield.Alfalfa (Medicago sativa) and awnless brome (Bromus inermis) are forage crops widely grown and utilized in agricultural and pastoral areas of China, with significant advantages, such as high protein and mowing resistance [13,14].The mixed planting of alfalfa and awnless brome exhibit advantages in mixing legume and grass forage [15], such as occupying the different above-ground and below-ground spaces and maximizing the utilization of light resources and below-ground nutrients, thus significantly improving Agronomy 2023, 13, 2733 2 of 15 forage yield and quality [3].Therefore, alfalfa and awnless brome can complement each other at several levels when planted together.
Previous studies mostly used forage nutritional quality and yield as important indicators to assess the forage mixing patterns [16,17], but this approach has some limitations.Since functional trait characteristics reflect species differences in resource competition, resistance to external pressure, and adaptation strategies, studies based on the functional traits of plant organs can better elucidate the interactions and coexistence mechanisms among species [18,19].Therefore, this study analyzed the characteristics of stem and leaf functional traits of alfalfa-awnless brome in mixed planting to reveal the plasticity of these traits and their adaptive responses.We also analyzed the correlation between stem and leaf functional traits of legumes and grasses under mixed sowing, providing a theoretical basis for elucidating the mechanism of efficient production in legume-grass mixed sowing systems.

The Study Area
This study was conducted in Jiaozhou City, Shandong Province (36 • 26 22 N, 120 • 04 43 E), China.The study area has a warm temperate monsoon climate with a characteristic oceanic climate, an average annual temperature of 12.1 • C (Figure 1), an average annual precipitation of about 695.6 mm, an average annual air pressure of 1015.6 MPa, and an annual frost-free period of 205.5 d, the annual sunshine duration is 2573 h.And the soil type is sandy ginger black, the information of physical and chemical properties of soil at 0~30 cm is in Table 1.   the utilization of light resources and below-ground nutrients, thus significantly imp ing forage yield and quality [3].Therefore, alfalfa and awnless brome can complem each other at several levels when planted together.Previous studies mostly used forage nutritional quality and yield as important cators to assess the forage mixing patterns [16,17], but this approach has some limitat Since functional trait characteristics reflect species differences in resource competition sistance to external pressure, and adaptation strategies, studies based on the functi traits of plant organs can better elucidate the interactions and coexistence mechan among species [18,19].Therefore, this study analyzed the characteristics of stem and functional traits of alfalfa-awnless brome in mixed planting to reveal the plasticity of t traits and their adaptive responses.We also analyzed the correlation between stem leaf functional traits of legumes and grasses under mixed sowing, providing a theore basis for elucidating the mechanism of efficient production in legume-grass mixed ing systems.

The Study Area
This study was conducted in Jiaozhou City, Shandong Province (36°26′22 120°04′43″ E), China.The study area has a warm temperate monsoon climate with a c acteristic oceanic climate, an average annual temperature of 12.1 °C (Figure 1), an ave annual precipitation of about 695.6 mm, an average annual air pressure of 1015.6 M and an annual frost-free period of 205.5 d, the annual sunshine duration is 2573 h. the soil type is sandy ginger black, the information of physical and chemical properti soil at 0~30 cm is in Table 1.

Experimental Design
The perennial leguminous forage alfalfa variety WL525HQ and graminaceous a less brome variety HARANO were obtained from Beijing Zhengdao Seed Industry Ltd., Beijing, China, for this study.The two forages were planted on 22 March 2019,

Experimental Design
The perennial leguminous forage alfalfa variety WL525HQ and graminaceous awnless brome variety HARANO were obtained from Beijing Zhengdao Seed Industry Co., Ltd., Beijing, China, for this study.The two forages were planted on 22 March 2019, in a mixed legume-grass seeding system, and the forage samples were collected at the early flowering stage of alfalfa on 4 June, 10 July, 14 August, and 21 September 2019 and on 3 June, 7 July, and 11 August 2020.A 1 m 2 sampling area was set up in each plot, and all plants within Agronomy 2023, 13, 2733 3 of 15 the sampling area were harvested manually after mowing at a 5 cm stable height.The fourth sampling in the second year (2020) was affected by an epidemic, and the alfalfa was already at the pod stage at the time of collection; thus, only seven crops were harvested across the two years.The awnless brome was not sampled after the first two sampling points to prevent thinning because the growth of the awnless brome was reduced after the first two sampling points each year.The alfalfa and awnless brome were sown individually as controls, while the sowing combination involved alfalfa-awnless brome peer mixing and planting in different rows.There were four treatments: alfalfa single sowing, awnless brome single sowing, alfalfa-awnless brome mixed sowing in alternating rows, and alfalfaawnless brome mixed sowing in the same rows, all of which were sown in strips.
The experimental plots were designed in a randomized block design with three replications of each treatment, with a total of 12 plots.Each plot was 12 m 2 (3 m × 4 m), and all plots were manually furrowed with a row spacing of 20 cm and a sowing depth of 3-4 cm.Based on the seed germination and thousand seed weight of awnless brome and alfalfa, the sowing rate was 37.5 kg/ha and 18.75 kg/ha for awnless brome and alfalfa, respectively, in order to ensure that all the mixed treatments were sown at a plant ratio of 1:1.

Sample Collection 2.3.1. Sample Collection
Fresh forage samples were sampled on 4 June 2019, for the first crop and on 10 July, for the second crop, when alfalfa was in the early flowering stage.Five sample points were selected in each plot for the W-shaped sampling [20,21].Two plants in good growth conditions were selected from each sample point by category, mowed and placed in ice boxes (YMJ-A from Hangzhou Lvbo Instrument Co., Ltd., Hangzhou, China) to be transported to the laboratory for storage at 4 • C.After sampling, three 1 m × 1 m plots were randomly selected in each plot and mowed with 4-5 cm of stubble height, with three rows of awnless brome and two rows of alfalfa in the mixed pattern.The yield of each forage was measured after sorting the mixed forages by species.After that, 500 g of mixed samples were randomly taken from each plot according to the quartering method and transported to the laboratory for further analysis.All plots were mown as a whole after sample collection.

Determination Indices and Methods
Measurement of the fresh samples: Leaf length, leaf width, leaf thickness, and leaf area were measured using a leaf area meter (YMJ-A from Hangzhou Lvbo Instrument Co., Ltd., Hangzhou, China), while the stem diameter was measured using a diameter meter (MA-YH from Nanjing Ming'ao Instrument Co., Ltd., Nanjing, China), and stem length was measured using a steel ruler (MA-YH from Nanjing Ming'ao Instrument Co., Ltd., Nanjing, China).Leaf and stem fresh weights were weighed on an electronic scale (YMJ-A from Hangzhou Lvbo Instrument Co., Ltd., Hangzhou, China).
Dry matter yield measurement: The forage was weighed in bags after mowing and oven-dried at 105 • C for 1 h and then 85 • C to constant weight to calculate the dry matter yield.
Determination of the crude protein yield: The forage samples were dried and crushed with an ultra-micro crusher, after which the total nitrogen content of the pasture was determined with an automatic Kjeldahl nitrogen tester.The crude protein yield of the mixed sown forage was calculated as forage dry matter yield (kg/ha) × forage protein content (%).

Principal Component and Gray Correlation Analyses
Principal component and gray correlation analyses were used to analyze the dry matter yield.Since dry matter yield is a more accurate indicator of grass production capacity [22], the dry matter yield of all treated alfalfa and awnless brome materials was selected as the reference column, which was recorded as {X 0 (k) (k = 1, 2, 3, 4, . .., n).Each index was used as the evaluation index in the comparison column, i.e., the set of observed values of the participating indices, and the indices were recorded as {(k)} (i = 1, 2, 3, 4, . .., m); k = (1, 2, 3, 4, . .., n).The traits used for the analysis included leaf length, leaf width, leaf area, leaf thickness, stem length, stem diameter, crude protein yield, and dry matter yield.The mowed stubble was denoted as X, and the traits were denoted as k; thus, the values of mowed stubble X at trait k formed the comparison series, X i .X 0 was used to construct the reference series.The correlation coefficients were calculated according to the formula (k), while the equal-weighted correlations were calculated using the formula r i .ρ is the resolution coefficient, ρ ∈ [0, 1], where 0.5 represents a moderate positive correlation.Here, we utilized the value 0.5.
Correlation coefficients: Equal-weighted correlations: where n is the number of samples.Weighting factor ω 1 = r i ∑ r i .According to the principal of correlation analysis, the larger the weighted coefficient, the greater the influence of the trait on the evaluation index, and the smaller the weighted coefficient, the smaller the influence on the evaluation index.

ANOVA and LSD Methods
Microsoft Excel 2019 (Redmond, WA, USA) was used for preliminary organization of the data, and one-way analysis of variance (ANOVA) and the least significant difference (LSD) method were used to compare and test the significance of differences among the indicators using SPSS 26.0 (BM, Inc, Armonk, NY, USA) statistical software, and Sigma Plot 12.5 (Systat Software, San lose, CA, USA) was used to generate the graphs.

Effects of Cutting Times and Sowing Patterns on Stem and Leaf Characteristics and Production Characteristics
According to the variance analysis of stem and leaf traits and nutritional quality indices of herbage in the alfalfa and awnless brome seeding model, it was found that cutting stubble times and seeding mode had significant effects on leaf length, leaf width, leaf thickness, leaf area, stem length, stem diameter, leaf fresh weight, protein yield, dry matter yield, and fresh matter yield (p < 0.05).Leaf length, leaf width, leaf thickness, leaf area, stem length, stem diameter, protein yield, and dry matter yield were all affected by mowing and seeding pattern (Table 2).During the first crop sampling in the first year, the leaf length, leaf width, and leaf thickness of alfalfa under the alfalfa-awnless brome peer mixture were significantly higher than those under the hetero-row mixture and monoculture treatments, with the values of 17.89 mm, 9.28 mm, and 0.19 mm, respectively.Leaf area showed significant differences (p < 0.05) under all three sowing patterns, with the ranking of heterocomplex > peer mix > monoculture, thus indicating that the monoculture sowing pattern had the lowest leaf area (Table 3).The leaf length of the first sampled awnless brome (147.38 mm) was significantly lower under the heterogeneous row mix treatment than under the monoculture treatment.Similarly, the leaf width of the first sampled awnless brome (5.43 mm) was significantly lower under the peer mix treatment than the monoculture treatment.The leaf length, leaf width, leaf thickness, and leaf area of the second sampled alfalfa were significantly lower than the first sampled crop under the alfalfa-awnless peer mixture treatment.The leaf width (6.79 mm) and leaf area (0.61 cm 2 ) of alfalfa were significantly lower under the peer row mix treatment than under the heterogeneous mixture and the monoculture treatments (Table 3).The leaf length of awnless brome (16.65 mm) under the heterogeneous row mixture treatment was significantly higher than that of the peer mixture treatment.Moreover, the leaf thickness under the monoculture treatment was significantly lower than that of the other two treatments.There was no significant difference in the fresh weight of leaves among the three treatments.However, the leaf length of the second sampled awnless brome was significantly higher compared to the first sampled crop under the hetero-row mixture and monoculture treatments, with values of 186.62 mm and 169.88 mm, respectively.The leaf width of the second sampled awnless brome (5.21 mm) was significantly lower compared to the first sampled crop under the monoculture treatment, and the leaf area of the monoculture treatment (6.95 cm 2 ) was significantly lower than that of the hetero-row mixture treatment.
For the first crop sampling in the second year, the leaf length of alfalfa under the alfalfa-awnless brome hetero-row mixture was significantly higher than that under the peer mixture and monoculture treatments (Table 4), with a value of 20.15 mm.The leaf thickness under alfalfa-awnless brome peer and hetero-row mixture treatments was significantly lower than that under the other treatments, with values of 0.14 mm and 0.15 mm, respectively.Moreover, the stem length was significantly higher (78.68 cm) in the alfalfaawnless brome mixed treatment than in the monoculture treatment, while stem diameter was significantly higher (3.21 mm) in the alfalfa-awnless brome allopathic mixed treatment than in the monoculture treatment.
The leaf length of the second sampled alfalfa crop was significantly higher in alfalfaawnless brome peer (24.87 mm) and hetero-row mixtures (24.67 mm) than in the other treatments.The leaf width under alfalfa-awnless brome peer treatment (9.47 mm) was significantly higher than that under monoculture treatments.Similarly, the leaf area under alfalfa-awnless brome peer and hetero-row mixture treatments was significantly higher than that under the monoculture, with a value of 1.77 cm 2 , and there were no significant differences in stem length and diameter between the treatments (p > 0.05).

Effects of Different Sowing Patterns and Stubble on Stem Traits of Alfalfa and Awnless Brome
For the first crop sampling in the first year, alfalfa stem length was the highest under the alfalfa-awnless brome monoculture treatment at 36.32 cm.The stem length of awnless brome was significantly lower in the single sowing treatment (5.30 cm) than in the other treatments (p < 0.05); however, there were no significant differences in stem diameter and stem fresh weight among the three treatments.For the second crop sampled in the same year, alfalfa stem length and diameter were higher under the alfalfa-awnless brome peer mixture than under the other two treatments, and the highest alfalfa stem fresh weight was 15.52 g under the heterogeneous row mixture treatment.The stem diameter of awnless brome significantly differed from its stem length and stem fresh weight, with the most significant differences occurring under the heterogeneous row mixture.However, there was a non-significant difference in stem length and stem fresh weight among the three sowing pattern treatments (p > 0.05).
For the first crop sampled in the second year, the stem diameter (3.21 mm) and stem fresh weight (1.40 g) of alfalfa were higher under the alfalfa-awnless brome hetero-row mix treatment than under the other treatments.The stem length under the peer-row mix treatment was the longest (78.68 cm), and the stem length and diameter under the monoculture treatment significantly differed from those of the other treatments (p < 0.05).
The alfalfa stem diameter of the second crop sampled was significantly higher (with an increase of 0.49 mm) compared to the first crop sampled in the same year under the alfalfaawnless brome monoculture treatment.Similarly, the stem length of the second sampled alfalfa was significantly higher than that of the first sampled crop under the heterogeneous mixed and monoculture treatments by 77.21 cm and 70.29 cm, respectively.The stem fresh weight of the second sampled awnless brome was significantly lower under all treatments compared to the first sampled crop, considering that awnless brome growth was affected by climate.However, there were no significant differences in the stem length, stem diameter, and stem fresh weight among the three treatments (Table 5).

Effect of Different Sowing Patterns and Stubble on the Stem-to-Leaf Ratio of Alfalfa and Awnless Brome
The stem-to-leaf mass ratio indicates the distribution of plant material in the aboveground tissues [23,24], and a smaller stem-to-leaf ratio indicates a more favorable distribution of nutrients to the leaves, better palatability, and higher protein content [25,26].Furthermore, the average leaf aspect ratio and stem length/diameter ratio indicate the geometric fineness of leaves and stems, and larger values indicate that the leaf and stem growth tend to be more vertical, while smaller values indicate that leaves and stems tend to grow laterally.
In the first year, the alfalfa stem length/diameter ratio of the first crop to be sampled was significantly higher under the alfalfa-awnless brome heterogeneous mixture than under the peer mix treatment, indicating that its stem growth tended to be more vertical.The leaf aspect ratio (27.32), stem length-to-diameter ratio (2.93), and stem fresh weight ratio (0.62) of awnless brome under the first monoculture treatment were the smallest among all treatments, indicating that awnless brome had shorter and wider leaves, thicker and shorter stems, and relatively higher leaf biomass and smaller stem biomass under the monoculture treatment (Table 6).
The alfalfa leaf aspect ratio (2.24) of the second crop sampled under the alfalfa-awnless brome peer mixture treatment was the highest of all treatments, indicating a greater tendency for vertical leaf growth.There was a significant decrease in the leaf aspect ratio of the second crop under the alfalfa monoculture treatment and a significant increase in stem length-to-diameter ratio in the peer mix versus monoculture treatment.A significant decrease also existed in the stem-to-leaf fresh weight ratio in the monoculture treatment.The stem length-to-diameter ratio of the second crop of awnless brome under the monoculture treatment was the smallest among the treatments, and the leaf length-to-width ratio was significantly increased in the monoculture treatment.All three sowing patterns of awnless brome had significantly higher stem length-to-diameter ratios in the second crop, indicating that the stems tended to grow vertically (Table 6).

Effects of Different Sowing Patterns and Stubble on Yield Traits of Alfalfa and Awnless Brome
In the first year, the crude protein, dry matter, and fresh matter yield of the first sampled alfalfa were significantly higher in the alfalfa-awnless brome hetero-row mix, with the values of 0.81 t/ha, 4.51 t/ha, and 14.72 t/ha, respectively, than those of the other two sowing patterns.Similarly, the crude protein yield (1.19 t/ha) and dry matter yield (5.57 t/ha) of the second sampled alfalfa were significantly higher in the alfalfa-awnless brome hetero-row mix than in the single sowing.There was no significant difference in fresh matter yield between the three sowing patterns (Table 7).
In the second year, the crude protein (1.82 t/ha) and dry matter (10.25 t/ha) yields of the first sampled alfalfa crop in the second year were significantly higher in the alfalfaawnless brome hetero-row mix than in the other sowing methods.Similarly, the crude protein (2.99 t/ha) and dry matter (14.81 t/ha) yields of the second sampled alfalfa crop were significantly higher under alfalfa-awnless brome peer-row mix than under the other sowing methods, but there was no significant difference in fresh matter yield among the treatments (p > 0.05).The crude protein, dry matter, and fresh matter yields of the first sampled awnless brome under the alfalfa-awnless brome mix were the lowest among the three sowing patterns, with the values of 0.12 t/ha, 0.52 t/ha, and 1.06 t/ha, respectively.The second sampled awnless brome under the monoculture pattern had the best performance, with the crude protein, dry matter, and fresh matter yields of 1.2 t/ha, 4.92 t/ha, and 11.90 t/ha, respectively.Furthermore, there was no significant variation between stubbles for alfalfa and awnless brome in the same sowing patterns (Table 7).

Correlation Coefficients of Alfalfa and Awnless Brome Stem and Leaf Traits with Their Dry Matter Yields under Different Treatments
The equal-weight correlation degree was derived using forage dry matter yield as the reference column; however, the equal-weight correlation degree can only evaluate the significance of different varieties under the condition that the traits are equally important.The data of correlation coefficients among different cuts and treatments of alfalfa is in Table 8.The data of correlation coefficients among different cuts and treatments of awnless brome is in Table 9.
Practically, the importance of different traits in alfalfa is different, and their weights are expressed according to the magnitude of the correlation degree.The values obtained were ω1 = 0.1674, ω2 = 0.1609, ω3 = 0.1742, ω4 = 0.1582, ω5 = 0.1654, and ω6 = 0.1739.A higher value means that the trait contributes more to the dry matter yield.Therefore, the order of weight of each trait in the evaluation index was leaf area > stem diameter > leaf length > stem length > leaf width > leaf thickness (Figure 2).The weighted index of each trait in awnless brome was ω1 = 0.1703, ω2 = 0.1668, ω3 = 0.1760, ω4 = 0.1613, ω5 = 0.1672, and ω6 = 0.1584, with the order of evaluation index being leaf area > leaf length > stem length > leaf width > leaf thickness > stem diameter (Figure 3).Spatial distribution changes of plant populations (in the same row and alternate rows) can directly affect their access to light energy, interspecific competition above and below ground, and yield potential [27].Stem and leaf traits are intrinsic physiological and external morphological measures of adaptive responses developed by plants in response to environmental changes.The traits are closely related to plant biomass, plant access to and utilization efficiency of resources [28], and the survival strategies adopted by plants to obtain maximum material harvest [29].In this study, we found that alfalfa stem and leaf traits of the first sampled crop in the first year were significantly higher under the alfalfa-awnless brome mixed treatment than under the hetero-row mixed and monoculture treatments.However, the stem and leaf traits of the second sampled alfalfa crop were significantly lower than those of the first sampled crop in the monoculture and alfalfaawnless brome hetero-row mixed treatments.This phenomenon may be because awnless brome has poor regeneration and growth rate compared to alfalfa after mowing, and awnless brome is a low plant, while alfalfa is a high plant in the mixed population [30].After mowing, awnless brome was subjected to stronger shading by alfalfa under the same row mix treatment compared to the different row mix treatment, inhibiting the growth and development of awnless brome, eventually resulting in poor development of awnless brome stems and leaves.The leaf width and diameter of awnless brome in the first year were significantly lower under the alfalfa-awnless brome treatment compared to the other treatments.In addition, alfalfa stem and leaf traits of the second sampled crop in the second year were significantly improved under the alfalfa-awnless brome hetero-row mix treatment than under the peer mix and monoculture treatments.Moreover, the stem and leaf traits were significantly higher in the second sample alfalfa crop compared to the first sampled crop, indicating that the growth of awnless brome under the hetero-row mix treatment was promoted with increasing planting years and crop frequency, but suppressed under the other two sowing methods.This could be because mowing stimulates tillering of grasses and branching of leguminous plants, and the more adequate exposure of awnless brome to light after mowing under the inter-row mix treatment reduces the interspecific competition and promotes the compensatory growth of awnless brome [31,32].

Effects of Different Sowing Patterns on the Quantitative Traits of Forage Grasses
Forage yield is among the most important indicators for forage evaluation.Different plant species differ significantly due to their interspecific competitiveness, which affects the population size and the yield of mixed planted grass [33].DU Jun-ying [30] showed that alfalfa yield was significantly higher under alfalfa-awnless brome mixed sowing treatment than under the monoculture system.In this study, the protein, dry matter, and fresh matter yields of the first and second sampled alfalfa crops under the alfalfa-awnless brome peer mixed sowing pattern were not significantly different from those under single sowing.However, the yield of alfalfa crops under the alfalfa-awnless brome heterogeneous mix was significantly higher than that of the peer mix and monoculture.This was consistent with the results of Yuqiang Tian [34], which showed that the heterogeneous row mix treatment significantly increased the relative density of grass forage compared to the peer mix treatment, significantly increasing the yield.
The protein, dry matter, and fresh matter yields of the first and second sampled awnless brome crops were significantly higher under the monoculture treatment than under the peer and heterogeneous row mix treatments.There was no significant difference in all yields of awnless brome under the peer and heterogeneous row mix treatments.The difference between the awnless brome yield reported in our study and the results reported by Wang Bin [35] may be because awnless brome was mowed at the same time as alfalfa (at the early flowering stage of alfalfa) when it was in the gestation stage.The effect of light intensity and wind speed within the grass population during the mowing period ranked as monosowed awnless brome > peer mixed and inter-row mixed > monosowed alfalfa.For the temperature, the pattern was monosowed awnless brome > monosowed alfalfa > peer mix, inter-row mix.The microclimatic conditions of awnless brome under monoculture were better than those under the other treatments.Moreover, the mixed sowing changed the biomass accumulation pattern of awnless brome [36], eventually reducing awnless brome yield in the mixed sowing system.

Effects of Different Stem and Leaf Traits on Forage Yield Traits
Stem and leaf traits are important indicators of plant growth, and as the main photosynthesis organs, the leaves play a critical role in light energy acquisition, water and nutrient uptake, and organic matter synthesis [37,38].The stem length and diameter of forage grasses reflect the growth and development of forage grasses and grass productivity [39], indicating their important role as indicators of the agronomic traits of forage grasses.The gray correlation analysis and the dry matter yield of forage stem and leaf traits showed leaf area, stem diameter, and leaf length were the top three stem and leaf traits that contributed more to the dry matter yield of alfalfa [40].The slower-growing species in the mixed sowing system were affected by the shade effect because shade strongly affects the development of forages by changing their biomass allocation strategy and reducing the growth height of the forage [41][42][43].The stem length of awnless brome contributed more to its dry matter yield than stem diameter, suggesting that adequate light resources increase the dry matter yield of awnless brome.

Conclusions
Mixing alfalfa with awnless brome in different rows can improve the stem and leaf functional traits, dry matter yield, and crude protein yield of both forages.The inter-row mixing significantly improved the stem and leaf functional traits of alfalfa but limited the growth and development of awnless brome in the late growth stage.The weighting order of the evaluation indices was dry matter yield > protein yield > leaf area > leaf length > leaf width > stem diameter > leaf thickness > stem length for alfalfa and dry matter yield > protein yield > leaf area > leaf length > leaf width > stem length > leaf thickness > stem diameter for awnless brome.In summary, heterogeneous row mixing is more suitable than the other treatments considered for alfalfa and awnless brome mixed planting.

Figure 1 .
Figure 1.Annual average temperature of test field in Jiaozhou area.Data from China Meteorolo Data Service Centre (http://data.cma.cn 8 February 2020).

Figure 1 .
Figure 1.Annual average temperature of test field in Jiaozhou area.Data from China Meteorological Data Service Centre (http://data.cma.cn 8 February 2020).

Figure 2 .
Figure 2. Correlation of Medicago sativa stem and leaf indicators under different sowing m LL is leaf length.LW is leaf width.LT is leaf thickness.LA is leaf area.SL is stem length.SD diameter.Peer mix (P), heterocomplex (H), monoculture (M).Significance levels are as foll p < 0.001, ** p < 0.01, * p < 0.05.

Figure 2 .
Figure 2. Correlation of Medicago sativa stem and leaf indicators under different sowing m LL is leaf length.LW is leaf width.LT is leaf thickness.LA is leaf area.SL is stem length.SD diameter.Peer mix (P), heterocomplex (H), monoculture (M).Significance levels are as follo p < 0.001, ** p < 0.01, * p < 0.05.

Figure 3 .
Figure 3. Correlation of Bromus inermis stem and leaf indicators under different sowing m LL is leaf length.LW is leaf width.LT is leaf thickness.LA is leaf area.SL is stem length.SD diameter.Peer mix (P), heterocomplex (H), monoculture (M).Significance levels are as follo p < 0.001, ** p < 0.01, * p < 0.05.

Author Contributions:
Conceptualization, X.Y. and F.M.; methodology, F.M.; software, X.T.; validation, X.L.; formal analysis, X.Y.; investigation, W.T.; resources, C.Y.; data curation, G.Y.; writingoriginal draft preparation, Y.Z.; writing-review and editing, J.S.; visualization, X.Y.; supervision, Y.X.; project administration, J.S.; funding acquisition, J.S.All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by the China Agriculture Research System (CARS-34), the Natural Science Foundation of Shandong Province (ZR2020QC189), the National Key Research and Development Program of China (No. 2022YFD1300802), the Doctoral Scientific Research Startup of Qingdao Agricultural University (6631116024), the Shandong Forage Research System (SDAIT-23-01), and the First-Class Grassland Science Discipline Program of Shandong Province (1619002), the "Youth Innovation Team Plan" of Universities in Shandong Province (2022KJ166).Data Availability Statement: Not applicable.

Table 1 .
The physical and chemical properties of soil at 0~30 cm. )

Table 1 .
The physical and chemical properties of soil at 0~30 cm.

Table 2 .
Effects of cutting times and sowing patterns on stem and leaf characteristics and production characteristics.
3.2.Effects of Different Sowing Patterns and Stubble on Stem and Leaf Traits of Alfalfa and Awnless Brome 3.2.1.Effect of Different Sowing Patterns and Stubble on Leaf Traits of Alfalfa and Awnless Brome

Table 3 .
The dynamics of the leaf traits of alfalfa.

Table 4 .
The second-year Alfalfa stem and leaf trait dynamics under different seeding methods.
Peer mix (P), heterocomplex (H), monoculture (M).Capital letters indicate significant differences in the same planting pattern at different stubble times at the 0.05 level (p < 0.05).Lowercase letters indicate significant differences in the different planting patterns at the same stubble times at the 0.05 level (p < 0.05).

Table 5 .
The dynamics of the stem traits of alfalfa and awnless brome.

Table 6 .
The dynamics of the stem-leaf ratio of alfalfa and awnless brome.
Peer mix (P), heterocomplex (H), monoculture (M).Capital letters indicate significant differences in the same planting pattern at different stubble times at the 0.05 level (p < 0.05).Lowercase letters indicate significant differences in the different planting patterns at the same stubble times at the 0.05 level (p < 0.05).

Table 7 .
The dynamics of the yield of alfalfa and awnless brome.
Peer mix (P), heterocomplex (H), monoculture (M).Capital letters indicate significant differences in the same planting pattern at different stubble times at the 0.05 level (p < 0.05).Lowercase letters indicate significant differences in the different planting patterns at the same stubble times at the 0.05 level (p < 0.05).

Table 8 .
Correlation coefficients among different cuts and treatments of alfalfa.

Table 9 .
Correlation coefficients among different cuts and treatments of awnless brome.