Next Article in Journal
AquaCrop-Simulated Response of Sorghum Biomass and Grain Yield to Biochar Amendment in South Sudan
Previous Article in Journal
Identification of Cabbage Seedling Defects in a Fast Automatic Transplanter Based on the maxIOU Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soybean (Glycine max (L.) Merr.) Yield Reduction due to Late Sowing as a Function of Radiation Interception and Use in a Cool Region of Northern Japan

by
Etsushi Kumagai
1,* and
Tomoki Takahashi
2
1
Agricultural Meteorology Group, Division of Agro-Environmental Research, Tohoku Agricultural Research Center, National Agriculture and Food Research Organization, 4 Akahira, Shimokuriyagawa, Morioka, Iwate 020-0198, Japan
2
Soil and Plant Health Management Group, Division of Lowland Farming Research, Tohoku Agricultural Research Center, National Agriculture and Food Research Organization, 3 Furukomichi, Yotsuya, Daisen, Akita 014-0102, Japan
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(1), 66; https://doi.org/10.3390/agronomy10010066
Submission received: 2 December 2019 / Revised: 26 December 2019 / Accepted: 27 December 2019 / Published: 2 January 2020

Abstract

:
Yields decrease when soybean is sown later than recommended in the cool climate of the Tohoku region of Japan. However, the factors responsible for this decrease are not fully understood. We investigated the effects of late sowing on growth, phenological development, yield, yield components, and radiation interception of three soybean cultivars in two consecutive years and analyzed the relationships of those variables with temperature and soil volumetric moisture content (SMC). Averaged across years and cultivars, yields decreased significantly when plants were sown approximately three weeks late. Yield reductions were partially due to reductions in node number per plant, dry matter production, and capture of cumulative irradiance, resulting from slowed canopy development during vegetative and early reproductive stages. The number of seeds per pod was one of the major determinants of the variation in yield. Owing to the delay in sowing date, the reduction in seeds per pod was likely due to low temperatures during the 20 days after seed filling began. Occasional lower SMC during reproductive stages did not affect yield, yield components, and growth parameters. However, these results were obtained from the two years’ experiments. Therefore, further investigations of the relationship of yield with temperature and SMC under different years and sites are needed.

1. Introduction

Soybean (Glycine max [L.] Merr.) is an important crop for food, protein, and oil. In Japan, soybean is used in traditional foods such as tofu, natto, and miso. However, the soybean yield in Japan is much lower than that in the other major producing countries such as the USA, Brazil, and China, and has stagnated over the past 30 years; the average (1982–2012) soybean yield is only 1.65 t ha−1 and has increased by less than 0.3 t ha−1 since 1982 [1]. Therefore, to increase and stabilize soybean yields, a strategy is needed that includes both practical agronomic management and genetics.
The optimization of sowing dates is the most important and least expensive agronomic practice that affects soybean yield [2]. Several studies show that recent yield increases in the USA were achieved by optimizing the sowing date with the cultivar maturity group (MG) [3,4]. In the Tohoku region of northern Japan, paddy rice monoculture is the dominant crop system, and the optimal transplanting time for paddy rice is mid-May to late May. The currently recommended sowing window for soybeans is late May to early June after the rice transplanting. Recently, the population of farmers in Japan has declined. Consequently, the scale at which farmers must operate has expanded [1]. Therefore, to avoid the increase in competition associated with rice transplanting and soybean sowing for large-scale production, soybean farmers are using a relatively longer period for sowing in the Tohoku region [5]. Reductions in soybean yields with delayed sowing are well documented in the Tohoku region, and the simple explanations for the yield reductions are a shorter period from sowing to maturity and cooler temperatures during the reproductive stage [6].
The effects of sowing date on soybeans depend on the meteorological conditions at the site, particularly the amount and the timing of precipitation [7,8]. Waterlogging is a serious problem in soybean production in converted paddy fields, which sometimes have poor drainage and a high water table in Japan [9]. The rainy season is June to July in the Tohoku region, so waterlogging occurs during the early vegetative stage. Excess soil water during the vegetative stage inhibits root growth of soybean. Plants with poor root systems could not absorb soil water and then seed yield decreased with water deficits in the summer [10]. Recently, at the Tohoku Agricultural Research Center (TARC) in Akita Prefecture, Japan, we used the FAO 56 evapotranspiration model [11] in a long-term continuous performance test and analyzed the correlations between soybean yield and soil volumetric moisture content (SMC) [12]. We found that seed yield over 33 years was significantly negatively correlated with average SMC during August. This result provided evidence that soybean yield decreased in drier years with lower SMC during August.
Yield potential is defined as the maximum yield achieved when a crop is grown in the absence of water and nutrient limitations and of damage by insects and pests [13]. It is parameterized by different efficiencies in the following equation suggested by Monteith [14]:
Yield potential = IR × RIE × RUE × HI
where incident radiation (IR) is the total incident solar radiation during the growing season; radiation interception efficiency (RIE) is determined by the rate of development and the duration of canopy coverage; radiation use efficiency (RUE) is determined by the amount of solar radiation that is transformed into biomass; the harvest index (HI) is determined by the amount of biomass allocated to vegetative versus reproductive organs. This equation can provide insight into how practical agronomic management and genetics determine the yield potential.
In the Tohoku region, late sowing may shorten the growth period from sowing to maturity and thereby may decrease the cumulative solar radiation (CumIR) intercepted. Moreover, it may affect RIE by affecting the development of the canopy and leaf expansion, which depend on climate factors such as temperature, precipitation, and soil water content [15,16]. Although RUE is constant for a given crop species in a given environment, large variability in RUE under optimal conditions has been reported in soybean [17]. The HI also tends to be constant and genotype-dependent; however, some studies show that sowing dates affect biomass partitioning into seed yield [18,19]. Despite many previous investigations of the effect of late sowing on soybean yield in the Tohoku region [6], the importance of radiation interception and use in explaining the lower yields associated with late planting has not been determined. In the USA, yield reductions with late plantings are partially explained by reductions in CumIR around flowering and early pod set [20,21].
In this study, we tested two hypotheses: (1) The yield reduction with late sowing date is explained by the reduction in CumIR during reproductive stages in this region; (2) the yield variations among the combinations of years, sowing dates, and cultivars are explained by variations in SMC. To test the hypotheses, we compared the changes in CumIR and RUE from flowering to late reproductive stages and in HI, yield, and yield components among three leading soybean cultivars between normal (early June) and late (late June) sowing in two consecutive years. Furthermore, we analyzed the effects of temperature and SMC on yields, yield components, and growth parameters of the three cultivars. Our results supported the first hypothesis but did not support the second hypothesis.

2. Materials and Methods

The experiments were conducted at the National Agriculture and Food Research Organization (NARO), TARC in Morioka, Iwate Prefecture, Japan (39°44’ N, 141°7’ E), in 2015 and 2016. In each year, different converted paddy fields were used. The soil characteristics of the two fields are summarized in Table S1. The soil in each field is an Andosol. Each field received 60 g m−2 fused phosphate fertilizer and 100 g m−2 magnesium lime a month before sowing in each year. In addition, the fields received 3 g m−2 N, 12.5 g m−2 P (P2O5 equivalent), and 5 g m−2 K (K2O equivalent) in the form of a compound fertilizer one day before sowing each year. Fertilizers were applied and incorporated to a depth of approximately 14 cm by using rotary tillers.
We used three determinate soybean cultivars: Ryuhou, Enrei, and Satonohohoemi. Within the Tohoku region, Ryuhou is widely grown in Akita and Iwate prefectures, and Enrei is widely grown in Yamagata Prefecture [22,23]. Satonohohoemi is a new cultivar developed at NARO TARC that was released as a newly recommended cultivar for Yamagata Prefecture in 2009 [22]. Ryuhou is a mid-season maturing cultivar and Enrei and Satonohohoemi are late-season maturing cultivars according to flowering dates at NARO TARC [22,23]. All three cultivars are equivalent to MG IV. Plants were grown to maturity. Seeds were treated with a combined insecticide and fungicide (CruiserMaxx, Syngenta Co., Tokyo, Japan) at the manufacturer’s recommended dose and then sown by hand at three seeds per hill on 5 June (normal sowing) and 25 June (late sowing) in 2015 and on 6 June (normal sowing) and 27 June (late sowing) in 2016. The planting density was 28.5 seeds m−2 (9.5 hills m−2) with 0.15 m between hills and 0.70 m between rows. This density corresponded to the recommended local planting densities in the Tohoku region (seven to 15 hills m−2). Plants were thinned to one per hill after establishment. A pre-emergence herbicide (Ecotop, Maruwa Biochemical Co., Tokyo, Japan), which contained 1.0% dimethenamide-p and 1.4% linuron, was applied immediately after sowing for weed control at the manufacturer’s recommended dose. A month after sowing, the area between the rows was tilled, and ridging was used to control weeds and lodging. In each year, the experimental design was a split-plot arrangement; the main plot was sowing date, and the subplot was cultivar, with three replicates. A subplot consisted of seven rows, each 4.5 m in length.
We surveyed the phenology of 10 plants in the center row of each plot at intervals of two or three days and recorded the dates of growth stages based on the staging system of Fehr and Caviness [24]: VE, emergence; R2, full flowering; R5, the beginning of seed filling; R7, the beginning of maturity. We determined the dry weight of aboveground parts per m2 by periodic sampling at R2, R5, and 20 days after R5. The aboveground parts of 0.84 m2 of the center four rows were harvested on each date. Their dry weights were measured after oven drying at 80 °C for three days. The differences in the dry weights of aboveground parts (increase in aboveground biomass—ΔAGB) between R2 and R5 and between R5 and R5 + 20 days were calculated. We estimated the intercepted solar radiation by using digital imaging techniques (GACS1; Kimura Ouyou-Kougei Co., Ltd., Saitama, Japan) according to Kumagai [8]: The fractional canopy cover was determined from digital images taken above the canopy at one-week intervals from the early vegetative to late reproductive stages. In each year, total daily solar radiation (MJ m−2 day−1) was measured by a pyranometer (CMP11, Campbell Scientific Inc., Logan, UT, USA) and recorded by a data logger (CR1000, Campbell Scientific) at a weather station near the study field, and daily incident radiation above the canopy was computed as the product of daily solar radiation and the fractional canopy cover. The CumIR from R2 to R5 and from R5 to R5 + 20 days was computed by summing the values of daily incident radiation above the canopy. The mean radiation use efficiency (RUE) for biomass production was determined by dividing the ΔAGB by the CumIR.
At maturity, we manually harvested the aboveground parts from 4.62 m2 of the center four rows of each plot. After the plants were completely air-dried, we removed the few remaining leaves and petioles and then weighed the remaining aboveground parts. The aboveground biomass equaled the sum of the weights of the stems, pod shells, and seeds. The seed yield per m2 was adjusted to 15% moisture content, and the harvest index (HI) was determined as the seed yield divided by the aboveground biomass. The yield components (numbers of nodes, pods, seeds per pod, pods per node, and 100-seed weight) were estimated from five representative plants in each plot. The 100-seed weight was also adjusted to 15% moisture.
In addition to solar radiation, daily means of temperature and precipitation during the growing season were recorded at the nearby weather station. The soil volumetric moisture content (SMC, mm3 mm−3) was monitored with time-domain transmission (TDR) sensors (CS625; Campbell Scientific) at approximately seven cm below the soil surface (the midpoint depth of the plow layer) in a plot with the normal sowing date at the two sites. SMC was measured and recorded every 30 min on a data logger (CR200X; Campbell Scientific).
At the end of June (approximately three weeks after sowing in the plot with the normal sowing date), we sampled the soil of the plow layer in 100 mL sampling tubes. The SMC of the samples at field capacity (θFC, pF = 1.5) was determined by the sand column method (DIK-3521; Daiki Rika Kogyo Co., Ltd., Saitama, Japan) and that at a permanent wilting point (θWP, pF = 4.2) by a dew point potentiometer (WP4C; Meter Group, Inc., Pullman, WA, USA). Because θFC and θWP differed somewhat between the two years (Table S1), the fraction of available soil water (FASW) was calculated as (SMC − θWP)/(θFC − θWP) in each year.
We first conducted ANOVA with a split-split-plot design, where treating year as the main factor, sowing date as the split factor, and cultivar as the split-split factor. Year, sowing date, cultivar, and their interaction were each considered as fixed effects, and replicate (block) was considered as a random effect. When the ANOVA produced a significant result (p < 0.05), Fisher’s LSD test was used to detect significant differences between means. We used multiple regression analysis to examine the relative importance of yield components in the two years of yield data. We used the natural logarithm transformation because the relationships between yield and its components were multiplicative. We also examined the relations between yield and yield components and between temperature and FASW during the reproductive stage by simple correlation (Pearson’s) and partial correlation analyses. All procedures were performed in SPSS v. 23.0 software (IBM, Tokyo, Japan).

3. Results and Discussion

3.1. Meteorological Conditions and Fraction of Available Soil Water

The two years had summer with higher temperatures and abundant solar radiation. The amount and pattern of precipitation during summer differed between years. Overall, 2015 was drier and 2016 was wetter than the long-term averages. In 2015 and 2016, the five-month mean temperature (from June to October) was 0.4 to 0.7 °C higher than the 30-year mean (Table 1). The largest departure of 2.1 °C above the 30-year mean was observed in July 2015 and September 2016. In addition, mean temperatures in August 2016 were 1.2 °C higher than the 30-year mean. The five-month mean daily solar radiation in 2015 and 2016 was 1.6 and 1.3 MJ m−2 day−1 respectively, higher than the 30-year mean. In 2015, the solar radiation in June and July was higher than normal, and in 2016, that in August was 3.1 MJ m−2 day−1 higher than the 30-year mean. The five-month mean precipitation in 2015 was lower than the 30-year mean (121.2 vs. 147.0 mm, respectively), whereas that in 2016 was higher (179.0 mm). In 2015, the monthly total precipitation was lower than the 30-year mean by 87.9 mm in July, 36.4 mm in August, and 51.4 mm in September. In 2016, the monthly total precipitation in June, August, and October was higher than the 30-year mean but that in July was lower.
The changes in the FASW during the two growing seasons were different as shown in Figure 1. In 2015, the FASW decreased to low levels at the end of July (day of the year—DOY, 203) and in the middle of August (DOY 223), which corresponded to the vegetative stages and the early reproductive stage (around R2), respectively of the normal sowing. However, it was high and relatively stable afterward. In 2016, the FASW fluctuated and decreased to relatively low levels several times from late July to early September (DOY 206, 211, 218, 226, and 247). DOY 218 and 226 corresponded to the early reproductive stage (around R2) in normal and late sowing, respectively. DOY 247 corresponded to the late-reproductive stage (after R5) in both sowing treatments. According to Allen et al. [11], soybean transpiration is inhibited when the FASW is depleted below 0.5. Therefore, drought stress might have sometimes limited soybean growth during the two growing seasons.

3.2. Aboveground Biomass, Seed Yield, Harvest Index, and Yield Components

The effects of year, sowing date, and cultivar on aboveground biomass, seed yield, harvest index, and yield components determined by ANOVA are shown in Table 2. The main finding of this study was the negative effects of late sowing date on aboveground biomass, seed yield, HI, number of pods, seeds per pod, and number of nodes as discussed below.
Although the effect of sowing date on aboveground biomass and seed yield was significant (p < 0.01), these variables were not affected by year or cultivar. Aboveground biomass and seed yield were not significantly affected by the interaction of year × sowing date × cultivar. When averaged across years and cultivars, late sowing led to significant decreases in aboveground biomass (23%, p < 0.01) and seed yield (29%) (p < 0.001). Cultivar and sowing date significantly affected the HI (p < 0.01). Averaged across years and sowing dates, Ryuhou had the highest HI, followed by Enrei and Satonohohoemi. Averaged across years and cultivars, late sowing decreased HI significantly from 0.55 to 0.51. However, the decrease in the HI with late sowing differed between years and among cultivars, as indicated by the significant interactions between year and sowing date (p < 0.05), between sowing date and cultivar (p < 0.01), and among year, sowing date, and cultivar (p < 0.001). The largest decrease in HI occurred in Enrei in 2015 and in Satonohohoemi in 2016 (Table S2).
Year significantly affected numbers of pods, seeds per pod, nodes, and pods per node (p < 0.01) but had no effect on 100-seed weight (Table 2). Averaged across sowing dates and cultivars, numbers of pods, nodes, and pods per node were higher and the number of seeds per pod was lower in 2016 than in 2015. Although we detected no difference in the number of pods, the numbers of seeds per pod, nodes, and pods per node and 100-seed weight differed significantly among the cultivars (p < 0.05). The number of seeds per pod was highest in Ryuhou. The number of nodes was higher in Enrei and Satonohohoemi. The number of pods per node was higher in Ryuhou and Satonohohoemi, and the 100-seed weight was highest in Satonohohoemi. Averaged across years and cultivars, late sowing significantly decreased the number of pods by 11% (p < 0.05), the number of nodes by 15% (p < 0.01), and 100-seed weight by 9% (p < 0.01). Late sowing significantly decreased the number of seeds per pod (p < 0.001), for which the interaction year × sowing date × cultivar was significant (p < 0.05). The decrease in seeds per pod differed among the combinations of years and cultivars (Table S2); that of Satonohohoemi in 2016 was the greatest among all combinations. The main effect of sowing date on pods per node was not significant, but the effect on pods per node differed among cultivars, as indicated by a significant interaction between sowing date and cultivar (p < 0.001). The number of pods per node in Enrei and Satonohohoemi decreased with late sowing, whereas that in Ryuhou did not change.
A complex relationship among yield components determined the seed yield. Therefore, we used multiple regression to determine the relative contributions of yield components (numbers of nodes, pods per node, and seeds per node, and 100-seed weight) to seed yield. All four components significantly positively affected yield (Figure 2). The number of seeds per pod contributed the most to yield (standardized multiple regression coefficient β = 0.661, p < 0.001), followed by the number of nodes (β = 0.598, p < 0.001). The value of β was significant (p < 0.01) but low for pods per node and 100-seed weight. This result indicates that the yield reduction due to late sowing was explained primarily by decreases in numbers of seeds per pod and nodes in the Tohoku region. In a review of previous experiments conducted in the Midwestern USA, which has the same latitude and a similar climate as the Tohoku region, soybean yield declined rapidly when seeds were planted after late May [7]. The yield decline following late sowing in north-central USA was due to fewer nodes [25]. De Bruin and Pedersen [26] show that reductions in the number of seeds per m−2 accounted for the yield decline associated with delayed planting. They also show that 100-seed weight is not associated with the yield decline due to late planting. The linear decline in yield per day of planting delay they observed was similar to what we observed (1.7–4.3 g m−2 d−1 vs. 4.2 g m−2 d−1).

3.3. Phenological Development

The other finding of this study was the negative effects of late sowing on the phenology from emergence (VE) to full flowering (R2), from R2 to the beginning of seed filling (R5), from R5 to the beginning of maturity (R7), and from VE to R7 (total growth, Table 3). The main factors of year, sowing date, and cultivar significantly affected the phenology from VE to R2 (p < 0.01), from R2 to R5 (p < 0.01), from R5 to R7 (p < 0.001), and from VE to R7 (p < 0.001), with the exception of the effect of sowing date on R5 to R7. Compared with 2015, the periods from VE to R2, from R5 to R7, and from VE to R7 in 2016 were longer but that from R2 to R5 was shorter. When averaged across years and cultivars, late sowing significantly decreased the time from VE to R2 by 12 days (p < 0.001), from R2 to R5 by one day (p < 0.05), and from VE to R7 by 13 days (p < 0.001). Late sowing affected the vegetative phase (from VE to R2) more strongly than the reproductive phase (from R2 to R5). The cultivar affected the duration of each phase: Satonohohoemi had the longest total (VE to R7), vegetative (VE to R2), and reproductive (R2 to R7) stages when averaged across years and sowing dates. The effect of late sowing on these durations differed substantially among the cultivars, as indicated by significant interactions between sowing date and cultivar (p < 0.01). In both years, the decrease in the duration of total growth (VE to R7) was greatest in Enrei (16 days in 2015, 15 days in 2016; Table S3). Previous studies found that a reduction in yield with late sowing is due to a shortened duration from sowing to flowering and to maturity [6,25], which is consistent with our results.

3.4. Canopy Cover, Dry Matter Production, Radiation Interception, and Use Efficiency

Indices of soybean growth, such as canopy cover [27], crop growth rate (CGR), RIE, and RUE [20], are correlated with seed yield. The changes in percentage canopy cover are shown in Figure S1. The cover was significantly lower in late sowing plots than in normal sowing plots for most of the growth period (DOY 200 to 240) in all cultivars in both years.
ΔAGB, CumIR, and RUE during the reproductive stages (from R2 to R5 and from R5 to R5 + 20 days) are summarized in Table 4. One of the important findings of this study was the reductions of ΔAGB and CumIR from R2 to R5 due to late sowing. When averaged across years, late sowing significantly decreased ΔAGB from R2 to R5 by 30% and CumIR from R2 to R5 by 31% (p < 0.001), and for both parameters, the interaction between sowing date and cultivar in the period was significant (p < 0.05). The decreases in ΔAGB and CumIR were larger for Satonohohoemi and Enrei than for Ryuhou (Table S4). By contrast, from R5 to R5 + 20 days, the effects of year, sowing date, and cultivar on ΔAGB were not significant, nor were the interactions among them. However, the effects of year, sowing date, and cultivar on CumIR were significant, as were all interactions among them (p < 0.01). As shown in Figure S2, the solar radiation after R5 differed between the two years. Therefore, these effects and interactions were due primarily to yearly changes in solar radiation during this period rather than to changes in canopy cover (which is related to RIE), because canopy cover reached approximately 100% at R5 (Figure S1), and there was no difference among all combinations. The effects of year, sowing date, and cultivar on RUE were not significant, nor were the interactions among them from R2 to R5 and from R5 to R5 + 20 days, except for a significant interaction between year and sowing date from R2 to R5 (p < 0.05). Our results indicate that the decrease in ΔAGB from R2 to R5 with late sowing was due to the decrease in CumIR rather than that in RUE. The decrease in CumIR that resulted from late sowing in this study was due to both slowed canopy development and a shortened period from R2 to R5.
Pearson’s correlation coefficients were determined for the relationships between ΔAGB and seed yield and yield components when data were combined across years, sowing dates, and cultivars (Table 5). ΔAGB from R2 to R5 was highly correlated with seed yield (Pearson’s r = 0.639, p < 0.001) and number of nodes (r = 0.606, p < 0.001) and was correlated with 100-seed weight (p < 0.05). However, from R5 to R5 + 20 days, the correlations between ΔAGB and yield components were not significant. Thus, the results support the first hypothesis that the yield reduction with late sowing date is explained by the reduction in CumIR during reproductive stages. Our results reveal that yield reduction with late sowing was due primarily to the reduction in the number of nodes, which in turn resulted from the decrease in the ΔAGB associated with the decrease in CumIR during the early to middle reproductive stage (R2 to R5). Although the appearance of main stem nodes stops soon after R2, branch nodes continue to increase from R1 to R5 in determinate soybean cultivars [28]. Therefore, by stage R5 in determinate soybean cultivars, the total number of nodes is determined. High CumIR, the ΔAGB, and the associated total number of nodes attained by R5 were important factors in preventing yield reduction of soybeans with the late sowing. Although an increase in the sowing rate is widely accepted to prevent yield decline with late sowing [19], the optimal combination of the sowing date with the sowing rate has not been determined. The combined effect of late sowing and dense planting on radiation interception and use in soybean will need to be investigated in future research. Our results are consistent with the well-known premise that the number of pods (seeds) and seed yield depend on assimilate supply during the period before R5 [17,21,29]. However, Shiraiwa et al. [30] found that soybean yield was significantly and positively correlated with CGR (ΔAGB/days) during the 20-day period after R5 in multiple environments (four seasons) in a warm region of Japan (Kyoto and Shiga prefectures), although they did not observe any positive correlations between seed yield and CGR before or after this period that were consistent among the environments. The weather conditions at the experimental sites and the cultivars used are possible explanations for the differences between the two studies. Thus, in future research, more comprehensive investigations of the relationship between yield and dry matter accumulation in many genotypes under different environmental conditions are needed.

3.5. Correlations of Seed Yield, Yield Components, and Growth Parameters with Temperature and Fraction of Available Soil Water

Cool temperatures during the reproductive stage are one of the limiting factors for seed yield in the northern region of Japan [6,31]. Both flooding and drought stress also influence soybean production in Japan [9,10,11]. As shown in Table 1 and Figure 1 and Figure S2, weather conditions and FASW during the reproductive stages varied across years, sowing dates, and cultivars in this study. Although the dataset was limited from two years at one site, we attempted to elucidate how temperature and FASW influenced seed yields, yield components, and growth parameters in the three cultivars sown on two sowing dates and in two years. The values of Pearson’s r were significant and positive for the relationships between mean temperature from R2 to R5 and pods per node (p < 0.05) and between mean temperature from R5 to R5 + 20 days and seed yield (p < 0.05; Table S5), presumably indicating that cooler temperatures during the reproductive stage due to late sowing (Figure S1) reduced seed yield and yield components. Furthermore, significant and negative correlation coefficients were found for the relationships between mean FASW during 20 days after R5 and both node number and pods per node (p < 0.05, Table S5). This result might indicate that excessive soil water reduced yield components. However, significant and negative correlation coefficients were detected for the relationships between mean temperature and mean FASW during the two periods (p < 0.001), indicating that excessive soil water because of the high amount of precipitation was associated with cold weather.
The partial correlation coefficients, when mean FASW was held constant, were determined for the relationships between yield components, and growth parameters and mean temperature during the reproductive stage (Table 6). The partial correlation coefficient between mean temperature and seeds per pod with a fixed mean FASW during the 20 days after R5 was 0.864 (p < 0.001). Although the number of seeds per pod remains stable under a wide range of temperatures [32], Kurosaki and Yumoto [33] report that low temperature (15.5 °C mean temperature) for four weeks in the flowering season reduced the number of seeds per pod in soybean cultivars in Hokkaido (the northernmost prefecture in Japan). In this study, the mean temperature was 22.6 °C during the period from R2 to R5 and 19.3 °C for the 20 days after R5 in the late-sown soybeans. Those temperatures might not negatively affect the number of seeds per pod. Kurosaki and Yumoto [33] also reported that the response of seeds per pod to low temperature was different among soybean cultivars. Therefore, one plausible explanation for the variation in the number of seeds per pod across years, sowing dates, and cultivars was both temperature variability during the initial seed filling stage and the differential responses of cultivars to low temperatures. Further investigation of the effect of low temperature on yield components in late-sown soybean in different years and sites will be required to optimize the combination between sowing dates and cultivars in cool regions of Japan.
The effect of water stress on soybean yield depends not only on the duration and intensity of the stress but also on the timing [34]. Drought stress during the reproductive stage increases the rate of soybean flower and pod abortion [35] and reduces seed size [36]. Frederick et al. [37] report that drought stress during the period from R1 to R7 decreases seed number and yield on branches, because most branch growth occurs between R1 and R5 in determinate soybean cultivars. Brevedan and Egli [38] found a lower photosynthetic rate in water-stressed soybean, which contributed to earlier maturity. However, according to Wani et al. [39], soybean yields are not substantially affected until root zone soil has been depleted below 0.6 FASW and leaf area development has not been hampered by prolonged mild water stress during the vegetative phase. Soltani and Sinclair [14] show that the threshold FASW for leaf area development is 0.31 and that for growth is 0.25. Although FASW occasionally fell below 0.5, at which transpiration in soybean is inhibited [24], around R2 (Figure 1), no positive linear relationship was detected between mean FASW during reproductive stages and yield, yield components, or growth parameters (Table S5), indicating that lower FASW in the plow layer at least did not affect soybean growth and yield. We did not detect the threshold FASW for inhibiting growth and yield in this study. These results did not support our second hypothesis that the yield variations among the combinations of years, sowing dates, and cultivars are explained by variations in SMC. However, it is possible that using the dataset from two years and two sowing dates prevented the present study from detecting the effects of drought stress. Moreover, soil physic characteristic was site-specific within the Tohoku region [40]. Therefore, further investigation under different years (weather) and sites (soil) will be needed to test this hypothesis in the future. Another possible explanation for the absence of a drought stress effect in this study could be that water flowed from the subsoil below the plow layer. Generally, the taproot of soybean cannot penetrate a hardpan (a compacted layer) immediately below the plow layer, and therefore the effective root zone is limited in converted paddy fields in Japan [40]. We found that the taproot length was approximately 14 cm in the two years (data not shown), which was equivalent to the plow depth. Therefore, uptake of water from the subsoil below the hardpan by taproots seems unlikely. However, capillary rise (upward water flow from subsoil) is a feature of Andosols [41], and thus taproots near the hardpan might avoid drought stress. To better understand water flow in soil and the effects on soybean growth, further investigations are also necessary.

4. Conclusions

The results from two years of field experiments in northern Japan show that sowing approximately three weeks late decreased soybean yields. The yield reduction was partially caused by reductions in the number of nodes, ΔAGB, and CumIR due to slower canopy development during vegetative and early reproductive stages. The number of seeds per pod was one of the major determinants of the variation in yield. Low temperatures caused by the delay in sowing date during the 20 days after R5 likely explained the reduction in the number of seeds per pod. Our study also showed that a lower FASW during reproductive stages did not affect seed yield, yield components, or growth parameters.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/10/1/66/s1: Table S1: Soil characteristics of the experimental soybean fields at TARC, NARO, Morioka, Iwate Prefecture, Japan, in 2015 and 2016; Table S2: Aboveground biomass (g m−2), seed yield (g m−2), harvest index (HI), and yield components (numbers of pods per plant, seeds per pod, nodes per plant, and pods per node, and 100-seed weight, g) of three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan; Table S3: Days from emergence (VE) to full flowering (R2), from R2 to the beginning of seed filling (R5), from R5 to the beginning of maturity (R7), and from VE to R7 for three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan; Table S4: Increase in aboveground biomass (ΔAGB, g m−2), cumulative intercepted irradiation (CumIR, MJ m−2), and radiation use efficiency for dry matter production (RUE, g MJ−1) during reproductive stages (from R2 to R5 and from R5 to 20 days after R5) of three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan; Table S5: Pearson’s coefficients of correlation of the relationships between mean temperature and the FASW (fraction of available soil water) during reproductive stages (from R2 to R5 and from R5 to R5 + 20 days) and seed yield (g m−2), yield components (numbers of pods per plant, seeds per pod, nodes per plant, pods per node, and 100-seed weight, g), increase in aboveground biomass (ΔAGB, g m−2), cumulative intercepted irradiation (CumIR, MJ m−2), and radiation use efficiency for dry matter production (RUE, g MJ−1) across three soybean cultivars and two sowing dates in 2015 and 2016 (n = 12) in the Tohoku region of northern Japan; Figure S1: Changes in the canopy cover (%) of three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) growing after normal and late sowing during the (A) 2015 and (B) 2016 seasons in the Tohoku region of northern Japan. DOY, day of the year. Values represent the mean (n = 3); Figure S2. Changes in daily mean temperature and solar radiation from June to October during the soybean growing season in (A and C) 2015 and (B and D) 2016 in the Tohoku region of northern Japan. DOY, day of the year. R2, beginning of flowering; R5, beginning of seed filling stage; | normal sowing date; ⁞ late sowing date; arrows indicate means of the three cultivars.

Author Contributions

Conceptualization, E.K. and T.T.; methodology, E.K. and T.T.; investigation, E.K.; writing—original draft preparation, E.K.; writing—review and editing, T.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Council for Science, Technology, and Innovation (CSTI), the Cross-ministerial Strategic Innovation Promotion Program (SIP), and the “Technologies for creating next-generation agriculture, forestry and fisheries” (funding agency: Bio-oriented Technology Research Advancement Institution, NARO, Japan).

Acknowledgments

The authors are grateful to Satoshi Nakano (Institute for Agro-Environmental Sciences, NARO), Naoki Matsuo (Kyushu Okinawa Agricultural Research Center, NARO), and Toshihiro Hasegawa (TARC, NARO) for their advice on conducting these experiments and analyzing the results. We thank Hisashi Tamura, Eisaku Kumagai, Fumihiko Saitou, Akihiko Umihata, and Kafumi Segawa for their technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. MAFF. The Statistical Yearbook of Ministry of Agriculture, Forestry and Fisheries; Ministry of Agriculture, Forestry and Fisheries: Tokyo, Japan, 2018. Available online: http://ww.maff.go.jp/j/tokei/kouhyou/sakumotu/sakkyou_kome/index.html (accessed on 2 December 2019).
  2. Robinson, A.P.; Conley, S.P.; Volenec, J.J.; Santini, J.B. Analysis of high yielding, early-planted soybean in Indiana. Agron. J. 2009, 101, 131–139. [Google Scholar] [CrossRef]
  3. Chen, G.; Wiatrak, P. Soybean development and yield are influenced by planting date and environmental conditions in the southeastern coastal plain, United States. Agron. J. 2010, 102, 1731–1737. [Google Scholar] [CrossRef]
  4. Salmerόn, M.; Gbur, E.E.; Bourland, F.M.; Golden, B.R.; Earnest, L.; Purcell, L.C. Soybean maturity group choices for maximizing light interception across planting dates in the U.S. Midsouth. Agron. J. 2015, 107, 2132–2142. [Google Scholar] [CrossRef]
  5. Mochida, H. Effects of sowing time and plating pattern on soybean yield and seed quality in locations along the Sea of Japan at Tohoku region. Bull. Tohoku Agric. Res. Center 2016, 118, 69–77. [Google Scholar]
  6. Tada, H.; Anamizu, K.; Kimura, S. The growth response of soybean plant to meteorological factors in Aomori prefecture. I. The effect of seeding date on the growth and yield. Tohoku J. Crop Sci. 1984, 27, 96–97. [Google Scholar]
  7. Egli, D.B.; Cornelius, P.L. A regional analysis of the response of soybean yield to planting date. Agron. J. 2009, 101, 330–335. [Google Scholar] [CrossRef]
  8. Kumagai, E. Effect of early sowing on growth and yield of determinate and indeterminate soybean (Glycine max (L.) Merr.) cultivars in a cool region of northern Japan. J. Agric. Meteorol. 2018, 74, 18–28. [Google Scholar] [CrossRef] [Green Version]
  9. Bajgain, R.; Kawasaki, Y.; Akamatsu, Y.; Tanaka, Y.; Kawamura, H.; Katsura, K.; Shiraiwa, T. Biomass production and yield of soybean grown under converted paddy fields with excess water during the early growth stage. Field Crops Res. 2015, 180, 221–227. [Google Scholar] [CrossRef]
  10. Hirasawa, T.; Nakahara, M.; Izumi, T.; Iwamoto, Y.; Ishihara, K. Effects of pre-flowering soil moisture deficits on dry matter production and ecophysiological characteristics in soybean plants under well irrigated conditions during grain filling. Plant Prod. Sci. 1998, 1, 8–17. [Google Scholar] [CrossRef]
  11. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines For Computing Crop Water Requirements; FAO Irrigation and Drainage Paper No 56; Food and Agriculture Organization, Land and Water: Rome, Italy, 1998. [Google Scholar]
  12. Kumagai, E.; Takahashi, T.; Nakano, S.; Matsuo, N. Relationship between soybean yield and drought in long-term continuous performance test at Tohoku agricultural research center, NARO: Analysis by using agrometeorological grid square data and FAO56 evapotranspiration model. Jpn. J. Crop Sci. 2018, 87, 233–241. [Google Scholar] [CrossRef] [Green Version]
  13. Evans, L.T.; Fischer, R.A. Yield potential: Its definition, measurement, and significance. Crop Sci. 1999, 39, 1544–1551. [Google Scholar] [CrossRef]
  14. Monteith, J.L. Climate and the efficiency of crop production in Britain. Phil. Trans. R. Soc. B 1997, 281, 277–294. [Google Scholar] [CrossRef]
  15. Soltani, A.; Sinclair, T.R. Phenology—Temperature and photoperiod. In Modeling Physiology of Crop Developments, Growth and Yield; Soltani, A., Sinclair, T.R., Eds.; CABI: Cambridge, MA, USA, 2012; pp. 73–87. [Google Scholar]
  16. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. AquaCrop—The FAO Crop Model to Simulate Yield Response to Water: I. Concepts and Underlying Principles. Agron. J. 2009, 101, 426–437. [Google Scholar] [CrossRef] [Green Version]
  17. Van Roekel, R.J.; Purcell, L.C.; Salmerón, M. Physiological and management factors contributing to soybean potential yield. Field Crops Res. 2015, 182, 86–97. [Google Scholar] [CrossRef]
  18. Calviño, P.A.; Sadras, V.O.; Andrade, F.H. Development, growth and yield of late-sown soybean in the southern Pampas. Eur. J. Agron. 2003, 19, 265–275. [Google Scholar] [CrossRef]
  19. Kawasaki, Y.; Yamazaki, R.; Katayama, K. Effects of late sowing on soybean yields and yield components in southwestern Japan. Plant Prod. Sci. 2018, 21, 339–348. [Google Scholar] [CrossRef] [Green Version]
  20. Egli, D.B.; Guffy, R.D.; Heitholt, J.J. Factors associated with reduced yields of delayed planings of soybean. J. Agron. Crop Sci. 1987, 159, 176–185. [Google Scholar] [CrossRef]
  21. Egli, D.B.; Bruening, W.P. Potential of early-maturing soybean cultivars in late plantings. Agron. J. 2000, 92, 532–537. [Google Scholar] [CrossRef]
  22. Kikuchi, A.; Kono, Y.; Kato, S.; Yumoto, S.; Takada, Y.; Shimada, S.; Sakai, T.; Shimada, H.; Takahashi, K.; Adachi, T.; et al. A new soybean cultivar, “Satonohohoemi”, with high protein content, lodging resistance, large seed size and good quality. Bull. Natl. Agric. Res. Cent. Tohoku Reg. 2011, 113, 1–15. [Google Scholar]
  23. Sato, Y.; Inoue, K.; Suzuki, M.; Igarashi, H.; Fujimoto, J.; Okada, K. Studies on the new recommended soybean variety ‘RYUHOU’. Res. Rep. Akita Agric. Res. Cent. 1998, 38, 80–93. [Google Scholar]
  24. Fehr, W.R.; Caviness, C.E. Stages of Soybean Development; Iowa Agriculture and Home Economics Experiment Station. Cooperative Extension Service; Iowa State University Digital Press: Ames, IA, USA, 1977; p. 80. [Google Scholar]
  25. Bastidas, A.M.; Setiyono, T.D.; Dobermann, A.; Cassman, K.G.; Elmore, R.W.; Graef, G.L.; Specht, J.E. Soybean Sowing Date: The Vegetative, Reproductive, and Agronomic Impacts. Crop Sci. 2008, 48, 727–740. [Google Scholar] [CrossRef] [Green Version]
  26. De Bruin, J.L.; Pedersen, P. Soybean seed yield response to planting date and seeding rate in the upper Midwest. Agron. J. 2008, 100, 696–703. [Google Scholar] [CrossRef]
  27. Steele, C.C.; Grabau, L.J. Planting dates for early-maturing soybean cultivars. Agron. J. 1997, 89, 449–453. [Google Scholar] [CrossRef]
  28. Board, J.E.; Harville, B.G. Soybean yield component responses to a light interception gradient during the reproductive period. Crop Sci. 1993, 33, 772–777. [Google Scholar] [CrossRef]
  29. Jiang, H.; Egli, D.B. Soybean seed number and crop growth rate during flowering. Agron. J. 1995, 87, 264–267. [Google Scholar] [CrossRef]
  30. Shiraiwa, T.; Ueno, N.; Shimada, S.; Horie, T. Correlation between yielding ability and dry matter productivity during initial seed filling stage in various soybean genotypes. Plant Prod. Sci. 2004, 7, 138–142. [Google Scholar] [CrossRef]
  31. Ohnishi, S.; Miyoshi, T.; Shirai, S. Low temperature stress at different flower developmental stages affects pollen development, pollination, and pod set in soybean. Environ. Exp. Bot. 2010, 69, 56–62. [Google Scholar] [CrossRef]
  32. Kumagai, E.; Yamada, T.; Hasegawa, T. Is the yield change due to warming affected by photoperiod sensitivity? Effects of the soybean E4 locus. Food Energy Secur. 2019, e186. [Google Scholar] [CrossRef]
  33. Kurosaki, H.; Yumoto, S. Effects of low temperature and shading during flowering on the yield components in soybeans. Plant Prod. Sci. 2006, 6, 17–23. [Google Scholar] [CrossRef]
  34. Desclaux, D.; Roumet, P. Impact of drought stress on the phenology of two soybean (Glycine max L. Merr.) cultivars. Field Crops Res. 1996, 46, 61–70. [Google Scholar] [CrossRef]
  35. Westgate, M.E.; Schussler, J.R.; Reicosky, D.C.; Brenner, M.L. Effect of water deficits on seed development in soybean. 2. Conservation of seed growth-rate. Plant Physiol. 1989, 91, 980–985. [Google Scholar] [CrossRef] [PubMed]
  36. DeSouza, P.I.; Egli, D.B.; Bruening, W.P. Water stress during seed filling and leaf senescence in soybean. Agron. J. 1997, 89, 807–812. [Google Scholar] [CrossRef]
  37. Frederick, J.R.; Camp, C.R.; Bauer, P.J. Drought-stress effects on branch and mainstem seed yield and yield components of determinate soybean. Crop Sci. 2001, 41, 759–763. [Google Scholar] [CrossRef]
  38. Brevedan, R.E.; Egli, D.B. Short periods of water stress during seed filling, leaf senescence, and yield of soybean. Crop Sci. 2003, 43, 2083–2088. [Google Scholar] [CrossRef]
  39. Wani, S.P.; Heng, L.K. Crop Yield Response to Water. FAO Irrigation and Drainage Paper No 66; Food and Agriculture Organization of the United Nation, Land and Water: Rome, Italy, 2013; pp. 24–131. [Google Scholar]
  40. Takahashi, T.; Mochida, H.; Sakakibara, M.; Morimoto, S.; Kobayashi, H.; Aiba, S. Investigation of the Factors Reducing Soybean Productivity in the Tohoku Region of Japan. Bull. Tohoku Agric. Res. Cent. 2014, 116, 89–118. [Google Scholar]
  41. Takahashi, T. Soil Moisture Model for Irrigation of Soybean. 2020; in preparation. [Google Scholar]
Figure 1. Changes in the fraction of available soil water (FASW) from July to September during the soybean growing season in (A) 2015 and (B) 2016 in the Tohoku region of northern Japan. DOY, day of the year. R2, beginning of flowering; R5, beginning of seed filling stage; | normal sowing date; ⁞ late sowing date; arrows indicate means of the three cultivars.
Figure 1. Changes in the fraction of available soil water (FASW) from July to September during the soybean growing season in (A) 2015 and (B) 2016 in the Tohoku region of northern Japan. DOY, day of the year. R2, beginning of flowering; R5, beginning of seed filling stage; | normal sowing date; ⁞ late sowing date; arrows indicate means of the three cultivars.
Agronomy 10 00066 g001
Figure 2. Diagram showing the regression coefficients of yield components in the multiple regression on seed yield and the coefficients of correlations between variables obtained across two years and three cultivars in the Tohoku region of northern Japan. Because the nature of the relationship between yield and yield components was multiplicative, the data of all cultivars in each block (n = 36) were natural-log transformed. Values with the single-headed arrows are the standardized multiple regression coefficients, and those with double-headed arrows are Pearson’s correlation coefficients. *** p < 0.001; ** p < 0.01; ns, not significant.
Figure 2. Diagram showing the regression coefficients of yield components in the multiple regression on seed yield and the coefficients of correlations between variables obtained across two years and three cultivars in the Tohoku region of northern Japan. Because the nature of the relationship between yield and yield components was multiplicative, the data of all cultivars in each block (n = 36) were natural-log transformed. Values with the single-headed arrows are the standardized multiple regression coefficients, and those with double-headed arrows are Pearson’s correlation coefficients. *** p < 0.001; ** p < 0.01; ns, not significant.
Agronomy 10 00066 g002
Table 1. Monthly and five-month (June–October) mean air temperature (°C), solar radiation (MJ m2 day1), and total precipitation (mm) in soybean fields during the 2015 and 2016 growing seasons compared with the 30-year means in the Tohoku region of northern Japan.
Table 1. Monthly and five-month (June–October) mean air temperature (°C), solar radiation (MJ m2 day1), and total precipitation (mm) in soybean fields during the 2015 and 2016 growing seasons compared with the 30-year means in the Tohoku region of northern Japan.
JuneJulyAugustSeptemberOctoberFive-Month Mean
Temperature (°C)201518.523.422.717.910.818.7
201618.121.624.120.311.219.0
30-year mean17.821.322.918.211.518.3
Solar radiation (MJ m−2 day−1)201519.217.815.613.211.115.4
201616.416.918.212.511.615.1
30-year mean17.115.115.111.99.913.8
Precipitation (mm)A2015123.5108.5145.5104.5124.0121.2
2016160.0131.0258.5160.5185.0179.0
30-year mean108.1196.4181.9155.992.5147.0
Table 2. Aboveground biomass (g m2), seed yield (g m2), harvest index (HI), and yield components (numbers of pods per plant, seeds per pod, nodes per plant, and pods per node, and 100-seed weight, g) of three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan.
Table 2. Aboveground biomass (g m2), seed yield (g m2), harvest index (HI), and yield components (numbers of pods per plant, seeds per pod, nodes per plant, and pods per node, and 100-seed weight, g) of three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan.
Above-Ground Biomass (g m−2)Seed Yield (g m−2)HINumber of Pods (Per Plant)Seeds Per PodNumber of Nodes (Per Plant)Number of Pods Per Node100-Seed Weight (g)
Year (Y)20154472370.5353.41.4240.01.3431.0
20165122770.5476.91.2349.21.5730.8
Sowing date (S)Normal5423000.5569.01.4448.31.4332.3
Late4162140.5161.31.2040.91.4829.4
Cultivar (Cv)Ryuhou4592740.59a64.91.51a41.7b1.54a30.0b
Enrei4602460.52b61.91.41b45.6a1.34b28.6b
Satonohohoemi5192520.47c68.61.04c46.3a1.48a34.0a
ANOVAYnsnsns********ns
S*************ns**
Cvnsns***ns**********
Y × Snsns*ns**nsnsns
Y × Cvnsns***nsnsns**ns
S × Cvnsns**nsnsns***ns
Y × S × Cvnsns***ns*nsnsns
Values of a parameter followed by the same letter did not differ significantly between cultivars (ANOVA followed by Fisher’s LSD test, p < 0.05). ANOVA results: *** p < 0.001; ** p < 0.01; * p < 0.05; ns, not significant.
Table 3. Days from emergence (VE) to full flowering (R2), from R2 to the beginning of seed filling (R5), from R5 to the beginning of maturity (R7), and from VE to R7 for three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan.
Table 3. Days from emergence (VE) to full flowering (R2), from R2 to the beginning of seed filling (R5), from R5 to the beginning of maturity (R7), and from VE to R7 for three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan.
Days from VE to R2Days from R2 to R5Days from R5 to R7Days from VE to R7
Year (Y)2015462046112
2016481950116
Sowing date (S)Normal532048121
Late411948108
Cultivar (Cv)Ryuhou45c18b45b107c
Enrei47b20a45b113b
Satonohohoemi49a20a54a123a
ANOVAY**********
S****ns***
Cv************
Y × S*nsns**
Y × Cvns**ns**
S × Cv**********
Y × S × Cv**ns*****
Values of a parameter followed by the same letter did not differ significantly between cultivars (ANOVA followed by Fisher’s LSD test, p < 0.05). ANOVA results: *** p < 0.001; ** p < 0.01; * p < 0.05; ns, not significant.
Table 4. Increase in aboveground biomass (ΔAGB, g m2), cumulative intercepted irradiation (CumIR, MJ m2), and radiation use efficiency for dry matter production (RUE, g MJ1) during the reproductive stages (from R2 to R5 and from R5 to 20 days after R5) of three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan.
Table 4. Increase in aboveground biomass (ΔAGB, g m2), cumulative intercepted irradiation (CumIR, MJ m2), and radiation use efficiency for dry matter production (RUE, g MJ1) during the reproductive stages (from R2 to R5 and from R5 to 20 days after R5) of three soybean cultivars (Ryuhou, Enrei, Satonohohoemi) with normal and late sowing dates in 2015 and 2016 in the Tohoku region of northern Japan.
ΔAGBCumIRRUE
(g m2)(MJ m2)(g MJ1)
R2–R5R5–R5 + 20 dR2–R5R5–R5 + 20 dR2–R5R5–R5 + 20 d
Year (Y)20152172122552030.861.06
20162312712562730.890.99
Sowing date (S)Normal2642473032340.871.08
Late1842352072420.890.97
Cultivar (Cv)Ryuhou191b252235b247a0.821.03
Enrei241a239269a232b0.911.06
Satonohohoemi240a233261a233b0.900.99
ANOVAYnsnsns***nsns
S***ns*****nsns
Cv**ns******nsns
Y × Snsns******ns
Y × Cvnsnsns***nsns
S × Cv*ns******nsns
Y × S × Cvnsnsns***nsns
Values of a parameter followed by the same letter did not differ significantly between cultivars (ANOVA followed by Fisher’s LSD test, p < 0.05). ANOVA results: *** p < 0.001; ** p < 0.01; * p < 0.05; ns, not significant.
Table 5. Pearson’s coefficients of correlation of the relationships between the increase in aboveground biomass (ΔAGB, g m−2) and seed yield (g m−2) and yield components (numbers of nodes per plant, pods per node, and seeds per pod, and 100-seed weight, g) during the reproductive stages (from R2 to R5 and from R5 to 20 days after R5) across three soybean cultivars and normal and late sowing dates in 2015 and 2016 (n = 36) in the Tohoku region of northern Japan.
Table 5. Pearson’s coefficients of correlation of the relationships between the increase in aboveground biomass (ΔAGB, g m−2) and seed yield (g m−2) and yield components (numbers of nodes per plant, pods per node, and seeds per pod, and 100-seed weight, g) during the reproductive stages (from R2 to R5 and from R5 to 20 days after R5) across three soybean cultivars and normal and late sowing dates in 2015 and 2016 (n = 36) in the Tohoku region of northern Japan.
Seed YieldNumber of NodesPods Per NodeSeeds Per Pod100-Seed Weight
ΔAGB R2–R50.6390.606−0.2760.2820.415
Significance******nsns*
ΔAGB R5–R5 + 20 d0.1740.2390.214−0.0350.043
Significancensnsnsnsns
*** p < 0.001; * p < 0.05; ns, not significant.
Table 6. Partial correlation coefficients, when mean fraction of available soil water is held constant, for the relationships between mean temperature during the reproductive stage (from R2 to R5 and from R5 to R5 + 20 days) and yield components (numbers of nodes, pods per node, and seeds per pod, and 100-seed weight), change in aboveground biomass (ΔAGB), cumulative intercepted irradiation (CumIR), and radiation use efficiency for dry matter production (RUE) across three soybean cultivars and two sowing dates in 2015 and 2016 (n = 12) in the Tohoku region of northern Japan.
Table 6. Partial correlation coefficients, when mean fraction of available soil water is held constant, for the relationships between mean temperature during the reproductive stage (from R2 to R5 and from R5 to R5 + 20 days) and yield components (numbers of nodes, pods per node, and seeds per pod, and 100-seed weight), change in aboveground biomass (ΔAGB), cumulative intercepted irradiation (CumIR), and radiation use efficiency for dry matter production (RUE) across three soybean cultivars and two sowing dates in 2015 and 2016 (n = 12) in the Tohoku region of northern Japan.
Number of NodesPods Per NodeSeeds Per Pod100-Seed WeightΔAGBRUECumIR
Temperature R2–R50.5440.557−0.5300.2070.1080.259−0.094
Significancensnsnsnsnsnsns
Temperature R5–R5 + 20 d−0.401−0.3290.864−0.2900.101−0.1330.226
Significancensns***nsnsnsns
*** p < 0.001; ns, not significant.

Share and Cite

MDPI and ACS Style

Kumagai, E.; Takahashi, T. Soybean (Glycine max (L.) Merr.) Yield Reduction due to Late Sowing as a Function of Radiation Interception and Use in a Cool Region of Northern Japan. Agronomy 2020, 10, 66. https://doi.org/10.3390/agronomy10010066

AMA Style

Kumagai E, Takahashi T. Soybean (Glycine max (L.) Merr.) Yield Reduction due to Late Sowing as a Function of Radiation Interception and Use in a Cool Region of Northern Japan. Agronomy. 2020; 10(1):66. https://doi.org/10.3390/agronomy10010066

Chicago/Turabian Style

Kumagai, Etsushi, and Tomoki Takahashi. 2020. "Soybean (Glycine max (L.) Merr.) Yield Reduction due to Late Sowing as a Function of Radiation Interception and Use in a Cool Region of Northern Japan" Agronomy 10, no. 1: 66. https://doi.org/10.3390/agronomy10010066

APA Style

Kumagai, E., & Takahashi, T. (2020). Soybean (Glycine max (L.) Merr.) Yield Reduction due to Late Sowing as a Function of Radiation Interception and Use in a Cool Region of Northern Japan. Agronomy, 10(1), 66. https://doi.org/10.3390/agronomy10010066

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