Next Article in Journal
Enhancing Drought Tolerance and Striga hermonthica Resistance in Maize Using Newly Derived Inbred Lines from the Wild Maize Relative, Zea diploperennis
Previous Article in Journal
Identifying Cotton Fields from Remote Sensing Images Using Multiple Deep Learning Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Competition for Light Interception in Cotton Populations of Different Densities

1
Zhengzhou Research Base, State Key Laboratory of Cotton Biology, Zhengzhou University, Zhengzhou 450001, China
2
State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Anyang 455000, China
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(1), 176; https://doi.org/10.3390/agronomy11010176
Submission received: 14 December 2020 / Revised: 13 January 2021 / Accepted: 15 January 2021 / Published: 18 January 2021

Abstract

:
Modification of the cotton canopy results in shade avoidance and competition for light, which shows that density and spatial arrangement of cotton have a great impact on light interception. This experiment was conducted in 2018 and 2019 in the experimental field at the Institute of Cotton Research of Chinese Academy of Agricultural Science in Anyang city, Henan Province, China. Six plant densities of cotton variety SCRC28 were used to assess spatial competition for light in cotton populations during the whole growing period. Light interception data were collected and analyzed according to the spatial grid method and the extension of Simpson’s 3/8 rule. The results showed that at the bottom of the canopy, greater light interception was observed at high densities than at low densities, while in the external part of the layer of the canopy in the horizontal direction, low light interception was recorded at low densities. Leaf area, aboveground biomass and plant height were obviously correlated with light interception, and the cotton population with a higher density (8.7 plants m−2) performed best at the light interception competition, and with the highest yield. The results will provide guidance on light management through the optimization of the structure of the canopy to provide more solar radiation and a significant basis by which to improve the management of light and canopy architecture.

1. Introduction

Cotton (Gossypium hirsutum L.) is an important cash crop in China [1] and plays a role in national defense, medicine, and industry [2]. However, as China’s total arable land area decreases each year [3], the competition between grain and cotton is increasing, and cotton production decreases each year. Keeping the total area of cotton cultivation unchanged or slightly reduced while achieving an increase in total cotton production and cotton yield requires understanding of light competition among plants. Traditional cotton planting technology is inefficient because labor costs are high and operation processes are cumbersome, and rapid development of the cotton industry is not applicable [4], hindering cotton production [5]. To solve these problems, a highly productive agricultural production system is needed that combines maximum resource utilization with automation provided by agricultural machinery [6]. Cotton is a crop with indeterminate growth habit, and lighting is a key factor that can directly and indirectly determine its quality and yield [7,8].
Photosynthetically active radiation (PAR) at wavelengths of 400–700 nm is used in photosynthesis by green plants [9,10,11]. It is important to identify the characteristics of light interception that improve the potential photosynthetic activity of green plants [12,13], as high efficiency in the absorption and conversion of PAR can affect crop biomass and improve crop economic yield [14].
Cotton is a species that is suitable for growth in high-light conditions [7]. The duration of sunlight has a direct impact on growth and development, and lack of light directly affects the structure of cotton stems and leaves and the flowering time. Studies have shown that cotton yields are affected by the total amount of solar radiation that is intercepted by cotton and radiation use efficiency (RUE) [15]. The amount of light intercepted by cotton is influenced by factors such as cotton canopy structure, while RUE is influenced by factors such as photosynthetic characteristics of different varieties and the environmental conditions during cultivation [16].
Density can be increased by using narrow rows [17]. Densities and spatial arrangement of cotton plants have a great impact on light interception because modification of the cotton canopy results in shade avoidance and competition for light. A number of studies have shown that the leaf area index (LAI) and canopy structure affect the accumulation of intercepted PAR [18,19,20,21]. The architecture of the cotton canopy has significant influence on the photosynthetic process and the ability of light to enter the canopy [22]. In cotton plants, the canopy is directly linked to leaf shape and structure, and directly affects light interception and overall yield, while the LAI determines the interception of solar radiation and light capture efficiency [23]. Previous studies have shown that the relationships between canopy structure and light interception can be measured by using different methods. Zhi’s study showed that the cotton canopy structure affected the spatial distribution of PAR, which was highly correlated with LAI and biomass [24].
A variety of plant models have been developed and used to describe the relationships between different plant growth. The 3-D virtual plant model, also called the functional-structural plant model, has been used to describe the physiological processes involved in plant growth, such as the coordination of organ expansion dynamics and geometric variables [25]. Herve Rey [26] presented a method to build a 3-D virtual flower model by using AMAP and MMR software and simulated light interception during different periods of growth and for different organs of the plant. Gu’s study [27] involved designing a functional-structural plant model named CottonXL to simulate the growth of cotton leaves and fruits and showed good accuracy in simulating the LAI. This model was used by Mao to estimate light interception by cotton in two growing systems and to determine row distance and plant densities [28].
This present study was carried out to identify the spatial distribution of light interception in cotton populations of different densities throughout the whole growing period, with a particular aim to evaluate (i) the trends of and differences in canopy light interception by these cotton populations during the process of cotton growth and (ii) the changes in LAI and aboveground biomass with canopy light interception among these cotton populations. The results provide guidance on light management of cotton through the optimization of canopy structure to increase plant acquisition of solar radiation.

2. Materials and Methods

2.1. Experimental Site

The experiment was conducted in 2018 and 2019 in an experimental field at the Institute of Cotton Research of the Chinese Academy of Agricultural Science in Anyang city, Henan Province, China (36°06′ N, 114°21′ E). The experimental field is a long-term cotton cropping field with medium loam soil and total N, P, and K contents in the 0–40 cm layer of 0.65, 0.01, and 0.11 g kg−1, respectively. Monthly average temperature data for the two years and the 5-year average (2013–2017) during the growing season are presented in Figure 1. The average air temperature during the growing season was 23 and 22 °C in 2018 and 2019, respectively, while the 5-year average was 22 °C. Annual rainfall in 2018 and 2019 was 402 and 419 mm, respectively, while the 5-year average was 422 mm. The monthly average PAR and total sunshine hour data for the two years during the cotton growing season are presented in Table 1. The annual sunshine hours were 2165 and 2143 h in 2018 and 2019, respectively, while the 5-year average was 1927 h. The average available PAR during the seedling and reproductive stages was higher than that during the other stages, while the average air temperature was lower. The whole cotton growing season was slightly cooler in 2019 than in 2018, and more rainfall occurred in 2019 than in 2018.

2.2. Experimental Design

Three replicates each of six plant densities (D1: 1.5, D2: 3.3, D3: 5.1, D4: 6.9, D5: 8.7, D6: 10.5 plants m−2) of cotton variety SCRC28 were arranged in a randomized complete block design (RCBD). Each plot had an area of 64 m2: 8 m in both length and width. Within each plot, the cotton was grown in 10 rows with row spacing of 0.8 m. The growth periods of the cotton are shown in Table 2. Sowing was performed on 20 April 2018, and 18 April 2019. Before sowing, the field was irrigated with an approximate volume of water of 290 m3 ha−1, and basal fertilizer was applied to the field (225 kg ha−1 N, 150 kg ha−1 P2O5, and 225 kg ha−1 K2O). Medium-term management involved mechanical cultivation methods for weeding and spraying pesticides and plant growth regulators, and chemical control methods were used to control pests and diseases.

2.3. Data Collection

Data on light interception, leaf area, biomass accumulation and plant height and yield were collected and recorded at specific stages of cotton growth during 2018 and 2019 on different days after emergence (DAE).

2.3.1. Light Interception

Light interception data were collected during the growing season at an interval of 15 days at 30, 45, 60, 75, 90, 105 and 120 DAE. Data were collected on a sunny day with no clouds, and the measurements began at 10 a.m. and finished half an hour later. Two specific rows of cotton were selected to measure light interception in each plot. In this study, the spatial grid method was used to identify stable measurement points in two cotton rows. In terms of horizontal distance, 0, 20, 40, 60 and 80 cm were selected as the measurement points. For vertical distance, the distance from the ground to the top of the canopy was divided into layers considering an interval of 20 cm. In this way, a spatial grid perpendicular to the cotton plants was obtained, and the measurement points in the two cotton rows are shown in Figure 2.
For each measurement point, a 1 m-long portable linear light quantum sensor (LI-191SA, LI-COR, Lincoln, NE, USA) was used to measure the intercepted PAR intensity and the reflected PAR intensity, and a LI-1400 instrument (LI-COR, Lincoln, NE, USA) was used to collect and record the data. When the measurements were performed at the various points, the sensor was placed parallel to the row direction and faceup to measure the intercepted PAR intensity and then reversed and placed facedown to measure the reflected PAR intensity. At the same time, another sensor was placed 20 cm above the cotton canopy to record the total intensity of the real-time PAR. Each sensor automatically recorded the total intensity of real-time incident PAR every 5 s.
A coordinate system was created in the vertical direction for the testing area, and the locations of the points measured and to be measured in a particular section of the testing area were used to create a Grid(i,j) file. Each element of the file recorded the position of the point measured and to be measured in the vertical direction. i and j represented the abscissas and ordinates of the section of the testing area, respectively. To evaluate the differences in light interception among the different densities of cotton, the grid math operation in Surfer software (Golden Software Inc., Golden, CO, USA) was used:
Grid C(i, j) = Grid A(i, j) − Grid B(i, j)
The fraction of intercepted PAR (iPAR), PAR reflectance (rPAR) and PAR transmittance (tPAR) at each measurement point in the cotton canopy were calculated according to Tang’s method [29]. tPAR was determined by dividing transmitted PAR intensity by total incident PAR intensity, and rPAR was calculated by dividing reflected PAR intensity by total incident PAR intensity.
According to a study by Zhi [23], under a planting density of 60,000 plants/ha, light reflectance of the cotton plants is only 3–5%, which has little effect on the growth and development of cotton. Therefore, in the current study, the influence of reflected light on the growth and development of cotton was ignored. Therefore, the iPAR was calculated according to the formula below:
iPAR = 1 − tPAR − rPAR
The cumulative interception rate and transmittance in the canopy of the entire population can be calculated in Surfer software with the extension of Simpson’s 3/8 rule. The formula is as follows:
A i = 3 Δ X 8 ( G i , 1 + 3 G i , 2 + 3 G i , 3 + 2 G i , 4 + + 2 G i , n c o l 1 + G i , n c o l )
V 3 Δ y 8 ( A 1 + 3 A 2 + 3 A 3 + 2 A 4 + + 2 A n c o l 1 + A n c o l )
where Gij is the ith row and jth column of the grid point excavation depth; Ai is the ith cross-sectional area; and ΔX and Δy are the grid data file column and row distances, respectively. The coefficients for Equations (3) and (4) are [1,3,3,2,……3,3,2,1].

2.3.2. Leaf Area

The LAI of the cotton plants was measured by collecting all leaves from two randomly selected plants growing in each plot and using a flatbed scanner (Phantom p800xl, MiCROTEK, Shanghai, China) to take photos of the leaves. The leaf images were processed with Image-Pro Plus 7.0 software. (Media Cybernetics, Rockville, MD, USA). Then, the total leaf area of a single plant was calculated. The LAI was determined by dividing total leaf area per plant by corresponding ground area.

2.3.3. Plant Height

Two plants were randomly selected and uprooted from the middle rows in each plot, and their height was measured using a ruler. The length from the cotyledon to the growth point of the plant was considered the height of a single plant. Average plant height was calculated with the data from all the selected plants from each treatment.

2.3.4. Biomass

The dry weight of the cotton plants was determined by collecting samples at an interval of approximately 15 days during the growing season at 30, 45, 60, 79, 87, 102 and 117 DAE in 2018 and at 30, 45, 60, 74, 88, 99, 116 and 129 DAE in 2019. Two plants from each plot were uprooted and separated into underground parts (roots), leaves, stems and reproductive parts. The samples were placed into an electric fan-assisted oven at 105 °C for 30 min to stop metabolism. The samples were then dried at 80 °C for 48 h until a constant weight was achieved. A logistic regression formula was used to calculate biomass accumulation [30]:
Y = A 1 + b e k t
In Formula (5), where Y (kg) represents biomass, A (kg) represents the maximum biomass, t (d) is the DAE, and b and k are constants.
The following formulas were calculated from Formula (5):
t 0 = lnb k ( t 0 = t )
t 1 = ln b ln ( 2 3 ) k
t 2 = ln b ln ( 2 + 3 ) k
V M = Ak 4
V A = Y 2 Y 1 Δ t
In the formulas above, VM represents the highest rate of biomass accumulation, t represents the fast accumulation period of maximum biomass, Y1 and Y2 represent the current biomass at t1 and t2, VA is the average biomass accumulation from t1 to t2, and ∆t is the total period of biomass accumulation.

2.3.5. Yield

Seed cotton yield and lint yield were measured after harvest, and the plants in each treatment were harvested three times each year. Lint percentage was calculated by dividing lint yield by seed cotton yield of 100 bolls.

2.3.6. Statistical Analysis

Microsoft Excel 2019 (Microsoft, Bothell, WA, USA) and Stata 16.0 (Stata Corporation, College Station, TX, USA) were used to process the data. SPSS 19.0 (SPSS Inc., Chicago, IL, USA) and Origin 2016 (OriginLab Corporation, Northampton, MA, USA) were used to analyze the yield and biomass data. Image-Pro Plus 6.0 was used to calculate the LAI. Light interception data were processed and analyzed in Stata 16.0, Surfer 18 (Golden Software Inc., Golden, CO, USA) and GS + 9.0 (Gamma Design Software Inc, Plainwell, MI, USA). Duncan’s multiple range test at the 5% probability level was used to test the differences among the mean values.

3. Results

3.1. Light Interception

3.1.1. Light Interception at Different Points

The light interception data were collected and analyzed the different positions of iPAR and the variation in light interception during the flowering stage to assess the competition for light among the different densities of cotton plants. According to Figure 3, light interception varied at different canopy positions. In 2018, the trend of iPAR ranged from 0.45 to 0.80 (low densities), 0.50 to 0.83 (medium densities), and 0.50 to 0.92 (high densities) at 50 cm in the vertical direction. iPAR decreased from bottom to the top of the cotton canopy, and there was no significant difference at the top of cotton canopies among different plant densities.
At each layer of the cotton canopies, light interception near the cotton rows was significantly higher than at the inter-row of cotton, with a value of 0.6–0.85 near the cotton rows and 0.45–0.75 at the inter-row of cotton. In addition, there was significant difference of light interception near cotton rows. However, there was no significant difference of light interception at the inter-row of cotton canopies, as the iPAR ranged from 0 to 0.85 in all plant densities at 50 cm in the horizontal direction. Additionally, light interception increased with increasing plant populations, with the trend of iPAR ranges from 0 to 0.85 (low densities), 0 to 0.90 (medium densities), and 0 to 0.95 (high densities). The result of 2019 showed the same trend of 2018, with the iPAR ranging from 0.42 to 0.85 (low densities), 0.50 to 0.97 (medium densities), and 0.44 to 0.98 (high densities) near the cotton rows at 30–80 cm in the vertical direction.

3.1.2. Variation of Light Interception in Different Positions

iPAR varied at different cotton plant densities, with the highest iPAR observed at high densities. Figure 4 and Figure 5 shows the variation in light interception by the whole cotton canopy associated with different plant densities for the two years, and each figure shows the iPAR differences between adjacent densities. According to the contour plots for the two years, the high densities of cotton plants showed a higher fraction of light interception than the low populations when adjacent densities were considered, with the value of iPAR variation ranging from −0.005 to 0.150. On the basis of vertical distance, most variation in the iPAR ranged from 0 to 0.02 during the two evaluated years, except for D5-D4 and D6-D5 in 2019, for which the range was 0–0.045. At the level of 40–80 cm, iPAR variation was 0–0.14 for D2-D1, 0–0.055 for D3-D2, D4-D3 and D5-D4 and 0–0.125 for D6-D5. At a height below 40 cm, iPAR variations were 0–0.05 for all comparisons. In addition, the results show that there were no significant differences at a height above 80 cm, while at a height below 80 cm, high-densities cotton plants showed a higher fraction light interception than the plants growing at a lower density. In terms of the horizontal position, most of iPAR variation values at a position between 30 and 50 cm ranged from −0.005 to 0.05 in 2018 and from −0.05 to 0.02 in 2019, showing that the variation in light interception by cotton plants did not obviously differ at the midpoint of the rows. However, most of the variation in iPAR at positions below 30 cm and above 80 cm ranged from −0.005 to 0.15 during the two evaluated years, indicating that there was little difference in light interception among the different plant densities, as the high densities of cotton plants were associated with high light interception.

3.2. Leaf Area

The total leaf area of a single plant and the LAI changed linearly on different DAE during the cotton-growing period. Total leaf area values of a single plant at different planting densities are shown in Figure 6.
Total leaf area of a single plant first increased, reached a peak at approximately 100 DAE, and then decreased over the course of the growth period in both years. Before 80 DAE in both years there was no significant difference between different plant densities in terms of leaf area. However, after 80 DAE, leaf area at D1 was higher than that at any other density during both years, and leaf area decreased as plant density increased.
The plants growing at the highest density reached a maximum leaf area of 4292.93 cm2 in 2018 and 4515.84 cm2 in 2019, while those growing at the lowest density reached a maximum leaf area of 6288.16 cm2 in 2018 and 6745.95 cm2 in 2019.
The LAI at D6 was higher than that at any other density during both years over the course of the whole growing period, and the LAI increased as plant density increased. The LAI first increased with the growth of the plants in both years, reached a peak at approximately 102 DAE and then decreased (Figure 7).
The plants growing at the highest density reached a maximum LAI of 4.18 in 2018 and 4.25 in 2019, while those growing at the lowest density reached a maximum LAI of 1.20 in 2018 and 1.31 in 2019.

3.3. Plant Height

Cotton plant height increased during the whole growing period, with plant heights in the order of D6 > D5 > D4 > D3 > D2 > D1 during both years during the whole growing period. The results for the two years showed the same trends in plant height (Figure 8).
The plants growing at the highest density reached a maximum plant height of 95.33 cm at 87 DAE in 2018 and 121.83 cm at 107 DAE in 2019, while those growing at the lowest density reached a maximum height of 77.83 cm at 87 DAE in 2018 and 100.33 cm at 107 DAE in 2019.

3.4. Biomass

3.4.1. Aboveground Biomass of Individual Plants

Cotton plant aboveground biomass associated with each density followed a normal logistic model, and that at the lowest density (D1) was higher than that at any other density during both years over the whole growing period (Figure 9).
The plants growing at the lowest density reached a maximum aboveground biomass of 345.06 g at 117 DAE in 2018 and 464.21 g at 129 DAE in 2019, while those growing at the highest density reached an aboveground biomass of 87.35 g at 117 DAE in 2018 and 110.85 g at 129 DAE in 2019.

3.4.2. Aboveground Biomass Accumulation

The cotton aboveground biomass significantly followed a normal logistic model on the basis of DAE (Table 3 and Table 4), and the aboveground biomass of the D6 population was higher than that at any other density during both years over the whole growing period (Figure 10). The simulation of biomass accumulation was based on Formula (5). Table 3 illustrates the starting and ending days of the fast accumulation period during 2018 and 2019, and the maximum aboveground biomass of all plants (1149.106 gm m−2 d−1) can be calculated. The average highest speeds of aboveground biomass accumulation for all plants were 55 and 102 DAE in 2018 and 61 and 93 DAE in 2019, with the highest average (VA = 11.5 and 17.0 gm m−2 d−1) and maximum rates (VM = 13.1 and 19.4 gm m−2 d−1) (Table 3).

3.5. The Relationships and Fitted Models

In this study, light interception data of each plot were calculated to observe the relationships with leaf area, aboveground biomass and plant height. As shown in Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15, across the different densities of cotton, the aboveground biomass of single plant was linearly related to cumulative iPAR in both years, with R2 values of 0.95 and 0.91, respectively, there was also a linear relationship between iPAR and plant height in this study, with R2 values of 0.93 and 0.91, respectively. The leaf area of single plant had a highly significant logarithmic correlation with iPAR, with R2 values of 0.89 and 0.87, respectively. For cotton populations, the LAI had a highly significant logarithmic correlation with iPAR, with R2 values of 0.83 and 0.74, respectively. The aboveground biomass were linearly related to cumulative iPAR, with R2 values of 0.93 and 0.92.

3.6. Yield

In both years, the D5 plant density resulted in the highest production of seed cotton, followed by the highest plant density (D6). The highest seed cotton yields were 4232.2 and 4351.5 kg ha−1 in 2018 and 2019, respectively, which were 14.9%, 15.0%, 25.6%, 37.4% and 54.5% higher than those in D6, D4, D3, D2 and D1, respectively, in 2018, and 4.3%, 6.2%, 7.3%, 15.5% and 31.9% higher than those in D6, D4, D3, D2 and D1, respectively, in 2019. The highest lint yields were 1545.3 and 1789.5 kg ha−1, produced at D5, in 2018 and 2019, respectively, which were 12.8%, 17.3%, 19.7%, 23.6% and 30.3% higher than those produced at D4, D6, D3, D2 and D1, respectively, in 2018, and 4.5%, 8.5%, 15.3%, 23.2% and 37.0% higher than those produced at D6, D4, D3, D2 and D1, respectively, in 2019. The highest lint percentage (43.6%) was recorded in D1, followed by D2, D3, D4 and D5 (Table 5).
The spatial distribution of cotton yield varied according to the different plant densities and at the different positions (Figure 16). With the exception of D6, a high plant density resulted in a high number of nodes with cotton bolls. As cotton plant density increased, the order of the number of nodes with cotton bolls was D1 > D2 > D3 > D6 > D4 > D5, and the distribution of cotton bolls was concentrated at the vertical center of the cotton plants. In terms of the horizontal position, as cotton density increased, the distribution of cotton bolls was concentrated toward the inner side of the cotton plants.

4. Discussion

Spatial distribution of light interception by cotton growing at different densities is of great importance, as light interception directly affects photosynthesis and influences biomass and yield [31]. The traditional method to estimate the vertical direction of light interception in the crop canopy is Beer’s law [32], however, leaf angle and leaf curve are needed to simulate the model to calculate PAR distribution of different height of the canopies and contains a series of mathematical functions which are difficult to measure and calculate [33]. The 3-D structure models are also widely used to calculate the light interceptions of crop canopies, which need a variety of measurement and data analysis [34]. This study presents the spatial grid method to identify stable measurement points in the cotton canopy between two cotton rows, which established a spatial grid perpendicular to the cotton plants, and focused on the comparison of light interception by cotton plants. By calculating the fraction of intercepted PAR (iPAR), PAR reflectance (rPAR) and PAR transmittance (tPAR) at each measurement point in the cotton canopy according to Tang’s method and the extension of Simpson’s 3/8 rule, multiple grid files representing different densities can be combined into one file to analyze differences in light interception by different cotton plant populations.
During the vigorous period of cotton growth and development, one of the most important factors of cotton growth was density [35], which showed obvious regional patterns and a high degree of spatial autocorrelation [36]. The plants growing at higher densities experienced greater competition than those at the lower densities with smaller cotton plants, fewer buds and bolls and lower productivity [37], which demonstrated the effect of spatial correlation between crop growth and yield. In both years, the D5 plant density with the highest seed cotton yield and lint yield, followed by D6, D4, D3, D2 and D1. The number of fruit branches per plant and the number of bolls per plant decreased with the increase of density. With the increase of density, the proportion of upper boll, middle boll and inner boll increased, and the proportion of lower boll and outer boll decreased [38].
The distribution of light interception at different canopy positions is spatially variable. Additionally, when canopy height increased, light interception decreased, it might be affected by the elevation of the sun or the diffuse radiation received at the lower layers of canopy [39,40,41]. In the horizontal direction iPAR values were greater near the cotton rows than at the inter-row of the canopy. According to Xue’s research, light interception of horizonal layers was not significantly different at the late stage of growth period, and in the early stage showed obvious difference of light interception among different densities, which increased light-use efficiency [42]. In this study, the distribution and variation in light interception during the flowering stage were analyzed, and the result showed that light interception among six different cotton densities showed significant differences in vertical direction, which obviously differed at canopy heights below 80 cm near the side of the cotton row and was mainly affected by the upper position of the canopy [43].
The competition between plants at different densities influenced the number and size of the leaves and influenced the canopy size, which affected the iPAR and yield [44]. In the whole cotton population, LAI increased when plant densities increased, however, the leaf area of a single plant was significantly higher in the low densities, which may be a reason why the cotton yield of D6 was high but slightly lower than D5, for the excessive leaf area may cause serious shedding in the field and leaf abscission, which affect light interception. If leaf area is small, though the light conditions are satisfied, the total photosynthetic products cannot meet the needs of boll growth, both situations lead to low yield of cotton [45]. A highly significant relationship was found between iPAR and leaf area, which is in strong agreement with previous studies [14,24,31].
Light is one of the original sources of energy which directly affect crop growth. The accumulation of biomass was significantly related to radiation-use efficiency in different crops [46,47,48,49]. In this study, higher densities showed higher aboveground biomass and yield in population, which may be because a higher density of plants with greater light interception presents a higher speed and longer period of fast biomass accumulation, for single plant the D1 density of plant was significantly higher than others, the reason was that the plants in low densities have less competition which each other, and each plant can receive enough light interception for growth, but the whole plant population was less competitive than the high-density ones. The trend of aboveground biomass accumulation followed the same rules as in Khan’s study [30], and the aboveground biomass of a single plant and the whole population growing at different densities were both linearly related to cumulative iPAR [14,41]. However, low plant densities present lower cotton productivity, that may be because early leaf senescence may cause low light-use efficiency in a population, which furthermore reduces boll weight [50].
The plant can also increase light interception by increasing plant height [51]. It was shown that plant height had a linear relationship with iPAR, which was because plant height may directly affect light interception and changes in the wind field within the cotton population, and, furthermore, affect the shedding area and photosynthesis at the bottom of crop canopies [52].
Canopy structure may also play an important role by directly affecting PAR and photosynthesis, with a further influence on light interception by the canopy and plant yield [53,54,55]. With increasing cotton plant density, cotton leaf area increased, while plant height decreased, further improving the canopy structure and increasing biomass accumulation [56].
Conclusively, the cotton population with a higher density (8.7 plants m−2) performed best at light interception, with the highest light-use efficiency and the highest yield.

5. Conclusions

Light is of great importance during the cotton growing period, and for the different densities of cotton plants, it was necessary to study the spatial distribution of light interception. The results of this study showed that (i) iPAR of the whole cotton canopy was higher at a high cotton plant density than at a lower density, and in the horizontal direction, differences in light interception were obvious at canopy heights below 80 cm near the side of the cotton row, (ii) LAI and aboveground biomass were obviously correlated with light interception. (iii) The cotton population with a higher density (8.7 plants m−2) performed best at the light interception, with the highest yield. The results will provide a significant basis by which to improve the management of light and canopy architecture.

Author Contributions

Y.L. (Yabing Li) and Y.H. designed the experiment. W.D., B.Y., S.X., M.X. and J.A. conducted the experiment. Y.L. (Yaping Lei), F.X., L.F., X.L. and G.W. helped and provided useful suggestions during the experiment. H.C. and X.Z. processed and analyzed data and wrote the first draft. Y.L. (Yabing Li), L.F. and Z.W. revised and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

We are thankful for the financial support by National Natural Foundation of China (31371561). The funder had no role in design, data collection and decision to publish or preparation of manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We highly acknowledge the help of technicians of the research station from the institute of cotton research Chinese Academy of Agricultural Sciences.

Conflicts of Interest

The authors declare no conflict of interests.

References

  1. Lu, X.; Jia, X.; Niu, J. The present situation and prospects of cotton industry development in China. Sci. Agric. Sin. 2018, 51, 26–36. [Google Scholar]
  2. Chen, M. Research on issues related to the development of China’s cotton industry--report at the 2014 annual meeting of China cotton association. China Cotton 2014, 41, 38–43. [Google Scholar]
  3. Mao, S.; Li, Y. Report on China’s Cotton Production Prosperity 2017; China Agriculture Press: Beijing, China, 2018. [Google Scholar]
  4. Zhang, B.; Bian, X. Evaluation of compresive benefits of smart agriculture in Heilongjiang province based on AHP. Chin. J. Agric. Resour. Reg. Plan 2019, 40, 109–115. [Google Scholar]
  5. Mao, S. The modernization of cotton cultivation technology in China also focuses on the research related to cotton cultivation in the 12th Five-Year Plan. China Cotton 2010, 37, 2–5. [Google Scholar]
  6. Dai, J.; Dong, H. Field Crops Research Intensive cotton farming technologies in China: Achievements, challenges and countermeasures. Field Crops Res. 2014, 155, 99–110. [Google Scholar] [CrossRef] [Green Version]
  7. Mao, S.; Li, Y.; Dong, H.; Bie, S.; Lin, Y.; Dong, H.; Tian, L.; Chen, D.; Ma, Q.; Ma, X.; et al. Institute of Cotton Research of CAAS. In Cotton Cultivation in China; Shanghai Scientific and Technical Publisher: Shanghai, China, 2013; pp. 123–124. [Google Scholar]
  8. Zhang, W.; Gou, L.; Li, M. Studies on The Relationship Between Canopy Apparent Photosynthesis Rate and Yield in Cotton in North Xinjiang. Cotton Sci. 1999, 11, 185–190. [Google Scholar]
  9. McCree, K.J. Photosynthetically Active Radiation. Physiological Plant Ecology I; Springer: Berlin/Heidelberg, Germany, 1981. [Google Scholar]
  10. McCree, K.J. Test of current definitions of photosynthetically active radiation against leaf photosynthesis data. Agric. Meteorol. 1972, 10, 443–453. [Google Scholar] [CrossRef]
  11. Gent, M.P.N. Canopy light interception, gas exchange, and biomass in reduced height isolines of winter wheat. Crop. Sci. 1995, 35, 1636–1642. [Google Scholar] [CrossRef]
  12. Zhang, L.Z.; Van der Werf, W.; Bastiaans, L.; Zhang, S.; Li, B.; Spiertz, J.H.J. Light interception and utilization in relay intercrops of wheat and cotton. Field Crops Res. 2008, 107, 29–42. [Google Scholar] [CrossRef]
  13. Trouwborst, G.; Oosterkamp, J.; Hogewoning, S.W.; Harbinson, J.; Van Ieperen, W. The responses of light interception, photosynthesis and fruit yield of cucumber to LED-lighting within the canopy. Physiol. Plant. 2010, 138, 289–300. [Google Scholar] [CrossRef]
  14. Zhi, X.; Han, Y.; Xing, F.; Lei, Y.; Wang, G.; Feng, L.; Yang, B.; Wang, Z.; Li, X.; Xiong, S. How do cotton light interception and carbohydrate partitioning respond to cropping systems including monoculture, intercropping with wheat, and direct-seeding after wheat? PLoS ONE 2019, 14, e0217243. [Google Scholar] [CrossRef] [PubMed]
  15. Louarn, G.; Lecoeur, J.; Lebon, E. A three-dimensional statistical reconstruction model of grapevine (Vitis vinifera) simulating canopy structure variability within and between cultivar/training system pairs. Ann. Bot. 2008, 101, 1167–1184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Vargas, L.A.; Andersen, M.N.; Jensen, C.R.; Jorgensen, U. Estimation of leaf area index, light interception and biomass accumulation of Miscanthus sinensis ‘Goliath’ from radiation measurements. Biomass Bioenergy 2002, 22, 1–14. [Google Scholar] [CrossRef]
  17. Heitholt, J.J.; Pettigrew, W.T.; Meredith, W.R. Light Interception and Lint Yield of Narrow-Row Cotton. Crop. Sci. 1992, 32, 728–733. [Google Scholar] [CrossRef]
  18. Wells, R. Soybean growth response to plant density: Relationships among canopy photosynthesis, leaf area, and light interception. Crop. Sci. 1991, 31, 755–761. [Google Scholar] [CrossRef]
  19. Goudriaan, J.; Monteith, J.L. A mathematical function for crop growth based on light interception and leaf area expansion. Ann. Bot. 1990, 66, 695–701. [Google Scholar] [CrossRef] [Green Version]
  20. Maddonni, G.A.; Otegui, M.E. Leaf area, light interception, and crop development in maize. Field Crops Res. 1996, 48, 81–87. [Google Scholar] [CrossRef]
  21. Du, X.; Chen, B.; Shen, T.; Zhang, Y.; Zhou, Z. Effect of cropping system on radiation use efficiency in double-cropped wheat–cotton. Field Crops Res. 2015, 170, 21–31. [Google Scholar] [CrossRef]
  22. Chapepa, B.; Mudada, N.; Mapuranga, R. The impact of plant density and spatial arrangement on light interception on cotton crop and seed cotton yield: An overview. J. Cotton Res. 2020, 3, 1–6. [Google Scholar] [CrossRef]
  23. Bhattacharya, A. Changing Climate and Resource Use Efficiency in Plants; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  24. Zhi, X.; Han, Y.; Mao, S.; Wang, G.; Feng, L.; Yang, B.; Fan, Z.; Du, W.; Lu, J.; Li, Y. Light spatial distribution in the canopy and crop development in cotton. PLoS ONE 2014, 9, e113409. [Google Scholar] [CrossRef]
  25. Vos, J.; Evers, J.B.; Buck-Sorlin, G.H.; Andrieu, B.; Chelle, M.; de Visser, P.H.B. Functional-structural plant modelling: A new versatile tool in crop science. J. Exp. Bot. 2010, 61, 2101–2115. [Google Scholar] [PubMed]
  26. Rey, H.; Dauzat, J.; Chenu, K.; Barczi, J.F.; Dosio, G.A.A.; Lecoeur, J. Using a 3-D virtual sunflower to simulate light capture at organ, plant and plot levels: Contribution of organ interception, impact of heliotropism and analysis of genotypic differences. Ann. Bot. 2008, 101, 1139–1151. [Google Scholar] [CrossRef] [PubMed]
  27. Gu, S.; Evers, J.B.; Zhang, L.; Mao, L.; Zhang, S.; Zhao, X.; Liu, S.; van der Werf, W.; Li, Z. Modelling the structural response of cotton plants to mepiquat chloride and population density. Ann. Bot. 2014, 114, 877–887. [Google Scholar] [CrossRef] [Green Version]
  28. Mao, L.; Zhang, L.; Evers, J.B.; Henke, M.; van der Werf, W.; Liu, S.; Zhang, S.; Zhao, X.; Wang, B.; Li, Z. Identification of plant configurations maximizing radiation capture in relay strip cotton using a functional–structural plant model. Field Crops Res. 2016, 187, 1–11. [Google Scholar] [CrossRef]
  29. Tang, Q. Studys on Canopy Photosynthesis Characteristic of Cotton with Different Leaf-Type; Shihezi University: Shihezi, China, 2014. [Google Scholar]
  30. Khan, N.; Xing, F.; Feng, L.; Wang, Z.; Xin, M.; Xiong, S.; Wang, G.; Chen, H.; Du, W.; Li, Y. Comparative Yield, Fiber Quality and Dry Matter Production of Cotton Planted at Various Densities under Equidistant Row Arrangement. Agronomy 2020, 10, 232. [Google Scholar] [CrossRef] [Green Version]
  31. Xing, F.; Han, Y.; Feng, L.; Zhi, X.; Wang, G.; Yang, B.; Fan, Z.; Lei, Y.; Du, W.; Wang, Z.; et al. Genotypic variation in spatiotemporal distribution of canopy light interception in relation to yield formation in cotton. J. Cotton Res. 2018, 1, 13. [Google Scholar] [CrossRef]
  32. Monsi, M.; Saeki, T. Über den Lichtfaktor in den Pflanzengesellschaften und seine Bedeutung für die Stoffproduktion. Jpn. J. Bot. 1953, 14, 22–52. [Google Scholar]
  33. Zhang, W.Y.; Tang, L.; Yang, X.; Liu, L.L.; Cao, W.X.; Zhu, Y. A simulation model for predicting canopy structure and light distribution in wheat. Eur. J. Agron. 2015, 67, 1–11. [Google Scholar] [CrossRef]
  34. Cao, W.X. Digital Farming Technology; Science Press: Beijing, China, 2018. [Google Scholar]
  35. Li, Y.; Mao, S.; Wang, X.; Han, Y.; Wang, G.; Fan, Z.; Zhu, Q. Cotton Growth Condition Spatial Variance Analysis Based on Geo-statistics. Cotton Sci. 2007, 19, 214–219. [Google Scholar]
  36. Liu, Y.; Pan, X. The Analysis on Dominant Provinces of Cotton Production in China. Chin. Agric. Sci. Bull. 2006, 22, 360–364. [Google Scholar]
  37. Wang, Y.; Cui, X.; Pan, X.; Pan, S. Study on the effect of competition among cotton plants on light energy utilization and growth and development. China Cotton 1999, 08, 3–5. [Google Scholar]
  38. Li, F. Study on the spatial distribution of cotton bells of different planting densities in SCRC28. Cotton Sci. 2012, 34, 25–27. [Google Scholar]
  39. Leuchner, M.; Hertel, C.; Menzel, A. Spatial variability of photosyntheticallyactive radiation in European beech and Norway spruce. Agric. For. Meteorol. 2011, 151, 1226–1232. [Google Scholar] [CrossRef]
  40. Sarlikioti, V.; de Visser, P.H.; Marcelis, L.F. Exploring the spatial distributionof light interception and photosynthesis of canopies by means of afunctional-structural plant model. Ann. Bot. 2011, 107, 875–883. [Google Scholar] [CrossRef] [Green Version]
  41. Zhi, X.; Han, Y.; Wang, G.; Feng, L.; Yang, B.; Fan, Z.; Du, W.; Lei, Y.; Mao, S.; Li, Y. Changes to the PAR Spatial Distribution, Biomass, and Fiber Quality in Response to Plant Densities. Cotton Sci. 2017, 29, 365–373. [Google Scholar]
  42. Xue, H.; Han, Y.; Li, Y.; Wang, G.; Feng, L.; Fan, Z.; Du, W.; Yang, B.; Cao, C.; Mao, S. Spatial distribution of light interception by different plant population densities and its relationship with yield. Field Crops Res. 2015, 184, 17–27. [Google Scholar] [CrossRef]
  43. Han, C.; Wei, S. The study of canopy characteristics of high-yielding groups in rice and their relationship with yield. Liaoning Agric. Sci. 1996, 3, 24–27. [Google Scholar]
  44. Bai, Z. Study on light interception and biomass production of different cotton cultivars. PLoS ONE 2016, 11, e0156335. [Google Scholar] [CrossRef]
  45. Shi, W. Effects of Plant Densities in Cotton (Gossypium hirsutum L.) Dry Matter Accumulation and Distribution, Yield and Fiber Quality; Nanjing Agricultural University: Nanjing, China, 2012. [Google Scholar]
  46. Kiniry, J.R.; Jones, C.A.; O’toole, J.C.; Blanchet, R.; Cabelguenne, M.; Spanel, D.A. Radiation-use efficiency in biomass accumulation prior to grain-filling for five grain-crop species. Field Crops Res. 1989, 20, 51–64. [Google Scholar] [CrossRef]
  47. Flénet, F.; Kiniry, J.R.; Board, J.E.; Westgate, M.E.; Reicosky, D.C. Rows pacing effects on light extinction coefficients of corn, sorghum soybean and sunflower. Agron. J. 1996, 88, 185–190. [Google Scholar] [CrossRef]
  48. Kiniry, J.R.; Williams, J.R.; Vanderlip, R.L.; Atwood, J.D.; Reicosky, D.C.; Mulliken, J.; Cox, W.J.; Mascagni, H.J., Jr.; Hollinger, S.E.; Wiebold, W.J. Evaluation of two maize models for nine U.S. locations. Agron. J. 1997, 89, 421–426. [Google Scholar] [CrossRef]
  49. Connor, D.J.; Fereres, E. A dynamic model of crop growth and partitioning of biomass. Field Crops Res. 1999, 63, 139–157. [Google Scholar] [CrossRef]
  50. Dai, J.; Li, W.; Tang, W.; Zhang, D.; Li, Z.; Lu, H.; Eneji, A.E.; Dong, H. Manipulation of dry matter accumulation and partitioning with plant density in relation to yield stability of cotton under intensive management. Field Crops Res. 2015, 180, 207–215. [Google Scholar] [CrossRef]
  51. Francescangeli, N.; Sangiacomo, M.A.; Martí, H. Effects of plant density in broccoli on yield and radiation use efficiency. Sci. Hortic. 2006, 110, 135–143. [Google Scholar] [CrossRef]
  52. Cheng, G.; Sun, H.; Zhang, Y.; Liu, L.; Wang, M.; Li, C. Effects of canopy structure on photosynthesis and yield in cotton. J. Agric. Univ. Hebei 2015, 38, 1–7. [Google Scholar]
  53. Lv, X.; Wang, K.; Liu, J.; Chu, G.; Yang, Y.; Dun, J.; Yao, C.; Zhou, F. An Analysis of Characters on Wheat/Corn Compound Canopy Structure. J. Shihezi Univ. (Nat. Sci.) 2001, 5, 120–122. [Google Scholar]
  54. Khan, N.; Han, Y.; Xing, F.; Feng, L.; Wang, Z.; Wang, G.; Yang, B.; Fan, Z.; Lei, Y.; Xiong, S.; et al. Plant Density Influences Reproductive Growth, Lint Yield and Boll Spatial Distribution of Cotton. Agronomy 2020, 10, 14. [Google Scholar] [CrossRef] [Green Version]
  55. Liu, F. Experimental and Mathematical Model of Rice Canopy Structure and Dry Matter Production under Controlled Irrigation; Hehai University: Nanjing, China, 2005. [Google Scholar]
  56. Zhang, D. Effect of Cotton Plant Densities on Light Distribution and LUE of Jujube-Cotton Intercropped Systems; Shihezi University: Shihezi, China, 2014. [Google Scholar]
Figure 1. Temperature during the cotton growing season in 2018 and 2019 and the average for the preceding five years.
Figure 1. Temperature during the cotton growing season in 2018 and 2019 and the average for the preceding five years.
Agronomy 11 00176 g001
Figure 2. The measurement points in the two cotton rows of cotton plot.
Figure 2. The measurement points in the two cotton rows of cotton plot.
Agronomy 11 00176 g002
Figure 3. Light interception during the flowering stage at different cotton plant densities in 2018 and 2019. (A) D1; (B) D2; (C) D3; (D) D4; (E) D5; (F) D6.
Figure 3. Light interception during the flowering stage at different cotton plant densities in 2018 and 2019. (A) D1; (B) D2; (C) D3; (D) D4; (E) D5; (F) D6.
Agronomy 11 00176 g003aAgronomy 11 00176 g003b
Figure 4. Variation in light interception at different positions in the canopy of plants growing at different densities in 2018 and 2019. (a) Near the cotton rows; (b) midpoint; (A) D2-D1; (B) D3-D2; (C) D4-D3; (D) D5-D4; (E) D6-D5.
Figure 4. Variation in light interception at different positions in the canopy of plants growing at different densities in 2018 and 2019. (a) Near the cotton rows; (b) midpoint; (A) D2-D1; (B) D3-D2; (C) D4-D3; (D) D5-D4; (E) D6-D5.
Agronomy 11 00176 g004
Figure 5. Variation in light interception at different levels in the canopy of plants growing at different densities in 2018 and 2019. (a) Bottom of the canopy; (b) middle of the canopy; (c) top of the canopy; (A) D2-D1; (B) D3-D2; (C) D4-D3; (D) D5-D4; (E) D6-D5.
Figure 5. Variation in light interception at different levels in the canopy of plants growing at different densities in 2018 and 2019. (a) Bottom of the canopy; (b) middle of the canopy; (c) top of the canopy; (A) D2-D1; (B) D3-D2; (C) D4-D3; (D) D5-D4; (E) D6-D5.
Agronomy 11 00176 g005
Figure 6. Leaf area per cotton plant growing at different planting densities in 2018 and 2019.
Figure 6. Leaf area per cotton plant growing at different planting densities in 2018 and 2019.
Agronomy 11 00176 g006
Figure 7. Leaf area index of cotton plants growing at different planting densities in 2018 and 2019.
Figure 7. Leaf area index of cotton plants growing at different planting densities in 2018 and 2019.
Agronomy 11 00176 g007aAgronomy 11 00176 g007b
Figure 8. Cotton plant height at different planting densities in 2018 and 2019.
Figure 8. Cotton plant height at different planting densities in 2018 and 2019.
Agronomy 11 00176 g008
Figure 9. The aboveground biomass of a single cotton plant grown at different planting densities in 2018 and 2019.
Figure 9. The aboveground biomass of a single cotton plant grown at different planting densities in 2018 and 2019.
Agronomy 11 00176 g009
Figure 10. The aboveground biomass of cotton grown at different planting densities in 2018 and 2019.
Figure 10. The aboveground biomass of cotton grown at different planting densities in 2018 and 2019.
Agronomy 11 00176 g010
Figure 11. Relationships and fitted models of intercepted PAR (iPAR) and leaf area per plant in 2018 and 2019.
Figure 11. Relationships and fitted models of intercepted PAR (iPAR) and leaf area per plant in 2018 and 2019.
Agronomy 11 00176 g011
Figure 12. Relationships and fitted models for the cumulative PAR interception and aboveground biomass per plant in 2018 and 2019.
Figure 12. Relationships and fitted models for the cumulative PAR interception and aboveground biomass per plant in 2018 and 2019.
Agronomy 11 00176 g012
Figure 13. Relationships and fitted models for the iPAR and leaf area index (LAI) per plant in 2018 and 2019.
Figure 13. Relationships and fitted models for the iPAR and leaf area index (LAI) per plant in 2018 and 2019.
Agronomy 11 00176 g013
Figure 14. Relationships and fitted models for the cumulative PAR interception and cotton plant aboveground biomass in 2018 and 2019.
Figure 14. Relationships and fitted models for the cumulative PAR interception and cotton plant aboveground biomass in 2018 and 2019.
Agronomy 11 00176 g014
Figure 15. Relationships and fitted models for iPAR and plant height in 2018 and 2019.
Figure 15. Relationships and fitted models for iPAR and plant height in 2018 and 2019.
Agronomy 11 00176 g015
Figure 16. Spatial distribution of average cotton yield at different planting densities averaged across in 2018 and 2019. (A) D1; (B) D2; (C) D3; (D) D4; (E) D5; (F) D6.
Figure 16. Spatial distribution of average cotton yield at different planting densities averaged across in 2018 and 2019. (A) D1; (B) D2; (C) D3; (D) D4; (E) D5; (F) D6.
Agronomy 11 00176 g016
Table 1. Sunshine hours and daily photosynthetically active radiation (PAR) during the cotton growing season in 2018 and 2019.
Table 1. Sunshine hours and daily photosynthetically active radiation (PAR) during the cotton growing season in 2018 and 2019.
Meteorological DataYearAprilMayJuneJulyAugustSeptemberOctober
Sunshine hours (h)2018198.59224.98248.91213.64234.48188.62205.45
2019297.67256.73260.16186.56213.74139.9199.93
Daily PAR (μmol·m−2·s−1)2018348.58398.81469.74418.51388.82310.04234.28
2019313.95490.95487.20448.48377.63304.86203.47
Table 2. Cotton growth stages according to the different plant densities in 2018 and 2019.
Table 2. Cotton growth stages according to the different plant densities in 2018 and 2019.
TreatmentEmergence/DAE 1Seedling/DAESquaring/DAEFlowering/DAEBoll Opening/DAE
Year 2018
D1April 28/0May 17/19May 29/31June 24/57August 6/100
D2April 28/0May 17/19May 30/32June 25/58August 8/102
D3April 28/0May 17/19May 30/32June 25/58August 8/102
D4April 28/0May 17/19May 30/32June 25/58August 10/104
D5April 28/0May 17/19May 30/32June 25/58August 10/104
D6April 28/0May 17/19May 31/33June 27/60August 13/107
Year 2019
D1April 28/0May 20/22May 28/30June 24/57August 14/108
D2April 28/0May 20/22May 28/30June 24/57August 15/109
D3April 28/0May 20/22May 30/32June 24/57August 15/109
D4April 28/0May 20/22May 30/32June 24/57August 16/110
D5April 28/0May 20/22May 31/33June 26/59August 17/111
D6April 28/0May 20/22May 31/33June 26/59August 19/112
1 DAE: Day after emergence.
Table 3. Eigen values of cotton biomass accumulation in the 2018 and 2019 growing season.
Table 3. Eigen values of cotton biomass accumulation in the 2018 and 2019 growing season.
TreatmentA (gm m−2 d−1)Fast Accumulation PeriodFastest Accumulation Point
t1 (DAE)t2 (DAE)Δt (d)VA (gm m−2 d−1)VM (gm m−2 d−1)at DAE
Year 2018
D1609.10564.0108.044.18.09.186.0
D2784.54258.4100.341.910.812.379.4
D3936.00354.9104.649.610.912.479.7
D41025.80554.5101.747.212.614.378.1
D51130.50551.3100.349.013.315.275.8
D61149.10648.999.450.513.115.074.2
Average55.3102.447.011.513.178.9
Year 2019
D1716.27464.7100.535.811.113.282.6
D2843.62663.094.131.115.617.878.6
D3918.05060.293.233.016.018.376.7
D4982.47059.390.731.418.820.675.0
D51091.45758.691.032.419.322.274.8
D61177.60658.290.432.221.024.174.3
Average60.793.332.617.019.477.0
t1 is the starting point and t2 is the ending point of the fast accumulation period (FAP). Dt is the total duration of the FAP. VA is the average and VM is the maximum rate of biomass accumulation during the FAP. DAE represents days after emergence.
Table 4. Simulation of plant biomass accumulation according to plant density during the cotton growing seasons of 2018 and 2019.
Table 4. Simulation of plant biomass accumulation according to plant density during the cotton growing seasons of 2018 and 2019.
TreatmentRegression EquationR2
Year 2018
D1Y = 6091.05095/(1 + 170.72e−0.05978t)0.985 ***
D2Y = 7845.41809/(1 + 147.37e−0.0629t)0.993 ***
D3Y = 9360.03681/(1 + 68.70e−0.05305t)0.992 ***
D4Y = 10,258.05161/(1 + 78.27e−0.05585t)0.991 ***
D5Y = 11,305.05257/(1 + 59.01e−0.05377t)0.989 ***
D6Y = 11,491.06419/(1 + 47.72e−0.05211t)0.989 ***
Year 2019
D1Y = 7162.73984/(1 + 433.56e−0.07351t)0.998 ***
D2Y = 8436.2559/(1 + 772.58e−0.08463t)0.997 ***
D3Y = 9180.49668/(1 + 455.31e−0.07978t)0.997 ***
D4Y = 9824.69778/(1 + 544.36e−0.08399t)0.996 ***
D5Y = 10,914.569/(1 + 439.83e−0.0814t)0.997 ***
D6Y = 11,776.0588/(1 + 436.68e−0.08186t)0.996 ***
*** Significant at the 0.001 probability level.
Table 5. Comparison of seed cotton and lint yields at various planting densities during the cotton growing seasons in 2018 and 2019.
Table 5. Comparison of seed cotton and lint yields at various planting densities during the cotton growing seasons in 2018 and 2019.
TreatmentSeed Cotton Yield (kg ha−1)Lint Yield (kg ha−1)Lint Percentage (%)
Year 2018
D12730.9 ± 96.0d1185.5 ± 32.0c43.6 ± 0.44a
D23080.5 ± 66.0cd1249.8 ± 27.4bc41.3 ± 0.58a
D33369.3 ± 84.5bc1291.0 ± 23.2bc38.4 ± 0.79a
D43679.5 ± 94.9b1369.4 ± 33.1b37.3 ± 0.28a
D54232.2 ± 116.1a1545.3 ± 28.2a36.5 ± 0.37a
D63682.2 ± 67.5b1316.8 ± 26.8bc35.8 ± 0.62a
Year 2019
D13299.1 ± 128.2d1306.4 ± 41.8d39.3 ± 1.07a
D23767.0 ± 68.4c1451.9 ± 23.3cd38.6 ± 1.15a
D34054.6 ± 20.2b1552.7 ± 10.4bc39.0 ± 1.10a
D44095.8 ± 72.8ab1649.4 ± 28.2ab39.1 ± 1.24a
D54351.5 ± 159.0a1789.5 ± 32.5a39.1 ± 1.31a
D64172.3 ± 105.4ab1712.8 ± 17.3ab38.7 ± 1.14a
Means followed by the same letters within the same category are statistically similar according to Duncan’s multiple range test at p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, H.; Zhao, X.; Han, Y.; Xing, F.; Feng, L.; Wang, Z.; Wang, G.; Yang, B.; Lei, Y.; Xiong, S.; et al. Competition for Light Interception in Cotton Populations of Different Densities. Agronomy 2021, 11, 176. https://doi.org/10.3390/agronomy11010176

AMA Style

Chen H, Zhao X, Han Y, Xing F, Feng L, Wang Z, Wang G, Yang B, Lei Y, Xiong S, et al. Competition for Light Interception in Cotton Populations of Different Densities. Agronomy. 2021; 11(1):176. https://doi.org/10.3390/agronomy11010176

Chicago/Turabian Style

Chen, Huanxuan, Xinxin Zhao, Yingchun Han, Fangfang Xing, Lu Feng, Zhanbiao Wang, Guoping Wang, Beifang Yang, Yaping Lei, Shiwu Xiong, and et al. 2021. "Competition for Light Interception in Cotton Populations of Different Densities" Agronomy 11, no. 1: 176. https://doi.org/10.3390/agronomy11010176

APA Style

Chen, H., Zhao, X., Han, Y., Xing, F., Feng, L., Wang, Z., Wang, G., Yang, B., Lei, Y., Xiong, S., Li, X., Xin, M., An, J., Du, W., & Li, Y. (2021). Competition for Light Interception in Cotton Populations of Different Densities. Agronomy, 11(1), 176. https://doi.org/10.3390/agronomy11010176

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