A method for determining a watershed vegetation and precipitation characteristic model based on evaporation force, an electronic device, and a storage medium

Abstract

The application belongs to the technical field of hydrology, and provides a watershed vegetation and precipitation characteristic model determination method based on evaporation force, an electronic device and a storage medium, which comprise data acquisition, evaporation force calculation model construction, evaporation force calculation model optimization, Gaussian interpolation processing, and evaporation force attenuation area and attenuation zone identification; the application can accurately invert the spatial distribution of evaporation force from reanalysis data such as temperature, humidity and air pressure by introducing the static force balance formula and the Clausius-Clapeyron equation, has good physical consistency and applicability; by fusing real-time leaf area index, the evaporation force model can respond to different underlying surface ecological characteristics, and the explanatory power of the model on the vertical water vapor transport mechanism is improved; by constructing the evaporation force field through Gaussian interpolation and Kriging algorithm, the evaporation force attenuation area is identified along the gradient direction, and then the precipitation spatial attenuation model is established, so that the response mechanism of precipitation on the watershed scale is effectively described.

Description

Technical Field

[0001] This invention relates to the field of hydrology technology, and in particular to a method for determining watershed vegetation and precipitation characteristics based on evaporation force, an electronic device, and a storage medium. Background Technology

[0002] Leaf area index (LAI) is a key indicator of climate change. In recent years, with the increasing frequency of global climate change and extreme weather events, a deeper understanding of the coupling mechanisms between vegetation distribution, water vapor transport, and regional precipitation is crucial for watershed ecological protection, water resource allocation, and drought early warning. Evapotranspiration, as an atmospheric dynamic variable driven by the vertical vapor pressure gradient caused by surface evaporation, effectively describes the energy and water vapor exchange processes between the underlying surface and the near-surface layer, reflecting the interaction characteristics between the Earth and the atmosphere.

[0003] However, traditional precipitation simulation and watershed hydrological feature identification methods largely rely on empirical parameters and statistical models, often failing to reveal the physical relationship between atmospheric thermodynamic processes and underlying vegetation characteristics. This is particularly problematic in areas with complex topography or significant variations, where accuracy is insufficient and adaptability is poor. Existing research mostly uses evapotranspiration as a theoretical variable, lacking systematic methods for its inversion, spatial interpolation, and quantitative coupling modeling with precipitation processes, thus hindering the full exploration of its application potential in vegetation and hydrological feature identification. Furthermore, although LAI, as an important parameter for measuring vegetation cover, has been widely used in ecological and evapotranspiration simulations, there is still a lack of methods to incorporate LAI into evapotranspiration-driven water vapor transport models to establish a quantitative mechanism for the regulatory role of vegetation on precipitation processes. Summary of the Invention

[0004] To overcome the shortcomings of existing technologies, the purpose of this invention is to provide a method, electronic device, and storage medium for determining watershed vegetation and precipitation characteristics based on evaporation force, which solves the problems of poor physical consistency and applicability of relying on empirical parameters and statistical models, inability to handle the physical relationship between atmospheric thermal processes and underlying vegetation characteristics, and lack of precipitation response mechanism at the watershed scale.

[0005] To achieve the above objectives, the present invention provides the following solution:

[0006] A method for determining watershed vegetation and precipitation characteristics based on evaporation force includes:

[0007] Data on temperature and humidity fields of various global pressure layers, as well as surface and sea surface temperature fields, are obtained to produce raw data.

[0008] An evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation.

[0009] The evaporation force calculation model was modified using the leaf area index to obtain an optimized evaporation force model.

[0010] The evaporation force optimization model is processed using Gaussian interpolation based on the original data to obtain an equal evaporation force field;

[0011] By identifying the evaporation force attenuation region and attenuation zone along the evaporation force gradient direction of the isoevaporation force field, a watershed precipitation attenuation characteristic model is obtained.

[0012] Preferably, data on the temperature and humidity fields of various global pressure layers, as well as surface and sea surface temperature fields, are acquired to obtain raw data, including:

[0013] Several pressure layers are defined; the pressure range of the pressure layers is from 100 hPa to 1000 hPa.

[0014] The air temperature and specific humidity data of each of the aforementioned pressure layers are collected to obtain the temperature and humidity field;

[0015] Collect surface temperature data and ocean surface temperature data to obtain the surface and ocean surface temperature field data;

[0016] The raw data is obtained by fusing the temperature and humidity field data with the surface and sea surface temperature field data.

[0017] Preferably, an evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation, including:

[0018] Extract the water vapor partial pressure and water vapor density from the raw data;

[0019] The static equilibrium formula is constructed based on the water vapor partial pressure and the water vapor density; the expression of the static equilibrium formula is: Where f is the output of the static equilibrium formula; z represents the partial pressure of the water vapor; z represents the altitude. The density of the water vapor is given. This is the gravity coefficient.

[0020] Preferably, an evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation, including:

[0021] Extract the air temperature from the raw data;

[0022] The Clausius-Clapeyron equation is constructed based on the air temperature; the expression of the Clausius-Clapeyron equation includes: ;in, The scale represents the vertical distribution of water vapor. The scale height is for dry air; T is the air temperature. R is the temperature constant; R is the gas constant; The molar mass of dry air.

[0023] Preferably, an evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation, including:

[0024] By combining the static equilibrium formula and the Clausius-Clapeyron equation, the evaporation force calculation model is obtained; the expression of the evaporation force calculation model is: ;in, ; This is the output of the evaporation force calculation model; This is the ratio of the dry air height to the water vapor height on the gauge. It is the molar mass ratio; The percentage of water vapor partial pressure; ρ is the molar mass of water vapor; p is the total atmospheric pressure.

[0025] Preferably, the evaporation force calculation model is modified using the leaf area index to obtain an optimized evaporation force model, including:

[0026] Extract the leaf area index from the original data;

[0027] The evaporation force optimization model is constructed based on the leaf area index and the evaporation force calculation model; the expression of the evaporation force optimization model is: ;in, The output of the evaporation force optimization model is denoted as ; LAI is the leaf area index.

[0028] Preferably, the evaporation force optimization model is processed using Gaussian interpolation based on the original data to obtain an iso-evaporation force field, including:

[0029] Define a Gaussian semivariogram model;

[0030] Set the search radius and generation precision;

[0031] Based on the search radius and the generation accuracy, the output of the evaporation force optimization model is processed using the Kriging interpolation algorithm and the Gaussian semi-variogram model to obtain the iso-evaporation force field.

[0032] Preferably, the evaporation force attenuation region and attenuation zone are identified along the evaporation force gradient direction of the isoevaporation force field to obtain a watershed precipitation attenuation characteristic model, including:

[0033] Identify the evaporation force attenuation region along the evaporation force gradient direction of the equal evaporation force field, and obtain the attenuation distance;

[0034] The attenuation zone is identified along the evaporation force gradient direction of the isoevaporation force field, thus obtaining the watershed precipitation attenuation characteristic model; the expression of the watershed precipitation attenuation characteristic model is: ;in, This is the output of the watershed precipitation attenuation characteristic model; denoted as precipitation at the starting point of the attenuation zone; x represents the distance from the starting point of the attenuation zone; and l represents the attenuation distance.

[0035] Preferably, an electronic device includes: at least one processor and a memory communicatively connected to the processor; wherein the memory stores instructions executable by the processor, the instructions being executed by the processor to enable the processor to perform the aforementioned method for determining a watershed vegetation and precipitation characteristic model based on evaporation force.

[0036] Preferably, a non-transitory computer-readable storage medium stores computer instructions for causing a computer to execute the aforementioned method for determining a watershed vegetation and precipitation characteristic model based on evaporation force.

[0037] The present invention discloses the following technical effects:

[0038] This invention provides a method, electronic device, and storage medium for determining watershed vegetation and precipitation characteristics based on evaporation force. By introducing the static equilibrium formula and the Clausius-Clapeyron equation, it solves the problems of poor physical consistency and applicability of models that rely heavily on empirical parameters and statistical models, and realizes the construction of an evaporation force calculation model based on multiple meteorological elements. By integrating real-time leaf area index, it solves the deficiency of existing methods in handling the physical relationship between atmospheric thermodynamic processes and underlying vegetation characteristics, and realizes nonlinear correction of vegetation regulation capabilities. By constructing an isoevaporation force field through Gaussian interpolation and Kriging algorithm, it solves the problem of existing methods lacking a precipitation response mechanism at the watershed scale, and realizes the identification of evaporation force attenuation regions along the gradient direction. Attached Figure Description

[0039] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0040] Figure 1 This is a schematic diagram illustrating the process of determining the watershed vegetation and precipitation characteristic model based on evaporation force, as provided in an embodiment of the present invention. Detailed Implementation

[0041] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0042] The purpose of this invention is to provide a method, electronic device, and storage medium for determining watershed vegetation and precipitation characteristics based on evaporation force, which solves the problems of poor physical consistency and applicability of relying on empirical parameters and statistical models, inability to handle the physical relationship between atmospheric thermodynamic processes and underlying vegetation characteristics, and lack of precipitation response mechanism at the watershed scale.

[0043] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0044] Figure 1 This is a schematic diagram illustrating the process of determining the watershed vegetation and precipitation characteristic model based on evaporation force, as provided in an embodiment of the present invention. Figure 1 As shown, this invention provides a method for determining watershed vegetation and precipitation characteristic models based on evaporation force, including:

[0045] Step 100: Obtain temperature and humidity fields of various global pressure layers, as well as surface and sea surface temperature fields, to obtain raw data;

[0046] Step 200: Construct an evaporation force calculation model based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation;

[0047] Step 300: The evaporation force calculation model is modified using the leaf area index to obtain the optimized evaporation force model;

[0048] Step 400: Process the evaporation force optimization model using the Gaussian interpolation method based on the original data to obtain an equal evaporation force field;

[0049] Step 500: Identify the evaporation force attenuation region and attenuation zone along the evaporation force gradient direction of the isoevaporation force field to obtain the watershed precipitation attenuation characteristic model.

[0050] Furthermore, data on the temperature and humidity fields of various global pressure layers, as well as surface and sea surface temperature fields, were acquired to obtain raw data, including:

[0051] Several pressure layers are defined; the pressure range of the pressure layers is from 100 hPa to 1000 hPa.

[0052] The air temperature and specific humidity data of each of the aforementioned pressure layers are collected to obtain the temperature and humidity field;

[0053] Collect surface temperature data and ocean surface temperature data to obtain the surface and ocean surface temperature field data;

[0054] The raw data is obtained by fusing the temperature and humidity field data with the surface and sea surface temperature field data.

[0055] Specifically, an evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation, including:

[0056] Extract the water vapor partial pressure and water vapor density from the raw data;

[0057] The static equilibrium formula is constructed based on the water vapor partial pressure and the water vapor density; the expression of the static equilibrium formula is: Where f is the output of the static equilibrium formula; z represents the partial pressure of the water vapor; z represents the altitude. The density of the water vapor is given. This is the gravity coefficient.

[0058] Furthermore, an evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation, including:

[0059] Extract the air temperature from the raw data;

[0060] The Clausius-Clapeyron equation is constructed based on the air temperature; the expression of the Clausius-Clapeyron equation includes: ;in, The scale represents the vertical distribution of water vapor. The scale height is for dry air; T is the air temperature. R is the temperature constant; R is the gas constant; The molar mass of dry air.

[0061] Specifically, an evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation, including:

[0062] By combining the static equilibrium formula and the Clausius-Clapeyron equation, the evaporation force calculation model is obtained; the expression of the evaporation force calculation model is: ;in, ; This is the output of the evaporation force calculation model; This is the ratio of the dry air height to the water vapor height on the gauge. It is the molar mass ratio; The percentage of water vapor partial pressure; ρ is the molar mass of water vapor; p is the total atmospheric pressure.

[0063] Furthermore, the evaporation force calculation model is modified using the leaf area index to obtain an optimized evaporation force model, including:

[0064] Extract the leaf area index from the original data;

[0065] The evaporation force optimization model is constructed based on the leaf area index and the evaporation force calculation model; the expression of the evaporation force optimization model is: ;in, The output of the evaporation force optimization model is denoted as ; LAI is the leaf area index.

[0066] Specifically, the evaporation force optimization model is processed using Gaussian interpolation based on the original data to obtain an iso-evaporation force field, including:

[0067] Define a Gaussian semivariogram model;

[0068] Set the search radius and generation precision;

[0069] Based on the search radius and the generation accuracy, the output of the evaporation force optimization model is processed using the Kriging interpolation algorithm and the Gaussian semi-variogram model to obtain the iso-evaporation force field.

[0070] Furthermore, by identifying evaporation force attenuation regions and attenuation zones along the evaporation force gradient direction of the isoevaporation force field, a watershed precipitation attenuation characteristic model is obtained, including:

[0071] Identify the evaporation force attenuation region along the evaporation force gradient direction of the equal evaporation force field, and obtain the attenuation distance;

[0072] The attenuation zone is identified along the evaporation force gradient direction of the isoevaporation force field, thus obtaining the watershed precipitation attenuation characteristic model; the expression of the watershed precipitation attenuation characteristic model is: ;in, This is the output of the watershed precipitation attenuation characteristic model; denoted as precipitation at the starting point of the attenuation zone; x represents the distance from the starting point of the attenuation zone; and l represents the attenuation distance.

[0073] Optionally, the global pressure layers refer to the 37 pressure layers from near the Earth's surface to the upper atmosphere, i.e., 100hPa-1000hPa; the global pressure layer temperature and humidity field data refer to the temperature and humidity data of each pressure layer; the surface and sea surface temperature field data refer to near-surface temperature data and ocean surface temperature data.

[0074] Specifically, the static equilibrium formula is: ;in, q represents the specific humidity; Let be the gas constant of water.

[0075] Furthermore, the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation can be expressed as: ; Furthermore, the constructed evaporation force calculation model is as follows: The total atmospheric pressure is taken from each pressure layer.

[0076] Optionally, leaf area index can be used for correction to ensure that the evaporation force model incorporates watershed vegetation characteristics, specifically: Wherein, LAI stands for Leaf Area Index (dimensionless). For evaporation force, It is the evaporative force corrected by the leaf area index.

[0077] Specifically, the isoevaporation force field is generated using Gaussian interpolation technology. This includes generating isoevaporation force lines using the Kriging interpolation algorithm, setting the semi-variogram model to Gaussian, setting the search radius to 50 km, and generating a 0.1°×0.1° grid. field.

[0078] Furthermore, based on the evaporation force gradient direction of the isoevaporation force field, the evaporation force attenuation region is identified, and the attenuation distance is determined. The attenuation zone is identified along the direction of maximum evaporation force gradient, resulting in a spatial attenuation characteristic model of precipitation. ;in, For precipitation; l is the attenuation length. denoted as precipitation at the starting point of the attenuation zone, and x as the distance from the starting point of the attenuation zone.

[0079] As an optional implementation, this embodiment also provides an electronic device, including: at least one processor and a memory communicatively connected to the processor; wherein the memory stores instructions executable by the processor, the instructions being executed by the processor to enable the processor to execute the aforementioned method for determining watershed vegetation and precipitation characteristics based on evaporation force.

[0080] As an optional implementation, this embodiment also provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the aforementioned method for determining watershed vegetation and precipitation characteristics based on evaporation force.

[0081] The beneficial effects of this invention are as follows:

[0082] This invention provides a method for determining watershed vegetation and hydrological characteristics based on evapotranspiration, overcoming the limitations of traditional empirical models in describing the spatial distribution of precipitation and enhancing the physical characterization of land-atmosphere interactions. By introducing the hydrostatic equilibrium formula and the Clausius-Clapeyron equation, an evapotranspiration calculation model based on multi-layered meteorological elements is constructed. This model can accurately invert the spatial distribution of evapotranspiration from reanalysis data such as temperature, humidity, and air pressure, exhibiting good physical consistency and applicability. Furthermore, by integrating the real-time leaf area index (LAI), a nonlinear correction to the vegetation regulation capacity is achieved, enabling the evapotranspiration model to respond to different underlying surface ecological characteristics and improving the model's explanatory power for vertical water vapor transport mechanisms. By constructing an isoevapotranspiration field using geostatistical Gaussian interpolation and the Kriging algorithm, evapotranspiration attenuation regions are identified along the gradient direction, thereby establishing a spatial precipitation attenuation model that effectively characterizes the response mechanism of precipitation at the watershed scale. This invention features high spatial resolution and a high degree of automation, making it suitable for identifying hydrological and ecological characteristics across multiple regions and time scales, and improving the accuracy of vegetation-hydrological coupling analysis.

[0083] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.

[0084] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for determining watershed vegetation and precipitation characteristics based on evaporation force, characterized in that, include: Data on temperature and humidity fields of various global pressure layers, as well as surface and sea surface temperature fields, are obtained to produce raw data. An evaporation force calculation model is constructed based on the static equilibrium formula and the atmospheric water vapor pressure gradient model of the Clausius-Clapeyron equation. The evaporation force calculation model was modified using the leaf area index to obtain an optimized evaporation force model. The evaporation force optimization model is processed using Gaussian interpolation based on the original data to obtain an equal evaporation force field; By identifying the evaporation force attenuation region and attenuation zone along the evaporation force gradient direction of the isoevaporation force field, a watershed precipitation attenuation characteristic model is obtained. An evaporation force calculation model is constructed based on the static equilibrium formula, the Clausius-Clapeyron equation, and the atmospheric water vapor pressure gradient model, including: Extract the water vapor partial pressure and water vapor density from the raw data; The static equilibrium formula is constructed based on the water vapor partial pressure and the water vapor density; the expression of the static equilibrium formula is: Where f is the output of the static equilibrium formula; z represents the partial pressure of the water vapor; z represents the altitude. The density of the water vapor is given. This is the gravity coefficient; An evaporation force calculation model is constructed based on the static equilibrium formula, the Clausius-Clapeyron equation, and the atmospheric water vapor pressure gradient model, including: Extract the air temperature from the raw data; The Clausius-Clapeyron equation is constructed based on the air temperature; the expression of the Clausius-Clapeyron equation includes: ;in, The scale represents the vertical distribution of water vapor. The scale height is for dry air; T is the air temperature. R is the temperature constant; R is the gas constant; Molar mass of dry air; An evaporation force calculation model is constructed based on the static equilibrium formula, the Clausius-Clapeyron equation, and the atmospheric water vapor pressure gradient model, including: By combining the static equilibrium formula and the Clausius-Clapeyron equation, the evaporation force calculation model is obtained; the expression of the evaporation force calculation model is: ;in, ; This is the output of the evaporation force calculation model; This is the ratio of the dry air height to the water vapor height on the gauge. It is the molar mass ratio; The percentage of water vapor partial pressure; ρ is the molar mass of water vapor; p is the total atmospheric pressure.

2. The method for determining the watershed vegetation and precipitation characteristics model based on evaporation force according to claim 1, characterized in that, Data on temperature and humidity fields of various global pressure layers, as well as surface and sea surface temperature fields, were obtained to produce raw data, including: Several pressure layers are defined; the pressure range of the pressure layers is from 100 hPa to 1000 hPa. The air temperature and specific humidity data of each of the aforementioned pressure layers are collected to obtain the temperature and humidity field; Collect surface temperature data of the Earth's surface and surface temperature data of the sea surface to obtain the surface and sea surface temperature field data; The raw data is obtained by fusing the temperature and humidity field data with the surface and sea surface temperature field data.

3. The method for determining watershed vegetation and precipitation characteristics based on evaporation force according to claim 1, characterized in that, The evaporation force calculation model is modified using the leaf area index to obtain an optimized evaporation force model, including: Extract the leaf area index from the original data; The evaporation force optimization model is constructed based on the leaf area index and the evaporation force calculation model; the expression of the evaporation force optimization model is: ;in, The output of the evaporation force optimization model; The leaf area index is mentioned.

4. The method for determining the watershed vegetation and precipitation characteristic model based on evaporation force according to claim 3, characterized in that, The evaporation force optimization model is processed using Gaussian interpolation based on the original data to obtain an iso-evaporation force field, including: Define a Gaussian semivariogram model; Set the search radius and generation precision; Based on the search radius and the generation accuracy, the output of the evaporation force optimization model is processed using the Kriging interpolation algorithm and the Gaussian semi-variogram model to obtain the iso-evaporation force field.

5. The method for determining the watershed vegetation and precipitation characteristic model based on evaporation force according to claim 4, characterized in that, Identifying evaporation force attenuation regions and attenuation zones along the evaporation force gradient direction of the isoevaporation force field yields a watershed precipitation attenuation characteristic model, including: Identify the evaporation force attenuation region along the evaporation force gradient direction of the equal evaporation force field, and obtain the attenuation distance; The attenuation zone is identified along the evaporation force gradient direction of the isoevaporation force field, thus obtaining the watershed precipitation attenuation characteristic model; the expression of the watershed precipitation attenuation characteristic model is: ;in, This is the output of the watershed precipitation attenuation characteristic model; denoted as precipitation at the starting point of the attenuation zone; x represents the distance from the starting point of the attenuation zone; and l represents the attenuation distance.

6. An electronic device, characterized in that, include: At least one processor, and a memory communicatively connected to the processor; wherein the memory stores instructions executable by the processor, the instructions being executed by the processor to enable the processor to perform a method for determining a watershed vegetation and precipitation characteristic model based on evaporation force, as described in any one of claims 1 to 5.

7. A non-transitory computer-readable storage medium storing computer instructions, characterized in that, The computer instructions are used to cause the computer to execute the method for determining the watershed vegetation and precipitation characteristics model based on evaporation force, as described in any one of claims 1 to 5.

Images