Growth parameters
Plant height that reflects growth condition among vegetative traits was significantly different under various conditions (Shadakshari et al., 2014) and these differences persisted even in the later stages, showing the order T0 > T1 > T2 > T3. Aboveground biomass was higher in T1 and T2 than in T0 in the early period (13 DAT) and even in the late period (42 DAT), biomass in T1 was greater than that of T0 (Table 2). Despite greater plant height in T0 compared to T1, the shoot weight was lower in T0, which was similar to reduced pod weight. The lower pod number observed in T2 and T3 is due to the occurrence of pod abortion after 30 DAT. This result is similar to the study by Torres et al., (2023), which found that end-trimming treatments resulted in reduced plant height, with the height difference not recovered thereafter. End-trimming treatments reduced plant height by 6-27%, implying a reduced risk of lodging; indeed, the incidence of lodging was lower in end-trimmed plants under field conditions (data not shown).
Pod weight per plant in T
1 was consistently higher than in T
0, whereas T
3 showed the lowest values. In the case of T
1, slight trimming removed leaves that previously functioned as sinks, allowing assimilates to be used for making pods. In contrast, the heavy trimming in T
3 led to assimilate translocation toward new leaf formation rather than pod filling, resulting in lower pod weight. As shown in Table 3, the pod weight showed significant differences between years, possibly due to variation in meteorological conditions (Fig 1). Seed size distribution based on 100-seed weight (Fig 2) revealed that T
1 produced seeds with greater size and less variability.
Mild trimming treatments promoted aboveground biomass accumulation. Aboveground biomass was higher in T
1 and T
2 than in T
0 during the early period and remained higher in T
1 than in T
0 in the late period. This result suggests that mild trimming treatments may promote aboveground biomass accumulation.
Seed weight is highly sensitive to leaf or pod removal (
McAlister and Krober, 1958). In our study, T
1 had a higher number and weight of seeds. T
1 had a higher pod number but fewer seeds than T
0, which resulted in a lower overall yield, but the 100-seed weight was greater. While pod removal generally increases seed size by enhancing assimilate translocation to the remaining pods (
McAlister and Krober, 1958), this is not always the case (
Munier-Jolain et al., 1998). Seed number is primarily influenced by canopy photosynthesis and assimilate availability during the early reproductive stages (R
1-R
3). In contrast, seed size is affected by assimilate supply during the seed growth stage (R
5) (
Egli and Bruening, 2003). The end-trimmed plants produced larger and more uniform seeds than the control (Fig 2). This result implies that end-trimming treatments promoted assimilate translocation.
In Korea, soybean yields exhibit significant year-to-year variation, largely driven by substantial annual fluctuations in weather conditions during the vegetative growth stage following flowering. Notably, significant differences were observed in aerial biomass rather than in seed number and 100-seed weight. In this condition, the end-trimming proved effective in preventing yield loss by mitigating the inhibition of reproductive growth caused by excessive vegetative growth.
The changes in canopy structure occurred due to altered leaf distribution on the main stem in response to trimming intensities. At 30 DAT, the control exhibited the most extended petiole lengths across all layers (Fig 3). When comparing treatments, petiole lengths in the middle layer were the longest across all treatments at 24 DAT. The petiole lengths of the upper layer increased, resulting in both upper and middle layers being longer than the lower layer. In contrast, T
3 did not show a notable increase in petiole length compared to other treatments and the middle layer had longer petioles than the upper layer. Petiole length in the upper layer increases to facilitate morphological adaptation and overcome light competition by accumulating leaves in the middle and upper layers. At 30 DAT, the total number of leaves per plant followed the order: T
1 > T
3 > T
0 > T
2, whereas leaf weight followed the order: T
1 > T
0 > T
2 > T
3 (Fig 3). The difference in these results is due to the thinner, smaller leaves of T
3, as indicated by the low leaf weight. This study showed that light reaching the lowest leaves was higher in treatments T
2 and T
3 (Fig 4).
Petiole length can affect canopy architecture and is a key factor in assimilate transport efficiency
(Zhu et al., 2010). The elongation of petioles in the upper canopy is a plant’s morphological adaptation to overcome light competition from the accumulation of leaves in the middle and upper layers. While 80% of assimilates are produced in the top 20% of the canopy (
Hatfield and Carlson, 1978), most pods are located in the lower two-thirds of the plant
(Hansen et al., 1981) and longer petioles therefore increase the distance for assimilate translocation (
Gent, 1982). The changes in petiole length in this research may explain the lower pod number in T
0 (Table 2 and 3), as the elongation of petioles in the middle layer could have hindered assimilate transport to the pods located there.
Excessive leaf area can limit light penetration and negatively affect yield
(Liu et al., 2008) by lowering the leaf photosynthesis
(Gholamhoseini et al., 2017). Conversely, a smaller leaf area can improve light penetration and air circulation in densely planted crops, positively impacting yield
(Jin et al., 2009). Our findings are consistent with previous reports that a dense leaf distribution in the upper canopy can severely limit light penetration (
Ogbuehi and Brandle, 1982). In our study, the light reaching the lowest leaves was higher in the T
2 and T
3 treatments (Fig 4). At 44 DAT, the leaf area index (LAI) for T
1, where end-trimming intensity was low, was similar to the control, whereas T3, subjected to higher end-trimming intensity, exhibited the lowest LAI. Because reductions in petiole length could modify plant architecture and improve canopy structure, positively influencing yield (
Egli, 1988) and plant architecture influences the soybean yield
(Zhu et al., 2024), the differences in PAR levels within the plant canopy were due to the altered canopy architecture resulting from the changed petiole length. In this study, T
3 maintained elevated chlorophyll levels in the lower leaves throughout the later growth stages.
Photosynthesis-related characteristics
At 44 DAT, leaf area index for T
1, where end-trimming intensity was low, was similar to the control, whereas T
3, subjected to higher end-trimming intensity, exhibited the lowest LAI. T
3 had shorter petioles and smaller newly developed leaves than the untreated control, resulting in a lower LAI. Chlorophyll content at 13 DAT showed an initial decrease with increasing end-trimming intensity, but at later stages, T
3 showed the highest chlorophyll content (Fig 5). When analyzing chlorophyll level by leaf position, treatments T
1 and T
2 showed similar or lower values in the upper leaves than the control, but higher values in the middle layer. T
3 maintained elevated chlorophyll levels in the lower leaves throughout the later growth stages.
We observed that stomatal conductance in trimmed plants remained higher than in the control. Because leaf color influences light distribution within the photosynthetic tissues, improving photosynthetic efficiency
(Adegoye et al., 2023), newly emerged leaves after trimming may exhibit strong photosynthetic activity. The plants that have new leaves emerging after trimming may grow vigorously. Overall, stomatal conductance was greater in end-trimmed plants, with T
3 showing uniformly high stomatal conductance across upper, middle and lower layers. At 30 DAT, photosynthetic rates in the upper layer increased with higher end-trimming intensity. T
3 showed the highest rate in the middle layer, while T
1 showed the lowest.
Strong end-trimming (T
3) increased photosynthesis in the upper layer and mild trimming (T
1) showed a smaller increase in photosynthetic rate (Fig 6). These results suggest that defoliation reduces resource competition between leaves, enhances the efficiency of photosynthesis in the remaining leaves and improves light penetration into the inside of the canopy, thereby accelerating the photosynthetic activity of individual plants. Generally, soybean leaves absorb 85-90% of incoming light, with the upper 25% of the canopy accounting for about 75% of total light interception (
Campbell and Norman, 1998). However, the uppermost leaves often reach light saturation, causing the excess irradiance to be dissipated as heat and reducing light-use efficiency. In contrast, the lower canopy receives less than half the light saturation point, resulting in a structurally inefficient light distribution, excess light at the upper layer and deficiency at the lower layer. Reducing chlorophyll content in the upper leaves limits their light absorption, allowing more light to penetrate the lower layer and resulting in a more balanced irradiation on the total leaves. Consistent with these findings, this study also observed that treatments T
1 and T
2 had higher chlorophyll content in the middle layer than in the upper layer (Fig 5). It was suggested that petiole length in the upper layer increases to facilitate morphological adaptation and overcome light competition by accumulating leaves in the middle and upper layers. SPAD values, stomatal conductance and photosynthetic rates all showed similar trends in response to end-trimming; however, stomatal conductance was the trait that particularly well-reflected the treatment effects. Stomatal conductance clearly demonstrated the effects of end-trimming at early after-treatments and, when combined with LAI, enabled verification of treatment efficiency.
Hyperspectral reflectances and vegetation indices
Hyperspectral reflectance showed noticeable differences among treatments at wavelengths above 700 nm during the early period after end-trimming (13 DAT). In contrast, no significant differences occurred during the later period (42 DAT) (Fig 7). Based on hyperspectral reflectance, there were notable differences among treatments in the early period after end-trimming (13 DAT), whereas by the later period (42 DAT), the differences between treatments weakened. T
2 exhibited the lowest vegetation index values at 13 DAT but showed recovered reflectance to intermediate levels by 42 DAT. In the early period, the control generally showed the highest index values, while vegetation indices tended to decrease with greater end-trimming intensity. However, as the plants grew after end-trimming from 13 to 42 DAT, T
3 recorded the highest values (Table 4). Vegetation indices such as Vog REI1, PSSRa, PSSRb, GNDVI and GCI were lowest in T
2 and T
1 at 13 and 42 DAT, respectively and maintained higher values in T
3 throughout the later period. PSRI (Plant Senescence Reflectance Index) produced negative values close to zero when leaf color was darker and shifted to positive values with senescence. PRI increased with end-trimming intensity at 13 DAT, then decreased at 42 DAT. Across treatments, the control consistently showed the highest value in most vegetation indices. T
1 showed vigorous initial growth at 13 DAT, followed by the occurrence of senescence symptoms. T
2 exhibited initial stress by mechanical trimming but recovered in the later period. T
3 showed better growth during the early period after end-trimming than T
2 and maintained better growth through mid-period before declining.
Typical plant leaves exhibit strong light absorption near 450 nm and 670 nm and a remarkable increase in reflectance beyond 690 nm
(Slattery et al., 2017). Hyperspectral reflectance indicated that all end-trimmed treatments showed lower NIR reflectance than the control (Fig 7). This result suggests that increased light penetration into the canopy following end-trimming resulted in greater NIR absorption. Especially, at 30 DAT, hyperspectral reflectance beyond 700 nm decreased as end-trimming intensity increased.
Plants respond to shading conditions by increasing specific leaf area, enhancing antenna complexes, decreasing the chlorophyll a:b ratio and allocating nitrogen between pigment-protein complexes to maximize light absorption and photosynthetic efficiency (
Evans and Poorter, 2001). In T
0, higher leaf density in the upper layer likely led to increased shading within the canopy, initiating nitrogen redistribution to enhance light absorption in the lower layer. Conversely, T
3 showed consistently lower reflectance over 450-700 nm than other treatments, suggesting stronger absorption of blue and red lights.
Vegetation indices calculated using hyperspectral data and multivariate statistical methods that enable the interpretation of entire reflectance spectra are helpful in field research due to instant and nondestructive
(Sellami et al., 2022). Specifically, vegetation index analysis is one of the most widely studied and applied methods for estimating crop conditions from spectral reflectance data (
Morcillo-Pallarés et al., 2019). Vegetation indices such as Vog REI1, PSSRa, PSSRb, GNDVI and GCI are related to chlorophyll content
(Vogelmann et al., 1993; Mahlein et al., 2010). T
2 showed lower values at 13 DAT than at 42 DAT, suggesting a recovery from inhibition by strong trimming as new leaves are emerging. PSRI (Plant senescence reflectance index), which is associated with leaf aging, produced negative values close to zero when leaves were darker and shifted to positive values with senescence
(Merzlyak et al., 1999). The vegetation indices used to estimate light use efficiency include PRI (Photochemical reflectance index) and SIPI (Structure insensitive pigment index), which reflect the radiation use ability of plants (
Kohzuma and Hikosaka, 2018). PRI increased with end-trimming intensity at 13 DAT, then decreased at 42 DAT.
Vegetation indices also provide insights into canopy function. Chlorophyll-related indices such as Vog REI1, PSSRa and GNDVI decreased with increasing trimming intensity in early stages but later increased in T
3, suggesting delayed senescence. PSRI values further confirmed this, as T
3 maintained lower values, indicating slower leaf aging. These results agree with the research of
Merzlyak et al., (1999), who developed PSRI as a reliable indicator of senescence progression. PRI and SIPI with higher values reflected improved light-use efficiency in trimmed plants, especially in T
3, which decreased as leaves aged at 42 DAT. According to
Sims and Gamon (2002), PRI is linked to xanthophyll cycle activity and photosynthetic radiation-use efficiency, supporting the interpretation that canopy restructuring improved radiation balance. Overall, hyperspectral traits confirmed that trimming modulated not only structural but also functional responses of the soybean canopy.