In our previous study (Islam et al. 2021) drainage treatments showed up to 50% increase in achene yield in sunflower relative to the undrained treatment under waterlogged saline conditions in 2018–19 and 2019–20. In this paper, we were seeking physiological causes for these effects. Here, there was a consistency between years in response to drainage treatments in waterlogged saline soil for all parameters measured (ion concentrations in leaves, leaf water potential, shoot growth rate and soil solute potential), but values in the first (2018–19) and second year (2019–20) differed slightly due to variations in rainfall, temperature, soil salinity, and dates of planting, waterlogging and harvesting.
Concentrations of Na+ and K+ in leaves and their ratio (Na+/K+)
In both years, Na+ concentration and Na+/K+ ratio in all leaves at 7 DAFI and in the older leaves at 10 DASI and FL were decreased by drainage treatments compared with the undrained treatment, while the K+ concentration was increased (Fig. 1). At 7 DAFI, the lowest leaf Na+ concentration was found in the most drained (SSD + SD) treatment (838 mmol kg− 1 in 2018–19 and 485 mmol kg− 1 in 2019–20), whereas the highest was with undrained treatment (1374 mmol kg− 1 in 2018–19 and 825 mmol kg− 1 in 2019–20). The SD and SSD treatments had Na+ concentrations between the SSD + SD treatment and the undrained treatment. Similar trends were observed at 10 DASI and FL (Fig. 1). By contrast, leaf K+ concentrations at 7 DAFI, 10 DASI and FL were 894–972 mmol kg− 1 in 2018–19 and 980–1023 mmol kg− 1 in 2019–20 with most drained treatment and 764–926 mmol kg− 1 in 2018–19 and 788–938 mmol kg− 1 in 2019–20 with the undrained treatment (Fig. 1). In both years, the most drained treatment had the lowest Na+/K+ ratio at 7 DAFI (all leaves), 10 DASI (older leaves) and FL (older leaves), which was 9–37, 19–32 and 32–52% lower than the SSD, SD and undrained treatments, respectively (Fig. 1). However, younger leaves did not show significant changes in leaf ion concentration or the Na+/K+ ratio (Supplementary Material, Fig. S1).
Leaf Na+ concentration and Na+/K+ ratio during the season were significantly and negatively associated with achene yield, while leaf K+ positively correlated with the achene yield (Figs. 2 and 3; yield data reported in Islam et al. 2021). Leaf Na+ concentration at different times explained 46–81% and 40–81% of the yield variation in 2018–19 (Fig. 2) and 2019–20 (Fig. 3), respectively. The achene yield variation explained by K+ concentration ranged from 53 to 70% in 2018–19 and 50 to 75% in 2019–20. The Na+/K+ ratio explained 50–79% of the yield variation in 2018–19 and 64–86% in 2019–20. The younger leaves showed weak relationships between the ion parameters and achene yield compared with older leaves.
In both years, increased SEW30 and soil EC1:5 at 0–15 cm, and decreased Ψs at 0–15 cm were significantly correlated with increased Na+, decreased K+ and increased Na+/K+ in leaves at different times during the crop growing season (Table 1). However, younger leaves showed either no relationship or weak relationships, particularly at flowering. The SEW30 explained 60–76, 52–54 and 56–73% of the variation in 2018–19 and 47–84, 43–79 and 34–89% of the variation in 2019–20 in leaf Na+, K+ and Na+/K+ respectively. The strongest correlations with Na+ (r2 = 0.76), K+ (r2 = 0.63) and Na+/K+ (r2 = 0.73) were observed in the older leaves at 10 DASI in 2018–19. In 2019–20, leaf Na+ (r2 = 0.84) and Na+/K+ (r2 = 0.89) showed similar response but leaf K+ (r2 = 0.79) showed strongest correlation at 7 DAFI. In the case of soil EC1:5, the r2 values for Na+, K+ and Na+/K+ were 0.49–0.72, 0.39–0.84 and 0.54–0.74, respectively in 2018–19 and 0.46–0.86, 0.43–0.78 and 0.72–0.90 in 2019–20. The strongest correlations with Na+ (r2 = 0.72 in 2018–19 and 0.86 in 2019–20) and Na+/K+ (r2 = 0.74 in 2018–19 and 0.90 in 2019–20) were observed in the older leaves at FL in both years. The strongest correlation with K+ (r2 = 0.84) were observed in the young leaves at 10 DASI in 2018–19, but in 2019–20 the correlation was strongest (r2 = 0.78) in the older leaves at 10 DASI.
The soil Ψs also showed significant linear relationship with leaf Na+ (r2 values of 0.58–0.75 in 2018–19 and 0.44–0.85 in 2019–20), K+ (r2 values of 0.46–0.77 and 0.47–0.70 in 2018–19 and 2019–20, respectively) and Na+/K+ (r2 values of 0.42–0.78 in 2018–19 and 0.37–0.92 in 2019–20) in both years although no correlation was found at 7 DAFI in 2018–19 (Table 1). The strongest correlation with leaf Na+ and Na+/K+ occurred in the older leaves at FL in both years. However, the strongest correlation with leaf K+ differed between years. In 2018–19, it was highest in the young leaves at 10 DASI, while in 2019–20, it was highest in the older leaves at 10 DASI.
Table 1. Significance of effects of SEW30, soil EC1:5 at 0–15 cm and Ψs at 0–15 cm on Na+, K+ and Na+/K+ in leaves at different times during the growing season in 2018–19 and 2019–20.
|
Significance level with r2 values and direction of the slope (in brackets)
|
2018–19
|
2019–20
|
SEW30
(cm days)
|
EC1:5
(dS m–1)
|
Ys
(kPa)
|
SEW30
(cm days)
|
EC1:5
(dS m–1)
|
Ys
(kPa)
|
At 7 DAFI
|
Na+ (AL)
|
(+) 0.74***
|
(+) 0.57**
|
NS
|
(+) 0.81***
|
(+) 0.72***
|
(-) 0.63**
|
K+ (AL)
|
(-) 0.52**
|
(-) 0.41*
|
NS
|
(-) 0.79**
|
(-) 0.59**
|
(+) 0.47*
|
Na+/K+ (AL)
|
(+) 0.73***
|
(+) 0.60**
|
NS
|
(+) 0.88***
|
(+) 0.73***
|
(-) 0.61**
|
At 10 DASI
|
Na+ (YL)
|
NS
|
(+) 0.49*
|
(-) 0.58***
|
(+) 0.47*
|
(+) 0.46*
|
(-) 0.44*
|
K+ (YL)
|
(-) 0.54**
|
(-) 0.84***
|
(+) 0.77***
|
(-) 0.43*
|
(-) 0.72***
|
(+) 0.57**
|
Na+/K+ (YL)
|
NS
|
(+) 0.54**
|
(-) 0.62**
|
(+) 0.61**
|
(+) 0.73***
|
(-) 0.65**
|
Na+ (OL)
|
(+) 0.76***
|
(+) 0.62**
|
NS
|
(+) 0.84***
|
(+) 0.60**
|
(-) 0.48*
|
K+ (OL)
|
(-) 0.63***
|
(-) 0.74***
|
(+) 0.53**
|
(-) 0.71***
|
(-) 0.78***
|
(+) 0.70***
|
Na+/K+ (OL)
|
(+) 0.73***
|
(+) 0.69***
|
(-) 0.42*
|
(+) 0.89***
|
(+) 0.72***
|
(-) 0.60**
|
At FL
|
Na+ (YL)
|
NS
|
NS
|
NS
|
NS
|
NS
|
NS
|
K+ (YL)
|
NS
|
(-) 0.49*
|
(+) 0.46*
|
NS
|
(-) 0.43*
|
(+) 0.48*
|
Na+/K+ (YL)
|
NS
|
NS
|
NS
|
(+) 0.34*
|
NS
|
(-) 0.37*
|
Na+ (OL)
|
(+) 0.60**
|
(+) 0.72***
|
(-) 0.75***
|
(+) 0.70***
|
(+) 0.86***
|
(-) 0.85***
|
K+ (OL)
|
NS
|
(-) 0.39*
|
(+) 0.42*
|
NS
|
(-) 0.57**
|
(+) 0.60**
|
Na+/K+ (OL)
|
(+) 0.56**
|
(+) 0.74***
|
(-) 0.78***
|
(+) 0.68***
|
(+) 0.90***
|
(-) 0.92***
|
Abbreviations: DAFI = days after first inundation, DASI = days after second inundation, FL = flowering, AL = all leaves, YL = younger leaves, OL = older leaves, Ys = solute potential, * = P <0.05, ** = P <0.01, *** = P <0.001, NS = non-significant. For all relationships n = 12.
Stomatal conductance, leaf water potential and leaf chlorophyll content
There was a significant difference (P <0.001) in stomatal conductance (gs) among the treatments at 3 DASI but no difference at 10 and 17 DASI (Table 2). The combined drain treatment (SSD+SD) showed the highest gs (633 mmol m–2 s–1), while the lowest gs was with the undrained treatment (401 mmol m–2 s–1) at 3 DASI, compared with the gs of 552 mmol m–2 s–1 for SSD and 562 mmol m–2 s–1 for SD treatments. Drainage treatments also influenced the leaf water potential (Yleaf) at different times in both years (Table 3). In 2018-19, the lowest Yleaf was with the SSD+SD treatment (-1.29 to -1.40 MPa) during the season, while the highest Yleaf was with the undrained treatment (-1.18 to -1.25 MPa). There was no difference between SD and SSD treatments. Similar trend was observed in 2019–20, showing the lowest Yleaf of -1.27 to -1.58 MPa with the most drained treatment and the highest Yleaf of -1.16 to -1.33 MPa with the undrained treatment) throughout the cropping season.
Measurements of leaf chlorophyll content (LCC) are reported in the Supplementary Materials (Table S1) with the units of chlorophyll content index (CCI). In 2018–19, the most drained treatment had the highest LCC (12.0–17.2 CCI) throughout the season, whereas the undrained treatment had the lowest LCC (9.0–14.8 CCI) (Supplementary Material, Table S1). The variation in LCC between SD and SSD was not significant, although the values were higher than the undrained treatment but lower than the SSD+SD treatment. In 2019–20, the response of LCC to drainage treatments was similar to the previous year. The ranges of LCC during the crop growing season in SSD+SD, SSD, SD and undrained treatments were 13.7–23.1, 12.1–21.0, 12.2–21.4 and 11.3–19.2 CCI, respectively.
The achene yield was significantly associated with gs (positive correlation; one year of data only), Yleaf (negative correlation) and LCC (positive correlation) in both years (Supplementary Material, Table S2). The gs, Yleaf and LCC accounted for 17–68, 45–67 and 47–69% of the variation in achene yield, respectively. In most cases, the greatest variation was observed at 3 DASI at P <0.001. In addition, soil Ys, EC1:5 at 0–15 cm soil depth and SEW30 all showed significant relationships with gs, LCC and Yleaf at different times during the cropping season (Supplementary Material, Table S3). The Ys in soil was positively correlated with gs (r2 = 0.45 at 17 DASI) and LCC (r2 values of 0.33–0.67), and negatively correlated with Yleaf (r2 values of 0.44–0.63). The strongest correlations with gs, LCC and Yleaf were at 17 DASI, FL and FL, respectively. In contrast, the soil EC1:5 gave a negative correlation with gs (r2 = 0.34 at 17 DASI) and LCC (r2 values of 0.59– 0.74), and a positive correlation with Yleaf (r2 values of 0.42–0.63). The strongest relationships with gs, LCC and Yleaf were at 17 DASI, FL and FL, respectively. The SEW30 also showed negative correlation with gs (r2 values of 0.19–0.80) and LCC (r2 values of 0.57–0.90), and a positive correlation with Yleaf (r2 values of 0.25–0.68). The highest r2 values were observed at 3 DASI. The results also showed that plant height (Supplementary Material, Fig. S2) and leaf area (Supplementary Material, Fig. S3) were negatively correlated with Yleaf at 3 DASI and at FL in both years. Plant height and leaf area explained 58–64% and 56–72% of the variation in Yleaf, respectively.
Table 2. Effects of drains on stomatal conductance at different times in 2018–19.
Treatment
|
Stomatal conductance
(mmol m–2 s–1)
|
3 DASI
|
10 DASI
|
17 DASI
|
SSD+SD
|
633 a
|
944
|
955
|
SSD
|
552 b
|
946
|
952
|
SD
|
562 b
|
937
|
947
|
Undrained
|
401 c
|
900
|
898
|
P-value
|
<0.001
|
NS
|
NS
|
Means with the same letter are not significantly different at 5% level of significance. Abbreviations: SSD = subsoil drain, SD = surface drain, DASI = days after second inundation.
Table 3. Effects of drains on leaf water potential at different times in 2018–19 and 2019–20.
Treatment
|
Leaf water potential (MPa)
|
2018–19
|
2019–20
|
3 DASI
|
10 DASI
|
17 DASI
|
FL
|
3 DASI
|
10 DASI
|
17 DASI
|
FL
|
SSD+SD
|
-1.29 c
|
-
|
-
|
-1.27 b
|
-1.29 b
|
-1.35 b
|
-1.58 c
|
-1.27 b
|
SSD
|
-1.22 ab
|
-
|
-
|
-1.19 ab
|
-1.22 a
|
-1.27 a
|
-1.44 b
|
-1.19 ab
|
SD
|
-1.24 b
|
-
|
-
|
-1.22 b
|
-1.25 ab
|
-1.29 ab
|
-1.47 b
|
-1.22 b
|
Undrained
|
-1.18 a
|
-
|
-
|
-1.16 a
|
-1.20 a
|
-1.25 a
|
-1.33 a
|
-1.16 a
|
P-value
|
<0.01
|
-
|
-
|
<0.001
|
<0.01
|
<0.05
|
<0.001
|
<0.001
|
Means with the same letter are not significantly different at 5% level of significance. Abbreviations: SSD = subsoil drain, SD = surface drain, DASI = days after second inundation, FL = flowering.
Solute potential
The drainage treatments and soil depths significantly influenced solute potential (Ys) in the soil at different times during the crop growing season in 2018–19 (Fig. 4) and 2019–20 (Fig. 5). In both years, the treatments and soil depths showed interaction effects throughout the sunflower growing period, but no interaction was found between treatments and position in the plot, or between depth and position. In general, the lowest Ys was measured in the undrained treatment (-148 to -614 kPa in 2018–19 and -75 to -555 kPa in 2019–20), while in most cases, the highest Ys was with the SSD+SD treatment (-79 to -482 kPa in 2018–19 and -51 to -300 kPa in 2019–20) throughout the cropping season. Late in the sunflower growing season (FL–HRV), the Ys in the undrained treatment was either lower or similar to the most draining treatment (SSD+SD). The Ys was lowest at the upper soil (average depth 7.5 cm) and increased sharply at average depth 22.5 cm. The values from average depths of 22.5 to 52.5 cm remained almost stable or slightly increased in SD and undrained treatments, while the SSD and SSD+SD treatments showed decreasing trends. Finally, at the deeper soil (average depth 52.5 cm), there was a little variation of Ys among the treatments.
At average depth 7.5 cm, the higher Ys was with SSD+SD treatment (-87 to -482 kPa) and the lower values with the undrained treatment (-177 to -614 kPa) throughout the season in both years. The SSD and SD treatments had a similar Ys in the topsoil (average depth 7.5 cm), but at average depth 22.5 cm, the SSD treatment had higher Ys than the SD treatment. Increasing Ys at average depth 7.5 cm during the cropping season was significantly correlated with increasing sunflower yield with r2 values of 0.29–0.92 in 2018–19 and 0.51–0.77 in 2019–20 (Fig. 6). In both years, the strongest relationship was observed at FL.
Relative growth rate and dry weight of shoot
The dry weight data used to calculate the relative growth rate (RGR) of the shoots are reported in Supplementary Materials (Table S4). In 2018–19, the RGR of the shoot during the period of T1 (before inundation, BI, to 7 DAFI), T2 (7 DAFI to 14 DAFI) and T3 (10 DASI to 17 DASI) were highest in the most drained treatment (0.07–0.17 g g–1 d–1) (Fig. 7). In contrast, the lowest RGR occurred in the undrained treatment (0.03–0.11 g g–1 d–1) during T1 and T2, but the values were higher than both SD and SSD treatments at T3. There was no variation among the treatments at T4. Similarly, in 2019–20, the undrained treatment had the lowest RGR during T1–T3 (0.06–0.11 g g–1 d–1), and at T4, the values were similar to other treatments (Fig. 7). In contrast, the highest RGR was in the SSD+SD treatments, although the values were not different from SD and SSD treatments at T1, T3 and T4.
In both years, the shoot dry weight at different times during the study period showed strong negative correlations with SEW30 (r2 values of 0.85–0.95 in 2018–19, and 0.82–0.94 in 2019–20) and soil EC1:5 (r2 values of 0.66–0.93 in 2018–19, and 0.69–0.84 in 2019–20) (Supplementary Material, Figs. S4 and S5). In contrast, there were positive linear relationships with Ys in soil, which explained 41–90 and 52–86% of the variation in shoot dry weight in 2018–19 and 2019–20, respectively (Supplementary Material, Fig. S6). In addition, the higher shoot dry weights at different times throughout the season were significantly correlated with higher leaf K+ (r2 values between 0.46 and 0.73 in 2018–19, and 0.47 and 0.80 in 2019–20) and lower leaf Na+ (r2 values between 0.78 and 0.80 in 2018–19, and 0.47 and 0.82 in 2019–20) and Na+/K+ ratio (r2 values ranged 0.34–0.78 in 2018–19, and 0.65–0.90 in 2019–20) (Supplementary Material, Figs. S7 and S8). The LCC also showed positive and linear relationships with shoot dry weight, which accounted for 75–89 and 62–80% of the variation in 2018–19 and 2019–20, respectively (Supplementary Material, Fig. S9).