- 1China Meteorological Administration Training Center, Beijing, China
- 2College of Resources and Environmental Sciences, China Agricultural University, Beijing, China
- 3China Meteorological Administration, Henan Key Laboratory of Agro-meteorological Support and Applied Technique, Zhengzhou, China
- 4College of Ecology, Lishui University, Lishui, China
- 5Tianjin Climate Center, Tianjin, China
Winter wheat (Triticum aestivum L.) is a crucial crop that guarantees food supply in the North China Plain (NCP). As the frequency of extreme cold events increases, it is necessary to explore the freezing resistance of different wheat varieties in order to clarify planting boundaries and help with risk assessment. In this study, 2-year controlled experiments were conducted to explore the effect of freezing temperatures (Tair) and freezing durations on three winterness types. A set of indexes were used to characterize the subfreezing stress on wheat tiller, leaf, and final yield. Logistical regressions were used to quantify the temperature threshold for 10%, 30%, and 50% of freezing injury. The results showed that the lower temperature threshold of tiller (LT) varied from −9.6 to −15.9°C, −10.7 to −19.1°C and −11.4 to −21.2°C for LT10, LT30, and LT50, respectively. The difference between LT and yield loss (YL) indexes reduced with decreased winterness types and was −0.1 to 3.4°C, −0.7 to 2.1°C, and 0.3 to 0.9°C higher compared with YL thresholds for winterness, semi-winterness, and weak-winterness types, respectively. The average minimum soil temperature was 7.5, 4.8, and 4.2°C higher than Tair for 1-, 2-, and 3-day treatment, respectively. Soil effective negative accumulated temperature hours (TSEh) ranged from 6.9 to 12.0, 48.4 to 6.9, and 84.7 to 106.9°C·h for 10%, 30%, and 50% tiller mortality, respectively. Freezing treatment with Tair < −12, −9, and −8°C obviously decreased leaf Fv/Fm for the three varieties and Fv/Fm declined obviously after 5 days of recovery under field conditions. Our results provided multiple indexes for quantifying subfreezing damage in practical wheat production and could shed light on future risk assessment.
1 Introduction
Winter sowing crops have been found to possess a competitive advantage over spring sowing crops, primarily due to their early development in spring. This early development allows for the full utilization of spring soil moisture and leads to earlier maturation (Pirjo et al., 2011; Rizza et al., 2011). This advantage is particularly significant in the North China Plain (NCP), where the double-cropping system of winter wheat (Triticum aestivum L.) and summer maize is prevalent (Wang et al., 2012; Gao et al., 2019). The NCP is responsible for approximately 68% of China’s wheat production, with its planting area accounting for 58% of the nation’s total (Sheng, 2016; Hong et al., 2023). Hence, ensuring food security in this region is of utmost importance. However, the success of winter crop growing is threatened by negative climatic conditions that can result in a decrease in overwinter survival rate, crop vigor, and final yields (Licker et al., 2013; Hong et al., 2023). To counteract the potential negative effects of climate change, it has been suggested that weaker winter varieties could be planted to achieve higher yields (Liu et al., 2009; Macholdt and Honermeier, 2017), given the rising temperatures in winter and spring (You et al., 2010; Sommer and Lengfellner, 2010). However, climate change also brings about increased climatic uncertainty and a higher frequency of extreme weather events (Judy et al., 2020). For example, China has experienced extreme cold events during the past decade, in 2008, 2010, 2016, and 2017, all of which have had a significant impact on wheat cultivation. The higher spring and autumn temperatures associated with climate change can delay plant hardening and accelerate dehardening, thereby making the crop more vulnerable to cold spells for an extended period (Rapacz et al., 2014). Considering these factors, it becomes clear that the benefits of overwintering crops become more complex to assess.
Winter survival of plants in the field is a complex process influenced by multiple factors, not solely dependent on temperature (Bergjord Olsen et al., 2018; Beil et al., 2019). The ability of winter wheat to survive winter conditions is referred to as winter hardiness. Winter hardiness is a complex trait that reflects tolerance to various combinations of freezing, desiccation, anoxia, ice encasement, and resistance to disease (Armonienė et al., 2013; Zheng et al., 2018). Research has shown that winter wheat is vulnerable to low-temperature damage, from both short exposures to extreme cold and prolonged exposure to milder sub-zero temperatures (Skinner and Mackey, 2009). Furthermore, repeated cycles of freezing and thawing can exacerbate the damage (Vico et al., 2014). It is important to recognize that the degree of injury is determined by both temperature and exposure time (Skinner and Mackey, 2009; Skinner, 2014). The sensitivity to low temperatures varies among different plant organs (Gusta and Wisniewski, 2013). In order to evaluate the freezing tolerance of winter wheat, most experimental studies have focused on the LT50 (temperature at which 50% of the plants are killed) threshold (Armonienė et al., 2013; Fowler et al., 2014; Bergjord Olsen et al., 2018), with particular emphasis on crown tissue as it is crucial for winter survival and the regeneration of other damaged plant organs (Li et al., 2015). However, it is also important to consider yield loss (YL) as an index for evaluating the severity of cold damage (Min et al., 2020; Li et al., 2021). Freezing injury can result in significant YL, outweighing the impact of other environmental factors. To estimate freezing tolerance, chlorophyll fluorescence-based techniques have been developed as reliable, non-invasive, and easy-to-use tools (Rapacz and Woźniczka, 2009; Rapacz et al., 2015; Jaimez et al., 2023). These techniques provide researchers with indirect information about the structure and function of the photosynthetic apparatus, aiding in the assessment of freezing tolerance (Rapacz et al., 2015). It is crucial to consider the plant as a whole and its surrounding environment to gain a comprehensive understanding of the factors affecting wheat’s cold tolerance. By addressing these aspects and employing appropriate protocols, researchers can study and identify sources of improvement in cold hardiness.
Numerous investigations have been conducted to assess the cold tolerance of winter wheat varieties (Mu et al., 2015; Zheng et al., 2018). However, most of these studies have focused on the ability of the plants’ tillers/crown to withstand extremely low temperatures. Considering the changes in prewinter growing conditions due to climate change, such as the hardening process, we hypothesize that the freezing tolerance of winter wheat may have been altered. Additionally, the indexes obtained from the freezing chamber and greenhouse may not accurately reflect the field conditions, which limits their practical application (Mu et al., 2015; Zheng et al., 2018). Moreover, there have been very few studies that have evaluated the freezing tolerance of different organs of winter wheat, especially in a non-destructive manner. Therefore, the objectives of this study are twofold: (1) to assess the freezing tolerance of various organs of winter wheat using indexes that accurately correspond to field conditions, and (2) to evaluate the impact of different lengths of freezing days on the survival of winter wheat.
2 Material and methods
2.1 Experimental design
Three winter wheat varieties with different winterness types were sown in round plastic pots (diameter = 25 to 30 cm), at 2 and 1 October in 2015 and 2016 at Shangzhuang and Baodi, located in the NCP. The winterness variety “Nongda211” (ND211) is typically planted at the north region of the NCP, which can resist mild cold freezing temperatures. The semi-winterness variety “Zhengmai366” (ZM366) is usually planted at south-middle regions of the NCP. The weak-winterness variety “Yanzhan4110” (YZ4110) is normally planted at south part of the NCP, with weakest resistance to cold freezing temperatures. The under part of plastic pots was removed and then filled with perforated and well-drained sandy loam soils to allow the supply of soil water from deeper soil layers. Initially, 150 kg N/ha of compound fertilizer (N:P:K = 12:12:12) was applied. Each pot contained 25 seeds and was thinned to 20 seeds per pot after emergence, giving a plant population of approximately 280 to 400 plants/m2. The experiment materials were watered every 10 days to prevent from drought stress before winter. Additionally, chemical fungicides were used during the tillering stage in line with the conventional field production of winter wheat to inhibit disease or pests.
All pots were placed outdoors under field growing conditions right after sowing until the day of freezing tests so that the hardening process was similar for all the plants. The pots were placed in pits in the soil so that the plants stayed at the same level as the grain field stubble surrounding them. To prevent the pots from being firmly frozen into the surrounding soil and being impossible to remove during winter, the pits were lined with a layer of plastic mat. The low temperature always struck by late December to early February at Baodi. A layer of straw was placed on top of the field crop when minimum temperature was lower than −10°C according to the weather forecast in the 2016–2017 experiment season to prevent the crops from accumulating freezing injury due to multi-freezing events.
2.2 Freezing tolerance experiments
Freezing tolerance for three varieties was tested by late December and early January, when the hardening process was finished and winter wheat had gained its full freezing tolerance, with three replications for each freezing test. Wheat crown was normally located in between 2 and 3 cm following normal seed drilling (Persson et al., 2017); thus, soil temperature at 2.5 cm depth was recorded every 1 h by automatic weather stations using temperature loggers placed in the plot with an expected accuracy of ±0.1°C. The occurrence of long-lasting cold event decreased since climate warming during winter in the NCP while the short-duration low-temperature events still existed, especially in the northern part of the NCP. The frequency of cold event at our experiment site is shown in Table 1. Specifically, the frequency of −10°C low temperature was higher than 90% (Table 1). Thus, the plants were first held in the freezer at 0°C for 1 h, and temperatures were reduced at a rate of 0.8 to 2.0°C/h (targeted temperatures from −8 to −20°C for different varieties). The temperatures were kept at the target temperatures for 3 h and then gradually raised to 0°C at the same rate. After the first duration of 1 day, 1/3 of the plants were moved out. The temperature change was repeated for 2- and 3-day treatments (Table 2). After the freezing tests, all plants were transferred to the field and kept for regrowth until maturity. The plastic mat was removed to allow access of root to the deeper soil layer. The survival proportion (mean of three replications) was counted after 1 week of green-up stage, and plants were kept in the field until maturity. Other management practices were kept the same as the field-grown wheat, irrigated twice at the jointing stage and at the flowering stage with 60 mm each time; 75 kg N/ha of urea was spread to the surface at the jointing stage (Zadock = 31).
2.3 Freezing tolerance indexes
To assess the damage degree of subfreezing temperatures, four different indexes were used in this paper: (i) the mortality rate of tillers (LT) was expressed as a numerical scale from 0 (no tillers were killed) to 1 (all tillers were killed) (Bergjord Olsen et al., 2018; Zheng et al., 2018) (Equation 1); (ii) the YL was the proportion of yield decreased compared to control treatment (Equation 2) (Wu et al., 2019; Dahal et al., 2021); (iii) soil degree hours as the sum of effective subfreezing temperatures at the crown depth during the freeze tests (Equations 3, 4); and (iv) the chlorophyll fluorescence index (see section 2.4) (Rapacz et al., 2015; Bai et al., 2021). The logistic function (Equation 5) was applied as it could better reproduce the effect of low temperature on crop growth (Tudela and Santibáñez, 2016; Zheng et al., 2018; Wang et al., 2021). During the 2015–2016 growing season, the three varieties experience another cold stress after the freezing test (Tmin = −15.3°C for 2 days), which caused a tiller death of 0.19, 0.39, and 0.48 for the CK treatment of the three varieties, which was excluded when calculating the LT caused by the controlled freezing chamber.
where F1 represents the tiller number before freezing experiment treatment and F2 is the tiller number after reviving for each treatment.
where Yck represents the harvested yield of control treatment (plants grown under field conditions) and Yi is the harvested yield for each freezing treatment.
where TSEh is the soil effective negative accumulated temperature hours (°C•h) and TSEd is the soil effective negative accumulated temperature days (°C•d). Tc is the threshold temperature of crown to withstand freezing injury (°C), Ti is the actual soil temperature at crown depth at ith h, Tsmini is the minimum soil temperature at ith d. When Ti > Tc and Tsmin > Tc, TSEh and TSEd would not accumulate.
where y represents the dependent change (LT or YL). The temperature (–b/a) corresponded to the temperature that causes 50% death or YL.
2.4 Chlorophyll fluorescence index
The chlorophyll fluorescence induction parameters were determined on the second leaf of the main leaf using OS-30P (Opti-Science inc., USA) using the Fv/Fm kinetics option. Leaves were dark-adapted for 20 min with a leaf clip holder placed on the central part. By exposure to saturated red light with a photon flux density of approximately 10,000 µmol/m2•s to estimate the initial (F0) and maximum (Fm) fluorescence values. Fv was calculated as (Fm − F0). The parameter Fv/Fm was used to assess the efficiency of excitation energy capture by an open PSII reaction center, which is the maximum capacity of light-dependent charge separation of PSII. It was supposed to decrease after freezing injury. Firstly, Fv/Fm was measured before freezing tests, and the second measure was performed after the freezing tests when all the samples were held at 0°C for 2 h (Bai et al., 2021). It was then recorded after 1, 5, and 10 days of recovery under field conditions. Five and three replications were performed for 2015–2016 and 2016–2017, respectively.
3 Results
3.1 Evaluation of overwintering freezing tolerance of winter wheat by tiller
Figure 1 shows the indexes of overwinter mortality, calculated in Section 2 for the 2-year experiments with different freezing temperatures and freezing days for three types of winter wheat varieties. Experiments were carried out at the end of December and at the start of January when the acclimation has completed under natural growing conditions. Winter damage increased with lower temperatures and extended freezing days. Mortality in individual trials ranged from 0 to 1 for “ZM366” and “YZ4110” during the 2-year experiments, indicating that the test conditions used were sufficient enough to result in significant different levels of survival because of variation in minimum temperature and the time held at the minimum temperature. The tiller mortality was small regardless of freezing days when Tmin was higher than −14, −12, and −10°C for ND211, ZM366, and YZ4110, respectively (Figures 1A-C). Specifically, for “ND211”, cold stress of approximately −12°C caused 2%–6% mortality, after which tiller mortality increased from 5% to 62% with increased cold stress and prolonged freezing days (Figure 2A). The effect of freezing days becomes obvious when Tmin ≤−16°C, where average tiller mortality increased by 12% with additional one freezing day (Figure 2A). Tiller mortality increased significantly with minimum temperature reduced from −12 to −16°C, while tiller mortality was ≤6% for freezing under −11°C. For cold stress lower than −14°C, average tiller mortality increased by 13.8% with an additional one freezing day for ZM366 (Figure 2B). It caused a tiller mortality of 3% to 100% when exposed to −9 to −14°C freezing stress for Yanzhan4110. The tiller mortality was within 25% regardless of freezing temperature or freezing days when Tmin ≥−11°C while tiller mortality increased significantly when exposed to −12°C stress and reached 100% death if freezing for 3 days under −14°C treatment (Figure 2C).
Figure 1. Relationship between minimum air temperature and tiller mortality during middle winter for different winterness types. (A) Nongda211, (B) Zhengmai366, (C) Yanzhan4110. F1 (hollow blue diamond), F2 (hollow green diamond), and F3 (hollow red diamond) indicated freezing for 1, 2, and 3 days, respectively, during the 2016–2017 growing season, F3-15 (solid black diamond) indicated freezing for 3 days during the 2015–2016 growing season.
Figure 2. Relationship between minimum air temperature and yield loss in middle winter for winterness types. (A) Nongda211, (B) Zhengmai366, (C) Yanzhan4110, F1 (hollow blue diamond), F2 (hollow green diamond), and F3 (hollow red diamond) indicated freezing for 1, 2, and 3 days, respectively, during the 2016–2017 growing season, F3-15 (solid black diamond) indicated freezing for 3 days during the 2015–2016 growing season.
To note, the winterness variety “Nongda211” showed greater tolerance for lower freezing damage (−13.1 to −21.2°C), followed by the semi-winterness variety “Zhengmai366” (−11.0 to −16.3°C) and the weak-winterness variety “Yanzhan4110” (−9.6 to −13.9°C). The freezing threshold temperatures of LT10, LT30, and LT50 ranged from −9.6 to −15.9°C, −10.7 to −19.1°C, and −11.4 to −21.2°C, respectively. As for increased freezing days, the threshold temperatures ranged from −10.9 to −21.2°C, −10.3 to −19.5°C, and −9.6 to −18.3°C, respectively. The differences between LT30 and LT10 were approximately 2.9–3.2°C, 1.4–2.4°C, and 1.4–1.8°C for the three varieties, while the differences between LT50 and LT30 were smaller, indicating that the lower the degree of low temperature, the greater the contribution of unit low temperature to mortality rate (Table 3).
3.2 Evaluation of overwintering freezing tolerance of winter wheat by accumulated soil effective negative accumulated temperature
Wheat crown tissue has the strongest tolerance for freezing injury, which was usually within 2–3 cm in the soil while temperatures fluctuated a lot in the top layer. The soil temperature at 2.5 cm was used to quantify the freezing injury on wheat plant. The average realized minimum soil temperatures were −2.2, −2.7, −4.6, −7.0, −8.0, −9.5, and −12.4°C in trials with target temperatures of −9, −10, −12, −14, −16, −18, and −20°C when freezing for 1 day. With extended freezing days, crown sensed minimum temperatures decreased, Tsmin during the third freezing day would be 2.7 to 4.0°C lower than was for the first freezing day, with little difference for Tsmin between the second day and the third day (Figure 3). The average realized minimum soil temperature was 7.5, 4.8, and 4.2°C higher than the designed minimum air temperature for the 1-, 2-, and 3-day treatment, respectively. Thus, the actual crown sensed soil temperature was used for quantifying tiller mortality instead of designed air minimum temperature irrespective of freezing days. According to experiment data, plants were killed when observed crown sensed temperature was approximately −3, −6, and −8°C for “YZ4110”, “ZM366”, and “ND211”, respectively, and thus were treated as effective freezing soil temperature for analysis later. Overwinter mortality increased linearly with increased effective negative accumulated temperature hours or days (Figure 4).
Figure 4. Relationship between soil effective negative accumulated temperature and tiller mortality for different variety types at middle winter. (A) Nongda211 (°C•h), (B) Nongda211 (°C•h), (C) Zhengmai366 (°C•h), (D) Zhengmai366 (°C•d), (E) Yanzhan4110 (°C•d), (F) Yanzhan4110 (°C•d).
Specifically, soil effective negative accumulated temperature hours (TSEh) for “YZ4110” would be 12.0 to 84.7°C·h, which is equal to a soil effective negative accumulated temperature days (TSEd) of 2.0 to 11.5°C·d, resulting in 10% to 50% mortality. For “ZM366”, TSEh would be 8.0 to 86.5°C·h, which is equal to a TSEd of 1.8 to 12.3°C·d, resulting in 10% to 50% mortality. “Nongda211” has a stronger resistance for freezing injury with a higher value of TSEh (6.9 to 106.9°C·h) and TSEd (2.0 to 11.5°C·d) (Table 4).
Table 4. The threshold of freezing injury based on soil effective negative accumulated temperature hours (TSEh) and soil effective negative accumulated temperature days (TSEd).
3.3 Evaluation of overwintering freezing injury on winter wheat yield loss
Winter wheat yield was more sensitive to freezing temperatures than crown tissue, which responded negatively to decreased temperatures for all three varieties. Statistical analysis showed that minimum temperature, time held at minimum temperature, and the interaction of these two factors had a significant effect on final YL (p = 0.05). YL increased with prolonged freezing days, especially under lowered temperatures (Figure 2).
Averaged YL was <10% when Tmin were higher than −12°C, 10%–20% when suffering from cold stress of −14 to −16°C, while a freezing stress of −20°C resulted in over 50% of YL for ND211. Specifically, a cold stress of approximately −12°C caused 10% YL, after which YL increased from 12.6% to 57.9% with increased cold stress and prolonged freezing days (Figure 2A). The effect of freezing days becomes obvious when Tmin ≤14°C, where average YL increased by 8% with one additional freezing day (Figure 2A). YL increased significantly with minimum temperature reduced from −11 to −16°C, while YL was approximately 10% for freezing under −10°C for ZM366. For cold stress lower than −12°C, average YL increased by 7% with one additional freezing day (Figure 2B). It caused a YL of 4% to 100% when exposed to −9 to −14°C freezing stress for Yanzhan4110. The YL was within 60% regardless of freezing temperature or freezing days while YL increased significantly when exposed to −14°C stress and reached 100% death if freezing for 2 and 3 days (Figure 2C).
Specifically, the freezing threshold temperatures of YL10, YL30, and YL50 ranged from −11.9 to −21.3°C, −8.9 to −15.9°C, and −8.7 to −13.4°C for ND211, ZM366, and YZ4110, respectively. As for increased freezing days, the threshold temperatures ranged from −8.7 to −12.5°C, −10.1 to −17.9°C, and −10.9 to −21.3°C for freezing 1, 2, and 3 days, respectively. The differences between YL30 and YL10 were approximately 3.8–5.4°C (ND211), 3.1–3.2°C (ZM366), and 1.3–2.0°C (YZ4110) for the three varieties, while the differences between YL50 and YL30 were smaller (2.5–3.4°C, 2.0–2.1°C, and 0.8–1.3°C), indicating that the lower the degree of low temperature, the greater contribution of unit low temperature to YL (Table 5).
3.4 Evaluation of overwintering freezing injury on winter wheat by Fv/Fm
Figure 5 shows the changes of Fv/Fm under diverse freezing treatments. Winterness variety exhibited evident advantage over the other two varieties. For most of the treatments, no obvious trend was found after freezing tests, but it caused a dramatic decrease of Fv/Fm after 5 days of recovery in the field for most of the cases, while the Fv/Fm response was different among different freezing days.
Figure 5. Effect of different freezing treatments on Fv/Fm ratio for different variety types in middle winter. The first three rows indicated the 2015–2016 growing season. (A–C) indicated freezing for 3 days in 2015 for “Nongda211”, “Zhengmai366”, and “Yanzhan4110”, respectively. (D–F) indicated freezing for 1 day for “Nongda211”, “Zhengmai366”, and “Yanzhan4110” during the 2016–2017 growing season, respectively. (G–I) indicated freezing for 2 days for “Nongda211”, “Zhengmai366”, and “Yanzhan4110” during the 2016–2017 growing season, respectively. (J–L) indicated freezing for 3 days for “Nongda211”, “Zhengmai366”, and “Yanzhan4110” during the 2016–2017 growing season, respectively. B means before treatment; A means after treatment; 1, 5, and 10 mean recovery after 1, 5, and 10 days, respectively.
For “ND211”, cold stress of 1 day did not have a significant influence on Fv/Fm except for −20°C. The Fv/Fm value fluctuated for other temperature treatments, while it decreased after 5 days of recovery when freezing at −20°C. Although the Fv/Fm value fluctuated after 1 day of recovery for the treatment of freezing for 2 days, a significant drop was detected for −18 and −20°C treatments after 5 days of recovery, and for −16°C after 10 days of recovery. Freezing for 3 days with a temperature below −14°C resulted in a major fall of Fv/Fm after 1 day of recovery. It seemed that the effect of low temperature was reversible when freezing at −10°C regardless of freezing days (Figures 5A, D, G, J).
The Fv/Fm values showed a slight increase after freezing treatment except freezing at −16°C for 3 days for “ZM366” for the 2016–2017 growing season. No obvious difference was observed for each low-temperature treatment after 1 day of recovery for freezing by 1 or 2 days, but after 5 days of recovery, the difference was apparent at −14°C and −12°C for 1-day and 2-day treatments, respectively. The Fv/Fm value was half of that before freezing after 10 days of recovery. When suffering at a temperature higher than −10°C for 3 days, Fv/Fm did not decrease much after 5 days of recovery for both experimental years. The injury was obvious after 10 days of recovery, with Fv/Fm valued from 0.38 to 0.45. It revealed that the effect of low temperature was irreversible with low-temperature treatment from −10 to −16°C regardless of freezing days (Figures 5B, E, H, K).
The Fv/Fm did not change much after freezing test both for 1-day and 2-day treatment for “Yanzhan4110”. The lowering of Fv/Fm was already evident after 5 days of recovery for temperature <−10°C. The injury was distinct after freezing for 3 days for both growing seasons, which continued to decrease under field recovery. The harm of cold stress was clear when Tmin < −12°C (Figures 5C, F, I, L).
4 Discussion
In this study, a set of freezing indexes was applied to evaluate the sub-freezing injury on three types of field-grown winter wheat in the NCP. Freezing tolerance was determined by LT, accumulated soil effective negative accumulated temperature hours/days (TSEh and TSEd), YL, and chlorophyll fluorescence Fv/Fm of plant samples, which were cold acclimated under natural conditions. As concerns had been stated that current experiments did not reflect true complexity of cold acclimation (Gusta and Wisniewski, 2013), plants cold acclimated under natural conditions seem to be a better way to reflect the plants suffering under field conditions (Rapacz et al., 2015). Previous studies showed that cold tolerance gained under controlled environment was stronger than that gained under field conditions (Rapacz et al., 2015; Short and Wolyn, 2022), which limit the prevailing of these indexes in reality. Actually, existing studies always put the after-treatment plants in the greenhouse, where the circumstantial environment was totally different from the field. Thus, in our experiment, we put the plants in the field after freezing treatment (with no additional freezing threat) to obtain disaster indicators that are more in line with practical production applications. Additionally, the existing studies mainly emphasize the temperature indicator while we tried to develop an index system that focuses on varied organs and quantifying methods.
4.1 The employment of multi-index thresholds for quantifying freezing injury
When evaluating the potential risks of winter injury, the LT50 parameter has commonly been utilized to quantify the likelihood of damage occurring when temperatures reach detrimental levels. This approach offers a straightforward and easily understandable method for characterizing the frequency at which harm may be expected to manifest. However, it has been argued that the 50% threshold for injury may be excessive for agronomic planting areas (Waalen et al., 2011; Vico et al., 2014). Furthermore, LT50 does not account for the severity of the damage, as it solely indicates whether or not damage has occurred. Consequently, this approach may yield similar risk levels across different locations where the damaging temperature has been reached (Bergjord Olsen et al., 2018). To address these limitations, a novel assessment method has been proposed in this study. The S curve was employed to assign slight injury, moderate injury, and severe injury ratings of 10%, 30%, and 50%, respectively, to evaluate a specific subfreezing damage for each cultivar. Accordingly, it has been suggested that considering 5% of kill as slight injury would align with the established standards for agrometeorological disasters in China. However, it should be noted that both overestimation and underestimation may occur at the extremities of the response curve (Tudela and Santibáñez, 2016).
Actually, the LT50 promoted by this study was a lot higher compared with the existing studies (Mu et al., 2015; Zheng et al., 2018). In our study, all wheat varieties underwent cold resistance training under natural conditions. After the low-temperature treatment, the wheat plants were placed in the field to reach their final maturity. However, wheat growth in the field is a complex process that involves various factors, including exposure to low temperatures (Gusta and Wisniewski, 2013). When wheat seedlings are exposed to freezing damage and left in the field, external factors such as strong winds and sunny days can exacerbate dehydration and lead to plant death. As a result, the freezing resistance index obtained from field experiments generally appeared weaker compared to those from greenhouse experiments (Mu et al., 2015; Zheng et al., 2018). For instance, Mu et al. (2015) reported that the LT50 values for winter wheat varieties were −21.1°C, −18.5°C, and −14.7°C for soil temperature, whereas we reported the similar threshold range using air temperature but not soil temperature. The reason may be that the freezing resistance of tillers is likely to be weaker compared to that of crown nodes/whole plants. Furthermore, climate change has significantly impacted China’s weather patterns since 1980, resulting in the absence of stable snow cover in the NCP during winter (Zhang et al., 2015). Many freezing damage indexes in China are based on research conducted before the 1980s when winter wheat varieties exhibited stronger resistance (Zheng, 1981; Gong et al., 1982). Moreover, in field production, the cumulative effects of repeated mild low-temperature threats throughout the long winter, along with temperature fluctuations and nutrient consumption, gradually weaken and relieve resistance over the entire overwintering period (Ferrante et al., 2024). Dry soil surfaces in winter, combined with strong winds and soil cracking, can lead to water loss and the natural death of overwintering seedlings. Therefore, even if the temperature in the field is slightly higher than the critical lethal temperature identified in the greenhouse, severe freezing damage and the death of seedlings can still occur.
4.2 The promotion of crown perceived soil temperature as freezing index
The soil temperature was less utilized when quantifying crop growth, since freezing tolerance is often revealed using LT50 of minimum air temperature in most winter damage studies. However, recent research has highlighted the significant impact of the duration of exposure to the minimum temperature on freezing tolerance (Skinner and Mackey, 2009; Skinner and Bellinger, 2011; Skinner, 2014). Furthermore, variations in soil temperature under different crop growing conditions have limited the applicability of air temperature measurements in diverse production areas (Yang et al., 2021). To address these limitations, researchers have explored the use of the sum of soil temperature below a threshold (often 0°C) to assess prolonged exposure to relatively milder temperatures. Nonetheless, the absence of cultivar-specific parameters has hindered the practical application of this approach (Vico et al., 2014). In this study, we propose the use of TSE, which incorporates cultivar-specific threshold temperatures derived from experimental results. This index provides a more accurate and comprehensive assessment of freezing tolerance. The TSE indexes are recommended for locations where shallow soil temperatures can be obtained, as the translation of air temperature to soil temperature can vary under different weather conditions (e.g., sunny, cloudy, rainy, or with snow cover). Considering the significant geophysical differences in winter wheat cultivation in the NCP, the varying seeding dates pose an additional concern. The density of plant cover before winter can greatly affect soil temperature dynamics, with denser covers resulting in less drastic changes compared to plants sown later with fewer leaves and a colder crown sensed temperature (Luo et al., 2018; Islam et al., 2022). The freezing of soil is a rapid process, and the minimum soil temperature realized became colder with prolonged freezing time (Figure 3). Occasionally, abnormal points may be observed in Figure 4 for all cultivars. It is worth noting that smaller accumulated TSE values resulted in higher death rates, while some larger TSE values led to lower death rates. This phenomenon can be explained by the freezing tolerance mechanism in wheat, which is activated only after exposure to a sufficiently low temperature for a sufficient duration (Skinner, 2014). Assessing the risks based on TSE provides a more stable comparison of freezing conditions across diverse climatic conditions. Therefore, incorporating the TSE indexes into winter damage studies can enhance our understanding of freezing tolerance in crops and inform decision-making regarding cultivation practices.
4.3 The application of yield loss index to measure low-temperature damage
Currently, the focus of most indexes lies in tiller death, disregarding the growth of after-treatment. However, in reality, YL plays a significant role in ensuring food security. The index of YL exhibits a similar trend to LT (low temperature), although the threshold temperature differs a lot. The reason may that invisible injury may occur in the growth cone, resulting in adverse effects on green-up growth. This effect is particularly severe in later-seeded plants in actual production, as excessive biomass consumption due to freezing injury leaves insufficient biomass for re-tillering, even if the crown has not been completely freezing to death. Shoot regrowth reflects both the number of surviving plants and their vigor (Waalen et al., 2011; Kinugasa and Gantsetseg, 2023); thus, excessive consumption during winter can have a considerable impact on later yields. Existing studies pointed out that freezing damage was composed of two components: the proportion of plants destroyed by freezing and the YL observed in damaged but surviving plants (Veisz et al., 2001). The side spikes, rather than the main spike, tend to be more susceptible to freezing, which cannot be accurately reflected solely by LT50. Additionally, the root systems of cereals are more vulnerable to freezing than the shoots, and cold freezing temperatures between −5 and −9°C can cause their death. The lack of root initiation in the field rather than just the death of a proportion of crowns are the cause of winter kill, and thus, differences existed between LT and YL.
4.4 The utilization of the non-destructive index for assessing low-temperature threats
The complexity of the response to cold damage is often underestimated (Gusta and Wisniewski, 2013), highlighting the need for considering the effect of cold damage on the entire plant. To evaluate the cold tolerance of leaves, this study employed the chlorophyll fluorescence Fv/Fm as a reliable, non-invasive, and easy-to-use tool (Rapacz et al., 2015; Jaimez et al., 2023). Our results revealed no distinct difference after the freezing tests, but a noticeable change in Fv/Fm was observed after a 5-day recovery period. The selected value of 5 days of recovery as a freezing indicator was supported by existing studies (Rapacz, 2007; Rapacz and Woźniczka, 2009), which stated some concerns on using Fv/Fm as a wheat freezing indicator when measures were made directly after freezing. Additionally, Rapacz et al. (2015) found that the highest correlations were obtained when chlorophyll fluorescence measurements were taken in mid-winter, which coincides with the time frame of our experiment. Although this indicator represents photo-inhibitory damage with a relatively short time lag, it provides farmers with the advantage of early detection of freezing injury before the next spring, enabling them to implement remedial measures to mitigate possible significant damage (Rizza et al., 2011; Wang et al., 2016) and adoption of agronomic management practices to promote re-tiller growth and compensate for YL.
The freezing stress experiments described in this paper were conducted in a controlled freezing chamber without snow cover on the surface. In field conditions, however, snow cover can sometimes be present, although it is not stable in the NCP, particularly under the influence of climate change. In the NCP, snow cover is typically shallow and does not significantly contribute to soil warming when the snow depth is less than 5 cm (Persson et al., 2017). Under the impact of climate change, the occurrence of freeze–thaw overwinter damage is on the rise (Bergjord Olsen et al., 2018); we confronted this kind of injury in the 2015–2016 experiment, and the mechanism of this type of damage should be further studied in the future. Moreover, it would be ideal to treat all three varieties at the same temperature, even though this would have increased the number of treatment combinations in the study; thus, we suggest future studies to pay attention to this kind of shortcoming.
5 Conclusion
This study aimed to investigate the tolerance thresholds of various winter wheat organs to freezing injury using a controlled freezing chamber. Following exposure to sub-freezing temperatures, the treated plants were grown in the field until maturity. The relationship between air temperature, soil temperature, and plant damage (10%, 30%, and 50%) was established, and the effects of temperature threats on leaf fluorescence (Fv/Fm) were quantified. The freezing tolerance varied among different plant organs. The lower-temperature threshold of tiller (LT) ranged from −9.6 to −15.9°C, −10.7 to −19.1°C, and −11.4 to −21.2°C for LT10, LT30, and LT50, respectively. The difference between LT and YL indexes decreased with less winterness, with values of −0.1 to 3.4°C, −0.7 to 2.1°C, and 0.3 to 0.9°C higher for winterness, semi-winterness, and weak-winterness types, respectively. The minimum soil temperature decreased with additional freezing days under the same target air temperature. The hourly soil effective negative accumulated temperature (TSEh) ranged from 12.0 to 84.7°C·h (YZ4110), 8.0 to 86.5°C·h (ZM366), and 6.9 to 106.9°C·h (ND211) for 10% to 50% tiller mortality. In most treatments, no significant trend was observed in Fv/Fm following freezing tests, but a dramatic decrease was observed after 5 days of recovery in the field. Overall, this study provides easily recorded indexes for assessing freezing injuries in winter wheat organs and demonstrates the non-destructive use of Fv/Fm in advance for freezing quantification and enable them to implement remedial measures to mitigate the possible significant damage.
Data availability statement
The data presented in this study are available on request from the corresponding author. The data are not publicly available because they need to be used in future work.
Author contributions
LW: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing. WY: Conceptualization, Data curation, Investigation, Methodology, Resources, Writing – review & editing. CC: Funding acquisition, Methodology, Project administration, Resources, Writing – original draft, Writing – review & editing. LF: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing. JY: Data curation, Investigation, Resources, Writing – original draft. DZ: Conceptualization, Data curation, Investigation, Methodology, Resources, Supervision, Writing – original draft. FY: Funding acquisition, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work is supported by National Natural Science Foundation of China (31371528), and the Lishui University Talent Launch Fund Project (6604CC01Z).
Acknowledgments
The authors thank the staff at the Tianjin climate center for providing help in conducting the experiment.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Armonienė, R., Liatukas, Ž, Brazauskas, G. (2013). Evaluation of freezing tolerance of winter wheat (Triticum aestivum L.) under controlled conditions and in the field. Zemdirbyste-Agriculture 100, 417–424. doi: 10.13080/z-a.2013.100.053
Bai, R., Lin, Y., Jiang, Y. (2021). Diverse genotypic variations of photosynthetic capacity, transpiration and antioxidant enzymes of lily hybrids to increasing salinity stress. Scientia Hortic. 280, 109939. doi: 10.1016/j.scienta.2021.109939
Beil, C. T., Anderson, V. A., Morgounov, A., Haley, S. D. (2019). Genomic selection for winter survival ability among a diverse collection of facultative and winter wheat genotypes. Mol. Breed. 39, 29. doi: 10.1007/s11032-018-0925-8
Bergjord Olsen, A. K., Persson, T., de Wit, A., Nkurunziza, L., Sindhøj, E., Eckersten, H. (2018). Estimating winter survival of winter wheat by simulations of plant frost tolerance. J. Agron. Crop Sci. 204, 62–73. doi: 10.1111/jac.12238
Dahal, N., Xiong, D. H., Neupane, N., Zhang, S., Yuan, Y., Zhang, B. J., et al. (2021). Estimating and analyzing the spatiotemporal characteristics of crop yield loss in response to drought in the Koshi river basin, Nepal. Theor. Appl. Climatol. 152, 1053–1073. doi: 10.1007/s00704-023-04447-8
Ferrante, A., Cossani, C. M., Able, J. A., Sadras, V. O. (2024). Yield response to frost in a set of historic wheat varieties. Field Crops Res. 310, 109336. doi: 10.1016/j.fcr.2024.109336
Fowler, D. B., Byrns, B. M., Greer, K. J. (2014). Overwinter low-temperature responses of cereals: analyses and simulation. Crop Sci. 54, 2395–2405. doi: 10.2135/cropsci2014.03.0196
Gao, F., Li, B., Ren, B., Zhao, B., Liu, P., Zhang, J. W. (2019). Effects of residue management strategies on greenhouse gases and yield under double cropping of winter wheat and summer maize. Sci. Total Environ. 112, 4287–4302. doi: 10.1002/agj2.20283
Gong, S., Zhang, L., Gu, Y. (1982). Simulation study of winter wheat freezing injury during winter. Meteorological 11, 32–34. doi: 10.7519/j.issn.1000-0526.1982.11.014
Gusta, L. V., Wisniewski, M. (2013). Understanding plant cold hardiness: an opinion. Physiologia plantarum 147, 4–14. doi: 10.1111/j.1399-3054.2012.01611.x
Hong, S. Z., Jiao, F. L., Li, M., Zhang, Q. F., Ma, Y. Z., Li, Q. (2023). Effects of delayed irrigation and different planting patterns on carbon emissions and yield of winter wheat in the North China Plain. Irrigation Sci. 41, 365–377. doi: 10.1007/s00271-022-00838-3
Islam, A. T., Islam, A. S., Islam, G. T., Bala, S. K., Salehin, M., Choudhury, A. K., et al. (2022). Adaptation strategies to increase water productivity of wheat under changing climate. Agric. Water Manage. 264, 107499. doi: 10.1016/j.agwat.2022.107499
Jaimez, R. E., Pena, G., Barragan, L., Chica, E., Arteaga, F., Cedeno, G. (2023). The effect of water deficit on leaf stomatal conductance, water relations, chlorophyll fluorescence and growth of rootstock-scion combinations of cacao. Scientia Hortic. 321, 112335. doi: 10.1016/j.scienta.2023.112335
Judy, L., Marjolijn, H., Robert, L. (2020). Climate change: making decisions in the face of deep uncertainty. Nature 580, 456. doi: 10.1038/d41586-020-01147-5
Kinugasa, T., Gantsetseg, B. (2023). Factors causing different CO2 responses in shoot regrowth after defoliation between two grassland species in Mongolia. J. arid environments 210, 1–7. doi: 10.1016/j.jaridenv.2022.104916
Li, H. W., Li, Y. P., Huang, G. H., Sun, J. (2021). Probabilistic assessment of crop yield loss to drought time: cales in Xinjiang, China. Int. J. Climatol. 41, 4077–4094. doi: 10.1002/joc.7059
Li, X., Pu, H., Liu, F., Zhou, Q., Cai, J., Dai, T. B., et al. (2015). Winter wheat photosynthesis and grain yield responses to spring freeze. Agron. J. 107, 1002–1010. doi: 10.2134/agronj14.0460
Licker, R., Kucharik, C. J., Doré, T., Lindeman, M. J., Makowski, D. (2013). Climatic impacts on winter wheat yields in Picardy, France and Rostov, Russia: 1973-2010. Agric. For. Meteorol. 176, 25–37. doi: 10.1016/j.agrformet.2013.02.010
Liu, Y., Wang, E., Yang, X., Wang, J. (2009). Contributions of climatic and crop varietal changes to crop production in the North China Plain, since 1980s. Global Change Biol. 16, 2287–2299. doi: 10.1111/j.1365-2486.2009.02077.x
Luo, Q., O’Leary, G., Cleverly, J., Eamus, D. (2018). Effectiveness of time of sowing and cultivar choice for managing climate change: wheat crop phenology and water use efficiency. Int. J. biometeorol. 62, 1049–1061. doi: 10.1007/s00484-018-1508-4
Macholdt, J., Honermeier, B. (2017). Importance of variety choice: Adapting to climate change in organic and conventional farming systems in Germany. Outlook Agric. 46, 178–184. doi: 10.1177/0030727017722420
Min, K., Chen, K., Arora, R. (2020). Short versus prolonged freezing differentially impacts freeze-thaw injury in spinach leaves: mechanistic insights through metabolite profiling. Physiologia plantarum 168, 777–789. doi: 10.1111/ppl.13033
Mu, C., Yang, X., Yang, J., Li, K. N., Zheng, D. X. (2015). Freezing resistance and injury indices for different cultivars of winter-spring wheat in Huang-Huai-Hai Plain. I. Comparison of freezing resistance for different cultivars of winter-spring wheat during mid-winter period. Chin. J. Appl. Ecol. 26, 3119–3125. doi: 10.13287/j.1001-9332.20150921.011
Persson, T., Bergjord, O. A. K., Nkurunziza, L., Sindhoj, E., Eckersten, H. (2017). Estimation of crown temperature of winter wheat and the effect on simulation of frost tolerance. J. Agron. Crop Sci. 203, 161–176. doi: 10.1111/jac.12187
Pirjo, P. S., Kaija, H., Lauri, J. (2011). Climate-induced overwintering challenges for wheat and rye in northern agriculture. Acta Agriculturae Scandinavica 61, 75–83. doi: 10.1080/09064710903535977
Rapacz, M. (2007). Chlorophyll a fluorescence transient during freezing and recovery in winter wheat. Photosynthetica 45, 409–418. doi: 10.1007/s11099-007-0069-2
Rapacz, M., Ergon, A., Hoglind, M., Jorgensen, M., Jurczyk, B., Ostrem, L., et al. (2014). Overwintering of herbaceous plants in a changing climate. Still more questions than answers. Plant sci.: an Int. J. Exp. Plant Biol. 225, 34–44. doi: 10.1016/j.plantsci.2014.05.009
Rapacz, M., Sasal, M., Wójcik-Jagła, M. (2015). Direct and indirect measurements of freezing tolerance: advantages and limitations. Acta Physiologiae Plantarum 37, 157. doi: 10.1007/s11738-015-1907-7
Rapacz, M., Woźniczka, A. (2009). A selection tool for freezing tolerance in common wheat using the fast chlorophyll a fluorescence transient. Plant Breed. 128, 227–234. doi: 10.1111/j.1439-0523.2008.01566.x
Rizza, F., Pagani, D., Gut, M., Prasil, I. T., Lago, C., Tondelli, A., et al. (2011). Diversity in the response to low temperature in representative barley genotypes cultivated in Europe. Crop Sci. 51, 2759–2779. doi: 10.2135/cropsci2011.01.0005
Short, W., Wolyn, D. (2022). The role of drought stress on the acquisition of freezing tolerance in asparagus under controlled conditions. Can. J. Plant Sci. 102, 864–874. doi: 10.1139/cjps-2021-0266
Skinner, D. Z. (2014). Time and temperature interactions in freezing tolerance of winter wheat. Crop Sci. 54, 1523–1529. doi: 10.2135/cropsci2013.09.0623
Skinner, D. Z., Bellinger, B. S. (2011). Differential response of wheat cultivars to components of the freezing process in saturated soil. Crop Sci. 51, 69–74. doi: 10.2135/cropsci2010.05.0286
Skinner, D. Z., Mackey, B. (2009). Freezing tolerance of winter wheat plants frozen in saturated soil. Field Crops Res. 113, 335–341. doi: 10.1016/j.fcr.2009.06.014
Sommer, U., Lengfellner, K. (2010). Climate change and the timing, magnitude, and composition of the phytoplankton spring bloom. Global Change Biol. 14, 1199–1208. doi: 10.1111/j.1365-2486.2008.01571.x
Tudela, V., Santibáñez, F. (2016). Modelling impact of freezing temperatures on reproductive organs of deciduous fruit trees. Agric. For. Meteorol. 226-227, 28–36. doi: 10.1016/j.agrformet.2016.05.002
Veisz, O., Braun, H. J., Bed, Z. (2001). Plant damage after freezing, and the frost resistance of varieties from the facultative and winter wheat observation nurseries. Euphytica 119, 179–183. doi: 10.1023/A:1017528624476
Vico, G., Hurry, V., Weih, M. (2014). Snowed in for survival: Quantifying the risk of winter damage to overwintering field crops in northern temperate latitudes. Agric. For. Meteorol. 197, 65–75. doi: 10.1016/j.agrformet.2014.06.003
Waalen, W. M., Tanino, K. K., Olsen, J. E., Eltun, R., Rognli, O. A., Gusta, L. V. (2011). Freezing tolerance of winter canola cultivars is best revealed by a prolonged freeze test. Crop Sci. 51, 1988–1996. doi: 10.2135/cropsci2011.02.0098
Wang, H. F., Huo, Z. G., Zhou, G. S., Liao, Q. H., Feng, H. K., Wu, L. (2016). Estimating leaf SPAD values of freeze-damaged winter wheat using continuous wavelet analysis. Plant Physiol. Biochem. 98, 39–45. doi: 10.1016/j.plaphy.2015.10.032
Wang, K., Su, L., Wang, Q. (2021). Cotton growth model under drip irrigation with film mulching: A case study of Xinjiang, China. Agron. J. 113, 2417–2436. doi: 10.1002/agj2.20667
Wang, J., Wang, E., Yang, X., Zhang, F. S., Yin, H. (2012). Increased yield potential of wheat-maize cropping system in the North China Plain by climate change adaptation. Climatic Change 113, 825–840. doi: 10.1007/s10584-011-0385-1
Wu, L., Feng, L., Li, Y., Wang, J., Wu, L. H. (2019). A yield-related agricultural drought index reveals spatio-temporal characteristics of droughts in Southwestern China. Sustainability 11, 1–13. doi: 10.3390/su11030714
Yang, X. M., Reynolds, W. D., Drury, C. F., Reeb, M. D. (2021). Cover crop effects on soil temperature in a clay loam soil in Southwestern Ontario. Can. J. Soil Sci. 101, 761–770. doi: 10.1139/cjss-2021-0070
You, Q., Kang, S., Aguilar, E., Pepin, N., Flügel, W. A., Yan, Y. P., et al. (2010). Changes in daily climate extremes in China and their connection to the large scale atmospheric circulation during 1961-2003. Climate Dynamics 36, 2399–2417. doi: 10.1007/s00382-009-0735-0
Zhang, H. L., Zhao, X., Yin, X. G., Liu, S. L., Xue, J. F., Wang, M., et al. (2015). Challenges and adaptations of farming to climate change in the North China Plain. Climatic Change 129, 213–224. doi: 10.1007/s10584-015-1337-y
Zheng, W. (1981). Mathematical model of winter wheat freezing injury during winter. Chin. J. Agrometeorol. 2, 35–44.
Keywords: freezing tolerance, lower temperature threshold of tiller, yield loss, soil effective negative accumulated temperature, Fv/Fm, winter wheat
Citation: Wu L, Yu W, Cheng C, Feng L, Yan J, Zheng D and Yang F (2024) Integrative indexes reveal the tolerance of winter wheat to different overwinter freezing injury. Front. Plant Sci. 15:1419381. doi: 10.3389/fpls.2024.1419381
Received: 18 April 2024; Accepted: 19 August 2024;
Published: 17 October 2024.
Edited by:
Balpreet Kaur Dhatt, Bayer Crop Science, United StatesReviewed by:
Ildikó Jócsák, Hungarian University of Agricultural and Life Sciences, HungaryMuhammad Ahmad Hassan, Anhui Academy of Agricultural Sciences, China
Copyright © 2024 Wu, Yu, Cheng, Feng, Yan, Zheng and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chen Cheng, chengsir1993@lsu.edu.cn; Liping Feng, fenglp@cau.edu.cn