CN115144556B - Method and device for quantifying rainfall vapor source based on isotope and leaf area index - Google Patents
Method and device for quantifying rainfall vapor source based on isotope and leaf area index Download PDFInfo
- Publication number
- CN115144556B CN115144556B CN202210793231.7A CN202210793231A CN115144556B CN 115144556 B CN115144556 B CN 115144556B CN 202210793231 A CN202210793231 A CN 202210793231A CN 115144556 B CN115144556 B CN 115144556B
- Authority
- CN
- China
- Prior art keywords
- rainfall
- isotope
- vapor
- data
- water vapor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 88
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 173
- 230000005068 transpiration Effects 0.000 claims abstract description 82
- 238000001704 evaporation Methods 0.000 claims abstract description 63
- 230000008020 evaporation Effects 0.000 claims abstract description 62
- 230000008569 process Effects 0.000 claims abstract description 37
- 238000002156 mixing Methods 0.000 claims abstract description 25
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 90
- 239000001301 oxygen Substances 0.000 claims description 90
- 229910052760 oxygen Inorganic materials 0.000 claims description 90
- 238000004364 calculation method Methods 0.000 claims description 37
- 238000005194 fractionation Methods 0.000 claims description 31
- 239000002689 soil Substances 0.000 claims description 19
- 238000012545 processing Methods 0.000 claims description 13
- SQMPFRFYKISAHD-YFKPBYRVSA-N ethyl (2s)-2,5-diamino-5-oxopentanoate Chemical compound CCOC(=O)[C@@H](N)CCC(N)=O SQMPFRFYKISAHD-YFKPBYRVSA-N 0.000 claims description 12
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 8
- 239000001257 hydrogen Substances 0.000 claims description 8
- 229910052739 hydrogen Inorganic materials 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 7
- 210000000476 body water Anatomy 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 230000000087 stabilizing effect Effects 0.000 claims description 3
- 239000002352 surface water Substances 0.000 claims description 3
- 230000001970 hydrokinetic effect Effects 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000013461 design Methods 0.000 description 20
- 230000000694 effects Effects 0.000 description 12
- 238000004590 computer program Methods 0.000 description 7
- 230000000875 corresponding effect Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 2
- 230000005494 condensation Effects 0.000 description 2
- 238000009833 condensation Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 241000218631 Coniferophyta Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/28—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring areas
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Computational Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Evolutionary Biology (AREA)
- Food Science & Technology (AREA)
- Probability & Statistics with Applications (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention discloses a method and a device for quantifying rainfall vapor sources based on isotopes and leaf area indexes, wherein the method comprises the following steps: acquiring hydrological data and rainfall isotope data in a sample area for preprocessing, and respectively calculating the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor by utilizing the hydrological data and the rainfall isotope data; determining a fitting relation between the leaf area index and the evaporation and transpiration ratio; model constraints are set and a water source mixing model is constructed so that contributions of different water vapor sources to rainfall are calculated through the water source mixing model. According to the invention, three sources of rainfall vapor can be divided through one stable isotope, the dependence of quantitative hydrologic cycle on isotope data is reduced, the isotope detection cost required for distinguishing the rainfall vapor sources is greatly reduced, and meanwhile, the contribution of each vapor source to rainfall vapor is calculated, so that the time-space change characteristics of the regional hydrologic cycle process can be accurately reflected.
Description
Technical Field
The invention belongs to the technical field of isotope hydrology, and particularly relates to a method and a device for quantifying rainfall vapor sources based on isotopes and leaf area indexes.
Background
Hydrologic cycle is an important theoretical basis for researching water resource space-time distribution and change rule. Water vapor transport is the most important component of hydrologic cycle for regulating and controlling precipitation and water resource balance in different areas. The traditional method for quantifying the water vapor source in the atmosphere mainly comprises two methods: an analytical model and a numerical model. The analytical model is similar to a black box method, the calculation process is relatively simple, but the simulation process is too simple, so that the physical information provided by the analytical model is limited. The numerical model tracks the path of the water vapor by adopting a numerical method so as to simulate the related atmospheric movement process including convection effect, evapotranspiration effect, condensation effect and the like, but the method has the defects of huge calculation amount and uncertainty of physical process parameters.
Because of certain limitations of the two methods, a physical tracing method based on stable isotopes is provided at present, and mutual verification and supplementary explanation are carried out on the method from the perspective of experiments and the conclusions of an analytical model and a numerical model; the physical tracing method based on stable isotopes is based on water balance and isotope mass balance, and the stable isotopes are used as tracers to construct a water source mixing model so as to quantify the contribution of water vapor from different sources to rainfall. Because the isotope content in the rainfall is the comprehensive product of the historical motion path of the air mass and the local specific hydrological meteorological conditions in the process of generating the rainfall by the water vapor condensation, the difference of the isotope content in different water bodies reflects the physical processes of water circulation water vapor transportation, phase change and the like.
Physical tracing methods based on stable isotopes are often divided into three-terminal mixed models and two-terminal mixed models according to the number of sources of rainfall. The three-terminal mixed model generally assumes that rainfall vapor includes three sources of local vapor with evaporation, local vapor generated by transpiration and convection vapor, and two stable isotopes are needed as tracers to divide contributions of the three vapor sources. The two-end member mixed model only adopts one stable isotope, but can only distinguish two rainfall vapor sources. The traditional two-end member mixed model is mainly used for arid areas neglecting evaporation and is used for measuring the influence of evaporation of water bodies such as lakes on water vapor recycling and the cloud secondary evaporation effect. Because isotope data are very deficient, researchers have proposed an improved two-end member mixing model, and by introducing evaporation and transpiration ratio parameters, three rainfall vapor sources are divided by only adopting one stable isotope, and the dependence on the isotope data is greatly reduced. However, the ratio parameters in the model are the average results of the area for many years established in the prior study, and the time-space variation characteristics of the area water circulation process cannot be described in detail.
Therefore, how to provide an effective solution to describe the regional hydrologic cycle process in detail has become a technical problem to be solved in the prior art.
Disclosure of Invention
The invention aims to provide a method and a device for quantifying rainfall vapor sources based on isotopes and leaf area indexes, which are used for solving the technical problem that the regional hydrologic cycle process cannot be described in detail in the prior art.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
in a first aspect, the present invention provides a method for quantifying a rainfall vapor source based on isotopes and leaf area index, comprising:
acquiring hydrological data and rainfall isotope data in a sample area and preprocessing the data;
based on the characteristic that isotopes participate in the fractionation process in hydrologic cycle, respectively calculating the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor by utilizing the pretreated hydrologic data and rainfall isotope data, wherein the isotope comprises oxygen isotope or hydrogen isotope;
determining a fitting relation between the leaf area index and the evaporation and transpiration ratio according to vegetation characteristics in the sample area;
based on the water body balance, the isotope mass balance and the fitting relation between the leaf area index and the evaporation and transpiration ratio, model constraint conditions are set, and a water source mixing model is constructed so as to calculate the contribution of different water vapor sources to rainfall through the water source mixing model.
In one possible design, the hydrographic data includes at least rainfall capacity, isotope content in rainfall, temperature, humidity, rainfall capacity, site name, corresponding longitude and latitude, and time.
In one possible design, a method of preprocessing hydro-meteorological data and rainfall isotope data includes:
all the hydrological data and the rainfall isotope data are stored into a unified format, and interpolation processing is carried out on the data with different precision by utilizing an interpolation method, so that the hydrological data and the rainfall isotope data with the same grid precision are obtained.
In one possible design, when the isotope is an oxygen isotope, the content of the oxygen isotope in the rainfall vapor, the evaporation vapor, the transpiration vapor and the convection vapor is calculated respectively by using the pretreated hydrological data and the rainfall isotope data based on the characteristic that the isotope participates in the fractionation process in the hydrologic cycle, and the method comprises the following steps:
based on the characteristic that the oxygen isotope fractionation is equivalent to the equilibrium fractionation in the rainfall process, the content of oxygen isotopes in rainfall vapor is calculated, and the calculation formula is as follows:
wherein delta 18 O PV Represents the content of oxygen isotopes in rainfall vapor, PV represents rainfall vapor, delta 18 O P Represents the isotope content of rainfall, P represents rainfall,represents an equilibrium fractionation factor, w represents liquid water, v represents gaseous water;
calculating the content delta of oxygen isotopes in the evaporated water vapor based on Craig-Gordon model 18 O Ev The calculation formula is as follows:
wherein delta 18 O S Represents the oxygen isotope content, delta in surface water 18 O A Represents the oxygen isotope content in free atmosphere, h represents relative humidity, epsilon represents an enrichment factor determined by equilibrium fractionation and dynamic fractionation processes;
based on the characteristic that the oxygen isotope content in the transpiration water vapor is equivalent to the oxygen isotope content in the vegetation stem water, and the oxygen isotope content in the vegetation stem water is equivalent to the oxygen isotope content in soil, the content delta of the oxygen isotope in the transpiration water vapor is calculated 18 O Tr The calculation formula is as follows:
wherein delta 18 O Plant Represents the oxygen isotope content delta in vegetation stem water 18 O Soil The method comprises the steps of representing the oxygen isotope content in soil, wherein k and b respectively represent two parameters obtained by linear fitting of the oxygen isotope content in soil water and the oxygen isotope content in rainfall;
based on the Rayleigh fractionation characteristic, the oxygen isotope content in the convection water vapor is calculated according to the following calculation formula:
wherein delta 18 O i Adv Represents the oxygen isotope content, delta, of convective water vapor located in upwind zone 18 O i-1 PV Representing the oxygen isotope content of rainfall vapor in the downwind zone, F represents the rainfall ratio of the upwind zone to the downwind zone.
In one possible design, determining a fit between the leaf area index and the transpiration ratio from the vegetation signatures in the sample area includes:
according to vegetation types distributed in a sample area, determining a fitting relation between leaf area indexes and evaporation and transpiration ratio values under each vegetation type, wherein the fitting relation expression is as follows:
R T ＝aLAI c ； (5)
wherein R is T The ratio of the evaporation effect and the transpiration effect of vegetation is represented, a and c are coefficients of a fitting relation expression corresponding to the current vegetation type, and LAI represents leaf area index.
In one possible design, the model constraints include:
A. the rainfall water vapor comprises convection water vapor and local water vapor generated by the action of transpiration, wherein the local water vapor comprises water vapor generated by soil evaporation, water surface evaporation and vegetation transpiration;
B. the oxygen isotope content of each water vapor source is different;
C. within each sub-sample region, the ratio of transpiration to transpiration is significantly related to the leaf area index.
In one possible design, setting model constraints and constructing a water source mixing model based on water balance, isotope mass balance, and fitting relationships between leaf area index and evaporation transpiration ratio, comprising:
according to constraint condition A, a water balance formula is constructed as follows:
Q PV ＝Q ET +Q Adv ； (6)
wherein Q is PV Represents rainfall, Q ET Represents the local water vapor quantity, Q Adv Represents the amount of convective water vapor;
according to constraint condition B, a solute balance formula for stabilizing the oxygen isotope is constructed as follows:
δ 18 O PV ×Q PV ＝δ 18 O ET ×Q ET +δ 18 O Adv ×Q Adv ； (7)
wherein delta 18 O PV Represents the oxygen isotope content, delta in rainfall vapor 18 O ET Represents the oxygen isotope content, delta, in the local water vapor generated by the evapotranspiration 18 O Adv Representing the oxygen isotope content in convection water vapor;
according to formulas (6) and (7), a calculation formula of contribution of convection water vapor to rainfall water vapor is constructed as follows:
wherein F is Adv Representing the duty cycle contribution of convective water vapor volume in rainfall;
according to constraint condition C, a calculation formula of contribution of body water vapor generated by transpiration to rainfall water vapor is constructed, wherein the calculation formula is as follows:
δ 18 O ET ＝δ 18 O Ev (1-f(LAI))+δ 18 O Tr f(LAI)； (9)
where f (LAI) represents a fitted relationship of leaf area indexes.
In a second aspect, the present invention provides an apparatus for quantifying a source of rainfall vapor based on isotopes and leaf area index, comprising:
the data processing module is used for acquiring hydrological meteorological data and rainfall isotope data in the sample area and preprocessing the data;
the content calculation module is used for calculating the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor respectively by utilizing the preprocessed hydrological data and rainfall isotope data based on the characteristic that the isotope participates in the fractionation process in hydrologic cycle, wherein the isotope comprises oxygen isotope or hydrogen isotope;
the fitting module is used for determining a fitting relation between the leaf area index and the evaporation and transpiration action ratio according to vegetation characteristics in the sample area;
the model construction module is used for setting model constraint conditions and constructing a water source mixing model based on water body balance, isotope mass balance and fitting relation between leaf area index and evaporation and transpiration ratio so as to calculate contribution of different water vapor sources to rainfall through the water source mixing model.
In one possible design, the hydrographic data includes at least rainfall capacity, isotope content in rainfall, temperature, humidity, rainfall capacity, site name, corresponding longitude and latitude, and time.
In one possible design, the data processing module is specifically configured to:
all the hydrological data and the rainfall isotope data are stored into a unified format, and interpolation processing is carried out on the data with different precision by utilizing an interpolation method, so that the hydrological data and the rainfall isotope data with the same grid precision are obtained.
In a third aspect, the present invention provides a computer device comprising a memory, a processor and a transceiver in sequential communication, wherein the memory is adapted to store a computer program, the transceiver is adapted to receive and transmit messages, and the processor is adapted to read the computer program and perform a method of quantifying a source of rainfall vapor based on isotope and leaf area index as described in any one of the possible designs of the first aspect.
In a fourth aspect, the invention provides a computer readable storage medium having instructions stored thereon which, when executed on a computer, perform a method of quantifying a rainfall vapor source based on isotope and leaf area index as described in any one of the possible designs of the first aspect.
In a fifth aspect, the invention provides a computer program product comprising instructions which, when run on a computer, cause the computer to perform a method of quantifying a source of rainfall vapor based on isotope and leaf area index as described in any one of the possible designs of the first aspect.
Compared with the prior art, the invention has the beneficial effects that:
according to the invention, the hydrological weather data and the rainfall isotope data in the sample area are obtained and preprocessed, and the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor is calculated respectively by utilizing the preprocessed hydrological weather data and rainfall isotope data; then, according to vegetation characteristics in the sample area, determining a fitting relation between the leaf area index and the evaporation and transpiration ratio; and finally, setting model constraint conditions and constructing a water source mixing model based on the water body balance, the isotope mass balance and the fitting relation between the leaf area index and the evaporation and transpiration ratio so as to calculate the contribution of different water vapor sources to rainfall through the water source mixing model. The invention can divide three rainfall vapor sources through one stable isotope, reduces the dependence of quantitative hydrologic cycle process on isotope data, greatly reduces isotope detection cost required by distinguishing rainfall vapor sources, and simultaneously can accurately reflect the space-time variation characteristics of regional water cycle process by calculating the contribution of each vapor source to rainfall vapor.
Drawings
FIG. 1 is a flow chart of a method for quantifying a rainfall vapor source based on isotope and leaf area index provided by the present invention;
fig. 2 is a graph showing a fitted curve between leaf area index and transpiration by evaporation ratio provided by the invention.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the present invention will be briefly described below with reference to the accompanying drawings and the description of the embodiments or the prior art, and it is obvious that the following description of the structure of the drawings is only some embodiments of the present invention, and other drawings can be obtained according to these drawings without inventive effort to a person skilled in the art. It should be noted that the description of these examples is for aiding in understanding the present invention, but is not intended to limit the present invention.
Examples
In order to accurately reflect the space-time variation characteristics of the regional water circulation process so as to overcome the defects of the prior art, the embodiment of the application provides a method for quantifying rainfall vapor sources based on isotopes and leaf area indexes.
Methods for quantifying rainfall vapor sources based on isotopes and leaf area indices provided in embodiments of the present application are described in detail below.
It should be noted that, the method for quantifying rainfall vapor sources based on isotope and leaf area indexes provided in the embodiments of the present application may be applied to any terminal device for performing data calculation, where the terminal device includes, but is not limited to, an industrial computer, a personal computer, a mobile terminal, etc., and for convenience of description, unless otherwise specified, the embodiments of the present application all use the industrial computer as an execution body. It will be appreciated that the execution subject is not limited to the embodiments of the present application, and in other embodiments, other types of terminal devices may be used as the execution subject.
1-2, are flowcharts of methods for quantifying rainfall vapor sources based on isotopes and leaf area index provided in embodiments of the present application, including, but not limited to, by steps S1-S4:
s1, acquiring hydrological meteorological data and rainfall isotope data in a sample area and preprocessing the data;
in step S1, the hydrometeorological data at least includes a rainfall, an isotope content in rainfall, a temperature, a humidity, a rainfall, a site name, a corresponding longitude and latitude, and a time. The data sources comprise a global precipitation stable isotope monitoring network, a Chinese meteorological data network and the like, a plurality of pieces of hydrological meteorological data can be obtained from each data website, and then required data are extracted from different file formats by MATLAB software.
Because the data sources of different data may be different, and thus the specific data format and data accuracy are different, the data needs to be uniformly processed before being used for subsequent calculation, and therefore, in step S1, preferably, the method for preprocessing the hydrographic meteorological data and rainfall isotope data includes:
all the hydrological data and the rainfall isotope data are stored into a unified format, and interpolation processing is carried out on the data with different precision by utilizing an interpolation method, so that the hydrological data and the rainfall isotope data with the same grid precision are obtained.
Specifically, required data can be extracted from different file formats by MATLAB software and stored as a mat format, and then, interpolation processing is carried out on data with different line precision by using a cubic spline interpolation method to obtain hydrological meteorological data and rainfall isotope data with the same grid precision.
S2, based on the characteristic that isotopes participate in the fractionation process in hydrologic cycle, respectively calculating the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor by utilizing the preprocessed hydrologic vapor data and rainfall isotope data, wherein the isotope comprises oxygen isotopes or hydrogen isotopes;
it should be noted that, in this embodiment, the contribution of each vapor source to the rainfall vapor can be calculated by using the oxygen isotope content in the rainfall vapor, the evaporation vapor, the transpiration vapor and the convection vapor, and the contribution of each vapor source to the rainfall vapor can be calculated by using the hydrogen isotope content in the rainfall vapor, the evaporation vapor, the transpiration vapor and the convection vapor, the specific working principle is basically the same, but the calculation of the relevant parameters has slight differences, but it can be understood that the calculation method using the hydrogen isotope can be easily deduced according to the calculation method using the oxygen isotope, therefore, for simplifying the description, the embodiment is described in principle by using the oxygen isotope content hereinafter, and will not be repeated.
In step S2, when the isotope is an oxygen isotope, based on the characteristics of the isotope participating in the fractionation process in the hydrologic cycle, the content of the oxygen isotope in the rainfall vapor, the evaporation vapor, the transpiration vapor and the convection vapor is calculated respectively by using the preprocessed hydrologic data and the rainfall isotope data, and the method comprises the following steps:
s21, calculating the content of oxygen isotopes in rainfall vapor based on the characteristic that the oxygen isotope fractionation is equivalent to equilibrium fractionation in the rainfall process, wherein the calculation formula is as follows:
wherein delta 18 O PV Represents the content of oxygen isotopes in rainfall vapor, PV represents rainfall vapor, delta 18 O P Represents the isotope content of rainfall, P represents rainfall,represents an equilibrium fractionation factor, w represents liquid water, v represents gaseous water;
calculating the content delta of oxygen isotopes in the evaporated water vapor based on Craig-Gordon model 18 O Ev The calculation formula is as follows:
wherein delta 18 O S Represents the oxygen isotope content, delta in surface water 18 O A Represents the oxygen isotope content in free atmosphere, h represents relative humidity, epsilon represents an enrichment factor determined by equilibrium fractionation and dynamic fractionation processes;
the Craig-Gordon model divides isotope fraction occurring in the process of evaporating the open water body to the atmosphere into 3 stages, including equilibrium fractionation close to the boundary layer, dynamic fractionation of the middle viscous thin layer, and non-fractionation of the turbulent mixing thin layer. Based on this, a simplified formula (2) can be obtained.
S22, calculating the content delta of the oxygen isotope in the transpiration water vapor based on the characteristic that the oxygen isotope content in the transpiration water vapor is equivalent to the oxygen isotope content in the vegetation stem water and the oxygen isotope content in the soil 18 O Tr The calculation formula is as follows:
wherein delta 18 O Plant Represents the oxygen isotope content delta in vegetation stem water 18 O Soil The method comprises the steps of representing the oxygen isotope content in soil, wherein k and b respectively represent two parameters obtained by linear fitting of the oxygen isotope content in soil water and the oxygen isotope content in rainfall;
the isotope content in the local vapor of the transpiration vapor is the same as the isotope content in the vegetation stem water, and can be approximately regarded as the isotope content in the soil water, and the isotope content in the soil water is mainly subjected to rainfall (delta) 18 O P ) And long-term shallow groundwater, so that there is a certain linear correlation between the isotope content in the soil water and the isotope content in rainfall in the area. Based on which equation (3) can be obtained.
S23, calculating the oxygen isotope content in the convection water vapor based on the Rayleigh fractionation characteristic, wherein the calculation formula is as follows:
wherein delta 18 O i Adv Represents the oxygen isotope content, delta, of convective water vapor located in upwind zone 18 O i-1 PV Representing the oxygen isotope content of rainfall vapor in the downwind zone, F represents the rainfall ratio of the upwind zone to the downwind zone.
It should be noted that the oxygen isotope content in the convection vapor may be referred to as the rayleigh fractionation process, that is, the process in which the product generated by the reaction is separated from the system immediately after being formed in an open system, and thus the formula (4) may be obtained.
S3, as shown in fig. 2, determining a fitting relation between the leaf area index and the evaporation and transpiration ratio according to vegetation characteristics in the sample area;
in step S3, determining a fitting relationship between the leaf area index and the evaporation and transpiration ratio according to the vegetation features in the sample area, including:
according to vegetation types distributed in a sample area, determining a fitting relation between leaf area indexes and evaporation and transpiration ratio values under each vegetation type, wherein the fitting relation expression is as follows:
R T ＝aLAI c ； (5)
wherein R is T The ratio of the evaporation effect and the transpiration effect of vegetation is represented, a and c are coefficients of a fitting relation expression corresponding to the current vegetation type, and LAI represents leaf area index.
The specific vegetation type can be determined by performing field investigation or data collection on the vegetation data in the sample area. Because the vegetation transpiration is influenced by the canopy conductivity, the larger the canopy conductivity is, the larger the vegetation transpiration is; whereas canopy conductance is positively correlated with leaf area, therefore, localized transpiration by evaporation is significantly correlated with also leaf area. Based on this, according to different types of vegetation, in combination with historical data, MATLAB software can be used to establish a fit relationship between leaf area index and evaporation and transpiration ratio, such as the relationship between leaf area index and vegetation evaporation and transpiration ratio in the case of conifer, broadleaf, hybrid and shrub coverage shown in fig. 2. The leaf area index is derived from MCD15 leaf area data, the data is based on Terra-Aqua MODIS remote sensing earth surface reflectivity data, a lookup table is simulated and generated by using a three-dimensional radiation transmission model algorithm, and real leaf area data is obtained through inversion.
S4, setting model constraint conditions and constructing a water source mixing model based on water body balance, isotope mass balance and fitting relation between leaf area index and evaporation and transpiration ratio, so as to calculate contribution of different water vapor sources to rainfall through the water source mixing model.
In step s4, the model constraint condition includes:
A. the rainfall water vapor comprises convection water vapor and local water vapor generated by the action of transpiration, wherein the local water vapor comprises water vapor generated by soil evaporation, water surface evaporation and vegetation transpiration;
B. the oxygen isotope content of each water vapor source is different;
C. within each sub-sample region, the ratio of transpiration to transpiration is significantly related to the leaf area index.
In step S4, model constraints are set and a water source mixing model is constructed based on the water body balance, the isotope mass balance, and the fitting relationship between the leaf area index and the evaporation and transpiration ratio, including:
s41, constructing a water body balance formula according to the constraint condition A, wherein the water body balance formula is as follows:
Q PV ＝Q ET +Q Adv ； (6)
wherein Q is PV Represents rainfall water vapor quantity, Q ET Represents the local water vapor quantity, Q Adv Represents the amount of convective water vapor;
s42, constructing a solute balance formula of the stable oxygen isotope according to the constraint condition B, wherein the solute balance formula is as follows:
δ 18 O PV ×Q PV ＝δ 18 O ET ×Q ET +δ 18 O Adv ×Q Adv ； (7)
wherein delta 18 O PV Represents the oxygen isotope content, delta in rainfall vapor 18 O ET Represents the oxygen isotope content, delta, in the local water vapor generated by the evapotranspiration 18 O Adv Representing the oxygen isotope content in convection water vapor;
s43, constructing a calculation formula of contribution of convection water vapor to rainfall water vapor according to formulas (6) and (7), wherein the calculation formula is as follows:
wherein F is Adv Representing the duty cycle contribution of convective water vapor volume in rainfall;
wherein, in order to obtain the contribution of each water vapor source to rainfall water vapor, the method comprises the following steps ofFormulas (6) and (7) are simultaneously divided by Q PV Yield 1=f ET +F Adv And delta 18 O PV ＝F ET δ 18 O ET +F Adv δ 18 O Adv The two formulas are combined, the formula (8) can be obtained through deformation, and the oxygen isotope content in the rainfall vapor, the transpiration vapor and the convection vapor obtained through calculation in the step S2 is substituted into the formula (8), so that the solution can be carried out.
S44, constructing a calculation formula of contribution of body water vapor generated by transpiration to rainfall water vapor according to constraint condition C, wherein the calculation formula is as follows:
δ 18 O ET ＝δ 18 O Ev (1-f(LAI))+δ 18 O Tr f(LAI)； (9)
wherein f (LAI) represents a leaf area index relationship, i.e., formula (5).
The contributions of different rainfall sources can be obtained by putting the calculated oxygen isotope content in rainfall, evaporation, transpiration and convection vapor and the leaf area index fitting relations into calculation formulas (8) and (9) established based on the assumed conditions.
Based on the disclosure, the embodiment of the application calculates the content of the same isotope in the rainfall vapor, the evaporation vapor, the transpiration vapor and the convection vapor respectively by acquiring the hydrological data and the rainfall isotope data in the sample area and performing pretreatment and using the hydrological data and the rainfall isotope data after pretreatment; then, according to vegetation characteristics in the sample area, determining a fitting relation between the leaf area index and the evaporation and transpiration ratio; and finally, setting model constraint conditions and constructing a water source mixing model based on the water body balance, the isotope mass balance and the fitting relation between the leaf area index and the evaporation and transpiration ratio so as to calculate the contribution of different water vapor sources to rainfall through the water source mixing model. The invention can divide three rainfall vapor sources through one stable isotope, reduces the dependence of quantitative hydrologic cycle process on isotope data, greatly reduces isotope detection cost required by distinguishing rainfall vapor sources, and simultaneously can accurately reflect the space-time variation characteristics of regional water cycle process by calculating the contribution of each vapor source to rainfall vapor.
In a second aspect, the present invention provides an apparatus for quantifying a source of rainfall vapor based on isotopes and leaf area index, comprising:
the data processing module is used for acquiring hydrological meteorological data and rainfall isotope data in the sample area and preprocessing the data;
the content calculation module is used for calculating the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor respectively by utilizing the preprocessed hydrological data and rainfall isotope data based on the characteristic that the isotope participates in the fractionation process in hydrologic cycle, wherein the isotope comprises oxygen isotope or hydrogen isotope;
the fitting module is used for determining a fitting relation between the leaf area index and the evaporation and transpiration action ratio according to vegetation characteristics in the sample area;
the model construction module is used for setting model constraint conditions and constructing a water source mixing model based on water body balance, isotope mass balance and fitting relation between leaf area index and evaporation and transpiration ratio so as to calculate contribution of different water vapor sources to rainfall through the water source mixing model.
In one possible design, the hydrographic data includes at least rainfall capacity, isotope content in rainfall, temperature, humidity, rainfall capacity, site name, corresponding longitude and latitude, and time.
In one possible design, the data processing module is specifically configured to:
all the hydrological data and the rainfall isotope data are stored into a unified format, and interpolation processing is carried out on the data with different precision by utilizing an interpolation method, so that the hydrological data and the rainfall isotope data with the same grid precision are obtained.
The working process, working details and technical effects of the foregoing apparatus provided in the second aspect of the present embodiment may be referred to as the method described in the foregoing first aspect or any one of the possible designs of the first aspect, which are not described herein again.
In a third aspect, the present invention provides a computer device comprising a memory, a processor and a transceiver in sequential communication, wherein the memory is adapted to store a computer program, the transceiver is adapted to receive and transmit messages, and the processor is adapted to read the computer program and perform a method of quantifying a source of rainfall vapor based on isotope and leaf area index as described in any one of the possible designs of the first aspect.
By way of specific example, the Memory may include, but is not limited to, random-Access Memory (RAM), read-Only Memory (ROM), flash Memory (Flash Memory), first-in first-out Memory (First Input First Output, FIFO), and/or first-in last-out Memory (First Input Last Output, FILO), etc.; the processor may not be limited to use with a microprocessor of the STM32F105 family; the transceiver may be, but is not limited to, a WiFi (wireless fidelity) wireless transceiver, a bluetooth wireless transceiver, a GPRS (General Packet Radio Service, general packet radio service technology) wireless transceiver, and/or a ZigBee (ZigBee protocol, low power local area network protocol based on the ieee802.15.4 standard), etc. In addition, the computer device may include, but is not limited to, a power module, a display screen, and other necessary components.
The working process, working details and technical effects of the foregoing computer device provided in the third aspect of the present embodiment may be referred to the above first aspect or any one of the possible designs of the first aspect, which are not described herein.
In a fourth aspect, the invention provides a computer readable storage medium having instructions stored thereon which, when executed on a computer, perform a method of quantifying a rainfall vapor source based on isotope and leaf area index as described in any one of the possible designs of the first aspect.
The computer readable storage medium refers to a carrier for storing data, and may include, but is not limited to, a floppy disk, an optical disk, a hard disk, a flash Memory, and/or a Memory Stick (Memory Stick), etc., where the computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable devices.
The working process, working details and technical effects of the foregoing computer readable storage medium provided in the fourth aspect of the present embodiment may refer to the method as described in the foregoing first aspect or any one of the possible designs of the first aspect, which are not repeated herein.
In a fifth aspect, the invention provides a computer program product comprising instructions which, when run on a computer, cause the computer to perform a method of quantifying a source of rainfall vapor based on isotope and leaf area index as described in any one of the possible designs of the first aspect.
The working process, working details and technical effects of the foregoing computer program product containing instructions provided in the fifth aspect of the present embodiment may be referred to as the method described in the foregoing first aspect or any one of the possible designs of the first aspect, which are not repeated herein.
Finally, it should be noted that: the foregoing description is only of the preferred embodiments of the invention and is not intended to limit the scope of the invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (7)
1. A method for quantifying a source of rain water vapor based on isotopes and leaf area indices, comprising:
acquiring hydrological data and rainfall isotope data in a sample area and preprocessing the data;
based on the characteristic that isotopes participate in the fractionation process in hydrologic cycle, respectively calculating the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor by utilizing the pretreated hydrologic data and rainfall isotope data, wherein the isotope comprises oxygen isotope or hydrogen isotope;
determining a fitting relation between the leaf area index and the evaporation and transpiration ratio according to vegetation characteristics in the sample area;
setting model constraint conditions and constructing a water source mixing model based on water body balance, isotope mass balance and fitting relation between leaf area index and evaporation and transpiration ratio so as to calculate contribution of different water vapor sources to rainfall through the water source mixing model;
determining a fit relationship between leaf area index and evaporation and transpiration ratio according to vegetation features in a sample area, wherein the fit relationship comprises the following steps:
according to vegetation types distributed in a sample area, determining a fitting relation between leaf area indexes and evaporation and transpiration ratio values under each vegetation type, wherein the fitting relation expression is as follows:
R T ＝aLAI c ； (5)
wherein R is T A and c are coefficients of a fitting relation expression corresponding to the current vegetation type, and LAI represents leaf area index;
the model constraints include:
A. the rainfall water vapor comprises convection water vapor and local water vapor generated by the action of transpiration, wherein the local water vapor comprises water vapor generated by soil evaporation, water surface evaporation and vegetation transpiration;
B. the oxygen isotope content of each water vapor source is different;
C. within each sub-sample region, there is a significant correlation between the ratio of transpiration to transpiration and the leaf area index;
based on the water body balance, isotope mass balance and fitting relation between leaf area index and evaporation and transpiration ratio, setting model constraint conditions and constructing a water source mixing model, comprising:
according to constraint condition A, a water balance formula is constructed as follows:
Q PV ＝Q ET +Q Adv ； (6)
wherein Q is PV Represents rainfall, Q ET Represents the local water vapor quantity, Q Adv Represents the amount of convective water vapor;
according to constraint condition B, a solute balance formula for stabilizing the oxygen isotope is constructed as follows:
δ 18 O PV ×Q PV ＝δ 18 O ET ×Q ET +δ 18 O Adv ×Q Adv ； (7)
wherein delta 18 O PV Represents the oxygen isotope content, delta in rainfall vapor 18 O ET Represents the oxygen isotope content, delta, in the local water vapor generated by the evapotranspiration 18 O Adv Representing the oxygen isotope content in convection water vapor;
according to formulas (6) and (7), a calculation formula of contribution of convection water vapor to rainfall water vapor is constructed as follows:
wherein F is Adv Representing the duty cycle contribution of convective water vapor volume in rainfall;
according to constraint condition C, a calculation formula of contribution of body water vapor generated by transpiration to rainfall water vapor is constructed, wherein the calculation formula is as follows:
δ 18 O ET ＝δ 18 O Ev (1-f(LAI))+δ 18 O Tr f(LAI)；(9)
where f (LAI) represents a fitted relationship of leaf area indexes.
2. The method of quantifying a rainfall vapor source based on isotope and leaf area index of claim 1 wherein the hydrokinetic data comprises at least rainfall capacity, isotope content in rainfall, temperature, humidity, rainfall capacity, site name, corresponding longitude and latitude and time.
3. The method of quantifying a rain water vapor source based on isotopes and leaf area indices of claim 1, wherein the method of preprocessing the hydro-meteorological data and rain isotope data comprises:
all the hydrological data and the rainfall isotope data are stored into a unified format, and interpolation processing is carried out on the data with different precision by utilizing an interpolation method, so that the hydrological data and the rainfall isotope data with the same grid precision are obtained.
4. The method for quantifying rainfall vapor sources based on isotopes and leaf area indexes according to claim 1, wherein when the isotopes are oxygen isotopes, the contents of oxygen isotopes in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor are calculated respectively using the pre-processed hydrological data and rainfall isotope data based on the characteristics of the isotopes participating in the fractionation process in the hydrologic cycle, comprising:
based on the characteristic that the oxygen isotope fractionation is equivalent to the equilibrium fractionation in the rainfall process, the content of oxygen isotopes in rainfall vapor is calculated, and the calculation formula is as follows:
wherein delta 18 O PV Represents the content of oxygen isotopes in rainfall vapor, PV represents rainfall vapor, delta 18 O P Represents the isotope content of rainfall, P represents rainfall,represents an equilibrium fractionation factor, w represents liquid water, v represents gaseous water;
calculating the content delta of oxygen isotopes in the evaporated water vapor based on Craig-Gordon model 18 O Ev The calculation formula is as follows:
wherein delta 18 O S Represents the oxygen isotope content, delta in surface water 18 O A Represents the oxygen isotope content in free atmosphere, h represents relative humidity, epsilon represents an enrichment factor determined by equilibrium fractionation and dynamic fractionation processes;
the content of oxygen isotope in the vapor based on transpiration is equivalent to vegetationThe oxygen isotope content in the stem water is equivalent to the characteristic of the oxygen isotope content in the soil in the vegetation stem water, and the oxygen isotope content delta in the transpiration water vapor is calculated 18 O Tr The calculation formula is as follows:
wherein delta 18 O Plant Represents the oxygen isotope content delta in vegetation stem water 18 O Soil The method comprises the steps of representing the oxygen isotope content in soil, wherein k and b respectively represent two parameters obtained by linear fitting of the oxygen isotope content in soil water and the oxygen isotope content in rainfall;
based on the Rayleigh fractionation characteristic, the oxygen isotope content in the convection water vapor is calculated according to the following calculation formula:
wherein delta 18 O i Adv Represents the oxygen isotope content, delta, of convective water vapor located in upwind zone 18 O i-1 PV Representing the oxygen isotope content of rainfall vapor in the downwind zone, F represents the rainfall ratio of the upwind zone to the downwind zone.
5. A device for quantifying a source of rainfall vapor based on isotopes and leaf area indices, comprising:
the data processing module is used for acquiring hydrological meteorological data and rainfall isotope data in the sample area and preprocessing the data;
the content calculation module is used for calculating the content of the same isotope in rainfall vapor, evaporation vapor, transpiration vapor and convection vapor respectively by utilizing the preprocessed hydrological data and rainfall isotope data based on the characteristic that the isotope participates in the fractionation process in hydrologic cycle, wherein the isotope comprises oxygen isotope or hydrogen isotope;
the fitting module is used for determining a fitting relation between the leaf area index and the evaporation and transpiration action ratio according to vegetation characteristics in the sample area;
the model construction module is used for setting model constraint conditions and constructing a water source mixing model based on water body balance, isotope mass balance and fitting relation between leaf area index and evaporation and transpiration ratio so as to calculate contribution of different water vapor sources to rainfall through the water source mixing model;
determining a fit relationship between leaf area index and evaporation and transpiration ratio according to vegetation features in a sample area, wherein the fit relationship comprises the following steps:
according to vegetation types distributed in a sample area, determining a fitting relation between leaf area indexes and evaporation and transpiration ratio values under each vegetation type, wherein the fitting relation expression is as follows:
R T ＝aLAI c ； (5)
wherein R is T A and c are coefficients of a fitting relation expression corresponding to the current vegetation type, and LAI represents leaf area index;
the model constraints include:
A. the rainfall water vapor comprises convection water vapor and local water vapor generated by the action of transpiration, wherein the local water vapor comprises water vapor generated by soil evaporation, water surface evaporation and vegetation transpiration;
B. the oxygen isotope content of each water vapor source is different;
C. within each sub-sample region, there is a significant correlation between the ratio of transpiration to transpiration and the leaf area index;
based on the water body balance, isotope mass balance and fitting relation between leaf area index and evaporation and transpiration ratio, setting model constraint conditions and constructing a water source mixing model, comprising:
according to constraint condition A, a water balance formula is constructed as follows:
Q PV ＝Q ET +Q Adv ； (6)
wherein Q is PV Represents rainfall, Q ET Represents the local water vapor quantity, Q Adv Represents the amount of convective water vapor;
according to constraint condition B, a solute balance formula for stabilizing the oxygen isotope is constructed as follows:
δ 18 O PV ×Q PV ＝δ 18 O ET ×Q ET +δ 18 O Adv ×Q Adv ； (7)
wherein delta 18 O PV Represents the oxygen isotope content, delta in rainfall vapor 18 O ET Represents the oxygen isotope content, delta, in the local water vapor generated by the evapotranspiration 18 O Adv Representing the oxygen isotope content in convection water vapor;
according to formulas (6) and (7), a calculation formula of contribution of convection water vapor to rainfall water vapor is constructed as follows:
wherein F is Adv Representing the duty cycle contribution of convective water vapor volume in rainfall;
according to constraint condition C, a calculation formula of contribution of body water vapor generated by transpiration to rainfall water vapor is constructed, wherein the calculation formula is as follows:
δ 18 O ET ＝δ 18 O Ev (1-f(LAI))+δ 18 O Tr f(LAI)；(9)
where f (LAI) represents a fitted relationship of leaf area indexes.
6. The apparatus for quantifying a rain water vapor source based on isotope and leaf area index of claim 5 wherein the hydrographic data comprises at least rainfall capacity, isotope content in rainfall, temperature, humidity, rainfall capacity, site name, corresponding latitude and longitude, and time.
7. The apparatus for quantifying a source of rain water vapor based on isotopes and leaf area indices of claim 5, wherein said data processing module is specifically configured to:
all the hydrological data and the rainfall isotope data are stored into a unified format, and interpolation processing is carried out on the data with different precision by utilizing an interpolation method, so that the hydrological data and the rainfall isotope data with the same grid precision are obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210793231.7A CN115144556B (en) | 2022-07-05 | 2022-07-05 | Method and device for quantifying rainfall vapor source based on isotope and leaf area index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210793231.7A CN115144556B (en) | 2022-07-05 | 2022-07-05 | Method and device for quantifying rainfall vapor source based on isotope and leaf area index |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115144556A CN115144556A (en) | 2022-10-04 |
| CN115144556B true CN115144556B (en) | 2024-03-26 |
Family
ID=83411400
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210793231.7A Active CN115144556B (en) | 2022-07-05 | 2022-07-05 | Method and device for quantifying rainfall vapor source based on isotope and leaf area index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115144556B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116628414B (en) * | 2023-07-25 | 2023-09-26 | 中国科学院地理科学与资源研究所 | Simulation prediction method for water vapor sources inside and outside drought and waterlogging disaster frequent areas |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106885892A (en) * | 2017-02-22 | 2017-06-23 | 北京林业大学 | The method and device split to forest ecosystem evapotranspiration |
| CN109598082A (en) * | 2018-12-17 | 2019-04-09 | 河海大学 | The calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope |
| CN111783027A (en) * | 2020-06-16 | 2020-10-16 | 湖北省水利水电规划勘测设计院 | Method for calculating replenishment rate of rainfall and underground water to surface water by utilizing isotopes |
-
2022
- 2022-07-05 CN CN202210793231.7A patent/CN115144556B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106885892A (en) * | 2017-02-22 | 2017-06-23 | 北京林业大学 | The method and device split to forest ecosystem evapotranspiration |
| CN109598082A (en) * | 2018-12-17 | 2019-04-09 | 河海大学 | The calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope |
| CN111783027A (en) * | 2020-06-16 | 2020-10-16 | 湖北省水利水电规划勘测设计院 | Method for calculating replenishment rate of rainfall and underground water to surface water by utilizing isotopes |
Non-Patent Citations (2)
| Title |
|---|
| 基于LMDZ 模型的西南地区水汽来源及水汽再循环率的分析;武茜茜等;地球与环境;第49卷(第4期);400-408 * |
| 基于氢氧同位素技术的滇东南峰林湖盆区不同群落降水分配过程;刘珏杉等;四川农业大学学报;第35卷(第3期);333-338 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115144556A (en) | 2022-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zheng et al. | Improved estimate of global gross primary production for reproducing its long-term variation, 1982–2017 | |
| Ciais et al. | Five decades of northern land carbon uptake revealed by the interhemispheric CO2 gradient | |
| Mishra et al. | Moist heat stress extremes in India enhanced by irrigation | |
| Yeung et al. | Isotopic constraint on the twentieth-century increase in tropospheric ozone | |
| Zhang et al. | Multi-decadal trends in global terrestrial evapotranspiration and its components | |
| Foken | The energy balance closure problem: An overview | |
| Zhu et al. | Greening of the Earth and its drivers | |
| Pan et al. | Responses of global terrestrial evapotranspiration to climate change and increasing atmospheric CO2 in the 21st century | |
| Jasechko et al. | Terrestrial water fluxes dominated by transpiration | |
| Risi et al. | Water‐stable isotopes in the LMDZ4 general circulation model: Model evaluation for present‐day and past climates and applications to climatic interpretations of tropical isotopic records | |
| CN113297528B (en) | NO based on multi-source big data2High-resolution space-time distribution calculation method | |
| Feinberg et al. | Improved tropospheric and stratospheric sulfur cycle in the aerosol–chemistry–climate model SOCOL-AERv2 | |
| Monerie et al. | Expected future changes in the African monsoon between 2030 and 2070 using some CMIP3 and CMIP5 models under a medium‐low RCP scenario | |
| Cuntz et al. | A comprehensive global three‐dimensional model of δ18O in atmospheric CO2: 1. Validation of surface processes | |
| CN112884079A (en) | Method for estimating near-surface nitrogen dioxide concentration based on Stacking integrated model | |
| Huang et al. | Improving the global MODIS GPP model by optimizing parameters with FLUXNET data | |
| Feng et al. | Precipitation characteristic changes due to global warming in a high‐resolution (16 km) ECMWF simulation | |
| Júnior et al. | EVAPO: A smartphone application to estimate potential evapotranspiration using cloud gridded meteorological data from NASA-POWER system | |
| Liu et al. | Mapping the FTS SWIR L2 product of XCO2 and XCH4 data from the GOSAT by the Kriging method–a case study in East Asia | |
| Maksyutov et al. | NIES/FRCGC global atmospheric tracer transport model: Description, validation, and surface sources and sinks inversion | |
| Fujita et al. | Observed behaviours of precipitable water vapour and precipitation intensity in response to upper air profiles estimated from surface air temperature | |
| Keenan et al. | A constraint on historic growth in global photosynthesis due to rising CO2 | |
| CN115144556B (en) | Method and device for quantifying rainfall vapor source based on isotope and leaf area index | |
| Mohawesh et al. | Comparison of Hargreaves and FAO56 equations for estimating monthly evapotranspiration for semi-arid and arid environments | |
| Voigt et al. | Arctic soil methane sink increases with drier conditions and higher ecosystem respiration |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |