Calculation method, system and device for simulating crop yield loss risk based on stomatal ozone absorption flux

Abstract

The invention provides a calculation method, system and device for simulating crop yield loss risk based on stomatal ozone absorption flux. The calculation method for simulating the crop yield loss risk based on the stomatal ozone absorption flux comprises the following steps: S1, obtaining test data related to stomatal conductance and yield in a rice growth season; s2, parameterizing a rice stomatal conductance model, and calculating rice stomatal conductance per hour by utilizing hour meteorological data; and S3, calculating the absorption flux of the crop pores to O3 under different threshold values. According to the calculation method, system and device for simulating the crop yield loss risk based on the stomatal ozone absorption flux provided by the invention, the risk that the influence of O3 on the crop yield cannot be accurately quantified in the prior art can be solved, and the influence of environmental factors on stomatal conductance and the difference of O3 damage threshold values of different crops are brought into crop yield risk quantification; quantitative accuracy can be improved, and method support is provided for accurately predicting crop yield loss caused by O3 and guaranteeing grain safety.

Description

Technical Field

[0001] The present invention belongs to the technical fields of agricultural meteorology and agricultural informatization, and particularly relates to a calculation method, system and device for quantitatively quantifying the risk of crop yield loss based on the ozone absorption flux of stomata. Background Art

[0002] With the rapid progress of industry and urbanization, the emissions of ozone (O 3 ) precursors are continuously increasing, resulting in varying degrees of increase in O 3 concentration globally. O 3 is a strongly oxidizing gas that is toxic to crop growth. It enters the crop through the leaf stomata, thereby damaging the photosynthetic rate of the leaves and reducing crop yields. Currently, in regions such as Asia and North America, the environmental O 3 concentration has exceeded the damage threshold of plants. O 3 pollution has become an important factor threatening China's food security. Reasonably quantifying the impact of O 3 pollution on crop yields is the key to ensuring food security;

[0003] Currently, the methods for quantifying the impact of O 3 pollution on crop yields are usually based on O 3 concentration and exposure time, ignoring the differences in the O 3 absorption amount due to the differences in crop leaf stomatal conductance and the differences in the sensitivity of different crops to O 3 . For example, the invention patent with the application number 202310246982.1 discloses a method and device for evaluating the risk of crop yield reduction caused by ozone pollution, which uses the cumulative index (AOT40) of the hourly O 3 concentration exceeding the threshold of 40 ppb to evaluate the risk of yield reduction caused by O 3 pollution. The existing technologies have the following disadvantages: (1) Under different climate conditions, the stomatal conductance of crop leaves will change. Only using the crop O 3 exposure dose AOT40 to evaluate the risk of yield reduction caused by O 3 pollution ignores the differences in the O 3 absorption flux of crop stomata under different climates and cannot accurately reflect the O 3 flux entering the crop through the stomata; (2) Different crops have different sensitivities to O 3 . The O 3 exposure dose AOT40 uses 40 ppb as the damage threshold for different crops, which is difficult to represent the differences in the sensitivity of different crops to O 3 and is difficult to be used to quantify the impact of O 3 on the yields of different crops; (3) Currently, constructing O 3The method of exposure dose and crop yield response requires prior normalization of the yield before integrating data from different experiments to construct the O 3 Exposure dose-yield response equation. In the normalization process, a linear equation fitting is mostly used, which ignores the O 3 Possible non-linear response relationship between the exposure dose and crop yield, resulting in an inaccurate dose response equation.

[0004] Therefore, the existing technology ignores the response changes of crop stomatal conductance to the environment in quantitatively O 3 The risk of crop yield impact, and also does not consider the differences in O 3 Sensitivity of different crops, making it difficult to provide effective guidance for mitigating O 3 The risk of crop yield loss. Therefore, it is necessary to provide a new method for quantitatively assessing the risk of crop yield loss based on stomatal ozone uptake flux to solve the above technical problems. Summary of the Invention

[0005] The technical problem solved by the present invention is to provide a method that can solve the problem that the existing technology cannot accurately quantify O 3 The risk of crop yield impact. By incorporating the influence of environmental factors on stomatal conductance and the differences in O 3 Injury thresholds of different crops into the quantification of crop yield risk, the accuracy of quantification can be improved, providing a method support for accurately predicting O 3 The risk of crop yield loss and ensuring food security, a method for quantitatively assessing the risk of crop yield loss based on stomatal ozone uptake flux.

[0006] To solve the above technical problems, the calculation method for simulating the risk of crop yield loss based on stomatal ozone uptake flux provided by the present invention includes:

[0007] S1. Obtain the experimental data related to stomatal conductance and yield during the rice growing season;

[0008] S2. Parameterize the rice stomatal conductance model and calculate the hourly rice stomatal conductance using hourly meteorological data;

[0009] S3. Calculate the absorption flux of crop stomata to O 3 at different thresholds;

[0010] S4. Construct the response relationship between rice yield and stomatal O 3 Absorption flux;

[0011] S5. Obtain the hourly meteorological data and O 3 Concentration data during the rice growth cycle in the region, calculate the absorption flux of rice stomata to O 3 and calculate the yield loss caused by O 3 ​

[0012] As a further solution of the present invention, the step S1 obtains experimental data related to stomatal conductance and yield during the rice growth season, and the step S2 parameterizes the rice stomatal conductance model, and calculates the hourly stomatal conductance of rice during the growth period according to the maximum stomatal conductance of the crop, photosynthetically active radiation, growing degree days, air temperature, vapor pressure deficit, and soil water content.

[0013] As a further solution of the present invention, the model calculation of the stomatal conductance adopts the following formula:

[0014] g s = g max ·f(GDD)·f(PAR)·f(T a )·f(VPD)·f(PAW)

[0015] where g max is the maximum stomatal conductance of the crop during the growth period, GDD represents the accumulated temperature above 0 degrees, with the flowering period of the crop as 0-degree accumulated temperature, the accumulated temperature before flowering is negative, and the accumulated temperature after flowering is positive, PAR represents photosynthetically active radiation, Ta represents air temperature, VPD represents vapor pressure deficit, and PAW represents soil water content.

[0016] As a further solution of the present invention, the phenological function f(GDD) is expressed as follows:

[0017] When A start ≤GDD < A start + phen 1

[0018]

[0019] When A start + phen 1 ≤GDD ≤ A end - phen 2

[0020] f(GDD)= 1

[0021] When A end - phen 2 <GDD ≤ A end

[0022]

[0023] where A start , A end respectively represent the accumulated temperature at the beginning and end of the crop growth period, phen 1 , phen 2 , phen a and phen eare model parameters;

[0024] The light intensity function f(PAR) is expressed as follows:

[0025] f(PAR) = 1 - exp(-light a ·PAR)

[0026] where light a are model parameters;

[0027] The temperature function f(T a ) is expressed as follows:

[0028]

[0029] where T min and T max are model parameters;

[0030] The vapor pressure deficit function f(VPD) is expressed as follows:

[0031]

[0032] where VPD min and VPD max are model parameters.

[0033] As a further aspect of the present invention, the soil water function f(PAW) is expressed as follows:

[0034]

[0035] where PAW t are model parameters;

[0036] By actually measuring the leaf stomatal conductance at different time points during the rice growth season, the stomatal conductance under different accumulated temperatures, light intensities, temperatures, and vapor pressure deficits was obtained. A total of 1625 pairs of data for four main cultivated conventional japonica rice varieties (Nanjing 5055, Wuyunjing 3, Wuyujing 27, and Huaidao 5) were obtained for the estimation of stomatal conductance model parameters;

[0037] Using the above-obtained stomatal conductance model parameters, and then inputting the hourly meteorological data during the rice growth period, the hourly actual stomatal conductance value of the rice can be calculated. By comparing the simulated value and the measured value, the model of this study has good accuracy. This step effectively simulates the stomatal conductance of conventional japonica rice during the growth cycle.

[0038] As a further aspect of the present invention, in step S3, the absorption flux of stomata to O 3 under different thresholds during the rice growth cycle is calculated, and according to the monitored hourly O 3Calculate the absorption flux of rice stomata to O based on the concentration and stomatal conductance calculated in step S2. The hourly absorption amount of O by crop stomata is calculated using the following formula: 3 The hourly absorption amount of O by crop stomata is calculated using the following formula: 3 F

[0039] F st =[O 3 ·g s ·0.663

[0040] Where g s is the hourly stomatal conductance calculated in step S2, and [O 3 is the corresponding O 3 concentration per hour. Then, the cumulative value of the absorption flux of stomata to O 3 exceeding the threshold y during the crop growth period is calculated as follows:

[0041] POD y =∫(F st -y)dt

[0042] Based on the above formula, calculate the stomatal O 3 absorption flux POD y at different thresholds y.

[0043] As a further solution of the present invention, step S4 is to construct the response relationship between rice yield and the absorption flux of stomata to O 3 Since different crops have different sensitivities to O 3 , the yield will be affected only when the absorption flux of stomata to O 3 exceeds a certain threshold y. In order to estimate the threshold y using the observed data and the absorption flux of stomata to O 3 , first, a non-linear response relationship between POD y and POD 0 is constructed. The function formula is as follows:

[0044] POD y =POD 0 ·exp(-a1·y)·(1-exp(-((POD 0 / (b1(exp(-b2·y)

[0045] -exp(-b3·y))))^(c1+c2·y))))

[0046] Where y, a1, b1, b2, b3, c1, and c2 are equation parameters. The relevant parameters of rice are obtained by fitting the function. In this study, the response relationship between POD 0 and POD y ;

[0047] Rice yield loss and stomatal O3 The cumulative absorption flux can be calculated using the following formula:

[0048] RY = 1 - S * POD y

[0049] where RY represents the relative yield loss of rice, S represents the sensitivity coefficient of rice, and POD y represents the stomatal O with a threshold of y 3 cumulative absorption flux. By combining the PODy conversion equation and the yield loss response to fit the observed data, the response relationship between rice yield loss and POD0 is obtained. At the same time, parameters such as y, a1, b1, b2, b3, c1, c2, and S are fitted, and the non-linear response equation of relative yield and POD0 is obtained.

[0050] As a further aspect of the present invention, the step S5 obtains the hourly meteorological data and O 3 concentration data during the rice growth period in the region, calculates the stomatal O 3 absorption flux of rice, calculates the yield loss caused by O 3 By obtaining the data of 217 stations in the middle and lower reaches of the Yangtze River and inputting them into the stomatal conductance model, the stomatal absorption flux POD0 is calculated. Using the stomatal O 3 absorption flux with a threshold of 14.40, the risk of O 3 causing rice yield loss in the Yangtze River Delta is calculated.

[0051] The present invention also provides a system for simulating the risk of crop yield loss based on stomatal ozone absorption flux, including:

[0052] Data acquisition module: used to acquire experimental data related to stomatal conductance and yield during the rice growth season;

[0053] Stomatal conductance calculation module: used to parameterize the rice stomatal conductance model and calculate the hourly stomatal conductance of rice using hourly meteorological data;

[0054] Stomatal O 3 Absorption flux calculation module: used to calculate the absorption flux of rice stomata to O 3 at different thresholds;

[0055] Dose-response relationship fitting module: used to calculate the dose-response relationship between rice yield and stomatal O 3 absorption flux at the optimal threshold;

[0056] Crop yield loss risk calculation module: used to obtain the hourly meteorological data and O 3 concentration data during the rice growth period in the region, calculate the stomatal O 3 absorption flux of rice, and quantify the risk of rice yield loss caused by O 3 ​

[0057] The present invention also provides a device for simulating the risk of crop yield loss based on the stomatal ozone absorption flux, including a processor and a memory. The memory stores a computer program, and the processor executes the computer program to implement the steps of the above method.

[0058] Compared with the related art, the calculation method, system and device for simulating the risk of crop yield loss based on the stomatal ozone absorption flux provided by the present invention have the following beneficial effects:

[0059] 1. The present invention integrates meteorological conditions and the physiological characteristics of crops into the calculation of stomatal conductance, and can accurately calculate the O 3 flux entering the plant body through the stomata, and more accurately calculates the impact of O 3 on the crop yields in different regions, providing a feasible method for regional assessment;

[0060] 2. The method of the present invention constructs a conversion equation between the stomatal O 3 absorption fluxes under different thresholds, providing favorable conditions for estimating the stomatal O 3 absorption flux thresholds of different crops;

[0061] 3. The present invention constructs a method of integrating the yields of different varieties of rice, and constructs a non-linear response relationship between the stomatal O 3 absorption flux and the rice yield through a one-step fitting method. The fitting accuracy of the dose-response equation is improved, which is beneficial to improving the accuracy of quantitatively predicting the risk of rice yield loss. It provides a method support for accurately predicting the crop yield loss caused by O 3 and ensuring food security. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] For the convenience of those skilled in the art to understand, the present invention will be further described below with reference to the accompanying drawings.

[0063] Figure 1 is a schematic flow chart of the present invention;

[0064] Figure 2 is a schematic diagram comparing the stomatal conductance simulated by the stomatal conductance model of the present invention with the observed stomatal conductance;

[0065] Figure 3 is a schematic diagram of the non-linear response relationship of POD y and POD 0 of the present invention;

[0066] Figure 4 is a schematic diagram of the non-linear response relationship between the relative rice yield loss and the stomatal O 3 absorption flux of the present invention;

[0067] Figure 5 Stomatal O in the Yangtze River Delta region of the present invention 3 Absorption flux distribution diagram. Specific implementation mode

[0068] Next, in combination with the accompanying drawings and specific implementation modes, the present invention will be further described. It should be noted that, on the premise of no conflict, any combination of the following-described embodiments or technical features can form a new embodiment.

[0069] A calculation method for simulating the risk of crop yield loss based on stomatal ozone absorption flux takes Chinese rice as the research object, and specifically elaborates the method of the present invention, including the following steps:

[0070] S1. Obtain experimental data related to stomatal conductance and yield during the rice growing season;

[0071] S2. Parameterize the rice stomatal conductance model and calculate the hourly stomatal conductance of rice using hourly meteorological data;

[0072] S3. Calculate the absorption flux of crop stomata to O 3 at different thresholds;

[0073] S4. Construct the response relationship between rice yield and stomatal O 3 absorption flux;

[0074] S5. Obtain hourly meteorological data and O 3 concentration data during the rice growth cycle in the region, calculate the absorption flux of rice stomata to O 3 , and calculate the yield loss caused by O 3 .

[0075] The step S1 obtains experimental data related to stomatal conductance and yield during the rice growing season.

[0076] The step S2 parameterizes the rice stomatal conductance model and calculates the hourly stomatal conductance of rice during the growth cycle according to the maximum stomatal conductance of the crop, photosynthetically active radiation, growing degree days, air temperature, vapor pressure deficit, and soil water content.

[0077] The model calculation of the stomatal conductance adopts the following formula:

[0078] g s = g max ·f(GDD)·f(PAR)·f(T a )·f(VPD)·f(PAW)

[0079] Among them, g maxis the maximum stomatal conductance of the crop during the growth period. GDD represents the accumulated temperature above 0 °C, with the crop flowering stage as the 0 °C accumulated temperature. The accumulated temperature before flowering is negative, and the accumulated temperature after flowering is positive. PAR represents photosynthetically active radiation, Ta represents air temperature, VPD represents the saturation vapor pressure deficit, and PAW represents soil water content.

[0080] The phenological function f(GDD) is expressed as follows:

[0081] When A start ≤ GDD < A start + phen 1

[0082]

[0083] When A start + phen 1 ≤ GDD ≤ A end - phen 2

[0084] f(GDD) = 1

[0085] When A end - phen 2 < GDD ≤ A end

[0086]

[0087] Wherein, A start , A end respectively represent the accumulated temperatures at the start and end of the crop growth period, phen 1 , phen 2 , phen a and phen e are model parameters.

[0088] The light intensity function f(PAR) is expressed as follows:

[0089] f(PAR) = 1 - exp(-light a · PAR)

[0090] Wherein, light a is a model parameter;

[0091] The temperature function f(T a ) is expressed as follows:

[0092]

[0093] Wherein, T min and T max are model parameters.

[0094] The saturated water vapor pressure difference function f(VPD) is expressed as follows:

[0095]

[0096] where VPD min and VPD max are model parameters.

[0097] The soil moisture function f(PAW) is expressed as follows:

[0098]

[0099] where PAW t is a model parameter;

[0100] By actually measuring the leaf stomatal conductance at different time points during the rice growth season, the stomatal conductance under different accumulated temperatures, light intensities, temperatures, and saturated water vapor pressure differences was obtained. In this case, a total of 1625 pairs of data for four major conventionally cultivated japonica rice varieties in China (Nanjing 5055, Wuyunjing 3, Wuyujing 27, and Huaidao 5) were obtained for the estimation of stomatal conductance model parameters, as Figure 2 shown.

[0101] Using the above-obtained stomatal conductance model parameters, and then inputting the hourly meteorological data during the rice growth period, the hourly actual stomatal conductance value of the rice can be calculated. By comparing the simulated value and the measured value, the model of this study has good accuracy (as Figure 3 shown), and this step effectively simulates the stomatal conductance of conventionally cultivated japonica rice during the growth cycle.

[0102] In step S3, the absorption flux of stomata to O 3 under different thresholds during the rice growth cycle is calculated. According to the monitored hourly O 3 concentration and the stomatal conductance calculated in step S2, the absorption flux of rice stomata to O 3 is calculated. The calculation formula for the hourly O 3 absorption amount by crop stomata is as follows:

[0103] F st =[O 3 ·g s ·0.663

[0104] where g s is the hourly stomatal conductance calculated in step S2, and [O 3 is the corresponding hourly O 3 concentration. Then, the cumulative value of the absorption flux of stomata to O 3 exceeding the threshold y during the crop growth cycle is calculated as follows:

[0105] POD y = ∫(F st - y)dt.

[0106] Based on the above formula, calculate the stomatal O 3 absorption flux POD y .

[0107] Step S4 is to construct the response relationship between rice yield and stomatal O 3 absorption flux. Since different crops have different sensitivities to O 3 , when the stomatal O 3 absorption flux exceeds a certain threshold y, it will affect the yield. In order to estimate the threshold y using the observed data and stomatal O 3 absorption flux, first, a non - linear response relationship between POD y and POD 0 was constructed (see Figure 4 ), and the function formula is as follows:

[0108] POD y = POD 0 · exp(-a1·y)·(1 - exp(-((POD 0 / (b1(exp(-b2·y)-exp(-b3·y))))^(c1 + c2·y))))

[0109] where y, a1, b1, b2, b3, c1, and c2 are equation parameters. The relevant parameters of rice were obtained by fitting the function. In this study, the response relationship between POD 0 and POD y can be seen in Figure 4 .

[0110] The rice yield loss and the cumulative absorption flux of stomatal O 3 can be calculated using the following formula:

[0111] RY = 1 - S * POD y

[0112] where RY represents the relative yield loss of rice, S represents the sensitivity coefficient of rice, and POD y represents the cumulative absorption flux of stomatal O 3 with a threshold of y. By combining the PODy conversion equation and the yield loss response to fit the observed data, the response relationship between the rice yield loss and POD0 was obtained, and at the same time, parameters such as y, a1, b1, b2, b3, c1, c2, and S were fitted. The non - linear response equation between the relative yield and POD0 can be seen in Figure 5 .

[0113] S5. Obtain hourly meteorological data and O 3 concentration data during the rice growth period in the area, calculate the O 3 absorption flux of rice stomata, and calculate the O 3 -induced yield loss. In this case, by obtaining data from 217 stations in the middle and lower reaches of the Yangtze River and inputting them into the stomatal conductance model, the stomatal absorption flux POD0 was calculated, and the risk of O 3 -induced yield loss of rice in the Yangtze River Delta was calculated using the stomatal O 3 absorption flux with a threshold of 14.40.

[0114] The present invention also discloses a system for simulating the risk of crop yield loss based on the stomatal ozone absorption flux, including:

[0115] Data acquisition module: used to obtain experimental data related to stomatal conductance and yield during the rice growth season;

[0116] Stomatal conductance calculation module: used to parameterize the rice stomatal conductance model and calculate the hourly stomatal conductance of rice using hourly meteorological data;

[0117] Stomatal O 3 absorption flux calculation module: used to calculate the absorption flux of rice stomata to O 3 under different thresholds;

[0118] Dose-response relationship fitting module: used to calculate the dose-response relationship between rice yield and stomatal O 3 absorption flux under the optimal threshold;

[0119] Crop yield loss risk calculation module: used to obtain hourly meteorological data and O 3 concentration data during the rice growth period in the area, calculate the O 3 absorption flux of rice stomata, and quantify the risk of O 3 -induced rice yield loss.

[0120] Based on the same inventive concept, the present invention also provides a device for simulating the risk of crop yield loss based on the stomatal ozone absorption flux, including a processor and a memory. The memory stores a computer program, and the processor executes the computer program to implement the steps of the above method.

[0121] Based on the above description, it can be known that the beneficial effects of the present invention compared with the prior art are:

[0122] 1. The present invention integrates meteorological conditions and the physiological characteristics of crops into the calculation of stomatal conductance, can accurately calculate the O 3 flux entering the plant through the stomata, and more accurately calculates the impact of O 3 on crop yields in different regions, providing a feasible method for regional assessment.

[0123] 2. The method of the present invention constructs a conversion equation between stomatal O 3 absorption fluxes at different thresholds, providing favorable conditions for estimating the stomatal O 3 absorption flux thresholds of different crops.

[0124] 3. The present invention constructs a method of integrating the yields of different rice varieties and constructing a non-linear response relationship between stomatal O 3 absorption flux and rice yield by a one-step fitting method. It improves the fitting accuracy of the dose-response equation, which is beneficial to improving the accuracy of quantitatively predicting the risk of rice yield loss. It provides a method support for accurately predicting O 3 causing crop yield loss and ensuring food security.

[0125] The above is only a preferred specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any person skilled in the art within the technical scope disclosed by the present invention, according to the technical solution and inventive concept of the present invention, makes equivalent substitutions or changes, and should be covered by the protection scope of the present invention.

Claims

1. A calculation method for simulating crop yield loss risk based on stomatal ozone absorption flux, characterized in that: The following steps are involved: S1. Obtain experimental data related to stomatal conductance and yield during the rice growing season; S2. Parameterized rice stomatal conductance model, using hourly meteorological data to calculate rice hourly stomatal conductance; S3. Calculate the absorption flux of O3 by crop stomata at different thresholds; S4. Construct the response relationship between rice yield and stomatal O3 absorption flux; S5. Obtain hourly meteorological data and O3 concentration data during the rice growth period in the region, calculate the O3 absorption flux through rice pores, and calculate the yield loss caused by O3.

2. The method for calculating the risk of crop yield loss based on stomatal ozone absorption flux simulation according to claim 1, characterized in that: The step S1 obtains experimental data related to stomatal conductance and yield during the rice growing season, and the step S2 parameterizes the rice stomatal conductance model to calculate the hourly stomatal conductance of rice during the growth cycle according to the maximum stomatal conductance of the crop, photosynthetically active radiation, growth accumulated temperature, air temperature, saturated water vapor pressure deficit and soil moisture content.

3. The method for calculating the risk of crop yield loss based on stomatal ozone absorption flux simulation according to claim 2, characterized in that: The model calculation of stomatal conductance adopts the following formula: g s =g max ·f(GDD)·f(PAR)·f(T a )·f(VPD)·f(PAW) Among them, g max It is the maximum stomatal conductance of the crop during its growth period, GDD indicates the accumulated temperature above 0 degrees, with the flowering period of the crop as 0 degrees accumulated temperature, the accumulated temperature before flowering is negative, and the accumulated temperature after flowering is positive, PAR indicates photosynthetically active radiation, Ta is air temperature, VPD indicates saturated water vapor pressure difference, and PAW indicates soil moisture content.

4. The method for calculating the risk of crop yield loss based on stomatal ozone absorption flux simulation according to claim 3, characterized in that: The phenological function f(GDD) is expressed as follows: When A start ≤ GDD < A start + phen1 When A start +phen1≤GDD≤A end -phen2 f(GDD)=1 When A end -phen2<GDD≤A end Among them, A start , A end Respectively represent the cumulative temperature at the beginning and end of the crop growth period, phen1, phen2, phen a and phen e are model parameters; The light intensity function f(PAR) is expressed as follows: f(PAR)=1-exp(-light a ·BY) Among them, light a are model parameters; The temperature function f(T a ) is expressed as follows: Among them, T min and T max are model parameters; The saturated vapor pressure difference function f(VPD) is expressed as follows: Among them, VPD min and VPD max is the model parameter.

5. The calculation method, system and device for simulating crop yield loss risk based on stomatal ozone absorption flux according to claim 3, characterized in that: The soil moisture function f(PAW) is expressed as follows: Among them, PAW t are model parameters; By measuring the stomatal conductance of leaves at different time points during the rice growing season, we obtained the stomatal conductance under different accumulated temperature, light intensity, temperature and saturated water vapor pressure difference. A total of 1625 pairs of data of four main cultivated conventional japonica rice varieties (Nanjing 5055, Wuyunjing 3, Wuyujing 27 and Huaidao 5) were obtained for the estimation of stomatal conductance model parameters. By using the stomatal conductance model parameters obtained above and then inputting the hourly meteorological data during the rice growth period, the actual hourly stomatal conductance value of rice can be calculated. By comparing the simulated value and the measured value, the model of this study has good accuracy. This step effectively simulates the stomatal conductance of conventional japonica rice during the growth cycle.

6. The method for calculating the risk of crop yield loss based on stomatal ozone absorption flux simulation according to claim 1, characterized in that: The step S3 calculates the O3 absorption flux of the pores at different thresholds during the rice growth cycle, and calculates the O3 absorption flux of the rice pores according to the monitored hourly O3 concentration and the stomatal conductance calculated in step S2. The O3 absorption of the crop pores per hour is calculated using the following formula: F st =[O3]·g s ·0.663 Among them, g s is the stomatal conductance per hour calculated in step S2, [O3] is the corresponding O3 concentration per hour, then the cumulative value of the stomatal O3 absorption flux exceeding the threshold y during the crop growth cycle is calculated as follows: Based on the above formula, the stomatal O3 absorption flux POD under different thresholds y is calculated y .

7. The method for calculating the risk of crop yield loss based on stomatal ozone absorption flux simulation according to claim 1, characterized in that: Step S4 is to construct the response relationship between rice yield and stomatal O3 absorption flux. Due to the different sensitivities of different crops to O3, the yield will be affected only when the stomatal O3 absorption flux exceeds a certain threshold y. In order to estimate the threshold y using the observed data and stomatal O3 absorption flux, a POD y The nonlinear response relationship with POD0 is as follows: POD y =POD0·exp(-a1·y)·(1-exp(-((POD0 / (b1(exp(-b2·y) -exp(-b3·y))))^(c1+c2·y)))) Among them, y, a1, b1, b2, b3, c1 and c2 are equation parameters. The relevant parameters of rice were obtained by fitting the function. In this study, POD0 and POD y Response relationship; Rice yield loss and stomatal O3 cumulative absorption flux can be calculated using the following formula: RY=1-S*POD y Among them, RY represents the relative yield loss of rice, S represents the sensitivity coefficient of rice, and POD y It represents the stomatal O3 cumulative absorption flux with a threshold value of y. By fitting the observed data with the combined PODy conversion equation and the yield loss response, the response relationship between rice yield loss and POD0 was obtained. At the same time, parameters such as y, a1, b1, b2, b3, c1, c2 and S, and the nonlinear response equation of relative yield and POD0 were obtained.

8. The method for calculating the risk of crop yield loss based on stomatal ozone absorption flux simulation according to claim 1, characterized in that: The step S5 obtains hourly meteorological data and O3 concentration data during the rice growth cycle in the region, calculates the rice stomatal O3 absorption flux, and calculates the yield loss caused by O3. The data of 217 stations in the middle and lower reaches of the Yangtze River are obtained and input into the stomatal conductance model to calculate the stomatal absorption flux POD0. The risk of rice yield loss caused by O3 in the Yangtze River Delta is calculated using the stomatal O3 absorption flux with a threshold of 14.

40.

9. A system for simulating crop yield loss risk based on stomatal ozone absorption flux, characterized in that: include: Data acquisition module: used to obtain experimental data related to stomatal conductance and yield during the rice growing season; Stomatal conductance calculation module: used to parameterize the rice stomatal conductance model and calculate the hourly stomatal conductance of rice using hourly meteorological data; Stomatal O3 absorption flux calculation module: used to calculate the absorption flux of O3 by rice stomata under different thresholds; Dose response relationship fitting module: used to calculate the dose response relationship between rice yield and stomatal O3 absorption flux under the optimal threshold; Crop yield loss risk calculation module: used to obtain hourly meteorological data and O3 concentration data during the rice growth cycle in the region, calculate the rice stomatal O3 absorption flux, and quantify the rice yield loss risk caused by O3.

10. A device for simulating crop yield loss risk based on stomatal ozone absorption flux, characterized in that: It comprises a processor and a memory, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of a method for calculating the risk of crop yield loss based on quantification of stomatal ozone absorption flux as described in any one of claims 1 to 8.

Images