Improving maize grain yield by matching maize growth and solar radiation

Matching of maize growth with solar radiation is of great importance for achieving high yield. We conducted experiments using different maize cultivars and planting densities under different solar radiations during grain filling to quantitatively analyze the relationships among these factors. We found that a decrease in solar radiation after silking caused a drop in maize grain yield and biomass, with lower solar radiation intensities leading to worse grain yields and biomass. Cultivar ZD958 was more sensitive to solar radiation changes than cultivar XY335; slight decreases in solar radiation (i.e., 15% shading) caused significant declines in ZD958 grain yield. When total solar radiation during grain filling was less than 486.9 MJ m−2 for XY335 and less than 510.9 MJ m−2 for ZD958, the two cultivars demonstrated high yields at lower planting density of 7.5 × 104 plants ha−1; average yields were 13.36 and 11.09 Mg ha−1, respectively. When radiation intensities were higher than 549.5 MJ m−2 for XY335 and higher than 605.8 MJ m−2 for ZD958, yields were higher at a higher planting density of 12 × 104 plants ha−1, with average yields of 20.58 Mg ha−1 for XY335 and 19.65 Mg ha−1 for ZD958.

Food security is fundamental for human survival and national security. By 2050, the global population is expected to increase to 9.2 billion from 6.5 billion today 1 . Population growth and expanding consumption will drive global food demand up by 50% by 2030, or double that of 2000 levels 2 . A primary strategy for maintaining food security in China and the rest of the world is to ensure adequate food supplies. Studies have shown that China's future arable land is expected to decrease, whereas food consumption and demand are expected to continue increasing 3 . In the absence of major breakthroughs in farming systems and agricultural production, the most viable strategy to ensure food security is to continuously increase crop yields 4 . However, maize yield in China is lower than that in the Europe and the U.S. 5 . Adapting maize growth to the environment will be an important innovation which will match maize growth and the climate (solar radiation, temperature etc.) and increase total biomass and then grain yield. Solar radiation drives crop photosynthesis, the formation and development of plant organs, and crop yields. Solar radiation is a primary determinant of potential crop productivity [6][7][8][9] , so increases in crop yields can be achieved by improving radiation use efficiency. Maize is a high light-use efficiency crop, with more than 95% of its dry matter formed through photosynthesis 10 . Modern maize production is based on high planting densities 11 . Climatic conditions, especially solar radiation, significantly affect maize growth and optimal planting densities [12][13][14] . Maize cultivation is widespread in the world and solar radiation and planting densities vary significantly among different maize planting countries and regions [14][15][16] . For example, in China solar radiation in the west region is more abundant than that in eastern region. As a result, grain yields and planting densities in the west are higher than in the east 14,17 . In recent years, due to climate warming, global solar radiation has shown a decreasing trend, with average declines of 1.4-2.7% per decade [17][18][19][20] . China is one country that has experienced severe declines in solar radiation, with total solar radiation decreasing by 4.5 W m −2 per decade and effective sunshine hours decreasing by 1.28% per decade on average 17 . These declines are most severe in eastern China where total solar radiation has decreased by more than 6% per decade [18][19][20] . Declining solar radiation significantly affects maize growth at all growing stages, especially during grain filling and thus decreases grain yield 21,22 .
Most researches have shown the effects of solar radiation on maize growth by conducting shading experiments 13,[21][22][23] . Furthermore, different degrees and periods of shading have different effects on maize growth and development 23 . Previous studies had shown that shading caused maize grain yields to decrease significantly. In general, reductions in maize grain yields are greater with post-silking shading than before-silking shading 21,24,25 . Moreover, different maize genotypes show different responses to shading in terms of their growth and development. For maize varieties that are sensitive to weak light, the photosynthetic rate and fluorescence parameters of their leaves decrease greatly with shading, whereas varieties that are insensitive to weak light show only slight decreases in these two parameters. However, most previous studies were conducted in low solar radiation areas 13 , so any shading experiments performed in these areas would result in radiation levels below those that exist in any maize planting regions in China and other countries and therefore have little significance in real situations [21][22][23] . Furthermore, few studies have investigated the interactive and quantitative relationships between solar radiation, planting density, and genotype. Therefore, Qitai and Yinchuan where the solar radiation resources were the most abundant in China were chosen as the experimental sites in this study which would be more meaningful for studying the effects of different solar radiation on maize growth by conducting different shading level experiments. Through experiments of different solar radiation levels during grain filling, planting densities, and genotypes, we studied the interactions and quantitative relationships between these factors and explored optimal planting densities and genotypes under different solar radiation conditions. Our hypotheses are: (1) as solar radiation increases, optimal planting density also increases and (2) different genotypes respond differently to changes in solar radiation. This study can inform best practices for achieving high-yield maize cultivation in regions with different solar radiation intensities in China and other countries in the world as well as in the face of climate change.

Results
Effects of different shading levels on maize grain yield. The results of the current study suggest that shading leads to accelerated maize senescence. The duration of maize crop growth was longer at Qitai than at Yinchuan (Table 1) and was longer for ZD958 than for XY335, and the positions of two experimental sites were shown in Fig. 1. As shown in Table 2, the average temperature at Yinchuan was higher than at Qitai, while solar radiation was lower. Table 3 presents the specific experimental parameters. Experimental results at the two sites were averaged for 2016 and 2017. Under natural light condition (CK), XY335 grain yield was higher than that of ZD958 (Fig. 2). The grain yields of both maize cultivars decreased as shading increased, but shading affected maize grain yields differently depending on planting density. The S-15% shading treatment produced no significant difference in XY335 grain yield compared with CK under all planting density treatments. We also found no significant difference in XY335 grain yield between the S-30% and CK treatments at the D1 planting density. However, XY335 yield at S-50% and D1 was 34.7% less than under the CK condition. At the D2 planting density, XY335 yields at S-30% and S-50% were 19.2% and 41.8% lower than yields under the CK condition, respectively. At the D3 planting density, XY335 yields were 23.6% and 51.3% lower than CK yields at S-30% and S-50%, respectively. In contrast, ZD958 yields under all shading treatments were significantly lower than for CK at all planting densities. ZD958 yields under S-15%, S-30%, and S-50% were 13.3%, 19.3%, and 50.1% lower, respectively, than the CK condition at the D1 planting density. These values were 14.8%, 25.2%, and 54.7% lower than the CK condition under the D2 planting density and 13.7%, 28.7%, and 63.5% lower under the D3 planting density. Our findings that XY335 displayed no significant yield differences between CK, S-15%, and S-30% at low planting densities, but ZD958 displayed significantly lower yields at all shading treatments at all planting densities indicate that ZD958 is more sensitive to reductions in solar radiation.
In 2016, XY335 grain yields in plots planted at the D1 density and under the CK, S-15%, and S-30% conditions were significantly higher than the grain yield under S-50% at Qitai; however, when planted at the D2 and D3 densities, XY335 yields under the CK and S-15% conditions were significantly higher than under S-30% and S-50%. For ZD958, yields under the CK and S-15% conditions were significantly higher than under S-30% and S-50% at the D1 planting density; however, yields under the different shading treatments at the D2 and D3 planting densities decreased in the order CK > S-15% > S-30% > S-50%. Shading produced similar changes in yields for the two cultivars, although ZD958 yield reductions were more significant. At the same planting densities, ZD958 grain yields at Yinchuan under different shading treatments decreased in the order: CK > S-15% > S-30% > S-50%.
In 2017, grain yields for the two maize cultivars at both sites in 2017 showed similar trends. Yields at Qitai in 2017 were slightly higher than in 2016 for the same treatments, whereas the opposite trend was observed at Yinchuan.
Effects of different shading levels on maize biomass and leaf area index. The biomass of both maize cultivars decreased as shading increased (Table 4). Dry matter accumulations at the same planting densities and cultivar followed the order of S-CK > S-15% > S-30% > S-50%. As planting density increased, biomass increased for all shading treatments except S-50%. There was a significant difference in dry matter production for XY335 between the CK and S-50% treatments at the D1 planting density; at D2 and D3, dry matter production under the CK and S-15% treatments was significantly higher than under S-30% and S-50%. For ZD958, dry matter production at the D1 planting density under CK and S-15% was significantly higher than under the S-30% and S-50% treatments, while dry matter production at D2 and D3 under the different shading treatments decreased significantly in the order CK > S-15% > S-30% > S-50%. Dry matter production for the two cultivars showed a similar response to increased shading, although the dry matter reductions for ZD958 were more significant. The trend in dry matter production was similar to that of yield, which indicates that reductions in dry matter production with increased shading was the primary cause of reduced grain yields.   www.nature.com/scientificreports www.nature.com/scientificreports/ The LAI of the two cultivars showed similar drecreasing trends with increased shading, although the trend was more significant for ZD958 (Table 4). There were no significant differences in the LAIs of the two cultivars at the different shade treatments at the D1 planting density. When planted at the D2 density, the LAIs of the two cultivars at S-50% were significantly lower than those of the other treatments. At the D3 density, the LAIs of the two cultivars at S-30% and S-50% were significantly lower than those at CK and S-15%. This indicates that leaves are more sensitive to changes in solar radiation at high planting densities and that low solar radiation accelerates leaf senescence at high planting densities.   www.nature.com/scientificreports www.nature.com/scientificreports/ Genotype, planting density, solar radiation, and their interactive effects on maize yield. Figure 3 shows the interactive effects of planting density and solar radiation on maize grain yields. Grain yields of the two cultivars increased as solar radiation increased during grain filling. Under the same amount of solar radiation, XY335 yields were higher than ZD958 yields. Furthermore, ZD958 required more solar radiation to produce the same grain yield as XY335. This indicates that the solar radiation use efficiency of XY335 was much higher than that of ZD958. When solar radiation levels during grain filling were in the ranges of 350-370, 370-410, 410-460, 460-555, and 555-615 MJ m −2 , the corresponding XY335 yields were about 10, 12, 14, 16, and 18 Mg ha −1 , respectively. When the solar radiation during grain filling was higher than 615 MJ m −2 and the planting density was 7.5 × 10 4 -9.7 × 10 4 plants ha − We conclude that when the solar radiation during grain filling is less than 623 MJ m −2 , increases in maize yield are primarily related to the amount of solar radiation during grain filling. In other words, the greater the solar radiation during grain filling, the greater the grain yield. When the solar radiation during grain filling is higher than 623 MJ m −2 , increases in maize yield are primarily related to planting density. Therefore, low planting densities are more suitable for areas where the solar radiation during grain filling is less than 623 MJ m −2 . However, in areas receiving greater than 623 MJ m −2 solar radiation during grain filling, radiation use efficiency can be improved by increasing planting density. XY335 yields gradually increased with increasing planting density, and did so more significantly than ZD958 yields. When the solar radiation during grain filling was higher than 623 MJ m −2 , XY335 grain yield was higher than that of ZD958 at the same planting density. The maximum grain yield for ZD958 occurred at the planting density of 9.5 × 10 4 -10.5 × 10 4 plants ha −1 . However, when the planting density increased, yield decreased. This indicates that XY335 has a higher yield potential than ZD958 under sufficient solar radiation conditions. When solar radiation is not a constraint, XY335 grain yield could be improved by increasing planting density.    www.nature.com/scientificreports www.nature.com/scientificreports/ addition, we found that yield increases with increasing planting density until a maximum yield is reached. This maximum yield and corresponding density are determined by the amount of radiation. The higher the radiation, the higher the density and yield. In addition, XY335 required less solar radiation to produce the same grain yield as ZD959, which indicates that the radiation use efficiency of XY335 was higher than that of ZD958 (Fig. 4). . This demonstrates that under different solar radiation conditions, appropriate planting densities are a prerequisite for high yields. In summary, for XY335 and ZD958, when the solar radiation during grain filling was less than 486.9 MJ m −2 and 510.9 MJ m −2 , respectively, the appropriate planting density was 7.5 × 10 4 plants ha −1 ; when the soalr radiation during grain filling was 486.9-502.4 MJ m −2 and 510.9-552.6 MJ m −2 , respectively, the appropriate planting density was 10.5 × 10 4 plants ha −1 ; and when the solar radiation during grain filling was higher than 549.5 MJ m −2 and 605.8 MJ m −2 , respectively, the appropriate planting density was 12 × 10 4 plants ha −1 . Comparing the two cultivars, we found that under the same planting densities and solar radiation levels, ZD958 yields were lower than those of XY335 and the yield differences between different planting densities were greater for ZD958 than for XY335.

Discussion
Most studies on the effects of solar radiation on maize growth have examined shading treatments 15,[22][23][24][25] . These studies have shown that reductions in solar radiation resulted in reduced maize yields and that the maize yield reductions increased gradually with decreasing solar radiation 26 . The largest yield reductions in summer maize yields due to shading occurred during silking [27][28][29] . In this study, maize yield decreased with decreasing solar radiation. Under the S-50% treatment, grain yield was the lowest. In addition, high density planting is an effective way to increase maize yields in modern maize production, but excessively high planting densities can exacerbate inter-plant competition for limited resources and eventually reduce maize yields. This study shows that reductions in solar radiation at different planting densities leads to different degrees of reduced yield, and that yield reduction increased with increasing planting density. For example, grain yield at the D3 planting density was the lowest under the lowest solar radiation condition of S-50% treatment. This is mainly because of that after shading and increasing planting densities, solar radiation in the canopy decreased sharply and inter-plant competition for limited solar radiation resources was exacerbated, and photosynthetic capacity of leaves decreased especially the lower leaves in the canopy [14][15][16] .
While rice and wheat yields have improved largely because of an increasing harvest index since the Green Revolution, maize yield increases are mainly attributed to large increases in total aboveground dry matter accumulation 2,4,9 . Previous studies have shown that maximum total biomass decreases under shade and that yield and dry matter accumulation decreases to different degrees correspondingly 30 . This study demonstrates that shading caused dry matter production of the two cultivars to decrease, but the trend was more significant for ZD958. The decreasing trend in dry matter production was similar to that of grain yield, which indicates that dry matter production is a primary cause of changes in grain yield. LAI is an agronomic variable that summarizes the complex combination of internal and external factors that influence the structures of plant canopies, which are strongly influenced by interactions between plant genotype and environmental factors 31 . Previous studies have shown that maintaining a green canopy and delaying leaf senescence prolongs photosynthesis for grain filling, and that high planting densities can accelerate leaf senescence 32 . LAI and chlorophyll content have been shown to decrease after shading 28 . This study shows that LAI decreased under low light conditions. This indicates that leaves are more sensitive to changes in solar radiation especially at high planting densities and that low solar radiation at high planting densities can accelerate leaf senescence. LAI reduction of ZD958 was more significant than for XY335, which may have caused the earlier leaf senescence and lower yields observed for ZD958.
Practices that increase planting densities can effectively improve maize yields regionally and achieve high yields on a worldwide scale 33,34 . Maize planting densities are being increased year by year in many maize producing countries, including the United States and China 35 . However, some researchers have discovered that excessively high planting densities negatively affect the light and temperature environment within the maize canopy 14,15,20 . Improved maize yields are best achieved by instituting reasonable planting densities and improving yield per plant 11,35 . Climate conditions in different maize planting areas in China differ significantly and influence optimal maize planting densities in these regions. Maize planting densities in northwestern China are significantly higher than in other regions, which is consistent with the distribution of solar radiation in China 14,36 . In this study, maize plots were subjected to different degrees of shading that simulated different degrees of solar radiation in different maize planting regions. Our results showed that for different genotypes and different levels of solar radiation during grain filling, the planting densities that achieved the highest yields were different. In areas where solar radiation during grain filling is less than 615 MJ m −2 for XY335 and 630 MJ m −2 for ZD958, increasing planting density does not increase yields. When solar radiation during grain filling was less than 486.9 MJ m −2 for XY335 and 510.9 MJ m −2 for ZD958, yields followed the order of D1 > D2 > D3 and the best planting density was 75,000 plants ha −1 , with average yields of 13.36 Mg ha −1 and 11.09 Mg ha −1 , respectively; when solar radiation during grain filling was higher than 549.5 MJ m −2 and 605.8 MJ m −2 , respectively, yields decreased in the order D3 > D2 > D1 and the appropriate planting density was 120,000 plants ha −1 , with average yields of 20.58 Mg ha −1 and 19.65 Mg ha −1 . With sufficient solar radiation, maize yields at the D3 planting density continued to increase. In contrast, yields at the D1 and D2 planting densities decreased as solar radiation increased. Several reasons could explain this phenomenon. On the one hand, when solar radiation is not a limiting factor, a high planting density can improve light interception; while at low planting densities, excessive solar radiation can result in photoinhibition and decreasing yields 37 . On the other hand, under conditions of limited solar radiation, assimilation at the D1 and D2 planting densities primarily achieves grain filling, whereas assimilation at the D3 density primarily achieves plant maintenance. Therefore, within a narrow range, low planting densities can produce higher yields than high planting densities under low solar radiation conditions 14,36 .
Previous studies have shown that different hybrids and their parents perform differently under weak-light stress. However, most weak-light resistant hybrids tend to be more tolerant of high planting densities 25,[38][39][40][41][42] . Some researchers have found that flat-type hybrids are more sensitive to shading than compact hybrids 20,43 . The current study showed that different maize hybrids responded differently to different levels of solar radiation. Under the same solar radiation conditions, XY335 yields were higher than ZD958 27 . Yield differences between the two cultivars were mainly influenced by solar radiation and not by growth duration, as evidenced by the fact that ZD958 produced lower yields than XY335 even though ZD958 has a longer growth duration. ZD958's yield under S-30% at low planting density, as well as under S-15% at medium planting density, was significantly lower than the CK condition. However, yield reductions for XY335 under the same conditions were not significant 27 . In the current study, maize yield decreased significantly with shading and the degree of decrease differed by cultivar. Yields of the two hybrids increased as solar radiation increased during grain filling. XY335 utilized weak light more efficiently than ZD958 and showed a greater potential to increase yield under high solar radiation conditions. The findings in this study can be useful not only to maize producers in regions with different solar radiation conditions but also to plant breeders.
In recent years, global warming caused many problems, such as drought, high temperature, and decreased solar radiation etc [17][18][19][20]44,45 . Global total solar radiation has shown a downward trend as global warming increases, and climate change has reduced maize yields. In China, for every 10% decrease in solar radiation, the average reduction in maize yields is 9.1% 46 . When solar radiation is not a constraint, grain yields can be improved by increasing planting densities. This study shows that it is possible to mitigate the effects of reduced solar radiation by selecting low-light tolerant cultivars and moderately reducing planting densities. These findings provide theoretical knowledge for achieving high-yield and high-efficiency cultivation of maize in response to future climate change conditions.  Table 1. After shading, the quality of incident light in the maize canopy has been found to remain unchanged 47,48 . The geographical locations of the sites and their climatic conditions during the two maize growing seasons are listed in Table 2. Meteorological data for the 2016 and 2017 maize growing seasons were obtained from meteorological stations located near the experimental sites.

Materials and Methods
We used scaffolding and shading nets to build temporary shading sheds. Shading nets of different shade strength were designed and fabricated. The amount of solar radiation reaching a distance of 1.5 m underneath www.nature.com/scientificreports www.nature.com/scientificreports/ each shading net was measured with a Sunscan (SUNSCAN, Delta-T, UK) and the nets that achieved shading of 15%, 30%, and 50% of natural light were chosen for the experiments. The shading sheds were fixed in place about 1.5 m above the maize canopy in order to maintain the same microclimate conditions (except for solar radiation) observed in the unshaded portions of the field. The shading treatment period began at the silking stage and lasted until maturity because this is the key growth period that determines maize yield and this period is rainy season and decreases in solar radiation often occur during this period 47,48 . The silking stage was defined as the date when 60% of the ears showed silk emergence. Physiological maturity was defined as the date when the black layer appeared.
Different planting densities can provide a better understanding of the interactions between solar radiation, density, and cultivars. The experimental plots were 11 × 10 m 2 in size and each plot included 22 rows that were 9 m long. A 1 m wide walkway was set up between adjacent plots. Borders around each plot consisted of three additional rows on each of two sides of the plot, plus a 1 m long extension at the beginning and end of every row 50 . Sufficient water was applied to prevent water stress. Base fertilizers were applied at rates of 75, 150 kg ha −1 N as urea, 188, 225 kg ha −1 P 2 O 5 (super phosphate), and 53, 75 kg ha −1 K 2 O (potassium sulfate) prior to sowing. An additional 300 kg ha −1 N was applied at both Yinchuan and Qitai during the growing season to prevent nutrient stress 50 . Weeds, diseases, and pests were controlled in the plots. sampling and measurements. At the physiological maturity stage, an area of 8 m 2 was harvested manually from the center of each plot (central four rows, 2 m long) and grain weight was measured. The grain moisture content was determined using a portable moisture meter (PM8188, Kett Electric Laboratory, Tokyo, Japan). Grain yield was determined at 14% moisture content. Aboveground dry matter at physiological maturity was measured by sampling five successive plants in the center row of each plot. Plant samples were dried at 80 °C to a constant weight. Leaf area was determined after 40 days of shading using a Portable Laser Leaf Area Analyser (CI-203; CID Bio-Science, US). statistical analysis. Statistical calculations were performed using Microsoft Excel (Office 2010). Analysis of variance (ANOVA) was used to test the differences between grain yield, biomass, and leaf area index (LAI) under the different shade conditions, planting densities, cultivars, sites, and years using SPSS ver. 21.0 (SPSS Institute Inc., US). Data analysis steps were as follows: (1) the treatments and the corresponding yield data were input into SPSS data analysis view and (2) "One-Way ANOVA" analysis at a 0.05 level of probability followed by Duncan test method was chosen to compare these differences.