<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2023.1262001</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Crop sensitivity to waterlogging mediated by soil temperature and growth stage</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Xu</surname>
<given-names>Fu-Li</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2380257"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Hu</surname>
<given-names>Pei-Min</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wan</surname>
<given-names>Xiao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Harrison</surname>
<given-names>Matthew Tom</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Ke</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xiong</surname>
<given-names>Qin-Xue</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2496908"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Agriculture, Yangtze University</institution>, <addr-line>Jingzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Meteorological Service Center, Jingzhou Meteorological Bureau</institution>, <addr-line>Jingzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Tasmanian Institute of Agriculture, University of Tasmania</institution>, <addr-line>Launceston, TAS</addr-line>, <country>Australia</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Najeeb Ullah, Qatar University, Qatar</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Junjun Ni, Southeast University, China; Ujjal Jyoti Phukan, University of Arizona, United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Pei-Min Hu, <email xlink:href="mailto:hupm@21cn.com">hupm@21cn.com</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>26</day>
<month>10</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1262001</elocation-id>
<history>
<date date-type="received">
<day>19</day>
<month>07</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>09</day>
<month>10</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Xu, Hu, Wan, Harrison, Liu and Xiong</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Xu, Hu, Wan, Harrison, Liu and Xiong</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<p>Waterlogging constrains crop yields in many regions around the world. Despite this, key drivers of crop sensitivity to waterlogging have received little attention. Here, we compare the ability of the SWAGMAN Destiny and CERES models in simulating soil aeration index, a variable contemporaneously used to compute three distinct waterlogging indices, denoted hereafter as WI <sub>Destiny</sub>, WI<sub>ASD1</sub>, and WI<sub>ASD2</sub>. We then account for effects of crop growth stage and soil temperature on waterlogging impact by introducing waterlogging severity indices, WI <sub>Growth</sub>, which accommodates growth stage tolerance, and WI <sub>Plus</sub>, which accounts for both soil temperature and growth stage. We evaluate these indices using data collected in pot experiments with genotypes &#x201c;Yang mai 11&#x201d; and &#x201c;Zheng mai 7698&#x201d; that were exposed to both single and double waterlogging events. We found that WI <sub>Plus</sub> exhibited the highest correlation with yield (-0.82 to -0.86) suggesting that waterlogging indices which integrate effects of temperature and growth stage may improve projections of yield penalty elicited by waterlogging. Importantly, WI <sub>Plus</sub> not only allows insight into physiological determinants, but also lends itself to remote computation through satellite imagery. As such, this index holds promise in scalable monitoring and forecasting of crop waterlogging.</p>
</abstract>
<kwd-group>
<kwd>waterlogging</kwd>
<kwd>SWAGMAN Destiny</kwd>
<kwd>hypoxia</kwd>
<kwd>anoxia</kwd>
<kwd>abiotic stress</kwd>
<kwd>soil aeration</kwd>
<kwd>crop physiology</kwd>
<kwd>model</kwd>
</kwd-group>
<counts>
<fig-count count="6"/>
<table-count count="2"/>
<equation-count count="16"/>
<ref-count count="41"/>
<page-count count="13"/>
<word-count count="6665"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Plant Abiotic Stress</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>In theory, optimal plant water balance may be realized through harmony between root water absorption and leaf transpiration (<xref ref-type="bibr" rid="B15">Ibrahim et&#xa0;al., 2018</xref>). Adequate water supply is essential for healthy plant growth, but when the soil becomes saturated, potentially resulting in plant waterlogging stress (<xref ref-type="bibr" rid="B35">Singhal et&#xa0;al., 2022</xref>). Soil waterlogging can manifest in agricultural fields for various reasons, such as excessive rainfall or irrigation, inadequate soil drainage, rising or perched water tables and lateral surface or subsurface flows (<xref ref-type="bibr" rid="B21">Liu et&#xa0;al., 2021</xref>). This can deplete oxygen levels within soil pores, leading to diminished growth, senescence and, in severe cases, crop mortality (<xref ref-type="bibr" rid="B10">Githui et&#xa0;al., 2022</xref>). Waterlogging can also indirectly impact on growth via soil mineral nitrogen, organic matter and carbon, as well as soil microbiota (<xref ref-type="bibr" rid="B27">Phelan et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B30">Rawnsley et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B31">S&#xe1;ndor et&#xa0;al., 2020</xref>).</p>
<p>The middle and lower reaches of the Yangtze River region are among China&#x2019;s primary wheat production zones, but they are also severely afflicted by waterlogging (<xref ref-type="bibr" rid="B37">Wan et&#xa0;al., 2022</xref>). Around 41% of total arable land area in this region grapples with waterlogging, significantly hampering the consistency and stability of crop yields (<xref ref-type="bibr" rid="B14">Hu, 2023</xref>). Investigating the impact of waterlogging, particularly hypoxia stress, on wheat growth and yield is crucial, together with how waterlogging may be amplified or interplay with other crop stresses, such as extreme heat, that may occur later in the growing season (<xref ref-type="bibr" rid="B18">Langworthy et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B12">Harrison, 2021</xref>).</p>
<p>Two approaches are typically employed to assess the extent of crop waterlogging. The first relies on soil groundwater depth, often computed in crop models using algorithms pertaining to cumulative excess groundwater depth (SEW30) (<xref ref-type="bibr" rid="B17">Kanwar et&#xa0;al., 1988</xref>; <xref ref-type="bibr" rid="B28">Qian et&#xa0;al., 2015</xref>) and consecutive suppression days of waterlogging (CSDI) (<xref ref-type="bibr" rid="B7">Evans and Skaggs, 1993</xref>). This method primarily elucidates the impact of groundwater and surface water on crop waterlogging stress, yet does not allow other waterlogging pathways to be captured, such as waterlogging triggered by excessive irrigation (<xref ref-type="bibr" rid="B22">Liu et&#xa0;al., 2020</xref>). Conversely, the second category is grounded upon occurrence of climatic elements, which then allows computation of indices such as the waterlogged day index (<xref ref-type="bibr" rid="B11">Guo et&#xa0;al., 2016</xref>) and wetland day index (<xref ref-type="bibr" rid="B41">Zhang et&#xa0;al., 2019</xref>). Such parameters ignore the multifaceted and integrated effects of soil type, topography and hydrology. Although remote sensing may potentially furnish insights into such factors at the regional scale (<xref ref-type="bibr" rid="B32">Shahpari et&#xa0;al., 2021</xref>), their availability to practitioners at present is limited.</p>
<p>The period from March to April in the middle and lower reaches of the Yangtze River region typically witnesses total precipitation from 300 to 400&#xa0;mm, constituting 30% to 40% of the annual rainfall. This timeframe coincides with critical nutritional and reproductive growth stages of wheat crops, and is often punctuated by frequent flood events (<xref ref-type="bibr" rid="B8">Gao et&#xa0;al., 2020</xref>). Although wheat has three distinct waterlogging stress responses (tolerant, inhibitory, and adaptive phases) (<xref ref-type="bibr" rid="B33">Shaw and Meyer, 2015</xref>), recurring episodes of waterlogging stress can significantly inhibit growth and yields, especially when high-intensity waterlogging coincides with extreme weather fluctuations. This explains why many models perform well with waterlogging stress under experimental conditions, but fall short of accurately predicting outcomes under field conditions (<xref ref-type="bibr" rid="B34">Shaw et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B13">Harrison et&#xa0;al., 2019</xref>).</p>
<p>Broadly, two approaches are employed to simulate waterlogging in crop models. The first calculates the Stress Day Index (SDI) based on groundwater table depth, as exemplified by DRAINMOD (<xref ref-type="bibr" rid="B36">Skaggs, 2018</xref>). However, this approach overlooks genotypic differences in crops, considering only groundwater and surface water, rendering it unsuitable for addressing waterlogging stemming from overirrigation. The second method computes hypoxic stress as a function of soil moisture. Notable examples include CROPR (<xref ref-type="bibr" rid="B29">Qian et&#xa0;al., 2017</xref>), SWAGMAN Destiny (<xref ref-type="bibr" rid="B26">Meyer et&#xa0;al., 1996</xref>), and the Agricultural Production Systems Simulator (APSIM) (<xref ref-type="bibr" rid="B2">Asseng et&#xa0;al., 1998</xref>; <xref ref-type="bibr" rid="B21">Liu et&#xa0;al., 2021</xref>). <xref ref-type="bibr" rid="B33">Shaw and Meyer (2015)</xref> showed that yield reductions were 6% higher at 65% air-filled pore space compared to 10%, leading to their proposal of three stages for simulating plant responses to waterlogging damage. Their hypoxic stress factor accounts for waterlogging duration and crop tolerance, making it applicable to damage caused by excessive irrigation and facilitating more precise quantitative damage analysis (<xref ref-type="bibr" rid="B34">Shaw et&#xa0;al., 2013</xref>). The Crop Estimation through Resource and Environment Synthesis (CERES) model, commonly used for predicting crop growth and yield, has been enhanced by Lizaso et&#xa0;al., who successfully integrated low oxygen (anaerobic) stress factors into the CERES-Wheat model (<xref ref-type="bibr" rid="B24">Lizaso and Ritchie, 1997</xref>).</p>
<p>The objectives of this study were thus to quantify waterlogging by integrating computations from both the SWAGMAN Destiny and CERES models. Specifically, we develop five wheat waterlogging injured indices: WI<sub>ASD2</sub> and WI<sub>ASD1</sub> were conceptualized as adaptations of the CERES model (<xref ref-type="bibr" rid="B24">Lizaso and Ritchie, 1997</xref>), WI <sub>Destiny</sub> was derived from the SWAGMAN Destiny model, WI <sub>Growth</sub> introduces differences in waterlogging tolerance across crop growth stages (following WI <sub>Destiny</sub>), while WI <sub>Plus</sub> considers the influence of soil temperature in addition to WI <sub>Growth</sub>. To assess the accuracy of these indicators, we conducted pot experiments with winter wheat (involving single or double flooding) using local genotypes (&#x201c;Yang mai 11&#x201d; and &#x201c;Zheng mai 7698&#x201d;). By analyzing relationships between our WIs and yield, we provide insight into how integrated physiological determinants impact on crop yield.</p>
</sec>
<sec id="s2" sec-type="material|methods">
<label>2</label>
<title>Material and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Plant material and experimental design</title>
<p>The study region is located in the Jianghan Plain, Hubei Province, China, characterized by a subtropical monsoon climate. The region&#x2019;s primary crops include rice, wheat, rapeseed, and cotton. Two wheat varieties were employed in the experiments: the waterlog-resistant Yang mai 11 (referred to as Y) and the waterlog-unresistant Zheng mai 7698 (referred to as Z) (<xref ref-type="bibr" rid="B40">Yu et&#xa0;al., 2014</xref>). The soil used in the pot experiment was selected from paddy soil, which is a typical cultivated soil in the middle and lower reaches of the Yangtze River (<xref ref-type="bibr" rid="B8">Gao et&#xa0;al., 2020</xref>).</p>
<p>The winter wheat pot experiment with waterlogging spanned from November 2020 to May 2022 and was conducted at the open field experimental site of the College of Agriculture, Yangtze University (longitude 112&#xb0;08&#x2019;, latitude 30&#xb0;21&#x2019;). This two-year experiment utilized storage boxes measuring 60&#xa0;cm in length, 45&#xa0;cm in width, and 35&#xa0;cm in height, arranged in a randomized complete block design.</p>
<p>The soil&#x2019;s basic physicochemical properties were as follows: pH value of 7.85, total nitrogen content of 1.19 g/kg, total phosphorus content of 0.77 g/kg, total potassium content of 10.45 g/kg, organic matter content of 15.84 g/kg, alkaline nitrogen content of 58.44 mg/kg, effective phosphorus content of 31.12 mg/kg, and available potassium content of 106.09 mg/kg. Wheat was sown on November 3, 2020, with each box receiving a mixture of 17.78g/kg of compound fertilizer (N: P<sub>2</sub>O<sub>5</sub>:K<sub>2</sub>O=15:15:15), 0.89g/kg of KCl, and 3.94g/kg of urea before sowing. At the tillering stage on January 25, 2022, urea was applied at a rate of 4.17 g/kg.</p>
<p>During the experiment, single waterlogging treatments lasting for 5 days, 12 days, and 20 days were administered on March 21 and April 8 in 2021 (or 2022). Following a 10-day interval, a second round of waterlogging treatments for the same durations (i.e., double waterlogging treatments) was applied to the potted plants in multiple boxes. Waterlogging was manually induced by irrigating with marked buckets to maintain soil moisture content at over 90% of the maximum field water-holding capacity. The criterion for waterlogging stress was maintaining a water level 5 centimeters above the soil surface. Daily irrigation volumes (in liters) for each box were meticulously recorded.</p>
<p>The experiment adopted a randomized block design involving continuous irrigation treatments during two critical periods: from tillering to flowering (March 21st to April 7th, denoted as B) and from flowering to maturity (April 8th to May 8th, denoted as Y), which corresponds to the crucial water requirement period for wheat. The experiment included three different durations of irrigation (5 days, 12 days, 20 days) and two types of irrigation treatments (single continuous waterlogging and double continuous waterlogging), resulting in a total of 72 treatments. Each treatment was replicated three times, with an additional six control boxes. Protective rows were established around the experiment using boxes without irrigation, featuring either the &#x201c;Yang mai 11&#x201d; variety or the &#x201c;Zheng mai 7698&#x201d; variety. The control boxes were positioned at the center. In total, the experiment encompassed 120 boxes.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Measurement indexes and methods</title>
<p>Soil Moisture: soil volumetric water content was assessed using an EM50 soil moisture measurement device connected to an EC-5 soil moisture sensor. The water content probe was inserted into the soil at a depth of 5-6cm. Observations were conducted daily from 8:00 to 20:00, with measurements taken at hourly intervals. The daily soil volumetric moisture content was calculated as the average of the 24-hour observation period.</p>
<p>Soil Temperature: Temperature data recorded using a fiber optic oxygen meter (PreSens Microx4, Germany Regensburg) with the measuring part inserted 5-6 cm into the soil. Then take the average as the daily average soil temperature.</p>
<p>Yield Measurement: At maturity (May 8), winter wheat from all boxes was harvested. The harvested wheat was threshed, sun-dried, and subsequently winnowed to remove impurities and empty grains. The total weight of the dried grains was then measured.</p>
<p>Meteorological Elements: The experimental station was equipped with a HOBO automatic weather station (HOBO Micro-Weather Station) that automatically recorded daily rainfall and other meteorological data.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Calculation method of wheat waterlogging injured index</title>
<p>
<italic>WI<sub>Destiny</sub>
</italic> takes into account the combined effects of waterlogging duration and soil moisture on the root system. <italic>WI<sub>ASD</sub>
</italic>
<sub>1</sub> and <italic>WI<sub>ASD</sub>
</italic>
<sub>2</sub> are improvements made by Lizaso, J. I. et&#xa0;al. to the CERES model, where the growth stage is divided into 1-9 phases (<xref ref-type="bibr" rid="B24">Lizaso and Ritchie, 1997</xref>). However, this method has only been validated in research on the impact of waterlogging on maize growth and development, and there is no relevant verification for wheat.</p>
<p>Crop sensitivity to waterlogging stress is heavily dependent on the developmental stage at which waterlogging occurs (<xref ref-type="bibr" rid="B23">Liu et&#xa0;al., 2023</xref>). As such, considering the lack of research on waterlogging tolerance during different growth stages in <italic>WI<sub>Destiny</sub>
</italic>, we have made modifications by using a sigmoid function to simulate the differences in waterlogging tolerance. This has led us to construct a new index called <italic>WI</italic>
<sub>Growth</sub>.</p>
<p>Constant critical values (0.65) have frequently been applied to represent oxygen stress (<xref ref-type="bibr" rid="B3">Barber et&#xa0;al., 2004</xref>; <xref ref-type="bibr" rid="B19">Le&#xe3;o et&#xa0;al., 2006</xref>), but constants are unlikely to be sufficient for any proxy. In the SWAGMAN Destiny model, the critical values for oxygen stress are likely to be inappropriate and must change with different soil types, varied temperature, organic matter content, soil depth and plant characteristics. Among them, the main factor is soil temperature (<xref ref-type="bibr" rid="B4">Bartholomeus et&#xa0;al., 2008</xref>). Therefore, based on <italic>WI</italic>
<sub>Growth</sub>, we have incorporated soil temperature and made improvements to the critical water filled pore space (CritWFPS=0.65), resulting in a revised index called <italic>WI<sub>Plus</sub>
</italic>. The improvement of the SWAGMAN Destiny model is shown in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Sequence diagram of the components required to compute WI Destiny in SWAGMAN destiny (Old WI) and its refinement (Improved WI). WI <sub>Growth</sub> primarily involved the addition of the weight coefficient for the growth period (Coef<sub>i</sub>). WI <sub>Plus</sub> incorporated the Soil Temperature (T<sub>soil</sub>) and the Weight coefficient of growth period (Coef<sub>i</sub>).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1262001-g001.tif"/>
</fig>
<sec id="s2_3_1">
<label>2.3.1</label>
<title>Calculating waterlogging injured index from SWAGMAN Destiny model (<italic>WI<sub>Destiny</sub>
</italic>)</title>
<p>In the SWAGMAN Destiny model, the calculation process for total influence factors of soil aeration index on roots involves several steps (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). Firstly, it is calculated the soil pore water content (<italic>WFPS</italic>) by obtaining the daily water content of the soil surfaces (<italic>S<sub>W</sub>
</italic>). Then, the characteristic parameter (<italic>A<sub>erf</sub>
</italic>) of the effect of low oxygen on the root system was calculated. Finally, combining the characteristic parameters (<italic>A<sub>erf</sub>
</italic>) with the duration of continuous waterlogging (<italic>D<sub>time</sub>
</italic>) to calculate the overall impact factor of low oxygen on root systems (<italic>L<sub>af</sub>
</italic>), providing a comprehensive assessment of the degree of impact of low oxygen stress on roots. The calculation process is as follows:</p>
<p>(1) Characteristic Parameters of hypoxia&#x2019;s impact on root function (<italic>A<sub>erf</sub>
</italic>)</p>
<p>The first step is to calculate the relative amount of pore space filled with water (Aerf) within the soil. This is a zero to unity factor where zero indicates the soil pores are completely full of water (the soil is saturated) and one indicates there is no water in the soil pores (the soil is oven dry). Aerf is calculated (<xref ref-type="bibr" rid="B34">Shaw et&#xa0;al., 2013</xref>) as:</p>
<disp-formula>
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable equalrows="true" equalcolumns="true">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>&gt;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where</p>
<disp-formula>
<label>(2)</label>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>WFPS</italic> is water filled pore space at soil water content; <italic>S<sub>W</sub>
</italic> is Volumetric soil water content (cm<sup>3</sup>
<sub>water</sub> cm<sup>-3</sup>
<sub>soil</sub>); <italic>B<sub>D</sub>
</italic> is Bulk density of dry soil (g/cm<sup>-3</sup>), and in accordance with the actual soil bulk density measured in pot experiments, it is taken as 1.4 g&#xb7;cm<sup>-3</sup>; The general soil particle density falls within the range of 2.6 to 2.8 g/cm<sup>3</sup>, and here it is taken as 2.7 g/cm<sup>3</sup>.</p>
<p>The critical water filled pore space (<italic>CritWFPS</italic>) is set to 0.65. Doran et&#xa0;al. found that respiratory activity of microbial function in soils decreased when the water filled pore space increased to a value above 0.65 (<xref ref-type="bibr" rid="B6">Doran et al., 1990</xref>).</p>
<p>(2) Number of days waterlogged (<italic>D<sub>time</sub>
</italic>)</p>
<p>When the average daily soil pore water content (<italic>CritWFPS</italic>) is greater than 0.65, it is defined that wheat is affected by waterlogging on that day. If the average daily soil pore water content (<italic>CritWFPS</italic>) is less than 0.65, then wheat is in a normal state with a continuous waterlogging duration of 0 days. Due to the delayed response of the root system to waterlogging damage, it is set that the damage will only affect the crop&#x2019;s root system after 3 days, and the impact remains unchanged after 60 days (<xref ref-type="bibr" rid="B33">Shaw and Meyer, 2015</xref>). The calculation formula is as follows (<xref ref-type="bibr" rid="B14">Hu, 2023</xref>):</p>
<disp-formula>
<label>(3)</label>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mo>{</mml:mo>
<mml:mtable equalrows="true" equalcolumns="true">
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0.65</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0.65</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0.65</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtext>&#x2004;&#x2004;&#x2004;&#x2004;&#x2004;&#x2004;&#x2004;</mml:mtext>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&lt;</mml:mo>
<mml:mn>0.65</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0.65</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>60</mml:mn>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>60</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0.65</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <italic>D<sub>time,i</sub>
</italic> represents the duration of waterlogging on the i-th day; <italic>WFPS<sub>,i</sub>
</italic> is the WFPS of the i-th day;</p>
<p>(3) Overall impact factor of hypoxia on root function (<italic>L<sub>af</sub>
</italic>)</p>
<p>The overall impact function of waterlogging on the root system is a comprehensive consideration of the effects of soil moisture and duration of waterlogging. The <italic>L<sub>af</sub>
</italic> value is a characteristic factor ranging from 0 to 1.&#xa0;A value closer to 1 indicates a greater impact of waterlogging on the root system of winter wheat, while a value closer to 0 indicates a lesser impact of waterlogging on the root system of winter wheat. its calculation formula is represented as (<xref ref-type="bibr" rid="B34">Shaw et&#xa0;al., 2013</xref>):</p>
<disp-formula>
<label>(4)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mfenced>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mo>&#x0020;</mml:mo>
<mml:mn>0.167</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>(4) Calculation of Waterlogging Index (<italic>WI<sub>Destiny</sub>
</italic>).</p>
<p>In the SWAGMAN Destiny model, <italic>L<sub>af</sub>
</italic> represents the simulated daily impact of waterlogging on winter wheat. Considering the small magnitude of values, it is multiplied by 1000 for scaling purposes. The daily impact function of wheat is averaged throughout the entire growing season to calculate the Waterlogging Index (<italic>WI<sub>Destiny</sub>
</italic>) for the entire crop cycle.</p>
<disp-formula>
<label>(5)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>n</mml:mi>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
<mml:mfenced>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>1000</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</sec>
<sec id="s2_3_2">
<label>2.3.2</label>
<title>Calculating waterlogging injured index from CERES model (<italic>WI<sub>ASD</sub>
</italic>
<sub>1</sub> and <italic>WI<sub>ASD</sub>
</italic>
<sub>2</sub>)</title>
<p>Jones et&#xa0;al. proposed a model to predict root growth that includes soil constraints affecting root growth (<xref ref-type="bibr" rid="B16">Jones et al., 1991</xref>). The model uses the fraction of water filled pore space as an aeration index to affect, root attributes such as rooting depth, branching and senescence.</p>
<p>(1) Soil aeration index</p>
<p>Aeration indexes are calculated in terms of soil porosity since oxygen diffusion depends on the air-filled pore space. According to formula (2), the water-filled porosity (WFPS) is calculated.</p>
<p>The water-filled pore space at saturation is assumed to have a maximum value (<italic>XWWFPS</italic>) of 0.93 to account for trapped air in the soil profile. Whenever the water-filled pore space (WFPS<sub>i</sub>) in a soil layer is larger than a critical value (<italic>CWFPS</italic>), a layer aeration factor (<italic>LAF<sub>i</sub>
</italic>) was calculated as (<xref ref-type="bibr" rid="B24">Lizaso and Ritchie, 1997</xref>):</p>
<disp-formula>
<label>(6)</label>
<mml:math display="block" id="M6">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mtext>LAF</mml:mtext>
</mml:mrow>
<mml:mtext>i</mml:mtext>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfenced>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>WFPS</mml:mtext>
</mml:mrow>
<mml:mtext>i</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo stretchy="false">/</mml:mo>
<mml:mfenced>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>CWFPS</mml:mtext>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>CWFPS</italic> is equal to 0.45 and <italic>LAF<sub>i</sub>
</italic> cannot be larger than 1. When a layer with <italic>WFPS</italic> of 0.9 or larger limits the aeration of lower layers. The Whole Rhizosphere Aeration Factor (WRAF) for the entire root zone is calculated by integrating all soil layers using root density as the weighting factor. The formula is as follows (<xref ref-type="bibr" rid="B24">Lizaso and Ritchie, 1997</xref>):</p>
<disp-formula>
<label>(7)</label>
<mml:math display="block" id="M7">
<mml:mrow>
<mml:mtext>WRAF</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mtext>i</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mtext>n</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
<mml:msub>
<mml:mrow>
<mml:mtext>LAF</mml:mtext>
</mml:mrow>
<mml:mtext>i</mml:mtext>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mo>*</mml:mo>
<mml:mtext>RLV</mml:mtext>
</mml:mrow>
<mml:mtext>i</mml:mtext>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo stretchy="false">/</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mtext>i</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mtext>n</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
<mml:msub>
<mml:mrow>
<mml:mtext>RLV</mml:mtext>
</mml:mrow>
<mml:mtext>i</mml:mtext>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where RKV<sub>i</sub> is the root length density for each soil layer. The soil in the root layer is divided into <italic>n</italic> layers. <italic>WRAF</italic> is a dimensionless unit ranging from 0 to 1, where a value closer to 1 indicates good soil aeration, while a value closer to 0 indicates severe root oxygen stress. Two dimensionless aeration stress indexes, ASD and ASD2, were calculated using WRAF to simulate the cumulative aeration status experienced by the plant.</p>
<p>The two dimensionless aeration stress indexes (ASD1 and ASD2) are calculated using the following formulas (<xref ref-type="bibr" rid="B24">Lizaso and Ritchie, 1997</xref>):</p>
<disp-formula>
<label>(8)</label>
<mml:math display="block" id="M8">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mo>{</mml:mo>
<mml:mtable equalrows="true" equalcolumns="true">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>WRAF</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>WRAF</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>WRAF</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>WRAF</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>0.855</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mtext>ISTAGE</mml:mtext>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
</disp-formula>
<p>
<inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mtext>ISTAGE</mml:mtext>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is an integer variable (1 to 9) identifying the phenological stage of the crop.</p>
<disp-formula>
<label>(9)</label>
<mml:math display="block" id="M9">
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mo>{</mml:mo>
<mml:mtable equalrows="true" equalcolumns="true">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>WRAF</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mtext>WRAF</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mtext>ASD</mml:mtext>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mtext>day</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
</disp-formula>
<p>
<italic>k</italic> is the coefficient of growth stage. During the vegetative growth stage,<italic>k</italic> =0.8 when <inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> &gt; 4.0 and <italic>k</italic>=0.1 when <inline-formula>
<mml:math display="inline" id="im3">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>&lt;4.0. During the remainder of the season, <italic>k</italic>=0.5.</p>
<p>(2) Calculation of <inline-formula>
<mml:math display="inline" id="im4">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im5">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>The average ASD value for each day of the entire growth season of winter wheat is calculated to obtain the waterlogging index of the whole reproductive period. The calculation formula is as follows:</p>
<disp-formula>
<label>(10)</label>
<mml:math display="block" id="M10">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>n</mml:mi>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:msub>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(11)</label>
<mml:math display="block" id="M11">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>n</mml:mi>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
<mml:mi>A</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>D</mml:mi>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
</sec>
<sec id="s2_3_3">
<label>2.3.3</label>
<title>Calculation of waterlogging injured index for growth stage (<italic>WI</italic>
<sub>Growth</sub>)</title>
<p>Modify Equation (4) by introducing the impact function of the growth stage. The calculation process is as follows:</p>
<disp-formula>
<label>(12)</label>
<mml:math display="block" id="M12">
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>L</mml:mi>
<mml:mi>a</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mo>[</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mo>&#x0020;</mml:mo>
<mml:mn>0.167</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msup>
<mml:mo>]</mml:mo>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(13)</label>
<mml:math display="block" id="M13">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:mtext>Growth</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>n</mml:mi>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>G</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>1000</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <inline-formula>
<mml:math display="inline" id="im6">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> refers to formula (1). <inline-formula>
<mml:math display="inline" id="im7">
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> refers to formula (3). <inline-formula>
<mml:math display="inline" id="im8">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the reaction weight coefficient of winter wheat to waterlogging during different growth stages, which ranges from 0 to 1. The impact of winter flooding on winter wheat is generally not significant during the overwintering period, with a small <inline-formula>
<mml:math display="inline" id="im9">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> value. The coefficient increases during the vegetative growth period, and reaches its maximum during the reproductive growth period, following an &#x201c;S&#x201d;-shaped curve pattern (<xref ref-type="bibr" rid="B33">Shaw and Meyer, 2015</xref>). Therefore, to simulate this characteristic using the sigmoid function, the formula is as follows (<xref ref-type="bibr" rid="B14">Hu, 2023</xref>):</p>
<disp-formula>
<label>(14)</label>
<mml:math display="block" id="M14">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mfenced>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.06</mml:mn>
<mml:mo>*</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>5.0</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mtext>&#xa0;</mml:mtext>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <italic>i</italic> represents the number of days from November 30th of the previous year. During the period from winter wheat sowing to November 30th, winter wheat is in its seedling stage, and it is assumed that waterlogging has no impact on its growth, therefore the <inline-formula>
<mml:math display="inline" id="im10">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> value is 0. The curve showing the variation of <inline-formula>
<mml:math display="inline" id="im11">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> with date can be obtained by calculating according to equation (14), as shown in <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>The weighting coefficient (Coef) for winter wheat waterlogging damage.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1262001-g002.tif"/>
</fig>
</sec>
<sec id="s2_3_4">
<label>2.3.4</label>
<title>New waterlogging injured index (<italic>WI<sub>Plus</sub>
</italic>)</title>
<p>We modify the value of CritWFPS in Eq. (1), keeping the other processes unchanged, and calculate the new waterlogging injured index (<italic>WI<sub>Plus</sub>
</italic>).</p>
<p>The format of the relative increase ratio D<sub>soil</sub> of root respiration at soil temperature T<sub>soil</sub> (&#xb0;C) is following (<xref ref-type="bibr" rid="B1">Amthor, 2000</xref>):</p>
<disp-formula>
<label>(15)</label>
<mml:math display="block" id="M15">
<mml:mrow>
<mml:msub>
<mml:mtext>D</mml:mtext>
<mml:mrow>
<mml:mtext>soil</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msubsup>
<mml:mtext>D</mml:mtext>
<mml:mn>10</mml:mn>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mtext>T</mml:mtext>
<mml:mrow>
<mml:mtext>soil</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>10</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>10</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>D<sub>10</sub> is the relative increase ratio at a temperature increase of 10&#xb0;C, equals 2.0 (<xref ref-type="bibr" rid="B25">Lloyd and Taylor, 1994</xref>).</p>
<p>According to the above formula, the critical value of the water filled pore space (<inline-formula>
<mml:math display="inline" id="im12">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>) at which <inline-formula>
<mml:math display="inline" id="im13">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> = 0.65 can be modified as follows:</p>
<disp-formula>
<label>(16)</label>
<mml:math display="block" id="M16">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mtext>soil</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>*</mml:mo>
<mml:mtext>&#x2004;</mml:mtext>
<mml:mo>(</mml:mo>
<mml:mn>0.8</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>0.2</mml:mn>
<mml:mo stretchy="false">/</mml:mo>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>After replacing <inline-formula>
<mml:math display="inline" id="im14">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> in equation (1) with <inline-formula>
<mml:math display="inline" id="im15">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mtext>soil</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and combining it with Eq. (12) to calculate <inline-formula>
<mml:math display="inline" id="im16">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>P</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results and analysis</title>
<sec id="s3_1">
<label>3.1</label>
<title>The impact of waterlogging duration on the rate of increase or decrease in wheat yield</title>
<p>The study analyzed the impact of single or double continuous waterlogging during two different growth stages (jointing stage to heading stage, B period; or heading stage to maturity stage, Y period) on the yield of two wheat varieties, &#x201c;Yang mai 11&#x201d; and &#x201c;Zheng mai 7698&#x201d;. The duration of waterlogging was calculated using formula (3), and the relationship between waterlogging duration and wheat yield rate is shown in <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>. It can be observed from the figure that there is a highly significant correlation (p&lt;0.01) between the number of days of waterlogging and the wheat yield rate. The correlation coefficients were -0.85 (&#x201c;Yang mai 11&#x201d;, single continuous waterlogging), -0.79 (&#x201c;Yang mai 11&#x201d;, double continuous waterlogging), -0.93 (&#x201c;Zheng mai 7698&#x201d;, single continuous waterlogging), and -0.67 (&#x201c;Zheng mai 7698&#x201d;, double continuous waterlogging). This means that the longer the duration of waterlogging, the lower the wheat yield rate.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>The relationship between the length of waterlogging and the rate of increase or decrease in wheat yield under different waterlogging treatments [ &#x201c;a.Yang mai 11&#x201d; and &#x201c;b.Zheng mai 7698&#x201d;]. **represented extremely significant correlation (<italic>P</italic>&lt;0.01).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1262001-g003.tif"/>
</fig>
<p>The linear simulation curves of the rate of increase or decrease in yield in relation to the duration of waterlogging show significant differences between single continuous waterlogging and double continuous waterlogging treatments. The slope of the curve for the single continuous waterlogging treatment is greater than that of the double continuous waterlogging treatment, indicating that the rate of yield increase/decrease of wheat is also influenced by the number of waterlogging events. When the duration of waterlogging is the same, the rate of yield increase/decrease of wheat during the Y period is significantly lower than that during the B period. The impact of single continuous waterlogging treatment on the rate of increase or decrease in wheat yield during the Y period is more pronounced for both &#x201c;Yang mai 11&#x201d; and &#x201c;Zheng mai 7698&#x201d; varieties compared to the double continuous waterlogging treatment.</p>
<p>Although there is a highly significant negative linear correlation between the duration of waterlogging and the rate of increase or decrease in wheat yield, using only the duration of waterlogging is not sufficient to represent the extent of crop waterlogging. The parameters of the linear simulation curve are evidently related to the duration of continuous waterlogging, the number of waterlogging events, and the varieties affected by waterlogging.</p>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>The relationship between five WI indices and wheat yield</title>
<p>Under all treatments,the relationship between the five Wl indices and the relative yield of two wheat varieties is shown in <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref> (&#x201c;Yang mai 11&#x201d;) and <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref> (&#x201c;Zheng mai 7698&#x201d;). From the graph, we can observe that the fitted curves between the WI indices and relative yield of wheat exhibit a clear upward trend initially, reaching a certain point before showing a downward trend. We refer to the ascending stage as the first stage and the descending stage as the second stage, with the WI value at the inflection point referred to as the threshold.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Fitted curves of five WI indices with the relative yields of &#x201c;Yang mai 11&#x201d;.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1262001-g004.tif"/>
</fig>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Fitted curves of five WI indices with the relative yields of &#x201c;Zheng mai 7698&#x201d;.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1262001-g005.tif"/>
</fig>
<p>The fitted curves of the five WI indices with wheat yield can reflect the tolerance of wheat varieties to waterlogging. The variety &#x201c;Yang mai 11&#x201d; exhibits strong tolerance to waterlogging, while &#x201c;Zheng mai 7698&#x201d; is not tolerant. In the second stage, the fitted curves between the five WI indices and the yield of &#x201c;Zheng mai 7698&#x201d; show a steeper slope compared to the fitted curves between the same indices and the yield of &#x201c;Yang mai 11&#x201d;. This means that under the same intensity of waterlogging stress, the decrease in yield for &#x201c;Zheng mai 7698&#x201d; is significantly greater than that for &#x201c;Yang mai 11&#x201d;.</p>
<p>(1) In the First Stage</p>
<p>From <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>, it can be seen that in both wheat varieties, the five indices show a positive correlation with the relative yield of wheat, indicating that larger WI values correspond to higher yields. The order of correlation coefficients is as follows: 0.71&gt; 0.70&gt; 0.66&gt; 0.53&gt; 0.50, which corresponds to WI <sub>Plus</sub> &gt; WI <sub>Growth</sub>&gt; WI <sub>Destiny</sub>&gt; WI<sub>ASD2</sub>&gt; WI<sub>ASD1</sub> (&#x201c;Yang mai11&#x201d;); 0.88&gt; 0.87&gt; 0.86&gt; 0.63&gt; 0.55, which corresponds to WI <sub>Growth</sub>&gt;WI <sub>Destiny</sub> &gt;WI <sub>Plus</sub> &gt;WI<sub>ASD2</sub>&gt;WI<sub>ASD1</sub>(&#x201c;Zheng mai 7698&#x201d;). Additionally, the results of fitting WI Growth, WI Destiny, and WI Plus with the yield of &#x201c;Zheng mai 7698&#x201d; all reach a highly significant level (p&lt; 0.01), with respective &#x2018;r&#x2019; values of 0.88, 0.87, and 0.86. WI Growth, WI Plus, and WI Destiny show significant fits with the yield of &#x201c;Yang mai 11&#x201d; (p&lt; 0.05) with respective &#x2018;r&#x2019; values of 0.70, 0.66, and 0.66.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>The simulation of the five WI and the relative yields of Wheat in the First Stage.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">variety</th>
<th valign="middle" align="center">Five Waterlogging Index</th>
<th valign="middle" align="center">Simulation curves</th>
<th valign="middle" align="center">
<italic>R</italic>&#xb2;</th>
<th valign="middle" align="center">
<italic>r</italic>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" rowspan="5" align="center">Yang mai 11</td>
<td valign="middle" align="center">WI <sub>Plus</sub>
</td>
<td valign="middle" align="center">y=0.0274x+0.9947</td>
<td valign="middle" align="center">0.51</td>
<td valign="middle" align="center">0.71*</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Growth</sub>
</td>
<td valign="middle" align="center">y=0.0455x+1.0334</td>
<td valign="middle" align="center">0.48</td>
<td valign="middle" align="center">0.70*</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Destiny</sub>
</td>
<td valign="middle" align="center">y=0.0392x+1.043</td>
<td valign="middle" align="center">0.44</td>
<td valign="middle" align="center">0.66*</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD2</sub>
</td>
<td valign="middle" align="center">y=0.021x+0.8548</td>
<td valign="middle" align="center">0.29</td>
<td valign="middle" align="center">0.53</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD1</sub>
</td>
<td valign="middle" align="center">y=0.0341x+0.9077</td>
<td valign="middle" align="center">0.25</td>
<td valign="middle" align="center">0.50</td>
</tr>
<tr>
<td valign="middle" rowspan="5" align="center">Zheng mai 7698</td>
<td valign="middle" align="center">WI <sub>Growth</sub>
</td>
<td valign="middle" align="center">y=0.0388x+0.9879</td>
<td valign="middle" align="center">0.77</td>
<td valign="middle" align="center">0.88**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Destiny</sub>
</td>
<td valign="middle" align="center">y=0.0361x+0.9913</td>
<td valign="middle" align="center">0.75</td>
<td valign="middle" align="center">0.87**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Plus</sub>
</td>
<td valign="middle" align="center">y=0.0204x+0.954</td>
<td valign="middle" align="center">0.73</td>
<td valign="middle" align="center">0.86**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD2</sub>
</td>
<td valign="middle" align="center">y=0.0267x+0.7294</td>
<td valign="middle" align="center">0.39</td>
<td valign="middle" align="center">0.63</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD1</sub>
</td>
<td valign="middle" align="center">y=0.0361x+0.8468</td>
<td valign="middle" align="center">0.31</td>
<td valign="middle" align="center">0.55</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>**represented extremely significant correlation (P&lt;0.01). *represented significant correlation (P&lt;0.05). The same as follow.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>In the first stage, the effect of the three indexes of WI <sub>Plus</sub>, WI <sub>Growth</sub> and WI <sub>Destiny</sub> were better than that of WI<sub>ASD1</sub> and WI<sub>ASD2</sub>.</p>
<p>(2) In the Second Stage</p>
<p>According to <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>, the fitting results of the five WI indices with the relative yield of both wheat varieties reach a highly significant level of 0.01.</p>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>The simulation of the five WI and the relative yields of Wheat in the Second Stage.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">variety</th>
<th valign="middle" align="center">Five Waterlogging Index</th>
<th valign="middle" align="center">Simulation curves</th>
<th valign="middle" align="center">
<italic>R</italic>&#xb2;</th>
<th valign="middle" align="center">
<italic>r</italic>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" rowspan="5" align="center">Yang mai 11</td>
<td valign="middle" align="center">WI <sub>Plus</sub>
</td>
<td valign="middle" align="center">y=-0.0035x+1.1863</td>
<td valign="middle" align="center">0.67</td>
<td valign="middle" align="center">-0.82**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Growth</sub>
</td>
<td valign="middle" align="center">y=-0.0094x+1.1814</td>
<td valign="middle" align="center">0.59</td>
<td valign="middle" align="center">-0.77**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Destiny</sub>
</td>
<td valign="middle" align="center">y=-0.0081x+1.1547</td>
<td valign="middle" align="center">0.53</td>
<td valign="middle" align="center">-0.73**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD2</sub>
</td>
<td valign="middle" align="center">y=-0.0024x+1.1046</td>
<td valign="middle" align="center">0.31</td>
<td valign="middle" align="center">-0.56**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD1</sub>
</td>
<td valign="middle" align="center">y=-0.0024x+1.0663</td>
<td valign="middle" align="center">0.27</td>
<td valign="middle" align="center">-0.52**</td>
</tr>
<tr>
<td valign="middle" rowspan="5" align="center">Zheng mai 7698</td>
<td valign="middle" align="center">WI <sub>Plus</sub>
</td>
<td valign="middle" align="center">y=-0.0053x+1.3051</td>
<td valign="middle" align="center">0.74</td>
<td valign="middle" align="center">-0.86**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD2</sub>
</td>
<td valign="middle" align="center">y=-0.0073x+1.3356</td>
<td valign="middle" align="center">0.72</td>
<td valign="middle" align="center">-0.85**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>ASD1</sub>
</td>
<td valign="middle" align="center">y=-0.0081x+1.2543</td>
<td valign="middle" align="center">0.68</td>
<td valign="middle" align="center">-0.82**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Growth</sub>
</td>
<td valign="middle" align="center">y=-0.0141x+1.3165</td>
<td valign="middle" align="center">0.53</td>
<td valign="middle" align="center">-0.73**</td>
</tr>
<tr>
<td valign="middle" align="center">WI <sub>Destiny</sub>
</td>
<td valign="middle" align="center">y=-0.015x+1.3356</td>
<td valign="middle" align="center">0.52</td>
<td valign="middle" align="center">-0.72**</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>**represented extremely significant correlation(P&lt;0.01). *represented significant correlation(P&lt;0.05). The same as follow.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>The five indices show a negative correlation with the relative yield of both wheat varieties, indicating that larger WI values are associated with more severe waterlogging and lower wheat yields. The order of the absolute correlation coefficients (|r|) is as follows: 0.82 &gt; 0.77 &gt; 0.73 &gt; 0.56 &gt; 0.52, corresponding to WI <sub>Plus</sub> &gt; WI <sub>Growth</sub> &gt; WI <sub>Destiny</sub> &gt;WI<sub>ASD2</sub>&gt;WI<sub>ASD1</sub> (&#x201c;Yang mai 11&#x201d;); 0.86 &gt; 0.85 &gt; 0.82 &gt; 0.73 &gt; 0.72, corresponding to WI <sub>Plus</sub> &gt;WI<sub>ASD2</sub>&gt;WI<sub>ASD1</sub>&gt; WI <sub>Growth</sub> &gt; WI <sub>Destiny</sub> (&#x201c;Zheng mai 7698&#x201d;). WI <sub>Plus</sub> reflects the waterlogging severity of wheat better, and it has the highest R and |r| values, regardless of whether it is for &#x201c;Zheng mai 7698&#x201d; or &#x201c;Yang mai 11&#x201d;.</p>
<p>WI<sub>ASD1</sub> and WI<sub>ASD2</sub> show a good correlation with the yield curve of &#x201c;Zheng mai 7698&#x201d;, second only to WI <sub>Plus</sub>, while their correlation with the yield curve of &#x201c;Yang mai 11&#x201d; is the lowest. Upon comparison, it is found that the order of correlation between WI <sub>Plus</sub>, WI <sub>Growth</sub> and WI <sub>Destiny</sub> with the yield curve of both varieties is WI <sub>Plus</sub> &gt; WI <sub>Growth</sub> &gt; WI <sub>Destiny</sub>.</p>
<p>In conclusion, the use of low oxygen stress as the waterlogging index (WI) can quantitatively characterize the degree of waterlogging in winter wheat. The improved WI <sub>Plus</sub>, which takes into account soil temperature and growth stage, better reflects the waterlogging severity of wheat.</p>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Validation of the accuracy of the improved index (WI <sub>Plus</sub> and WI <sub>Growth</sub>)</title>
<p>We used the experimental data of soil moisture, temperature, and growth period in 2021 to calculate the values of WI <sub>Destiny</sub>, WI <sub>Growth</sub>, and WI <sub>Plus</sub> for 2021. Then we calculated the predicted yields by bringing the three WI values into the equations in <xref ref-type="table" rid="T1">
<bold>Tables&#xa0;1</bold>
</xref>, <xref ref-type="table" rid="T2">
<bold>2</bold>
</xref>. Compared with the observed yield, it was found that the improved WI <sub>Plus</sub> index provided more accurate predictions for the yield of two wheat varieties (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Comparison of observed (Observed) and simulated (Sim). The left is &#x201c;Yang mai 11&#x201d;, and the right is &#x201c;Zheng mai 7698&#x201d;. R<sup>2</sup> represents the coefficient of determination, and RMSE represents the Root Mean Squared Error.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1262001-g006.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>The middle and lower reaches of the Yangtze River region are the main wheat-growing areas in China. They are greatly influenced by the subtropical monsoon climate. In the months of March and April each year, extensive waterlogging occurs, and water stress can lead to a decrease in wheat yield or even complete crop failure (<xref ref-type="bibr" rid="B9">Ghobadi et&#xa0;al., 2017</xref>). This study has established five WI indices and analyzed them through pot experiments on two wheat varieties, &#x201c;Yang mai 11&#x201d; and &#x201c;Zheng mai 7698,&#x201d; under different waterlogging treatments.</p>
<p>Among the four waterlogging treatments, there is a strong negative correlation between the duration of waterlogging and the wheat yield increment/decrement rate. In other words, the longer the waterlogging period, the lower the wheat yield increment/decrement rate. When comparing the effects of waterlogging during the jointing stage and the booting stage on wheat yield increment/decrement rate, it was found that waterlogging during the booting stage had a greater impact on wheat yield compared to the jointing stage, which is consistent with the findings of Ghobadi et&#xa0;al. Additionally, the study discovered that, under the same duration of waterlogging, the slope of the yield increment/decrement curve for a single waterlogging treatment was greater than that of the curve for two waterlogging treatments. This could be attributed to the adaptive phase that occurs when the waterlogging duration reaches a critical point (5 days) in wheat growth, where the crop gradually adapts and recovers normal growth. This includes processes such as stomatal regulation, physiological metabolism adjustment, root adaptation, and nutrient regulation to cope with the impact of waterlogged conditions. During this phase, wheat relies on self-regulation and adaptive capacity to withstand the stress caused by waterlogging and maintain growth and yield (<xref ref-type="bibr" rid="B33">Shaw and Meyer, 2015</xref>). Therefore, although there is a strong negative correlation between the duration of waterlogging and the wheat yield increment/decrement rate, the linear simulation curve parameters of the waterlogging duration and wheat yield increment/decrement rate are significantly influenced by factors such as the waterlogging period, frequency, and wheat variety. Using only the duration of waterlogging as a universal characteristic to represent the degree of waterlogging in wheat is not highly applicable.</p>
<p>It has been found through research that there is a correlation between the five WI indices and wheat yield. Before reaching the threshold, the fitting curve of the five WI indices and wheat yield shows an upward trend, indicating that a higher WI value corresponds to a higher wheat yield. However, once the threshold is exceeded, the curve starts to decline, and an increase in WI value will lead to a decrease in wheat yield. In the second stage of the fitting curve, the five WI indices show the following characteristics in relation to wheat yield: for the &#x201c;Zheng mai 7698&#x201d; variety, the curve&#x2019;s slope relating to the WI index and yield is greater than that of the &#x201c;Yang mai 11&#x201d; variety. This indicates that the &#x201c;Yang mai 11&#x201d; variety has stronger tolerance to waterlogging, while the &#x201c;Zheng mai 7698&#x201d; variety is less capable of withstanding waterlogged conditions, which is consistent with the findings of Yu J J and others (<xref ref-type="bibr" rid="B40">Yu et&#xa0;al., 2014</xref>). Therefore, it can be seen that the WI indices are influenced by the tolerance of the variety to waterlogging, although the present study did not consider the tolerance of different varieties. However, the WI indices can reflect the level of waterlogging tolerance among different wheat varieties.</p>
<p>Many studies have shown a positive correlation between soil respiration and temperature (<xref ref-type="bibr" rid="B25">Lloyd and Taylor, 1994</xref>; <xref ref-type="bibr" rid="B5">Das and Paul, 2015</xref>). Considering the influence of temperature on soil oxygen content under waterlogged conditions where oxygen diffusion is restricted, the commonly used fixed oxygen threshold may not be suitable. Experimental results on wheat also indicate that the improved WI <sub>Plus</sub> is more suitable for quantitatively analyzing the impact of low oxygen stress on wheat yield. The threshold for WI <sub>Plus</sub> is 14.9, where WI <sub>Plus</sub> &#x2264; 14.9 (&#x201c;Zheng mai 7698&#x201d;) and &#x2264; 10.8 (&#x201c;Yang mai 11&#x201d;) show a positive correlation between WI <sub>Plus</sub> and wheat yield, indicating that the wheat is not waterlogged. When WI <sub>Plus</sub> &gt; 14.9 (&#x201c;Zheng mai 7698&#x201d;) and &gt; 10.8 (&#x201c;Yang mai 11&#x201d;), there is a negative correlation between WI <sub>Plus</sub> and wheat yield, indicating that wheat is affected by waterlogging, and a higher WI value indicates more severe waterlogging.</p>
<p>This study calculates the soil aeration index based on soil volumetric water content and constructs the WI index using characteristic values of wheat root responses to low oxygen (anoxic) stress, making it applicable in general. Additionally, advancements in remote sensing technology, particularly satellite remote sensing, have shown that soil moisture can be measured using various remote sensing techniques (<xref ref-type="bibr" rid="B38">Wang et&#xa0;al., 2022</xref>). Firstly, spectral pictures of wheat planting regions may be collected using high-resolution remote sensing equipment carried by satellites or drones. Secondly, we can preprocess the collected images, including radiation calibration, atmospheric correction, geocoding, and other steps, to eliminate noise and distortion in the images and ensure the accuracy of the data. Thirdly, using image processing technology, extract features related to wheat waterlogging, such as soil temperature and soil moisture content. Finally, the relevant feature quantities were brought into the exponential model of this study, and evaluated and verified through yield data. This can achieve a wider range of wheat stain monitoring. However, its practical application is limited by factors such as insufficient temporal resolution, difficulty in extracting characteristic quantities of waterlogging, and large computational load. Therefore, in the application process, these factors need to be fully considered, and these effects can be reduced through optimization methods.</p>
<p>Currently, there are three remote sensing methods commonly used for soil moisture retrieval: optical remote sensing, thermal infrared remote sensing, and microwave remote sensing. Among them, microwave remote sensing is unaffected by various meteorological conditions and has strong cloud penetration capability. Radar backscatter coefficients are sensitive to changes in soil moisture (<xref ref-type="bibr" rid="B20">Li et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B38">Wang et&#xa0;al., 2022</xref>). For example, Xiong Q X., et&#xa0;al. used Sentinel-1A data to extract soil surface relative moisture content data. Then, combined with the precipitation index in the previous period, they used the Kalman filter interpolation method to obtain daily soil moisture information in the observation area. In addition, remote sensing technology can also obtain soil temperature information (<xref ref-type="bibr" rid="B39">Xiong et al., 2021</xref>). This provides more opportunities for the widespread application of the WI index in monitoring waterlogging. However, since this experiment was conducted in pots, there is minimal vertical variation in soil volumetric water content, which significantly differs from field moisture characteristics. If the results are further applied in field conditions, it would be an area for future improvement.</p>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusions</title>
<p>Using only the duration of waterlogging cannot quantitatively reflect the severity of wheat waterlogging. It is also related to factors such as the tolerance of wheat varieties and the frequency of waterlogging events. Furthermore, the WI index can be applied to assess waterlogging damage under excessive irrigation conditions and can also reflect the impact of multiple waterlogging events on wheat yield. Comparing WI <sub>Destiny</sub>, WI <sub>Growth</sub>, WI <sub>ASD1</sub> and WI <sub>ASD2</sub>, the improved WI <sub>plus</sub> shows better performance. When WI <sub>plus</sub> is &#x2264;14.9 (&#x201c;Zheng mai 7698&#x201d;) or &#x2264;10.8 (&#x201c;Yang mai 11&#x201d;), wheat does not experience waterlogging. When WI <sub>plus</sub> is &gt;14.9 (&#x201c;Zheng mai 7698&#x201d;) or &gt;10.8 (&#x201c;Yang mai 11&#x201d;), wheat experiences waterlogging.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Material</bold>
</xref>. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>F-LX: Writing-original draft, Formal analysis, Methodology. P-MH: Writing &#x2013; review &amp; editing, Funding acquisition. XW: Writing &#x2013; review &amp; editing, Date curation, Investigation. MH: Writing &#x2013; review &amp; editing, Revise the paper. KL: Writing &#x2013; review &amp; editing, Revise the paper. Q-XX: Writing &#x2013; review &amp; editing, Conceptualization, Methodology. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. National Natural Science Foundation of China (No. 31871516) Grain Research Development Corporation (UOT2306-001RTX)</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>The authors would like to thank Xiao-Qin Yang, Qian Xiang, Yu-Jin Huang and other members of the research group for their assistance.</p>
</ack>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2023.1262001/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2023.1262001/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Table_1.xlsx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
<supplementary-material xlink:href="Table_2.xlsx" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
<supplementary-material xlink:href="Table_3.xlsx" id="SM3" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
<supplementary-material xlink:href="Table_4.xlsx" id="SM4" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
<supplementary-material xlink:href="Table_5.xlsx" id="SM5" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amthor</surname> <given-names>J. S.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>The McCree&#x2013;de Wit&#x2013;Penning de Vries&#x2013;Thornley respiration paradigms: 30 years later</article-title>. <source>Ann. Bot.</source> <volume>86</volume> (<issue>1</issue>), <fpage>1</fpage>&#x2013;<lpage>20</lpage>. doi: <pub-id pub-id-type="doi">10.1006/anbo.2000.1175</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Asseng</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Keating</surname> <given-names>B. A.</given-names>
</name>
<name>
<surname>Fillery</surname> <given-names>I. P. R.</given-names>
</name>
<name>
<surname>Gregory</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Bowden</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Turner</surname> <given-names>N. C.</given-names>
</name>
<etal/>
</person-group>. (<year>1998</year>). <article-title>Performance of the APSIM-wheat module in Western Australia</article-title>. <source>Field Crops Res.</source> <volume>57</volume>, <fpage>163</fpage>&#x2013;<lpage>179</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0378-4290(97)00117-2</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barber</surname> <given-names>K. R.</given-names>
</name>
<name>
<surname>Leeds-Harrison</surname> <given-names>P. B.</given-names>
</name>
<name>
<surname>Lawson</surname> <given-names>C. S.</given-names>
</name>
<name>
<surname>Gowing</surname> <given-names>D. J. G.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Soil aeration status in a lowland wet grassland</article-title>. <source>Hydrol. Processes</source> <volume>18</volume> (<issue>2</issue>), <fpage>329</fpage>&#x2013;<lpage>341</lpage>. doi: <pub-id pub-id-type="doi">10.1002/hyp.1378</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bartholomeus</surname> <given-names>R. P.</given-names>
</name>
<name>
<surname>Witte</surname> <given-names>J. P. M.</given-names>
</name>
<name>
<surname>van Bodegom</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>van Dam</surname> <given-names>J. C.</given-names>
</name>
<name>
<surname>Aerts</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Critical soil conditions for oxygen stress to plant roots: Substituting the Feddes-function by a process-based model</article-title>. <source>J. Hydrol.</source> <volume>360</volume> (<issue>1-4</issue>), <fpage>147</fpage>&#x2013;<lpage>165</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jhydrol.2008.07.029</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Das</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Paul</surname> <given-names>P. K.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Present status of soil moisture estimation by microwave remote sensing</article-title>. <source>Cogent Geosci.</source> <volume>1</volume> (<issue>1</issue>), <fpage>1084669</fpage>. doi: <pub-id pub-id-type="doi">10.1080/23312041.2015.1084669</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Doran</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Mielke</surname> <given-names>L. N.</given-names>
</name>
<name>
<surname>Power</surname> <given-names>J. F.</given-names>
</name>
</person-group> (<year>1990</year>). &#x201c;<article-title>Microbial activity as regulated by soil water-filled pore space</article-title>,&#x201d; in <conf-name>Transactions 14th International Congress of Soil Science</conf-name>, <conf-loc>Kyoto, Japan</conf-loc>, Vol. <volume>III</volume>. <fpage>94</fpage>&#x2013;<lpage>99</lpage>.</citation>
</ref>
<ref id="B7">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Evans</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Skaggs</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>1993</year>). &#x201c;<article-title>Stress day index models to predict corn and soybean yield response to water table management</article-title>,&#x201d; in <conf-name>15th congress on irrigation and drainage &#x2018;Water management in the next century &#x2018;Workshop on subsurface drainage simulation</conf-name>. <fpage>15</fpage> p. <publisher-name>IRSTEA</publisher-name>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Su</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>A. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Research progress of the response mechanism of wheat growth to waterlogging stress and the related regulating managements</article-title>. <source>Chin. J. Appl. Ecol.</source> <volume>31</volume> (<issue>12</issue>), <fpage>4321</fpage>&#x2013;<lpage>4330</lpage>. doi: <pub-id pub-id-type="doi">10.13287/J.1001-9332.202012.028</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ghobadi</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Ghobadi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Zebarjadi</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Effect of waterlogging at different growth stages on some morphological traits of wheat varieties</article-title>. <source>Int. J. Biometeorol.</source> <volume>61</volume> (<issue>4</issue>), <fpage>635</fpage>&#x2013;<lpage>645</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00484-016-1240-x</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Githui</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Beverly</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Aiad</surname> <given-names>M.</given-names>
</name>
<name>
<surname>McCaskill</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Modelling waterlogging impacts on crop growth: a review of aeration stress definition in crop models and sensitivity analysis of APSIM</article-title>. <source>Int. J. Plant Biol.</source> <volume>13</volume> (<issue>3</issue>), <fpage>180</fpage>&#x2013;<lpage>200</lpage>. doi: <pub-id pub-id-type="doi">10.3390/ijpb13030017</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Si</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Dynamic risk assessment of waterlogging disaster for maize based on CERES-Maize model in Midwest of Jilin Province, China</article-title>. <source>Natural Hazards</source> <volume>83</volume>, <fpage>1747</fpage>&#x2013;<lpage>1761</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s11069-016-2391-0</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Climate change benefits negated by extreme heat</article-title>. <source>Nat. Food</source> <volume>2</volume>, <fpage>855</fpage>&#x2013;<lpage>856</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s43016-021-00387-6</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Roggero</surname> <given-names>P. P.</given-names>
</name>
<name>
<surname>Zavattaro</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Simple, efficient and robust techniques for automatic multi-objective function parameterization: case studies of local and global optimization using APSIM</article-title>. <source>Environ. Model. Softw.</source> <volume>117</volume>, <fpage>109</fpage>&#x2013;<lpage>133</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.envsoft.2019.03.010</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname> <given-names>P. M.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Constructing the wheat waterlogging damage index for risk assessment in the Middle and Lower Reaches of the Yangtze River</article-title>. <source>Trans. Chin. Soc. Agric. Eng. (Transactions CSAE)</source> <volume>02)</volume>, <fpage>100</fpage>&#x2013;<lpage>106</lpage>. doi: <pub-id pub-id-type="doi">10.11975/j.issn.1002-6819.202207005</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ibrahim</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Meinke</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Fan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A regulator of early flowering in barley (Hordeum vulgare L.)</article-title>. <source>Plops One</source> <volume>13</volume> (<issue>7</issue>), <fpage>1</fpage>&#x2013;<lpage>16</lpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0200722</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jones</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Bland</surname> <given-names>W. L.</given-names>
</name>
<name>
<surname>Ritchie</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Williams</surname> <given-names>J. R.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Simulation of root growth. Modeling plant and soil systems</article-title>. <source>American Society of Agronomy, Inc.</source> <fpage>91</fpage>&#x2013;<lpage>123</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kanwar</surname> <given-names>R. S.</given-names>
</name>
<name>
<surname>Baker</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Mukhtar</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>Excessive soil water effects at various stages of development on the growth and yield of corn</article-title>. <source>Trans. ASAE</source> <volume>31</volume> (<issue>1</issue>), <fpage>133</fpage>&#x2013;<lpage>0141</lpage>. doi: <pub-id pub-id-type="doi">10.13031/2013.30678</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langworthy</surname> <given-names>A. D.</given-names>
</name>
<name>
<surname>Rawnsley</surname> <given-names>R. P.</given-names>
</name>
<name>
<surname>Freeman</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Pembleton</surname> <given-names>K. G.</given-names>
</name>
<name>
<surname>Corkery</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Potential of summer-active temperate (C3) perennial forages to mitigate the detrimental effects of supraspinal temperatures on summer home-grown feed production in south-eastern Australian dairying regions</article-title>. <source>Crop Pasture Sci.</source> <volume>69</volume> (<issue>8</issue>), <fpage>808</fpage>&#x2013;<lpage>820</lpage>. doi: <pub-id pub-id-type="doi">10.1071/CP17291</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Le&#xe3;o</surname> <given-names>T. P.</given-names>
</name>
<name>
<surname>Da Silva</surname> <given-names>A. P.</given-names>
</name>
<name>
<surname>Macedo</surname> <given-names>M. C. M.</given-names>
</name>
<name>
<surname>Imhoff</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Euclides</surname> <given-names>V. P. B.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Least limiting water range: A potential indicator of changes in near-surface soil physical quality after the conversion of Brazilian Savanna into pasture</article-title>. <source>Soil Tillage Res.</source> <volume>88</volume> (<issue>1-2</issue>), <fpage>279</fpage>&#x2013;<lpage>285</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.still.2005.06.014</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>Z. J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y. M.</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>X. L.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>J. S.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Soil moisture retrieval from remote sensing</article-title>. <source>J. Beijing Normal Univ. (Natural Science)</source> <volume>56</volume> (<issue>3</issue>), <fpage>474</fpage>&#x2013;<lpage>481</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.12202/j.0476-0301.2020231</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Archontoulis</surname> <given-names>S. V.</given-names>
</name>
<name>
<surname>Huth</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>D. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Climate change shifts forward flowering and reduces crop waterlogging stress</article-title>. <source>Environ. Res. Lett.</source> <volume>16</volume> (<issue>9</issue>), <fpage>094017</fpage>. doi: <pub-id pub-id-type="doi">10.1088/1748-9326/ac1b5a</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Ibrahim</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Nuruzzaman Manik</surname> <given-names>S. M.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Tian</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Genetic factors increasing barley grain yields under soil waterlogging</article-title>. <source>Food Energy Secur.</source> <volume>9</volume> (<issue>4</issue>), <elocation-id>e238</elocation-id>. doi: <pub-id pub-id-type="doi">10.1002/fes3.238</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>D. L.</given-names>
</name>
<name>
<surname>Meinke</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Hoogenboom</surname> <given-names>G.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Silver lining to a climate crisis in multiple prospects for alleviating crop waterlogging under future climates</article-title>. <source>Nat. Commun.</source> <volume>14</volume> (<issue>1</issue>), <fpage>765</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41467-023-36129-4</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Lizaso</surname> <given-names>J. I.</given-names>
</name>
<name>
<surname>Ritchie</surname> <given-names>J. T.</given-names>
</name>
</person-group> (<year>1997</year>). &#x201c;<article-title>A modified version of CERES to predict the impact of soil water excess on maize crop growth and development</article-title>,&#x201d; in <conf-name>Applications of Systems Approaches at the Field Level: Volume 2 Proceedings of the Second International Symposium on Systems Approaches for Agricultural Development</conf-name>, <conf-loc>IRRI, Los Ba&#xf1;os, Philippines</conf-loc>, <conf-date>6&#x2013;8 December 1995</conf-date>. <fpage>153</fpage>&#x2013;<lpage>167</lpage>. (<publisher-name>Springer Netherlands</publisher-name>).</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lloyd</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Taylor</surname> <given-names>J. A.</given-names>
</name>
</person-group> (<year>1994</year>). <article-title>On the temperature dependence of soil respiration</article-title>. <source>Funct. Ecol.</source> <fpage>315</fpage>&#x2013;<lpage>323</lpage>. doi: <pub-id pub-id-type="doi">10.2307/2389824</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Meyer</surname> <given-names>W. S.</given-names>
</name>
<name>
<surname>Godwin</surname> <given-names>D. C.</given-names>
</name>
<name>
<surname>White</surname> <given-names>R. J. G.</given-names>
</name>
</person-group> (<year>1996</year>). &#x201c;<article-title>SWAGMAN&#xae; Destiny:A tool to predict productivity change due to salinity, waterlogging and irrigation management</article-title>,&#x201d; in <conf-name>&#x2018;Proceedings of the 8th Australian Agronomy Conference&#x2019;</conf-name>, <conf-loc>Toowomba, Qld</conf-loc>. <fpage>425</fpage>&#x2013;<lpage>428</lpage>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Phelan</surname> <given-names>D. C.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Kemmerer</surname> <given-names>E. P.</given-names>
</name>
<name>
<surname>Parsons</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Management opportunities for boosting productivity of cool-temperate grazed dairy farms under climate change</article-title>. <source>Agric. Syst.</source> <volume>138</volume>, <fpage>46</fpage>&#x2013;<lpage>54</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.agsy.2015.05.005</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qian</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. G.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>W. B.</given-names>
</name>
<name>
<surname>Jia</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Yield reduction analysis and determination of drainage index in cotton under waterlogging followed by submergence</article-title>. <source>Trans. Chin. Soc. Agric. Eng.</source> <volume>13)</volume>, <fpage>89</fpage>&#x2013;<lpage>97</lpage>. doi: <pub-id pub-id-type="doi">10.11975/j.issn.1002-6819.2015.13.013</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qian</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. G.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>W. B.</given-names>
</name>
<name>
<surname>Qi</surname> <given-names>Z. M.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>H. W.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>Y. Y.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>An improved CROPR model for estimating cotton yield under soil aeration stress</article-title>. <source>Crop Pasture Sci.</source> <volume>68</volume> (<issue>4</issue>), <fpage>366</fpage>&#x2013;<lpage>377</lpage>. doi: <pub-id pub-id-type="doi">10.1071/CP16426</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rawnsley</surname> <given-names>R. P.</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>A. P.</given-names>
</name>
<name>
<surname>Christie</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Eckard</surname> <given-names>R. J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Current and future direction of nitrogen fertiliser use in Australian grazing systems</article-title>. <source>Crop Pasture Sci.</source> <volume>70</volume> (<issue>12</issue>), <fpage>1034</fpage>&#x2013;<lpage>1043</lpage>. doi: <pub-id pub-id-type="doi">10.1071/CP18566</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>S&#xe1;ndor</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Ehrhardt</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Grace</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Recous</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Snow</surname> <given-names>V.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Ensemble modelling of carbon fluxes in grasslands and croplands</article-title>. <source>Field Crops Res.</source> <volume>252</volume>, <fpage>107791</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fcr.2020.107791</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shahpari</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Allison</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Stanley</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>An integrated economic, environmental and social approach to agricultural land-use planning</article-title>. <source>Land</source> <volume>10</volume> (<issue>4</issue>), <fpage>1</fpage>&#x2013;<lpage>18</lpage>. doi: <pub-id pub-id-type="doi">10.3390/land10040364</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shaw</surname> <given-names>R. E.</given-names>
</name>
<name>
<surname>Meyer</surname> <given-names>W. S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Improved empirical representation of plant responses to waterlogging for simulating crop yield</article-title>. <source>Agron. J.</source> <volume>107</volume> (<issue>5</issue>), <fpage>1711</fpage>&#x2013;<lpage>1723</lpage>. doi: <pub-id pub-id-type="doi">10.2134/agronj14.0625</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shaw</surname> <given-names>R. E.</given-names>
</name>
<name>
<surname>Meyer</surname> <given-names>W. S.</given-names>
</name>
<name>
<surname>McNeill</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Tyerman</surname> <given-names>S. D.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Waterlogging in Australian agricultural landscapes: a review of plant responses and crop models</article-title>. <source>Crop Pasture Sci.</source> <volume>64</volume> (<issue>6</issue>), <fpage>549</fpage>&#x2013;<lpage>562</lpage>. doi: <pub-id pub-id-type="doi">10.1071/CP13080</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Singhal</surname> <given-names>R. K.</given-names>
</name>
<name>
<surname>Fahad</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Kumar</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Choyal</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Javed</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Jinger</surname> <given-names>D.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Beneficial elements: New Players in improving nutrient use efficiency and abiotic stress tolerance</article-title>. <source>Plant Growth Regul.</source>, (<issue>2</issue>) <fpage>237</fpage>&#x2013;<lpage>265</lpage>. doi: <pub-id pub-id-type="doi">10.1007/S10725-022-00843-8</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="other">
<person-group person-group-type="author">
<name>
<surname>Skaggs</surname> <given-names>R. W.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>DRAINMOD: A simulation model for shallow water table soils</article-title>. Available at: <ext-link ext-link-type="uri" xlink:href="https://tigerprints.clemson.edu/cgi/viewcontent.cgi?article=1210&amp;context=scwrc">https://tigerprints.clemson.edu/cgi/viewcontent.cgi?article=1210&amp;context=scwrc</ext-link>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wan</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>X. Q.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y. T.</given-names>
</name>
<name>
<surname>Xiang</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Yuan</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>S. Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Change characteristics of soil oxygen content of waterlogged wheat and waterlog damage index determination</article-title>. <source>J. South. Agric.</source> <volume>(11)</volume>, <fpage>3058</fpage>&#x2013;<lpage>3067</lpage>. doi: <pub-id pub-id-type="doi">10.3969/j.issn.2095-1191.2022.11.006</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Soil moisture retrieval by integrating SAR and optical data over winter wheat fields</article-title>. <source>Appl. Sci.</source> <volume>12</volume> (<issue>23</issue>), <fpage>12057</fpage>. doi: <pub-id pub-id-type="doi">10.3390/app122312057</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiong</surname> <given-names>Q. X.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Extracting the spatial distribution information of crop sub-surface waterlogging using antecedent precipitation index and sentinel-1A SAR data</article-title>. <source>Trans. Chin. Soc Agric. Eng.</source> <volume>37</volume>, <fpage>175</fpage>&#x2013;<lpage>183</lpage>. doi: <pub-id pub-id-type="doi">10.11975/j.issn.1002-6819.2021.05.020</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. Y.</given-names>
</name>
<name>
<surname>Duan</surname> <given-names>Y. Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. L.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Studies on the waterlogging resistance of main wheat varieties in Jianghan Plain</article-title>. <source>Hubei Agric. Sci.</source> <volume>53</volume> (<issue>4</issue>), <fpage>760</fpage>&#x2013;<lpage>764</lpage>. doi: <pub-id pub-id-type="doi">10.14088/j.cnki.issn0439-8114.2014.04.025</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>H. Y.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>H. D.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Risk assessment of rapeseed spring waterlogging in the middle and lower reaches of the Yangtze River</article-title>. <source>Meteorol. Environ. Sci.</source> <volume>42</volume> (<issue>1</issue>), <fpage>11</fpage>&#x2013;<lpage>17</lpage>. doi: <pub-id pub-id-type="doi">10.16765/j.cnki.1673-7148.2019.01.002</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>