Abstract

Population structure determines sugarcane yield, of which canopy structure is a key component. To fully understand the relations between sugarcane yield and parameters of the canopy structure, 17 sugarcane varieties were investigated at five growth stages. The results indicated that there were significant differences between characterized parameters among sugarcane populations at different growth stages. During sugarcane growth after planting, leaf area index (LAI) and leaf distribution (LD) increased, while transmission coefficient for diffuse radiation (TD), mean foliage inclination angle (MFIA), transmission coefficient for solar beam radiation penetration (TR), and extinction coefficient ( ) decreased. Significant negative correlations were found between sugarcane yield and MFIA, TD, TR, and at the early elongation stage, while a significant positive correlation between sugarcane yield and LD was found at the same stage. A regression for sugarcane yield, with relative error of yield fitting less than 10%, was successfully established: sugarcane yield = 2380.12 + 46.25 × LD − 491.82 × LAI + 1.36 × MFIA + 614.91 × TD − 1908.05 × TR − 182.53 ×   + 1281.75 × LD − 1.35 × MFIA + 831.2 × TR − 407.8 ×   + 8.21 × MFIA − 834.50 × TD − 1695.49 ×   .

1. Introduction

Sugarcane (Saccharum spp. hybrids), with economical importance, has by far the highest yield in crops. Since vegetative stalks are targeted in sugarcane harvest, population structure forms the basis of sugarcane yield. The population structure refers to the dynamics of the distribution and arrangement of each single plant, total leaf area, total plant number, and total root weight in time and space. The crop canopy structure plays a key role in the population structure because it directly affects not only the interception of sunlight but also the photosynthetic efficiency and crop yield of the population through the influence on microenvironment of water, heat, and atmosphere on the canopy [13]. Therefore, the light use efficiency of a crop population is directly related to its canopy structure. The biological yield of a crop and its organ distributions are the ultimate results of the canopy structure, which will be reflected in photosynthetic characteristics. There is vital significance for the exploration of appropriate population structure and the cultivation of high-yield and high-quality crops through the study on canopy structure [4]. The leaf morphologies and population structures of other crops such as rice, wheat, and cotton have already been investigated [511].

The cultivation of an optimized population structure to improve optical radiation distribution within the canopy and to increase energy utilization is the basis for high yield in crops. The main indices for light radiation within the crop canopy are leaf area index (LAI), mean foliage inclination angle (MFIA), transmission coefficient for diffuse radiation (TD), transmission coefficient for solar beam radiation penetration (TR) in different altitude angles and azimuth angles, extinction coefficient ( ), and leaf distribution (LD). It is reported there is a close relationship among these indices [511]. Canopy light interception can be influenced by the LAI, and it rises along with the increasing LAI. Light interception peaks when LAI is at its optimum. The photosynthetic rate can be affected by the light interception and an appropriate increase in light interception rate of a population can improve the photosynthetic capacity and thus increase production [610]. Canopy value at all levels reflects vertical distribution of the leaf area and leaf angle and also the vertical diminishing status of canopy light. In recent years, the investigation on the canopy structure of crops has been greatly facilitated due to the emergence of canopy image analysis techniques. Moreover, the population structure characteristic values such as LAI, MFIA, TD, TR in different altitude angles and azimuth angles, , and LD were successfully determined [9].

Sugarcane is not only a main sugar crop but also an important energy one. The selection of yield and quality-related traits is a key for breeding superior sugarcane varieties. The relationships between leaf morphologies and yields or sucrose contents in sugarcane have been explored [1215]. However, ecophysiological studies on relationships between characteristics of the canopy structure and sugarcane growth as well as yields have been rarely reported [1214]. Currently, there is no report on the effects of sugarcane leaf morphology and its spatial distribution on radiation characteristics of canopy.

From all the above, the prediction of sugarcane yield-related traits in population structure at different growth stages is significant in theory and practice for breeding high-yield sugarcane varieties. The canopy spaces and light distributions of different sugarcane populations were investigated in the present study. The relationships between structural characteristics and their yield-related traits were also revealed among different varieties. This study aims to provide references for breeding of high-yield sugarcane varieties on the basis of high-quality population structure.

2. Materials and Methods

2.1. Materials

The tested sugarcane varieties were FN94-0403, FN94-0744, FN95-1726, FN96-0907, FN98-10100, ROC10, GT94-116, GT95-118, GT96-211, GT96-44, GT97-18, MT70-611, YT92-1287, YT96-107, YT9-6794, YT96-835, and YT96-86. These varieties were planted in the farm of Sugarcane Research Institute, Fujian Agriculture and Forestry University (longitude: 119.23 E, latitude: 26.08 N), as a randomized complete block design. There were three rows for each variety, with the amount of 45,000 sugarcane two-bud sets per hm−2. The plot area was 31.35 m2, with the row length of 9.5 meter and the row space of 1.1 meter. Field management level was slightly higher than that of the local production fields, with timely cultivator earth, fertilization, irrigation, and pest control. During field management, the same pilot and the same technical measure should be completed on the same day.

2.2. Methods

The canopy structure parameters of sugarcane were determined by CI-100 digital plant canopy imager (Washington, USA) at the stages of seedling (May), tillering (June), early elongation (July), rapid elongation (August), and late elongation (September), respectively. The measured indices include LAI, MFIA, TD, TR, , and LD. Fifteen plants of each variety were randomly selected for measurement on canopy structure parameters with three replicates. In mid-November, plant height, stalk diameter, and stalk number for all tested varieties were measured. Single stem weight and cane yield were calculated according to the following formulas: single stem weight = 0.785 × plant height × stalk diameter2, cane yield = single stem weight × effective stalk number per unit area.

The measurement was conducted according to established methods [13]. During evening time when sunshine was not particularly strong and the sky remained cloudless, observation stick with a fish-eye probe was installed centrally between the rows and was adjusted to exclude any influence from the outside including shadows. Five images were taken in each district. Software for plant canopy analysis, provided by U.S. CID, Inc. (Washington, USA), was applied to calculate the canopy structure parameters. During measurement, LD was represented by distribution frequency of the leaf within each azimuth. The leaf blade azimuth is the angle between the normal direction of leaf surface projection and the direction North; this was then divided clockwise into four directions 90° between each other. and TR were derived according to zenith angle, which was the angle between the target direction and the zenith direction. In this study, CI100 fish-eye lens were set as the observation point, and four angles, 27°, 45°, 63°, and 81°, were determined from the beginning of 9° to 18° zenith angle.

2.3. Statistical Analyses

The variance analysis, factor analysis, and regression analysis were performed on the means of data by DPS statistical software [16].

3. Result and Analysis

3.1. Seasonal Changes of Population Structure Characteristics of Sugarcane Varieties
3.1.1. Leaf Area Index (LAI)

Sugarcane LAIs varied significantly at different growth stages. As seasons change, the LAIs of all varieties increased (Table 1). There was no significant difference in LAI between sugarcane varieties at seedling (May) and tillering stages (June), but at the early elongation stage (July) there was a significant difference. The average LAIs of FN94-0403, FN98-10100, GT96-44, GT94-116, FN96-0907 and MT70-611 at the early elongation stage (July) were larger than those of YT96-794, YT96-835, YT96-107, and GT95-118. There was a relatively small difference in LAI between sugarcane varieties at the rapid elongation stage (August) and the late elongation stage (September). At the rapid elongation stage, the largest average LAI was MT70-611, while the smallest average LAIs were YT96-794 and ROC10. At the late elongation stage (September), the largest average LAIs were GT95-118 and FN98-10100, while the smallest average LAI was YT96-794.

3.1.2. Mean Foliage Inclination Angle (MFIA)

There were significant differences in the sugarcane MFIAs at different growth stages (Table 2). As seasons changed, the maximum MFIA appeared at the seedling stage (May), while the minimum one appeared at the rapid elongation stage (August). There was no significant difference in MFIA between sugarcane varieties at seedling stage (May), tillering stage (June), and the rapid elongation stage (August), but there was a significant difference at the early elongation stage (July). Compared to those of GT95-118, ROC10, and FN95-1726 at the early elongation stage (July), the MEIAs of FN94-0403 and MT70-611 were smaller. There were large differences in the average MFIA among sugarcane varieties at the late elongation stage (September). The MFIAs of three sugarcane varieties, FN96-0907, GT97-18, and FN94-0403, were larger than those of another two sugarcane varieties, GT96-44 and YT96-107.

3.1.3. Transmission Coefficient for Diffuse Radiation (TD)

There were significant differences in the sugarcane TD at different growth stages (Table 3). As seasons changed, the maximum TD appeared at the seedling stage (May), while the minimum TD was found at the late elongation stage (September). There was no significant difference in TD between sugarcane varieties at the seedling stage, the tillering stage (June), the rapid elongation stage (August), and the late elongation stage (September), but there was a significant difference at the early elongation stage (July). The TD of GT95-118 at the early elongation stage (July) was larger than other sugarcane varieties including FN98-10100.

3.1.4. Transmission Coefficient for Solar Beam Radiation Penetration (TR)

The incident light, reflecting solar radiation, can be considered as parallel light beam reaching ground, which is the major light component of plant photosynthesis. The distribution of incident light in a population is determined by the population structure. It is generally considered that light decreases from top to bottom layers of the population, and even the lowest layer of the population would receive this light. The TR reflects the distribution of incident light in the crop population, relating to the transmittance of crop population and the interception of light. Effects were extremely significant on TR from different sugarcane varieties, growth stages, zenith angles, variety × zenith angle, and growth stage × zenith angle. TR decreased linearly with increasing the zenith angle from top to bottom in the vertical direction (Table 4). When the zenith angles were 9° and 27°, respectively, there was little difference in TR. However, when the zenith angle reached 63°, TR decreased sharply. When the zenith angle increased from 9° to 81°, TR decreased from 0.23 to 0.04. The results indicated that the distribution and transmission of light in sugarcane population decreased gradually from top to bottom, while the interception of light increases from top to bottom. There were significant differences in TRs of sugarcane varieties at different growth stages. With a change of seasons, the maximum TR appeared at the seedling stage (May), while the minimum TR appeared at the late elongation stage (September). There were significant differences in the TR between sugarcane varieties at each growth stage. Among all tested sugarcane varieties, GT97-18, FN95-1726, and YT92-1287 had the largest TRs at the seedling stage (May), and GT96-211 had the smallest. Similarly, FN98-10100, MT70-611, and YT96-794 had the largest TRs and YT92-1287 had the smallest at the tillering stage (June). GT95-118 and ROC10 had the largest TRs, while FN94-0403 and FN98-10100 had the smallest at the early elongation stage (July). ROC10 and YT96-794 had the largest TRs, while FN94-0744 and MT70-611 had the smallest at the rapid elongation stage (August). FN94-0403 and YT96-794 had larger TRs than those of FN98-10100 and GT95-118 at the late elongation stage (September).

3.1.5. Extinction Coefficient ( )

The extinction coefficient ( ) represents the attenuation characteristic parameters of the light intensity in the crop population in the vertical direction [2]. to represent different values of canopy at different altitude angles (0~90°). As indicated in Table 5, increased with an increasing altitude angle from top to bottom in the vertical direction. When the altitude angles were 9° and 27°, respectively, the difference of was small. When the altitude angle reached 45°, increased sharply. It should be noted that there were significant differences in values at different growth stages. As seasons changed, the maximum appeared at the seedling stage (May) and the minimum at the late elongation stage (September). However, there was no significant difference in the between sugarcane varieties at different growth stages.

3.1.6. Leaf Distribution (LD)

Leaf distribution (LD) refers to the scattering density in leaves at each direction. This index characterizes the spatial distribution characteristics of the canopy leaves of the crop population. There were significant differences in LD at different azimuth angles (Table 6). At the seedling stage (May), the tillering stage (June), the early elongation stage (July), and the rapid elongation stage (August), the LDs at 270~360° azimuth angle were the highest, followed by those at 180~270°, while those at 90~180° were the lowest. However, at the late elongation stage (September), the LDs at 0~90° azimuth angle were the highest, while those at 180~270° were the lowest. There were significant differences in LDs at different growth stages. As seasons changed, the highest LDs appeared at the seedling stage (May) and the lowest at the late elongation stage (September). Within some stages like the seedling (May) and the early elongation (July), there were significant differences between some sugarcane varieties. For example, FN98-10100 and GT96-211 had significantly higher LDs than that of GT97-18 at the seedling stage (May), and FN98-10100 had a higher LD than that of GT95-118 at the early elongation stage (July). There was no significant difference in LD between sugarcane varieties at the tillering stage (June), the rapid elongation stage (August), and the late elongation stage (September).

3.2. Factor and Regression Analyses of Canopy Structure Parameters and Their Correlation with Yield-Related Traits in Sugarcane

From thirty canopy parameters of sugarcane at five stages, the first six principal factors related to four yield-related traits were characterized, the cumulative contribution rates were generated, and the initial factor loading matrix is presented in Table 7. The cumulative contribution rate of the characteristic values of the first six principal factors reached 90%. The main information contained in the first six principal factors represent most of the information from the thirty canopy parameters and the initial characteristic parameters of the four yield-related traits. The selection of these 6 factors not only preserved the major part of the information but also achieved the goal of dimensionality reduction. The communalities of the traits were listed in Tables 8 and 9. The communality of the th trait is the contribution of the total principal factors to the total variance of the th trait. The larger the communality is, the better the representativeness of this principal factor to the variable is.

The communalities of 34 traits were large, indicating that the six principal factors were well representative of these 34 traits after performing the varimax rotation, which maximized the total variance of each factor in the loading matrix. The factor loading matrix by varimax rotation of the 6 principal factors was shown in Table 8. The factor loading matrix after varimax rotation had significantly increased loading values of the important variables represented by the principal factors (refer to the initial loading values in Table 7). The principal factors become more biologically meaningful after varimax rotation. As shown in Table 8, the traits with the major roles in the first principal factor were the MFIA and at the early elongation stage; LAI, MFIA, TD, TR, and LD at the rapid elongation stage; and TR at the late elongation stage. Their loading values were large. The loading values of LAI, M FIA, TD, TR, , and LD at the seedling stage were larger than others in the second principal factor. The loading values of LAI, M FIA, TD, TR, , and LD and stalk number per hectare at the early elongation stage were larger than others in the third principal factor. The loading values of LAI, MFIA, TD, TR, , and LD at the tillering stage were larger than others in the fourth principal factor. The loading values of at the rapid elongation stage, LAI and TD at the late elongation stage, and the plant height were larger in the fifth principal factor. The loading values of MFIA and at the late elongation stage, stem diameter, single-stalk weight, and sugarcane yield were larger in the sixth principal factor. There were significant negative correlations between the plant height and the MFIA, between the plant height and the at the rapid elongation stage, and between the plant height and the TR at the late elongation stage, while there were significant positive correlations between the MFIA and TD at the late elongation stage. There were significant negative correlations between the stalk number per hectare and TD or TR at the early elongation stage, while there were significant positive correlations between stalk number per hectare and the LAI or LD at the early elongation stage. There was a significant negative correlation between the stem diameter and LD at tillering stage. There were no significant correlations between stem diameter and the MFIA or at the late elongation stage. There was a significant negative correlation between the single-stalk weight and LD at tillering stage, while there were significant positive correlations between single-stalk weight and the MFIA or at the late elongation stage. There were significant negative correlations between the sugarcane yield and the MFIA, TD, TR, or at early elongation stage, while there was a significant positive correlation of the sugarcane yield with the LD of the early elongation stage. There were correlations of sugarcane yield with the MFIA and at the late elongation stage, but they were not significant. A regression analysis was also conducted for the canopy structure parameters and their correlation with yield-related traits in sugarcane. The sugarcane yield was predicted by the canopy structure characteristic parameters. The regression equation was as follows: sugarcane yield = 2380.12 + 46.25 × LD (tillering stage) − 491.82 × LAI (early elongation stage) + 1.36 × MFIA (early elongation stage) + 614.91 × TD (early elongation stage) − 1908.05 × TR (early elongation stage) − 182.53 ×   (early elongation stage) + 1281.75 × LD (early elongation stage) − 1.35 × MFIA (rapid elongation stage) + 831.2 × TR (rapid elongation stage) − 407.8 ×   (rapid elongation stage) + 8.21 × MFIA (late elongation stage) − 834.50 × TD (late elongation stage) − 1695.49 ×   (late elongation stage) ( ). The relative error of yield fitting for all these 17 test sugarcane varieties was less than 10%.

4. Discussion and Conclusions

Crop yield mainly depends on the light receiving ability and light distribution in the crop population. The leaves are the main organ for photosynthesis. Therefore, the characteristics of leaves are important in light distribution and photosynthetic traits of a crop’s population [14]. The shape, size, quantity, and spatial distribution of leaves determine light use efficiency and influence light environment and vice versa. Sugarcane varieties with short, narrow, thick, and straight canopy have high-yield potentials [12]. Leaf morphology and canopy structure directly affect light interception of canopy and light distribution at different layers [610]. Leaf morphology and spatial distribution influence radiation characteristics of the canopy [1719]. In the present study, large radiation parameters in canopy, especially in lower layers of the canopy, were found in sugarcane genotypes with short, narrow, and straight leaves, which is beneficial to light absorption, resulting in favorable conditions for sugarcane yield and sugar accumulation during the later growth period.

Appropriate population structure is essential to improve energy use efficiency and thus the high crop yield. A positive correlation was found between proper LAI and dry matter accumulation [20]. Enhancement of population LAI, especially at the later development stages, can increase dry matter accumulation. MFIA influences the plant ability to receive solar radiation and the canopy distribution, especially on values. By breeding varieties with proper leaf angles and distribution, an ideal population structure can be obtained with the most effective leaf area. Zhang et al. [21] revealed that the canopy MFIA of winter wheat decreases with crop growth along with increasing nitrogen application. Even with the same LAI on two populations, different light distributions, canopy photosynthetic rates, and dry matter accumulations can be observed if spatial distribution of leaves and average leaf angles are different [20]. Numeral researches have proved that crop productions tightly correlate with the functions of canopy populations [2022]. Maddonni et al. [22] pointed out that light attenuation of maize, a row crop, was influenced by canopy architecture in terms of size, shape, and orientation of shoot components. Monsi and Saeki [23] showed that, under full light conditions, canopy photosynthesis was maximized at a high LAI with vertically inclined leaves, while, under low light conditions, it was at a low LAI with horizontal leaves. They suggested that plants have evolved standing stem with different orientations of branches and leaves to maximize canopy photosynthesis. Canopy photosynthesis models are useful to analyze the size structure of populations in plant communities and to predict the structure and function of future terrestrial ecosystems [24]. Stewart et al. [25] found that plants with very upright leaves can have both the smallest and largest daily canopy photosynthesis, depending on row widths. Reta-Sánchez and Fowler [26] revealed that early canopy modifications simulating plant characteristics such as reduced plant height, short branches, and modified leaf shape increased light availability at the medium and upper part of the canopy. Liu et al. [27] demonstrated that narrow-row planting patterns improved the canopy structure, allowed more IPAR (light intercepted by a canopy) to reach the middle-row strata of the canopy, and enhanced the leaf photosynthetic characteristics of maize crops at silking stage compared with the control. As illustrated in the present study, sugarcane canopy structure can be affected by planting pattern and other physiological characters such as light interception and radiation use efficiency. Besides, a large crop photosynthetic area at the late growth stage is one of the most important features of population structure on productive sugarcane farmland, which are similar to the previous reports [20–27].

Canopy structure, reflected by vertical and spatial distribution, orientation, and density of foliage and its supporting structures, determines the pattern of light attenuation, the distribution of photosynthesis, respiration, transpiration, and nutrient cycling in the canopy and therefore the yield and quantity of crops [28]. The present study indicated that different influences of canopy structure on sugarcane yield component factors at different stages. The mutual influences on each factor were reflected by the correlations of the parameters in the canopy structure. Understanding these influences helps the breeders to select varieties in optimization of the canopy structure at each developmental stage and thus helps the sugarcane farmers in their cultivation and yield prediction. The correlation and regression models established in this study demonstrated that the canopy structure at the early elongation stage is closely related with the formation of sugarcane yield because the sugarcane varieties with small MFIA, TD, TR, , and large LD at the early elongation stage have the greatest number of leaves in unit space and the spatial distributions of leaves are uniform. Ecologically, this kind of parameter combination is in favor of the interception of light by the sugarcane population and the improvement of light energy utilization. Physiologically, more photosynthetic products could be accumulated in this condition, which would generate high-yield sugarcane population based on a solid foundation of a high-quality canopy structure for sugarcane growth and development [17, 28]. What should also be stressed is that all the 17 sugarcane varieties are popular currently in China and they do cover a relatively wide range in population structure for making regression, which should help to distinguish between different types of sugarcane varieties and canopy structure characteristics and thereby reduce the test result errors and ensure the reliability of the regression for sugarcane yield established in this study.

Given that the organic substances required for sugarcane growth are produced by sugarcane leaves through photosynthesis, the canopy structure and light distribution of sugarcane at the elongation stage thus play a crucial role in the structure formation, the accumulation and partitioning of photosynthetic products, and thus the yield formation. Our study clearly demonstrates this phenomenon by the correlations between the main canopy structure parameters at the early elongation stage and sugarcane yield. As seasons changed, the LAI and LD increase, while the TD, MFIA, TR, and decrease, which is in accordance with previous reports [1719, 28]. Attention to some traits, which are negatively correlated with yields at specific growth stages, is required in breeding and cultivation. When selecting the yield-related traits based on the canopy structure parameters of sugarcane, the values for the above mentioned traits are still required to be optimized. For example, the mutual shading between leaves should be avoided. In short, it is significant in breeding and cultivation of high-yield and high-quality sugarcane varieties through the rapid and accurate confirmation of characteristics of canopy structure and timely optimization of the population size.

Conflict of Interests

The authors declare no conflict of interests.

Acknowledgments

This work was supported by the earmarked fund for the Modern Agriculture Technology of China (CARS-20) and Research Funds for Distinguished Young Scientists in Fujian Agriculture and Forestry University (xjq201202). The authors especially thank Luguang Wu in School of Agriculture and Food Sciences, University of Queensland, Brisbane, Australia, for his critical revision and valuable comments on this paper.